[发明专利]乳腺X光图像中胸肌区域的分割方法

权利要求书

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的胸肌壁分界线。