Thursday, March 12, 2009

Laplacion filter

% Read a gray scale image and apply a 3 x 3 Laplacion filter
close all; clear all; clc;
a = imread('rice.png');
a=im2double(a);
%a=imcrop(a);
w = [-1 -1 -1; -1 8 -1; -1 -1 -1];
[m n] = size(a);
si=1;
cs=si+1;
for i=1:m
for j=1:n
con=0;s=0;
for k= -si:si
for p= -si:si
if (i+k > 0 && j+p > 0) && (i+k < m && j+p < n )
we=w(cs+k,cs+p);
if we < 0
w1=we*-1;
va= -1;
else
w1=we;
va= 1;
end
for wi=1:w1
% con=con+1;
s=s+va*a(i+k,j+p);
end
end
end
end
b(i,j)=s;
% clear s;
end
end
figure,imshow(a),title('Original Image');
figure,imshow(b),title('Image After Sharpaning');

No comments: