一种改进心脏核磁共振图像多步配准方法.docVIP

一种改进的心脏核磁共振图像多步配准方法 　　摘要：提高核磁共振成像速度的研究方兴未艾。基于k-t的成像技术比如k-t SENSE和k-t PCA，利用成像序列在时间维度上的冗余性，能够将采样加速比提高到8倍以上。由于它们利用成像序列的时间冗余性，对成像对象的运动特别敏感，不能移除运动导致的伪影，从而造成重构图像在空间域的模糊。为了使基于k-t的成像技术能够运用于具有较大运动场合的成像对象，研究了一种改进的图像配准方法，通过多步配准来完成仿射变换。实验结果表明，该方法比传统的仿射变换配准精度更高。 　　关键词：核磁共振成像；图像配准；仿射变换 　　DOIDOI：10.11907/rjdk.161091 　　中图分类号：TP317.4 　　文献标识码：A 文章编号：1672-7800（2016）005-0207-03 　　0 引言 　　核磁共振成像（Magnetic Resonance Imaging，MRI） 是一种傅立叶成像技术，与传统的基于X射线的成像技术（如计算机断层扫描，Computed tomography ，CT）不同，MRI利用射频脉冲作用下水分子中的氢元素在磁场中产生共振的现象来成像[1]。其特点是多参数和高对比度成像，并且可以在任意方位断层，提供丰富的诊断信息。由于人体生理和物理局限性，在同一时刻只能采集到图像的一个傅立叶编码系数，所以成像速度较慢。 　　并行成像（Parallel Imaging） [2-3]是已经获得临床应用的快速成像方法，PI并不采集所有重构需要的k-space数据，而是根据预先设计的采样模式，采集整个k-space 的一部分，从而减少采样时间。比如SENSE、GRAPPA和SPIRiT，利用k-space数据的空间冗余特性，使用多个感应线圈以减少成像需要采集的k-space值。由于PI以每帧独立的方式重构，所以PI对成像对象的运动不敏感，加速率一般只有2-3倍[4]。基于k-t的成像技术利用成像序列在时间维度上的冗余性，进一步提高数据采集的加速比。k-t SENSE和k-t PCA利用低分辨率图像训练时间基函数，并使用这些时间基函数调整重构[5-6]；压缩传感（Compressed Sensing， CS），比如Sparse MRI和k-t SLR，利用重构图像的变换稀疏性调整重构，这些方法能够在8倍以上采样加速比下移除图像伪影[7-8]。由于是利用成像序列的时间冗余，所以对成像对象的运动特别敏感，不能移除运动导致的伪影，造成重构图像模糊。 　　图像配准技术能够将成像对象中的关注区域在空间上对齐[9]，从而提高重构图像的变换稀疏性，使得基于k-t 的成像技术能够运用于具有较大运动场合的成像对象，提高了MRI技术的适应性。然而对于一般配准问题而言，模型的灵活性与结果的精确性是相互矛盾的。模型选取越灵活，自由度越大，结果的精度就越低。本文采用多步配准的方法来提高配准的精确度。首先进行灵活性较差的粗模型配准，得到的解大致接近真实解。然后通过限制粗模型中的参数，进行灵活性较大的精配准，得到更加接近真实的解。仿真数据和真实数据的实验结果表明，这种方法的精确度比传统的一步配准精确度更高。 　　1 基本原理与算法流程 　　1.1 基本原理 　　配准过程通常涉及两幅图像，一副图像为参考图像，另一幅图像为浮动图像。配准的目的是根据给定的相似性度量求得一个变换，通过这个变换，可以将浮动图像变换后与参考图像在空间上对齐。配准算法常常将图像转换到物理空间内执行，而通常被表述为一个优化过程，即通过搜索求得最优的空间变换参数值，可用下式表达： 　　在配准过程中，变换模型选取越灵活，自由度会越大，但却导致结果精确度不高。心脏运动可以近似看成仿射运动，其自由度较大，若采用传统的仿射变换模型来处理，配准精确度就不高。为了保证模型的灵活性并提高配准的精度，本文对仿射变换模型进行了改进，采用多步配准方法完成配准过程，即先进行平移变换，再进行刚体变换，最后进行仿射变换。从理论上分析，先经过平移变换模型配准，可以使平移分量的参数大致接近真实解。再用平移变换模型的配准结果进行刚体变换，由于平移分量与真实解接近，在搜索最优解的过程中可以限制平移分量的取值，使解的搜索范围缩小，得到的配准结果比直接进行刚体变换结果要精确。最后在刚体变换配准的基础上进行灵活性更大的仿射变换配准，既可以保证模型的灵活性，又可以提高配准精度。对于仿射运动配准可采用如下公式表示： 　　其中，φt、φr、φa分别表示平移变换、刚体变换、仿射变换，S表示经过多步变换后的图像。在具体配准时，将平移变换后的结果作为刚体变换的输入，然后将刚体变换后的结果作为仿射变换的输入，最后完成仿射运动纠正。 　　1.2 算法流程 　　配准的目的是使