原文见https://blog.csdn.net/qq_23100787/article/details/51036396
I=imread('d:\\测试2.bmp');%读取原图像
figure(1);imshow(I);%显示原图像
I_gray=rgb2gray(I);%原图像变为灰度图像
level=graythresh(I_gray);%计算图像I_gray的全局阈值,level为标准化灰度值,其范围为[0 1]
[height,width]=size(I_gray);%计算灰度图像的长宽
I_bw=im2bw(I_gray,level);%im2bw使用阈值level将灰度图像转换为二值图像.
figure(2);imshow(I_bw);%显示二值图像(背景为白色)
for i=1:height %%循环中进行反色
for j=1:width
if I_bw(i,j)==1
I_bw(i,j)=0;
else I_bw(i,j)=1;
end
end
end
figure(3);imshow(I_bw);%显示取反后的二值图像(背景为黑色)
[L,num]=bwlabel(I_bw,8);%bwlabel标注二值图像I_bw中的目标物体,返回标识矩阵L和I_bw中目标物体的数量num,8表示连通数.
plot_x=zeros(1,num);%%zeros(m,n)产生m×n的全0矩阵.用于记录质心位置的横坐标
plot_y=zeros(1,num);%zeros(m,n)产生m×n的全0矩阵.用于记录质心位置的纵坐标
for k=1:num %%num个区域依次统计质心位置
sum_x=0; sum_y=0; area=0; %初始化
for i=1:height
for j=1:width
if L(i,j)==k
sum_x=sum_x+i; %计算第K区域的横坐标总和
sum_y=sum_y+j; %计算第K区域的纵坐标总和
area=area+1; %计算第K区域的由多少个坐标点表示
end
end
end
plot_x(k)=fix(sum_x/area); %计算第K区域的质心横坐标
plot_y(k)=fix(sum_y/area);%计算第K区域的质心纵坐标
end
figure(4);imshow(I_bw);%显示取反后的二值图像(背景为黑色),并在图上标注质心点位置
for i=1:num
hold on
plot(plot_y(i) ,plot_x(i), '*')
end
clear;clc;close all
%%读入图像
I_gray=imread('130.bmp');
level=graythresh(I_gray); %%求二值化的阈值
[height,width]=size(I_gray);
bw=im2bw(I_gray,level); %%二值化图像
figure(1),imshow(bw); %%显示二值化图像
[L,num]=bwlabel(bw,8); %%标注二进制图像中已连接的部分
plot_x=zeros(1,1); %%用于记录质心位置的坐标
plot_y=zeros(1,1);
%%求质心
sum_x=0;sum_y=0;area=0;
[height,width]=size(bw);
for i=1:height
for j=1:width
if L(i,j)==1
sum_x=sum_x+i;
sum_y=sum_y+j;
area=area+1;
end
end
end
%%质心坐标
plot_x(1)=fix(sum_x/area);
plot_y(1)=fix(sum_y/area);
figure(2);imshow(bw);
%%标记质心点
hold on
plot(plot_y(1) ,plot_x(1), '*')