[發明專利]一種乳腺圖像中分割胸肌的方法

有效

• 申請號: 201410028053.4
• 申請日: 2014-01-22
• 公開（公告）號: CN104182965A
• 公開（公告）日: 2014-12-03
• 發明（設計）人: 李華;石海林;王瀟崧
• 申請（專利權）人: 上海聯影醫療科技有限公司
• 主分類號: G06T7/00
• 分類號: G06T7/00
• 代理公司: 無
• 代理人: 無

1.一種乳腺圖像中分割胸肌的方法，其特征在于，包括如下步驟：

對乳腺圖像進行形態學限制的區域增長，獲得第一胸肌區域；

對所述第一胸肌區域進行第一校正，獲得分割胸肌區域的乳腺圖像。

2.如權利要求1所述的乳腺圖像中分割胸肌的方法，其特征在于，所述形態學限制的區域增長包括如下步驟：

(1)定義形態學限制；

(2)進行形態學限制的自適應區域增長，所述區域增長的閾值下限為：th_L=(1-t)*(T_max+T_min)／2

其中T_max和T_min分別為當前已生長區域的最大值和最小值；

(3)判斷：若已生長區域的面積S′與形態學限制區域的面積S之差△S≤a*S，則獲得所述第一胸肌區域；

否則，調整參數t＝t*b，并重復所述步驟(2)和(3)；

其中a、b、t均為介于(0，1)之間的常數。

3.如權利要求2所述的乳腺圖像中分割胸肌的方法，其特征在于，所述定義形態學限制包括如下步驟：在所述乳腺圖像上

(1)定義點p(1，y)：取乳腺圖像第一行上胸肌與乳房的分界點為點p；

(2)定義點q(x，1)：取乳腺圖像中乳房的列寬最大值的點所在的行數m，并下移c行，令點q所在行x=m+c；

(3)取所述點p和點q的連線pq，所述連線pq將圖像分割為第一部分和第二部分，所述第一部分為所述形態學限制的區域。

4.如權利要求1所述的乳腺圖像中分割胸肌的方法，其特征在于，所述第一校正包括如下步驟：

(1)對所述第一胸肌區域進行去除突兀結構操作，獲得第二胸肌區域；

(2)對所述第二胸肌區域在n種尺度的乳腺圖像上進行精調，獲得分割胸肌區域的乳腺圖像；

其中，所述n為大于或等于1的正整數。

5.如權利要求4所述的乳腺圖像中分割胸肌的方法，其特征在于，所述去除突兀結構操作包括如下步驟：

遍歷所述第一胸肌區域邊緣上所有的點，取任意兩個端點A和B，連接A、B兩端點形成直線AB，若直線AB同時滿足：

a)直線AB距離

b)直線AB位于所述第一胸肌區域內部；

c)直線AB距離d_AB<d*C；

則去除突兀結構；

若否，則不作處理；

其中，所述L為所述第一胸肌區域中胸肌與乳腺相鄰邊界的總長；所述C為所述第一胸肌區域中連接A、B兩點的胸肌邊緣總長；所述d為大于1的常數。

6.如權利要求4所述的乳腺圖像中分割胸肌的方法，其特征在于，所述精調包括如下步驟：

(1)取所述第二胸肌區域中胸肌的邊界點，取任意一個邊界點，作所述邊界點的法線，取邊界點和邊界點在法線方向上相鄰的至少2個像素點；

(2)計算所述至少3個像素點的灰度值，定義相鄰兩點的灰度值之差為梯度；

(3)確定梯度最大點，定義該點對應的圖像坐標為所述邊界點的更新坐標；

(4)采用Snake算法在n種尺度的乳腺圖像中第二胸肌區域內精調胸肌的邊界，獲得所述分割胸肌區域的乳腺圖像。

其中，所述Snake算法的狀態方程為：

狀態方程E_snake＝E_int+E_ext；

所述更新坐標為所述Snake狀態方程中E_ext的外部能量；所述胸肌邊界的平滑性為所述Snake狀態方程中的內部能量E_int。

7.如權利要求4或6所述的乳腺圖像中分割胸肌的方法，其特征在于，所述精調依據所述乳腺圖像的分辨率由低到高依次進行迭代精調。

8.如權利要求4所述的乳腺圖像中分割胸肌的方法，其特征在于，所述精調前，對所述第二胸肌區域進行邊緣平滑操作。

9.如權利要求1所述的乳腺圖像中分割胸肌的方法，其特征在于，獲得所述第一胸肌區域前，首先所述判斷圖像中乳腺的方向。

10.如權利要求1所述的乳腺圖像中分割胸肌的方法，其特征在于，所述判定乳腺方向的方法包括如下步驟：

取所述乳腺圖像中含有乳腺部分的任意一行1和中間列m，所述第1行和第m列將所述第一圖像的上半部分對稱分割為第三部分和第四部分；

計算所述第三和第四兩部分圖像的灰度平均值T_right和T_left；

若T_right>T_left，則乳腺位于所述第一圖像的右側；

若T_right<T_left，則乳腺位于所述第一圖像的左側。