Top-hat变换

clear all;
close all;%‘Fig2.22（a）.jpg‘
a=imread（‘C:\Documents and Settings\Administrator\桌面\49.bmp‘）;
[nmp]=size（a）;
figuresubplot（121）imshow（a）;title（‘原图像‘）;
a=double（a）;
b=zeros（19）;
e=zeros（nm）;
f=zeros（nm）;

for i=2:（n-1）
    for j=2:（m-1）
        b（11）=-1+a（i-1j-1）;
        b（12）=-2+a（i-1j）;
        b（13）=-3+a（i-1j+1）;
        b（14）=-4+a（ij-1）;
        b（15）=-5+a（ij）;
        b（16）=-6+a（ij+1）;
        b（17）=-7+a（i+1j-1）;
        b（18）=-8+a（i+1j）;
        b（19）=-9+a（i+1j+1）;
        min=255;
        for k=1:9         
            if b（1k）                min=b（1k）;
            end
        end
        e（ij）=min;
    end
end
for i=2:（n-1）
    for j=2:（m-1）
        b（11）=1+e（i-1j-1）;
        b（12）=2+e（i-1j）;
        b（1