資源描述:
《Matlab醫(yī)學(xué)圖像分割區(qū)域生長》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫。
1、實(shí)驗(yàn)十醫(yī)學(xué)圖像分割(二)實(shí)驗(yàn)?zāi)康模?.了解圖像分割的基本理論和方法;2.掌握閾值分割的方法和閾值的選擇;3.掌握基于分水嶺分割的原理和應(yīng)用;實(shí)驗(yàn)內(nèi)容:1.區(qū)域生長法利用圖像像素間的相似性進(jìn)行分割.調(diào)用regiongrow函數(shù)對(duì)圖像weld.tif進(jìn)行處理.注意參數(shù)中S(種子值).T(閾值)的選擇對(duì)分割效果的影響。S=255.T=65和S=255T=150和S=150.T=65三組值進(jìn)行處理.理解在區(qū)域生長法的原理。同時(shí)對(duì)liver.bmp.自己選擇合適的S和T.以較好得分割出肝臟。代碼1:closeall;f=imread('weld.tif');figure(1),im
2、show(f),title('原圖');[g1,NR]=regiongrow(f,255,65);[g2,NR]=regiongrow(f,255,150);[g3,NR]=regiongrow(f,150,65);figure(2),imshow(g1),title('S=255,T=65');figure(3),imshow(g2),title('S=255,T=150');figure(4),imshow(g3),title('S=150,T=65');..[y,x]=ginput();x=floor(x);y=floor(y);G=zeros(size(g));G
3、(g3==g3(x,y))=255;結(jié)果:....代碼2:closeall;f=imread('liver1.bmp');figure(1),imshow(f),title('原圖');g=regiongrow(f,255,135);%肝臟分離figure(2),imshow(g),title('S=255,T=135');[y,x]=ginput();x=floor(x);y=floor(y);..G=zeros(size(g));G(g==g(x,y))=255;figure(3),imshow(G);結(jié)果:....1.編寫m文件.實(shí)現(xiàn)單一型或質(zhì)心型區(qū)域生長算法。由用
4、戶通過ginput函數(shù)選擇種子點(diǎn).然后根據(jù)用戶輸入的閾值T進(jìn)行區(qū)域生長。輸入為:待分割的圖片f和閾值T.輸出為分割后的二值圖片g。代碼1:(單一型)function[a]=regiongrow2(f,T)imshow(f);se=[0,1,0;1,1,1;0,1,0];[y,x]=ginput();x=floor(x);y=floor(y);a=im2bw(zeros(size(f)));a(x,y)=1;while1b=imdilate(a,se);P=b-a;find(P);index=find(P);n=size(index);fori=1:nt=abs(f(x,y
5、)-f(index(i)));..ift>=Tb(index(i))=0;P(index(i))=0;endenda=b;ifsum(sum(P))==0break;endend應(yīng)用:f=imread('weld.tif');figure(1),imshow(f)g=regiongrow3(f,30);figure(2),imshow(g)代碼2:(質(zhì)心型)function[a]=regiongrow4(f,T)imshow(f);se=[0,1,0;1,1,1;0,1,0];[y,x]=ginput();x=floor(x);y=floor(y);..a=im2bw(z
6、eros(size(f)));a(x,y)=1;average=f(x,y);while1b=imdilate(a,se);P=b-a;index=find(P);n=size(index);fori=1:nt=abs(average-double(f(index(i))));ift>=Tb(index(i))=0;P(index(i))=0;endenda=b;indexb=find(b);N=length(indexb);All=0;fori=1:NAll=All+double(f(indexb(i)));end..average=All/N;ifsum(sum(P)
7、)==0break;endend結(jié)果:..1.用課本中分水嶺的距離變換分水嶺分割、梯度分水嶺分割和控制標(biāo)記符分水嶺分割對(duì)liver1和liver2圖像進(jìn)行分割.分析上述算法的特點(diǎn)和存在的問題.考慮如何進(jìn)行改進(jìn)?代碼1:closeall;f=imread('liver1.bmp');figure(1),imshow(f);..%g=im2bw(f,graythresh(f));[g,NR]=regiongrow(f,255,135);%肝臟分離figure(2),imshow(g);gc=~g;D=bwdist(gc);L=wa