资源描述:
《图像处理中的标记分水岭分割算法》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、.图像处理中的标记分水岭分割算法如果图像中的目标物体是连接在一起的,则分割起来会更困难,分水岭分割算法经常用于处理这类问题,通常会取得比较好的效果。分水岭分割算法把图像看成一幅“地形图”,其中亮度比较强的区域像素值较大,而比较暗的区域像素值较小,通过寻找“汇水盆地”和“分水岭界限”,对图像进行分割。直接应用分水岭分割算法的效果往往并不好,如果在图像中对前景对象和背景对象进行标注区别,再应用分水岭算法会取得较好的分割效果。有很多图像处理工具箱函数可以用到,如fspecial、imfilter、watershed、lable2rgb、imopen
2、、imclose、imreconstruct、imcomplement、imregionalmax、bwareaopen、graythresh、和imimposemin函数等。下面进行一个例子,步骤如下。1、读取图像并求其边界,代码如下。rgb=imread('pears.png');%读取原图像I=rgb2gray(rgb);%转化为灰度图像figure;subplot(121)%显示灰度图像imshow(I)text(732,501,'ImagecourtesyofCorel',...'FontSize',7,'HorizontalAli
3、gnment','right')hy=fspecial('sobel');%sobel算子hx=hy';Iy=imfilter(double(I),hy,'replicate');%滤波求Y方向边缘Ix=imfilter(double(I),hx,'replicate');%滤波求X方向边缘gradmag=sqrt(Ix.^2+Iy.^2);%求模subplot(122);imshow(gradmag,[]),%显示梯度title('Gradientmagnitude(gradmag)')在这一步骤中,首先读取一套真彩色图像,然后把真色图像转
4、化为灰度图像,结果如图所示:..图1原图和梯度图像使用sobel边缘算子对图像进行水平和垂直方向的滤波,然后求取模值,sobel算子滤波后的图像在边缘处会显示比较大的值,在没有边界处的值会很小,如上图右图所示。1、直接使用梯度模值进行分水岭算法,代码如下。L=watershed(gradmag);%直接应用分水岭算法Lrgb=label2rgb(L);%转化为彩色图像figure;imshow(Lrgb),%显示分割后的图像title('Watershedtransformofgradientmagnitude(Lrgb)')直接使用梯度模值
5、图像进行分水岭算法得到的结果往往会存在过度分割的现象,如下图所示。因此通常需要对前景对象和背景对象进行标记,以获得更好的分割效果。..图2直接对梯度图像进行分水岭分割1、分别对前景和背景进行标记,代码如下。se=strel('disk',20);%圆形结构元素Io=imopen(I,se);%形态学开操作figure;subplot(121)imshow(Io),%显示执行开操作后的图像title('Opening(Io)')Ie=imerode(I,se);%对图像进行腐蚀Iobr=imreconstruct(Ie,I);%形态学重建sub
6、plot(122);imshow(Iobr),%显示重建后的图像title('Opening-by-reconstruction(Iobr)')Ioc=imclose(Io,se);%形态学关操作figure;subplot(121)imshow(Ioc),%显示关操作后的图像title('Opening-closing(Ioc)')Iobrd=imdilate(Iobr,se);%对图像进行膨胀Iobrcbr=imreconstruct(imcomplement(Iobrd),...imcomplement(Iobr));%形态学重建Iob
7、rcbr=imcomplement(Iobrcbr);%图像求反subplot(122);imshow(Iobrcbr),%显示重建求反后的图像title('Opening-closingbyreconstruction(Iobrcbr)')fgm=imregionalmax(Iobrcbr);%局部极大值figure;imshow(fgm),%显示重建后局部极大值图像title('Regionalmaximaofopening-closingbyreconstruction(fgm)')..I2=I;I2(fgm)=255;%局部极大值处像
8、素设置为255figure;imshow(I2),%在原图上显示极大值区域title('Regionalmaximasuperimposedonoriginalim