乳腺x光图像中胸肌区域的分割方法
【专利摘要】本发明公开一种乳腺X光图像中胸肌区域的分割方法，其包括：在乳腺X光图像中确定感兴趣区域ROI；将兴趣区域ROI二值化分割获得二值化图像ROIBW；提取二值化图像ROIBW的边界图像ROIedge；检测边界图像ROIedge中的直线，确定感兴趣区域ROI的胸肌壁分界线；将胸肌壁分界线延伸至与边界图像ROIedge的边界相交，获得乳腺X光图像中腺体组织与胸肌壁的分界线；依据获得的腺体组织与胸肌壁的分界线对乳腺X光图像进行分割处理。本发明利用乳腺X光图像中腺体组织与胸肌壁之间灰度特征的差异，确定出腺体组织与胸肌壁的分界线，能够准确定位胸肌区域进行图像分割。
【专利说明】乳腺X光图像中胸肌区域的分割方法
【技术领域】
[0001]本发明涉及一种医学图像处理技术，尤其是涉及一种乳腺X光图像中胸肌区域的分割方法。
【背景技术】
[0002]乳腺X光(mammography)图像原片清晰度和对比度的经常不足，为了帮助临床医生更好地诊断病情，需要对原片进行图像处理和增强。而在侧斜位(MLO，MediolateralOblique)图像中，乳腺组织往往含有胸肌部分，其性质与乳腺区域有较大的不同且不含钙化点。根据每个区域各自独特的图像特征进行信息检测与处理，是微钙化和肿块检测研究的前提和重要基础。胸肌区域的分割在乳腺X线图像的左右乳比较研究、肿块检测区域定界方面有着重要的作用。目前能够搜索到得近似相关专利有:
[0003]中国申请号200910100655.5的专利申请，公开一种用于乳腺钥靶X线图像中乳腺区域提取的预处理方法，此方法简单地使用中值滤波方法和图像中左右两侧像素密度分布来判定乳腺区域在图像中的位置，然后使用水平和垂直未曝光边框剥离得到最终的乳房区域。此方法无法区分出胸肌区域。
[0004]中国申请号200810217937.9的专利申请，一种确定乳腺外边缘轮廓的方法，此方法使用图像质心提取、图像二值化及区域生长等技术，确定乳腺外边缘轮廓。但是其方法无法定位胸肌区域，并且区域生长法产生的边缘容易丢失皮肤线。
[0005]现有技术一般采用简单的阈值分割和密度分布进行乳腺X光图像分割，只能得到单一的乳房轮廓区域，对于轮廓区域中含有大面积胸肌的MLO图像则无法进一步区分开来。这种分割方法给乳腺X光图像的后续处理和计算机辅助诊断带来了困难。
【发明内容】
[0006]本发明提出一种乳腺X光图像中胸肌区域的分割方法，以解决目前分割方法存在无法准确定位胸肌区域导致难于分割出胸肌区域的技术问题。
[0007]本发明采用如下技术方案实现:一种乳腺X光图像中胸肌区域的分割方法，其包括步骤:
[0008]B、在乳腺X光图像中选取宽度是乳腺X光图像宽度的1/N、高度是乳腺X光图像高度的1/N的矩形区域定义为感兴趣区域R0I，且感兴趣区域ROI中至少包含了部分的胸肌区域；
[0009]C、将兴趣区域ROI 二值化分割获得二值化图像ROIbw ；
[0010]D、提取二值化图像ROIbw的边界图像ROIedge ；
[0011]E、检测边界图像ROIedge中的直线，确定感兴趣区域ROI的胸肌壁分界线；
[0012]F、将胸肌壁分界线延伸至与边界图像ROIedge的边界相交，获得乳腺X光图像中腺体组织与胸肌壁的分界线；
[0013]G、依据获得的腺体组织与胸肌壁的分界线对乳腺X光图像进行分割处理。[0014]其中，在步骤B之前还包括步骤:A、对乳腺X光图像进行预处理，去除乳腺X光图像的背景和标签。
[0015]其中，步骤A具体包括:计算乳腺X光图像中所有像素点的总平均灰度u:u =w0*u0+wl*ul ;从最小灰度值到最大灰度值遍历t,按公式g = wO*(uO-u)2+wl*(ul-u)2计算每一个t值的g值，当g最大时的t值即为全局阈值T ;将像素点的值大于全局阈值T的图像区域全部置0，小于全局阈值T的像素点不做改变；其中，所有组织的像素点数在乳腺X光图像中占总像素点数的比例为w0，所有组织的像素点的平均灰度为uO，背景的像素点数在乳腺X光图像中占总像素点数的比例为wl，背景的像素点数的平均灰度为ul。
[0016]其中，Μ、N都为3。
[0017]其中，步骤C具体包括:Cl、令i = 0，1，…，GLmax为感兴趣区域ROI中的每一个灰度级，统计每个灰度级别的在感兴趣区域ROI中所占的比例p(i)，给定初始阈值
【权利要求】
1.一种乳腺X光图像中胸肌区域的分割方法，其特征在于，包括步骤: B、在乳腺X光图像中选取宽度是乳腺X光图像宽度的1/N、高度是乳腺X光图像高度的1/N的矩形区域定义为感兴趣区域ROI，且感兴趣区域ROI中至少包含了部分的胸肌区域； C、将兴趣区域ROI二值化分割获得二值化图像ROIbw ； D、提取二值化图像ROIbw的边界图像ROIedge； E、检测边界图像ROIedge中的直线，确定感兴趣区域ROI的胸肌壁分界线； F、将胸肌壁分界线延伸至与边界图像ROIedge的边界相交，获得乳腺X光图像中腺体组织与胸肌壁的分界线； G、依据获得的腺体组织与胸肌壁的分界线对乳腺X光图像进行分割处理。
2.根据权利要求1所述乳腺X光图像中胸肌区域的分割方法，其特征在于，在步骤B之前还包括步骤: A、对乳腺X光图像进行预处理，去除乳腺X光图像的背景和标签。
3.根据权利要求2所述乳腺X光图像中胸肌区域的分割方法，其特征在于，步骤A具体包括: 计算乳腺X光图像中所有像素点的总平均灰度u:u = w0*u0+wl*ul ； 从最小灰度值到最大灰度值遍历t,按公式g = w0*(u0-u)2+wl* (ul-u)2计算每一个t值的g值，当g最大时的t值即为全局阈值T ； 将像素点的值大于全局阈值T的图像区域全部置0，小于全局阈值T的像素点不做改变； 其中，所有组织的像素点数在乳腺X光图像中占总像素点数的比例为W0，所有组织的像素点的平均灰度为U0，背景的像素点数在乳腺X光图像中占总像素点数的比例为wl，背景的像素点数的平均灰度为Ul。
4.根据权利要求1所述乳腺X光图像中胸肌区域的分割方法，其特征在于，M、N都为3。
5.根据权利要求1所述乳腺X光图像中胸肌区域的分割方法，其特征在于，步骤C具体包括: Cl、令i = 0，1，…，GLmax为感兴趣区域ROI中的每一个灰度级，统计每个灰度级别的 在感兴趣区域ROI中所占的比例P⑴，给定初始阈值
6.根据权利要求1所述乳腺X光图像中胸肌区域的分割方法，其特征在于，步骤D包括: 使用高斯平滑模板对二值化图像ROIbw进行平滑处理； 使用Canny边缘算子提取平滑结果中所有的边缘，获得边界图像ROIedg-
7.根据权利要求1所述乳腺X光图像中胸肌区域的分割方法，其特征在于，步骤E包括: 使用霍夫变换在边界图像ROIedge中检测直线； 如果检测结果中包含多条直线，则选择斜率在1~3.7之间的且斜率最小的一条直线作为感兴趣区域ROI内的胸肌壁分界线； 如果霍夫变换仅仅检测出一条直线，那么就将这条直线作为感兴趣区域ROI的胸肌壁分界线。
【文档编号】G06T7/00GK103700085SQ201210367387
【公开日】2014年4月2日 申请日期:2012年9月28日 优先权日:2012年9月28日
【发明者】赵明, 康雨 申请人:深圳市蓝韵实业有限公司