% sharpening di immagine con filtro sobel

immagine = '43552_det.jpg';
filter = [1 2 1; 0 0 0; -1 -2 -1];

% oppure:
% filter = fspecial('sobel');

I = imread(immagine);
T = rgb2gray(I);

X1 = imfilter(double(T), filter); % derivata orizzontale
X2 = imfilter(double(T), filter'); % derivata verticale

Y = uint8(double(T)+(X1+X2)/4);

figure; imshow(Y); % sharpened
figure; imshow(T)  % original
