图像处理——周长、面积计算(Matlab)

题目：面积与周长的计算检测下图细胞数量，并计算面积与周长的比。处理过程首先将图像变换为灰度图像，并取其边缘二值图像。然后分别对边缘二值掩膜进行膨胀、填充空洞、清除边缘。最后对处理后的二值图像的连通区域进行检测，求其连通区域个数、周长和面积，并以此求出面积与周长的比。处理结果最后得到细胞数量a=88，面积m=24333，周长z=4956.4，面积与周长的比b=4.9094。程序清单clc;clearallI=imread('7.bmp');figure;imshow(I);I2=rgb2gray(I);[junk,threshold]=edge(I2,'sobel');fudgeFactor=.5;BWs=edge(I2,'sobel',threshold*fudgeFactor);figure;subplot(221),imshow(BWs),title('边缘梯度二值掩膜');se90=strel('line',3,90);se0=strel('line',3,0);BWsdil=imdilate(BWs,[se90se0]);subplot(222);imshow(BWsdil),title('膨胀梯度掩膜');BWdfill=imfill(BWsdil,'holes');subplot(223);imshow(BWdfill);title('填充空洞后的二值图像');BWnobord=imclearborder(BWdfill,4);subplot(224);imshow(BWnobord),title('清除边缘的二值图像');k1=bwlabel(BWnobord);I5=~BWnobord;figure;imshow(I5);a=max(max(k1))[labeled,numObjects]=bwlabel(BWnobord,4);celldata=regionprops(labeled,'all');fori=1:1:acelldata(i).Areacelldata(i).Perimeter;endallcellm=[celldata.Area];allcellz=[celldata.Perimeter];m=sum(allcellm)z=sum(allcellz)b=m/z