计算连通域的面积、周长、圆形度、矩形度以及长宽比

二值化

最小矩

结果：

Ar =

8220 924 17191 93 0 5020 3133

perimeter =

162 54 262 54 118 112

metric =

3.9360 3.9819 3.1471

4.0078 4.5305 3.1 386

Pwl =

1.0816 1.0303 1.2720 1.0303 1.1467 1.2

500

Pr =

0.7913 0.8235 0.8650 0.8289 0.7783 0.5 420

>>

代码

clear,clc

RGB = imread('pillsetc.png');

figure,imshow(RGB);

I = rgb2gray(RGB);

threshold = graythresh(I);

bw = im2bw(I,threshold);

figure,imshow(bw)

% remove all object containing fewer than 30 pixels

bw = bwareaopen(bw,30);

% fill a gap in the pen's cap

se = strel('disk',2);

bw = imclose(bw,se);

% fill any holes, so that regionprops can be used to estimate

% the area enclosed by each of the boundaries

bw = imfill(bw,'holes');

figure,imshow(bw)

ed=edge(bw);

figure,imshow(ed)

%%%%%以上是图像二值化上面制作二值化图像

L = bwlabel(bw);

L1 = bwlabel(ed);

Ar=zeros(1,max(L(:)));

perimeter=zeros(1,max(L1(:)));

metric=zeros(1,max(L1(:)));

Pwl=zeros(1,max(L1(:)));

Pr=zeros(1,max(L1(:)));

for i=1:max(L(:))

Ar(i)=sum(bw(L==i));

perimeter(i)=sum(ed(L==i));

metric(i) = 4*pi*Ar(i)/perimeter(i)^2;

[y,x]=find(L==i);

x0=min(x(:));

x1=max(x(:));

y0=min(y(:));

y1=max(y(:));

hold on

rectangle('Position',[x0,y0,x1-x0,y1-y0],'edgeColor','g','LineWidth', 1)

if x1-x0>=y1-y0

Pwl(i)=(x1-x0)/(y1-y0);

else

Pwl(i)=(y1-y0)/(x1-x0);

end

Pr(i)=Ar(i)/((y1-y0)*(x1-x0));

end

Ar

perimeter

metric

Pwl

Pr

clear,clc

RGB = imread('D:\360°2è??ˉàà?÷????\tupian\1.png');

I = rgb2gray(RGB);

threshold = graythresh(I);

bw = im2bw(I,threshold);

bw = bwareaopen(bw,30);

se = strel('disk',2);

bw = imclose(bw,se);

bw = imfill(bw,'holes');

ed=edge(bw);

L = bwlabel(bw);

L1 = bwlabel(ed);

perimeter=zeros(1,max(L1(:)));

for i=1:max(L(:))

perimeter(i)=sum(ed(L==i));

[y,x]=find(L==i);

x0=min(x(:));

x1=max(x(:));

y0=min(y(:));

y1=max(y(:));

hold on

rectangle('Position',[x0,y0,x1-x0,y1-y0],'edgeColor','g','LineWidth', 1)

end

perimeter