matlab练习程序（Top-Hat操作）

clear all;
close all;
clc;

img=imread('lena.jpg');
img=im2bw(img,graythresh(img));
imshow(img)

[height width]=size(img);

img_re=zeros(height,width);
temp=[];
for i=2:height-1                 %腐蚀
for j=2:width-1
temp=img(i-1:i+1,j-1:j+1);
img_re(i,j)=min(temp(:));
end
end

for i=2:height-1                %膨胀
for j=2:width-1
temp=img_re(i-1:i+1,j-1:j+1);
img_re(i,j)=max(temp(:));
end
end

img=double(img);
img_re=img-img_re;  %原图像减开操作后的图像,还有一种是闭操作后的图像减原图像
figure,imshow(mat2gray(img_re));