[發明專利]乳腺X光圖像中胸肌區域的分割方法

在審

• 申請號: 201210367387.5
• 申請日: 2012-09-28
• 公開（公告）號: CN103700085A
• 公開（公告）日: 2014-04-02
• 發明（設計）人: 趙明;康雨
• 申請（專利權）人: 深圳市藍韻實業有限公司
• 主分類號: G06T7/00
• 分類號: G06T7/00
• 代理公司: 深圳冠華專利事務所(普通合伙) 44267
• 代理人: 諸蘭芬

1.一種乳腺X光圖像中胸肌區域的分割方法，其特征在于，包括步驟：?

B、在乳腺X光圖像中選取寬度是乳腺X光圖像寬度的1/N、高度是乳腺X光圖像高度的1/N的矩形區域定義為感興趣區域ROI，且感興趣區域ROI中至少包含了部分的胸肌區域；?

C、將興趣區域ROI二值化分割獲得二值化圖像ROI_BW；?

D、提取二值化圖像ROI_BW的邊界圖像ROI_edge；?

E、檢測邊界圖像ROI_edge中的直線，確定感興趣區域ROI的胸肌壁分界線；?

F、將胸肌壁分界線延伸至與邊界圖像ROI_edge的邊界相交，獲得乳腺X光圖像中腺體組織與胸肌壁的分界線；?

G、依據獲得的腺體組織與胸肌壁的分界線對乳腺X光圖像進行分割處理。?

2.根據權利要求1所述乳腺X光圖像中胸肌區域的分割方法，其特征在于，在步驟B之前還包括步驟：?

A、對乳腺X光圖像進行預處理，去除乳腺X光圖像的背景和標簽。?

3.根據權利要求2所述乳腺X光圖像中胸肌區域的分割方法，其特征在于，步驟A具體包括：?

計算乳腺X光圖像中所有像素點的總平均灰度u：u＝w0*u0+w1*u1；?

從最小灰度值到最大灰度值遍歷t，按公式g＝w0*(u0-u)²+w1*(u1-u)²計算每一個t值的g值，當g最大時的t值即為全局閾值T；?

將像素點的值大于全局閾值T的圖像區域全部置0，小于全局閾值T的像素點不做改變；?

其中，所有組織的像素點數在乳腺X光圖像中占總像素點數的比例為w0，所有組織的像素點的平均灰度為u0，背景的像素點數在乳腺X光圖像中占總像素點數的比例為w1，背景的像素點數的平均灰度為u1。?

4.根據權利要求1所述乳腺X光圖像中胸肌區域的分割方法，其特征在于，M、N都為3。?

5.根據權利要求1所述乳腺X光圖像中胸肌區域的分割方法，其特征在于，步驟C具體包括：?

C1、令i＝0，1，…，GL_max為感興趣區域ROI中的每一個灰度級，統計每個灰度級別的在感興趣區域ROI中所占的比例p(i)，給定初始閾值

C2、計算感興趣區域ROI中高于閾值部分統計值均值和低于閾值部分的統計值均值

C3、根據m_b和m_a的值更新t_ROI，

C4、重復以上步驟C1-C3直到|t_ROI-m_b|＜1，這時得到的t_ROI為感興趣區域ROI的迭代閾值；?

C5、根據這個迭代閾值將感興趣區域ROI進行二值化：大于迭代閾值的像素為1，小于迭代閾值的像素點為0，得到的二值化圖像ROI_BW。

6.根據權利要求1所述乳腺X光圖像中胸肌區域的分割方法，其特征在于，步驟D包括：?

使用高斯平滑模板對二值化圖像ROI_BW進行平滑處理；?

使用Canny邊緣算子提取平滑結果中所有的邊緣，獲得邊界圖像ROI_edge。?

7.根據權利要求1所述乳腺X光圖像中胸肌區域的分割方法，其特征在于，步驟E包括：?

使用霍夫變換在邊界圖像ROI_edge中檢測直線；?

如果檢測結果中包含多條直線，則選擇斜率在1～3.7之間的且斜率最小的?一條直線作為感興趣區域ROI內的胸肌壁分界線；?

如果霍夫變換僅僅檢測出一條直線，那么就將這條直線作為感興趣區域ROI的胸肌壁分界線。?