clc
clear
open（‘截圖.fig‘）;
h=get（gca‘Children‘）;
I=get（h‘Cdata‘）;%導入fig圖片
Ig?=?rgb2gray（I）;
%T=?graythresh（Ig）;

Bimg?=?~（im2bw（Ig0.2））;?%二值化然后取反
L?=?bwlabel（Bimg）;???????%給每個區域加標簽
for?i?=?1?:?max（L（:））
????NrLi?=?sum（L==i）;?????%第i個區域的面積。
????if?NrLi?>1?????????%這個數你自己算，估計一下那些小孔最大能占多少個pixel。

????????L（L==i）?=?10;?????????%太大的不是你要的孔，去掉
????end
end
figure;?imshow（L）;
axis?on
[mn]=?size（L）;
%[r?c]?=?find（L?==?1）%c對應水平位置、r對應垂直位置
%[rowcol]?=?find（L）%查詢滿足一定條件的元素的行和列
%%%%%%%%%%%上
x=0;
y=0;
for?i=3:m
????for?j=3:n
????????if?（（Bimg（ij）==0）&&（Bimg（ij）==Bimg（i+1j））&&（Bimg（ij-1）==1）&&（Bimg（ij-1）==Bimg（ij+1））&&（Bimg（ij+1）==Bimg（i-1j））||（（Bimg（ij）==0）&&（Bimg（ij-1）==0）&&（Bimg（ij-2）==0）&&（Bimg（ij+1）==0）&&（Bimg（ij+2）==0）&&（Bimg（i+1j）==0）&&（Bimg（i-1j）==1）））
????????????x=x+1;
???????????a（1x）=5.82-（i+8）*3.04/369+0.118;
?????????????y=y+1;
???????????a（2y）=j;
????????end
????end
end
%%%%%%%%%%%%%%%中上
x1=0;
y1=0;
for?i=3:m
????for?j=3:n
????????