图像处理实验三--圆形物体的图像分割与分析

实验三圆形物体的图像分割与分析一、实验目的掌握在计算机上进行图像分割与形态学图像处理和分析的方法二、实验要求应用所学方法完成coins.tif图像的分割，用图像分析的方法估算整幅图像中钱币的总额。1)对指定图像（coins.tif）完成图像分割、将分割后的图像显示于屏幕上。2)对分割后图像中不同目标的尺寸（面积或半径）分布进行直方图显示。3)讨论你的方法的误差和局限。4)图像中有2元、1元、5角和1角的四种硬币，计算图像中钱币的总额。三、实验仪器设备及软件HPD538、MATLAB四、实验原理以自己的语言结合课堂笔记进行总结，要求过程推导清晰明了。五、实验步骤及程序实验步骤、程序流程、实验源程序和注释齐全实验源程序: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);gridA=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