一种非刚体的配准方法发明专利.docxVIP

一种非刚体的配准方法 技术领域 本发明涉及医学图像处理领域，尤其涉及一种非刚体的配准方法。 背景技术 医学图像配准是医学图像三维重建与可视化的先决条件。医学图像配准 包括刚体配准和非刚体配准，刚体配准主要用于脑部或全身脊柱等只存在全 局旋转和平移变化的扫描图像；而非刚体配准用于心脏、肺或乳腺等存在局 部位移的器官扫描图像。 基于自由形变的非刚体配准是一种传统并且比较通行的非刚体配准方 法，其对应的图像坐标点的位移是其周围网格控制点位移的B样条函数，所 述B样条函数的函数关系即为配准中的坐标变换。优化过程中求解的是控制 点的位移量。优化的代价函数如果结合弯曲能量可使坐标变换平滑，保持非 刚体拓扑变换的连续可导。采用由稀疏到稠密的多级控制点密度逐步地对浮 动图像进行校正，本身即符合非刚体形变存在的不同尺度特性，可相应提高 配准的精度和速度。 具体地，在基于自由形变(FFD，free?form?deformation)的非刚体配准 方法中，当采用3阶B样条时，像素点在X、Y和Z三个方向上的位移分量是 其周围的若干个数量的控制点的位移分量的B样条函数。如图1所示为像素 点与对应B样条控制点的关系示意图。其中，图中阵列排布的点为控制点。X 为待移动的像素点位置，位于所述图像上的控制点数目为5*5*5，图像外各有 1排或2排控制点。其中，在三维方向上的位移量分别为u、v、w的表达式如 下： x?′?y?′?z?′?=?T?(?x?y?z?)?=?x?y?z?+?Σ?l?=?0?3?Σ?m?=?0?3?Σ?n?=?0?3?B?l?(?u?)?B?m?(?v?)?B?n?(?w?)?dx?i?+?l?,?j?+?m?,?k?+?n?Σ?l?=?0?3?Σ?m?=?0?3?Σ?n?=?0?3?B?l?(?u?)?B?m?(?v?)?B?n?(?w?)?dy?i?+?l?,?j?+?m?,?k?+?n?Σ?l?=?0?3?Σ?m?=?0?3?Σ?n?=?0?3?B?l?(?u?)?B?m?(?v?)?B?n?(?w?)?dz?i?+?l?,?j?+?m?,?k?+?n?]] 式中是控制点的序号，如图1所示， B0(u)=(1-u)3／6，B1(u)=(3u3-6u2+4)/6，B2(u)=(-3u3+3u2+3u+1)/6， B3(u)=u3／6. dx，dy，dz分别为控制点在X、Y和Z三个方向上的位移分量，即配准过程 的待优化参数。 为了使位移场分布的比较平滑，有时候会将一个称作弯曲能量(bending? energy)的项添加到参考图像与校正的浮动图像之间的相似性测度上，使得 优化程序收敛时得到一个变化比较平滑的位移场。 具体可参见文献1：D.Rueckert，L.I.Sonoda，C.Hayes，D.L.G.Hill，M.O. Leach，and?D.J.Hawkes，‘Nonrigid?Registration?Using?Free-FormDeformations： Application?to?Breast?MR?Images’，IEEE?TRANSACTIONS?ON?MEDICAL? IMAGING，VOL.18,NO.8,AUGUST1999和文献2：Marius?Staring，Stefan? Klein，‘Spatial?derivatives?and?penalty?terms?in?ITK?and?elastix’。 所述弯曲能量的计算公式如下： 1?V?∫?0?X?∫?0?Y?∫?0?z?(?(???2?T???x?2?)?2?+?(???2?T???y?2?)?2?+?(???2?T???z?2?)?2?+?2?(???2?T???x???y?)?2?+?2?(???2?T???x???z?)?2?+?2?(???2?T???y???z?)?2?)?dxdydz?]] ，其中式中V为全部像素数。 如图2所示，对于这种基本的FFD非刚体配准方法，为了校正较细小的 浮动图像形变，采用比较稠密的控制点时，控制点位移向量会出现相互交叉 等拓扑畸变，导致校正后的图像局部出现怪异的畸变。 美国专利US8,184,129B2利用归一化互信息测度(NMI，Normalized? Mutual?Information)和FFD对图像进行非刚体配准。具体地对图像进行了子 区域的划分，进行对应子区域单独配准，然后将位移场进行插值平滑得到全 局的位移场。但是算法实现复杂，计算较费时。美国专利申请2009／0304252A1 在配准CT图像和超声图像用胸隔膜和血管的灰度与梯度等信息生成优化的代 价函数，属于基于特征而不是基于灰度的配准方法，算法适用范围较窄。 发明内容 本发明解决的问题是提供一种非刚体的配准方法，用以避免在较稠密控 制点密度