Matlab医学图像分割区域生长Word文档下载推荐.docx

1、）;figure（1）,imshow（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=65figure（3）,imshow（g2）,title（S=255,T=150figure（4）,imshow（g3）,title（S=150,T=65y,x=ginput（）;x=floor（x）;y=floor（y）;G=zeros（size（g）;G（g3=g3（x,y）=255;结果：代码2：

2、liver1.bmpg=regiongrow（f,255,135）;%肝脏分离figure（2）,imshow（g）,title（S=255,T=135G（g=g（x,y）=255;figure（3）,imshow（G）;2.编写m文件，实现单一型或质心型区域生长算法。由用户通过ginput函数选择种子点，然后根据用户输入的阈值T进行区域生长。输入为：待分割的图片f和阈值T，输出为分割后的二值图片g。（单一型）function a=regiongrow2（f,T）imshow（f）;se=0,1,0;1,1,1;0,1,0;a=im2bw（zeros（size（f）;a（x,y）=1;whil

3、e 1 b=imdilate（a,se）; P=b-a; find（P）; index=find（P）; n=size（index）; for i=1:n t=abs（f（x,y）-f（index（i）; if t=T b（index（i）=0; P（index（i）=0; end a=b; if sum（sum（P）=0 break; end end应用：figure（1）,imshow（f）g=regiongrow3（f,30）;figure（2）,imshow（g）（质心型）function a=regiongrow4（f,T）average=f（x,y）; t=abs（average-

4、double（f（index（i）; indexb=find（b）; N=length（indexb）; All=0;N All=All+double（f（indexb（i）; average=All/N;3.用课本中分水岭的距离变换分水岭分割、梯度分水岭分割和控制标记符分水岭分割对liver1和liver2图像进行分割，分析上述算法的特点和存在的问题，考虑如何进行改进？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=w

5、atershed（-D）;w=L=0;figure（3）,imshow（w）;g2= g & w;figure（5）,imshow（g2）;F=regiongrow（f,255,135）;%F=im2bw（f,graythresh（f）;figure（1）,imshow（F）;h=fspecial（sobelfd=double（f）;g=sqrt（imfilter（fd,h,replicate）.2+. imfilter（fd,h,）.2）;L=watershed（g）;wr=L=0;figure（3）,imshow（wr）;g2=imclose（imopen（g,ones（3,3）,ones

6、（3,3）;L2=watershed（g2）;wr2=L2=0;f2=F;f2（wr2）=255;figure（4）,imshow（f2）;代码3：rm=imregionalmin（g）;im=imextendedmin（f,2）;fim=f;fim（im）=175;Lim=watershed（bwdist（im）;em=Lim=0;figure（4）,imshow（em）;g2=imimposemin（g,im|em）;f2=f;f2（L2=0）=255;figure（5）,imshow（f2）;实验报告：1.实验目的2.实验内容（部分可略写）每题分开写，并标明题号3.实验小结和体会小结：1、第一题重在理解S和T的涵义，若S是一个标量，则它定义一个亮度值，f中有着该值的所有点都变成种子；若T是标量，则它会定义一个全局阈值，阈值用来测试图像中的像素是否与该种子或8连接种子足够相似。2、第二题比较注重算法。用膨胀的方法解决四连接生长点的确定，即先膨胀再进行判断，直到没有出现新的生长点为止。效果与重建很像，但算法不同，适合不同的场合。3、第三题重在理解分水岭的原理，即基于距离变换找出分水岭。但对于梯度分水岭和控制标记符的分水岭分割不理解，需要听老师和其他同学的讲解。还有第二张liver2.bmp有点问题，需要用以前的图片。