基于瓣膜对齐及概率图的心脏CT‑TEE配准方法与流程

本发明属于图像处理技术领域，特别涉及一种对心脏ct与tee图像的配准方法，可用于对心脏解剖结构的识别与跟踪。

背景技术：

随着医学成像技术以及计算机处理技术的快速发展，医学成像的模态越来越丰富，如ct、mr、pet、spect、超声图像等。不同的成像模态之间具有很大的差异性，多模态医学图像配准通过融合不同模态的成像信息，将多种图像进行结合，在同一幅图像上显示多幅图像的成像信息，其目的是将多种多样的信息准确地融合到同一幅图像中，以便能够更精确地从各个角度观察病灶和解剖结构。

根据不同的衡量准则，图像配准有不同的分类方式，例如，根据空间维数的分类、根据变换方式的分类、根据优化算法的分类、根据算法所基于的特征的分类、根据所采用的相似性测度的分类等；根据空间维度，图像配准可分为2d-2d、2d-3d、3d-3d的配准。根据空间变换方式的不同，可以划分为刚体变换和非刚体变换。图像配准中优化算法的选择也是多样的，通常有梯度下降法、牛顿法、powell法、遗传算法等。医学图像可以提取的特征非常丰富，通常包括特征点、表面纹理、像素强度以及外部特征等。基于特征点的配准是指通过选取在几何上有特别意义的可以定位的特征点集，比如不连续点，图形的转折点，线交叉点、在医学上具有解剖意义的点等，并进行坐标匹配；基于表面的配准是指通过分割的方式提取出图像感兴趣区域的轮廓作为配准的特征空间；基于像素值的配准是指利用整幅图像的像素值或体素值构成特征空间；对于医学图像，基于外部特征的配准是指通过在受检者体内固定标记物或向体内注入显影物质以获得在图像上的确定的标记点。根据所采用的相似性测度的不同，图像配准又可以分为基于互相关的配准、基于互信息的配准等。

由collignon等提出的互信息法是近几年研究的热点之一，特别是在多模态医学图像配准领域中。在国外，最先利用互信息法进行医学图像配准的是viola等。之后又有学者在互信息的基础上提出了许多改进算法，maes等提出了归一化互信息，降低了传统的互信息作为测度函数对两幅配准图像重叠区域的敏感度。josien提出了一种互信息结合梯度的相似性度量准则gmi，成功应用到了mr、ct、pet等异源图像的配准上。fan等人将小波变换与互信息结合，成功地应用到可见光图像和红外光图像的配准上。但在心脏ct与tee图像的配准上，由于两种成像模态在像素值上差异性显著，且心脏tee图像的边缘和纹理表现模糊，这些医学图像配准方法均难以对其进行精确有效的配准。

技术实现要素：

本发明的目的在于针对ct和tee图像成像模态的巨大差异性，提出一种基于瓣膜对齐及概率图的心脏ct-tee配准方法，通过融合这两种成像模态下的结构和纹理信息，实现对心脏ct-tee图像的精确配准。

为实现上述目的，本发明包括如下步骤：

(1)分别对心脏ct图像ir(x)及tee图像if(x)进行交互式分割，得到ct图像ir(x)的分割图，即感兴趣区域gr(x)，以及tee图像if(x)的分割图，即感兴趣区域gf(x)；

(2)在进行交互式分割的同时引入ct图像ir(x)的3个心脏瓣膜端点或中点位置坐标x1r＝(i1r,j1r),x2r＝(i2r,j2r),x3r＝(i3r,j3r)，作为ct图像的3个特征点，同时引入tee图像if(x)中与ct图像ir(x)的3个特征点对应的心脏瓣膜端点或中点位置坐标x1f＝(i1f,j1f),x2f＝(i2f,j2f),x3f＝(i3f,j3f)，作为tee图像的3个特征点，将ct图像的3个特征点和tee图像的3个特征点构成3对特征点对，将这3对特征点对作为基础配准的先验信息；

(3)将ct图像ir(x)作为参考图像，将tee图像if(x)作为浮动图像，基于瓣膜先验信息，对参考图像和浮动图像进行基础配准，得到基础配准的变换矩阵t1,并得到基础配准结果s1(x)；

(4)设定ct图像ir(x)的增强矩阵vr(x)，分别对ct图像ir(x)和感兴趣区域gr(x)进行区域增强，得到增强后的ct图像irh(x)和增强后的感兴趣区域grh(x)，设定tee图像if(x)的增强矩阵vf(x)，分别对tee图像if(x)和感兴趣区域gf(x)进行区域增强，得到增强后的tee图像ifh(x)和增强后的感兴趣区域gfh(x)；

(5)基于区域增强后的ct图像irh(x)和增强后的感兴趣区域grh(x)，生成ct图像ir(x)的概率图pr(x)，基于区域增强后的tee图像ifh(x)和增强后的感兴趣区域gfh(x)，生成tee图像if(x)的概率图pf(x)；

(6)对(5)中生成的两个概率图pr(x)和pf(x)进行相似性度量，并将(3)中得到的基础配准的变换矩阵t1作为最终配准中寻优算法的初始参数，基于概率图pr(x)和pf(x)的相似性对ct图像ir(x)和tee图像if(x)进行最终配准，求得变换矩阵t2，并得到最终配准结果s2(x)。

本发明与现有技术相比具有以下优点：

1、本发明基于实际应用中对瓣膜位置信息的需求，以及tee图像中瓣膜成像清晰且ct图像中瓣膜位置可见的现实分析，在基础配准中引入了瓣膜的空间位置信息，以此作为先验信息进行特征点匹配，极大地提高了ct与tee图像配准的有效性；

2、本发明将基础配准得到的变换矩阵作为最终配准中寻优算法powell的初始参数，避免了powell算法因为初始参数选择不合适而进入局部最优的问题，提升了局部寻优算法的效率及准确性；

3、本发明基于对感兴趣部位进行区域增强生成概率图，极大的简化了概率图的生成过程，并且将本发明中生成的概率图与归一化互信息相结合，极大的提升了相似性度量的准确性和可靠性；

附图说明

图1是本发明的实现总流程图；

图2是本发明所使用的心脏ct图像，即配准的参考图像；

图3是本发明所使用的心脏tee图像，即配准的浮动图像；

图4是本发明中基于图2生成的概率图像；

图5是本发明中基于图3生成的概率图像；

图6是本发明对图2与图3进行基础配准后的配准图像；

图7是本发明对图2与图6进行融合后的融合图像；

图8是本发明对图2与图3进行最终配准后的配准图像；

图9是本发明对图2与图8进行融合后的融合图像。

具体实施方式

以下结合附图，对本发明的具体实施方案和效果作进一步的解释和说明：

参照图1，本发明基于瓣膜对齐及概率图互信息的心脏ct-tee图像配准方法

其实现步骤如下：

步骤1：对输入的心脏ct图像ir(x)和tee图像if(x)进行交互式分割。

1a)输入心脏ct图像ir(x)，如图2所示，人为地选取ct图像ir(x)中右心房以及主动脉所在的区域，并选取该区域中少量几个目标像素点以及几个背景点，进行数次迭代，得到ct图像的感兴趣区域gr(x)；

1b)输入tee图像if(x)，如图3所示，人为地选取tee图像if(x)中右心房以及主动脉所在的区域，并选取该区域中少量几个目标像素点以及几个背景点，进行数次迭代，得到tee图像的感兴趣区域gf(x)。

步骤2：获得心脏瓣膜先验信息。

2a)在交互式分割的同时，引入ct图像ir(x)的3个心脏瓣膜端点或中点位置坐标x1r＝(i1r,j1r),x2r＝(i2r,j2r),x3r＝(i3r,j3r)，作为ct图像的3个特征点，同时引入tee图像if(x)中与ct图像ir(x)的3个特征点对应的心脏瓣膜端点或中点位置坐标x1f＝(i1f,j1f),x2f＝(i2f,j2f),x3f＝(i3f,j3f)，作为tee图像的3个特征点；

2b)根据心脏瓣膜可以分为二尖瓣、三尖瓣、主动脉瓣或肺动脉瓣的特性，用ct图像的3个特征点与tee图像的3个特征点构成3对特征点对，并将这3对特征点对作为瓣膜先验信息。

步骤3：对ct图像ir(x)和tee图像if(x)进行基础配准。

3a)将ct图像ir(x)作为参考图像，将tee图像if(x)作为浮动图像；

3b)基于瓣膜先验信息，即基于参考图像的3个特征点的位置坐标x1r＝(i1r,j1r),x2r＝(i2r,j2r),x3r＝(i3r,j3r)和浮动图像的3个特征点的位置坐标x1f＝(i1f,j1f),x2f＝(i2f,j2f),x3f＝(i3f,j3f)，按如下公式求得3对特征点对的坐标变换矩阵t1：

3c)基于变换矩阵t1，对浮动图像进行仿射变换，并通过三次插值得到基础配准结果s1(x)：s1(x)＝t1×if(x)，如图6所示，将图2与图6融合，得到融合图，如图7所示，从图7可以看出，经过基础配准后，参考图像图2与浮动图像图3在物理位置上得到了大致的匹配。

步骤4：对ct图像ir(x)和其感兴趣区域gr(x)进行区域增强。

4a)根据ct图像ir(x)的像素分布设定增强矩阵vr(x)：

设定vr(x)的大小与ct图像ir(x)的大小一致，并将vr(x)中与感兴趣区域坐标一致的点均设为一个大于0的固定值，本发明中设定的是80，其余位置的点均设为0；

4b)对ct图像ir(x)和感兴趣区域gr(x)进行灰度增强，得到增强后的ct图像irh(x)和增强后的感兴趣区域grh(x)：

irh(x)＝ir(x)+vr(x)

grh(x)＝gr(x)+vr(x)。

步骤5：对tee图像if(x)和其感兴趣区域gf(x)进行区域增强。

5a)根据tee图像if(x)的像素分布设定增强矩阵vf(x)：

设定vf(x)的大小与tee图像if(x)的大小一致，并将vf(x)中与感兴趣区域坐标一致的点均设为一个大于0的固定值，本发明中设定为80，其余位置的点均设为0；

5b)对tee图像if(x)和其感兴趣区域gf(x)进行灰度增强，得到增强后的tee图像ifh(x)和增强后的感兴趣区域gfh(x)：

ifh(x)＝if(x)+vf(x)

gfh(x)＝gf(x)+vf(x)。

步骤6：生成ct图像ir(x)的概率图pr(x)。

6a)基于步骤4得到的区域增强后的ct图像irh(x)和增强后的感兴趣区域grh(x)，生成经过区域增强后的ct图像irh(x)关于增强后感兴趣区域grh(x)的概率密度函数：

该概率密度函数相当于一个概率检索表，当输入某个像素点时，输出即为该像素点可能位于感兴趣区域中的概率值；

6b)将区域增强后的ct图像irh(x)作为概率密度函数fr(i)的输入，得到ct图像ir(x)的概率图pr(x)：pr(x)＝fr(x),x∈irh(x)，如图4所示。

步骤7：生成tee图像if(x)的概率图pf(x)。

7a)基于步骤5得到的区域增强后的tee图像ifh(x)和增强后的感兴趣区域gfh(x)，生成经过区域增强后的tee图像ifh(x)关于增强后感兴趣区域gfh(x)的概率密度函数：

7b)将区域增强后的tee图像ifh(x)作为概率密度函数ff(i)的输入，得到tee图像if(x)的概率图pf(x)：pf(x)＝ff(x),x∈ifh(x)，如图5所示。

步骤8：对ct图像ir(x)和tee图像if(x)进行最终配准。

8a)对概率图pr(x)和pf(x)进行相似性度量

相似性度量准则可以是互信息、归一化互信息、互相关、梯度互信息等，本发明采用的是归一化互信息，其求解过程如下：

8a1)求出图像ir(x)的概率图pr(x)的熵hr(x)和图像if(x)的概率图pf(x)的熵

hf(x)：

其中，p(pr(x))为概率图pr(x)的概率密度函数，p(pf(x))为概率图pf(x)的概率密度函数；

8a2)求出概率图pr(x)和pf(x)的联合熵hr,f(x)：

其中，p(pr(x),pf(x))为概率图pr(x)和pf(x)的联合概率密度函数；

8a3)得到概率图pr(x)和pf(x)的归一化互信息nmpi：

8b)对概率图pr(x)和pf(x)的归一化互信息nmpi进行寻优

寻优算法的种类很多，如粒子群算法、模拟退火算法、蚁群算法、梯度下降法、powell法，本发明中采用的是powell法，具体实现过程如下：

8b1)将基础配准的变换矩阵t1作为powell法的初始参数；

8b2)通过powell法寻优得到当概率图pr(x)和pf(x)的归一化互信息nmpi最大时的坐标变换矩阵t2，通过对tee图像if(x)进行坐标变换得到最终配准结果s2(x)：s2(x)＝t2×if(x)，如图8所示，将图2与图8进行融合，得到融合图，如图9所示。

从图9可以看出来，经过最终配准后，参考图像图2与浮动图像图3得到了精准的配准效果。

以上描述仅是本发明的一个具体实例，不构成对本发明的任何限制，显然对于本领域的专业人员来说，在了解了本发明内容和原理后，都可能在不背离本发明原理、结构的情况下，进行形式和细节上的各种修正和改变，但是这些基于本发明思想的修正和改变仍在本发明的权利要求保护范围之内。