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

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

图1                                               图2

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

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

matlap代码：

clear
clc
I = imread('IMG-0025-00001.jpg');%read pic
p=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中文论坛 )