CN103892824A - 除去标测图数据中的赝象 - Google Patents

Abstract

一种除去标测图数据中的赝象。本发明公开了一种用于标测的方法，所述方法包括：接收生理参数的初始的一组测量值，所述生理参数的测量值在身体器官中的相应位置处测量，以及接收包括一系列空间标测图元素的、所述器官的三维(3D)标测图。所述方法包括在测量所述测量值处的所述相应位置与所述标测图元素的子组之间形成对应关系，并且响应于所述对应关系将所述生理参数的相应元素值与除所述子组之外的标测图元素进行关联。所述方法还包括调整所述相应元素值，使得相邻各组的所述标测图元素形成短程线，以及显示所述器官的标测图，所述标测图示出了所调整的元素值。

Description

除去标测图数据中的赝象

技术领域

本发明整体涉及物理参数的测量，并具体地讲涉及与身体器官例如心脏相关联的参数的测量。

背景技术

在医疗手术(例如标测心脏的电活动)中，测量通常相对嘈杂。用以降低测量过程中的噪声影响的系统将是有益的。

发明内容

本发明实施例提供了一种用于标测的方法，所述方法包括：

接收生理参数的初始的一组测量值，所述测量值在身体器官中的相应位置处测量；

接收包括一系列空间标测图元素的、器官的三维(3D)标测图；

在测量所述测量值处的相应位置与标测图元素的子组之间形成对应关系；

响应于所述对应关系将生理参数的相应元素值与除所述子组之外的标测图元素进行关联；

调整相应元素值，使得相邻各组的标测图元素形成短程线；以及

显示器官的标测图，该标测图示出了调整过的元素值。

在一个实施例中，身体器官包括心脏，并且生理参数包括心脏的局部活化作用时间。

通常，空间标测图元素包括具有顶点的平面多边形，所述顶点对应于身体器官的壁上的位置。在本发明所公开的实施例中，平面多边形包括三角形，并且所述子组包括相应位置的、相应最近的三角形。除子组之外的标测图元素可包括一个或多个相邻三角形，所述一个或多个相邻三角形与相应最近的三角形中的每一者具有公共边。

该方法还可包括响应于与相应最近的三角形中的每一者相关联的慢度矢量而识别一个或多个相邻三角形。

调整相应元素值可包括调整给定的最近的三角形的质心的元素值和邻近给定的最近的三角形的一个或多个三角形的质心的元素值，以形成短程线。作为另外一种选择或除此之外，调整相应元素值可包括最小化与质心相关联的局部位移矢量。

在可供选择的实施例中，除子组之外的标测图元素包括一个或多个三角形，所述一个或多个三角形不与相应最近的三角形中的每一者具有公共边。

在另外的可供选择的实施例中，短程线包括空间短程线。

在另外的可供选择的实施例中，显示标测图包括将与器官相关联的等时线结合到标测图中。

在另外所公开的实施例中，短程线包括时间短程线。

接收初始的一组测量值可包括从相应位置处的多个相应电极同时接收生理参数的多个测量值。

根据本发明的一个实施例，还提供了用于标测的设备，所述设备包括：

探针，该探针被配置成：

生成生理参数的初始的一组测量值，所述测量值在身体器官中的相应位置处测量，以及

生成包括一系列空间标测图元素的、器官的三维(3D)标测图；以及

处理器，所述处理器被配置成：

在测量所述测量值处的相应位置与标测图元素的子组之间形成对应关系，

响应于所述对应关系将生理参数的相应元素值与除所述子组之外的标测图元素进行关联，

调整相应元素值，使得相邻各组的标测图元素形成短程线，以及

显示器官的标测图，该标测图示出了调整过的元素值。

结合附图，通过以下对本发明实施例的详细说明将更全面地理解本发明。

附图说明

图1为根据本发明实施例的电生理信号分析系统的示意图；

图2为根据本发明实施例的示出心脏内壁的三维标测图的示意图；

图3A和3B为根据本发明实施例的图2的标测图中三个三角形的示意图；

图4为根据本发明实施例的在计算心脏中的估计局部活化作用时间的过程中，由处理器所执行的步骤的流程图；并且

图5示出了根据本发明实施例的等时线的示意图。

具体实施方式

综述

本发明的实施例提供用于标测与身体器官相关联的物理参数的方法。通常，如本文所假设，身体器官是受试者的心脏，并且物理参数包括与心脏跳动相关联的局部活化作用时间(LAT)。LAT是通过心脏壁的电活动流的指示，并且本发明的实施例使用初始的一组LAT测量值。

在心脏壁的初始的一组LAT与空间标测图元素(通常为网格形式的多边形元素)的子组之间形成对应关系。除非另外指明，否则在以下描述中，假设多边形空间标测图元素包括三角形元素。包括子组在内的空间标测图元素可从心脏壁的测量位置生成。

一旦已执行对应关系，则估计与除子组中的那些之外的标测图元素(即，除子组中的三角形之外的三角形)相关联的LAT。可对估计的LAT值进行调整，使得与LAT相关联的三角形(通常包括子组中所包括的三角形)的质心形成短程线。短程线通常为空间短程线，其中距离是最小的。在一些实施例中，短程线包括时间短程线，在这种情况下使质心的LAT之间的测量时间最小化。估计过程通常形成多条短程线。

一旦LAT已被调整以形成多条短程线，则可将LAT排序以生成等时线。可显示出通常基于上述网格的心脏标测图，其示出了等时线形式的调整过的LAT。

调整LAT使得其形成心脏壁的标测图上的短程线，这使得能够相对于心脏壁的部分的位置对LAT进行平滑处理。本发明人已发现与通过现有技术方法的平滑处理相比，此类平滑处理产生更好的结果。

系统描述

现在参见图1，其为根据本发明实施例的电生理信号分析系统20的示意图。系统20可被配置成分析基本上任何生理参数或此类参数的组合，但在本文的描述中，以举例的方式，假设分析的信号为心内心电图(ECG)电势-时间关系。为了充分表征此类关系，信号需要及时彼此参考。

在本发明的实施例中，时间参考通过针对参考信号及时测量例子(本文称为参考例子)而实现。在本文中，以举例的方式，假设参考信号包括参考ECG电势与时间信号的关系。又如，假设参考例子是ECG参考信号的QRS复合波的起始。对于被标测的心脏中的任何给定位置，被称为位置的电活动的局部活化作用时间(LAT)的物理参数可依据满足预定条件的电活动来限定。

在以下描述中，假设预定条件包括该位置处的电势的最大快速变化发生的时间，并且假设LAT为从参考例子至该位置处的最大快速电势偏转开始的时间。LAT可以是正的或负的。用于确定最大快速电势偏转发生的时间的方法，以及确定LAT的其他定义和条件是本领域的技术人员所熟知的，并且假设所有此类方法、定义和条件均包括在本发明的范围内。

为简洁和清晰起见，除非另外指明，否则下述描述均假设是其中系统20使用探针24感测来自心脏34的电信号的研究过程。假设探针的远端32具有电极22。本领域的普通技术人员将能够修改对于具有一个或多个电极的多个探针，以及对于由除心脏之外的器官产生的信号的描述。

通常，探针24包括导管，所述导管在系统20的用户28执行标测过程期间插入到受试者26的体内。在本文的描述中，以举例的方式假设用户28为医疗专业人员。在该过程中，假设受试者26附接到接地电极23。此外，假设电极29在心脏34的区域内附接到受试者26的皮肤。

系统20可由系统处理器40控制，所述系统处理器包括与存储器44连通的处理单元42。处理器40通常安装在控制台46中，所述控制台包括操作控件38，所述操作控件通常包括专业人员28用来与处理器互动的指示装置39，例如鼠标或轨迹球。处理器使用存储在存储器44中的软件，包括探针跟踪器模块30和ECG模块36来操作系统20。将由处理器40执行的操作的结果在显示器48上呈现给专业人员，所述显示器48通常向用户呈现图形用户界面、由电极22感测的ECG信号的视觉表示和/或正在被研究的心脏34的图像或标测图。例如，该软件可以电子形式通过网络下载到处理器40，或者作为另外一种选择或除此之外，该软件可被提供和/或存储在非临时性有形介质(例如，磁性存储器、光学存储器或电子存储器)上。

联接ECG模块36，以便接收来自电极22和电极29的电信号。模块被配置用于分析所述信号并且可将分析结果以标准ECG格式(通常为随时间移动的图形表示)呈现在显示器48上。

探针跟踪模块30在探针位于受试者26体内时跟踪探针24的多个部分。跟踪模块通常在受试者26的心脏内跟踪探针24的远端32的位置和取向。在一些实施例中，模块30跟踪探针的其他部分。跟踪器模块可使用本领域中已知的用于跟踪探针的任何方法。例如，模块30可在受试者附近操作磁场发射器，使得来自发射器的磁场与被跟踪的探针部分中的跟踪线图相互作用。线圈与磁场相互作用生成信号，所述信号被传输到模块，而且模块对所述信号进行分析，以确定线圈的位置和取向。(为简洁起见，图1中并未示出此类线圈和发射器。)由Biosense Webster，Diamond Bar，(CA)生产的系统使用此类跟踪方法。作为另外一种选择或除此之外，跟踪器模块30可通过以下方式跟踪探针24：测量电极23、电极29与电极22之间的阻抗，以及可位于探针上的其他电极的阻抗。(在这种情况下，电极22和/或电极29可提供ECG以及跟踪信号。)由Biosense Webster生产的

系统采用磁场发射器以及阻抗测量两者以用于跟踪。

通过使用跟踪器模块30，处理器40能够测量远端32的位置。此外，通过使用跟踪器模块30和ECG模块36二者，处理器能够测量远端的位置，以及在这些具体的位置处检测到的信号的LAT。为清楚起见，在本公开和权利要求中，不具有相关的LAT测量值的远端的测量位置在本文中被称为非LAT位置，并且具有相应的LAT测量值的远端的测量位置被称为LAT位置。在本发明的实施例中，假设非LAT位置被用于生成心脏34的壁的三维(3D)解剖标测图。

图2为根据本发明实施例的示出心脏34的内壁的3D标测图62，以及LAT位置的示意图60。示意图60可呈现在显示器48上。为简洁起见，完整标测图的仅一部分显示在图2中。标测图62被制定为包括多个非LAT位置点64的网格，所述多个非LAT位置点64的位置已由跟踪器模块30进行了评估。心脏壁是运动的，但在评估非LAT位置点的位置的过程中，所述模块例如通过在心脏跳动期间诸如心房收缩开始时将所有测量点调整到参考时间而允许此类运动。通过本领域已知的方法，处理器40通过点问直线66连接点64，以形成连接的平面多边形的网格。所述平面多边形可具有任何便利数量的边，并且例如可包括五边形或六边形。为简洁起见，在以下的说明中，假设连接的平面多边形包括三角形70，并且本领域的普通技术人员将能够修改对具有除三条边之外的边的平面多边形的情况的说明。连接的三角形70形成近似于心脏内壁表面的表面。

示意图还显示出LAT位置68，每个LAT位置具有相关的LAT。通常，对LAT位置及其相关的LAT在某一时间段处进行评估，所述时间段不同于处理器40用以生成标测图62的时间段。(在显示器48上呈现示意图60的过程中，与给定LAT位置相关联的LAT的值可通过对代表LAT位置的点进行颜色编码而表示。)在本公开中，根据需要，特定的非LAT位置64、线66和LAT位置68通过将参考字母作为后缀和/或前缀添加到标识标号来加以区分。例如，在示意图60中，三个非LAT位置64D，64E，64F形成三角形70D，而LAT位置68D靠近三角形70D，但与其分离。

就非LAT位置而言，将LAT位置调整至参考时间。原则上，LAT位置68应与三角形70的表面配准，因为两种类型的位置(即，LAT位置和非LAT位置)应当位于心脏壁上。然而，在实施过程中，由于例如位置的测量过程中存在误差，以及调整测量位置的过程中存在误差，因此所述位置没有配准。误差通常至少部分地由于心脏的运动。本发明的实施例纠正了两种类型位置的失配。

图3A和3B为根据本发明实施例的标测图62中三个三角形70A，70B，70C的示意图。该示意图结合一组正交的xyz轴绘制，使得图3B为三个三角形的侧视图，并且图3A为顶视图。三角形70A具有作为顶点的非LAT位置64L64M，64N，所述顶点通过线66L，66M，66N连接。三角形70B具有作为顶点的非LAT位置64P，64M，64N，所述顶点通过线66Q，66M，66P连接。三角形70C具有作为顶点的非LAT位置64L，64M，64Q，所述顶点通过线66L，66Q，66S连接。顶点64M，64N和线66M是三角形70A，70B公共的，并且顶点64M，64L和线66L是三角形70A，70C公共的。以举例的方式，假设xyz正交轴具有平行于线66M的z轴，并且被配置为使得三角形70A位于平行于xz平面的平面72中。LAT位置68B靠近三角形70A，并且不在平面72中。

三角形70A具有几何质心C70A，三角形70B具有几何质心C70B，三角形70C具有几何质心C70C。应当理解，质心C70A，C70B，C70C可由相应三角形70A，70B，70C的顶点的已知值计算。图3A和3B的其他元素在下文有所描述。

心脏的电活动可以被看做是在每次心脏跳动开始时在窦房结处发起，且流经构成心脏的心肌和结缔组织的电势。在心脏的腔壁上的任何点处，该点处的LAT由经过该点的电势流引起。如上所述，网格62近似于心脏的壁。如下所述，本发明的实施例使用测量的LAT位置及其相关的LAT值，连同测量的非LAT位置一起，通过对经过网格62的电活动流以及网格的三角形质心处的LAT进行估计，而生成心脏腔壁上的点处的LAT估计值。

图4为根据本发明实施例的在计算心脏34中的估计LAT的过程中由处理器40所执行步骤的流程图100。

在标测步骤102中，处理器和跟踪器模块30接收和采集非LAT位置64的3D值。所述接收和采集可通过移动探针24的远端32直到该远端接触心脏壁而实现。

在生成网格的步骤104中，处理器使用直线连接非LAT位置64，以便形成网格62。网格62被形成为构成三角形70的一些列空间标测图元素。由一组3D点形成三角形元素网格的方法是本领域已知的，并通常包括将任何给定的点连接至其最近相邻点中的一个或多个。

在质心步骤106中，处理器计算每个网格三角形的相应质心的值。

在“原始”LAT步骤108中，处理器40使用跟踪器模块30和ECG模块36接收和采集各组LAT位置68及其相关的LAT。所述采集基本上类似于对非LAT位置64的采集(步骤102)，但通常远端32在每个LAT位置68处的驻留时间长于在非LAT位置处的该时间，以允许处理器采集所述位置的LAT。

步骤102和108彼此独立。因此，所述步骤可相继执行，或者它们可基本上同时执行。

在流程图100的剩余步骤中，处理器对网格62的基本上所有三角形70，以及基本上所有LAT位置68执行计算。为清楚起见，对步骤的解释是指图3A和3B的三角形70A，70B和70C。

在投射步骤110中，处理器确定最靠近的三角形70，即对应于每个LAT位置68的三角形70。此类三角形在本文中称为基础三角形。一旦已确定了给定LAT位置68的基础三角形，处理器便将LAT位置投射到LAT投射点处的基础三角形上。因此，参见图3A，处理器40确定三角形70A最靠近LAT位置68B，所以三角形70A是基础三角形。处理器将LAT位置68B投射到基础三角形70A上的LAT投射点P68B。

在相邻三角形步骤111中，处理器确定与位于步骤110中的基础三角形邻近(即，相邻)的三角形。

在γ和LAT分配步骤112中，对于网格中的每个基础三角形元素，即具有LAT投射点的每个三角形，处理器根据公式(1)分配标测图元素值(本文假设所述标测图元素值包括质心的局部活化作用时间)：

t_i-γ_ki·τ_k=ε_t    (1)

其中t为基础三角形的质心的已分配局部活化作用时间(LAT)，

i为基础三角形的标识符，

τ为LAT位置的LAT，

k为LAT位置的标识符，

ε_t为常数，其为三角形质心的所测量的LAT值(τ_k)与期望的LAT值(t_i)的相似性的量度。ε_t可被预设为较小的随机值，通常在±(～0.01-0.05)ms的范围内，并且

γki为参数。

该公式将在步骤119中使用。

参数γki通常在0-1的范围内，并且参数的值可如下文中所说明通过处理器40进行调整。通常，将参数的值设置为随着LAT位置与质心之间的距离减少而更接近1。

在一个实施例中，可使用指数函数根据公式(1a)制定参数γki的预设值：

γ_k，i=e^{-α·d(i，k)}

              (1a)

其中α为度量常数，并且

d(i，k)为第i个质心与第k个LAT位置之间的距离。

对于图3A和3B中示出的三角形，三角形70A是基础三角形，并且公式(1)变为：

t_70A=γ_68B,70B·τ_68B    (1b)

公式(1)可被视为衍生自通过探针24上的单个电极22生成的测量值。更普遍的情况(其中包括公式(1))是其中LAT通过探针24上的M多个电极22同时测量的情况，其中M为正整数。对于这种普遍情况，适用公式(1c)：

$t_{i+m}-{\mathrm{\γ}}_{\mathrm{ki}}^{\mathrm{\Σ}}\·({\mathrm{\τ}}_k+{\mathrm{\Δ}}_{\mathrm{km}})={\mathrm{\ϵ}}_t---\left(1c\right)$

其中t、i、τ、ε_t和k如上文针对公式(1)所定义，

m为电极数，1≤m≤M，并且

Δkm为在第k个LAT位置处的第m个电极的时延。

在公式(1c)中，参数

对于所有m个同时测量的点是通用的。在一个实施例中，参数

的值可被计算为由公式(1d)给定的各个参数值的平均数或加权平均数，所述公式(1d)为：

γ_k，i[m]=e^{-α·d(i+m，m)}    (1d)

其中α如上文针对公式(1a)所定义，并且

d(i+m，m)为第(i+m)个质心与第k次测量的第m个电极位置之间的距离。

为简洁起见，以下说明假设公式(1)和(1a)适用于由单个电极生成测量值的情况。本领域的普通技术人员将能够修改对由公式(1c)和(1d)所涵盖的情况(即，多个电极同时生成测量值的情况)的说明。

在短程线步骤117中，处理器40计算基础三角形质心和相邻三角形质心之间的空间短程线Δ。处理器对网格中的所有三角形执行该计算。空间短程线使两个质心之间的位移最小化，并且空间短程线Δ可根据公式(2)定义：

$\mathrm{\Δ}=\mathrm{Min}({\stackrel{\→}{d}}_i+{\stackrel{\→}{d}}_j)---\left(2\right)$

其中分别为从基础三角形和相邻三角形质心到所述两个三角形所公共的线的局部位移矢量；所述两个矢量具有在公共线上的公共顶点。

在本公开和权利要求中，元素之间的短程线应被理解为元素之间的最短路径，并且短程线转移是参数沿着此类路径的转移。短程线可以是空间短程线，在该情况下最短路径是元素之间的最短空问路径。公式(3)是使用空间短程线的例子。作为另外一种选择，短程线可以是时间短程线，在该情况下最短路径是元素之间的最短时间路径。

质心C70B和C70A之间的示例性局部位移矢量在图3A中示出，并且具有在线66M上的公共顶点76。(为清楚起见，在线66L上具有公共顶点的质心C70A和C70C之间的位移矢量未在图中示出。)

在LAT和慢度估计步骤119中，对于所有三角形(基础和相邻三角形，即步骤110，111中识别的那些)，处理器对被分配至所述三角形的相应质心的慢度矢量进行计算

慢度值可根据关系(3)定义：

$\stackrel{\→}{S_1}=\frac{\stackrel{\→}{V_j}}{{\left|V_i\right|}^2}---\left(3\right)$

其中

为估计的速度矢量，从而测量第i个基础三角形中电活动流的速度和方向。

如从公式(3)明显可见，慢度是大小为速度的倒数且方向与速度方向相同的矢量。

连接至质心C70A的示例性慢度矢量

在三角形70A中示出。

在步骤119中，处理器40使用如下公式(4)计算所有三角形的慢度矢量。

$\stackrel{\→}{S_1}-\stackrel{\→}{S_j}=\stackrel{\→}{{\mathrm{\ϵ}}_S}---\left(4\right)$

其中i为一个三角形的标识符，

j为相邻三角形i的三角形的标识符，并且

为常数矢量，该常数矢量是相邻位置的慢度值的相似性的量度；

通常，

为慢度预期值的大约0.01％-1％。

的值可由使用者28提供给处理器，并且可适用于所有网格62。在一些实施例中，ε_S的不同值可通常根据与心脏三角形i，j的元素有关的位置提供。

连接至质心C70B的示出了示例性的三角形慢度矢量。

同样在步骤119中，对于每个基础三角形，处理器40使用上述公式(1)计算基础三角形的质心的已分配局部活化作用时间(LAT)。

处理器根据如下公式(5)对分配至三角形质心的活化作用时间进行计算。

$t_j-t_i=\stackrel{\→}{S_i}\·{\stackrel{\→}{d}}_i+\stackrel{\→}{S_j}\·{\stackrel{\→}{d}}_j---\left(5\right)$

其中i为前标记三角形的标识符，

j为邻接前标记三角形i的三角形(即，后标记三角形)的标识符，

t_i、t_j为i和j三角形的质心的相应局部活化作用时间(LAT)，

为i和j三角形的慢度矢量，并且

由公式(2)生成。对公式(5)右侧的检查表明，第一项等于电活动在前标记三角形中的分量从公共三角形线到前标记三角形质心的行进时间。相似地，第二项等于电活动在后标记三角形中在公共三角形线与后标记三角形质心之间的行进时间。

对于前标记三角形质心和后标记三角形质心二者，处理器40应用公式(5)将t_j和S_j的值联系起来。公式(5)的应用针对所有基础三角形和共轭三角形实施，因此生成了由步骤110和111确定的所有三角形(在本文中被称为所关注区)的活化作用时间和慢度。

应当理解，应用公式(5)假设在所关注区的所有三角形质心之间的电活动发生短程线转移。

在延续步骤122中，计算

的过程通过执行对连续围绕所关注区的三角形的计算而扩展。每个前标记三角形通常具有与其相关联的一个或多个另外的前标记三角形，所述另外的前标记三角形形成前标记三角形的相邻序列。可通过分析与慢度矢量相交的边而对所述另外的前标记三角形进行识别。

相似地，每个后标记三角形可具有与其相关联的一个或多个另外的后标记三角形，所述另外的后标记三角形也形成后标记三角形的相邻序列。

例如，在图3A中，慢度矢量与边66P相交，使得存在与三角形70B具有公共边66P的另外的前标记三角形(未示出)。相似地，附接到质心C70C的慢度矢量(未示出)与边66R或66S相交，因此存在另外的后标记三角形，其具有这些边中的一条作为与三角形70C的公共边。

所述两个序列可视为从一个基础三角形开始，并且由所述两个序列形成的路径中的三角形的所有质心均位于空间短程线上。

在步骤122中，处理器实施步骤111的操作，以识别另外的相邻三角形。随着每个另外的相邻三角形被识别出，处理器实施步骤117的操作。处理器继续生成相邻三角形以及它们的活化作用时间，直到遇到网格62的边缘或另一条短程线上的三角形。

正在进行并包括步骤122的操作生成一组空间短程线。每条短程线构成一组相邻的三角形，所述相邻的三角形连接以形成被标测的电活动流所采取的路径。

通常，一旦步骤122已完成，则可能仍然有如下网格62的三角形，所述三角形不在识别出的短程线上，并因此不具有分配的活化作用时间。

在步骤124中，处理器计算执行步骤122后发现的相邻三角形的活化作用时间和慢度。通过求解这些三角形的公式(4)和(5)的系统，以及使用存在于步骤119中的相邻所关注区的三角形的已确定活化作用时间和慢度来计算活化作用时间和慢度。

在最终步骤128中，处理器在显示器48上显示网格62的三角形质心的活化作用时间和/或慢度矢量。通常，将平滑后的等时线结合到显示器中。

通过查看流程图100的说明，应当理解上文呈现的公式(公式(1)、(4)和(5))不是独立的，并且不能依次求解。

图5示出了根据本发明实施例的等时线的示意图。图150示出了由局部活化作用时间生成的等时线，该等时线如上文针对步骤108所述生成。仅使用由最小内推和外推法测量的LAT来绘制所述线。图152示出了根据流程图100，如应用于与图150所用相同的LAT而生成的等时线。为清楚起见，等时线通过指示与所述线相关联的不同时间的不同灰度分开。从图152显而易见的是，图152中等时线的平滑度大大超过图150中的那些。

虽然上述说明通常涉及使系统20分析和调整局部活化作用时间，但应当理解，上述的分析和调整可适用于与器官(例如心脏)相关联的其他生理参数。例如，除随时间运行之外，系统20可被配置成随横跨器官的电压运行。又如，在消融器官期间存在热流动，并且通过器官的热流动可自身表现为器官的温度变化。系统20可被配置成分析和调整器官的测量温度。本领域的普通技术人员将能够识别系统20适用的其他生理参数，例如器官发挥功能时从其发出的红外辐射，并假设所有此类参数包括在本发明的范围内。

此外，虽然上述说明已假设电活动的短程线转移包括通过空间短程线的转移，但本发明的实施例还包括通过时间短程线的转移。本领域的普通技术人员将能够调整说明(例如不使用公式(3)，而是最小化公式(5))，以适应时间短程线转移。

因此应意识到，上述实施例均以举例方式举出，并且本发明不受上文特别显示和描述的内容限制。相反，本发明的范围包括上述各种特征的组合和子组合以及它们的变化形式和修改形式，本领域的技术人员在阅读上述说明时将会想到所述变化形式和修改形式，并且所述变化形式和修改形式并未在现有技术中公开。

Claims (26)

1. 一种用于标测的方法，包括：

接收生理参数的初始的一组测量值，所述测量值在身体器官中的相应位置处测量；

接收包括一系列空间标测图元素的、所述器官的三维(3D)标测图；

在测量所述测量值处的所述相应位置与所述标测图元素的子组之间形成对应关系；

响应于所述对应关系将所述生理参数的相应元素值与除所述子组之外的标测图元素进行关联；

调整所述相应元素值，使得相邻各组的所述标测图元素形成短程线；以及

显示所述器官的标测图，所述标测图示出了所调整的元素值。

2. 根据权利要求1所述的方法，其中所述身体器官包括心脏，并且其中所述生理参数包括所述心脏的局部活化作用时间。

3. 根据权利要求1所述的方法，其中所述空间标测图元素包括具有顶点的平面多边形，所述顶点对应于所述身体器官的壁上的位置。

4. 根据权利要求3所述的方法，其中所述平面多边形包括三角形，并且其中所述子组包括所述相应位置的、相应最近的三角形。

5. 根据权利要求4所述的方法，其中除所述子组之外的所述标测图元素包括一个或多个相邻三角形，所述一个或多个相邻三角形与所述相应最近的三角形中的每一者具有公共边。

6. 根据权利要求5所述的方法，并且包括响应于与所述相应最近的三角形中的每一者相关联的慢度矢量而识别所述一个或多个相邻三角形。

7. 根据权利要求5所述的方法，其中调整所述相应元素值包括调整给定的最近的三角形的质心的元素值以及邻近所述给定的最近的三角形的所述一个或多个三角形的质心的元素值，以形成所述短程线。

8. 根据权利要求7所述的方法，其中调整所述相应元素值包括最小化与所述质心相关联的局部位移矢量。

9. 根据权利要求4所述的方法，其中除所述子组之外的所述标测图元素包括一个或多个三角形，所述一个或多个三角形不与所述相应最近的三角形中的每一者具有公共边。

10. 根据权利要求1所述的方法，其中所述短程线包括空间短程线。

11. 根据权利要求1所述的方法，其中显示所述标测图包括将与所述器官相关联的等时线结合到所述标测图中。

12. 根据权利要求1所述的方法，其中所述短程线包括时间短程线。

13. 根据权利要求1所述的方法，其中接收所述初始的一组测量值包括从所述相应位置处的多个相应电极同时接收所述生理参数的多个测量值。

14. 用于标测的设备，包括：

探针，所述探针被配置成：

生成生理参数的初始的一组测量值，所述测量值在身体器官中的相应位置处测量，以及

生成包括一系列空间标测图元素的、所述器官的三维(3D)标测图；以及

处理器，所述处理器被配置成：

在测量所述测量值处的所述相应位置与所述标测图元素的子组之间形成对应关系，

响应于所述对应关系将所述生理参数的相应元素值与除所述子组之外的标测图元素进行关联，

调整所述相应元素值，使得相邻各组的所述标测图元素形成短程线，以及

显示所述器官的标测图，所述标测图示出了所调整的元素值。

15. 根据权利要求14所述的设备，其中所述身体器官包括心脏，并且其中所述生理参数包括所述心脏的局部活化作用时间。

16. 根据权利要求14所述的设备，其中所述空间标测图元素包括具有顶点的平面多边形，所述顶点对应于所述身体器官的壁上的位置。

17. 根据权利要求16所述的设备，其中所述平面多边形包括三角形，并且其中所述子组包括所述相应位置的、相应最近的三角形。

18. 根据权利要求17所述的设备，其中除所述子组之外的所述标测图元素包括一个或多个相邻三角形，所述一个或多个相邻三角形与所述相应最近的三角形中的每一者具有公共边。

19. 根据权利要求18所述的设备，其中所述处理器被配置成响应于与所述相应最近的三角形中的每一者相关联的慢度矢量而识别所述一个或多个相邻三角形。

20. 根据权利要求18所述的设备，其中调整所述相应元素值包括调整给定的最近的三角形的质心的元素值以及邻近所述给定的最近的三角形的所述一个或多个三角形的质心的元素值，以形成所述短程线。

21. 根据权利要求20所述的设备，其中调整所述相应元素值包括最小化与所述质心相关联的局部位移矢量。

22. 根据权利要求17所述的设备，其中除所述子组之外的所述标测图元素包括一个或多个三角形，所述一个或多个三角形不与所述相应最近的三角形中的每一者具有公共边。

23. 根据权利要求14所述的设备，其中所述短程线包括空间短程线。

24. 根据权利要求14所述的设备，其中显示所述标测图包括将与所述器官相关联的等时线结合到所述标测图中。

25. 根据权利要求14所述的设备，其中所述短程线包括时间短程线。

26. 根据权利要求14所述的设备，其中所述探针包括多个电极，所述多个电极从所述相应位置同时接收所述生理参数的所述初始的一组测量值。

Images