%%?function?最終腐蝕函數
%?將二值圖像腐蝕成互不相連通的像素點
%%
function?[imgn]=dilation（img）
[m?n]=size（img）;
imgn=zeros（mn）;
preimg=imgn;
%腐蝕矩陣
se=[1?1?1;1?1?1;?1?1?1];
while?sum（sum（preimg-img））~=0?????%當腐蝕結果不變時跳出
????preimg=img;
????img=img>0;
???
????[imglabel]=bwlabel（img）;??????%標記不同區域，label是區域個數
????[imgn]=imerode（imgse）;
????
????%腐蝕之后是否有哪個被標記的區域消失了
????Hist=zeros（1label）;????????????
????for?i=1:m
????????for?j=1:n
????????????if?imgn（ij）~=0
????????????????Hist（imgn（ij））=imgn（ij）;??
????????????end
????????end
????end

????%統計消失區域的標號
????H=[];
????for?i=1:label
????????if?Hist（i）==0
????????????H=[H?i];???????
????????end
????end
????
????%如果這個區域消失了，那么再把這個區域恢復過來
????if?~isempty（H）
????????l=length（H）;
????????for?i=1:m
????????????for?j=1:n
??????????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1122??2017-06-04?18:55??dilation.m
?????文件??????455566??2017-06-04?05:26??Fig1039（a）（polymersomes）.tif
?????文件?????????345??2017-06-04?19:05??Region.m
?????文件????????1206??2017-06-04?19:01??RegionGrow.m
?????文件?????????861??2017-06-04?18:53??RegionGrowing.m
?????文件???????61925??2017-06-04?18:46??RegionGrowing.png