实验三

实验三圆形物体的图像分割与分析一、实验目的掌握在计算机上进行图像分割与形态学图像处理和分析的方法二、实验要求应用所学方法完成coins.tif图像的分割，用图像分析的方法估算整幅图像中钱币的总额。1)对指定图像（coins.tif）完成图像分割、将分割后的图像显示于屏幕上。2)对分割后图像中不同目标的尺寸（面积或半径）分布进行直方图显示。3)讨论你的方法的误差和局限。4)图像中有2元、1元、5角和1角的四种硬币，计算图像中钱币的总额。三、实验仪器设备及软件HPD538、MATLAB四、实验原理1）图像腐蚀：设A为目标图像，B为结构元素，则目标图像A被结构元素B腐蚀可定义为：其中，y是一个表示集合{ABx|(B)A平移的位移量。2）图像膨胀：设A为目标图像，B为结构元素，则目标图像A被结构元素B膨胀可定义为：其中，y是一个表示集合{AB{x|((B)A)}平移的位移量。五、实验步骤及程序I1=imread('LENA.BMP');I=double(I1);%转化为double型clearall;I=imread('coins.bmp');subplot(231);imshow(I);title('原图');I1=im2bw(I,graythresh(I));subplot(232);imshow(I1),title('原图二值化');H=medfilt2(I1);subplot(233);imshow(H),title('原图中值滤波');SE=strel('disk',4);%结构元素，创建半径为4的圆盘G=imdilate(H,SE);%膨胀D=bwdist(G);%计算二值图象欧式距离clearall;I=imread('coins.bmp');subplot(231);imshow(I);title('原图');I1=im2bw(I,graythresh(I));subplot(232);imshow(I1),title('原图二值化');H=medfilt2(I1);subplot(233);imshow(H),title('原图中值滤波');SE=strel('disk',4);%结构元素，创建半径为4的圆盘G=imdilate(H,SE);%膨胀D=bwdist(G);%计算二值图象欧式距离L=watershed(-D);%分水岭分割subplot(234);imshow(L),title('分水岭分割');R=(L==0);%得到目标区域Q=imerode(G,SE);%腐蚀S=((~R)&(~Q));subplot(235);imshow(S),title('分割后图像');[Ilabel,num]=bwlabel(S,8);%对图象标T=regionprops(Ilabel,'EquivDiameter');%测量图像的区域特性U=[T.EquivDiameter];subplot(236);hist(U,6);grid;title('直方图');A=length(find(U35));B=length(find(U35&U30));C=length(find(U30&U25.2));D=length(find(U25.2&U20));SUM=A*2+B*1+C*0.5+D*0.1六、结果与分析SUM=38.1000分析：通过对图像二值化以及分割处理，可以较为准确地实现硬币的分类与统计。