一种乳腺x光图像中乳头区域的分割算法
【技术领域】
[0001]本发明涉及一种医学图像中的分割算法，尤其是涉及一种乳腺X光图像中乳头区域的分割算法。
【背景技术】
[0002]在乳腺X光图像中，乳头区域因亮度与背景差别相对较小，在乳腺增强中经常被模糊。乳头分割可为乳头区域增强创造先决条件，帮助临床医生更好地诊断病情。目前的乳头分割算法比较少，一些乳腺相关的分割算法只能检测出乳头区域，但不具备乳头区域分割功能。在临床上，乳头区对于医生的诊断具有参考价值，因此乳头的分割是有必要的，且对分割精细度要求较高。
[0003]申请号为ZL95192508.3、公开号CN1145674A的中国专利申请，提出了一种数字乳房X射线照片中的自动分割及轮廓线和乳头的检测，通过检测轮廓线把数字乳房X射线图分割成乳腺区和背景区，较为有效的检测到乳头区域。然而，此方法不能够完成乳头区域的分割工作。
[0004]申请号为ZL201110245860.8、公开号为CN102956035A的中国专利申请，提出一种用于乳腺X线图像中提取乳腺区域的预处理方法及系统，此发明通过全局阈值和局部阈值法得到光滑的乳房轮廓，能够完整检测到乳头区域，并很好地去除背景杂物干扰。然而，此发明并没有公开乳头区域的分割算法。
【发明内容】
[0005]本发明是申请人于2011年8月25日提出的名称为《用于乳腺X线图像中提取乳腺区域的预处理方法及系统》、申请号为ZL201110245860.8的中国发明专利申请的衍生发明，是该发明的后续算法，需要在乳腺分割的基础上完成。
[0006]本发明在乳房区域分割完整的基础上，沿乳房边缘获得图像局部特性并确定乳头点，然后通过罗氏定理获得乳头与皮肤的交接点，最后通过曲线拟合的方式得到乳头与乳腺组织的分割线，并将乳头区分离出来。
[0007]本发明采用以下技术方案:一种乳腺X光图像中乳头区域的分割算法，包括如下步骤:
[0008](a)预处理:将图像旋转至乳房靠右且胸肌壁在右上角，对图像进行分割，获取包含完整的乳头的乳房区域分割图noBackgroundSrc ；
[0009](b)确定乳头点:沿着所述分割图noBackgroundSrc中的乳房边缘，以边缘点位置为中心逐个取方形小块，计算每个所述方形小块中的像素亮度总和，确定像素亮度最大的所述方形小块的中心为乳头点的位置；
[0010](C)截取nipblock:以所述乳头点为中心,在所述分割图noBackgroundSrc上截取一方形区域nipblock,并将其中组织取边缘,所述方形区域nipblock中包含乳头与乳房皮肤的两个交接点，两个所述交接点为上拐点和下拐点；
[0011](d)寻找乳头拐点:将所述乳头点分别与上下端点相连，得到两条直线，以所述乳头点为隔断分为上下两部分，上部分点计算与靠上直线的距离，得到最大距离的点即定义为所述上拐点；下部分点计算与靠下直线的距离，得到最大距离的点即定义为所述下拐占.
[0012](e)曲线拟合:分别取所述上拐点之前、所述下拐点之后的若干个点进行曲线拟合，得到所述上拐点与所述下拐点之间的曲线；
[0013](f)标记乳头区域:将所述曲线左侧标记为乳头区域，右侧标记为乳腺区域，并输出最后的图像。
[0014]优选地，所述步骤(C)中还包括将边缘线去噪，并使用边缘追踪法记录下所有边缘点位置。
[0015]优选地，所述步骤(d)中所述上部分点与所述下部分点均只取直线右边的点进行计算其与对应直线的距离。
[0016]优选地，所述步骤(e)中所述曲线拟合为二次多项式曲线拟合。
[0017]优选地，所述步骤(f)中输出最后的图像之前还包括将所述乳头区域填充为灰色，将所述乳腺区域填充为白色，将所述背景区域填充为黑色。
[0018]与现有技术相比，本发明具有如下有益效果:本发明利用乳头区域的像素特征定位乳头位置，然后通过寻找上下拐点和二次多项式曲线拟合的方式将乳头区域与乳腺区域分离。此算法为改善图像增强中乳头区的模糊化问题打下了基础，方便医生诊断，提高诊断效率，并有利于最终医学诊断。
【附图说明】
[0019]图1为本发明一种乳腺X光图像中乳头区域的分割算法的流程图。
【具体实施方式】
[0020]以下结合附图1对本发明一种乳腺X光图像中乳头区域的分割算法做进一步的详细描述，应该指出的是，所列举的实施例不应理解对发明的限制。
[0021]本发明提出一种乳腺X光图像中乳头区域的分割算法，包括如下步骤:
[0022](a)预处理:将图像旋转至标准位置，即乳房靠右且胸肌壁在右上角，对图像进行分割，获取包含完整的乳头的乳房区域分割图noBackgroundSrc ；
[0023](b)确定乳头点:沿着所述分割图noBackgroundSrc中的乳房边缘，以边缘点位置为中心逐个取一定宽度的方形小块，计算每个所述方形小块中的像素亮度总和，确定像素亮度最大的所述方形小块的中心为乳头点的位置；因为乳头区域在图中属于较为明亮的区域，因此像素亮度总和最大的方形小块中必定包含乳头区域，其中心即为乳头点的位置；
[0024](c)截取nipblock:以所述乳头点为中心,在所述分割图noBackgroundSrc上截取一方形区域nipblock,并将其中组织取边缘,所述方形区域nipblock中包含乳头与乳房皮肤的两个交接点，两个所述交接点为上拐点和下拐点；将边缘线去噪，并使用边缘追踪法记录下所有边缘点位置，可以防止遗漏拐点；
[0025]Cd)寻找乳头拐点:使用罗氏定理在所述方形区域nipblock中寻找拐点，将所述乳头点分别与上下端点相连，得到两条直线，以所述乳头点为隔断分为上下两部分，上部分点计算与靠上直线的距离，得到最大距离的点即定义为所述上拐点；下部分点计算与靠下直线的距离，得到最大距离的点即定义为所述下拐点；为了防止出现错误，所述上部分点与所述下部分点均只取直线右边的点进行计算其与对应直线的距离；
[0026](e)曲线拟合:找到拐点之后，分别取所述上拐点之前、所述下拐点之后的若干个点进行曲线拟合，得到所述上拐点与所述下拐点之间的曲线，算法中采取二次多项式曲线拟合，其目的是在所述上拐点与下拐点之间得到一条较为自然的连接线以分割乳头区域，而不是简单地用一条直线隔开；
[0027](f)标记乳头区域:将所述曲线左侧标记为乳头区域，右侧标记为乳腺区域，并将所述乳头区域填充为灰色，将所述乳腺区域填充为白色，将所述背景区域填充为黑色，最终输出一幅三色图。
[0028]显然，上述内容只是为了说明本发明的特点，而并非对发明的限制，有关技术领域的普通技术人员根据本发明在相应的技术领域做出的变化应属于本发明的保护范畴。
【主权项】
1.一种乳腺X光图像中乳头区域的分割算法，其特征在于，包括如下步骤: Ca)预处理:将图像旋转至乳房靠右且胸肌壁在右上角，对图像进行分割，获取包含完整的乳头的乳房区域分割图noBackgroundSrc ； (b)确定乳头点:沿着所述分割图noBackgroundSrc中的乳房边缘，以边缘点位置为中心逐个取方形小块，计算每个所述方形小块中的像素亮度总和，确定像素亮度最大的所述方形小块的中心为乳头点的位置； (c)截取nipblock:以所述乳头点为中心,在所述分割图noBackgroundSrc上截取一方形区域nipblock，并将其中组织取边缘，所述方形区域nipblock中包含乳头与乳房皮肤的两个交接点，两个所述交接点为上拐点和下拐点； Cd)寻找乳头拐点:将所述乳头点分别与上下端点相连，得到两条直线，以所述乳头点为隔断分为上下两部分，上部分点计算与靠上直线的距离，得到最大距离的点即定义为所述上拐点；下部分点计算与靠下直线的距离，得到最大距离的点即定义为所述下拐点； Ce)曲线拟合:分别取所述上拐点之前、所述下拐点之后的若干个点进行曲线拟合，得到所述上拐点与所述下拐点之间的曲线； Cf)标记乳头区域:将所述曲线左侧标记为乳头区域，右侧标记为乳腺区域，并输出最后的图像。2.根据权利要求1所述的乳腺X光图像中乳头区域的分割算法，其特征在于，所述步骤(c)中还包括将边缘线去噪，并使用边缘追踪法记录下所有边缘点位置。3.根据权利要求1所述的乳腺X光图像中乳头区域的分割算法，其特征在于，所述步骤Cd)中所述上部分点与所述下部分点均只取直线右边的点进行计算其与对应直线的距离。4.根据权利要求1所述的乳腺X光图像中乳头区域的分割算法，其特征在于，所述步骤Ce)中所述曲线拟合为二次多项式曲线拟合。5.根据权利要求1所述的乳腺X光图像中乳头区域的分割算法，其特征在于，所述步骤(f)中输出最后的图像之前还包括将所述乳头区域填充为灰色，将所述乳腺区域填充为白色，将所述背景区域填充为黑色。
【专利摘要】本发明公开一种乳腺X光图像中乳头区域的分割算法，包括如下步骤：（a）预处理；（b）确定乳头点；（c）截取nipblock；（d）寻找乳头拐点；（e）曲线拟合；（f）标记乳头区域；本发明利用乳头区域的像素特征定位乳头位置，然后通过寻找上下拐点和二次多项式曲线拟合的方式将乳头区域与乳腺区域分离。此算法为改善图像增强中乳头区的模糊化问题打下了基础，方便医生诊断，提高诊断效率，并有利于最终医学诊断。
【IPC分类】G06T7/60, G06T5/00
【公开号】CN104978753
【申请号】CN201410130621
【发明人】闫桂霞, 赵明
【申请人】深圳市蓝韵实业有限公司
【公开日】2015年10月14日
【申请日】2014年4月2日