基于MATLAB的数字图像处理_规则图形平行四边形面积测量

1、数字图像处理_规则图形平行四边形面积测量%-平行四边形图形测量源程序%-作者:ZH-%-单位:HX-%-日期:JULY-%-步骤1*图像预处理*-image = imread('image_name.bmp' %读入原始图像image1 = imresize(image,0.2; %调整原始图像大小缩小至它的0.2倍x y z = size(image1figure,imshow(image1;title('图1 原始图像'%显示调整后的原始图像(如图1所示%-步骤2*图像灰度变换*-image2=rgb2gray(image1; %将调整后的原始图像转换为灰度

2、图figureimshow(image2; %显示调整后的原始图像的灰度图(如图2所示title('图2 灰度图像'%-步骤3*图像分割*-image = imread('image_name.bmp' %输入原始图像image1 = imread('Background2.bmp' %输入背景图像grayscale = rgb2gray(image; %灰度转换grayscale1 = rgb2gray(image1;subtract=imsubtract(grayscale1,grayscale; %差分运算去除背景figureimshow(

3、subtract; %显示输出图像即差分图像(如图3所示title('图 3差分图像'%-%coon=4;%noborder=imclearborder(image,coon; %用imclearborder命令消除图像中的边缘对象%-R=imclearborder(subtract,8;figureimshow(R; %显示去除边缘对象的图像(如图4所示title('图4 去除边缘对象的图像'%-步骤4*调节图像对比度*-% image=imadjust(image,0.2 0.6,0 1;%image=imadjust(R,0.2 0.6,0 1;figur

4、eimshow(image; %显示调节对比度后的图像(如图5所示title('图5 调节对比度后的图像'%-步骤5*边缘检测*-image1=edge(image,'sobel'image2=edge(image,'roberts'image3=edge(image,'prewit'image4=edge(image,'canny'figure;subplot(2,2,1;imshow(image1;title('sobel图像边缘提取'subplot(2,2,2;imshow(image2;ti

5、tle('roberts图像边缘提取'subplot(2,2,3;imshow(image3;title('prewit图像边缘提取'subplot(2,2,4;imshow(image4;title('canny图像边缘提取' %图6 显示四种算子的边缘提取图%title('图6 四种算子的边缘提取图'%-步骤6*图像形态学处理*-%-1 填补缝隙-I=strel('line',8,90;H=strel('line',6,0;image=imdilate(image2,I H;figure;ims

6、how(image; %显示填补边缘缝隙后的图像(如图7所示title('图7 填补边缘缝隙后的图像'%-2 填充-filled = imfill(image,'holes'figure,imshow(filled; %显示填充后的图像(如图8所示title('图8 填充后的图像'%-3 平滑-SeD=strel('diamond',1;final=imerode(filled,SeD;final=imerode(final,SeD;figure,imshow(final; %显示平滑后的图像(如图9所示title('图9

7、 平滑后的图像'%-!图像中值滤波-%-1 中值滤波-final = medfilt2(final,5 5; %对图像进行中值滤波,滤波窗口默认值为5 5 final = medfilt2(final,5 5;final = medfilt2(final,5 5;final = medfilt2(final,5 5;final = medfilt2(final,5 5;final = medfilt2(final,5 5;final = medfilt2(final,5 5;final = medfilt2(final,5 5;final = medfilt2(final,5 5;fi

8、nal = medfilt2(final,5 5;final = medfilt2(final,5 5;figure;imshow(final; %显示滤波后图像(见图10所示title('图10 滤波后的图像'%-2 标记-L=bwlabel(final,8;s=regionprops(L,'Area','Centroid'k= ismember(L, find(s.Area >= 2200; %去除小于2200的区域number=size(s,1; %孔数area=s.Area; %面积centroid=cat(1,s.Centroid; % Centroid,为以后标记提供位置imshow(k;hold on; %在区域上标记面积(如图11所示%title('图11 区域上标记面积后的图像'for i=1:numbertext(centroid(i,1,centroid(i,2,num2str(area(i; %标记end%-计算平行四边形面积-%根据图形和参照物的像素个数,根据测量原理知道,%S0=4平方厘米,参照物像素个数M0=52458 ,图形像素个数M1=104951 。故可以得到图形面积%S0=4;M0=num2str(area(2;M1=num2str(area(1;disp(M0;