matlab函数提取遥感图像的边缘_内含原理说明文档_matlab_遥感图边缘检测

clc %分别通过Sobel边缘检测的计算公式及Edge函数中Sobel参数的设定来提取图片边缘 f0=imread('shida.jpg'); f1=rgb2gray(f0); f=double(f1); %转为double型 [m,n]=size(f); %Sobel检测的计算公式(dxs、dys分别是近似一阶导数值的掩模) for i=2:m-1 for j=2:n-1 ay(i,j)=f(i-1,j+1)+2*f(i,j+1)+f(i+1,j+1)-f(i-1,j-1)-2*f(i,j-1)-f(i+1,j-1); ax(i,j)=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); f2(i,j)=sqrt((ax(i,j)^2+ay(i,j)^2)); end end subplot(3,3,1) imshow(f0); title('原图像') subplot(3,3,2) imshow(f1); title('原图像的灰度图') subplot(3,3,3) imhist(f1); title('直方图') subplot(3,3,4) imshow(f2,[]); title('Sobel检测计算公式所得图') subplot(3,3,5) imshow(edge(f1,'sobel',0.105)); title('Edge函数Sobel') subplot(3,3,6) imshow(edge(f1,'canny')),title('Edge函数Canny'); subplot(3,3,7) imshow(edge(f1,'roberts')),title('Edge函数Roberts'); subplot(3,3,8) imshow(edge(f1,'prewitt')),title('Edge函数Prewitt');