基于空间和阻抗的carto电解剖图与ct图像的配准方法和装置的制造方法
基于空间和阻抗的carto电解剖图与ct图像的配准方法和装置的制造方法
【技术领域】
[0001] 本发明涉及手术导航中图像配准技术领域，尤其涉及心房颤动（房颤）射频消融 介入手术导航中的CART0电解剖图与CT图像的配准方法和装置。
【背景技术】
[0002] 房颤是临床上最常见的心律失常，有着较高的致残率和致死率。三维标测和消融 系统CART0是治疗房颤的首选手术平台。
[0003] CART0电解剖图是消融手术中实时模拟生成的心房解剖结构图，依据医生在心房 内壁标测的若干个标测点构建。医生可以直接在CART0电解剖图上操作导管实施消融，然 而，由于模拟构建CART0电解剖图的标测点较少，导致电解剖图分辨率低，肺静脉的数目、 具体开口部分及每支肺静脉各自的分支情况等等这些复杂的解剖结构无法在电解剖图上 清晰地显示出来，医生难以准确定位并有效隔离肺静脉，致使房颤复发率高，并可能伴随肺 静脉狭窄等严重并发症，从而极大地影响消融治疗效果。
[0004] 众所周知，CT图像分辨率高，能清晰地显示上述消融靶区的复杂结构，然而，CT图 像只能在术前获得，医生不能在上面直接实施导管操作。
[0005] 配准CART0电解剖图和CT图像，可帮助实时并清楚地显示心腔和肺静脉的解剖结 构。依据两者的配准图像，医生将导管头固定在发生房颤的肺静脉前庭，在上面设定围绕左 或右侧上、下肺静脉的环状消融线；导管头释放射频电能，逐点消融，把导致房颤的纤维组 织烧掉，直至肺静脉出现与左心房无关的自律性电位。因此，CART0电解剖图和CT图像的 配准图像能够让医生从多个角度观察消融靶区的结构，从而增加消融位点的精确性，减少 手术相关并发症，缩短手术时间和X线曝光时间。
[0006] 然而，现有的CART0电解剖图和CT图像配准技术达到的配准精度均难以满足消融 导管精确导航的要求。
【发明内容】
[0007] 有鉴于此，本发明提供了一种房颤射频消融介入手术导航中的CART0电解剖图与 CT图像的配准方法和装置，以提高CART0电解剖图与CT图像的配准精度，为房颤消融导管 提供精确导航。
[0008] 需要说明的是，CART0电生理解剖图包括有CART0电解剖图、CART0阻抗图等； CART0阻抗图是在CART0电解剖图上叠加阻抗信息生成的，其实质是叠加了阻抗信息的 CART0电解剖图，阻抗图与电解剖图有着完全相同的空间信息。因此，配准了CART0阻抗图 与CT图像，也就是配准了CART0电解剖图与CT图像。
[0009] 为了达到上述发明目的，本发明的第一方面提出了一种CART0电解剖图与CT图像 的配准方法，包括如下步骤：
[0010] 获取CT图像和CART0电生理解剖图，所述CART0电生理解剖图包括CART0电解剖 图和CARTO阻抗图；所述CT图像是对CT原始图像进行重建分割之后获得的心房内壁的曲 面图像；所述CART0阻抗图是在CART0电解剖图上叠加阻抗信息生成的，有着与CART0电解 剖图完全相同的空间信息，同时，CART0阻抗图中的阻抗值与所述CT图像的CT值相关；
[0011] 配准所述CART0电解剖图与CT图像，得到粗配准变换；
[0012] 以所述粗配准变换作为初始变换，配准所述CART0阻抗图与CT图像，得到精配准 变换；
[0013] 依据所述精配准变换对所述CART0电解剖图作变换，将变换后的CART0电解剖图 叠加在所述CT图像上，得到最终配准的CART0电解剖图与CT图像。
[0014] 基于本发明的第一方面提出的CART0电解剖图与CT图像的配准方法，本发明同时 还提出了一种CART0电解剖图与CT图像的配准装置，包括：
[0015] 获取图像数据单元，用于获取CT图像和CART0电生理解剖图，所述CART0电生理 解剖图包括CART0电解剖图和CART0阻抗图；所述CT图像是对CT原始图像进行重建分割 之后获得的心房内壁的曲面图像；所述CART0阻抗图是在CART0电解剖图上叠加阻抗信息 生成的，有着与CART0电解剖图完全相同的空间信息，同时，CART0阻抗图中的阻抗值与所 述CT图像的CT值相关；
[0016] 粗配准单元，用于配准所述CART0电解剖图与CT图像，得到粗配准变换；
[0017] 精配准单元，用于以所述粗配准变换作为初始变换，配准所述CART0阻抗图与CT 图像，得到精配准变换；
[0018] 输出配准图像单元，用于依据所述精配准变换对所述CART0电解剖图作变换，将 变换后的CART0电解剖图叠加在所述CT图像上，得到最终配准的CART0电解剖图与CT图 像。
[0019] 相较于现有技术，本发明具有以下有益效果：
[0020] 相对现有技术，本发明提供的CART0电解剖图和CT图像配准方法和装置不仅利用 了CART0电解剖图所提供的空间信息，同时，也充分利用了CART0阻抗图所提供的阻抗信 息，有效克服了CART0电解剖图中标测点少、信息量不足的缺陷，使得配准精度得到很大提 高。本发明最终获得的CART0电解剖图与CT图像配准精度能够满足房颤消融导管精确导 航的要求，从而保证了导管消融位点的准确性。
【附图说明】
[0021] 为了清楚地理解本发明的技术方案，下面对描述本发明【具体实施方式】时用到的附 图进行简要说明。显而易见，这些附图仅是本发明的一部分附图，本领域普通技术人员在不 付出创造性劳动的前提下，还可以获得其它附图。
[0022] 图1是本发明实施例提供的CART0电解剖图与CT图像配准方法的流程示意图；
[0023] 图2是本发明实施例提供的CART0电解剖图与CT图像粗配准过程的流程示意图；
[0024] 图3是本发明实施例提供的CART0阻抗图与CT图像精配准过程的流程示意图；
[0025] 图4是本发明实施例提供的CART0电解剖图与CT图像配准装置的结构示意图。
【具体实施方式】
[0026] 为使本发明的目的、技术手段和有益效果更加清楚完整，下面结合附图对本发明 的【具体实施方式】进行描述。
[0027] 在介绍本发明的【具体实施方式】之前，首先介绍CART0电生理解剖图和CT图像的生 成过程。
[0028] CART0电生理解剖图是在房颤消融手术中实时生成的。手术过程中，医生在X线电 视监控下，用一根电极导管从静脉插入左心房，导管在心腔内壁逐点标测，这样依次获得数 百个标测点，导管头采集的磁场信号返回有标测点的空间位置信息，依此构建出心腔的空 间三维结构模型，即CART0电解剖图。导管头采集的磁场信号同时还记录了心房不同部位 的阻抗大小、时间激动的早晚顺序等电生理信息，这些电生理信息以颜色变化叠加在CART0 电解剖图上，分别表现为显示阻抗信息的阻抗图、显示激动顺序信息的激动时间图等。这些 图均独立显示，统称为CART0电生理解剖图。换句话说，CART0电生理解剖图不仅包括显示 空间位置信息的CART0电解剖图，还包括显示阻抗信息的的CART0阻抗图，以及显示激动顺 序信息的CART0激动时间图等。
[0029] CT图像则是在消融手术前生成的。在进行消融手术之前，先对患者进行CT扫描， 得到CT原始图像；然后在CART0标测系统中，利用相应软件对CT原始图像进行重建分割， 得到心房内壁的曲面图像，即本发明所用CT图像。
[0030] 从上述生成CART0电生理解剖图的过程中，可以看到，CART0阻抗图和CART0电解 剖图各自独立存在，CART0阻抗图实质是叠加了阻抗信息的CART0电解剖图，两者有着完全 相同的空间解剖结构信息。因此，配准了CART0阻抗图与CT图像，也就是配准了CART0电 解剖图与CT图像，反之亦然。
[0031] 图1是本发明实施例提供的CART0电解剖图与CT图像配准方法的流程示意图。如 图1所示，该配准方法包括以下步骤：
[0032] S101、获取CT图像A和CART0电生理解剖图B， CART0电生理解剖图B包括CART0 电解剖图&和CART0阻抗图B2:
[0033] 首先，利用CART0三维标测系统，生成CART0电生理解剖图和CT图像；然后，将 CART0电生理解剖图和CT图像导入普通PC机内，在PC机上实施CART0电解剖图与CT图像 配准。
[0034] S102、配准所述CART0电解剖图&与CT图像A，得到粗配准变换：
[0035] 首先，从CART0电生理解剖图B中分离出CART0电解剖图&和CART0阻抗图B2。
[0036] 其次，采用刚体变换模型，以基于主轴的方法配准CART0电解剖图^与^!图像A， 得到粗配准变换。
[0037] CART0电解剖图与CT图像均是形状完整的封闭目标，适合采用基于主轴的配准算 法。同时，为了减少医生与病人在X线下的暴露时间，提高配准效率，本发明选择了刚体变 换模型。由于CT图像比CART0电解剖图分辨率高，细节信息更丰富，所以，本发明将CT图 像作为参考图像，CART0电解剖图作为浮动图像。
[0038] 如图2所示，粗配准过程具体具体包括以下步骤：
[0039] S201、计算CT图像A的质心CA、惯性矩阵IA和特征矩阵EA;计算CART0电解剖图 的质心CB1、惯性矩阵IB1和特征矩阵EB1。
[0040] S202、计算与CT图像具有主轴对应关系的CART0电解剖图B/ :
[0041] 惯性矩阵中的每一个特征向量对应一个主轴，坐标系中的任何坐标轴旋转 180度，并不会影响图像的惯性矩阵。因此，设定CT图像的特征矩阵￡4中的特征向量 排列固定不变，与之对应，构成CARTO电解剖图的特征矩阵EB1的特征向量排列方式有
种，即CARTO电解剖图的特征矩阵EB1共有48种。对每一个E8汾别计算对 应的CARTO电解剖图B/，计算公式如下：
[0043] 其中，私是特征矩阵EB1的转置矩阵；
[0044]S203、计算并比较48个CART0电解剖图B/到所述CT图像A的距离，与CT图像 距离最小的CART0电解剖图&'即是粗配准后的CART0电解剖图B#，其特征矩阵为EB1* :
[0045] 在本发明实施例中，CART0电解剖图B1'到所述CT图像A的距离可以定义为所有 CART0标测点到CT图像最近点欧式距离的均方根。因为仅需遍历48种主轴的对应关系，寻 找其中使得两图距离最小的一种作为最终的配准结果，所以粗配准算法中不涉及任何优化 搜索算法。
[0046]S204、根据特征矩阵EB1*计算刚性的粗配准变换I；，计算公式如下：
[0048] 其中，R为旋转变换，Tr为平移变换。
[0049]S103、以所述粗配准变换作为初始变换，配准所述CART0阻抗图与CT图像，得到精 配准变换：
[0050] 如前所述，CART0阻抗图实质是叠加了阻抗信息的CART0电解剖图，阻抗图与电解 剖图有着完全相同的空间信息。所以，配准了CART0阻抗图与CT图像，也就是配准了CART0 电解剖图与CT图像，反之，配准了CART0电解剖图与CT图像，也就是配准了CART0阻抗图 与CT图像。
[0051] 发明人经过研宄发现，在众多的电生理信息中，心房的阻抗值与CT值之间存在很 大的相关性，可以运用互信息度量CART0阻抗图与CT图像之间的相似度，以基于互信息的 方法配准CART0阻抗图和CT图像，进而实现CART0电解剖图与CT图像的配准。
[0052] 因为心脏是弹性体，考虑速度与精度的平衡，本步骤采用了不同于粗配准刚体变 换的仿射变换模型。
[0053] 如图3所示，精配准过程具体包括以下步骤：
[0054]S301、对CART0阻抗图B2作仿射变换，得到仿射变换后的CART0阻抗图。
[0055] S302、计算所述仿射变换后的CART0阻抗图与CT图像A之间的互信息。
[0056]S303、采用单纯形法，以所述粗配准变换作为优化搜索的初始变换，在仿射变换空 间域内搜索，寻找使得所述互信息达到最大值的仿射变换T*，即为精配准变换。
[0057]S104、依据上述精配准变换，对CART0电解剖图&作变换，将变换后的CART0电解 剖图&叠加在CT图像A上，得到最终配准的CART0电解剖图与CT图像。
[0058] 以上为本发明实施例提供的CART0电解剖图与CT图像配准方法的【具体实施方式】。
[0059] 相对现有技术仅利用空间信息的方法，本发明提供的CART0电解剖图和CT图像配 准方法不仅利用了CART0电解剖图所提供的空间信息，同时，也充分利用了CART0阻抗图所 提供的阻抗信息，有效克服了CART0电解剖图中标测点少、信息量不足的缺陷，使得配准精 度得到很大程度提高。本发明最终得到的CART0电解剖图与CT图像配准精度能够满足房 颤消融导管精确导航的要求，从而保证了导管消融位点的准确性。
[0060] 基于上述实施例提供的CART0电解剖图与CT图像的配准方法，本发明实施例还提 供了CART0电解剖图与CT图像的配准装置。具体参见以下实施例。
[0061] 图4是本发明实施例提供的CART0电解剖图与CT图像配准装置的结构示意图。如 图4所示，该配置装置包括以下单元：
[0062] 获取图像数据单元41，用于获取CT图像和CART0电生理解剖图，所述CART0电生 理解剖图包括CART0电解剖图和CART0阻抗图；所述CT图像是对CT原始图像进行重建分 割之后获得的心房内壁的曲面图像；所述CART0阻抗图是在CART0电解剖图上叠加阻抗信 息生成的，有着与CART0电解剖图完全相同的空间信息，同时，CART0阻抗图中的阻抗值与 所述CT图像的CT值相关；
[0063] 粗配准单元42,用于配准所述CART0电解剖图与CT图像，得到粗配准变换；
[0064] 精配准单元43,用于以所述粗配准变换作为初始变换，配准所述CART0阻抗图与 CT图像，得到精配准变换；
[0065] 输出配准图像单元44,用于依据所述精配准变换对所述CART0电解剖图作变换， 将变换后的CART0电解剖图叠加在所述CT图像上，得到最终配准的CART0电解剖图与CT 图像。
[0066] 通过上述实施例提供的配置装置，最终获得的配准后的CART0电解剖图与CT图像 配准精度能够满足消融导管精确导航要求，能够保证导管消融位点的准确性。
[0067] 进一步地，所述粗配准单元42包括采用刚体变换模型，以基于主轴的配准方法粗 配准CART0电解剖图与CT图像，得到粗配准变换的子单元。
[0068] 更进一步地，所述粗配准单元42包括：
[0069] 第一计算子单元421，用于计算CT图像A的质心CA、惯性矩阵IA和特征矩阵EA;计 算CART0电解剖图&的质心CB1、惯性矩阵IB1和特征矩阵EB1;
[0070] 第二计算子单元422,用于计算每一个EB1对应的CART0电解剖图B/;设定CT图 像A的特征矩阵￡4固定不变，每一个特征向量对应一个主轴，与CT图像主轴对应的CART0 电解剖图主轴有48种排列方式，即CART0电解剖图特征矩阵EB1的特征向量有48种排列； 对每一个EB1分别计算对应的CART0电解剖图B/，计算公式如下：
[0072] 其中，砥是特征矩阵EB1的转置矩阵；
[0073] 第三计算子单元423,用于计算并比较48个CART0电解剖图B/到所述CT图像A 的距离；与CT图像A距离最小的CART0电解剖图B/即是粗配准后的CART0电解剖图&*， 其特征矩阵为;
[0074] 第四计算子单元424,用于根据特征矩阵计算刚性的粗配准变换L，计算公式 如下：
[0076] 其中，R为旋转变换，I；为平移变换。
[0077] 再进一步地，所述精配准单元43包括采用仿射变换模型，以基于互信息的方法配 准所述CARTO阻抗图与CT图像，得到精配准变换的子单元。
[0078] 更进一步地，所述精配准单元43包括：
[0079] 仿射变换子单元431，用于对CART0阻抗图做仿射变换，得到仿射变换后的CART0 阻抗图；
[0080] 第五计算子单元432,用于计算所述仿射变换后的CART0阻抗图与所述CT图像之 间的互信息；
[0081] 优化搜索子单元433,用于采用单纯形法，以粗配准变换作为初始变换，在仿射变 换空间域内搜索，寻找使得所述互信息达到最大值的仿射变换T*，所述仿射变换T*即为精 配准变换。
[0082] 以上所述仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制。虽 然本发明以较佳实施例揭露如上，然而并非用以限定本发明。任何熟悉本领域的技术人员， 在不脱离本发明技术方案范围情况下，都可利用上述揭示的方法和技术内容对本发明技术 方案做出许多可能的变动和修饰，或修改为等同变化的等效实施例。因此，凡是未脱离本发 明技术方案的内容，依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化 及修饰，均仍属于本发明技术方案保护的范围内。
【主权项】
1. 一种CARTO电解剖图与CT图像的配准方法，其特征在于，包括： 获取CT图像和CARTO电生理解剖图，所述CARTO电生理解剖图包括CARTO电解剖图和 CARTO阻抗图；所述CT图像是对CT原始图像进行重建分割之后获得的心房内壁的曲面图 像；所述CARTO阻抗图是在CARTO电解剖图上叠加阻抗信息生成的，有着与CARTO电解剖图 完全相同的空间信息，同时，CARTO阻抗图中的阻抗值与所述CT图像的CT值相关； 配准所述CARTO电解剖图与CT图像，得到粗配准变换； 以所述粗配准变换作为初始变换，配准所述CARTO阻抗图与CT图像，得到精配准变 换； 依据所述精配准变换对所述CARTO电解剖图作变换，将变换后的CARTO电解剖图叠加 在所述CT图像上，得到最终配准的CARTO电解剖图与CT图像。2. 根据权利要求1所述的方法，其特征在于，配准所述CARTO电解剖图与CT图像，得到 粗配准变换，具体包括： 采用刚体变换模型，以基于主轴的配准方法粗配准CARTO电解剖图与CT图像，得到粗 配准变换。3. 根据权利要求2所述的方法，其特征在于，采用刚体变换模型，以基于主轴的方法粗 配准CARTO电解剖图与CT图像，得到粗配准变换，具体包括： A1、计算CT图像A的质心Ca、惯性矩阵Ia和特征矩阵E A;计算CARTO电解剖图B i的质 心Cbi、惯性矩阵Ibi和特征矩阵E B1; A2、设定CT图像A的特征矩阵EA固定不变，每一个特征向量对应一个主轴，与CT图像 主轴对应的CARTO电解剖图主轴有48种排列方式，即CARTO电解剖图特征矩阵Ebi的特征 向量有48种排列；对每一个Ebi分别计算对应的CARTO电解剖图B i'，计算公式如下： B1'= Ea *ETm *{Im-Cm) + CA ； 其中，矻是特征矩阵Ebi的转置矩阵； A3、计算并比较48个CARTO电解剖图B/到所述CT图像A的距离，与CT图像A距离 最小的CARTO电解剖图B/即是粗配准后CARTO电解剖图瓦，其特征矩阵为, A4、根据特征矩阵十算刚性的粗配准变换Ttl，计算公式如下： T0= [R,Tr] ； R-EA^ETBf； xr=CA-R*CB1； 其中，R为旋转变换，?；为平移变换。4. 根据权利要求1所述的方法，其特征在于，配准所述CARTO阻抗图与CT图像，得到精 配准变换，具体包括； 采用仿射变换模型，以基于互信息的方法配准所述CARTO阻抗图与CT图像，得到精配 准变换。5. 根据权利要求4所述的方法，其特征在于，采用仿射变换模型，以基于互信息的方法 配准所述CARTO阻抗图与CT图像，得到精配准变换，具体包括： Β1、对CARTO阻抗图作仿射变换，得到仿射变换后的CARTO阻抗图； Β2、计算仿射变换后的CARTO阻抗图氏与CT图像A之间的互信息； Β3、采用单纯形法，以粗配准变换作为初始变换，在仿射变换空间域内搜索，寻找使得 所述互信息达到最大值的仿射变换T*，所述仿射变换T*即为精配准变换。6. -种CARTO电解剖图与CT图像的配准装置，其特征在于，包括： 获取图像数据单元，用于获取CT图像和CARTO电生理解剖图，所述CARTO电生理解剖 图包括CARTO电解剖图和CARTO阻抗图；所述CT图像是对CT原始图像进行重建分割之后 获得的心房内壁的曲面图像；所述CARTO阻抗图是在CARTO电解剖图上叠加阻抗信息生成 的，有着与CARTO电解剖图完全相同的空间信息，同时，CARTO阻抗图中的阻抗值与所述CT 图像的CT值相关； 粗配准单元，用于配准所述CARTO电解剖图与CT图像，得到粗配准变换； 精配准单元，用于以所述粗配准变换作为初始变换，配准所述CARTO阻抗图与CT图像， 得到精配准变换； 输出配准图像单元，用于依据所述精配准变换对所述CARTO电解剖图作变换，将变换 后的CARTO电解剖图叠加在所述CT图像上，得到最终配准的CARTO电解剖图与CT图像。7. 根据权利要求6所述的装置，其特征在于，所述粗配准单元包括： 第一计算子单元，用于计算CT图像A的质心Ca、惯性矩阵Ia和特征矩阵Ea;计算CARTO 电解剖图B1的质心C B1、惯性矩阵Ibi和特征矩阵E Β1; 第二计算子单元，用于计算每一个Ebi对应的CARTO电解剖图B 设定CT图像A的特 征矩阵EA固定不变，每一个特征向量对应一个主轴，与CT图像主轴对应的CARTO电解剖图 主轴有48种排列方式，即CARTO电解剖图特征矩阵Ebi的特征向量有48种排列；对每一个 Ebi分别计算对应的CARTO电解剖图B /，计算公式如下： Βχ' = Εα*ΕΙχ*{Ιβι-€β,) + €α· 其中，砥是特征矩阵Ebi的转置矩阵； 第三计算子单元，用于计算并比较48个CARTO电解剖图B/到所述CT图像A的距离； 与CT图像A距离最小的CARTO电解剖图B1'即是粗配准后的CARTO电解剖图B1*,其特征 矩阵为; 第四计算子单元，用于根据特征矩阵计算刚性的粗配准变换Ttl，计算公式如下： T0= [R,Tr] ,R^E^ETm，] Tr=CA-R*CB1； 其中，R为旋转变换，?；为平移变换。8. 根据权利要求6或7所述的装置，其特征在于，所述精配准单元包括采用仿射变换模 型，以基于互信息的方法配准所述CARTO阻抗图与CT图像，得到精配准变换的子单元。9. 根据权利要求8所述的装置，其特征在于，所述精配准单元包括： 仿射变换子单元，用于对CARTO阻抗图做仿射变换，得到仿射变换后的CARTO阻抗图； 第五计算子单元，用于计算所述仿射变换后的CARTO阻抗图与所述CT图像之间的互信 息； 优化搜索子单元，用于采用单纯形法，以粗配准变换作为初始变换，在仿射变换空间域 内搜索，寻找使得所述互信息达到最大值的仿射变换Τ*，所述仿射变换Τ*即为精配准变 换。
【专利摘要】本发明公开了一种基于空间和阻抗的CARTO电解剖图与CT图像的配准方法和装置。CARTO电生理解剖图包括电解剖图、阻抗图等，电解剖图与阻抗图有着完全相同的空间信息。首先，利用空间信息，配准CARTO电解剖图与CT图像，得到粗配准变换；然后，以粗配准变换为初始变换,利用阻抗信息，配准CARTO阻抗图与CT图像，得到精配准变换；最后，依据所述精配准变换对CARTO电解剖图作变换，将变换后的电解剖图叠加在CT图像上，得到配准后的CARTO电解剖图与CT图像。相对现有方法，本发明既利用了空间信息，也利用了阻抗信息，配准后的CARTO电解剖图与CT图像的配准精度得到很大提高，可以为房颤射频消融手术中的导管提供精确导航。
【IPC分类】G06T7/00
【公开号】CN104899886
【申请号】CN201510311980
【发明人】舒丽霞, 蔺嫦燕
【申请人】首都医科大学附属北京安贞医院
【公开日】2015年9月9日
【申请日】2015年6月8日
【技术领域】
[0001] 本发明涉及手术导航中图像配准技术领域，尤其涉及心房颤动（房颤）射频消融 介入手术导航中的CART0电解剖图与CT图像的配准方法和装置。
【背景技术】
[0002] 房颤是临床上最常见的心律失常，有着较高的致残率和致死率。三维标测和消融 系统CART0是治疗房颤的首选手术平台。
[0003] CART0电解剖图是消融手术中实时模拟生成的心房解剖结构图，依据医生在心房 内壁标测的若干个标测点构建。医生可以直接在CART0电解剖图上操作导管实施消融，然 而，由于模拟构建CART0电解剖图的标测点较少，导致电解剖图分辨率低，肺静脉的数目、 具体开口部分及每支肺静脉各自的分支情况等等这些复杂的解剖结构无法在电解剖图上 清晰地显示出来，医生难以准确定位并有效隔离肺静脉，致使房颤复发率高，并可能伴随肺 静脉狭窄等严重并发症，从而极大地影响消融治疗效果。
[0004] 众所周知，CT图像分辨率高，能清晰地显示上述消融靶区的复杂结构，然而，CT图 像只能在术前获得，医生不能在上面直接实施导管操作。
[0005] 配准CART0电解剖图和CT图像，可帮助实时并清楚地显示心腔和肺静脉的解剖结 构。依据两者的配准图像，医生将导管头固定在发生房颤的肺静脉前庭，在上面设定围绕左 或右侧上、下肺静脉的环状消融线；导管头释放射频电能，逐点消融，把导致房颤的纤维组 织烧掉，直至肺静脉出现与左心房无关的自律性电位。因此，CART0电解剖图和CT图像的 配准图像能够让医生从多个角度观察消融靶区的结构，从而增加消融位点的精确性，减少 手术相关并发症，缩短手术时间和X线曝光时间。
[0006] 然而，现有的CART0电解剖图和CT图像配准技术达到的配准精度均难以满足消融 导管精确导航的要求。
【发明内容】
[0007] 有鉴于此，本发明提供了一种房颤射频消融介入手术导航中的CART0电解剖图与 CT图像的配准方法和装置，以提高CART0电解剖图与CT图像的配准精度，为房颤消融导管 提供精确导航。
[0008] 需要说明的是，CART0电生理解剖图包括有CART0电解剖图、CART0阻抗图等； CART0阻抗图是在CART0电解剖图上叠加阻抗信息生成的，其实质是叠加了阻抗信息的 CART0电解剖图，阻抗图与电解剖图有着完全相同的空间信息。因此，配准了CART0阻抗图 与CT图像，也就是配准了CART0电解剖图与CT图像。
[0009] 为了达到上述发明目的，本发明的第一方面提出了一种CART0电解剖图与CT图像 的配准方法，包括如下步骤：
[0010] 获取CT图像和CART0电生理解剖图，所述CART0电生理解剖图包括CART0电解剖 图和CARTO阻抗图；所述CT图像是对CT原始图像进行重建分割之后获得的心房内壁的曲 面图像；所述CART0阻抗图是在CART0电解剖图上叠加阻抗信息生成的，有着与CART0电解 剖图完全相同的空间信息，同时，CART0阻抗图中的阻抗值与所述CT图像的CT值相关；
[0011] 配准所述CART0电解剖图与CT图像，得到粗配准变换；
[0012] 以所述粗配准变换作为初始变换，配准所述CART0阻抗图与CT图像，得到精配准 变换；
[0013] 依据所述精配准变换对所述CART0电解剖图作变换，将变换后的CART0电解剖图 叠加在所述CT图像上，得到最终配准的CART0电解剖图与CT图像。
[0014] 基于本发明的第一方面提出的CART0电解剖图与CT图像的配准方法，本发明同时 还提出了一种CART0电解剖图与CT图像的配准装置，包括：
[0015] 获取图像数据单元，用于获取CT图像和CART0电生理解剖图，所述CART0电生理 解剖图包括CART0电解剖图和CART0阻抗图；所述CT图像是对CT原始图像进行重建分割 之后获得的心房内壁的曲面图像；所述CART0阻抗图是在CART0电解剖图上叠加阻抗信息 生成的，有着与CART0电解剖图完全相同的空间信息，同时，CART0阻抗图中的阻抗值与所 述CT图像的CT值相关；
[0016] 粗配准单元，用于配准所述CART0电解剖图与CT图像，得到粗配准变换；
[0017] 精配准单元，用于以所述粗配准变换作为初始变换，配准所述CART0阻抗图与CT 图像，得到精配准变换；
[0018] 输出配准图像单元，用于依据所述精配准变换对所述CART0电解剖图作变换，将 变换后的CART0电解剖图叠加在所述CT图像上，得到最终配准的CART0电解剖图与CT图 像。
[0019] 相较于现有技术，本发明具有以下有益效果：
[0020] 相对现有技术，本发明提供的CART0电解剖图和CT图像配准方法和装置不仅利用 了CART0电解剖图所提供的空间信息，同时，也充分利用了CART0阻抗图所提供的阻抗信 息，有效克服了CART0电解剖图中标测点少、信息量不足的缺陷，使得配准精度得到很大提 高。本发明最终获得的CART0电解剖图与CT图像配准精度能够满足房颤消融导管精确导 航的要求，从而保证了导管消融位点的准确性。
【附图说明】
[0021] 为了清楚地理解本发明的技术方案，下面对描述本发明【具体实施方式】时用到的附 图进行简要说明。显而易见，这些附图仅是本发明的一部分附图，本领域普通技术人员在不 付出创造性劳动的前提下，还可以获得其它附图。
[0022] 图1是本发明实施例提供的CART0电解剖图与CT图像配准方法的流程示意图；
[0023] 图2是本发明实施例提供的CART0电解剖图与CT图像粗配准过程的流程示意图；
[0024] 图3是本发明实施例提供的CART0阻抗图与CT图像精配准过程的流程示意图；
[0025] 图4是本发明实施例提供的CART0电解剖图与CT图像配准装置的结构示意图。
【具体实施方式】
[0026] 为使本发明的目的、技术手段和有益效果更加清楚完整，下面结合附图对本发明 的【具体实施方式】进行描述。
[0027] 在介绍本发明的【具体实施方式】之前，首先介绍CART0电生理解剖图和CT图像的生 成过程。
[0028] CART0电生理解剖图是在房颤消融手术中实时生成的。手术过程中，医生在X线电 视监控下，用一根电极导管从静脉插入左心房，导管在心腔内壁逐点标测，这样依次获得数 百个标测点，导管头采集的磁场信号返回有标测点的空间位置信息，依此构建出心腔的空 间三维结构模型，即CART0电解剖图。导管头采集的磁场信号同时还记录了心房不同部位 的阻抗大小、时间激动的早晚顺序等电生理信息，这些电生理信息以颜色变化叠加在CART0 电解剖图上，分别表现为显示阻抗信息的阻抗图、显示激动顺序信息的激动时间图等。这些 图均独立显示，统称为CART0电生理解剖图。换句话说，CART0电生理解剖图不仅包括显示 空间位置信息的CART0电解剖图，还包括显示阻抗信息的的CART0阻抗图，以及显示激动顺 序信息的CART0激动时间图等。
[0029] CT图像则是在消融手术前生成的。在进行消融手术之前，先对患者进行CT扫描， 得到CT原始图像；然后在CART0标测系统中，利用相应软件对CT原始图像进行重建分割， 得到心房内壁的曲面图像，即本发明所用CT图像。
[0030] 从上述生成CART0电生理解剖图的过程中，可以看到，CART0阻抗图和CART0电解 剖图各自独立存在，CART0阻抗图实质是叠加了阻抗信息的CART0电解剖图，两者有着完全 相同的空间解剖结构信息。因此，配准了CART0阻抗图与CT图像，也就是配准了CART0电 解剖图与CT图像，反之亦然。
[0031] 图1是本发明实施例提供的CART0电解剖图与CT图像配准方法的流程示意图。如 图1所示，该配准方法包括以下步骤：
[0032] S101、获取CT图像A和CART0电生理解剖图B， CART0电生理解剖图B包括CART0 电解剖图&和CART0阻抗图B2:
[0033] 首先，利用CART0三维标测系统，生成CART0电生理解剖图和CT图像；然后，将 CART0电生理解剖图和CT图像导入普通PC机内，在PC机上实施CART0电解剖图与CT图像 配准。
[0034] S102、配准所述CART0电解剖图&与CT图像A，得到粗配准变换：
[0035] 首先，从CART0电生理解剖图B中分离出CART0电解剖图&和CART0阻抗图B2。
[0036] 其次，采用刚体变换模型，以基于主轴的方法配准CART0电解剖图^与^!图像A， 得到粗配准变换。
[0037] CART0电解剖图与CT图像均是形状完整的封闭目标，适合采用基于主轴的配准算 法。同时，为了减少医生与病人在X线下的暴露时间，提高配准效率，本发明选择了刚体变 换模型。由于CT图像比CART0电解剖图分辨率高，细节信息更丰富，所以，本发明将CT图 像作为参考图像，CART0电解剖图作为浮动图像。
[0038] 如图2所示，粗配准过程具体具体包括以下步骤：
[0039] S201、计算CT图像A的质心CA、惯性矩阵IA和特征矩阵EA;计算CART0电解剖图 的质心CB1、惯性矩阵IB1和特征矩阵EB1。
[0040] S202、计算与CT图像具有主轴对应关系的CART0电解剖图B/ :
[0041] 惯性矩阵中的每一个特征向量对应一个主轴，坐标系中的任何坐标轴旋转 180度，并不会影响图像的惯性矩阵。因此，设定CT图像的特征矩阵￡4中的特征向量 排列固定不变，与之对应，构成CARTO电解剖图的特征矩阵EB1的特征向量排列方式有
种，即CARTO电解剖图的特征矩阵EB1共有48种。对每一个E8汾别计算对 应的CARTO电解剖图B/，计算公式如下：
[0043] 其中，私是特征矩阵EB1的转置矩阵；
[0044]S203、计算并比较48个CART0电解剖图B/到所述CT图像A的距离，与CT图像 距离最小的CART0电解剖图&'即是粗配准后的CART0电解剖图B#，其特征矩阵为EB1* :
[0045] 在本发明实施例中，CART0电解剖图B1'到所述CT图像A的距离可以定义为所有 CART0标测点到CT图像最近点欧式距离的均方根。因为仅需遍历48种主轴的对应关系，寻 找其中使得两图距离最小的一种作为最终的配准结果，所以粗配准算法中不涉及任何优化 搜索算法。
[0046]S204、根据特征矩阵EB1*计算刚性的粗配准变换I；，计算公式如下：
[0048] 其中，R为旋转变换，Tr为平移变换。
[0049]S103、以所述粗配准变换作为初始变换，配准所述CART0阻抗图与CT图像，得到精 配准变换：
[0050] 如前所述，CART0阻抗图实质是叠加了阻抗信息的CART0电解剖图，阻抗图与电解 剖图有着完全相同的空间信息。所以，配准了CART0阻抗图与CT图像，也就是配准了CART0 电解剖图与CT图像，反之，配准了CART0电解剖图与CT图像，也就是配准了CART0阻抗图 与CT图像。
[0051] 发明人经过研宄发现，在众多的电生理信息中，心房的阻抗值与CT值之间存在很 大的相关性，可以运用互信息度量CART0阻抗图与CT图像之间的相似度，以基于互信息的 方法配准CART0阻抗图和CT图像，进而实现CART0电解剖图与CT图像的配准。
[0052] 因为心脏是弹性体，考虑速度与精度的平衡，本步骤采用了不同于粗配准刚体变 换的仿射变换模型。
[0053] 如图3所示，精配准过程具体包括以下步骤：
[0054]S301、对CART0阻抗图B2作仿射变换，得到仿射变换后的CART0阻抗图。
[0055] S302、计算所述仿射变换后的CART0阻抗图与CT图像A之间的互信息。
[0056]S303、采用单纯形法，以所述粗配准变换作为优化搜索的初始变换，在仿射变换空 间域内搜索，寻找使得所述互信息达到最大值的仿射变换T*，即为精配准变换。
[0057]S104、依据上述精配准变换，对CART0电解剖图&作变换，将变换后的CART0电解 剖图&叠加在CT图像A上，得到最终配准的CART0电解剖图与CT图像。
[0058] 以上为本发明实施例提供的CART0电解剖图与CT图像配准方法的【具体实施方式】。
[0059] 相对现有技术仅利用空间信息的方法，本发明提供的CART0电解剖图和CT图像配 准方法不仅利用了CART0电解剖图所提供的空间信息，同时，也充分利用了CART0阻抗图所 提供的阻抗信息，有效克服了CART0电解剖图中标测点少、信息量不足的缺陷，使得配准精 度得到很大程度提高。本发明最终得到的CART0电解剖图与CT图像配准精度能够满足房 颤消融导管精确导航的要求，从而保证了导管消融位点的准确性。
[0060] 基于上述实施例提供的CART0电解剖图与CT图像的配准方法，本发明实施例还提 供了CART0电解剖图与CT图像的配准装置。具体参见以下实施例。
[0061] 图4是本发明实施例提供的CART0电解剖图与CT图像配准装置的结构示意图。如 图4所示，该配置装置包括以下单元：
[0062] 获取图像数据单元41，用于获取CT图像和CART0电生理解剖图，所述CART0电生 理解剖图包括CART0电解剖图和CART0阻抗图；所述CT图像是对CT原始图像进行重建分 割之后获得的心房内壁的曲面图像；所述CART0阻抗图是在CART0电解剖图上叠加阻抗信 息生成的，有着与CART0电解剖图完全相同的空间信息，同时，CART0阻抗图中的阻抗值与 所述CT图像的CT值相关；
[0063] 粗配准单元42,用于配准所述CART0电解剖图与CT图像，得到粗配准变换；
[0064] 精配准单元43,用于以所述粗配准变换作为初始变换，配准所述CART0阻抗图与 CT图像，得到精配准变换；
[0065] 输出配准图像单元44,用于依据所述精配准变换对所述CART0电解剖图作变换， 将变换后的CART0电解剖图叠加在所述CT图像上，得到最终配准的CART0电解剖图与CT 图像。
[0066] 通过上述实施例提供的配置装置，最终获得的配准后的CART0电解剖图与CT图像 配准精度能够满足消融导管精确导航要求，能够保证导管消融位点的准确性。
[0067] 进一步地，所述粗配准单元42包括采用刚体变换模型，以基于主轴的配准方法粗 配准CART0电解剖图与CT图像，得到粗配准变换的子单元。
[0068] 更进一步地，所述粗配准单元42包括：
[0069] 第一计算子单元421，用于计算CT图像A的质心CA、惯性矩阵IA和特征矩阵EA;计 算CART0电解剖图&的质心CB1、惯性矩阵IB1和特征矩阵EB1;
[0070] 第二计算子单元422,用于计算每一个EB1对应的CART0电解剖图B/;设定CT图 像A的特征矩阵￡4固定不变，每一个特征向量对应一个主轴，与CT图像主轴对应的CART0 电解剖图主轴有48种排列方式，即CART0电解剖图特征矩阵EB1的特征向量有48种排列； 对每一个EB1分别计算对应的CART0电解剖图B/，计算公式如下：
[0072] 其中，砥是特征矩阵EB1的转置矩阵；
[0073] 第三计算子单元423,用于计算并比较48个CART0电解剖图B/到所述CT图像A 的距离；与CT图像A距离最小的CART0电解剖图B/即是粗配准后的CART0电解剖图&*， 其特征矩阵为;
[0074] 第四计算子单元424,用于根据特征矩阵计算刚性的粗配准变换L，计算公式 如下：
[0076] 其中，R为旋转变换，I；为平移变换。
[0077] 再进一步地，所述精配准单元43包括采用仿射变换模型，以基于互信息的方法配 准所述CARTO阻抗图与CT图像，得到精配准变换的子单元。
[0078] 更进一步地，所述精配准单元43包括：
[0079] 仿射变换子单元431，用于对CART0阻抗图做仿射变换，得到仿射变换后的CART0 阻抗图；
[0080] 第五计算子单元432,用于计算所述仿射变换后的CART0阻抗图与所述CT图像之 间的互信息；
[0081] 优化搜索子单元433,用于采用单纯形法，以粗配准变换作为初始变换，在仿射变 换空间域内搜索，寻找使得所述互信息达到最大值的仿射变换T*，所述仿射变换T*即为精 配准变换。
[0082] 以上所述仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制。虽 然本发明以较佳实施例揭露如上，然而并非用以限定本发明。任何熟悉本领域的技术人员， 在不脱离本发明技术方案范围情况下，都可利用上述揭示的方法和技术内容对本发明技术 方案做出许多可能的变动和修饰，或修改为等同变化的等效实施例。因此，凡是未脱离本发 明技术方案的内容，依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化 及修饰，均仍属于本发明技术方案保护的范围内。
【主权项】
1. 一种CARTO电解剖图与CT图像的配准方法，其特征在于，包括： 获取CT图像和CARTO电生理解剖图，所述CARTO电生理解剖图包括CARTO电解剖图和 CARTO阻抗图；所述CT图像是对CT原始图像进行重建分割之后获得的心房内壁的曲面图 像；所述CARTO阻抗图是在CARTO电解剖图上叠加阻抗信息生成的，有着与CARTO电解剖图 完全相同的空间信息，同时，CARTO阻抗图中的阻抗值与所述CT图像的CT值相关； 配准所述CARTO电解剖图与CT图像，得到粗配准变换； 以所述粗配准变换作为初始变换，配准所述CARTO阻抗图与CT图像，得到精配准变 换； 依据所述精配准变换对所述CARTO电解剖图作变换，将变换后的CARTO电解剖图叠加 在所述CT图像上，得到最终配准的CARTO电解剖图与CT图像。2. 根据权利要求1所述的方法，其特征在于，配准所述CARTO电解剖图与CT图像，得到 粗配准变换，具体包括： 采用刚体变换模型，以基于主轴的配准方法粗配准CARTO电解剖图与CT图像，得到粗 配准变换。3. 根据权利要求2所述的方法，其特征在于，采用刚体变换模型，以基于主轴的方法粗 配准CARTO电解剖图与CT图像，得到粗配准变换，具体包括： A1、计算CT图像A的质心Ca、惯性矩阵Ia和特征矩阵E A;计算CARTO电解剖图B i的质 心Cbi、惯性矩阵Ibi和特征矩阵E B1; A2、设定CT图像A的特征矩阵EA固定不变，每一个特征向量对应一个主轴，与CT图像 主轴对应的CARTO电解剖图主轴有48种排列方式，即CARTO电解剖图特征矩阵Ebi的特征 向量有48种排列；对每一个Ebi分别计算对应的CARTO电解剖图B i'，计算公式如下： B1'= Ea *ETm *{Im-Cm) + CA ； 其中，矻是特征矩阵Ebi的转置矩阵； A3、计算并比较48个CARTO电解剖图B/到所述CT图像A的距离，与CT图像A距离 最小的CARTO电解剖图B/即是粗配准后CARTO电解剖图瓦，其特征矩阵为, A4、根据特征矩阵十算刚性的粗配准变换Ttl，计算公式如下： T0= [R,Tr] ； R-EA^ETBf； xr=CA-R*CB1； 其中，R为旋转变换，?；为平移变换。4. 根据权利要求1所述的方法，其特征在于，配准所述CARTO阻抗图与CT图像，得到精 配准变换，具体包括； 采用仿射变换模型，以基于互信息的方法配准所述CARTO阻抗图与CT图像，得到精配 准变换。5. 根据权利要求4所述的方法，其特征在于，采用仿射变换模型，以基于互信息的方法 配准所述CARTO阻抗图与CT图像，得到精配准变换，具体包括： Β1、对CARTO阻抗图作仿射变换，得到仿射变换后的CARTO阻抗图； Β2、计算仿射变换后的CARTO阻抗图氏与CT图像A之间的互信息； Β3、采用单纯形法，以粗配准变换作为初始变换，在仿射变换空间域内搜索，寻找使得 所述互信息达到最大值的仿射变换T*，所述仿射变换T*即为精配准变换。6. -种CARTO电解剖图与CT图像的配准装置，其特征在于，包括： 获取图像数据单元，用于获取CT图像和CARTO电生理解剖图，所述CARTO电生理解剖 图包括CARTO电解剖图和CARTO阻抗图；所述CT图像是对CT原始图像进行重建分割之后 获得的心房内壁的曲面图像；所述CARTO阻抗图是在CARTO电解剖图上叠加阻抗信息生成 的，有着与CARTO电解剖图完全相同的空间信息，同时，CARTO阻抗图中的阻抗值与所述CT 图像的CT值相关； 粗配准单元，用于配准所述CARTO电解剖图与CT图像，得到粗配准变换； 精配准单元，用于以所述粗配准变换作为初始变换，配准所述CARTO阻抗图与CT图像， 得到精配准变换； 输出配准图像单元，用于依据所述精配准变换对所述CARTO电解剖图作变换，将变换 后的CARTO电解剖图叠加在所述CT图像上，得到最终配准的CARTO电解剖图与CT图像。7. 根据权利要求6所述的装置，其特征在于，所述粗配准单元包括： 第一计算子单元，用于计算CT图像A的质心Ca、惯性矩阵Ia和特征矩阵Ea;计算CARTO 电解剖图B1的质心C B1、惯性矩阵Ibi和特征矩阵E Β1; 第二计算子单元，用于计算每一个Ebi对应的CARTO电解剖图B 设定CT图像A的特 征矩阵EA固定不变，每一个特征向量对应一个主轴，与CT图像主轴对应的CARTO电解剖图 主轴有48种排列方式，即CARTO电解剖图特征矩阵Ebi的特征向量有48种排列；对每一个 Ebi分别计算对应的CARTO电解剖图B /，计算公式如下： Βχ' = Εα*ΕΙχ*{Ιβι-€β,) + €α· 其中，砥是特征矩阵Ebi的转置矩阵； 第三计算子单元，用于计算并比较48个CARTO电解剖图B/到所述CT图像A的距离； 与CT图像A距离最小的CARTO电解剖图B1'即是粗配准后的CARTO电解剖图B1*,其特征 矩阵为; 第四计算子单元，用于根据特征矩阵计算刚性的粗配准变换Ttl，计算公式如下： T0= [R,Tr] ,R^E^ETm，] Tr=CA-R*CB1； 其中，R为旋转变换，?；为平移变换。8. 根据权利要求6或7所述的装置，其特征在于，所述精配准单元包括采用仿射变换模 型，以基于互信息的方法配准所述CARTO阻抗图与CT图像，得到精配准变换的子单元。9. 根据权利要求8所述的装置，其特征在于，所述精配准单元包括： 仿射变换子单元，用于对CARTO阻抗图做仿射变换，得到仿射变换后的CARTO阻抗图； 第五计算子单元，用于计算所述仿射变换后的CARTO阻抗图与所述CT图像之间的互信 息； 优化搜索子单元，用于采用单纯形法，以粗配准变换作为初始变换，在仿射变换空间域 内搜索，寻找使得所述互信息达到最大值的仿射变换Τ*，所述仿射变换Τ*即为精配准变 换。
【专利摘要】本发明公开了一种基于空间和阻抗的CARTO电解剖图与CT图像的配准方法和装置。CARTO电生理解剖图包括电解剖图、阻抗图等，电解剖图与阻抗图有着完全相同的空间信息。首先，利用空间信息，配准CARTO电解剖图与CT图像，得到粗配准变换；然后，以粗配准变换为初始变换,利用阻抗信息，配准CARTO阻抗图与CT图像，得到精配准变换；最后，依据所述精配准变换对CARTO电解剖图作变换，将变换后的电解剖图叠加在CT图像上，得到配准后的CARTO电解剖图与CT图像。相对现有方法，本发明既利用了空间信息，也利用了阻抗信息，配准后的CARTO电解剖图与CT图像的配准精度得到很大提高，可以为房颤射频消融手术中的导管提供精确导航。
【IPC分类】G06T7/00
【公开号】CN104899886
【申请号】CN201510311980
【发明人】舒丽霞, 蔺嫦燕
【申请人】首都医科大学附属北京安贞医院
【公开日】2015年9月9日
【申请日】2015年6月8日
最新回复(0)