图形图像处理

开发平台：

Matlab

priwittsem.m：源码内容

							I=imread('SEM.bmp');
I=im2double(I);
%figure;
%imshow(I);title('org img');
[height width R]=size(I);
for i=2:height-1
    for j=2:width-1
        Dx=[I(i+1,j-1)-I(i-1,j-1)]+[I(i+1,j)-I(i-1,j)]+[I(i+1,j+1)-I(i-1,j+1)];
        Dy=[I(i-1,j+1)-I(i-1,j-1)]+[I(i,j+1)-I(i,j-1)]+[I(i+1,j+1)-I(i+1,j-1)];
       P(i,j)=sqrt(Dx^2+Dy^2);
     
    end
end
%figure;
%imshow(P,[]);
for i=1:height-1
    for j=1:width-1
        if (P(i,j)<0.65)
            P(i,j)=1;
        else P(i,j)=0;
        end
    end
end
figure;
imshow(P,[]);

						