function?[plot_xplot_yheightwidth]=Centroid（I?）

SE1=[1?1?1?1?1?1?1?1?1];
SE2=[0?0?0?0?1?0?0?0?0;
????0?0?0?1?1?1?0?0?0;
????0?0?1?1?1?1?1?0?0;
????0?1?1?1?1?1?1?1?0;
????1?1?1?1?1?1?1?1?1;
????0?1?1?1?1?1?1?1?0;
????0?0?1?1?1?1?1?0?0;
????0?0?0?1?1?1?0?0?0;
????0?0?0?0?1?0?0?0?0];
I=imadjust（I）;
%?I=medfilt2（I[10?10]）?;???????????????%中值濾波
%??w2=fspecial（‘average‘[5?5]）;?%%?先定義一個濾波器
%?I=imfilter（Iw2‘replicate‘）;?%%讓圖像通過濾波器
fmax1=double（max（max（I）））;%I的最大值并輸出雙精度型
fmin1=double（min（min（I）））;%I的最小值并輸出雙精度型
level=（fmax1-（fmax1-fmin1）/3）/255;%獲得最佳閾值
I=im2bw（Ilevel）;
%?I=~I;
I=imdilate（ISE1）;%膨脹
%?I=imdilate（ISE1‘）;%膨脹
%?I=imopen（ISE2）;
I=medfilt2（I[10?10]）?;???????????????%中值濾波
I?=?bwareaopen（I?20）;?????????????????%刪除二值圖像中面積小于30的區域
%?imshow（I）;%顯示原圖像
[heightwidth]=size（I）;%計算灰度圖像的長寬
%?[Lnum]=bwlabel（I8）;%bwlabel標注二值圖像I_bw中的目標物體，返回標識矩陣Ｌ和I_bw中目標物體