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

计算连通域的面积、周长、圆形度、矩形度以及长宽比 1 计算连通域的面积、周长、圆形度、矩形度以及长宽比 二值化 最小矩 结果： 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 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 pens 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)); 3 计算连通域的面积、周长、圆形度、矩形度以及长宽比 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°²È«ä¯ÀÀÆ÷ÏÂÔØ\t