[發明專利]肺區CT序列自動分割方法

無效

• 申請號: 201310257482.4
• 申請日: 2013-06-26
• 公開（公告）號: CN103400365A
• 公開（公告）日: 2013-11-20
• 發明（設計）人: 曲建明
• 申請（專利權）人: 成都金盤電子科大多媒體技術有限公司
• 主分類號: G06T7/00
• 分類號: G06T7/00;G06K9/34;A61B6/03
• 代理公司: 成都中亞專利代理有限公司 51126
• 代理人: 王崗

1.一種肺區CT序列自動分割方法，其特征在于：包括如下步驟；

（1）輸入肺區CT序列圖像；?

（2）利用交互式區域生長對圖像I進行分割；

（3）據分割結果得到初始輪廓，計算相鄰影像的種子點坐標；

（4）基于I中種子點坐標利用步驟2中交互式區域生長對序列中的當前圖像II進行分割；

（5）重復2、3、4步驟，判斷CT序列中的圖像是否已全部分割完成，否則轉步驟（3），直到所有圖像分割完成。

2.根據權利要求1所述的肺區CT序列自動分割方法，其特征在于：

步驟（2）按照如下方式進行；

利用OTSU算法求得CT圖像最優閾值作為終止條件并利用四角置點區域生長算法得到背景區域，進而連通域標記合并小面積區域去除背景，得到Body-Mask；

其中采用四角置點區域生長的方法濾除背景流程如下：

a、閾值以上文得到的最優閾值作為迭代終止條件；

b、分別從影像的四角向中心區域掃描，若掃描的像素值小于閾值即認為是背景區域并將像素值置為0（黑色），不斷掃描直至所有背景像素點合并完畢；

c、生長過程：從種子點開始遍歷其鄰域像素，若像素灰度值小于閾值，則標記為背景，若大于閾值且其四鄰域有小于閾值的像素，則標記為邊界；

d、將標記為背景的像素置為白色；

標記并針對軀干模板數據，以步驟中閾值分割并取反，得到粗略肺實質模板；?

對二值的肺實質模板輪廓跟蹤，對得到的閉合輪廓對象進行面積過濾，去除主氣管，心臟血管等孤立小面積閉合區域；

構造圓形結構數據對肺實質模板邊界進行平滑修復。

3.根據權利要求1所述的肺區CT序列自動分割方法，其特征在于：步驟（3）具體流程如下：

?遍歷當前層影像的目標區域邊界，得到灰度均值A及灰度方差S；

據邊界點坐標由逼近算法迭代得到目標區域正放置最大內接矩形R，并得到形心坐標（I，J）；

將（I，J）投影到下一層影像的相應位置，并提取其灰度值Gray；

計算abs（Gray-A），若abs（Gray-A）小于S，則種子點有效；否則，種子點無效，則以該點為圓心，在半徑不大于R的范圍內遍歷其鄰域像素，重復（4）操作，從而最終得到種子點。

4.根據權利要求3所述肺區CT序列自動分割方法，其特征在于：其中用于種子點計算的最大內接矩形近似逼近算法：?

在得到的目標邊界上尋找X，Y方向的中軸線，并劃分邊界坐標為四段子邊界；

選定四段中任一段邊界，選定為矩形頂點（X，Y），構建對角方向射線，得到對角射線方向與邊界點的交點，構建矩形，判斷該矩形邊界是否與邊界點有交點，有則舍棄，無則記錄矩形面積Si，頂點坐標，方向角到數組；

以一定的角度旋轉對角方向射線，重復b步驟；

以M為步長，在該字段邊界選定頂點重復迭代a，b過程；若子段遍歷結束未找到則遍歷下一子段；

子段邊界遍歷結束，在面積數組中選取面積最大的數組項，作為近似的最大內接矩形；求的形心坐標（I，J），同時可得最大內切圓半徑R。