matlap实现肺实质区域初始分割,去除肺部气管及背景

来源:互联网 发布:keep设置网络 编辑:程序博客网 时间:2024/03/29 21:03


下图(图1)是通过直线扫描方法选取躯干种子点进行区域生长后的效果。为了提取肺实质区域(标记为1和2),我们先将图像进行反转,再进行8连通区域标记,以去除中间肺气管(标记为3)。

图1                                               图2

   但实验发现,无法去除。原因可由图像得知,气管和肺实质在连通标记后粘连在一块。  

为此,我们在进行连通标记前先进行腐蚀操作,扩大肺实质和气管间距离。再进行连通标记处理,去除面积较小的气管和面积较大的背景。

matlap代码:

clearclcI = imread('IMG-0025-00001.jpg');%read picp=imcomplement(I);se=strel('square',2);se1=strel('disk',2);bw=imerode(p,se);[L,N]=bwlabel(bw,4);s = regionprops(L,'Area');bw1=ismember(L,find([s.Area]>=10000 & [s.Area]<=50000  ));bw2=imdilate(bw1,se);figure;subplot(2,3,1);imshow(p);title('binaryinvert');subplot(2,3,2);imshow(bw);title('beforeimeroderegionpops');subplot(2,3,3);imshow(bw1);title('afterimerodese>10000<50000');subplot(2,3,4);imshow(bw2);title('imdilatese');
显示效果如下:


后续我们还需要对肺实质轮廓进行修补操作以及肺结节的分割,假阳性去除等操作

本文相关链接可见我的matlap论坛帖子:区域连通标记去除小面积区域
http://www.ilovematlab.cn/thread-326016-1-1.html
(出处: MATLAB中文论坛 )


0 0