实验4 医学图像分割(二).pdfVIP

实验4医学图像分割（二） 实验目的： 1. 掌握形态学图像处理的原理、函数； 2. 会灵活选取模板形态学运算的模板； 3. 理解图像分割的基本理论和方法； 4. 掌握阈值分割的方法和阈值的选择； 实验内容： 1. 对图片text.tif进行膨胀操作，选择结构元素为[010;111;0 1 0] [1 1 1;1 1 1;1 1 1] 和 观察其处理效果的异同，如果一次 膨胀的效果不明显，可多次膨胀再比较。 原图 f=imread(text.tif); 结构元素为[010;111;010] w=[010;111;010]; f1=imdilate(f,w); imshow(f);figure,imshow(f1); f2=imdilate(f1,w); imshow(f2); 二次膨胀 结构元素为[111;111;111] 2. 对图片shapes.tif进行开、闭运算，比较其处理的效果。 原图 开运算 闭运算 f=imread(shapes.tif); se=strel(square,20); fo=imopen(f,se); imshow(f);figure,imshow(fo); fc=imclose(f,se); imshow(fc); 3. 用 形 态 学 方 法 ， 对 存 在 杂 散 点 的 指 纹 图 片 noisy-fingerprint.tif进行处理。 f=imread(noisy-fingerprint.tif); se=strel(square,3); fo=imopen(f,se); foc=imclose(fo,se); imshow(f),figure,imshow(fo),figure,imshow(foc); 对图像做开运算来消除杂散点，再用闭运算补充开运算产生的缺口 1 lung.bmp im2bw 4. （ ）对图片 用 进行处理，其分割的阈值 0.05,0.15 0.7 分别设为 和 ，分析哪个效果好（结合图像的 直方图）。 f=imread(lung.bmp); f1=im2bw(f); imhist(f); 2 Otsu （ ）用全局阈值法或 算法计算出合适的阈值。 T=graythresh(f); T T= 0.1647 g1=im2bw(f,0.05); g2=im2bw(f,0.15); g3=im2bw(f,0.7); figure,imshow(g1) figure,imshow(g2) figure,imshow(g3) figure,imhist(f) k=graythresh(f); g=im2bw(f,k); figure,imshow(g) 3 （ ）提高题：对分割出来的二值图像运用所学的形态学知 识，得到肺部图像。 提示： 如 为原图像，f fillhole为处理后的二值图像（其中肺部为 白色，其他为黑色），则以下语句可以得到原图中的肺部图 像。 f(fillhole==0)=255;% f(fillhole==0)=0 或 也可以 f=imread(lung.bmp); k=graythresh(f); g=im2bw(f,k); imshow(g); h=1-g; a=imfill(h,holes); f(a==0)=255; imshow(f); 原图 Fillhole 最终图