CN105654450B - 基于局部和全局区域测地线模型的mr图像分割和偏移场矫正方法 - Google Patents
基于局部和全局区域测地线模型的mr图像分割和偏移场矫正方法 Download PDFInfo
- Publication number
- CN105654450B CN105654450B CN201410625664.7A CN201410625664A CN105654450B CN 105654450 B CN105654450 B CN 105654450B CN 201410625664 A CN201410625664 A CN 201410625664A CN 105654450 B CN105654450 B CN 105654450B
- Authority
- CN
- China
- Prior art keywords
- image
- local
- biased field
- function
- segmentation
- 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.)
- Active
Links
Abstract
本发明属MR图像分割及应用领域,涉及一种基于全局和局部区域的测地线模型进行MR图像分割和偏移场矫正的方法,本发明基于MR图像分割在医学图像分析中的至关重要性,以及MR图像灰度不均匀性使分割更为困难,噪声和偏移场导致灰度不均匀性等原因,采用基于全局和局部的符号压力函数提取图像的全局和局部信息处理灰度不均匀图像;在局部符号压力函数中加入偏移场矫正项,实现同时分割MR图像和偏移场矫正,克服偏移场造成的灰度不均匀性;将模型由二项水平集扩展到四项水平集,实现脑MR图像灰质、白质和脑脊液的准确分割。本方法应用在合成图像和MR图像中,分割结果表明其具有显著的准确性和高效性。
Description
技术领域
本发明属于图像分割(Magnetic Resonance,MR)及应用领域。涉及新的基于水平集的分割算法,具体涉及一种基于局部和全局区域的测地线模型的方法,该方法能实现同时MR图像分割和偏移场矫正。
背景技术
现有技术公开了水平集方法(level set method)是一种求解曲线演化的方法,其以一种隐含的方式来表达平面闭合曲线或立体闭合曲面,从而避免对闭合曲线演化过程中的跟踪,将曲线演化转化成一个纯粹的求解偏微分方程数值解问题。近年来,水平集方法逐渐成为研究的热点,并被应用于图像分割、图像平滑、运动分割即运动目标跟踪,甚至立体视觉及图像的修复。与传统的分割方法比较,水平集方法优势明显,它能够在统一的数学模型中约束图像数据、初始形状、先验知识、统计分析等信息,这种特性使它被广泛地应用于复杂的医学图像的分割。
在实际应用中,由于偏移场的存在,MR图像存在灰度不均匀性。图像灰度不均匀性现象为图像分割带来了严峻的挑战,也阻碍了医疗领域计算机辅助系统的发展。为了克服MR图像灰度不均匀现象的不利影响,国内外学者先后提出了多种基于纠正偏移场思想的解决方案,在一定程度上缓解了灰度不均匀对分割的影响。截至目前为止,MR图像分割方法很多,但分割精度仍然不是很高。
因此,本申请的发明人拟提供一种提高MR图像分割的准确性和高效性的新的基于水平集的分割算法,使能实现同时MR图像分割和偏移场矫正。
于本发明相关的现有技术有如下参考文献:
[1]Kass M,Witkin A,Terzopoulos D.Snakes:Active contour models[J].International journal of computer vision,1988,1(4):321-331.
[2]Caselles V,Kimmel R,Sapiro G.Geodesic active contours[J].International journal of computer vision,1997,22(1):61-79.
[3]Goldenberg R,Kimmel R,Rivlin E,et al.Fast geodesic active contours[J].Image Processing,IEEE Transactions on,2001,10(10):1467-1475.
[4]Chan T F,Vese L A.Active contours without edges[J].Imageprocessing,IEEE transactions on,2001,10(2):266-277.
[5]Li C,Kao C Y,Gore J C,et al.Implicit active contours driven bylocal binary fitting energy[C]//Computer Vision and Pattern Recognition,2007.CVPR'07.IEEE Conference on.IEEE,2007:1-7.
[6]Li C,Kao C Y,Gore J C,et al.Minimization of region-scalablefitting energy for image segmentation[J].Image Processing,IEEE Transactionson,2008,17(10):1940-1949.
[7]Li C,Huang R,Ding Z,et al.A level set method for imagesegmentation in the presence of intensity inhomogeneities with application toMRI[J].Image Processing,IEEE Transactions on,2011,20(7):2007-2016.
[8]Zhang K,Zhang L,Song H,et al.Active contours with selective localor global segmentation:a new formulation and level set method[J].Image andVision computing,2010,28(4):668-676.
[9]Wang L,Li C,Sun Q,et al.Active contours driven by local and globalintensity fitting energy with application to brain MR image segmentation[J].Computerized Medical Imaging and Graphics,2009,33(7):520-531.
[10]Aubert G,Kornprobst P.Mathematical problems in image processing:partial differential equations and the calculus of variations[M].Springer,2006.
[11]Xu C,Yezzi Jr A,Prince J L.On the relationship between parametricand geometric active contours[C]//Signals,Systems and Computers,2000.Conference Record of the Thirty-Fourth Asilomar Conference on.IEEE,2000,1:483-489.
[12]Lankton S,Tannenbaum A.Localizing region-based active contours[J].Image Processing,IEEE Transactions on,2008,17(11):2029-2039.
[13]Guillemaud R,Brady M.Estimating the bias field of MR images[J].Medical Imaging,IEEE Transactions on,1997,16(3):238-251.
[14]Vovk U,Pernus F,Likar B.A review of methods for correction ofintensity inhomogeneity in MRI[J].Medical Imaging,IEEE Transactions on,2007,26(3):405-421.
[15]Vese L A,Chan T F.A multiphase level set framework for imagesegmentation using the Mumford and Shah model[J].Internationaljournal ofcomputer vision,2002,50(3):271-293.
[16]Chung G,Dinov I D,Toga A W,et al.MRI tissue segmentation using avariational multilayer approach[M]//Computational Biomechanics forMedicine.Springer New York,2010:5-16.。
发明内容
本发明的目的在于克服现有技术的缺陷,提出了一种新的基于局部和全局区域的测地线模型的方法,该方法能实现同时的MR图像分割和偏移场矫正,高效精确的将脑MR图像分割为灰质、白质和脑脊液。
本发明基于MR图像分割在医学图像分析中的至关重要性,以及MR图像灰度不均匀性使分割更为困难,噪声和偏移场导致灰度不均匀性等原因,提出了一种新的基于局部和全局区域的测地线模型的方法,其中包括:采用基于全局和局部的符号压力函数提取图像的全局和局部信息处理灰度不均匀图像;在局部符号压力函数中加入偏移场矫正项,实现同时分割MR图像和偏移场矫正,克服偏 移场造成的灰度不均匀性;将模型由二项水平集扩展到四项水平集,实现脑MR图像灰质、白质和脑脊液的准确分割。
具体的,本发明主要针对由于偏移场而导致灰度不均匀的MR图像,进一步扩展到四项水平集,实现脑MR图像灰质、白质和脑脊液的准确分割;其中包括步骤:
首先,提出了一个新的符号压力函数,包括代表局部图像信息的局部符号压力函数和代表全局图像信息的全局符号压力函数;所述的局部符号压力函数可以使轮廓准确地停在目标边界,限制图像的灰度不均匀性;所述的全局符号压力函数可以使轮廓快速的靠近目标,因此,使初始轮廓曲线选择更加灵活和算法更加高效;
其次,对MR图像进行偏移场矫正,由于偏移场在局部范围内缓慢变化的,本发明中使用局部符号压力函数计算偏移场,水平集演化过程中不断地矫正偏移场,达到同时MR图像分割和偏移场矫正的目的;
最后,为了分割多个目标,将其中的由二项水平集扩展到四项,成功的分割出脑MR图像中的灰质、白质和脑脊液。
更具体的,本发明的基于局部和全局区域测地线模型的MR图像分割和偏移场矫正方法,其特征在于,其包括:
假设要处理的图像为I(x),x∈Ω,初始轮廓为C;通过下述步骤:
1)定义新的符号压力函数;
如图1所示,初始轮廓(绿圆)将图像分为前景和背景两个区域,同时它将局部轮廓
(黑圆)分为局部的内部(红色)区域和外部(绿色)区域,通过计算得到前景和背景的平均灰
度值C1和C2,局部拟合值f1(x)和f2(x),
本发明提出了一种新的符号压力函数来代替测地线模型中的边缘停止函数其中,Kσ是高斯核函数,具体见公式φ(x)是 水平集函数,H(φ(x))是Heaviside函数,具体见公式ε是一个常量,一般取ε=1;
SL为局部符号压力函数,
SG为全局符号压力函数,
新的符号压力函数为:
sNew(I(x))=(1-ω)·sL(I(x))+ω·sG(I(x)),x∈Ω, (3)
ω为高斯权重
ω=exp(-β(Mmax-Mmin)2), (4)
其中Mmax和Mmin为局部小窗口内的最大和最小灰度值,由此可见,在灰度变化大的地方ω的值越大,全局符号压力函数SG占的比重越大,快速将轮廓趋于目标边界。反之,在灰度变化小的地方ω的值越小,局部符号压力函数SL占的比重越大,准确的将轮廓停在目标边界。
将新的符号压力函数替代测地线模型中的边缘停止函数,最小化能量函数得到水平集演化方程为,
2)偏移场矫正:
偏移场是造成MR图像灰度不均匀的主要原因,其表现为脑MR图像上同一类别的像素灰度沿空间呈缓慢的变化,含偏移场的MR图像可建模为:I(x)=b(x)·J(x)+n(x)。式中,I(x)为观测图像,J(x)为真实图像,b(x)为偏移场,n(x)是噪声。实际上,真实图像J(x)在每个组织中可以用值c(x)来代替,并且在一个小的窗口内偏移场是缓慢变化的,可以看做是一个常量,因此,fi(x)≈b(x)·c(x),
将偏移场带入局部符号压力函数,
令M1(φ(x))=H(φ(x)),M2(φ(x))=1-H(φ(x)),ci和b为,
3)扩展到四项水平集来实现MR图像灰质、白质和脑脊液的分割:
如图2所示,初始轮廓(蓝圆和红圆)代表两个水平集函数,同时它将局部(黑圆)区域分为局部内部(红色)区域和局部外部(绿色)区域,两个水平集φ1和φ2将图像分为四个区域{φ1>0,φ2>0},{φ1>0,φ2<0},{φ1<0,φ2>0},{φ1<0,φ2<0}。让M1=H(φ1)H(φ2),M2=H(φ1)(1-H(φ2)),M3=(1-H(φ1))H(φ2),M4=(1-H(φ1))(1-H(φ2)),通过计算得到四个区域的平均灰度值C1、C2、C3、C4,局部拟合值f1(x)、f1(x)、f3(x)、f4(x),
定义整体能量函数为
其中局部压力函数定义为:
全局符号压力函数定义为:
b,c1...,c4,C1...,C4定义为:
最小化的能量函数(公式9)得到关于φ1和φ2的梯度下降流
其中,δ(φ)是Dirac函数,具体见公式ε是一个常量,其值为0.5。
本发明经合成图像和MR图像实验验证,结果表明,本方法的分割结果明显提高了准确性和高效性。
本发明方法具有如下优点:
(1)初始活动轮廓选择灵活,分割效率高;
(2)结果准确性非常好,基本上都类似于不同专家手工测量的结果;
(3)稳定性高,有可重复性。
(4)提高病灶检测的准确性和实用性,降低成本,辅助诊断。
附图说明
图1.二项水平集说明图。
图2四项水平集说明图。
图3MR图像分割示例,其中,(a)列是原始图像和初始轮廓,(b)列是分割结果,(c)列是偏移场,(d)列是偏移场矫正后的图像。
图4偏移场矫正前后的直方图,其中,(a)和(b)分别表示上述图3第一行原始图像和偏移场矫正后图像的直方图,(c)和(d)分别表示上述图3第二行原始图像和偏移场矫正后图像的直方图。
图5脑MR图像分割结果,其中,(a)是原始图像,(b)是Li’s方法分割结果,(c)本发明分割结果。
图6脑MR图像灰质、白质和脑脊液的分割结果显示,其中,(a)专家分割结果,(b)Li’s分割结果,(c)本发明分割结果。
具体实施方式
实施例1
本实施例采用基于局部和全局区域测地线模型的MR图像分割和偏移场矫正方法进行了试验验证,选取脑MR图像的轴状面、矢状面和冠状面各20张进行MR图像分割和偏移场矫正并统计学分析,将人工分割图像和本发明方法分割图像的重合率作为精度信息,并将Li’s方法和本发明方法进行比较,结果如表1所述,其中分别显示了白质和灰质(Tissue白质和灰质)在轴状面Axial、矢状面Sagittal和冠状面Coronal的分割精度(均值+方差),结果表明,Li‘s方法和本发明方法均得到了较高的分割结果,但本发明方法精度更高,稳定性更好。表1是对20张脑MR图像分割,并与专家分割结果比较,Li’s方法和本发明方法精度显示(均值+方差)。
表1
Claims (4)
1.基于局部和全局区域测地线模型的MR图像分割和偏移场矫正方法,其特征在于,主要针对由于偏移场而导致灰度不均匀的MR图像,进一步扩展到四项水平集,实现脑MR图像灰质、白质和脑脊液的准确分割;其包括步骤:
首先,提出一个新的符号压力函数,包括代表局部图像信息的局部符号压力函数和代表全局图像信息的全局符号压力函数;
其次,对MR图像进行偏移场矫正,其中使用局部符号压力函数计算偏移场,水平集演化过程中不断地矫正偏移场,同时对MR图像分割和偏移场矫正;
最后,将方法中的二项水平集扩展到四项,分割多个目标。
2.按权利要求1所述的方法,其特征在于,所述的方法中采用三个全局符号压力函数和三个局部符号压力函数,通过水平集演化,准确的分割多个目标。
3.按权利要求1所述的方法,其特征在于,其通过下述步骤同时对MR图像分割和偏移场矫正:
假设要处理的图像为I(x),x∈Ω,初始轮廓为C;按下述步骤:
1)定义新的符号压力函数;
所述初始轮廓中,将图像分为前景和背景两个区域,同时将局部轮廓分为局部的内部区域和外部区域,通过计算得到前景和背景的平均灰度值C1和C2,局部拟合值f1(x)和f2(x),
其中,采用新的符号压力函数代替测地线模型中的边缘停止函数
其中,Kσ是高斯核函数;
公式φ(x)是水平集函数,H(φ(x))是Heaviside函数,
公式取ε=1;
SL为局部符号压力函数,
SG为全局符号压力函数,
新的符号压力函数为:
sNew(I(x))=(1-ω)·sL(I(x))+ω·sG(I(x)),x∈Ω, (3)
ω为高斯权重
ω=exp(-β(Mmax-Mmin)2), (4)
其中Mmax和Mmin为局部小窗口内的最大和最小灰度值;
将新的符号压力函数替代测地线模型中的边缘停止函数,最小化能量函数得到水平集演化方程为,
2)偏移场矫正:
含偏移场的MR图像建模为:I(x)=b(x)·T(x)+n(x),式中,I(x)为观测图像,T(x)为真实图像,b(x)为偏移场,n(x)是噪声;实际上,真实图像T(x)在每个组织中用值c(x)代替,并且在一个小窗口内偏移场缓慢变化,看做是一个常量;因此,fi(x)≈b(x)·c(x);
将偏移场带入局部符号压力函数,
令M1(φ(x))=H(φ(x)),M2(φ(x))=1-H(φ(x)),ci和b为,
3)扩展到四项水平集来实现MR图像灰质、白质和脑脊液的分割:
初始轮廓代表两个水平集函数,同时它将局部区域分为局部内部区域和局部外部区域,两个水平集φ1和φ2将图像分为四个区域{φ1>0,φ2>0},{φ1>0,φ2<0},{φ1<0,φ2>0},{φ1<0,φ2<0};让M1=H(φ1)H(φ2),M2=H(φ1)(1-H(φ2)),M3=(1-H(φ1))H(φ2),M4=(1-H(φ1))(1-H(φ2)),通过计算得到四个区域的平均灰度值C1、C2、C3、C4,局部拟合值f1(x)、f2(x)、f3(x)、f4(x),
定义整体能量函数为
其中局部压力函数定义为:
全局符号压力函数定义为:
b,c1...,c4,C1...,C4定义为:
最小化能量函数公式(9)得到关于φ1和φ2的梯度下降流
其中,δ(φ)是Dirac函数,具体见公式ε是一个常量,其值为0.5。
4.按权利要求1所述的方法,其特征在于,所述的分割多个目标为分割出脑MR图像中的灰质、白质和脑脊液。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410625664.7A CN105654450B (zh) | 2014-11-09 | 2014-11-09 | 基于局部和全局区域测地线模型的mr图像分割和偏移场矫正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410625664.7A CN105654450B (zh) | 2014-11-09 | 2014-11-09 | 基于局部和全局区域测地线模型的mr图像分割和偏移场矫正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105654450A CN105654450A (zh) | 2016-06-08 |
CN105654450B true CN105654450B (zh) | 2019-01-15 |
Family
ID=56482832
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410625664.7A Active CN105654450B (zh) | 2014-11-09 | 2014-11-09 | 基于局部和全局区域测地线模型的mr图像分割和偏移场矫正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105654450B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106570837B (zh) * | 2016-10-28 | 2019-12-17 | 中国人民解放军第三军医大学 | 基于高斯多尺度空间的脑组织mri图像偏移场校正方法 |
CN106981096B (zh) * | 2017-02-13 | 2019-12-06 | 上海大学 | 一种vr环境中基于卷积曲面逼近的建模方法 |
CN107798684B (zh) * | 2017-11-07 | 2020-02-28 | 河南师范大学 | 一种基于符号压力函数的活动轮廓图像分割方法及装置 |
CN110047046B (zh) * | 2019-04-15 | 2022-09-06 | 南京邮电大学 | 图像分割和偏移场矫正方法、可读存储介质和终端 |
CN110390667B (zh) * | 2019-06-18 | 2023-10-20 | 平安科技(深圳)有限公司 | 基于眼底oct影像的病灶提取方法、装置、设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101714253A (zh) * | 2009-12-04 | 2010-05-26 | 西安电子科技大学 | 基于测地线活动区域模型的交互式图像分割修正方法 |
CN102063723A (zh) * | 2011-01-28 | 2011-05-18 | 南京理工大学 | 高通量成像系统下的斑马鱼间脑和中脑自动分割方法 |
CN102306373A (zh) * | 2011-08-17 | 2012-01-04 | 深圳市旭东数字医学影像技术有限公司 | 一种腹部器官的三维医学图像分割方法及系统 |
CN102737378A (zh) * | 2012-05-24 | 2012-10-17 | 中国科学院深圳先进技术研究院 | 测地线活动轮廓分割图像的方法和系统 |
CN103729843A (zh) * | 2013-12-20 | 2014-04-16 | 电子科技大学 | 基于马尔可夫的医学图像分割方法 |
CN103996193A (zh) * | 2014-05-16 | 2014-08-20 | 南京信息工程大学 | 结合加权邻域信息与偏移场恢复的脑mr图像分割方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8498481B2 (en) * | 2010-05-07 | 2013-07-30 | Microsoft Corporation | Image segmentation using star-convexity constraints |
AU2011265383A1 (en) * | 2011-12-20 | 2013-07-04 | Canon Kabushiki Kaisha | Geodesic superpixel segmentation |
-
2014
- 2014-11-09 CN CN201410625664.7A patent/CN105654450B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101714253A (zh) * | 2009-12-04 | 2010-05-26 | 西安电子科技大学 | 基于测地线活动区域模型的交互式图像分割修正方法 |
CN102063723A (zh) * | 2011-01-28 | 2011-05-18 | 南京理工大学 | 高通量成像系统下的斑马鱼间脑和中脑自动分割方法 |
CN102306373A (zh) * | 2011-08-17 | 2012-01-04 | 深圳市旭东数字医学影像技术有限公司 | 一种腹部器官的三维医学图像分割方法及系统 |
CN102737378A (zh) * | 2012-05-24 | 2012-10-17 | 中国科学院深圳先进技术研究院 | 测地线活动轮廓分割图像的方法和系统 |
CN103729843A (zh) * | 2013-12-20 | 2014-04-16 | 电子科技大学 | 基于马尔可夫的医学图像分割方法 |
CN103996193A (zh) * | 2014-05-16 | 2014-08-20 | 南京信息工程大学 | 结合加权邻域信息与偏移场恢复的脑mr图像分割方法 |
Non-Patent Citations (1)
Title |
---|
一种基于测地线模型图像分割的隐式算法;罗宏文等;《吉林大学学报(理学版)》;20070930;第45卷(第5期);819-821 |
Also Published As
Publication number | Publication date |
---|---|
CN105654450A (zh) | 2016-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ji et al. | Active contours driven by local likelihood image fitting energy for image segmentation | |
CN105654450B (zh) | 基于局部和全局区域测地线模型的mr图像分割和偏移场矫正方法 | |
Liu et al. | Distance regularized two level sets for segmentation of left and right ventricles from cine-MRI | |
Grosgeorge et al. | Graph cut segmentation with a statistical shape model in cardiac MRI | |
Peters et al. | Automatic whole heart segmentation in static magnetic resonance image volumes | |
Klodt et al. | A convex framework for image segmentation with moment constraints | |
Saritha et al. | A comprehensive review: Segmentation of MRI images—brain tumor | |
CN104268873A (zh) | 基于核磁共振图像的乳腺肿瘤分割方法 | |
Zhang et al. | A variational level set model for multiscale image segmentation | |
Min et al. | A novel dual minimization based level set method for image segmentation | |
Chen et al. | A generalized asymmetric dual-front model for active contours and image segmentation | |
CN109285176B (zh) | 一种基于正则化图割的大脑组织分割方法 | |
Pei et al. | Improved brain tumor growth prediction and segmentation in longitudinal brain MRI | |
Kama et al. | Segmentation of soft tissues and tumors from biomedical images using optimized k-means clustering via level set formulation | |
Min et al. | A polynomial piecewise constant approximation method based on dual constraint relaxation for segmenting images with intensity inhomogeneity | |
Platero et al. | A new label fusion method using graph cuts: application to hippocampus segmentation | |
Quan et al. | Segmentation of tumor ultrasound image via region-based Ncut method | |
Cho et al. | Level-set segmentation of brain tumors using a new hybrid speed function | |
Saikumar et al. | Improved fuzzy c-means clustering algorithm using watershed transform on level set method for image segmentation | |
Li et al. | Novel model using kernel function and local intensity information for noise image segmentation | |
Huang et al. | MR image segmentation based on fuzzy c-means clustering and the level set method | |
Saikumar et al. | MRI brain image segmentation algorithm using watershed transform and kernel fuzzy C-means clustering on level set method | |
Chang et al. | Subvoxel segmentation and representation of brain cortex using fuzzy clustering and gradient vector diffusion | |
Yang et al. | Joint prior models of neighboring objects for 3D image segmentation | |
Cao et al. | A variational level set SAR image segmentation approach based on statistical 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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |