计算图像的质心

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中目标物体