心电图与体表测量的相关性的制作方法

专利名称：心电图与体表测量的相关性的制作方法
技术领域：
本发明涉及心律不齐的评估和治疗。更具体地，本发明涉及改进心脏的电 子绘亂用于在评估4獐不齐及其消融治疗中i柳。
背景技术：
基于体表心电图(ECG)技术、用于心脏中的电位的非侵/i会图的方法是 已知的。这些方法将3维成像与ECG 结合，以,生产生心外膜^M上 以及心内膜表面上电位的3维图。
文献Modre等的"均匀节律数据的心房非侵入^^活绘图(Atrial Noninvasive Activation Mapping of Paced Rththm Data) "J. Cardiovasc. Electrophysiology 14:712-719， (2003年7月)，基于磁共振成像(MRI)和ECG绘图数据，描述 了表面心脏模型激活时间(AT)成像方法。心内膜和心外膜表面舰能以这种 方式绘图。AT模式与心房电位的CARTO，图比较。夕卜部的解剖iH己在内部绘 图后fflil将导管端移动到在体表的iH己位置，用于将CARTO数据偶合到MRI 定位系统。有人建议心房内的AT成像对于患有局灶心律失常患者的心房活性 的非侵入成像是有用的。

发明内容
在本发明的实施例中，心脏的电子图通过两种形式获取 一种相对较大侵 入而另一种较少侵入。在一个实施例中，较大侵入的图劍顿导管绘图系统获 取的心内膜的图。多个较少侵入的图可以基于体表ECG非侵Ai也获取。可替 换地或另外地，较少侵入的图可以4顿心外膜的接触探测器(如穿过胸壁插入 心包的导管)来获取。较少侵入的图典型地显示心外膜上的电位，尽管它们也
可以显示心内膜电位，如在上面提到的文献"均匀节律数据的心房非侵入激活
绘图"(Modre等人)中描述的，较少侵入的图基于的是心鹏卜部获取的娜，其 中主要感兴趣的是心外膜电位，尽管它们可以附带地包括心内膜信息。
心内膜和心外膜电子图的特征可以〗細解剖對斜己和/或各图的电Tt寺征相
互配准。配准的目的是M:两张图的电子特征之间的对应关系或相关性。 一旦
已经建立了相关性，就可以将电子变换函数应用到随后的、重复的心外膜图， 以生成心内膜图的新版本而不必重复侵入的心内膜绘图过程。
本发明的各方面使用小量的心内膜测量和相对大量的ECG通道，来产生 心内膜电位的病人特定的时空图，ECG M作为体表电位图而获取。此外，能 够估计时空的心内膜图。
例如，本发明的方法可以用于跟m过非侵入的或心包的方式进行心律不
齐的消融治疗之后。在治疗之后，重复心外膜图的获取以便确定治疗是否成功 可是足够的。在心外膜图中观察到的改变可以棚t到心内膜图上，以便^i正治 疗结果。
本发明的方法也可以用在能够以类似方式测量和绘希啲心脏的其它属性的 相关心内膜和心外膜图。
本发明的实施例提供一种用于生成活体对象的心脏电子图的方法，其通过 将探测器(如导管)插入心脏腔室中来执行。该探测器具有至少一个电极，并 可以配M第一位置传感器。电信号从来自心脏内的至少一个^M点的该至少 一个电极mt，并且在至少一个接收点接收。典型地，存在多个對寸点和接收 点。接收点可以W^的内部敬卜部。该方法iB131相对于鄉点定位接收点、 以及计算划寸的电信号和接收的电信号之间的函数关系(如线性矩阵关系)来 执行。该计算可以通过确定测量的导联场矩阵(lead field matrix)、以M31从 测量的导联场矩阵计算反转的导联场矩阵来执行，该测量的导联场矩阵定义发 射的电信号和接收的电信号之间的统性矩阵关系。可替换地，育^够确定反转的 导联场矩阵而不用明确地计算观懂的导联场矩阵。该方法CT1在相同或其它 接收点接收电生理学信号、并将反转的导联场矩阵应用到该电生理学信号以获 取心内膜电子图来执行。
在该方法的一个方面中，通过将接收点与第二位置传SI关联、并且读取 第一位置传 和第二位置传感器以确定其间的差别，来使接收点相对于发射 根据该方法的另一方面，导管具有至少两个电极，电信号从电极的不同子
集，。mt的信号可以是时分OT或频分复用的。
根据该方法的另一方面，电极是单极电极。可替换地，该电极可以是双极的。
在该方法的一个方面中，在接收点和^l寸点的子集之间确定阻抗。 在该方法的另一方面中，从在^lf点的子集中生成的电偶接收信号。 在该方法的另一方面中，相对于m的呼吸周期的预定相位，执行发射电
信号、接收发射的电信号、确定测量的导联场矩阵、以及计算反转的导联场矩阵。
在该方法的另一方面中，相对于皿的心动周期的预定相位，执行皿电 信号、接收划寸的电信号、确定测量的导联场矩阵、以及计算反转的导联场矩 阵。
该方法的一个方面包括获取胸部的解剖图像、使用该解剖图像来准备胸部 的有限元素模型、以及调整该有限元素模型的参数以使计算的导联场矩阵符合 测量的导联场矩阵。
该方法的另一方面，通过调节测量的导联场矩阵来计算反转的导联场矩 阵。通过移除反转导联场矩阵的零空间来实现规则化。
本发明的其它方法樹共用于实现该方法的装置。


为了更好的理解本发明，可参考通标例方式的本发明的详细描述，其要 结合以下附图阅读，其中相同的元件给定了相同的标号，并且在其中
图1是用于关联多个电子心脏图的系统的高级图，该电子心脏图根据本发 明公开的实施例来构造和实施；
图2是在图1示出的系统中i柳的导管的远端的平面图3是关联心内膜和心外膜电子图的系统的图示，该电子图根据本发明的 可替换的实施例来构造和操作；
图4是根据本发明公开实施例的简化的胸部咅舰图，示出躯干背心以及分 布在胸部附近的电极；
图5是根据本发明公开的实施例的示出图4中躯干背心细节的示意图6是根据本发明公开的实施例的用于关联心内膜和心外膜图的方法的流 程图7是根据本发明公开的实施例的图示关联心内膜和心外膜图的方法的各 方面的示意图8是根据本发明的公开的实施例的图示用于求解正向矩阵问题的技术
图9是根据本发明的公开的实施例的用于通过信号注入来确定导联场矩阵 的方法的流程图；以及
图10是根据本发明公开的实施例的用于开发3维病人特定的心脏电解剖 模型的功能滩图。
具体实施例方式
在下面的描述中，提出了大量具体的细节，以便彻底的理解本发明。然而， 对本领域技术人员明显的是在没有这些具体细节的情况下可实施本发明。在其 它实例中，^^f周知的电力、控制逻辑、以及传统算法和过程的计算机程序指 令的细节，没有详细地示出以便不必要地模糊本发明。
实施例1
现在转到附图，首先参照图1，其是根据本发明的公开实施例构造和操作 的系统20的图示。系统20用以确定探测器或导管22的位置、用于获職军剖 和电子数据、以及使用导管22用于组织消融。在获取心内膜电子图期间，使 用已知的血管内方式将导管22方M^ 26的心脏24的腔室中。为了获得心 外膜的电子图，导管22可以经^I^地插入包風L、脏24的心包腔。可替换地， 心外膜电子图可以非侵Ai也获得，如以下进一步详细描述的。在美国专利No. 5,471,982、 5,391,199、 6,226,542、 6,301,496和6,892,091中，以及在PCT专利 公开WO94/06349、 WO96/05768和W097/24981中描述了用于心脏绘图的示 例方法和设备，这些公开作为参考结合于此。例如，美国公开No. 5,391,199描 述了这样的导管，其包括用于领糧心脏电活动的传感器、以及用于确定导管相 对于外部施加磁场的^S的小型线圈。^ffl这种导管，會巨够通过确定在多个位 置处的电活动以及确定Ml置的空间坐标，来在短的时间段内从一乡M样的点
收集麵。
现在参照图2，其是导管22 (图l)的实施例的远端44的平面图。远端44 包括在它的顶端48的电极46，用于测量心脏组织的电属性。电极46可以是单 极的或双极的电极。电极46也可以用于将电信号发送到心脏用于诊断的目的 (例如，用于起博绘图)，和/鋼于治疗的目的(例如，用于消融有缺陷的心 脏组织)。导管22的远端44可慰也包括用于测量远场电信号的非接触电极54 的阵列52。阵列52是线性阵列，其中非接触电极54沿额端44的纵轴线性 排列。导管22的远端44还包括至少一个^a传自56，其生j^ffl于至少确定 在体内导管的位置的信号。在一些实施例中，还确定在体内导管的方向。位置 传感器56 ,地邻近顶端48，导管22的该实施例在上面提到的美国专利 No.6,892,091中进一步描述。
再次参照图l，导管22的远端44的电极和换能器iM:电缆58，经过导管 22的插入管而连接到控制处理器28 (图1),控制处理器28控制系统20的其 它元件，包括用于将信号传输到导管22的信号发生器29、图像处理器21、以 及EKG处理器27。为了方便，控制处理器28显示为单个单元。然而，它可 以实现为多个用以执行在此描述的不同的处理任务的处理设备。控制处理器28 确定导管22相对于心脏24的特定^H己或特征的^fi坐标。控制处理器28驱 动显示器40，显示器40显示身体内部的导管位置，并显示由该系统生成的功 能性图。控制处理器28还驱动消融换能器，该消融换能器一般位于导管22的 顶部。导管22用以生成解剖图或心内膜电子图。另外，导管上的各电极可以 用于消融。系统20會,是CARTO XP EP导航和消S^、统，从Biosense Webster 公司(Diamond Canyon路3333号，Diamond Bar， CA 91765，美国)可得， 适当地修改以实现以下进一步详细公开的本发明的某^f寺征。
在本发明的一些实施例中，^ffl多个体表电极31能够非侵Ai也获得心外 膜电子图，代表性ite出了三个体表电极31，该领域已知的是当使用非侵入技 术时，典型地需要更大的电极阵列以便获取精确的心外膜电子图。电极31可 以方便地安装在如任何下面的文献中描述的多电极胸面板中，这些文件的所有 内容作为参考结合于此Ransbuiy等人的美国专利申请公开No.2004/0015194; Sippensgroenewegen的美国专利申请公开No.2001/0056289;互联网上公开的 Ramanathan等人的"用于心脏电生理和心律不齐的非侵入心电图成像(Noninvasive electrocardiographic Imaging Electrophysiology and Arrhythmia) (Nature Medicine);以及上面提到的Modre等人的文献。电极31 ffijl电缆33 连接到控制处理器28，并链接到EKG处理器27。
电极31能够可替换地内部地置于对象内。例如，它们可以是暂时iik^7lc 久地在如冠状窦或心夕鹏的区域内布置的賴导联或电极。
上面提到的心包内技术能够用以产生心外膜电子图。该方制乃然比上面描 述的用于获得心内膜电子图的导管插入技术更少侵入。该技术使用心夕卜膜接触 探测器作为导管22， {顿已知的弓l入技术舰胸壁将导管22插入心刨莫。
在任一情况下，心外膜电子图典型地示出心外膜上的电位，尽管它也可以 示出心内膜的电位。无论如何，术语"心外膜电子图"在此使用作为从心脏外 部获得的首要感兴趣的 。
使用与显示器40链接的图像处理器21，基于解剖标己和/或该图的电子特 征，来将心内膜和心外膜电子图相互配准。该配准的目的是建立两幅图的电子 特征之间的变换函数，在此也称为变换式。 一旦己经建立了变换式，随后获得 的心外膜电子图能够变换为新版本的心内膜电子图。在一些实施例中，M将 心夕卜膜电子图简单地在视觉上搠寸到心内膜电子图，变换能够生效。可替换地， 数学变换可以施加至噺的心外膜电子图，以便创織版本的心内膜电子图。
在共同受让人的美国专利No. 6,650,927中公开了在一个实施例中使用的适
当的配准技术，并作为参考结合于此。简要地描述该技术
心外膜电子图和已、内膜电子图育濒是3维图。^f顿在共同受让人的申请No. 11/215,435、题目"使用生理数据的多模式图像的分割和配准(Segmentation and Registration of Multimodal Images using Physiological Data)"中公开的方法，肖g够
执《m些图的配7隹。
实施例2
现在参照图3，其是根据本发明的实施例构建和操作的系统106的图示。 系统106类似于系统20 (图l)。然而，膽26现在穿上了躯干背心108，躯 干背心108具有多个电极110，典型地大约在125到250个电极之间，这些电 极放置汩区干背心108内以W^26的背心的前面、后面和侧面樹共电位的 测量。电极110经由导线112和电缆33连接到控制处理器28。修改了控制处 理器28用于接收和处棘自躯干背心108的 。
现在控制处理器28包括用于阻抗检测的电子电路，如在2005年1月7日 提交的美国专利申请11/030,934中描述的，该申请转让给本专利申请的受让人， 并且其公开作为参考结合于此。修改该系统以基于在少量心脏内的点和电极110 之间的阻抗测量，产生其之间的函数关系。在一个实施例中，该M是线性多 维系数矩阵，在此称为导联场矩阵。然后粒矩阵的反转，例如在(YoramRudy 等人的)美国专利申请公开No. 2003/012163中描述的，其公开作为参考结合 于此。在该公开内容中，反辦巨阵对应于心外膜电位。然而在系统106中，矩 阵的反转对应于心内膜的电导，其对于现有技术是先进的。在过去，还不可能 可靠地评估外部测量和心内膜电位之间的传递函数。这是因为电场横跨心肌内 的纤维肌组织。如上面提到的，这样的组织的量和方向在个体中变化。可替换 地，在系统106的一些实施例中，导联场矩阵和它的反转可以涉及基于心外膜 的电导的图。以下，进一步详细讨论导联场矩阵的反转。
可以只《顿一个。、内的点。接收点或多个接收点旨嫩W^内部^卜部。 例如， 一个或更多，导联冠状窦电极、心外的或甚至心肌内电极能够用作各 接收点。
现在参照图4，其是根据本发明的公开实施例的显示躯干背心108的胸部 114以及分布在胸部周围的电极110的简化剖视图。图4也显示右心房116， 并且包括三个心内膜点118、 120、 122。如下面解释的，在位于心内膜点118、 120、 122的导管电极和电极110之间进行了P服测量。在一些实施例中，也在 心外放置的电极(未在图4中显示)和电极110之间测量阻抗。
现在参照图5，其是根据本发明的公开实施例的图示躯干背心108 (图3) 的细节的示意图。构憩区干背心108以包括分布^ffi力点，其可与电极110符 合。然而，这样的f袷是为了方便，并且不是基本的。该压力点124 M^软 的齿槽126连接，具有预定的自由度。齿槽126导劍区干背心108更紧密地符 合胸部114的几何开邻(图4)。躯干背心108包括至少一个^fi传自128， 该位置传自128是包括电极110的坐标系统中的参照点。参照共同受让的美 国专利申请公开No. 2004/0068178中的位置系统，教导了这样的^S传感器的 使用，其公开作为参考结合于此。iiS传感器128使得电极110的位量在医疗 过程期间能够!戯跟宗、并且舰差异计算能够被关联到心脏内的电极。传感器 128的位置不是基本的，只要电极110會的多相对于心脏内各点方爐。现在参照图6，其是根据本发明的可替换的实施例的用于关联心内膜和心
夕卜膜图的方法的流程图。在图6中描述的许多过程步骤柳醉是示例性的，并
且倉嫩改变，这将对本领,术人员来说是明显的。
在开始步骤130，膽穿±|区干背心108，并JJ^接至啦制处理器28 (图 3)。如社面提至啲美国专利申请公开No. 2003/0120163中描述的，获取了体 表心外膜电子图。
在步骤132，将导管插入心脏，并且相对于^顿躯干背心108 (图3)上的 电极110，获取了有限阻抗图。典型地，少量的点如心内膜点118、 120、 122 (图4) IM以获取该有限阻抗图。
在步骤134，获得了心脏的解剖图像。这會的多在执行了初始步骤130和步 骤132的相同时间段期间预获取或获得。实际上，如果病人在计^l几断层成像 (CT)、心肌灌注SPECT期间，穿着体表电极的"背心"，夷口么4顿已知的技术 例如超声成像，育^够可选地以接近实时地获得解剖图像的获取。然而，如上所 述获取的阻抗 可能经常足以发展出通常的身術莫型，并且还从通常的身体 模型发展出病人特定的模型。
在步骤136，应用在上面提到的美国专利申i青公开No. 2003/0120163中描 述的矩阵求解技术，在心动周期中的不同点确定了心内膜点118、 120、 122和 針电极110之间的电导，以生成卩鹏图。
在步骤138，在开始步骤130和步骤132中产生的图被组合，并与在步骤 134中获得的解剖图像配准。该步骤将体表图转换为具体的组合心内膜图。
步骤140是可选的。在一些应用中，获得关于心脏腔室的部分信息是重要 的。在步骤140中，图被分害似产生一个或更多区域图。用于分割图像并产生 部分 的图像处理技术在本领域中是^^周知的，并且ftf可适合的方法可以 应用于该步骤。
在步骤142，通常地开舰心脏的消融治疗，然后，在步骤144，通过重 复在开始步骤130和步骤132、 136、 138中描述的绘图和处理，M获得新的 详细心内膜图来銜正消融。用于评估消融的技术是已知的，并且在例如共同受 让的申请No. 11/357,512，题目为"M定步的损害评估(Lesion Assessment by pacing)"中描述，其公开作为参考结合于此。
在最后的步骤146，通过一次或更多次上面详述的过程，获得了在开始步 骤130中描述的新的体表心外膜电子图。然后，使用相同的解剖图像和在步骤
132和134中获得的有限卩鹏图，来产生了一幅顿多新的结合的心内膜亂 用于长期监视。如上面提到的，新的图能够可选地分割。虽然在准备第一电子 图时，背心的各电极通常将不与它们的位置f拾，但是不论如何该过程是有效 的，只要能够确定背心电极相对于心内点的相賴體。 一般操作
现在参照图7，其是根据本发明的公开实施例的进一步图示参照图6描述 的方法的实现细节的示意图。示出了赠148的胸部的咅舰图，穿着具有电极 152的多电极胸部面板150。心脏内导管154方^g在心脏腔室156内。导管154 具有位置传感器158和多个电极160。发生器162刺激电极160。信号在电极152 被检测并被传导至接收器164。链接至赎收器164的处理器166然后确定腔室 156和电极152之间的电导，并产生有限心内膜电导图。将意识到与电极152 的数量相比，相对少的电极160，并且因此相对少的心脏内的点被用来获得电 导或阻抗测量。
现在能够在从腔室156的心内膜表面上的点168划寸出来的信号和由电极 152接收的信号之间^lz:矩阵关系。在导管插入过程期间参照位置传感器158， 确定了点168的精确的相对位置。M反转矩阵，使用下面和在上面提到的文 献(Rudy等人的)中公开的各种技术， 一旦在电极152处的信号已知，就可 以在点168处计算腔室156内的电位。M:在心动周期和呼吸周期期间的不同 的时间以及再次在收回导管之后执行该计算，能够生成腔室156的时间变化的 心内膜电子图。通过湖糧多电极胸部面板中的信号以及将它们用作为相同调整 的矩阵的系数并重复矩阵反转、或直接将它们施加至IJ之前反转的矩阵，就能够 在将来的时间段重新生成该图。
该方法因此涉及来自导管的顶端的信息综合。首先，将测量的心内膜电位 添加到导联场矩阵(有时候称为"观懂的导联场矩阵")的各元素。其次，{顿 在导管的顶端和背心的电极之间得到的阻抗测量来更新导联场矩阵。当导管被 移动时，它的健被连续i腿跟宗，为心内膜电位和导联场矩阵积聚了更多测量。 这些测量IM以搏腿地M导联场矩阵的反转求解。
粒从腔室156的心内膜表面上的点168對寸的信号和由电极152接收的 信号之间的可靠线性矩阵关系，对获得不时旨,重新生成并与之前的实例进行
比较的适当心内膜(或心外膜)电子图是基本的。己知为"反演问题，的这个问 题在固有地存在空间不定性方面已知是复杂的。在过去，施加到心内膜电位的 数学过程趋于产生差的图像射jf率。在使用心外膜电位求解反演问题交付了合 理的结果时，应用到心内膜会导致涂抹的图像。该问题还被在病人中变化的心 脏肌肉的纤维结构而一一步复杂。处理这个需要张量阻抗模型，其不能容易地
从MRI或CT扫描估计。 正演问题
"正演问题"魏于给定的电源分布，在具有给定的电的和介电的特性(传 导率)的介质中发现电位的过程。该问题导致具有唯一解的统性矩阵方程
<formula>formula see original document page 15</formula> (i)
其中A是传递矩阵(导联场矩阵)，X是电流源或心内膜或心夕卜膜电位或跨膜
电位，以及5是在体表背心湖糧的电位的阵列，即体表电位图。 基于模型的方法
在一种方法中，特征化了普通的人体组织传导率和MR/CT模型的分割。 使用有限元素方法(FEM)解算器计算了导联场矢巨阵。机械或电子模型具有一 些自由度5tt示影响矩阵A的关键因素G)，艮口， <formula>formula see original document page 15</formula>(2)
关键因素(？)能够是器官的几何位置或大小，以^H察域内的各种组织 的比例。器官和组织的传导率在个体病人中不同。例如，光纤方向影响心肌中 的传导率张量的方向。通过特征化在导联场矩阵中反映的这些差异来改进并完 成病人特定的1iM。
通过注入的信号完成的测量被用以最优化关键因素(？)，使得FEM角軒去 变为
<formula>formula see original document page 15</formula> (3)
其中力是注入到导管中的电流，该导管放置在位置？处(在参照坐标系中，
在该时间接地)。；。)是点源、—的由FEM解算器计算的阻抗矩阵。对于一
组机械的或电子的模型参数G)，值&是指示接收位置组中所观糧的电压的
向量。测量可以是阻抗测量。可替换地，它们可以是基于来自电偶的各信号的
导联场矩阵的测量。这样的偶极可以通过在导管中M5:两个相邻的电极之间的
差别来产生。如果各电t腿当地排列，那么偶极會,在三个正交方向中定向。
由方程4给出了用于观懂的卩鹏("的最优化。 》=argM"Zl2訓丽广^^她,(到l2 (4).
用于该组参数(？)的选择包括器官大小、传导率、光纤方向以及各向异性比 率。已知的最优化搜索算法能够用以确定用于该组参数(？)的值，例如遗传 算法、仿真退火、以及神经网络或其混合。在D.Farina、 O.Skipa、 C.Kaltwasser、 0.D6ssel和W.RBauer的"基于最优化的心脏去极化重建(Optimazation-based reconstruction of depolarisation of the heart) ，， (Proc. Computer in Cardiology,芝 加哥，美国，2004， 31， 129-132)中公开了这样的最优化的示例。
当施加到导联场矩阵时，偶极产生的电^到与单电压源相同的作用。在 一些情况下，将心脏电活动描述为一组在心ULh传播的偶极是更方便的，在该 心肌中，密度和方向参数随着心动周期变化。
使用该组参数G)，會嫩校准FEM牛IM，以便精确地，病人的解剖结 构。从心脏内部产生信号给出关于心肌的属性的重要信息，即关于光纤方向的 假设的正确性。该信号也能够展现其它的诊断信息。缺血的、疤痕和扩张相关 (stent-related)组织可以具有原理正常传导率的显著偏差。该方法的增加的好 处是心肌阻抗的视觉图，其自身可以具有诊断价值。另外地或可替换地，导管 可以心外放置，并且信号注入的方法施加到心外导联。当以这种方式执行该过 程时，可以基于外部的读数，来生成心脏的心外膜的电子图。以下的讨论关注 于心内膜图生成，但已作必要的修正，可应用到心外膜图生成。如果期望，在
此公开的原理，可以应用到M:将涉及心夕卜膜电位的矩阵转换到心内膜信号。
可替换地，电流能够通过两个相邻的心内电极注入，因此模拟电流偶极。 以这种方式，导联场矩阵A的一磐巨阵元素f詢多直接测量。该过程肯嫩用以更 新原始的导联场矩阵，或它育嫩用以构建导联场矩阵A而不用倒可MRI或CT 组。在任一情况下，4顿这样的信号产生了体表电位图(BSPM)，该信号 4OT注入的电流而不使用从生物电源出现的电流偶极。
对反演问题的数学求角科皮呼吸劣化。这倉,Mil门f空在呼吸周期中的一个 点(例如结束呼气)处的数据收集而最小化。但是所期望的是，齡呼吸周期
中的数据，以及构建依赖呼吸的导联场矩阵，对于呼吸周期的相位来校正该导 联场矩阵。当完成这个时对反演问题的求解的质量改进相当多。
发明人已经公开了如何克服矩阵反转中固有的技术困难。M:将不同组合 的信号典型地注入到位于心脏腔室中的相对少的源中，并且测量在接收点的信 号，而能够一列一列精确地确定导联场矩阵。使用与矩阵中存在的未知系数一 样多的心脏内电极的不同组合，来重复该计算。当然，系数的数量依赖于心脏 内的源和外部的导联的数量。
如果4顿了多电极导管，例如PENTARAY，高密度绘图导管(从Biosense Webster公司可得)，当导管固定时，在许多点测量了心内膜电位。这允许躯干 背心测量快速地完成。
现在参照图8，其是根据本发明的公开实施例的图示用于展开正向矩阵的 技术的图示。在图8的上部中的矩阵表格182中图示了方程1。 ^f顿心脏184 的腔室内的两个或更多的源例如电极160 (图7)，注入了电子信号。如图8的 上部中图示的，该信号表示为向量186,其对应于图的上部中的向量188。表 示为向量190、并X寸应于图的上部中的向量192的多个卩m测量，在各源和身区 干194的多条外部导駒列如电极152 (图7)之间得到。这允许粒导联场矩 阵198的一列196。在图8的示例中，列196是导联场矩阵198的^&列200。 通过重复地改变源信号的位置，或在多元件绘图导管的情况中变化源信号的组 合，而创建并求解了矩阵方程。同时，以极大的精度确定了导联场矩阵198的 系数，噪声由相对大幅度的SA信号淹没。
现在参照图9，其是根据本发明的公开实施例的图示战用于M信号注 入来确定观懂的导联场矩阵的方法流程图。在初始步骤202，通过应用躯干背 心、和插入链接到位置传繊例如系统106 (图3)的心脏绘图导管，来使病 人做准备。期望绘图导管具有多个电极。当會嫩4OT足够的源时，当不用在心 脏内额外导航导管时这加速该过程。需要得到以便确定导联场矩阵的领糧的数 量被记录。
接下来，在步骤204，选择了至少两个心脏内源导联的组合。借助系统106 的位置处理工具，精确得知这些源相对于心鹏賴啲参照特征的j體。
如对本领嫩术人员明显的是，在步骤204中，也可以{顿单极导联，同 时在导联场矩阵的计算中进行适当调整。
接下来在步骤206， j顿电流源插入电子信号以创建电偶。用于该信号的 合适的值在l-100kHz为l-10mA。该信号可以是疸电压或恒电流信号。
接下来在步骤208，记录了在步骤206中选择的源导联和^^躯干背心的 导联之间的阻抗测量。可以使用共同受让的美国专利申请公开 No.2007/0060832，题目为"舰阻抗检测(Detection of skin Impedance)"的技 术来测量阻抗，在ltkil过弓间并入该申请公开。在i顿电偶的实施例中，在躯 干背心导联确定了偶极位置和方向。
现在控制进行到判定步骤210，其中确定是否需要更多的测量。如果在决 定步骤210的确定是肯定的，那么控帝腿回到步骤204 ^i^择另一个源。
如果在决定步骤210的确定是否定的，另P么控制进行到最后步骤212。矩 阵方程被求解并且导联场矩阵值l^艮告。
如上面提到的，图8和图9中描述的测量相对于呼吸周期被门控。此外， 它们相对于心动周期被门控。通过在周期中的不同点重复观糧，能够以倒可期 望的时空分辨率，来获得时间变化的病人特定的导联场矩阵组。
将从有限元素模型的考虑回想起参数最优化是必须的。在执行最后的步骤 212之后，并且一旦可靠的导联场矩阵可用，就能够应用最优化算纟妹在方程 3中建立符合实际导联场矩阵的参数G)，同时大大地减少了最终的病人特定 的有限元素模型中的模型错误。然后有更大的信心来能够展开反演问题的求 解。最后，肯的多从少量最初的心脏内测量产生描绘病人的心脏中的电生理学的 接近完美的图像，而没有由呼吸周期引起的伪像。會嫩不时地重复这些"4D"图 像例如来评估治疗。应当注意到，在随后的时间段中的接收点不需要与原始时 间段中的接收点相同，只有例如通过参照坐标系，它们相对于原始传输点或原 始接收点的相对位置是可识另啲。将反转的导联场矩阵应用到新的接收点保持 有效。
实施例3<formula>formula see original document page 18</formula>
(5)
方程5中的符号与方程4相同。由方程4描述的基于模型的最优化技术， 现在直接应用来确定反转的导联场矩阵；i ，而不用明确地计算导联场矩阵自身。
实施例4
现在参照图10，其是根据本发明的公开实施例的用于开发3维病人特定 的心脏电解剖模型的方法的功能方框图。该功能图能够舰规范化系统106 (图 3)并应用上面参照实施例1、 2、 3公开的技7(^，、以及^顿参照图9描 述的方^^实现。
最初，在功能组216中使用传统的3维解剖成像模态，来准备3维解剖的 病人模型214，这与步骤134 (图6)类似，但通常延伸到心脏自身之外的胸部 器官和组织。功能组218展开了在具有对象的身体中具有已知的坐标的电子信 号的矩阵，其在功能块219中与模型214合并。
在功能组220中，如在实施例l中描述的在心脏导管插入期间产生了电子 解剖图。该图可以j顿CARTO XP EP导航和消融系统(从Biosense Webster 公司，3333 Diamond Canyon Road Diamond Bar, CA 91765可得)获得。
在导管插入期间，在功能组222中应用了参照图9描述的方法。在功能块 224中计算了观糧阻抗阵列。该P车列在块226中应用。该阵列用以舰功能i央 226的模型214，因此影响模型228。模型228然后l細以在功育旨块230中求解 反演问题和最优化FEM参数(方程3、 4、 5)。典型地，以下描述的反演问题 可以在经由功能组222的数据的影响之前，最初在块230中求解，并且在功能 块232中产生了最初的电子解剖图像。在最初的求解中，在功能块232中可以 开发初步的心内膜电位以及可选的心外膜电位的电子解剖图像组。随后，在包 括功能组222的好处后，在功能块234中生成了舰的图像。由功能组222提 供的数据的另一效果是，在功能±央236中计算导联场矩阵。这可以重复地与在 功能块224中观糧的矩阵比较，并应用来4顿用来求解方程3、 4、 5的最优化 算絲鹏丽。
图10中描述的系统操作會,M确定呼吸参数《t)，用以确定在心动周期 cc(t)期间导联场矩阵对运动的依赖。该技术最终加速了导管插入过程，改进了 心律不齐和f赎形成的诊断，并改进了介入治疗的结果。
反演问题
{顿图6和图7中描述的方法和系统的主要目的在于，基于少量心脏内电 位和许多体表测量计算心肌中的电源，其已知为"反演问题'。在方程1中，测 量的信号S是己知的。如上戶，，导联场矩阵A可以M使用病人的MRI或 CT扫描求解几个正演计算来计算。然后描述源的向量；从下面的方程确定
A一1 . S = S (6)
数学地，发现导联场矩阵A的反转A-'是必要的。不幸的是，该问题是不 适定(ill-posed)的。矩阵A具有非零的零空间，艮卩能够发现不同的向量S， 其属于该零空间并导致小于噪声的电位S 。向量^和零空间外的向量的每个线 性组合导致相同的电位^ ,并且因此是方程1的解。此外，由于未知的个体的 传导率，矩阵A的系数具有一些不确定。
规则化的求解
规则化是广泛描述的用于反演问题的方法，其中丢弃了反转的零空间。最 通常使用的方法是具有为零的Tikhonov新因子的Tikhonov规则化，其选择具 有最小范数(最小可能的源)的解。
当根据，各方法确定导联场效应时，有两种主要因子改进反转的质量。 第一，由于注入信号的幅度而改进了信噪比。ECG跨膜电位在10mV的数量 级并且持续时间相对短，而生成的信号在伏的数量级并且持续时间长的多。因 此，平均技术會,用以进一步提高信噪比。幅度的两个数織的鹏容易实现。 这允许较少能量的特征值有意义地参与到解中，并且最后改进观察的心内膜电 位的可靠性。第二，当^顿更大量的心脏内的源时，舰了矩阵的质量。为此， 4顿多元件绘图导管是期望的，以及f顿具有大量接收电极糊区干背心。
本领域技术人员将意识到本发明不限于以上己经特定示出和描述的。相 反，本发明的范围包括不在现有技术中的、以上描述的各种特征的组合和子组 合及其变化和修改，其可由本领Jl^术人员fflil阅读之前的描述想到。
权利要求
1.一种用于生成活体对象的心脏的电子图的方法，包括以下各步骤将探测器插入心脏的腔室中，所述探测器具有至少一个电极；从来自心脏内的至少一个发射点的所述电极发射电信号；在至少一个接收点接收所述发射的电信号；相对于所述至少一个发射点定位所述接收点；确定所述发射的电信号和所述接收的电信号之间的函数关系；在新的接收点接收电生理学信号；以及将所述函数关系应用到所述电生理学信号，以获得心内膜电子图。
2. 如权权利要求的方法，其中所述函数关系是测量的反转导联场矩阵。
3. 如权禾腰求2戶脱的方法，还包括以下步骤 获取ffWm的胸部的解剖图像；4OT所述解剖图像准备所述具有各参数的胸部的有限元素模型，所述有限 元素模型具有计算的反转导联场矩阵；以及调整所述各参数以使戶脱计算的反转导联场矩阵符合所述测量的反转导联 场矩阵。
4. 如权利要求1所述的方法，还包括步骤在执行戶，接收电生理学信号和应用戶腿函数关系的步mt前，所述对象撤回所述探测器。
5. 如权利要求1所述的方法，其中所述的少一个接收点^0Mxm外部。
6. 如权利要求1所述的方法，其中所述的至少一个接收点在戶; 内部。
7. 如权利要求1所述的方法，其中所述的探测器具有至少两个电极，并且鹏{顿戶脱电极的不同子集时分飾戶腿电信号，来执行激寸电信号。
8. 如权利要求1所述的方法，其中所述探测器具有至少两个电极，并且通过〗OT戶，电极的不同子集频分OT戶;M电信号，来执行Mt电信号。
9. 如权利要求l所述的方法，其中戶，电极是单极电极。
10. 如权利要求l所述的方法，其中所述电极极电极。
11. 如权利要求i所述的方法，其中相对于所述的呼吸周期的预定相 {姊执行戶诚魁寸电信号、接脱鄉的电信号、以及确定函数l系的步骤。
12. 如权利要求i所述的方法，其中相对于戶;f^f象的心动周期的预定相 位来执行戶腿鄉电信号、接柳腿鄉的电信号、以及确定函数关系的步骤。
13. —种用于产生活体对象的心脏的电子图的方法，包括以下步骤将导管插入心脏的腔室中，所述导管具有第一定位传感器和至少一个电极；从在心脏内的多个，点的所述电极发射电信号； 在所述对象外部的多个接收点接收戶舰激寸的电信号；相对于所述发射点定位所述接收点；确定测量的导联场矩阵，以定义所述发射的电信号和所述接收的电信号之间的线性矩阵关系；根据所述测量的导联场矩阵来计算反转的导联场矩阵； 在所述接收点接收电生理学信号；以及 将所述反转的导联场矩阵应用到所述电生理学信号，以获取心内膜电子图。
14. 如权利要求所述腿的方法，其中所述定位所述接受点的步骤包括关联所述接收点与第二位置传感器；以及读取所述第一位置传感器合所述第二位置传感器，以切丁之间的差别。
15. 如权利要求13所述的方法，其中戶诚导管具有至少两个电极，并且用所述电极的不同子集执行发射电信号。
16. 如权利要求13所述的方法，其中所述电极是单极电极。
17. 如权利要求13所述的方法，其中所述电极单极电极。
18. 如权利要求13所述的方法，其中ffl31确定所述接收点和所述激寸点的子集之间的阻抗，来执行所述接受所述发射的电信号的步骤。
19. 如权利要求13所述的方法，其中通过测量由在戶舰发射点的子集中生成的电偶产生的信号，来执行所述接受所述发射信号的步骤。
20. 如权利要求1 所述的方法，其中相对于所述对象的呼吸周期的预定相位来执行发射电信号，接受所述发射的电信号，确定测量的导联场矩阵、 以及计算反转的导联场矩阵的步骤。
21. 如权利要求13所述的方法，其中相对于所述的心动周期的预定相位来执行所述发射电信号、接收所述发射电信号、确定测量的导联场矩阵、 以及计算反转的导联场矩阵的步骤。
22. 如权利要求13所述的方法，还包括以下步骤获取戶;iwm的胸部的解剖图像；使用所述解剖图像准备所述具有各参数的胸部的有限元素模型，0M有限元素模型具有计算的导联场矩阵；以及调整所述各参数以〗妙脱计算的导联场矩阵f始戶舰测量的反转导联场矩阵。
23. 如权利要求13所述的方法，其中所述计算反转的导联场矩阵的步骤 包括通过移除所述反转导联场矩阵的零空间来规则化所述测量的导联场矩阵。
24. —种用于成像活体膽中的心脏的系统，包括 成像设备；信号发生器；以及链接到穿在所^m上糊区干背心的处理器，所淑区干背心包括多个接收 器和第一位置传感器，所述处理器链接到所述成像设备、所述信号发生器以及 链接到适于插入戶，心脏的绘图导管，所述绘图导管具有绘图电极，所述处理 器操作用于读取第一位置传繊，以相对于戶服接收器定^0M绘图电极，所 述处理器操作用于{妙脱信号发生器将电信号顺序地发送到戶舰绘图电极，并 且所述绘图电极从所述心脏中的不同鄉点",地发生电信号，其中所述发射 的电信号经由所述接收器传递到所述处理器，以作为接收的电信号，所述处理 器还操作用于确定湖糧的导联场矩阵，该观懂的导联场矩阵定义在其各自的位 置处的所述對t的电信号和所述接收的电信号之间的线性矩阵关系，根据戶脱 观糧的导联场矩阵来计算反转的导联场矩阵，并应用戶舰反转的导联场矩阵来处理在戶，接收器接收的戶; 的电生理学信号，以根据所述电生理学信号来生成心内膜电子图，并且将戶，心内膜电子图显示在戶皿成像设备上。
25. 如权利要求24所述的系统，其中戶脱绘图导管具有至少两个绘图电极，并且戶腿信号发生器操作用于m将戶腿电信号顺序地传送给戶舰绘图电极的不同子集，来产^0M^I寸的电信号。
26. 如权利要求24所述的系统，处理器操作用于在所述对象的呼吸周期的预定相位处确定戶;M测量的导联场矩阵。
27. 如权利要求24戶腿的系统，处理器操作用于在戶;MX^的心动周期的预定相位处确定戶;M测量的导联场矩阵。
28. 如权利要求24所述的系统，其中所述处理器操作用于接^BW^ 的胸部的解剖图像，并且^ffl所述解剖图像来准备具有各参数的所述胸部的有 限元素模型，戶脱有限元素模型具有计算的导联场矩阵，以及调^)5脱各参数以^0M计算的导联场矩阵f始戶;M测量的反转导联场矩阵。
29. 如权利要求24戶服的系统，其中戶服处理器操作用于通过移除戶做 反转导联场矩阵的零空间规则化所述测量的导联场矩阵，来计算戶腿反转的导 联场矩阵。
全文摘要
通过构建少量心内膜点和使用多电极胸部面板的外部接收点之间的矩阵关系，来获得了一种可靠的心内膜图。矩阵的反转产生允许心内膜图被构建的信息。使用多电极胸部面板、将新的电信号应用到矩阵关系、以及再次反转矩阵以产生新的心内膜电子图，从而非侵入地获得随后的图。
文档编号A61B5/0402GK101199416SQ200710182138
公开日2008年6月18日 申请日期2007年9月6日 优先权日2006年9月6日
发明者J·波拉思, M·巴-塔尔, Y·施瓦茨 申请人:韦伯斯特生物官能公司