连通域的几何特征-程序

附录二

RGB = imread('F:/01054.png'); % 读入一幅彩色图像

figure,imshow(RGB); % 显示彩色图像

I = rgb2gray(RGB); % 转换为灰度图像

threshold = graythresh(I); %求阈值

bw 1= im2bw(I,threshold); %使用阈值变换法将灰度图转换成二值图像figure,imshow(bw1) %显示转化后的二值图像

Bw1 = bwareaopen(bw1,30); %删除二值图像bw中面积小于30的对象

se = strel('disk',2); %构造结构模板

bw 1= imclose(bw1,se); %闭运算，填补轮廓线

Bw1 = imfill(bw1,'holes'); %填充图像区域和“空洞”

figure,imshow(bw1)

bw=uint8(bw1) %数据类型转换

ed=edge(bw); % 边缘化

figure,imshow(ed) %ed为获得的边缘图像

L = bwlabel(bw); %贴标签

L1 = bwlabel(ed);

Ar=zeros(1,max(L(:))); %取L中所有元素,赋0值

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); %找出L中等于i的点的坐标

x0=min(x(:));

x1=max(x(:));

y0=min(y(:));

y1=max(y(:));

hold on %新画图像之后不想覆盖原图像就要加上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