乳房图像获取方法及装置、乳房摄影系统

申请号 CN201510640746.3 申请日 2015-09-30 公开(公告)号 CN105374025A 公开(公告)日 2016-03-02
申请人 上海联影医疗科技有限公司; 发明人 张娜; 周海华; 滕万里; 杨乐;
摘要 一种乳房图像获取方法及装置、乳房摄影系统。所述乳房图像获取方法包括:分割医学图像以获得乳房区域;基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘;基于所述乳房区域和所述乳房边缘确定所述乳房图像。本 发明 技术方案可以获取不同剂量的医学图像中的乳房图像,且最终获得的乳房图像的精确度较高,提高了乳房图像的 质量 ,进而也降低了误诊率,且在一定程度上也提高了医生对乳房区域病灶的诊断效率。
权利要求

1.一种乳房图像获取方法,其特征在于,包括:
分割医学图像以获得乳房区域;
基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘;
基于所述乳房区域和所述乳房边缘确定所述乳房图像。
2.如权利要求1所述的乳房图像获取方法,其特征在于,所述基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘包括:
对所述医学图像进行去噪以获得第一图像;
对所述第一图像进行预处理以使得待检测乳房边缘的灰度值最高,获得第二图像;
基于边缘检测方法对所述第二图像中的乳房边缘进行检测。
3.如权利要求2所述的乳房图像获取方法,其特征在于,对所述第一图像进行预处理包括:对所述第一图像进行梯度变换或差分运算。
4.如权利要求2所述的乳房图像获取方法,其特征在于,基于迭代法或大津法对所述第二图像中的乳房边缘进行检测。
5.如权利要求1所述的乳房图像获取方法,其特征在于,所述基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘包括:
对所述医学图像进行对数变换以获得第三图像;
对所述第三图像进行去噪以获得第四图像;
对所述第四图像进行预处理以使得待检测乳房边缘的灰度值最高,获得第五图像;
基于边缘检测方法对所述第五图像中的乳房边缘进行检测。
6.如权利要求5所述的乳房图像获取方法,其特征在于,对所述第四图像进行预处理包括:对所述第四图像进行梯度变换或差分运算。
7.如权利要求5所述的乳房图像获取方法,其特征在于,基于迭代法或大津法对所述第五图像中的乳房边缘进行检测。
8.如权利要求1所述的乳房图像获取方法,其特征在于,还包括去除所述乳房图像中存在的孤立像素点。
9.如权利要求1所述的乳房图像获取方法,其特征在于,基于大津法分割所述医学图像以获得乳房区域。
10.一种乳房图像获取装置,其特征在于,包括:
分割单元,用于分割医学图像以获得乳房区域;
检测单元,用于基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘;
确定单元,用于基于所述乳房区域和所述乳房边缘确定所述乳房图像。
11.如权利要求10所述的乳房图像获取装置,其特征在于,所述检测单元包括:
第一去噪单元,用于对所述医学图像进行去噪以获得第一图像;
第一预处理单元,用于对所述第一图像进行预处理以使得待检测乳房边缘的灰度值最高,获得第二图像;
第一检测子单元,用于基于边缘检测方法对所述第二图像中的乳房边缘进行检测。
12.如权利要求10所述的乳房图像获取装置,其特征在于,所述检测单元包括:
变换单元,用于对所述医学图像进行对数变换以获得第三图像;
第二去噪单元,用于对所述第三图像进行去噪以获得第四图像;
第二预处理单元,用于对所述第四图像进行预处理以使得待检测乳房边缘的灰度值最高,获得第五图像;
第二检测子单元,用于基于边缘检测方法对所述第五图像中的乳房边缘进行检测。
13.一种乳房摄影系统,其特征在于,包括权利要求10至12任一项所述的乳房图像获取装置。

说明书全文

乳房图像获取方法及装置、乳房摄影系统

技术领域

[0001] 本发明涉及图像处理技术领域,特别涉及一种乳房图像获取方法及装置、乳房摄影系统。

背景技术

[0002] 随着计算机科学和信息技术的发展,医学成像技术也得到了迅速的发展,各种医用影像系统不断涌现。全视野数字乳房X线摄影系统(FFDM,full-field digital mammography)作为乳腺癌的筛查和诊断得到了广泛的应用。
[0003] 对采用FFDM拍摄获得的医学图像而言,乳房图像的提取是基于FFDM进行摄片的计算机辅助诊断中的关键步骤。对乳房图像进行提取能够将病灶的检测范围限定在乳房区域之内,进而降低背景区域的干扰。
[0004] 目前,通常采用以下三种方法对乳房图像进行分割:
[0005] 第一种:大津法,利用阈值分割原图像中的前景图像和背景图像。采用大津法实现对乳房图像的分割较简单,但是分割后获得的乳房图像不精确,有所缺失。
[0006] 第二种:分算法,相对于大津法而言采用分水岭算法分割后获得的乳房图像比采用大津法分割获得的乳房图像较精确,但是对于剂量较高的医学图像,采用分水岭算法分割后获得的乳房图像仍不精确。
[0007] 第三种:Snake算法,采用Snake算法对医学图像进行分割,分割获得的乳房图像较精确,但其分割速度较慢,不适合用于实际的产品中。
[0008] 因此,如何能够提供一种精确度高且速度快的乳房图像获取方法,成为目前亟待解决的问题之一。

发明内容

[0009] 本发明要解决的问题是提供一种乳房图像获取方法及装置、乳房摄影系统,以使得获得的乳房图像精度高且获得所述乳房图像的速度快。
[0010] 为解决上述问题,本发明技术方案提供一种乳房图像获取方法,包括:
[0011] 分割医学图像以获得乳房区域;
[0012] 基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘;
[0013] 基于所述乳房区域和所述乳房边缘确定所述乳房图像。
[0014] 可选的,所述基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘包括:
[0015] 对所述医学图像进行去噪以获得第一图像;
[0016] 对所述第一图像进行预处理以使得待检测乳房边缘的灰度值最高,获得第二图像;
[0017] 基于边缘检测方法对所述第二图像中的乳房边缘进行检测。
[0018] 可选的,对所述第一图像进行预处理包括:对所述第一图像进行梯度变换或差分运算。
[0019] 可选的,基于迭代法或大津法对所述第二图像中的乳房边缘进行检测。
[0020] 可选的,所述基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘包括:
[0021] 对所述医学图像进行对数变换以获得第三图像;
[0022] 对所述第三图像进行去噪以获得第四图像;
[0023] 对所述第四图像进行预处理以使得待检测乳房边缘的灰度值最高,获得第五图像;
[0024] 基于边缘检测方法对所述第五图像中的乳房边缘进行检测。
[0025] 可选的,对所述第四图像进行预处理包括:对所述第四图像进行梯度变换或差分运算。
[0026] 可选的,基于迭代法或大津法对所述第五图像中的乳房边缘进行检测。
[0027] 可选的,所述乳房图像获取方法,还包括去除所述乳房图像中存在的孤立像素点。
[0028] 可选的,基于大津法分割所述医学图像以获得乳房区域。
[0029] 为解决上述问题,本发明技术方案还提供一种乳房图像获取装置,包括:
[0030] 分割单元,用于分割医学图像以获得乳房区域;
[0031] 检测单元,用于基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘;
[0032] 确定单元,用于基于所述乳房区域和所述乳房边缘确定所述乳房图像。
[0033] 可选的,所述检测单元包括:
[0034] 第一去噪单元,用于对所述医学图像进行去噪以获得第一图像;
[0035] 第一预处理单元,用于对所述第一图像进行预处理以使得待检测乳房边缘的灰度值最高,获得第二图像;
[0036] 第一检测子单元,用于基于边缘检测方法对所述第二图像中的乳房边缘进行检测。
[0037] 可选的,所述检测单元包括:
[0038] 变换单元,用于对所述医学图像进行对数变换以获得第三图像;
[0039] 第二去噪单元,用于对所述第三图像进行去噪以获得第四图像;
[0040] 第二预处理单元,用于对所述第四图像进行预处理以使得待检测乳房边缘的灰度值最高,获得第五图像;
[0041] 第二检测子单元,用于基于边缘检测方法对所述第五图像中的乳房边缘进行检测。
[0042] 为解决上述问题,本发明技术方案还提供一种乳房摄影系统,包括上述的乳房图像获取装置。
[0043] 与现有技术相比,本发明技术方案具有以下优点:
[0044] 先分割医学图像以获得乳房区域,然后基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘,进而基于所述乳房区域和所述乳房边缘确定所述乳房图像,由于先分割出乳房区域且考虑了实际医学图像中乳房边缘的灰度变化特性来检测所述乳房边缘,因此对于不同剂量的医学图像而言,最终获得的乳房图像的精确度较高,减少了背景区域的干扰,提高了乳房图像的质量,进而也降低了误诊率,在一定程度上也提高了医生对乳房区域病灶的诊断效率。
[0045] 进一步地,对所述医学图像进行去噪以获得第一图像,然后对所述第一图像进行预处理获得待检乳房边缘灰度值最高的第二图像,避免了噪声对所述第二图像的影响,实现了待检测乳房边缘的灰度值最高,乳房区域的灰度值次高、背景区域的灰度值最低。故使得后续采用边缘检测方法对所述第二图像中的乳房边缘进行检测时,减少将背景区域的像素点检测为乳房区域的像素点的概率,进一步的提高了最终获得的乳房图像的精确度,提高了获得的乳房图像的质量。
[0046] 进一步地,考虑到X射线的衰减特性服从指数分布,在对所述医学图像去噪前,对所述医学图像进行对数变换获得第三图像,便于后续的处理,在一定程度上降低了图像处理过程中计算的复杂度和计算量,进而也提高了获得乳房图像的速度。
[0047] 进一步地,采用大津法对所述医学图像进行分割以获得乳房区域,由于乳房图像的灰度直方图为双峰图,故采用大津法分割所述医学图像以获得乳房区域的速度快且精度高,同时由于考虑了实际医学图像中乳房边缘的灰度变化特性来检测所述乳房边缘,因此对于不同剂量的医学图像而言,最终获得的乳房图像的精确度较高,进而提高了乳房图像的质量,在降低了误诊率的同时也提高了医生对乳房区域的病灶的诊断效率。此外采用迭代法检测所述第二图像中的乳房边缘,计算量小,因此进一步的提高了获得乳房图像的速度。附图说明
[0048] 图1是本发明实施方式的乳房图像获取方法的流程示意图;
[0049] 图2是FFDM拍摄的医学图像;
[0050] 图3是本发明实施例的乳房图像获取方法的流程示意图;
[0051] 图4是采用本发明实施例的乳房图像获取方法获得的乳房图像;
[0052] 图5是本发明实施方式的乳房图像获取装置的结构示意图。

具体实施方式

[0053] 为使本发明的上述目的、特征和优点能够更为明显易懂,下面结合附图对本发明的具体实施方式做详细的说明。在以下描述中阐述了具体细节以便于充分理解本发明。但是本发明能够以多种不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广。因此本发明不受下面公开的具体实施方式的限制。
[0054] 正如现有技术中所描述的,目前所采用的对乳房图像进行分割的算法,分割后获得的乳房图像的精确度低,不符合实际的临床需求,若要获得精确度较高的乳房图像,则获得乳房图像的速度较慢,不能应用在实际的产品中。因此,现有的乳房图像分割算法无法实现以较快的速度获得精确度较高的乳房图像。
[0055] 发明人对FFDM获得的医学图像进行分析后得出,对于实际拍摄获得的医学图像而言,乳房边缘的灰度值与背景区域的灰度值比较接近,尤其是对于高剂量的医学图像,其乳房边缘和背景区域基本上难以区分,进而导致采用分割算法分割乳房图像(包括乳房区域和乳房边缘)时,乳房边缘会被分割到背景区域中,出现过分割的现象,分割获得的乳房图像的精确度较低。
[0056] 因此,发明人考虑先将医学图像中的乳房区域进行提取,然后基于乳房边缘的灰度变化特性来检测乳房边缘,将获得的乳房区域和乳房边缘进行叠加,以获得乳房图像。
[0057] 参见图1,图1是本发明实施方式的乳房图像获取方法的流程示意图,如图1所示,本发明实施方式的乳房图像获取方法包括:
[0058] S101:分割医学图像以获得乳房区域;
[0059] S102:基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘;
[0060] S103:基于所述乳房区域和所述乳房边缘确定所述乳房图像。
[0061] 以下结合具体的实施例对本发明实施方式的乳房图像获取方法进行详细的说明。
[0062] 执行S101,对通过FFDM采集到的医学图像进行分割,获得乳房区域,考虑到乳房图像的灰度直方图为双峰图,采用大津法分割在速度和精度上最为有效,因此本实施例中先获取所述医学图像的灰度直方图,然后采用大津法对所述医学图像进行分割,分割后获得的乳房区域精度高,且分割速度较快。在其他实施例中,也可以采用其他的图像分割算法对所述医学图像中的乳房区域进行分割,如:采用分水岭算法对所述医学图像进行分割。
[0063] 执行S102,基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘。由上述可知对于临床获得的医学图像而言,乳房边缘的灰度值和背景区域的灰度值非常接近,此外,乳房边缘的灰度值又是所述医学图像中灰度变化最剧烈的部分,因此,可以先抬高所述医学图像中待检测乳房边缘的灰度值,使得所述待检测乳房边缘在所述医学图像中凸显出来以便于后续对所述乳房边缘的检测。此外,考虑到通过FFDM采集到的医学图像,不同剂量下采集到的医学图像的噪声各不相同,且噪声的存在势必会对提高待检测乳房边缘的灰度值有所影响,因此在提高所述医学图像中待检测乳房边缘的灰度值之前对所述医学图像进行去噪,避免噪声对后续处理的影向。
[0064] 此外,考虑到X射线的衰减服从指数分布,为了降低计算复杂度和计算量进而进一步的提高检测乳房边缘的速度,本实施例中在对所述医学图像去噪前,对所述医学图像先进行对数变换以获得第三图像,然后对所述第三图像进行去噪以获得第四图像。
[0065] 本实施例中,采用高斯滤波的方式对所述医学图像去噪,在其他实施例中也可以采用均值滤波的方式对所述医学图像去噪。
[0066] 在对所述医学图像进行对数变换和去除了经过对数变换的医学图像中的噪声以获得所述第四图像后,对所述第四图像进行预处理以使得待检测乳房边缘的灰度值最高,也即抬高所述第四图像中待检测乳房边缘的灰度值,使得所述待检测乳房边缘在所述第四图像中凸显出来以便于后续对所述乳房边缘的检测。本实施例中所述预处理包括对所述第四图像进行梯度变换,具体的通过sobel梯度算子对所述第四图像进行梯度变换。在其他实施例中,也可以对所述第四图像进行差分运算,也即采用差分法来提高所述第四图像中待检测乳房边缘的灰度值。
[0067] 本领域技术人员知晓,对于采用FFDM拍摄获得的医学图像而言,其乳房区域的灰度值通常低于背景区域的灰度值,即为正片,如图2所示,正片中背景区域的灰度值>乳房边缘的灰度值>乳房区域的灰度值,对所述正片(医学图像)进行对数变换后获得第三图像,再对所述第三图像进行去噪可以使得所述第三图像中灰度值波动较大的背景区域的灰度波动变小,或者说可以使得所述第三图像的背景区域中所有像素点的灰度值趋于同一水平。因此,在对去噪后的第三图像即第四图像进行梯度变换时,可以使得获得的第五图像中待检测乳房边缘的灰度值>乳房区域的灰度值>背景区域的灰度值。故,可以避免出现在采用边缘检测方法提取所述第五图像中的乳房边缘时,若背景区域的灰度值较高,可能会将背景区域中的像素点检测为乳房边缘的像素点的现象,因此,对进行了对数变换的医学图像先进行去噪再进行梯度变换,在一定程度上为后续获取的乳房图像的精确度提供了一定的保证,有利于提高乳房图像的精确度。
[0068] 在对所述医学图像进行了对数变换、去噪和梯度变换获得第五图像后,基于边缘检测方法对所述第五图像中的乳房边缘进行检测。考虑到实际应用中获得的第五图像为浮点型的数据,而采用迭代法对乳房边缘进行检测时,不需要将浮点型的数据转换为整型数据故图像精度不会有所损失,此外迭代法计算量较小且速度快,因此本实施例中,采用迭代法检测所述第五图像中的乳房边缘。具体地,
[0069] 设置初始阈值T0,为了加快迭代速度,本实施例中以所述第五图像中最大灰度值和最小灰度值之和的一半为所述初始阈值T0,第一次迭代时,根据所述初始阈值T0,所述第五图像被分为两个区域,像素点灰度值大于所述初始阈值T0的第一区域和像素点灰度值小于等于所述初始阈值T0的第二区域,计算所述第一区域和第二区域的灰度均值(各区域像素点灰度值之和除以对应的像素点个数),以所述第一区域和第二区域的灰度均值之和的一半T1为新的阈值以更新初始阈值T0,判断abs(T1-T0)是否大于1,若abs(T1-T0)>1,则以T1为阈值继续对所述第五图像进行分割,直至abs(Ti+1-Ti)<1,迭代终止,Ti+1为最终的分割阈值,像素点灰度值大于Ti+1的像素点的集合为乳房边缘。
[0070] 在其他实施例中,也可以采用大津法或其他分割方法对第二图像中的乳房边缘进行提取。
[0071] 此外,在实际迭代过程中,还可以对迭代过程中获得的阈值乘以小于1的权重因子,所述权重因子可以根据对多幅医学图像的乳房边缘进行检测时的实际迭代过程中迭代收敛前阈值的变化规律而定。对迭代过程中获得的阈值乘以小于1的权重因子,可以进一步的提高获取的乳房边缘的精确度。
[0072] 执行S103,基于所述乳房区域和所述乳房边缘确定所述乳房图像。本实施例中,具体地就是将S101获得的乳房区域和S102获得的乳房边缘进行叠加以获得最终的乳房图像。考虑到实际处理过程中,可能会有背景区域的像素点被误判为乳房图像的像素点,因此需对叠加后获得的乳房图像中的孤立像素点进行去除,以提高乳房图像的精确度。
[0073] 以下结合实际的应用和图3对本发明实施例的乳房图像获取方法进行简单的说明。
[0074] 如图3所示,输入原始投影图像(RawImage,即FFDM获得的医学图像),计算RawImage的灰度直方图,采用大津法(Otsu)分割RawImage以获得乳房区域。对RawImage先进行log变换再进行高斯滤波得到Image_1,对Image_1进行梯度变换得到Image_2,采用迭代法检测所述Image_2中的乳房边缘,将检测得到的乳房区域和乳房边缘进行叠加以获得乳房图像,去除所述乳房图像中的孤立像素点以输出最终的乳房图像。
[0075] 参见图4,图4是采用本发明实施例的乳房图像获取方法获得的乳房图像,图中白色区域所示即为最终获得的乳房图像,从图4中可以看出采用本发明实施例的乳房图像获取方法获得的乳房边缘清晰,可以精确的将乳房图像从医学图像中分割出。
[0076] 本实施例中,采用大津法分割乳房区域,通过对医学图像进行对数变换、去噪和梯度变换获得所述第五图像,并通过迭代法对所述第五图像中的乳房边缘进行检测,对分割后的乳房区域和检测获得的乳房边缘进行叠加所获得的乳房图像的精确度较高,且采用大津法分割乳房区域,迭代法检测乳房边缘速度快,另外采用对数变换可以降低图像处理过程中的复杂度和计算量,进而也加快了检测乳房边缘的速度。此外,去除获得的乳房图像中的孤立像素点进一步地提高了最终获得的乳房图像的精确度,进而也提高了乳房图像的质量,在降低了误诊率的同时也提高了医生对乳房区域病灶的诊断效率。
[0077] 本发明实施方式还提供一种乳房获取装置,请参见图5,图5是本发明实施方式的乳房图像获取装置的结构示意图,如图5所示所述乳房获取装置1包括:
[0078] 分割单元10,用于分割医学图像以获得乳房区域;
[0079] 检测单元11,用于基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘;
[0080] 确定单元12,用于基于所述乳房区域和所述乳房边缘确定所述乳房图像。
[0081] 本实施例中,所述检测单元11包括:
[0082] 变换单元110,用于对所述医学图像进行对数变换以获得第三图像;
[0083] 第二去噪单元111,用于对所述第三图像进行去噪以获得第四图像;
[0084] 第二预处理单元112,用于对所述第四图像进行预处理以使得待检测乳房边缘的灰度值最高,获得第五图像;
[0085] 第二检测子单元113,用于基于边缘检测方法对所述第五图像中的乳房边缘进行检测。
[0086] 在其他实施例中,所述检测单元也可以包括:
[0087] 第一去噪单元,用于对所述医学图像进行去噪以获得第一图像;
[0088] 第一预处理单元,用于对所述第一图像进行预处理以使得待检测乳房边缘的灰度值最高,获得第二图像;
[0089] 第一检测子单元,用于基于边缘检测方法对所述第二图像中的乳房边缘进行检测。
[0090] 所述乳房图像获取装置的具体实施可参考所述乳房图像获取方法的实施,在此不再赘述。
[0091] 本实施例还提供一种乳房摄影系统,包括上述的乳房图像获取装置。
[0092] 综上所述,本发明实施方式提供的乳房图像获取方法,至少具有如下有益效果:
[0093] 先分割医学图像以获得乳房区域,然后基于乳房边缘的灰度变化特性检测所述医学图像中的乳房边缘,进而基于所述乳房区域和所述乳房边缘确定所述乳房图像,由于先分割出乳房区域且考虑了实际医学图像中乳房边缘的灰度变化特性来检测所述乳房边缘,因此对于不同剂量的医学图像而言,最终获得的乳房图像的精确度较高,减少了背景区域的干扰,提高了乳房图像的质量,进而也降低了误诊率,在一定程度上也提高了医生对乳房区域病灶的诊断效率。
[0094] 进一步地,对所述医学图像进行去噪以获得第一图像,然后对所述第一图像进行预处理获得待检乳房边缘灰度值最高的第二图像,避免了噪声对所述第二图像的影响,实现了待检测乳房边缘的灰度值最高,乳房区域的灰度值次高、背景区域的灰度值最低。故使得后续采用边缘检测方法对所述第二图像中的乳房边缘进行检测时,减少将背景区域的像素点检测为乳房区域的像素点的概率,进一步的提高了最终获得的乳房图像的精确度,提高了获得的乳房图像的质量。
[0095] 进一步地,考虑到X射线的衰减特性服从指数分布,在对所述医学图像去噪前,对所述医学图像进行对数变换获得第三图像,便于后续的处理,在一定程度上降低了图像处理过程中计算的复杂度和计算量,进而也提高了获得乳房图像的速度。
[0096] 进一步地,采用大津法对所述医学图像进行分割以获得乳房区域,由于乳房图像的灰度直方图为双峰图,故采用大津法分割所述医学图像以获得乳房区域的速度快且精度高,同时由于考虑了实际医学图像中乳房边缘的灰度变化特性来检测所述乳房边缘,因此对于不同剂量的医学图像而言,最终获得的乳房图像的精确度较高,进而提高了乳房图像的质量,在降低了误诊率的同时也提高了医生对乳房区域的病灶的诊断效率。此外采用迭代法检测所述第二图像中的乳房边缘,计算量小,因此进一步的提高了获得乳房图像的速度。
[0097] 本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。
QQ群二维码
意见反馈