bw=imread（‘223.bmp‘）;
bw=imcomplement（bw）;??????????????%反白
bw=im2bw（bwgraythresh（bw））;??????%轉化為灰度圖
imshow（bw）;

[BL]?=?bwboundaries（bw‘noholes‘）;??%檢測邊界
figureimshow（bw）;

%?Display?the?label?matrix?and?draw?each?boundary??%填充顏色
figureimshow（label2rgb（L?@jet?[.5?.5?.5]））
hold?on
for?k?=?1:length（B）
??boundary?=?B{k};
??plot（boundary（:2）?boundary（:1）?‘w‘?‘LineWidth‘?2）
end


stats?=?regionprops（L‘Area‘‘Centroid‘）;
threshold?=?0.8;????????????????????????????????%閾值設定

%?loop?over?the?boundaries????????????????????%計算面積與周長關系“metric
for?k?=?1:length（B）

??%?obtain?（XY）?boundary?coordinates?corresponding?to?label?‘k‘
??boundary?=?B{k};

??%?compute?a?simple?estim