CN104902808B - 在解剖映射系统中估计恢复曲线 - Google Patents

在解剖映射系统中估计恢复曲线 Download PDF

Info

Publication number
CN104902808B
CN104902808B CN201380069283.9A CN201380069283A CN104902808B CN 104902808 B CN104902808 B CN 104902808B CN 201380069283 A CN201380069283 A CN 201380069283A CN 104902808 B CN104902808 B CN 104902808B
Authority
CN
China
Prior art keywords
action potential
potential duration
estimated
cycle length
recovery
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201380069283.9A
Other languages
English (en)
Other versions
CN104902808A (zh
Inventor
普拉莫德辛格·H·塔库尔
巴伦·马斯卡巴
尚塔·阿科特-克里什纳穆尔第
希巴吉·肖梅
苏尼帕·萨哈
艾伦·C·舒罗斯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Boston Scientific Scimed Inc
Original Assignee
Boston Scientific Scimed Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Boston Scientific Scimed Inc filed Critical Boston Scientific Scimed Inc
Publication of CN104902808A publication Critical patent/CN104902808A/zh
Application granted granted Critical
Publication of CN104902808B publication Critical patent/CN104902808B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/25—Bioelectric electrodes therefor
    • A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
    • A61B5/28—Bioelectric electrodes therefor specially adapted for particular uses for electrocardiography [ECG]
    • A61B5/283—Invasive
    • A61B5/287—Holders for multiple electrodes, e.g. electrode catheters for electrophysiological study [EPS]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6846—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be brought in contact with an internal body part, i.e. invasive
    • A61B5/6847—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be brought in contact with an internal body part, i.e. invasive mounted on an invasive device
    • A61B5/6852—Catheters
    • A61B5/6858—Catheters with a distal basket, e.g. expandable basket
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B18/04—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by heating
    • A61B18/12—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by heating by passing a current through the tissue to be heated, e.g. high-frequency current
    • A61B18/14—Probes or electrodes therefor
    • A61B18/1492—Probes or electrodes therefor having a flexible, catheter-like structure, e.g. for heart ablation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B2018/00315—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body for treatment of particular body parts
    • A61B2018/00345—Vascular system
    • A61B2018/00351—Heart
    • A61B2018/00357—Endocardium
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B2018/00571—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body for achieving a particular surgical effect
    • A61B2018/00577—Ablation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods

Abstract

一种用于映射解剖结构的方法包括用设置在所述解剖结构中或者附近的多个电极感测生理活动的激活信号，每个激活信号具有相关联的周期长度，估计每个周期长度的动作电位时程和舒张间隔，基于所估计的动作电位时程和由在前的周期长度所估计的舒张间隔产生恢复曲线，迭代地优化每个所估计的动作电位时程和相应的舒张间隔以最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的函数关系，和基于所优化的动作电位时程和舒张间隔产生动作电位时程恢复曲线。

Description

在解剖映射系统中估计恢复曲线
对相关申请的交叉参考
本申请根据35U.S.C.第119(e)章要求于2013年1月2日提交的题为“在解剖映射系统中估计恢复曲线”的美国临时申请61/748181的权益，在此将其全部通过引用并入。
技术领域
本公开涉及心脏映射系统。更具体地，本公开涉及被配置为基于所记录的激活信号来估计恢复曲线的心脏映射系统。
背景
电恢复是动作电位时程的变化与第一心脏收缩和额外的收缩之间发生的变化的舒张间隔之间的关系。恢复(restitution)反映心脏组织相对于额外收缩开始时的复原性。可以通过对一系列舒张间隔测量动作电位时程来构建电恢复曲线。该曲线最初在短舒张间隔导致动作电位时程更大缩短的位置处非常陡峭。最初陡峭部分之后，随着动作电位时程(APD)在更长的舒张间隔处达到最大而达到平稳状态。整个范围的舒张间隔的电恢复曲线的斜率或者曲线的最陡峭部分的斜率可以用于APD变化对舒张间隔的变化的响应的度量。动作电位时程和不应性的分散的增加与心律失常的风险的增加相关联。由缩短的舒张间隔产生的缩短的动作电位时程改变心肌的不应性，其被认为是以极大地增强了对纤维性颤动的敏感性的方式为折返去极化(reentrant depolarization)建立了路径。
发明内容
本文公开了用于从心脏导管感测的内在心脏激活信号映射解剖结构的方法以及使用这种方法的解剖映射系统的各种实施方案。
在实施例1中，一种用于映射解剖结构的方法包括用设置在所述解剖结构中或者附近的多个电极感测生理活动的激活信号，每个激活信号具有相关联的周期长度，估计每个周期长度的动作电位时程和舒张间隔，基于所估计的动作电位时程和由在前的周期长度所估计的舒张间隔产生初始恢复曲线，迭代地优化每个所估计的动作电位时程和相应的舒张间隔以最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的函数关系，和基于所优化的动作电位时程和舒张间隔产生最终的动作电位时程恢复曲线。
在实施例2中，根据实施例1所述的方法，其中，估计动作电位时程和舒张间隔的步骤包括随机地划分(partition)每个周期长度以估计相应的动作电位时程和舒张间隔。
在实施例3中，根据实施例1或者2所述的方法，其中，迭代地优化每个所估计的动作电位时程和相应的舒张间隔的步骤进一步包括用第一分区分辨率迭代地优化，和在用所述第一分区分辨率优化之后用第二分区分辨率迭代地优化，其中，所述第二分区分辨率高于所述第一分区分辨率。
在实施例4中，根据实施例1-3中的任一个所述的方法，其中，产生所述恢复曲线的步骤包括产生使每个估计的动作电位时程与由在前的周期长度所估计的舒张间隔配对的二维数据集，以及绘制所述二维数据集以产生初始恢复曲线。
在实施例5中，根据实施例1-4中的任一个所述的方法，其中，最大化所述函数关系的步骤包括最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的交互信息，所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的最大信息系数得分，和最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的距离相关度量(a distance correlation measure)中的至少一项。
在实施例6中，根据实施例1-5中的任一个所述的方法，所述方法还包括最大化所产生的恢复曲线和预定的恢复模板之间的相似性。
在实施例7中，根据实施例1-6中的任一个所述的方法，其中，最大化所述相似性的步骤包括最大化所产生的恢复曲线和所述预定的恢复模板之间的相关系数，和最小化所产生的恢复曲线和所述预定的恢复模板之间的均方误差中的至少一项。
在实施例8中，根据实施例1-7中的任一个所述的方法，其中，所述预定的恢复模板对应于指数函数、S形函数和先前产生的经优化的恢复曲线中的任何一项。
在实施例9中，根据实施例1-8中的任一个所述的方法，其中，为所述多个电极中的每一个产生动作电位时程恢复曲线。
在实施例10中，根据实施例1-9中的任一个所述的方法，其中，响应于由所述多个电极中的至少一个制定的预定的起搏方案来感测生理活动的激活信号。
在实施例11中，根据实施例1-10中的任一个所述的方法，所述方法还包括：基于所产生的动作电位时程恢复曲线产生所述解剖结构的解剖映射。
在实施例12中，一种用于映射解剖结构的方法包括用设置在所述解剖结构中或者附近的多个电极感测心脏激活信号，每个激活信号与周期长度相关联，将每个周期长度划分成动作电位时程和舒张间隔，产生恢复数据集从而使得每个数据点包括动作电位时程和其前一个的舒张间隔，产生所述数据点的二维绘图，所述绘图具有预定义的离散网格位置的网格，确定所述恢复数据集内的函数关联，调整每个周期长度的划分使得每个数据点从一个网格位置移动到另一个网格位置，重复所述调整和确定步骤以最大化所述恢复数据集内的所述函数关联，以及基于对应于最大化的函数关联的恢复数据集产生动作电位时程恢复曲线。
在实施例13中，根据实施例12所述的方法，其中，所述函数关联包括最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的交互信息，最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的最大信息系数得分和最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的距离相关度量中的至少一项。
在实施例14中，根据实施例12和13中任一个所述的方法，其中，所述函数关联基于最大化所产生的恢复曲线和预定的恢复模板之间的相似性。
在实施例15中，根据实施例12-14中的任一个所述的方法，其中，所述相似性包括最大化所产生的恢复曲线和所述预定的恢复模板之间的相关系数，以及最小化所产生的恢复曲线和所述预定的恢复模板之间的均方误差中的至少一项。
在实施例16中，根据实施例12-15中任一个所述的方法，其中，所述预定的恢复模板对应于指数函数，S形函数，和先前产生的经优化的恢复曲线中的任何一项。
在实施例17中，根据实施例12-16中的任一个所述的方法，其中，在最大化所述函数关联之后，所述方法包括增加预定义网格的分辨率，和重复所述调整和确定步骤直到所述函数关联在整个恢复数据集进一步增加。
在实施例18中，根据实施例12-17中的任一个所述的方法，其中，响应于由所述多个电极中的至少一个制定的预定的起搏方案来感测所述心脏激活信号。
在实施例19中，根据实施例12-18中的任一个所述的方法，其中，所述预定的起搏方案包括从统计分布取样的带有脉冲间间隔的随机起搏脉冲的至少一个。
在实施例20中，根据实施例12-19中的任一个所述的方法，所述方法还包括基于所产生的动作电位时程恢复曲线产生心脏映射。
在实施例21中，一种用于映射解剖结构的方法包括用设置在所述解剖结构中或者附近的多个电极感测生理活动的激活信号，每个激活信号具有相关联的周期长度，产生动作电位时程和在前的周期长度的函数，基于所产生的函数将每个周期长度分割成动作电位时程和舒张间隔，以及基于所划分的动作电位时程和在前的舒张间隔产生动作电位时程恢复曲线。
在实施例22中，根据实施例21所述的方法，其中，产生函数的步骤包括基于所述在前的周期长度产生每个激活信号的周期长度的直方图，基于所产生的直方图确定周期长度的不应期，基于所确定的不应期确定动作电位时程，和对于多个在前的周期长度重复产生直方图，确定所述不应期，并估计所述动作电位时程的步骤。
在实施例23中，根据实施例20-21任一个所述的方法，其中，产生直方图的步骤包括确定所记录的周期长度的范围，将所确定的范围分割(segment)成等距竖条(equidistance bin)，识别三个连续的激活信号的所有三联体，和为共享第一周期长度竖条(bin)的三联体产生第二和第三激活信号的周期长度的直方图。
在实施例24中，根据实施例20-23中的任一个所述的方法，其中，为所有周期长度竖条确定所述不应期。
在实施例25中，根据实施例20-24中的任一个所述的方法，其中，确定所述不应期的步骤包括识别所产生的直方图的介于零和具有至少等于预定的阈值的频数(count)的第一竖条之间的区域，以及根据所识别的区域确定所述不应期。
在实施例26中，根据实施例20-25中的任一个所述的方法，其中，所估计的动作电位时程是所述不应期的持续时间。
在实施例27中，根据实施例20-26中的任一个所述的方法，其中，确定所述不应期的步骤包括拟合平滑函数至所产生的直方图，和将所产生的直方图的介于零和所述平滑函数下降到低于预定阈值的位置之间的区域识别成所述不应期。
在实施例28中，根据实施例20-27中的任一个所述的方法，其中，所述平滑函数是S形平滑函数。
在实施例29中，根据实施例20-28中的任一个所述的方法，其中，响应于通过所述多个电极中的至少一个制定的预定的起搏方案来感测生理活动的激活信号。
在实施例30中，根据实施例20-29中的任一个所述的方法，所述方法还包括基于所产生的动作电位时程恢复曲线产生心脏映射。
虽然已公开多个实施方案，通过示出和描述了本发明的说明性实施方案的下面的详细描述，仍然有本发明的其它实施方案对于本领域的技术人员变得明显。因此，附图和详细描述应被认为是说明性的而不是限制性的。
附图的简要说明
图1是用于进入(access)身体中的靶标组织区以便诊断和治疗的导管系统的实施方案的示意图。
图2是与图1的系统联合使用的具有筐功能元件承载结构的映射导管的实施方案的示意图。
图3是包括多个映射电极的筐功能元件的实施方案的示意侧视图。
图4是显示包括动作电位时程和舒张间隔的周期长度的动作电位曲线的实施例的示图。
图5是动作电位时程恢复曲线的实施例的示图。
图6A是初始动作电位恢复曲线的示图。
图6B-6C是分割之后的图6A的初始动作电位恢复曲线的示图。
图7A是起搏和窦性心律期间的周期长度直方图的实施例的示图。
图7B是房颤期间的周期长度直方图的实施例的示图。
虽然本发明可修改为各种修改和替代形式，但具体实施方案已在附图中通过实施例示出并在下面详细描述。然而，意图不是将本发明限制于所描述的特定实施例。与此相反，本发明意在覆盖落入由所附的权利要求所限定的本发明的范围内的所有修改，等同物和替代。
详细说明
图1是用于进入身体中的靶标组织区以便诊断和治疗的系统10的示意图。图1大致显示了部署在心脏的左心室中的系统10。可选地，系统10可以部署在心脏的其他区域中，如左心房、右心房或右心室。尽管所示的实施方案示出了系统10用于消融心肌组织，然而可以可选地将系统10(以及本文所描述的方法)配置用于其他组织消融应用，如用于消融前列腺、脑、胆囊、子宫和身体的其他区域中，包括在不一定基于导管的系统中的组织的手术(procedure)。
系统10包括映射探针14和消融探针16。在图1中，每个单独地经由合适的经皮通路通过静脉或动脉(例如股静脉或动脉)引入选定的心脏区域12。可选地，映射探针14和消融探针16可组装在整合结构中用于同时引入和部署在心脏区域12中。
映射探针14具有柔性导管主体18。导管主体18的远端带有三维的多电极结构20。在示出的实施方案中，结构20采取限定开放的内部空间22的筐形式(参见图2)，但是可以使用其他的多电极结构，其中电极结构的几何学和电极位置是已知的。多电极结构20承载多个映射电极24，每个具有电极位置和通道。将每个电极24配置为感测其上将要实施消融过程的解剖区域中的内在生理活动。在一些实施方案中，将所述电极配置成检测解剖结构内的内在生理活动的激活信号，例如心脏活动的激活时间。
电极24电连接到处理系统32。信号线(未示出)电连接到筐结构20上的每个电极24。所述线延伸穿过探针14的主体18并且将每个电极24电连接到处理系统32的输入端，这将在后面更详细地描述。电极24感测解剖区域，例如，心肌组织中的内在电活动。通过处理系统32处理所感测的活动(即，激活信号)以通过产生解剖映射(例如，动作电位时程(APD)恢复曲线或者传导速度(CV)恢复曲线)协助医生识别适于消融的心脏内的位点或者多个位点。
消融探针16包括带有一个或多个消融电极36的柔性导管主体34。该一个或多个消融电极36电连接到射频(RF)发生器37，其被配置为将消融能量传送到一个或多个消融电极36。消融探针16可相对于待处理的解剖特征以及结构20移动。消融探针16可定位在结构20的电极24之间或者可定位成与结构20的电极24相邻，因为一个或者多个消融电极36相对于待处理的组织定位。
处理系统32向显示器40输出所产生的APD图给医生。在示出的实施方案中，处理系统32包括输出显示装置40(例如，CRT、LED显示器或打印机)。装置40以对医生最有用的格式呈现APD图。此外，处理系统32可以产生用于在显示装置40上显示的识别位置的输出，其辅助医生指导(多个)消融电极36与经识别用于消融的位点处的组织接触。
图2示出适合用在图1所示的系统10中的在远端具有电极24的映射导管14的实施方案。映射导管14具有柔性导管主体18，其远端承载被配置为承载映射电极或传感器24的三维结构20。映射电极24感测心肌组织中的内在电活动(例如，激活信号)，然后通过处理系统32处理所感测的活动以通过产生的和显示的APD恢复曲线协助医生识别具有心律紊乱或其它心肌病理的位点或者多个位点。该过程通常被称为映射。然后该信息可以用于确定用于施加合适的治疗(例如消融)到所识别的位点的合适位置，以及用于导航(navigate)一个或者多个消融电极36到所识别的位点。
所示出的三维结构20包括基座部件41和端盖42，在这两者之间柔性花键44通常以圆周间隔的关系延伸。如上所述，三维结构20采用限定开放的内部空间22的筐的形式。在一些实施方案中，花键44由弹性惰性材料制成，例如镍钛诺金属或硅橡胶，并且以弹性的、预加张力(pretensed)的状态在基座部件41和端盖42之间连接，从而弯向和符合它们接触的组织表面。在示出的实施方案中，八个花键44形成三维结构20。附加的或更少的花键44可以在其他实施方案中使用。如图所示，每个花键44承载八个映射电极24。在三维结构20的其他实施方案中，额外的或者更少的映射电极24可以设置在每个花键44上。在示出的实施方案中，三维结构20相对较小(例如，直径为40mm或者更小)。在可选的实施方案中，三维结构20较大(例如，直径为40mm或者更大)。
可滑动的护套50可沿着导管主体30的主轴移动。向前移动护套50(即，朝着远端)导致护套50在三维结构20上移动，从而将结构20折叠成适于引入内部空间(诸如，例如心脏)内的紧凑的、低剖面(low profile)的状态。与此相反，向后移动护套50(即，朝着近端)暴露三维结构20，使结构20弹性地扩展，并呈现图2中示出的预加张力的姿态。三维结构20的实施方案的进一步的细节公开在美国专利第5647870号中，名称为“多电极支撑结构”，其在此通过引用以其整体并入本文。
信号线(未示出)电连接到每个映射电极26。该线延伸穿过映射导管20的主体30进入手柄54，在手柄中它们连接到外部连接器56，所述外部连接器可以是多插头连接器。连接器56将映射电极24电连接到处理系统32。用于处理由映射导管产生的处理信号的映射系统和方法的进一步的细节在名称为“用于引导多电极结构内的可移动电极元件的系统和方法”的美国专利第6070094号，名称为“心脏映射和消融系统”的美国专利第6233491号和名称为“用于改善体腔的已注册的映射的系统和工艺”的美国专利第6735465号中讨论，其公开内容通过引用并入本文。
应当注意，其他多电极结构可以部署在映射导管14的远端。应当进一步注意到，多映射电极24可以设置在多于一个的结构上，而不是例如，图2中示出的单个映射导管14上。例如，如果利用多个映射结构在左心房内映射，那么可以使用包括带有多个映射电极的冠状窦导管和位于左心房中的带有多个映射电极的筐导管的布置。作为另一个实施例，如果利用多个映射结构在右心房内映射，那么可以使用这样一种布置，该布置包括用于在冠状窦中定位的带有多个映射电极的十极电极(decapolar)导管，和用于围绕三尖瓣环定位的带有多个映射电极的环状导管。
虽然映射电极24已被描述为由专用映射探针承载，诸如映射导管14，但是映射电极可以在非映射专用探针或者多功能探针上承载。例如，消融导管，诸如消融导管16，可以被配置为包括设置在导管主体的远端上并且连接到信号处理系统32和引导系统38的一个或者多个映射电极24。作为另一个实施例，处于消融导管的远端的消融电极可以连接到信号处理系统32从而也作为映射电极操作。
为了示出系统10的操作，图3是包括多个映射电极24的筐结构20的实施方案的示意侧视图。在示出的实施方案中，筐结构包括64个映射电极24。映射电极24以八个电极的组(标示为1，2，3，4，5，6，7，和8)设置在八个花键(标示为A，B，C，D，E，F，G和H)中的每一个上。虽然显示64个映射电极24的布置被设置在筐结构20上，然而映射电极24可以可选地以不同的数量布置，布置在不同的结构上，和/或布置在不同的位置中。另外，多个筐结构可以部署在相同或不同的解剖结构中以同时从不同的解剖结构获得信号。
在将筐结构20定位于邻近待处理的解剖结构(例如，心脏的左心房，左心室，右心房，或右心室)之后，处理系统32被配置为记录与解剖结构的生理活动相关的来自每个电极24通道的激活信号，即，电极24测量解剖结构的生理学固有的电激活信号。响应于内在的生理活动或基于由多个电极24中的至少一个制定的预定的起搏方案可以感测生理活动的激活信号。处理系统32可以被配置成基于先前的舒张间隔(pDI)和动作电位时程之间的，以及类似地，先前的舒张间隔和解剖结构的传导速度之间的交互信息相关性来确定恢复曲线。可选地，处理系统32可以被配置为利用基于动作电位时程与先前的舒张间隔(pDI)之间的相关性的直方图来确定恢复曲线。当试图识别或表征心肌病理学(例如心律失常)时，恢复曲线对于医生是有用的。
图4示出了包括周期长度(CL)的动作电位，两次连续心跳的激活之间的持续时间，其由动作电位时程(APD)和舒张间隔(D1)组成，动作电位时程(APD)即同一次心跳的复原和激活之间的持续时间，舒张间隔(D1)即下一次心跳的激活和当前心跳的复原之间的持续时间。决定APD的一个因素是在前的CL。当CL恒定时，APD恒定。然而，当CL被缩短或延长时，在前的APD被分别缩短或延长。APD的变化和变化的先前DI(pDI)之间的关系被称为恢复。恢复反映了解剖结构(例如心脏组织)的复原性，并且对医生来说有诊断意义。
如图5所示，恢复曲线是APD相对于pDI(或CL)变化的绘图并且在理解心脏的病理，诸如异常信号传导和/或纤维性颤动方面具有意义。恢复曲线也可以用来引导动作电位波形插入到激活信号中以可视化激活信号传播。通常对一系列的舒张间隔测量APD恢复曲线并且将APD恢复曲线绘制成类似于或对应于指数曲线。以往，为测量每个对应的DI的APD，通过心脏起搏方案将DI设定为恒定从而人为地产生各种DI。
在一些实施方案中，处理系统32被配置为基于所获取的动作电位的APD和pDI之间的函数关系或关联确定APD恢复曲线。在记录激活信号之后，处理系统32确定每个激活信号的周期长度并估计每个周期长度的APD和DI。对每个周期长度的APD和DI的初始估计可以是随机的或基于先验信息，例如预定的恢复模板或先前确定的恢复曲线。处理系统32产生使每个估计的APD与估计的pDI(即，来自其前一个激活信号的估计的DI)配对的二维恢复数据集。在二维坐标平面上绘制恢复数据集以产生如图5所示的初始恢复曲线，其中，x轴是以毫秒计的pDI而y轴是以毫秒计的APD，因此，恢复曲线上的每个点是(pDI，APD)的形式。由于初始的恢复曲线是基于估计的和/或任意选择的APD的和DI的初始恢复数据集，初始恢复曲线可能不类似于图5中示出的指数曲线。然而，如果在为每个CL估计初始APD和DI时应用先验知识，那么初始恢复曲线可能类似于APD恢复曲线通常特有的指数曲线。
处理系统32确定恢复数据集的两个维度(即，APD和pDI)之间的函数关联。函数关联可包括最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的交互信息、最大信息系数和距离相关性中的至少一项。由于所记录的激活信号的每个CL的实际的APD和DI是未知的，处理系统32能够使用迭代函数关联估计APD和DI。例如，处理系统进行迭代变化以优化初始恢复数据集直到函数关联，即交互信息，最大信息系数，以及距离相关性被最大化，届时从优化的恢复数据集得到的APD恢复曲线应该类似于图5示出的指数曲线。
交互信息相关性量化两个变量，即APD和pDI之间的关联或相关性，而没有它们之间的关系的函数形式的任何假设。最大信息系数是两个变量X和Y，即APD和pDI之间的线性或非线性的关联的强度的量度。距离相关性是任意的，不一定相等维度的两个随机变量或两个随机向量，即APD和pDI之间的统计相关性的量度，并且距离相关性通常来源于距离方差，距离标准偏差和距离协方差。
为了根据函数关联优化初始恢复数据集，处理系统32将初始恢复曲线(图6A)分割成如图6B中所示的具有多个网格位置的网格并且确定初始函数关联。处理系统32选择网格位置并且将点从所选择的网格位置移动到六个相邻的网格位置中的每一个。因为恢复曲线的点由(pDI，APD)表示，为了移动恢复曲线上的点，处理系统32调整对应于该点的所估计的APD和pDI中的至少一个从而使得该点位于相邻的网格位置的其中一个中。在点被迁移到相邻网格位置之后，函数关联被确定。如果在先前的和当前的函数关联之间观察到函数关联的正向变化，那么在恢复数据集中更新该变化。以该方式，使每个网格位置内的点迁移到相邻的网格位置直到对于所有的网格迁移来说所有的网格位置处的函数关联被最大化。所得到的恢复数据集是由优化的APD和pDI值组成的优化的恢复数据集。
在一些实施方案中，为了进一步优化恢复数据集和相应的估计的APD和pDI值，如图6C所示，将网格的分辨率从较粗分辨率增加至较高分辨率，并且迁移点从一个网格位置到相邻的网格位置直到函数关联的迭代过程被进一步最大化。例如，在选择的网格位置处，处理系统32通过调整与位于所选择的网格位置内的恢复数据集的点关联的APD和pDI值将该点迁移至八个相邻的网格位置中的任何一个。确定函数关联，然后将其与先前的函数关联相比较。如果观察到函数关联的正向变化，那么在恢复数据集中更新该变化。对于每个相应的相邻网格位置在每个网格位置处重复该过程直到函数关联被最大化。
在一些实施方案中，恢复数据集从多个电极的所选择的电极组产生。如前所述地根据函数关联优化恢复数据集。将所得到的优化的APD恢复曲线设定为预定的恢复模板并且用作前面提及的用于为相邻的电极组产生优化的恢复数据集的先验知识。
该先验信息可用于最大化所产生的恢复数据集和预定恢复模板之间的相似性。相似性可包括最大化所产生的恢复曲线和预定的恢复模板之间的相关系数或最大化所产生的恢复曲线和预定的恢复模板之间的均方误差。恢复模板可以包括同一患者或不同患者的先前产生的恢复模板，在同一研究中从相邻的或者邻近的电极组产生的恢复曲线，上升指数函数，和S形函数中的任何一项。
在一些实施方案中，处理系统32被配置为基于动作电位时程和在前的周期长度的函数确定恢复曲线。例如，该函数可以是动作电位时程(APD)与在前的周期长度(pCL)之间的直方图相关性。起搏或窦性CL的典型直方图预计是如图7A所示的窄的。然而，在房颤或者随机的起搏方案期间的CL直方图可以是宽的，如图7B所示。如果有足够的激活信号记录，由于不应期60，直方图的下尾收敛为零。不应期是动作电位之后的时间间隔，在该时间间隔期间不能引发新的动作电位。该不应期60可以用来估计每个周期长度的APD，从该APD可以估计相应的DI和pDI以产生所记录的激活信号的APD恢复曲线。
处理系统32感测来自多个电极24的激活信号并为每个激活信号确定周期长度(CL)。响应于内在的生理活动或者基于通过多个电极24中的至少一个制定的预定的起搏方案可以感测生理活动的激活信号。基于在前的周期长度产生确定的CL的直方图。处理系统32确定所有的CL的范围并且将该范围分割成等距竖条。例如，如果所确定的范围是110毫秒(ms)到500毫秒，那么该处理系统可以产生5毫秒间隔的竖条，如110-115毫秒，115-120毫秒，120-125毫秒等。但是，此间隔长度可以更短或更长。如前所述，APD取决于其前一个的DI或pDI。该处理系统识别三个连续的激活信号的所有组并为共享相同的第一周期长度竖条的所有三联体产生后两个(the second two)激活信号的周期长度的直方图。例如，对于具有落入110ms-115ms竖条内的第一周期长度的所有三联体，第二和第三激活信号之间的周期长度被直方图化。该处理系统32被配置成用平滑函数(例如，S形平滑函数)从所产生的直方图过滤噪声和异常值。
基于所产生的直方图估计不应期60。所产生的直方图内介于零和具有至少等于预定阈值的频数的第一竖条之间的区域被识别为不应期60。预定的阈值频数可以是零或基于频数的百分比的相对阈值。确定所有周期长度竖条的不应期60。根据相应的不应期60估计每个竖条的动作电位时程并且可以将动作电位时程设定为相应的不应期的持续时间。
对于所有的先前的周期长度pCL竖条，即115-120毫秒，120-125毫秒等，重复产生直方图，从所产生的直方图确定不应期，以及从所确定的不应期估计APD的过程以获得多个(APD，pCL)数据点，即，对于所有的竖条。可以将数据点绘图以产生APD和先前的周期长度的函数。例如，根据对每个pCL竖条所确定的APD可以产生APD对pCL曲线。处理系统32根据该函数将每个所确定的CL划分为APD和DI。可以将所得到的APD相对于其相应的在前舒张间隔进行绘图从而为所划分的CL产生恢复曲线。
可以对所讨论的示例性实施方案进行各种改变和添加而不偏离本发明的范围。例如，虽然以上描述的实施方案是指特定的特征，但是本发明的范围还包括具有特征的不同组合的实施方案和不包括所有的描述的特征的实施方案。因此，本发明的范围旨在涵盖所有这样的替代，修改和变化落入权利要求的范围之内，以及它们的所有的等同物。

Claims (20)

1.一种用于映射解剖结构的装置，所述装置包括：
映射探针，其具有柔性导管主体和位于所述柔性导管主体的远端的多个电极，其中，所述多个电极被配置成用设置在所述解剖结构中或者附近的所述多个电极感测生理活动的激活信号，每个激活信号具有相关联的周期长度；
用于确定每个激活信号的周期长度的构件；
用于估计每个周期长度的动作电位时程和舒张间隔的构件；
用于基于所估计的动作电位时程和由在前的周期长度所估计的舒张间隔产生恢复数据集的构件；
用于迭代地优化每个所估计的动作电位时程和相应的舒张间隔以最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的函数关系的构件；和
用于基于对应于最大化的函数关系的恢复数据集产生动作电位时程恢复曲线的构件。
2.根据权利要求1所述的装置，其中，估计动作电位时程和舒张间隔包括随机地划分每个周期长度以估计相应的动作电位时程和舒张间隔。
3.根据权利要求2所述的装置，其中，迭代地优化每个所估计的动作电位时程和相应的舒张间隔进一步包括：
用第一分区分辨率迭代地优化；和
在用所述第一分区分辨率优化之后用第二分区分辨率迭代地优化，其中，所述第二分区分辨率高于所述第一分区分辨率。
4.根据权利要求2所述的装置，其中，产生所述恢复曲线包括：
产生使每个估计的动作电位时程与由在前的周期长度所估计的舒张间隔配对的二维数据集；
绘制所述二维数据集以产生初始恢复曲线。
5.根据权利要求1所述的装置，其中，最大化所述函数关系包括以下的至少一项：
最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的交互信息；
最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的最大信息系数得分；
最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的距离相关性。
6.根据权利要求1所述的装置，其还包括：用于最大化所产生的恢复曲线和预定的恢复模板之间的相似性的构件。
7.根据权利要求6所述的装置，其中，最大化所述相似性包括以下的至少一项：
最大化所产生的恢复曲线和所述预定的恢复模板之间的相关系数；和
最小化所产生的恢复曲线和所述预定的恢复模板之间的均方误差。
8.根据权利要求6所述的装置，其中，所述预定的恢复模板对应于指数函数、S形函数和先前产生的优化的恢复曲线中的任何一项。
9.根据权利要求1所述的装置，其中，为所述多个电极中的每一个产生动作电位时程恢复曲线。
10.根据权利要求1所述的装置，其中，响应于由所述多个电极中的至少一个制定的预定的起搏方案来感测生理活动的激活信号。
11.根据权利要求1所述的装置，所述装置还包括：
用于基于所产生的动作电位时程恢复曲线产生所述解剖结构的解剖映射的构件。
12.一种用于映射解剖结构的装置，所述装置包括：
映射探针，其具有柔性导管主体和位于所述柔性导管主体的远端的多个电极，其中，所述多个电极被配置成用设置在所述解剖结构中或者附近的所述多个电极感测心脏激活信号，每个激活信号与周期长度相关联；
用于将每个周期长度划分为动作电位时程和舒张间隔的构件；
用于产生恢复数据集从而使得每个数据点包括动作电位时程和其前一个的舒张间隔的构件；
用于产生所述数据点的二维绘图的构件，所述绘图具有预定义的离散网格位置的网格；
用于确定所述恢复数据集内的函数关联的构件；
用于调整每个周期长度的划分使得每个数据点从一个网格位置移动到另一个网格位置的构件；
用于重复所述调整和确定步骤以最大化所述恢复数据集内的所述函数关联的构件；以及
用于基于对应于最大化的函数关联的恢复数据集产生动作电位时程恢复曲线的构件。
13.根据权利要求12所述的装置，其中，所述函数关联包括以下的至少一项：
最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的交互信息；
最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的最大信息系数得分；以及
最大化所估计的动作电位时程和由在前的周期长度所估计的舒张间隔之间的距离相关度量。
14.根据权利要求12所述的装置，其中，所述函数关联是基于最大化所产生的恢复曲线和预定的恢复模板之间的相似性。
15.根据权利要求14所述的装置，其中，所述相似性包括以下的至少一项：
最大化所产生的恢复曲线和所述预定的恢复模板之间的相关系数；以及
最小化所产生的恢复曲线和所述预定的恢复模板之间的均方误差。
16.根据权利要求14所述的装置，其中，所述预定的恢复模板对应于指数函数，S形函数，和先前产生的优化的恢复曲线中的任何一项。
17.根据权利要求12所述的装置，其进一步包括：
用于在最大化所述函数关联之后增加预定义网格的分辨率的构件；和
用于在最大化所述函数关联之后重复所述调整和确定步骤直到所述函数关联在整个恢复数据集进一步增加的构件。
18.根据权利要求12所述的装置，其中，响应于由所述多个电极中的至少一个制定的预定的起搏方案来感测所述心脏激活信号。
19.根据权利要求18所述的装置，其中，所述预定的起搏方案包括从统计分布取样的带有脉冲间隔的随机起搏脉冲的至少一个。
20.根据权利要求12所述的装置，其还包括：
用于基于所产生的动作电位时程恢复曲线产生心脏映射的构件。
CN201380069283.9A 2013-01-02 2013-12-20 在解剖映射系统中估计恢复曲线 Active CN104902808B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201361748181P 2013-01-02 2013-01-02
US61/748,181 2013-01-02
PCT/US2013/077143 WO2014107355A1 (en) 2013-01-02 2013-12-20 Estimating restitution curves in an anatomical mapping system

Publications (2)

Publication Number Publication Date
CN104902808A CN104902808A (zh) 2015-09-09
CN104902808B true CN104902808B (zh) 2018-04-24

Family

ID=49943587

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201380069283.9A Active CN104902808B (zh) 2013-01-02 2013-12-20 在解剖映射系统中估计恢复曲线

Country Status (5)

Country Link
US (1) US9089272B2 (zh)
EP (1) EP2941181B1 (zh)
JP (1) JP6144361B2 (zh)
CN (1) CN104902808B (zh)
WO (1) WO2014107355A1 (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10695126B2 (en) 2008-10-06 2020-06-30 Santa Anna Tech Llc Catheter with a double balloon structure to generate and apply a heated ablative zone to tissue
US9259156B2 (en) * 2014-07-10 2016-02-16 General Electric Company Methods and systems to translate two dimensional mapping into a three dimensional derived model
US10398509B2 (en) 2015-09-18 2019-09-03 General Electric Company System and method for optimal catheter selection for individual patient anatomy
CN108348186B (zh) * 2015-10-07 2021-06-18 圣犹达医疗用品心脏病学部门有限公司 用于标测心脏恢复的方法和系统
EP3344136B1 (en) 2015-10-07 2020-07-01 St. Jude Medical, Cardiology Division, Inc. Methods and systems for mapping cardiac repolarization
US11331140B2 (en) 2016-05-19 2022-05-17 Aqua Heart, Inc. Heated vapor ablation systems and methods for treating cardiac conditions
EP3456246A1 (en) * 2017-09-19 2019-03-20 Tata Consultancy Services Limited A cascaded binary classifier for identifying rhythms in a single-lead electrocardiogram (ecg) signal
KR102257498B1 (ko) * 2019-05-29 2021-05-27 연세대학교 산학협력단 활동 전위 기간 상환 현상을 반영한 심장 모델 생성 방법 및 생성 장치

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4421114A (en) * 1978-10-30 1983-12-20 Berkovits Barouh V Tachycardia treatment
US8165666B1 (en) * 2011-05-02 2012-04-24 Topera, Inc. System and method for reconstructing cardiac activation information

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6094593A (en) 1991-05-17 2000-07-25 Cedars-Sinai Medical Center Method and apparatus for detecting prospenity for ventricular fibrillation using action potential curves
DE69425249T2 (de) 1993-03-16 2001-03-22 Ep Technologies Träger-anordnung für mehrfach-elektroden
US6233491B1 (en) 1993-03-16 2001-05-15 Ep Technologies, Inc. Cardiac mapping and ablation systems
US5876336A (en) 1994-10-11 1999-03-02 Ep Technologies, Inc. Systems and methods for guiding movable electrode elements within multiple-electrode structure
US6735465B2 (en) 2001-10-24 2004-05-11 Scimed Life Systems, Inc. Systems and processes for refining a registered map of a body cavity
US7123954B2 (en) * 2002-09-19 2006-10-17 Sanjiv Mathur Narayan Method for classifying and localizing heart arrhythmias
WO2005006946A2 (en) 2003-07-03 2005-01-27 New York Univeristy System and method for assessment of cardiac electrophysiologic stability and modulation of cardiac oscillations
US8676303B2 (en) 2008-05-13 2014-03-18 The Regents Of The University Of California Methods and systems for treating heart instability
WO2010042826A1 (en) 2008-10-09 2010-04-15 Sanjiv Narayan Machine and process for the automatic localization of sources of biological rhythm disorders
US9282910B2 (en) 2011-05-02 2016-03-15 The Regents Of The University Of California System and method for targeting heart rhythm disorders using shaped ablation

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4421114A (en) * 1978-10-30 1983-12-20 Berkovits Barouh V Tachycardia treatment
US8165666B1 (en) * 2011-05-02 2012-04-24 Topera, Inc. System and method for reconstructing cardiac activation information

Also Published As

Publication number Publication date
CN104902808A (zh) 2015-09-09
WO2014107355A1 (en) 2014-07-10
US20140187989A1 (en) 2014-07-03
US9089272B2 (en) 2015-07-28
JP6144361B2 (ja) 2017-06-07
JP2016508057A (ja) 2016-03-17
EP2941181B1 (en) 2020-01-22
EP2941181A1 (en) 2015-11-11

Similar Documents

Publication Publication Date Title
CN104902808B (zh) 在解剖映射系统中估计恢复曲线
CN108471975B (zh) 用于统计分析局部异常心室活动的电描记图并标测该局部异常心室活动的方法和系统
De Groot et al. Critical appraisal of technologies to assess electrical activity during atrial fibrillation: A position paper from the European heart rhythm association and European society of cardiology working group on eCardiology in collaboration with the heart rhythm society, asia pacific heart rhythm society, Latin American heart rhythm society and computing in cardiology
US9579034B2 (en) Enhanced activation onset time optimization by similarity based pattern matching
CN113749664A (zh) 使用自动编码器减小心内心电图噪声及使用深度学习训练损失函数以细化心内和体表心电图
KR20210027077A (ko) 다수의 전극 카테터에 의해 측정된 국소적 활성화 시간(lat)의 오차 추정
CN105708442B (zh) 心室远场衰减
JP2012509701A (ja) 心臓を画像化する画像化装置
US20220202345A1 (en) Method and apparatus to find abnormal activations in intra-cardiac electrocardiograms based on specificity and sensitivity
CN114587369A (zh) 基于学习数据分析来设置自动感兴趣窗口
CN107249446A (zh) 用于检测器官中电活动的导管和方法
CN114601468A (zh) 用于解剖结构标测数据的信号和校正处理
CN114795227A (zh) 通过导管的在线磁性校准来管理医疗装置装备
CN114052657A (zh) 使用多个成像模态来自动识别有机组织内的疤痕区域
CN114052891A (zh) 大面积环状消融点的解剖结构的自动分割
US20220192576A1 (en) Cardiac mapping via near-field identification from physiological signals
US20230181246A1 (en) Cardiac vein ablation visualization system and catheter
US20230181087A1 (en) Intracardiac unipolar far field cancelation using multiple electrode cathethers
EP4079223A1 (en) Visualization of ventricular tachycardia causing reentrant circuits via pseudo-activation maps
US20230146716A1 (en) Digital twin of atria for atrial fibrillation patients
WO2023105493A1 (en) Cardiac vein ablation visualization system and catheter
CN114681039A (zh) 经由视觉信号进行消融后验证
CN114642435A (zh) 自动起搏后间期测量和显示
CN114617560A (zh) 用于解剖结构的相干标测的信号流的速度流的信号处理

Legal Events

Date Code Title Description
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant