testedge.m

用拉普拉斯算子检测图像的阶梯状边缘,其中的模板可以进行修改

%***************************************************************************************************
% functionName:testedge.m
%     function: 用拉普拉斯算子检测阶梯状边缘
%         date:2006.05.08
%    progrmmer: wt
%***************************************************************

I=imread('e:\work\picture\ss.jpg');
I=rgb2gray(I);
%I=im2bw(I);
I=double(I);
[m n]=size(I);
row=zeros(1,n);
col=zeros(m+2,1);
I1=[row;I;row];
I1=[col,I1,col];
[m1 n1]=size(I1);
W=[ 0 -1  0;
   -1  4 -1;
    0 -1  0];
[p q]=size(W);
M=zeros(m1,n1);
for i=2:m1-2
    for j=2:n1-2
        I_block=I1(i-1:i+1,j-1:j+1);
        M_block=I_block.*W;
        R=sum(sum(M_block))/(p*q);
        M(i,j)=R;
    end
end
M=M(2:m1-1,2:n1-1);
%imshow(uint8(I),[]);
%figure,imshow(uint8(M),[]);
subplot(2,1,1),imshow(uint8(I),[]);
subplot(2,1,2),imshow(uint8(M),[]);