第三次数字图像处理实验.docVIP

%用Laplacian算子进行静态图像边缘检测 clear all; f=imread(e:\tupian\hui.jpg); imshow(f); title(原图); [m,n]=size(f); g=zeros(m,n); for i=2:m-1 for j=2:n-1 g(i,j)=4*f(i,j)-f(i-1,j)-f(i,j-1)-f(i,j+1)-f(i+1,j); end end figure, imshow(g); title(用Laplacian算子进行静态图像边缘检测后的图片) %用Sobel算子进行静态图像边缘检测 clear all; f=imread(E:\tupian\hui.jpg); imshow(f); title(原图); [m,n]=size(f); g=zeros(m,n); for i=2:m-1 for j=2:n-1 h=double(f(i+1,j-1)+2*f(i+1,j)+f(i+1,j+1)-f(i-1,j-1)-2*f(i-1,j)-f(i-1,j+1)); k=double(2*f(i,j+1)+f(i-1,j+1)+f(i+1,j+1)-f(i-1,j-1)-2*f(i,j-1)-f(i+1,j-1)); l=h*h;