matlab 边缘检测程序
=imread('rice.png');
subplot(2,2,1);
imshow(I);
J = imnoise(I,'gaussian');;
subplot(2,2,2);
imshow(J);
[m,n]=size(J);
K=J;
h1=[1,2,1;2,4,2;1,2,1]/16;
h2=[1,1,1;1,1,1;1,1,1]/9;
for i=2:m-1;
for j=2:n-1
J(i,j)=h1(1)*J(i-1,j-1)+h1(2)*J(i-1,j)+h1(3)*J(i-1,j+1)+h1(4)*J(i,j-1)+h1(5)*J(i,j)+h1(6)*J(i,j+1)+h1(7)*J(i+1,j-1)+h1(8)*J(i+1,j)+h1(9)*J(i+1,j+1);
if J(i,j)>255
J(i,j)=255;
elseif J(i,j)<0
J(i,j)=abs(J(i,j));
else
J(i,j)=double(J(i,j));
end
end
end
subplot(2,2,3);
imshow(J);
for i=2:m-1;