clear;
A=[111;111;111];
[A_heightA_width]=size（A）;
B=padarray（zeros（A_heightA_width）[1?1]1‘both‘）;???%結(jié)構(gòu)元素
I0=imread（‘bird2.jpg‘）;
I=rgb2gray（I0）;
imshow（I[]）;title（‘源圖像‘）;
I=double（I）;
[im_heightim_width]=size（I）;
[B_heightB_width]=size（B）;
A_halfheight=floor（A_height/2）;
A_halfwidth=floor（A_width/2）;
B_halfheight=floor（B_height/2）;
B_halfwidth=floor（B_width/2）;
[se_origin]=floor（（size（B）+1）/2）;
E=padarray（Ise_origin255‘both‘）;?%用于腐蝕
D=padarray（Ise_origin0‘both‘）;?%用于膨脹


for?m1=1:A_height*A_width
????E1（::m1）=E;
end
for?m2=1:B_height*B_width
????E2（::m2）=E;????????????????????%多層矩陣用于腐蝕
end


for?m1=1:A_height*A_width
????D1（::m1）=D;
end
for?m2=1:B_height*B_width
????D2（::m2）=D;????????????????????%多層矩陣用于膨脹
end


%腐蝕運(yùn)算
k1=1

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????11516??2012-05-10?17:10??bird.jpg

?????文件???????5310??2012-07-13?01:14??Soft_Morphology_Top_Hat.m

-----------?---------??----------?-----??----

????????????????16826????????????????????2