CN103337074A - 一种基于主动轮廓模型分割乳腺dce-mri病灶的方法 - Google Patents
一种基于主动轮廓模型分割乳腺dce-mri病灶的方法 Download PDFInfo
- Publication number
- CN103337074A CN103337074A CN2013102476179A CN201310247617A CN103337074A CN 103337074 A CN103337074 A CN 103337074A CN 2013102476179 A CN2013102476179 A CN 2013102476179A CN 201310247617 A CN201310247617 A CN 201310247617A CN 103337074 A CN103337074 A CN 103337074A
- Authority
- CN
- China
- Prior art keywords
- focus
- image
- interest
- mri
- mammary gland
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 210000005075 mammary gland Anatomy 0.000 title claims abstract description 39
- 238000013535 dynamic contrast enhanced MRI Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000011218 segmentation Effects 0.000 claims abstract description 32
- 238000002595 magnetic resonance imaging Methods 0.000 claims abstract description 24
- 238000009826 distribution Methods 0.000 claims abstract description 18
- 230000003902 lesion Effects 0.000 claims description 24
- 238000012800 visualization Methods 0.000 claims description 11
- 238000005520 cutting process Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000000638 solvent extraction Methods 0.000 claims description 5
- 230000001105 regulatory effect Effects 0.000 claims description 4
- 238000009877 rendering Methods 0.000 claims description 4
- 230000002708 enhancing effect Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 230000002452 interceptive effect Effects 0.000 claims description 3
- 238000005728 strengthening Methods 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 2
- 230000008676 import Effects 0.000 claims description 2
- 230000010429 evolutionary process Effects 0.000 abstract description 3
- 238000012545 processing Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000013500 data storage Methods 0.000 abstract description 2
- 238000003709 image segmentation Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 21
- 238000013461 design Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 208000026310 Breast neoplasm Diseases 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 206010006187 Breast cancer Diseases 0.000 description 4
- 239000002872 contrast media Substances 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 4
- 206010028980 Neoplasm Diseases 0.000 description 3
- 210000000481 breast Anatomy 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000001575 pathological effect Effects 0.000 description 3
- 230000000505 pernicious effect Effects 0.000 description 3
- 238000004195 computer-aided diagnosis Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 230000003211 malignant effect Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- LGMLJQFQKXPRGA-VPVMAENOSA-K gadopentetate dimeglumine Chemical compound [Gd+3].CNC[C@H](O)[C@@H](O)[C@H](O)[C@H](O)CO.CNC[C@H](O)[C@@H](O)[C@H](O)[C@H](O)CO.OC(=O)CN(CC([O-])=O)CCN(CC([O-])=O)CCN(CC(O)=O)CC([O-])=O LGMLJQFQKXPRGA-VPVMAENOSA-K 0.000 description 1
- 230000001435 haemodynamic effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000002601 intratumoral effect Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 210000004088 microvessel Anatomy 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000000819 phase cycle Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000004614 tumor growth Effects 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 210000005166 vasculature Anatomy 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
一种基于主动轮廓模型分割乳腺DCE-MRI病灶的方法,属于医学图像分割领域,包括由MRI扫描设备获取乳腺DCE-MRI影像序列数据,人工选取感兴趣区域、自动获取减影后感兴趣体积、主动轮廓分割病灶和可视化显示病灶。本发明充分利用乳腺DCE-MRI影像背景统计分布一致病灶内部分布多样这一特点设计主动轮廓模型的边缘停止函数,实现了病灶的可靠分割,有效地避免了边界泄漏现象。模型演化过程中无需重新初始化符号距离函数,使***的实时性更高。本发明的实施对人工操作要求低、智能化程度高、数据存储空间的需求量低、处理速度快,并且三维角度的分割可以有效获得病灶全面立体的空间信息,方便医生多角度观察分析病灶。
Description
技术领域
本发明涉及医学影像分割技术领域,特别涉及一种基于水平集理论的快速可靠分割乳腺动态对比增强磁共振影像序列中三维病灶的主动轮廓模型方法。
背景技术
近年来,乳腺DCE-MRI扫描已经成为一种检测乳腺疾病最具潜力的影像技术之一,在早期乳腺癌的临床诊断和医学研究中发挥着重要的作用。同时,基于DCE-MRI的计算机辅助诊断(computer aided diagnosis,CAD)***因其能够检出病灶继而生成有意义的诊断建议,能够辅助放射科医师提高诊断效率减轻诊断工作量,已被逐渐应用到临床实践中。精确可靠的病灶分割不仅标志着病灶的检出,而且制约其后形态学特征、纹理特征和血流动力学特征提取的精度进而影响诊断结论的正确性,因此病灶分割是CAD的基础环节。临床中,医师通常依据经验手工分割病灶,该过程繁琐、存在主观性、可重复性差,且不易实现三维分割也无法形成病灶完整的空间立体表示,所以实现高效的乳腺DCE-MRI三维病灶分割功能以提高CAD***的可操作性及诊断性能势在必行。
由于主动轮廓模型对具有噪声和边界间隙的图像具有较强的鲁棒性,可得到光滑连续的目标边界,且边界表达能够达到亚像素级精度等特点被广泛应用于解决乳腺DCE-MRI病灶分割问题。其中,Jiazheng Shi等人在他们乳腺肿瘤治疗反馈估计的研究中,将模糊C聚类和水平集方法结合来分割病灶,模糊C聚类获得演化初始轮廓并作为一项外力来控制水平集演化。L.Bahreini等人利用梯度矢量流方法分割二维乳腺病灶,梯度矢量流是一种参数主动轮廓模型,分割过程需调节大量参数,因此若应用于临床,将对人工理论知识要求较高。Shannon C.Agner等人提出了一种混合主动轮廓模型方法分割二维病灶,该算法结合流形学方法将图像空间转换到谱空间以改善病灶在原始空间中边缘模糊的问题,由于计算谱空间中的张量梯度耗时长,会限制其处理复杂影像的实时性。此外,多数基于主动轮廓模型的方法中为了保证算法数值计算的稳定性,曲线或曲面演化过程中需要重复进行符号距离函数初始化工作,耗时长,这很大程度上限制了算法的实时性应用。
总结目前动态对比增强乳腺磁共振影像病灶的分割研究,可发现存在以下几个难点:(1)病灶与正常组织毗邻,且病灶边界存在模糊,若不能有效区分,分割容易引起边界泄漏现象;(2)病灶内部灰度分布多样,这是进行病灶分期及诊断的有效信息,不具有一致的统计特征,常规建模往往不能涵盖病灶的全面情况;(3)病灶形状多样,无法利用基于先验形状的方法对其进行分割;(4)为获得病灶全面立体的空间信息,三维角度分割十分必要,然而DCE-MRI影像数据量大,不合理的分割过程会造成分割运算的时间长效率低。因此在考虑病灶多样性分布并有效建模正常组织基础上,进行病灶三维快速自动分割***的研究是非常必要的。
鉴于上述原因,本发明统计分析了影像中病灶周围正常组织的灰度分布和病灶内部的灰度分布,发现:病灶周围正常组织灰度分布具有一致性而病灶内部灰度分布多样。这一DCE-MRI影像表征的统计分析结果与病理学解释相吻合:乳腺癌DCE-MRI扫描的病理学基础研究指出,肿瘤的生长依赖于肿瘤内血管的生成,肿瘤内脉管***往往比正常组织更具有渗透性,当对比剂进入静脉后,会因为不同性质的肿瘤间血管的生成、微血管密度、肿瘤内毛细血管的通透性和组织细胞外间隙的不同而在影像中反映出不同的灰度分布,而正常组织则不会具有这种多样的变化。在上述病理学研究的有力支持下,我们发明了一种基于乳腺DCE-MRI背景信息的主动轮廓模型方法分割乳腺三维病灶。
发明内容
本发明提供了一种基于主动轮廓模型分割乳腺DCE-MRI病灶的方法,所要解决的技术问题是,提供一种对操作者要求低、分割病灶性能可靠、基于水平集理论能够快速识别乳腺DCE-MRI病灶边界的主动轮廓模型方法。该方法只需操作者在病灶增强明显的断层影像中进行简单的框选操作,之后便可由***自动快速分割出病灶,在二维断层影像上标记病灶和三维可视化显示病灶。
本发明的技术方案是:由MRI扫描设备获取乳腺DCE-MRI影像序列数据,影像预处理工作站对影像序列进行筛选,MRI影像存储服务器保存所有被检查者的影像序列,由MRI影像后处理工作站从影像存储服务器获取待处理的影像序列进行分析和计算得到病灶分割结果并显示在MRI影像后处理工作站的计算机显示屏上;在MRI影像后处理工作站中设置有人工选取感兴趣区域模块、自动获取减影后感兴趣体积模块、主动轮廓模型分割病灶模块和可视化显示病灶模块。
人工选取感兴趣区域模块,用于接收从MRI影像存储服务器传来的待分割乳腺DCE-MRI影像序列,操作者自由浏览该序列选择病灶增强明显的一帧断层影像框选包含病灶的感兴趣区域,并将所选断层影像和区域位置传送给自动获取减影后感兴趣体积模块。
自动获取减影后感兴趣体积模块,利用记录的人工选取的断层影像和感兴趣区域位置自动获取减影后的感兴趣体积,将感兴趣体积作为分割算法实际处理的体数据传送给主动轮廓模型分割病灶模块。
主动轮廓模型分割病灶模块,采用本发明提供的基于影像背景信息利用水平集理论设计的主动轮廓模型分割来自自动获取减影后感兴趣体积模块的体数据,并将分割结果传送给可视化显示病灶模块。
可视化显示病灶模块,用于将从主动轮廓模型分割模块获得的分割结果进行二维标记和三维可视化显示,并将结果显示在计算机屏幕上供放射科医师分析或存储到硬盘,供其它影像分析软件进行进一步分析。
具体方法如下:
步骤一、在人工选取感兴趣区域模块中,人工选取感兴趣区域模块接收从MRI影像存储服务器传来的待分割乳腺DCE-MRI影像序列,在可以自由浏览选择乳腺DCE-MRI影像序列的交互式环境下操作者在增强后的影像序列中选择病灶增强明显的一帧断层影像,利用可以调节大小的矩形窗选取包含目标病灶的感兴趣区域,由计算机自动记录所选断层影像在整个影像序列中的位置以及感兴趣区域在所选断层影像中的位置。
步骤二、在自动获取减影后感兴趣体积模块中,利用步骤一中记录的人工选取的断层影像位置,结合待分割乳腺DCE-MRI影像序列中每个时相的扫描层数自动确定在哪些断层影像序列上进行感兴趣区域的提取,然后,在所选的断层影像序列中,依据记录的感兴趣区域位置,自动获得感兴趣区域序列,获得感兴趣区域序列后,利用各时相增强后的感兴趣区域序列减增强前序列再取平均值得到减影后的感兴趣体积,即分割算法实际处理的体数据。
步骤三、在主动轮廓模型分割病灶模块中,对步骤二获得的减影后的感兴趣体积利用Gamma分布建模影像背景分布,采用水平集理论实现本发明提出的主动轮廓模型方法进行病灶分割,分割算法如下:
(1)本发明对初始轮廓的要求较低,采用二值水平集函数实现。初始轮廓面的位置可以根据减影后的感兴趣体积情况自适应调整,初始轮廓面形状可以任意选取。初始水平集函数为:
(2)考虑到影像背景的实际分布和乳腺癌DCE-MRI扫描的病理学基础,以及两个参数的分布如Gamma分布比一个参数的分布如Rayleigh分布或Poisson分布对来自不同设备的变化较大的影像数据有更好的适应性。本发明提出利用Gamma分布建模感兴趣体积背景的灰度分布。设计新的边缘停止阈值th的计算公式如下:
th=F-1(p|κ,θ)={th:F(th|κ,θ)=p} (2)
其中是Gamma累积分布函数,θ和κ分别是Gamma分布的尺度参数和形状参数,利用最大似然法结合Choi等人的理论估计Gamma分布的参数,得到的估计公式如下:
其中,I(x,y,z)表示减影后感兴趣体积函数,Ωbackground表示感兴趣体积的背景区域,A为Ωbackground区域内的体素个数。
(3)本发明提出的边缘停止函数bdspf(I)的设置公式如下:
bdspf(I)又可称为基于影像背景分布的符号压力函数(backgrounddistribution-based signed pressure function),该函数利用影像区域信息满足归一化的要求,并且符号可正可负。
(4)本发明提出的曲面演化函数是在测地线活动轮廓模型(geometric activecontour,GAC)框架下重新设计的,并结合KaihuaZhang等人的理论,通过将新的边缘停止函数取代GAC水平集演化式中的边缘停止函数形成新的演化式如下:
其中,φ为水平集函数,α为一常数用以加快演化速度。
本发明提出的主动轮廓模型算法可采用C语言或其它计算机编程语言实现,具体可按照如下步骤进行:
1)向计算机输入待分割体数据,水平集函数φ初始化成二值函数φ0如公式(1)所示,c通常设置为1;
2)设置迭代阈值ε即φt+1-φt<ε,ε取很小的数,设置高斯核Gσ的参数,置初始迭代次数t=0;
3)设置累积分布概率p,p的范围一般取0.95—1;
4)利用公式(3)估计Gamma分布的尺度参数θt和利用公式(4)估计Gamma分布的形状参数κt;
5)利用公式(2),计算边缘停止阈值tht;
6)利用公式(5)计算边缘停止函数bdspft(I);
7)根据公式(6),更新水平集函数;
8)判断是否满足收敛条件,如果满足则停止迭代,输出分割结果,否则转向步骤9),若演化曲面越界则转向步骤3)重新设置累积分布概率;
9)重设水平集函数φt,如果φt>0,φt=c,否则φt=-c;
10)使用高斯核Gσ光滑水平集函数φt,即Gσ*φt转向步骤4)。
步骤四、在可视化显示病灶模块中,利用步骤三得到的病灶分割结果在原始断层影像上标记病灶,利用表面绘制算法生成病灶的三维立体图像,并将二维和三维结果显示在计算机屏幕上供放射科医师分析、或存储到硬盘上供其它影像分析软件进行进一步分析。
本发明的有益效果是:
本发明充分考虑了病灶多样性分布并对正常组织进行有效建模,较好地解决了复杂、大数据量的乳腺DCE-MRI影像序列中三维病灶的分割问题。采用人工选取感兴趣区域的操作可以排除具有明显亮度变化的心脏等组织器官对后续病灶分割过程的影响,极大地降低计算机处理的数据量,加快计算机处理的速度。各时相影像序列间的减影操作可以简化影像,减弱有明显对比剂变化的大血管等结构的影响,使影像背景简单病灶目标明确,从而降低病灶分割的复杂度。本发明提供的主动轮廓模型采用水平集理论实现,水平集演化函数简单,轮廓演化过程中无需重复初始化符号距离函数降低了算法实现的难度,减少耗时,增加了算法的实时性。利用影像的背景信息设计驱动力可以排除由于病灶多样、内部分布不一致等因素造成的分割误差,进而使得本发明能够准确识别病灶模糊边界,避免边界泄漏现象的产生,分割可靠精度高。本发明的***框架完整、分割过程合理,对所进行的每一步处理都充分考虑了计算的复杂性和对数据存储空间的需求量,使分割快速稳定。本发明要求的人工操作简单、智能化程度高。
附图说明
图1是本发明从获取影像数据到输出分割结果的流程图
图2是本发明乳腺三维病灶分割***模块的流程图
图3是本发明主动轮廓模型算法实施的程序框图
图4是本发明用不同初始轮廓分割乳腺DCE-MRI病灶的分割效果图
图5是本发明用于乳腺DCE-MRI断层影像序列的分割效果图
图6是本发明用于乳腺DCE-MRI三维病灶分割过程中表面绘制效果图
具体实施方式
下面结合实例和附图对本发明的实现方法做进一步详细说明。
为实施本发明提出的乳腺DCE-MRI三维病灶分割方法需要按照如图1所示的过程,从MRI扫描设备获取被检查者的乳腺MRI影像序列数据经影像预处理工作站筛选,将筛选后的数据保存到用来保存所有被检查者影像序列的MRI影像存储服务器中,再由MRI影像后处理工作站从MRI影像存储服务器获取待处理的影像序列进行分析和计算。MRI影像后处理工作站中的乳腺三维病灶分割***构建包括人工选取感兴趣区域模块、自动获取减影后感兴趣体积模块、主动轮廓模型分割病灶模块和可视化显示病灶模块。这些模块的工作需要按照如图2所示的步骤进行:
第1步,接收从MRI影像存储服务器中传来的某位被检查者的待分割的乳腺DCE-MRI影像序列,根据实际情况,获取的乳腺DCE-MRI三维影像序列的扫描时相数一般为7-9个时相,每个时相可能包含76-184幅断层影像。也就是说至少会有500张左右的断层影像需要处理。
第2步,考虑到待处理影像巨大的数据量,以及存在大量的冗余信息,并且为了简化操作,降低处理对象的复杂度,采用人工选取感兴趣区域的方式。在交互式操作环境下,操作者浏览增强后的影像序列,选择病灶动态增强明显的断层影像,在此断层影像中框选出包含病灶的感兴趣区域。矩形窗大小可根据具体情况调节。在实例分割中我们选择80×80像素(该尺寸远大于病灶的大小)固定大小的矩形窗。计算机自动记录该断层影像在整个序列中的位置和人工选取的感兴趣区域在断层影像中的位置,如被选取的感兴趣区域位于第4增强时相第100帧影像从左上角坐标(301,301)到右下角坐标(380,380),在此采用图像坐标。
第3步,考虑到早期乳腺癌病灶比较小、MRI存在容积效应,使得病灶在同一时相仅可见于几幅断层影像,例如图5第二行良性病灶仅可见于6帧断层影像。计算机根据每个时相扫描的层数以及人工选取感兴趣区域影像所在的位置自动确定对哪些断层影像进行分割。计算机根据记录感兴趣区域位置,自动计算出所选断层影像序列所有该位置的感兴趣区域序列。对获得感兴趣序列进行减影操作。所采用的方式是将各时相增强影像序列取平均,再减去增强前影像序列,得到减影序列,即包含病灶的感兴趣体积。
第4步,利用主动轮廓模型算法分割包含病灶的感兴趣体积提取病灶,程序流程如图3所示。
第5步,在断层影像序列中标记病灶,可视化显示三维病灶,以供医师进行分析,或将最终分割的结果存储到硬盘中,供其它影像分析软件进行进一步处理。
本发明采用1.5T或3.0T的MRI扫描设备获取的DCE-MRI影像序列数据实验。被检查者采取俯卧位,双侧乳腺充分暴露,自然垂入扫描床中的专用线圈内。采用乳腺容积扫描成像序列进行轴位扫描,3.0T设备上的重复时间(TR)为4.0ms,回波时间(TE)为2.1ms,1.5T设备上的TR为4.2ms,TE为2.0ms,层扫无间隔,帧影像矩阵大小为512×512,层厚为2mm或2.2mm。对比剂采用钆喷酸葡胺(Gd-DTPA)。对比剂注入57秒后开始不间断增强扫描。扫描层数从76-184不等。为了泛化本发明中影像特征研究的意义,不限制DCE-MRI的时相数,如:病例数据可以是7时相、8时相或9时相的影像序列。因此,一个病例的总影像数据量为扫描层数与时相序列数的乘积。
人工选取感兴趣区域模块、自动获取减影后感兴趣体积模块、主动轮廓模型分割病灶模块和可视化显示病灶模块所在的MRI影像后处理工作站要求的最低配置为:硬件:CPU:intel i3530,主频2.93G,双核心四线程,DMI总线2.5GT/s,二级缓存2×256K,三级缓存4M,倍频22倍,外频133M,主板:英特尔H55芯片组,内存:2G,DDR3,1333MHz,硬盘:500GB以上,5400-7200转/分,SATA接口,数据传输率300M/秒,显卡:核心频率750MHz,GDDR3显存512M,位宽128bit以上,显示器:19.1英寸,分辨率1440×900,软件:Windows XP或Windows7操作***。
图4为本发明用不同初始轮廓分割乳腺DCE-MRI病灶的分割效果图,图中分别利用了不同初始轮廓来分割良恶性病灶。取p=0.98,ε=10-8,初始轮廓位于目标影像的不同位置,经过若干次迭代后均收敛到病灶边缘。前三列为恶性病灶,后三列为良性病灶。由图可见,采用本发明提出的主动轮廓模型进行分割时,初始轮廓的选择更加灵活。
图5是本发明用于乳腺DCE-MRI断层影像序列的分割效果图,取p=0.98,ε=10-8。其中,第一行表示恶性病灶,第二行表示良性病灶,由效果图可知,恶性病灶边缘模糊,形状不规则,良性病灶边缘光滑,形状规则。图6是本发明用于乳腺DCE-MRI三维病灶分割过程中表面绘制效果图。其中,第一列表示初始轮廓面,第二列表示中间迭代分割结果,第三列表示最终的分割结果。图5第一行与图6第一行对应是恶性病灶,图5第二行与图6第二行对应是良性病灶。由图5和图6可见,本发明提供的主动轮廓模型可以很好地分割出乳腺三维病灶。
Claims (1)
1.一种基于主动轮廓模型分割乳腺DCE-MRI病灶的方法,由MRI扫描设备获取乳腺DCE-MRI影像序列数据,影像预处理工作站对影像序列进行筛选,MRI影像存储服务器保存所有被检查者的影像序列,由MRI影像后处理工作站从影像存储服务器获取待处理的影像序列进行分析和计算得到病灶分割结果并显示在MRI影像后处理工作站的计算机显示屏上;其特征在于,在MRI影像后处理工作站中设置有人工选取感兴趣区域模块、自动获取减影后感兴趣体积模块、主动轮廓模型分割病灶模块和可视化显示病灶模块,具体方法如下:
步骤一、在人工选取感兴趣区域模块中,人工选取感兴趣区域模块接收从MRI影像存储服务器传来的待分割乳腺DCE-MRI影像序列,在可以自由浏览选择乳腺DCE-MRI影像序列的交互式环境下操作者在增强后的影像序列中选择病灶增强明显的一帧断层影像,利用可以调节大小的矩形窗选取包含目标病灶的感兴趣区域,由***自动记录所选断层影像在整个影像序列中的位置以及感兴趣区域在所选断层影像中的位置;
步骤二、在自动获取减影后感兴趣体积模块中,利用步骤一中记录的人工选取的断层影像位置,结合待分割乳腺DCE-MRI影像序列中每个时相的扫描层数自动确定在哪些断层影像序列上进行感兴趣区域的提取,然后,在所选的断层影像序列中,依据记录的感兴趣区域位置,自动获得感兴趣区域序列,获得感兴趣区域序列后,利用各时相增强后的感兴趣区域序列减增强前序列再取平均值得到减影后的感兴趣体积,即分割算法实际处理的体数据;
步骤三、在主动轮廓模型分割病灶模块中,对步骤二获得的减影后的感兴趣体积利用Gamma分布建模影像背景分布,采用水平集理论实现主动轮廓模型方法进行病灶分割,分割算法的流程如下:
1)向计算机输入待分割数据,水平集函数φ初始化成二值函数φ0如下所示:
2)设置迭代阈值ε即φt+1-φt<ε,ε取很小的数,高斯核Gσ的参数,置初始迭代次数t=0;
3)设置累积分布概率p,p的范围一般取0.95—1;
4)利用)估计Gamma分布的尺度参数θt和利用估)A计2lGoagmm}a分布的形状参数κt,其中,I(x,y,z)表示减影后感兴趣体积函数,Ωbackground表示感兴趣体积的背景区域,A为Ωbackground区域内的体素个数;
8)判断是否满足收敛条件,如果满足则停止迭代,输出分割结果,否则转向步骤9),若演化曲面越界则转向步骤3)重新设置累积分布概率;
9)重设水平集函数φt,如果φt>0,φt=c,否则φt=-c;
10)使用高斯核Gσ光滑水平集函数φt,即Gσ*φt转向步骤4);
步骤四、在可视化显示病灶模块中,利用步骤三得到的病灶分割结果在原始断层影像上标记病灶,利用表面绘制算法生成病灶的三维立体图像,并将二维和三维结果显示在计算机屏幕上供放射科医师分析、或存储到硬盘上供其它影像分析软件进行进一步分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310247617.9A CN103337074B (zh) | 2013-06-18 | 2013-06-18 | 一种基于主动轮廓模型分割乳腺dce-mri病灶的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310247617.9A CN103337074B (zh) | 2013-06-18 | 2013-06-18 | 一种基于主动轮廓模型分割乳腺dce-mri病灶的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103337074A true CN103337074A (zh) | 2013-10-02 |
CN103337074B CN103337074B (zh) | 2016-01-13 |
Family
ID=49245222
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310247617.9A Expired - Fee Related CN103337074B (zh) | 2013-06-18 | 2013-06-18 | 一种基于主动轮廓模型分割乳腺dce-mri病灶的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103337074B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103886564A (zh) * | 2013-11-29 | 2014-06-25 | 沈阳东软医疗***有限公司 | 一种pet心脏图像中心肌轮廓线分割方法和装置 |
CN104657984A (zh) * | 2015-01-28 | 2015-05-27 | 复旦大学 | 三维超声乳腺全容积图像感兴趣区域的自动提取方法 |
CN104867151A (zh) * | 2015-05-26 | 2015-08-26 | 沈阳市第四人民医院 | 一种基于延展性主动轮廓模型的眼底图像血管分割方法 |
CN104933288A (zh) * | 2014-03-18 | 2015-09-23 | 三星电子株式会社 | 用于对医学图像中的解剖元素进行可视化的设备和方法 |
CN105892959A (zh) * | 2016-04-28 | 2016-08-24 | 中山大学孙逸仙纪念医院 | 基于***mri影像数据定制个性化硅胶假体模具的方法 |
CN106780720A (zh) * | 2016-11-30 | 2017-05-31 | 上海联影医疗科技有限公司 | 医学图像显示方法及装置 |
CN107657622A (zh) * | 2017-09-05 | 2018-02-02 | 广州医科大学附属肿瘤医院 | 放射性核素断层影像与放射性浓集水平关系的分析方法 |
CN107748900A (zh) * | 2017-11-08 | 2018-03-02 | 山东财经大学 | 基于区分性卷积神经网络的乳腺肿瘤分类方法和装置 |
CN109492653A (zh) * | 2018-11-15 | 2019-03-19 | 深圳市比邻星精密技术有限公司 | 乳腺病灶体积测量方法、装置、计算机设备及存储介质 |
CN110009645A (zh) * | 2019-04-11 | 2019-07-12 | 东北大学 | 一种肝癌病灶影像的双层轮廓分割方法 |
CN112802040A (zh) * | 2021-01-28 | 2021-05-14 | 上海藤核智能科技有限公司 | 一种基于边缘感知的x光片气胸分割及评估方法 |
CN113034426A (zh) * | 2019-12-25 | 2021-06-25 | 飞依诺科技(苏州)有限公司 | 超声图像病灶描述方法、装置、计算机设备和存储介质 |
CN113112476A (zh) * | 2021-04-14 | 2021-07-13 | 中国人民解放军北部战区总医院 | 一种识别海马硬化所致颞叶癫痫的致痫灶和/或预测其病理分型的方法和*** |
CN113298770A (zh) * | 2021-05-20 | 2021-08-24 | 武汉工程大学 | 图像噪声水平估计方法、装置及计算机存储介质 |
CN114820591A (zh) * | 2022-06-06 | 2022-07-29 | 北京医准智能科技有限公司 | 图像处理方法、装置、电子设备及介质 |
CN116051553A (zh) * | 2023-03-30 | 2023-05-02 | 天津医科大学朱宪彝纪念医院(天津医科大学代谢病医院、天津代谢病防治中心) | 一种在三维医学模型内部进行标记的方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101334895A (zh) * | 2008-08-07 | 2008-12-31 | 清华大学 | 一种针对动态增强乳腺磁共振影像序列的影像分割方法 |
CN102663755A (zh) * | 2012-04-18 | 2012-09-12 | 北京理工大学 | 一种针对灰度不均匀的核磁共振图像的分割方法 |
CN102682449A (zh) * | 2012-04-25 | 2012-09-19 | 中国人民解放军军事医学科学院卫生装备研究所 | 软组织核磁图像自适应外力水平集自动分割及实现方法 |
US20130071000A1 (en) * | 2010-05-17 | 2013-03-21 | Tatsuya Takagi | Radiographic-image processing device |
-
2013
- 2013-06-18 CN CN201310247617.9A patent/CN103337074B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101334895A (zh) * | 2008-08-07 | 2008-12-31 | 清华大学 | 一种针对动态增强乳腺磁共振影像序列的影像分割方法 |
US20130071000A1 (en) * | 2010-05-17 | 2013-03-21 | Tatsuya Takagi | Radiographic-image processing device |
CN102663755A (zh) * | 2012-04-18 | 2012-09-12 | 北京理工大学 | 一种针对灰度不均匀的核磁共振图像的分割方法 |
CN102682449A (zh) * | 2012-04-25 | 2012-09-19 | 中国人民解放军军事医学科学院卫生装备研究所 | 软组织核磁图像自适应外力水平集自动分割及实现方法 |
Non-Patent Citations (2)
Title |
---|
CHUNMING LI 等: "A Level Set Method for Image Segmentation in the Presence of Intensity Inhomogeneities With Application to MRI", 《IMAGE PROCESSING, IEEE TRANSACTIONS ON》 * |
李珂 等: "乳腺计算机辅助诊断中DCE-MRI图像特征的选择与分析", 《北京生物医学工程》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103886564A (zh) * | 2013-11-29 | 2014-06-25 | 沈阳东软医疗***有限公司 | 一种pet心脏图像中心肌轮廓线分割方法和装置 |
CN104933288B (zh) * | 2014-03-18 | 2019-08-16 | 三星电子株式会社 | 用于对医学图像中的解剖元素进行可视化的设备和方法 |
CN104933288A (zh) * | 2014-03-18 | 2015-09-23 | 三星电子株式会社 | 用于对医学图像中的解剖元素进行可视化的设备和方法 |
US10383602B2 (en) | 2014-03-18 | 2019-08-20 | Samsung Electronics Co., Ltd. | Apparatus and method for visualizing anatomical elements in a medical image |
CN104657984A (zh) * | 2015-01-28 | 2015-05-27 | 复旦大学 | 三维超声乳腺全容积图像感兴趣区域的自动提取方法 |
CN104657984B (zh) * | 2015-01-28 | 2018-10-16 | 复旦大学 | 三维超声乳腺全容积图像感兴趣区域的自动提取方法 |
CN104867151A (zh) * | 2015-05-26 | 2015-08-26 | 沈阳市第四人民医院 | 一种基于延展性主动轮廓模型的眼底图像血管分割方法 |
CN104867151B (zh) * | 2015-05-26 | 2018-05-04 | 沈阳市第四人民医院 | 一种基于延展性主动轮廓模型的眼底图像血管分割方法 |
CN105892959A (zh) * | 2016-04-28 | 2016-08-24 | 中山大学孙逸仙纪念医院 | 基于***mri影像数据定制个性化硅胶假体模具的方法 |
CN106780720A (zh) * | 2016-11-30 | 2017-05-31 | 上海联影医疗科技有限公司 | 医学图像显示方法及装置 |
CN106780720B (zh) * | 2016-11-30 | 2020-06-16 | 上海联影医疗科技有限公司 | 医学图像显示方法及装置 |
CN107657622A (zh) * | 2017-09-05 | 2018-02-02 | 广州医科大学附属肿瘤医院 | 放射性核素断层影像与放射性浓集水平关系的分析方法 |
CN107657622B (zh) * | 2017-09-05 | 2021-02-05 | 广州医科大学附属肿瘤医院 | 放射性核素断层影像与放射性浓集水平关系的分析方法 |
CN107748900A (zh) * | 2017-11-08 | 2018-03-02 | 山东财经大学 | 基于区分性卷积神经网络的乳腺肿瘤分类方法和装置 |
CN109492653A (zh) * | 2018-11-15 | 2019-03-19 | 深圳市比邻星精密技术有限公司 | 乳腺病灶体积测量方法、装置、计算机设备及存储介质 |
CN109492653B (zh) * | 2018-11-15 | 2022-04-12 | 深圳市比邻星精密技术有限公司 | 乳腺病灶体积测量方法、装置、计算机设备及存储介质 |
CN110009645B (zh) * | 2019-04-11 | 2023-06-23 | 东北大学 | 一种肝癌病灶影像的双层轮廓分割方法 |
CN110009645A (zh) * | 2019-04-11 | 2019-07-12 | 东北大学 | 一种肝癌病灶影像的双层轮廓分割方法 |
CN113034426A (zh) * | 2019-12-25 | 2021-06-25 | 飞依诺科技(苏州)有限公司 | 超声图像病灶描述方法、装置、计算机设备和存储介质 |
CN113034426B (zh) * | 2019-12-25 | 2024-03-08 | 飞依诺科技股份有限公司 | 超声图像病灶描述方法、装置、计算机设备和存储介质 |
CN112802040A (zh) * | 2021-01-28 | 2021-05-14 | 上海藤核智能科技有限公司 | 一种基于边缘感知的x光片气胸分割及评估方法 |
CN112802040B (zh) * | 2021-01-28 | 2024-05-31 | 上海藤核智能科技有限公司 | 一种基于边缘感知的x光片气胸分割及评估方法 |
CN113112476A (zh) * | 2021-04-14 | 2021-07-13 | 中国人民解放军北部战区总医院 | 一种识别海马硬化所致颞叶癫痫的致痫灶和/或预测其病理分型的方法和*** |
CN113112476B (zh) * | 2021-04-14 | 2023-08-29 | 中国人民解放军北部战区总医院 | 一种识别致痫灶和/或预测其病理分型的方法和*** |
CN113298770A (zh) * | 2021-05-20 | 2021-08-24 | 武汉工程大学 | 图像噪声水平估计方法、装置及计算机存储介质 |
CN114820591A (zh) * | 2022-06-06 | 2022-07-29 | 北京医准智能科技有限公司 | 图像处理方法、装置、电子设备及介质 |
CN116051553A (zh) * | 2023-03-30 | 2023-05-02 | 天津医科大学朱宪彝纪念医院(天津医科大学代谢病医院、天津代谢病防治中心) | 一种在三维医学模型内部进行标记的方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103337074B (zh) | 2016-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103337074B (zh) | 一种基于主动轮廓模型分割乳腺dce-mri病灶的方法 | |
US11983799B2 (en) | System and method for synthesizing low-dimensional image data from high-dimensional image data using an object grid enhancement | |
Zhang et al. | Review of breast cancer pathologigcal image processing | |
US8155405B2 (en) | System and method for lesion segmentation in whole body magnetic resonance images | |
US9792703B2 (en) | Generating a synthetic two-dimensional mammogram | |
EP2916738B1 (en) | Lung, lobe, and fissure imaging systems and methods | |
US8401264B2 (en) | Solid modeling based on volumetric scans | |
CN103236058B (zh) | 获取四维心脏图像感兴趣体积的方法 | |
CN111563897B (zh) | 基于弱监督学习的乳腺核磁影像肿瘤分割的方法及装置 | |
US9940706B2 (en) | Image processing apparatus, image processing method and storage medium for extracting a lesion having a ground glass opactiy | |
CN107622492A (zh) | 肺裂分割方法及*** | |
Kalra | Developing fe human models from medical images | |
CN103501699A (zh) | 用于隔离成像数据中的潜在异常的方法和装置和其对医学图像的应用 | |
CN109191510B (zh) | 一种病理切片的3d重建方法及其装置 | |
Xiang et al. | Skeleton cuts—An efficient segmentation method for volume rendering | |
CN114972362A (zh) | 一种基于RMAU-Net网络的医学图像自动分割方法与*** | |
Bu et al. | Lung nodule detection based on YOLOv3 deep learning with limited datasets | |
CN103530902A (zh) | 医学影像二维处理及三维重建*** | |
Chen et al. | Sketch-based Volumetric Seeded Region Growing. | |
CN104331864A (zh) | 基于非下采样轮廓波和视觉显著模型的乳腺影像处理 | |
Hopp et al. | Breast tissue characterization by sound speed: correlation with mammograms using a 2D/3D image registration | |
CN110264428A (zh) | 一种基于3d卷积反卷积和生成对抗网络的医学图像降噪方法 | |
CN110533667A (zh) | 基于图像金字塔融合的肺部肿瘤ct影像3d分割方法 | |
Wu et al. | Automatic segmentation of ultrasound tomography image | |
Garcia et al. | Multimodal breast parenchymal patterns correlation using a patient-specific biomechanical model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160113 |
|
CF01 | Termination of patent right due to non-payment of annual fee |