I =imread('erciyuan.jpg');
img =rgb2gray(I);
img =im2double(img);[m,n]=size(img);for i=3:m-2for j=3:n-2% LoG算子粗提取图像边缘
log_margin(i,j)=-img(i-2,j)-img(i-1,j-1)-2*img(i-1,j)-img(i-1,j+1)-img(i,j+2)-2*img(i,j-1)+16*img(i,j)-2*img(i,j+1)-img(i,j+2)-img(i+1,j-1)-2*img(i+1,j)-img(i+1,j+1)-img(i+2,j);
end
end
subplot(1,2,1);imshow(I);title('原图');subplot(1,2,2);imshow(log_margin);title('显示图像边缘');