一种分割乳腺病灶的方法
技术领域
本发明涉及医学图像处理领域,尤其涉及一种分割乳腺病灶的方法。
背景技术
计算机辅助诊断(Computer Aided Diagnosis,CAD)是指通过影像学、医学图像处理技术以及其它可能的生理、生化手段,结合计算机的分析计算,辅助影像科医师发现病灶,提高诊断的准确率。随着技术的发展,以核磁共振成像、计算机断层扫描为代表的现代高清晰影像设备为临床疾病的诊断提供了极大方便。目前,CAD研究较为成熟的是在乳腺和肺部病变方面;
乳腺癌是女性最常见的癌症之一,是导致女性死亡的第二大癌症杀手。核磁共振成像因其高敏感度而成为辅助医生检测乳腺癌的有效手段。基于核磁共振成像的乳腺病灶分割技术对于医生了解病灶的形态,大小等重要特性扮演着非常重要的角色。
现有技术中,大量乳腺病灶分割的研究都是在用户介入(例如手动提供种子点)的情况下进行的半自动分割技术。有的是基于用户给出的灰度阈值来分析病灶的二维边缘和形态特征;有的在用户给出病灶内一个种子点后进行病灶的三维分割;有的是在用户定义的感兴趣区域内进行病灶分割。然而无论是给出感兴趣区域还是给出病灶内一个种子点,这些用户手动介入的过程都是非常耗时的。一种无需用户介入的全自动病灶分割技术可以节省大量的人力和时间,对乳腺病灶的分析有着非常重大的意义。一种不包括精确分割过程的全自动病灶的检测技术,利用细胞神经网络的方法先对乳腺区域进行划分,然后再对划分好的乳腺区域内利用模板的方法对病灶可进行全自动的检测。然而利用细胞神经网络的方法有时不能准确地划分乳腺区域,高亮度的心脏部分有时会被归入乳腺区域,而且此方法不适用于抑脂的乳腺图像。再者利用模板的检测病灶方法由于模板自身的局限性,可能不适用于边缘增强内部较暗的病灶类型以及其他特殊病灶模态。
发明内容
本发明解决的问题是提供一种分割乳腺病灶的方法,用以解决现有技术中必须在用户手动给出种子点的情况下,高效、准确地自动分割处理一个以上病灶的问题。
为了解决上述技术问题,本发明提供了一种分割乳腺病灶的方法,包括如下步骤:
(1)对注射造影剂前的磁共振图像基于椭圆模型进行胸部粗分割,获得第一图像;
(2)对所述注射造影剂前、后的磁共振图像进行剪影,在获得剪影图像上基于所述第一图像提取乳腺部分的剪影图像;
(3)对每层乳腺部分的剪影图像进行病灶检测,获得第二图像;
(4)在所述第二图像上依次找到三维连通域,依次对每个连通域以重心为种子点在所述剪影图像上进行自适应的区域增长,获得分割后的乳腺可疑病灶的磁共振图像。
进一步地,所述胸部粗分割前进行滤波器去噪处理和闭操作。
进一步地,所述滤波器去噪过程利用中值滤波器、高斯滤波器、高通滤波器、均值滤波器或低通滤波器。
进一步地,所述第一图像的取得包括如下步骤:
(1)对所述注射造影剂前的磁共振图像从上到下逐行进行检测,当检测到连通域的范围达到阈值M时,停止检测并记录行数N′;
(2)以行数N1=N′+n作为椭圆的起始行,以所述注射造影剂前的磁共振图像的最后一行作为椭圆的终止行为行数N2,进行椭圆分割处理;
(3)对椭圆分割后的图像从上到下逐行检测,当测得的灰度值为非零值时,停止检测,记录行数N′1作为乳腺部分的起始行,并记录所述非零值的像素点所在列数L;
(4)对椭圆分割后的图像从下到上逐行检测第L列,当测得的灰度值为非零值时,停止检测,记录行数N′2作为乳腺部分的终止行,以所述起始行N′1和终止行N′2提取所述乳腺部分,获得第一图像;
所述阈值M的取值范围为图像宽度的40~60%,所述行数n取自然数,取值范围为5~20行。
进一步地,所述剪影步骤前对所述注射造影剂前、后的磁共振图像进行配准操作。
进一步地,所述注射造影剂后的磁共振图像为造影剂初始增强的磁共振图像。
进一步地,所述病灶检测的方法包括如下步骤:
(1)对每层乳腺部分的剪影图像根据灰度直方图进行阈值分割;
(2)在阈值分割后的图像上进行检测,若三维连通区域的体积小于阈值b,则去除所述三维连通区域;
(3)在阈值分割后的图像上,依次得到三维连通域,对每个连通域进行二维投影,对所述二维投影进行检测,当满足时,则去除所述三维连通域,获得第二图像;
所述lmax为所述二维投影最长轴的长度,所述S为该连通域二维投影的面积,所述阈值b的取值范围为10~50立方毫米,所述c为倍数,其取值范围为0.9~1.2。
进一步地,所述阈值分割方法包括如下步骤:
(1)对每层乳腺部分的剪影图像计算灰度直方图,根据所述灰度直方图分别定义第k层的高阈值Tk1和低阈值Tk2;
(2)若第k层图像的高阈值Tk1<a×Tave,则以高阈值Tk1对该层进行阈值分割;
(3)若第k层图像的高阈值Tk1≥a×Tave,则以低阈值Tk2对该层进行阈值分割;
所述k取自然数,所述阈值Tave为各层剪影图的高阈值平均值,所述a为倍数,取值范围为0.9~1.2。
进一步地,其特征在于,对所述阈值分割后的图像进行滤波器去噪处理。
进一步地,所述连通域检测前进行闭操作。
与现有技术相比,本发明具有以下优点:(1)本发明提出一种分割乳腺病灶的方法,无须用户手动提供种子点,简单、有效且快速,可广泛应用于各种注射造影剂的核磁共振成像数据。(2)稳定性以及鲁棒性强,不仅适用于普通的乳腺核磁共振成像,对于抑脂的乳腺图像也有效。(3)不会因为成像脂肪部分不够饱满而出现乳腺区域分割不准确的情形。(4)对分割后的乳腺区域进行病灶检测不存在模板局限性的问题。(5)可以检测到一个以上的病灶区域。
附图说明
图1所示为本发明一个实施例分割乳腺病灶方法的流程示意图;
图2所示为本发明一个实施例分割乳腺病灶方法的胸部粗分割流程示意图;
图3所示为本发明一个实施例分割乳腺病灶方法的病灶检测的流程示意图;
图4a~图4j为图1实施例所示的各阶段的结果示意图。
具体实施方式
在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施的限制。
其次,本发明利用示意图进行详细描述,在详述本发明实施例时,为便于说明,所述示意图只是实例,其在此不应限制本发明保护的范围。
为了解决背景技术中的技术问题,本实施例提供了一种分割乳腺病灶的方法。图1所示为本发明一个实施例分割乳腺病灶方法的流程示意图,请参见图1:
首先,执行步骤S11:对注射造影剂前的磁共振图像基于椭圆模型进行胸部粗分割,获得第一图像,所述胸部粗分割步骤用于去除诸如心脏等高亮部分的影像,减少图像处理过程中的数据计算量,同时减少高亮影像部分所导致的假阳性肿块,提高分割乳腺病灶的准确性。
具体地,如图2本发明实施例分割乳腺病灶方法的胸部粗分割流程示意图所示,所述第一图像的取得包括如下步骤:
首先执行步骤S21:对所述注射造影剂前的磁共振图像从上到下逐行进行检测,当检测到连通域的范围达到阈值M时,停止检测并记录行数N′;接着执行步骤S22:以行数N1=N′+n作为椭圆的起始行,以所述注射造影剂前的磁共振图像的最后一行作为椭圆的终止行为行数N2,在所述起始行N1和终止行N2之间进行椭圆分割处理;需要说明的是,两处乳房并不连通,所以乳腺部分连通域的面积小于胸部连通域的面积,可以通过设定阈值M来检测乳腺部分和胸部部分,进行所述胸部的粗分割;所述阈值M的取值范围为所述磁共振图像宽度的40~60%,本实施例中优选阈值M为50%,所述行数n的取值范围为5~20行,以防遗漏乳房与胸部连接边缘的可疑病灶,本实施例优选n为10行。
继续执行步骤S23~S27:在所述椭圆分割处理后的磁共振图像上,从上到下逐行检测,当测得灰度值为非零值时,停止检测,记录行数N′1作为乳腺部分的起始行,同时记录非零值的像素点所在的列数L;接着在所述椭圆分割处理后的磁共振图像上从下到上逐行检测第L列,当测得的灰度值为非零值时,停止检测,记录行数N′2作为乳腺部分的终止行,最后在所述起始行N′1和终止行N′2之间提取所述乳腺部分,去除所述胸部椭圆分割操作后的背景,获得所述第一图像。
需要说明的是,如图4a所示的注射造影剂前的磁共振图像进行胸部粗分割时,首先进行滤波器的去噪处理,用于去除背景噪音,获得如图4b所示去噪后的磁共振图像,所述滤波器去噪处理可以采用中值滤波器、高斯滤波器、高通滤波器、低通滤波器、均值滤波器等对图像进行滤波,本实施例优选中值滤波器,用于去除点状的礁岩噪音;接着进行闭操作,对滤波后的图像进行边缘规则化的填充操作,获得如图4c所示的闭操作后的磁共振图像,以防止图像上出现断层和空洞现象。
接着,执行步骤S12:对所述注射造影剂前、后的磁共振图像进行剪影,获得如图4d所示的剪影图像,接着根据第一图像,在所述剪影图像中提取乳腺部分,获得如图4e所示乳腺部分的剪影图像;所述注射造影剂前、后的磁共振图像进行剪影前,首先进行配准操作。
需要说明的是,本实施例中注射造影剂后的磁共振图像为如图4f所示的造影剂初始增强的磁共振图像,所述造影剂初始增强的磁共振图像为注射造影剂后大约60~70毫秒左右采集到的磁共振图像,以区别其它增强图像造成的正常组织同样变亮而引起的误诊和漏诊。
继续执行步骤S13:对每层乳腺部分的剪影图像进行病灶检测,获得第二图像,用以检测可疑的乳腺病灶并去除血管、小体积噪点等假阳性肿块,提高分割乳腺病灶的准确率。
具体地,如图3本发明一个实施例分割乳腺病灶方法的病灶检测的流程示意图所示,首先执行步骤S31~S33:对每层乳腺部分的剪影图像根据灰度直方图进行阈值分割,首先对每层乳腺部分的剪影图像计算灰度直方图,根据所述灰度直方图分别定义每层图像的高、低阈值,例如第k层的高阈值Tk1和低阈值Tk2;若第k层图像的高阈值Tk1<a×Tave,则以高阈值Tk1对该层进行阈值分割;若第k层图像的高阈值Tk1≥a×Tave,则以低阈值Tk2对该层进行阈值分割,最终获得如图4g所示阈值分割后的磁共振图像。其中,所述k取自然数,阈值Tave为各层乳腺部分剪影图的高阈值平均值,所述a为倍数,其取值范围为0.9~1.2。
具体地,取两层乳腺部分的剪影图像,根据第一层乳腺部分剪影图像的灰度直方图选取保留至少最亮的300个像素点的高阈值T11,设定T11=1000;同样地根据第二层剪影图像的灰度直方图选取保留至少最亮的300个像素点的高阈值T21,设定T21=3000,两层剪影图像的高阈值平均值Tave=1500,取倍数a=1.2,则T11<1.2×Tave,说明第一层剪影图像中比较亮的像素点数较少,含有可疑病灶的可能性较小,因此,以高阈值T11=1000对该层进行阈值分割,而T21≥1.2×Tave,说明该层包含比较亮的像素点数较多,含有病灶的可能性高,因此以该层的低阈值T22对图像进行阈值分割,以尽可能保留潜在病灶的信息,例如第二层剪影图像的低阈值T22=2000进行阈值分割后,可以获得至少最亮的3000个像素点,包含更多的潜在病灶信息,提高检测分割病灶的准确性。
接着执行步骤S34:在阈值分割后的图像上首先进行滤波器去噪处理,获得如图4h所示去噪后的阈值分割图像,接着对所示去噪后的阈值分割图像进行检测,若三维像素点数的体积小于阈值b,则去除所述三维像素点;因病灶具有一定体积,可以通过本步骤排除体积过小的噪音等干扰因素,提高分割病灶的准确性。本实施例中阈值b为经验值,依据图像情况而定,在本实施例中阈值b取10~50立方毫米。
继续执行步骤S35:首先在阈值分割后的图像上,依次得到三维连通域,依次对每个连通域进行二维投影,以防止不同层的血管与肿瘤在二维投影时发生重叠或连接,影响分割结果。接着依次对所述二维投影进行检测,当满足时,则去除所述三维连通域。其中,所述lmax为所述二维投影最长轴的长度,面积S为该连通域的二维投影面积,所述c为倍数,其取值范围为0.9~1.2。本步骤用于去除非圆形假阳性肿块,例如血管等。最后执行步骤S36:通过进行病灶检测,获得第二图像;所述第二图像为如图4i所示的病灶检测后的磁共振图像。
本步骤通过阈值分割、三维像素点的体积大小以及连通域的二维投影面积与最长轴的比较,有效地去除不同形状、大小的假阳性肿块,提高本实施例分割乳腺病灶方法的准确性。
最后执行步骤S14:在所述第二图像上依次找到三维连通域,依次以每个连通域的重心为种子点在所述剪影图像上进行自适应的区域增长,获得如图4j所示的分割病灶后的磁共振图像,以供医生作为辅助医学诊断。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。