基于互信息的医学图像配准与拼接方法研究

发布时间：2020-04-20 04:09

【摘要】：随着医学成像技术的迅速发展,出现了不同的医学成像设备,如计算机断层成像(CT)、磁共振成像(MRI)、数字减影血管造影(DSA)、单光子发射断层成像(SPECT)、正电子发射断层成像(PET)等。不同的成像技术能从人体获取不同方面的信息,然而,由于成像原理的差异,每种成像设备并不能提供人体各个方面的信息。为了使多次成像或多种模式成像的信息可以得到综合利用,减少由于信息不完整等因素引起的临床诊断和治疗误差；研究人员开始寻找一项能够在一幅图像上同时表达来自人体不同信息的技术,如功能图像和解剖图像；而达到这一目的最基本的要求是使不同模式图像中的同一解剖结构达到空间位置上的一致,从而产生了医学图像配准技术。 医学图像配准是对一幅医学图像寻求一种空间变换,使它与另一幅医学图像上的对应点达到空间位置上的一致。配准过程本质上是一个多参数最优化问题,是图像融合的前提；在诊断医学、手术规划、模拟仿真以及放疗计划等方面都有重要应用。因此,对医学图像配准的研究,具有重要的学术意义和临床应用价值。 由Collignon和Maes、Viola和Wells几乎同时提出的互信息配准法不需要对被处理图像进行任何预处理,不需要对配准图像间的灰度关系进行限制；并且,互信息法几乎可以用于任何不同模态图像的配准,已广泛应用到多模态医学图像的配准中。图像配准方法包括特征空间、空间变换、空间搜索、相似度测量四个方面。互信息配准方法是将待配准图像的全部灰度信息作为特征空间,用互信息作为相似度来测量两幅图像的相关性。 在刚性配准过程中,计算两幅图像的互信息时,容易引入新的灰度值；为了避免这种问题的出现,一般采用部分体积插值统计两幅图像的联合直方图。但在该方法计算中,当图像平移整数点时,联合直方图统计会出现误差,使目标函数出现局部极值,从而造成误配准。针对这一问题,本文将高斯函数引入到直方图统计中,选取适当的邻域,用高斯函数计算邻域内各点像素对联合直方图的贡献。利用高斯函数的平滑性,避免了在互信息计算过程中统计图像联合直方图时出现误差。使用Powell优化方法,寻找最佳的优化参数,实现图像的最佳配准。采用CT-PET数据进行实验,证明该方法平滑了目标函数,有效地消除了局部极值,提高了多模态图像配准的精确性,并且,对噪音图像配准也产生很好的效果。该方法适用于多模态医学图像配准,克服了传统互信息计算时的不足,提高了配准的正确率和精确度。 目前提出的许多医学图像配准算法只是针对刚性变换,弹性配准虽然已经提出了一些方法,但同刚性配准算法相比还不成熟,不能满足临床实时需求；由于研究对象具有多样性和复杂性,尽管已提出了许多弹性配准算法,可并没有一种配准方法能在各个方面都达到临床需求。也就是说已提出的方法都有一定的局限性,如缺少实时性和有效的全自动性,这些不足在一定程度上制约着医学图像弹性配准算法在实际临床中的应用。因此医学图像弹性配准有着广泛的临床应用,也是医学图像处理领域的研究热点。 文章实现的非刚性配准算法,是基于B-样条的FFD变换模型,用互信息作为相似度测量,满足了配准计算速度和精度的要求。本文选用三次B-样条曲线做变换模型,具有良好的局部控制特性,每个控制点变换时,只影响周围4×4的邻域;控制面网格的疏密可以用来控制形变程度,当控制面网格较密时,变形模型趋向于描述局部形变；控制面网格较疏时,变形模型趋向于对全局进行形变；控制点的移动只影响其周围网格邻域内的邻近点,在计算互信息时,无需计算两幅图像整体的互信息来判断该控制点是否移动到形变的准确位置,只需计算两幅图像对应控制点周围邻域的互信息,从而达到减少计算量的目的。 图像拼接是在医学领域有着广泛的应用。随着医疗技术的发展,临床医学对x线影像有一些新的需求,在脊柱和长骨等的骨科诊断和手术治疗中,获取完整解剖结构影像、解剖径线测量和角度测量。对于术前诊断、术中监控和术后评价都具有十分重要的意义。因此,由图像拼接获取的全景图像具有重要的临床应用价值。 医学图像拼接分搜索重叠区域、图像配准、灰度映射、图像融合等几部分。文章中实验数据为医学图像采用基于区域的灰度相关方法进行重叠区域搜索,用灰度映射对图像灰度调整及融合。基于区域的拼接方法,一般直接利用待拼接图像和参考图像中像素的灰度相关性,确定相邻两幅图像的重叠区域,然后构造模板进行匹配,从而得到完整图像。图像拼接中相邻帧间灰度差较大时,找到重叠区域,直接进行图像简单的叠加会出现皮带现象,导致全景图像不连续。针对这一问题提出了灰度映射的方法,统计两幅图像重叠区域的联合直方图,计算出新的灰度映射值；将一幅图像映射到另一幅图像,从而减少灰度差,完成图像的色差调整。采用多组CR数据进行实验,该方法对数据进行色差调整,实现了图像的无缝拼接。实验证明,该方法用互信息作为相似度,提高了匹配的精确性,对医学图像拼接有一定的使用价值。
【图文】：

配准示意图

则两变量的互信息最大。图像配准过程即寻找最优的空间变换参数T0: T0=argmaxI(R，T(F))基于互信息配准算法流程图，，如图2一2。优优化化 插插值值 空空间变换换图2一2配准流程图Figure.2一 2RegistrationProeess
【学位授予单位】：南方医科大学
【学位级别】：硕士
【学位授予年份】：2011
【分类号】：R814;TP391.41

【引证文献】

中国硕士学位论文全文数据库 前4条

1 李勇;基于互信息的图像拼接算法研究[D];哈尔滨工业大学;2016年

2 冯骥;基于互信息和B样条FFD模型的图像配准算法研究[D];哈尔滨理工大学;2016年

3 袁茂恺;自吸离散式数字PCR芯片检测装置的关键技术研究与系统开发[D];浙江大学;2016年

4 王海风;机器人辅助外科视觉导航系统图像配准研究[D];吉林大学;2014年

本文编号：2634134