[發(fā)明專利]基于胸部CT圖像的肺部自動提取方法

有效

• 申請?zhí)枺?/td>: 201310745336.6
• 申請日: 2013-12-30
• 公開（公告）號: CN104751437A
• 公開（公告）日: 2015-07-01
• 發(fā)明（設(shè)計）人: 黃雪艷
• 申請（專利權(quán)）人: 藍網(wǎng)科技有限公司
• 主分類號: G06T7/00
• 分類號: G06T7/00
• 代理公司: 深圳冠華專利事務(wù)所(普通合伙) 44267
• 代理人: 諸蘭芬

1.一種基于胸部CT圖像的肺部自動提取方法，其特征在于，所述方法包括步驟：

A、將胸部CT圖像進行二值化，得到初始二值圖像；

B、求出初始二值圖像中的最大連通域，在初始二值圖像中去除最大連通域以外的其他區(qū)域；

C、利用歐氏距離變換提取帶狀區(qū)域，去除與胸腔壁粘連的干擾區(qū)；

D、求得初提取出的肺部二值圖，將初提取出的肺部二值圖像進行反色，得到初提取出的肺部二值反色圖，對初提取出的肺部二值圖像反色圖進行歐氏距離變換求得歐氏距離變換圖；

E、利用初提取出的肺部二值反色圖和歐氏距離變換圖，得到肺部缺失的區(qū)域二值圖；

F、將肺部缺失的區(qū)域二值圖與初提取出的肺部二值圖相或后，判斷左右肺是否粘連，若粘連，則將左右肺進行分離，得到完整的肺部區(qū)域二值圖模版。

2.根據(jù)權(quán)利要求1所述基于胸部CT圖像的肺部自動提取方法，其特征在于，步驟B中找出最大連通域的步驟包括：

對初始二值圖像進行八領(lǐng)域標記；

計算出每個標記區(qū)域的面積，最后找到面積最大的區(qū)域就是最大連通域。

3.根據(jù)權(quán)利要求2所述基于胸部CT圖像的肺部自動提取方法，其特征在于，對二值化圖像進行八領(lǐng)域標記的步驟包括：

從左向右，從上到下掃描初始二值圖像，對未標記的目標像素點進行標記，對該像素點的八領(lǐng)域點中沒有出界且不為背景點和已標記點的像素點用相同的標記號進行標記，對該標號的點數(shù)進行統(tǒng)計；

利用相同的方法對新標記的像素點的八領(lǐng)域進行相同的搜索與標記，直到該標記號的像素點的八領(lǐng)域點都為已標記點或為背景點或出界，則該連通區(qū)域標記完成，將標記號加一，對下個連通區(qū)域進行標記。

4.根據(jù)權(quán)利要求1所述基于胸部CT圖像的肺部自動提取方法，其特征在于，步驟C和步驟D中利用歐氏距離變換的步驟具體包括：

求取初提取出的肺部二值圖的邊界二值圖像和去除邊界后的二值圖像；

生成初始化歐氏距離變換圖，將得到的邊界二值圖像中為1的像素點定義為初始化歐氏距離變換圖的(0,0)，初始化歐氏距離變換圖中其他像素點填充為極大值(MaxDistance,MaxDistance)，

MaxDistance=max(1.414*ImageWidth,1.414*ImageHeight)，式中，ImageWidth和ImageHeight分別為胸部CT圖像圖像的寬度和高度；

使用歐氏距離變換模板逐行掃描初始化歐氏距離變換圖，將歐氏距離變換模板和歐氏距離變換模板覆蓋下的像素點對應(yīng)相加得（x_i,y_i），將計算(x_i*x_i+y_i*y_i)得到最小值所對應(yīng)的點（x_i,y_i）賦值給中心像素（pEDMap.x，pEDMap.y），然后將初始化歐氏距離變換圖中各點像素取值為sqrt(pEDMap.x*pEDMap.x+pEDMap.y*pEDMap.y)，得到歐氏距離變換圖；

其中，pEDMap.x和pEDMap.y分別表示初提取出的肺部二值圖的目標像素點到距離其最近的背景像素點，在x方向和y方向的距離。

5.根據(jù)權(quán)利要求1所述基于胸部CT圖像的肺部自動提取方法，其特征在于，步驟C是將經(jīng)過步驟B處理后的初始二值圖像與提取的帶狀區(qū)域二值圖相減，將與胸腔壁粘連的干擾區(qū)去除。

6.根據(jù)權(quán)利要求1所述基于胸部CT圖像的肺部自動提取方法，其特征在于，判斷左右肺是否粘連的步驟包括：對肺部缺失的區(qū)域二值圖與初提取出的肺部二值圖相或后的圖像進行連通區(qū)域個數(shù)的判斷，若只有一個連通區(qū)域，且連通區(qū)域的像素個數(shù)大于預(yù)設(shè)閾值，則判斷左右肺粘連。