一种乳腺钼靶图像的乳腺区域提取方法及系统与流程

1.本发明涉及医学图像处理技术领域，具体涉及一种乳腺钼靶图像的乳腺区域提取方法及系统。


背景技术：

2.乳腺钼靶成像，因为其成像分辨率高，易于操作，价格低廉且可以检查出医生触摸不到的肿块，对钙化灶也可以较好地显示，所以对于女性尤其是对致密性乳腺较多的女性来说，乳腺钼靶成像技术更适合进行早期的乳腺癌诊断工作。
3.目前乳腺钼靶图像检查主要还是由医生通过阅览乳腺钼靶图像得出诊断结果，不同的主治医生的经验丰富程度不同，而且长时间地阅片之后会产生一定程度的视觉疲劳，专注力相应减弱，因此会导致最终的检查结果出现一定的误差。尤其是面对早期乳腺癌的时候，肿块的特征比中晚期更难以观察到，也更容易存在漏检的现象。乳腺钼靶图像，由于其存在大量黑色背景，同时图像的对比度较差，这会对检测有较大的负面影响。
4.现阶段乳腺钼靶图像的乳腺区域的方法主要还是基于乳腺钼靶图像，利用深度学习技术构建的乳腺癌辅助诊断系统，虽然该方法已经取得不错的成效，其结果也获得了专业医生的认可。不过，这类系统需要大量的标注数据作为训练样本，而这些训练数据通常需要发射科医生来完成标注。医生标注需要凭借大量的临床经验，而且不同的医生会具备不同的评价标准，这一标注过程耗时耗力，当数据量特别大的时候，这项工作变成一个巨大的挑战。


技术实现要素：

5.为了解决上述技术问题，本发明提出了一种乳腺钼靶图像的乳腺区域提取方法，包括如下步骤：s1、获取乳腺钼靶图像，使用阈值分割算法将完整的乳腺区域从乳腺钼靶图像中分割出来；s2、使用截断归一化法对分割后的乳腺区域进行归一化处理，得到截断归一化图像；s3、使用自适应直方图均衡算法来增强截断归一化图像的边界，得到增强图像；s4、将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据，得到融合图像，并对融合图像进行镜像处理。
6.进一步地，步骤s2中，得到分割后的乳腺区域的灰度分布，对乳腺区域的灰度值进行排序，选择大于最小灰度值5%的灰度值作为灰度值下限，选择小于最大灰度值1%的灰度值作为灰度值上限，对乳腺区域中的每个灰度值m进行截断处理：；
对截断后的乳腺区域中的每个灰度值m进行归一化处理，得到归一化后的灰度值p：。
7.进一步地，步骤s3中，采用乳腺钼靶图像的原始图的累积分布函数作为变换函数，均衡化变换如下：；其中，为输出图像中的亮度值，为输入图像的亮度值，l为像素数，像素点序号j=1,2,
……
，l，为图像概率密度函数，。
8.进一步地，步骤s4包括如下步骤：s41、将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据；将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据r，融合算法如下式：
9.s42、对融合图像进行镜像处理；将融合后图像做水平翻转处理，将图像的左半部分和右半部分以图像垂直中轴线为中心进行镜像对换，设(x0，y0)为图像中像素点的坐标，(x0，y0)进行水平翻转后坐标将变为(w-x0，y0)，图像宽度为w，变换矩阵表达式如下式所示：。
10.进一步地，步骤s1中，阈值分割算法的具体步骤如下：设阈值为k，阈值k把乳腺钼靶图像中的所有像素灰度值分为两部分，两部分的像素灰度均值分别为,，图像的总体像素灰度均值为,两部分的像素数占总像素数的比例分别为,,则：；类间方差表示为：；求出使最大时对应的灰度均值，把灰度均值作为阈值k去分割乳腺钼靶图像。
11.本发明还提出了一种乳腺钼靶图像的乳腺区域提取系统，用于实现乳腺区域提取方法，包括：乳腺钼靶检查装置，乳腺区域提取装置，图像输出装置；所述乳腺钼靶检查装置，用于获取乳腺钼靶图像；所述乳腺区域提取装置，用于对乳腺钼靶图像中的乳腺区域进行提取并进行数据处理；图像输出装置，用于输出数据处理后的乳腺区域。
12.进一步地，乳腺区域提取装置包括：图像分割单元、归一化处理单元、边界增强单元和融合处理单元；所述图像分割单元，用于使用阈值分割算法将完整的乳腺区域从乳腺钼靶图像中分割出来；所述归一化处理单元，用于使用截断归一化法对分割后的乳腺区域进行归一化处理，得到截断归一化图像；所述边界增强单元，用于使用自适应直方图均衡算法来增强截断归一化图像的边界，得到增强图像；所述融合处理单元，用于将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据，得到融合图像，并对融合图像进行镜像处理。
13.相比于现有技术，本发明具有如下有益技术效果：本发明将乳腺区域和背景区域较好地分割开来，使用截断归一化法对分割后的乳腺区域进行归一化处理，得到截断归一化图像，改善了原始图像中乳腺边界模糊的问题，经过累积分布函数增强截断归一化图像的边界，使得图像中乳腺边界清晰度得到了进一步的提升，将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据，既保留了原始的数据分布信息，又可以提供乳腺边界信息，有利于取乳腺图像中特征数据的提取。
附图说明
14.为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
15.图1为本发明的乳腺钼靶图像的乳腺区域提取方法的流程示意图。
16.图2为本发明的乳腺钼靶图像的乳腺区域提取系统的结构示意图。
具体实施方式
17.为使本技术实施例的目的、技术方案和优点更加清楚，下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
18.在本发明的具体实施例附图中，为了更好、更清楚的描述系统中的各元件的工作原理，表现所述装置中各部分的连接关系，只是明显区分了各元件之间的相对位置关系，并不能构成对元件或结构内的信号传输方向、连接顺序及各部分结构大小、尺寸、形状的限定。
19.如图1所示，为本发明的乳腺钼靶图像的乳腺区域提取方法的流程示意图，包括如下步骤：s1、获取乳腺钼靶图像，使用阈值分割算法将完整的乳腺区域从乳腺钼靶图像中分割出来。
20.获取乳腺钼靶图像的方式为从每个患例中采集左乳cc位、左乳mlo位、右乳cc位和
右乳mlo位的乳腺钼靶图像。由于钼靶图像成像特性，图像除了乳腺部分外还有很大的黑色区域，其中带有成像信息，这部分区域是无用的，需要剔除掉这部分，保证乳腺部分充盈整副图像，乳腺钼靶图像尺度为672
×
448。
21.统计乳腺钼靶图像的像素灰度直方图，灰度分布用图的形式表现出来，像素灰度直方图中的两个峰，一个代表乳腺区域的像素分布，另一个代表除去乳腺区域的背景区域的灰度分布。在两个峰值之间必然存在一个谷值，把谷值设为图像分割的阈值，通过该谷值能够将乳腺区域和背景区域较好地分割开来。
22.在对乳腺钼靶图像进行阈值分割时，所选取的分割阈值应使乳腺区域的平均灰度、背景区域的平均灰度与整幅乳腺钼靶图像的平均灰度之间的差值最大，这种差异用类间方差来表示。本发明使用最大方差算法，来计算分割阈值。
23.阈值分割算法的具体步骤如下：设阈值为k，阈值k把乳腺钼靶图像中的所有像素灰度分为两部分，两部分的像素灰度均值分别为,，图像的总体像素灰度均值为,两部分的像素数占总像素数的比例分别为,,则：：；类间方差表示为：；简化得到：。
24.其中，求出使最大时对应的灰度均值，把灰度均值作为阈值k去分割乳腺钼靶图像。
25.s2、使用截断归一化法对分割后的乳腺区域进行归一化处理，得到截断归一化图像。
26.采用截断归一化法进行归一化处理的具体过程如下：s2.1、得到分割后的乳腺区域的灰度分布。
27.s2.2、对乳腺区域的灰度值进行排序，选择大于最小灰度值5%的灰度值作为灰度值下限，选择小于最大灰度值1%的灰度值作为灰度值上限，对乳腺区域中的每个灰度值mj进行截断处理，如下述公式：；其中，j为像素点序号，j=1,2,
……
，l。
28.截断处理避免了乳腺区域中的噪声和异常值产生的影响，最大程度保留有效的乳腺区域图像的灰度值分布。
29.s2.3、对截断后的乳腺区域中的每个灰度值mj进行如下公式所示的归一化处理，得到归一化后的灰度值p：
；截断归一化处理之后,可以最大程度的保留乳腺区域的像素差异性，使乳腺区域的图像更加清晰，病灶的边界信息也更明显，病灶和周围组织的对比度更高，更容易进行观察。
30.s3、使用自适应直方图均衡算法增强截断归一化图像的边界，得到增强图像。具体包括如下步骤：采用乳腺钼靶图像的原始图的累积分布函数作为变换函数，进行增强化变换如下：；其中，为输出亮度值，为输入的每个像素点处的亮度值，l为像素数，像素点序号j=1,2,
……
，l，为图像概率密度函数，。
31.经过增强处理，既保留了原始图像的数据分布信息，还增加了区域边界信息，为区域的定位和提取提供了充分的图像特征信息。
32.s4、将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据，得到融合图像，并对融合图像进行镜像处理。
33.为了给深度学习模型提供原始的图像数据分布信息，本发明将截断归一化图像和增强图像的像素值作为三个通道，进行数据处理，包括如下步骤：s41、将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据。
34.将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据r，融合算法如下式：
35.s42、对融合图像进行镜像处理。
36.通过将融合后图像做水平翻转处理，将图像的左半部分和右半部分以图像垂直中轴线为中心进行镜像对换。设 (x0，y0)为图像中像素点的坐标，(x0，y0)进行水平翻转后坐标将变为(w-x0，y0)，图像宽度为w，变换矩阵表达式如下式所示：；翻转前、后不影响其图像的融合数据。
37.如图2所示，为本发明的乳腺钼靶图像的乳腺区域提取系统的结构示意图，包括：乳腺钼靶检查装置，乳腺区域提取装置，图像输出装置。
38.乳腺钼靶检查装置，用于获取乳腺钼靶图像；乳腺区域提取装置，用于对乳腺钼靶图像中的乳腺区域进行提取并进行数据处理。乳腺区域提取装置包括：图像分割单元、归一化处理单元、边界增强单元和融合处理单元。
39.图像分割单元，用于使用阈值分割算法将完整的乳腺区域从乳腺钼靶图像中分割
出来；归一化处理单元，用于使用截断归一化法对分割后的乳腺区域进行归一化处理，得到截断归一化图像；边界增强单元，使用累积分布函数增强截断归一化图像的边界，得到增强图像；融合处理单元，用于将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据，得到融合图像，并对融合图像进行镜像处理。
40.图像输出装置，用于输出数据处理后的乳腺区域。
41.本发明将乳腺区域和背景区域较好地分割开来，使用截断归一化法对分割后的乳腺区域进行归一化处理，得到截断归一化图像，改善了原始图像中乳腺边界模糊的问题，经过累积分布函数增强截断归一化图像的边界，使得图像中乳腺边界清晰度得到了进一步的提升，将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据，既保留了原始的数据分布信息，又可以提供乳腺边界信息，有利于取乳腺图像中特征数据的提取。
42.在上述实施例中，可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用软件实现时，可以全部或部分地以计算机程序产品的形式实现。所述计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行所述计算机程序指令时，全部或部分地产生按照本技术实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中，或者通过所述计算机可读存储介质进行传输。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质，(例如，软盘、硬盘、磁带)、光介质(例如，dvd)、或者半导体介质(例如，固态硬盘(solid state disk，ssd))等。
43.以上所述，仅为本技术的具体实施方式，但本技术的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本技术揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本技术的保护范围之内。因此，本技术的保护范围应以权利要求的保护范围为准。

技术特征：
1.一种乳腺钼靶图像的乳腺区域提取方法，其特征在于，包括如下步骤：s1、获取乳腺钼靶图像，使用阈值分割算法将完整的乳腺区域从乳腺钼靶图像中分割出来；s2、使用截断归一化法对分割后的乳腺区域进行归一化处理，得到截断归一化图像；s3、使用自适应直方图均衡算法增强截断归一化图像的边界，得到增强图像；s4、将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据，得到融合图像，并对融合图像进行镜像处理。2.根据权利要求1所述的乳腺区域提取方法，其特征在于，步骤s2中，得到分割后的乳腺区域的灰度分布，对乳腺区域的灰度值进行排序，选择大于最小灰度值5%的灰度值作为灰度值下限，选择小于最大灰度值1%的灰度值作为灰度值上限，对乳腺区域中的每个灰度值m进行截断处理：；对截断后的乳腺区域中的每个灰度值m进行归一化处理，得到归一化后的灰度值p：。3.根据权利要求2所述的乳腺区域提取方法，其特征在于，步骤s3中，采用乳腺钼靶图像的原始图的累积分布函数作为变换函数，均衡化变换如下：；其中，为输出图像中的亮度值，为输入图像的亮度值，l为像素数，像素点序号j=1,2,
……
，l，为图像概率密度函数，。4.根据权利要求3所述的乳腺区域提取方法，其特征在于，步骤s4包括如下步骤：s41、将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据；将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据r，融合算法如下式：；s42、对融合图像进行镜像处理；将融合后图像做水平翻转处理，将图像的左半部分和右半部分以图像垂直中轴线为中心进行镜像对换，设(x0，y0)为图像中像素点的坐标，(x0，y0)进行水平翻转后坐标将变为(w-x0，y0)，图像宽度为w，变换矩阵表达式如下式所示：。5.根据权利要求1所述的乳腺区域提取方法，其特征在于，步骤s1中，阈值分割算法的具体步骤如下：设阈值为k，阈值k把乳腺钼靶图像中的所有像素灰度值分为两部分，两部分的像素灰度均值分别为,，图像的总体像素灰度均值为,两部分的像素数占总像素数的比例
分别为,,则：；类间方差表示为：；求出使最大时对应的灰度均值，把灰度均值作为阈值k去分割乳腺钼靶图像。6.一种乳腺钼靶图像的乳腺区域提取系统，用于实现如权利要求1-5任意一项所述的乳腺区域提取方法，其特征在于，包括：乳腺钼靶检查装置，乳腺区域提取装置，图像输出装置；所述乳腺钼靶检查装置，用于获取乳腺钼靶图像；所述乳腺区域提取装置，用于对乳腺钼靶图像中的乳腺区域进行提取并进行数据处理；图像输出装置，用于输出数据处理后的乳腺区域。7.根据权利要求6所述的乳腺区域提取系统，其特征在于，乳腺区域提取装置包括：图像分割单元、归一化处理单元、边界增强单元和融合处理单元；所述图像分割单元，用于使用阈值分割算法将完整的乳腺区域从乳腺钼靶图像中分割出来；所述归一化处理单元，用于使用截断归一化法对分割后的乳腺区域进行归一化处理，得到截断归一化图像；所述边界增强单元，用于使用自适应直方图均衡算法来增强截断归一化图像的边界，得到增强图像；所述融合处理单元，用于将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据，得到融合图像，并对融合图像进行镜像处理。

技术总结
本发明提出了一种乳腺钼靶图像的乳腺区域提取方法及系统，涉及医学图像处理技术领域，获取乳腺钼靶图像，使用阈值分割算法将完整的乳腺区域从乳腺钼靶图像中分割出来；使用截断归一化法对分割后的乳腺区域进行归一化处理，得到截断归一化图像；使用自适应直方图均衡算法增强截断归一化图像的边界，得到增强图像；将截断归一化图像的灰度值和增强图像的亮度值融合为融合数据，得到融合图像，并对融合图像进行镜像处理。既保留了原始的数据分布信息，又可以提供乳腺边界信息，有利于取乳腺图像中特征数据的提取。图像中特征数据的提取。图像中特征数据的提取。


技术研发人员：王健 陈贺凯 贾治华
受保护的技术使用者：天津市第五中心医院
技术研发日：2023.06.13
技术公布日：2023/7/17