论文部分内容阅读
核磁共振成像(Magnetic Resonance Imaging, MRI)技术已经成为脑疾病临床诊断的重要辅助手段。准确地分割脑MR图像对脑部解剖、脑疾病诸如阿尔茨海默氏病、帕金森氏病以及精神分裂症的分析与研究等具有重要的指导意义。然而在实际应用中,由于射频场的不均匀性等因素,导致脑MR图像的灰度均匀性变差,其表现为同一组织的像素灰度沿空间呈缓慢平滑的变化,直接导致传统的基于全局灰度值的分割模型分割失败;其次,脑MR图像成像过程中由于受仪器设备等物理原因影响,使得图像中经常含有噪声,影响分割精度。除此之外,图像中存在Partial Volume(PV)效应,即各个软组织之间边界比较模糊,不明确和不连续。这些都给脑MR图像的准确分割造成了很大困难。和成人脑MR图像相比,新生儿脑MR图像分辨率更低,PV效应更加显著;由于新生儿的大脑处在迅速发育中,导致图像中的灰度不均匀性的成因更加复杂;白质/灰质对比度的逆转,使得处在脑脊液和灰质交界处的像素具有和白质相似的灰度,使得分割更加困难。所有这些使得脑MR图像的分割仍然是一个充满挑战的课题。大部分脑MR图像分割方法是基于单个像素的,是离散的,容易受到初始化的影响。为此,本文对成人脑MR图像和新生儿脑MR图像分割进行了深入研究,首先采用连续的方法,提出了一些基于活动轮廓的分割模型来克服脑MR图像中存在的灰度不均匀性;然后,进一步将这些分割模型统一到一个凸优化框架中,同时分割脑MR图像和估计偏移场,特别地,由于采用了凸优化技术,可以很好地克服初始化的影响;最后我们进一步将该凸优化框架模型和活动轮廓模型应用到新生儿脑MR图像分割中,取得了较好的效果。所做的主要工作和研究成果如下:(1)为了克服灰度不均匀性,提出了一种基于图像局部均值的分割模型。利用多相位水平集方式来拟合局部灰度,并采用一种简单而有效的初始化方法,达到更快更准地分割脑MR图像。由于采用了多相位水平集方法,该模型可以同时得到白质、灰质和脑脊液。另外水平集方法很好地保证了分割结果的光滑性,可以得到光滑的边界/曲面。(2)为了克服对初始化的敏感性,提出了一种基于局部均值和全局均值的活动轮廓模型。该模型分别定义了一个局部拟合力和全局拟合力。其中,局部拟合力用于驱动曲线向目标靠近并在目标边界处停止,全局拟合力用于驱动远离目标边界的曲线向目标靠拢。在这两个拟合力的作用下,该模型对初始曲线较不敏感同时克服脑MR图像中存在的灰度不均匀性。(3)提出了一种基于局部高斯概率的活动轮廓分割模型。通过分析每个像素邻域的灰度分布,利用窗口函数定义了一个局部高斯拟合能量。局部高斯概率中的两个参数,局部均值和局部方差是随着空间变化而变化的,因此本文方法可以很好地克服脑图像中存在的灰度不均匀性。值得一提的是,本文方法的局部均值和局部方差不再是显式定义的,而是可以通过变分法原理严格推导出来的。3T和7T的脑MR图像分割实验证明本文方法不仅可以有效地克服灰度不均匀性,还可以克服低对比度。(4)提出了一种新颖的基于凸优化的脑MR图像分割模型,同时进行图像分割和估计偏移场,因此可以克服灰度不均匀性的影响。特别地,由于该模型关于图像分割变量是凸的,因此对初始化具有很强的鲁棒性。另外该模型也可以克服噪声的影响。和传统的基于梯度下降法求解相比,我们采用Split Bregman方法对该模型进行快速求解,使得对一幅2D的图像处理时间一般不超过1秒,3D的数据处理时间一般不超过100秒。大量的二维和三维实验证明本文方法无论在时间效率上还是准确度上都取得了满意的结果。而且,我们还进一步给出一个基于凸优化的脑MR图像分割框架模型,该框架综合利用了图像的局部统计信息。比较实验也证明了本文提出的框架模型是有效的。(5)提出基于凸优化和耦合水平集方法的新生儿脑MR图像分割框架。耦合水平集方法综合利用了图像的局部灰度统计信息、基于Atlas的先验信息、大脑皮层的厚度信息。为了克服初始化对分割结果的影响,我们采用凸模型进行分割,然后进行必要的校正,从而为耦合水平集方法提供一个可靠的初始化。该新生儿脑MR图像分割框架已经在10个新生儿数据上取得了满意的结果,和目前流行的新生儿MR图像分割方法相比,本文方法取得了较高的分割精度。