sobel.m

来自「自编边缘提取方法」· M 代码 · 共 35 行

x=imread('E:\磁共振小波作业\翻译\老闫作业\边缘提取\mri1.jpg');
X0=rgb2gray(x);
X=double(X0);
S=size(X);
G=zeros(S);

%计算X的水平和垂直梯度
for i=2:S(1)-1
    for j=2:S(2)-1
        Dx=(X(i+1,j-1)-X(i-1,j-1)+2*(X(i+1,j)-X(i-1,j))+X(i+1,j+1)-X(i-1,j+1))^2; %Sobel算子的Dx(垂直梯度)
        Dy=(X(i-1,j+1)-X(i-1,j-1)+2*(X(i,j+1)-X(i,j-1))+X(i+1,j+1)-X(i+1,j-1))^2; %Sobel算子的Dx(水平梯度)
        G(i,j)=round(sqrt(Dx+Dy)); %此为梯度模,round为取整函数
    end
end 

%对梯度进行阈值处理
for i=1:S(1)
    for j=1:S(2)
        if G(i,j)>100  %把梯度值和阈值比较 ，大于T的把图像的灰度变为255小于T的把图像的灰度变为0
            X(i,j)=255; 
        else
            X(i,j)=0;
        end
    end
end

subplot(1,2,1);
imshow(X0);
title('原始图像')

subplot(1,2,2);
imshow(X);
colormap(gray(256));
title('经过阈值处理的水平+垂直边缘');