14组实验六 实验6 图像分割要点.doc

实验6 图像分割 实验目的： 熟悉区域生长法； 2. 分水岭分割算法 实验内容： 区域生长法利用图像像素间的相似性进行分割，调用regiongrow函数对图像weld.tif进行处理，注意参数中S（种子值），T（阈值）的选择对分割效果的影响。S=255，T＝65和S＝255 T＝150和S＝150，T＝65三组值进行处理，理解在区域生长法的原理。同时对liver.bmp，自己选择合适的S和T，以较好的分割出肝脏。参考书本P309例10.8 答： f=imread(weld.tif); imshow(f) title(weld.tif); f1=regiongrow(f,255,65); figure,imshow(f1) title(S=255,T=65); f2=regiongrow(f,255,150); figure,imshow(f2) title(S=255,T=150); f3=regiongrow(f,150,65); figure,imshow(f3) title(S=150,T=65); f=imread(liver.bmp); imshow(f) title(liver.bmp); [g, NR, SI, TI] = regiongrow(f,175,20); figure,imshow(g) title(区域生长 S=175，T=20); gf=imfill(g); figure,imshow(gf) title(孔洞填充); se=strel(disk,10); ge=imerode(gf,se); figure,imshow(ge) title(腐蚀); ir=imreconstruct(ge,gf); figure,imshow(ir) title(重构); se=strel(disk,2); d=imdilate(ir,se); figure,imshow(d); title(膨胀 分离出的肝脏); 利用分水岭分割算法对图像rice.tif进行分割，并对出现过分割现象提出解决方案。(可以尝试tophat等学过的算法对图像进行预处理后再分割；使用不同的结构元素和参数；使用梯度分水岭算法与标记符控制算法等。比较得到最好的结果）。 答： f=imread(rice.tif); imshow(f) title(rice.tif); g=im2bw(f,graythresh(f)); gc=~g; D=bwdist(gc); L=watershed(-D); w=L==0; g2= g ~w; figure,imshow(g2) title(用距离和分水岭变换分割) f=imread(rice.tif); se=strel(disk,10); ft=imtophat(f,se); h=fspecial(sobel); fd=double(ft); g=sqrt(imfilter(fd,h,replicate) .^2+imfilter(fd,h,replicate) .^2); L=watershed(g); wr =L==0; g2=imclose(imopen(g,ones(3,3)),ones(3,3)); L2 = watershed(g2); wr2 = L2==0; f2=ft; f2(wr2)=255; figure,imshow(f2) title(用梯度的分水岭分割) f=imread(rice.tif); se=strel(disk,10); ft=imtophat(f,se); h=fspecial(sobel); fd=double(ft); g=sqrt(imfilter(fd,h,replicate) .^2+imfilter(fd,h,replicate) .^2); L=watershed(g); t=graythresh(ft); im = imextendedmin(ft,t*255); im=~im; Lim=watershed(bwdist(im)); em=Lim==0; g2=imimposemin(g,im|em); L2=watershed(g2); f2=ft; f2(L2==0)=255; figure,imshow(f2) title(控制标记符的分水岭分割) f=imread(rice.tif); se=strel(disk,20); f=imtophat(fse); hy = fspecial(sobel); hx = hy; fy = imfilter(double(f), hy, replicate); fx = i