www.gusucode.com > 图像变换及频域滤波 空间域图像增强 图像边缘检测源码程序 > digitalimageprocessprograms/ruihua.m

    clc;
fid=fopen('d:/img/lena.img','r');
f=fread(fid,[256,256],'uchar');
subplot(3,1,1);
imshow(f,[0,255]);
g1=f;g2=f;
a=1;
for x=2:255
    for y=2:255
        g1(x,y)=f(x,y)-(f((x-1),y)+f((x+1),y)+f(x,(y-1))+f(x,(y+1))-f(x,y).*4).*a;
        g2(x,y)=f(x,y).*4*a-(f((x-1),y)+f((x+1),y)+f(x,(y-1))+f(x,(y+1))).*a;
    end
end
subplot(3,1,2);
imshow(g1,[0,255]);
subplot(3,1,3);
imshow(g2,[0,255]);