代码是这个样子
A=imread('pic2.png'); %读入要处理的图像
figure(1)
imshow(A)
X=rgb2gray(A);
[m,n]=size(X);
S=size(X); %图像的大小(长和宽)
for i=2:m-1 %进行边界提取
for j=2:n-1
Dx=abs(X(i+1,j-1)-X(i-1,j-1)+2*(X(i+1,j)-2*X(i-1,j))+X(i+1,j+1)-X(i-1,j+1)); %Sobel算子的Dx(垂直梯度)
Dy=abs(X(i-1,j+1)-X(i-1,j-1)+2*(X(i,j+1)-2*X(i,j-1))+X(i+1,j+1)-X(i+1,j-1)); %Sobel算子的Dx(水平梯度)
g(i,j)=Dx+Dy;
end
end
figure(2)
imshow(uint8(g))