CN110859612A - 自动识别异常lat差异 - Google Patents
自动识别异常lat差异 Download PDFInfo
- Publication number
- CN110859612A CN110859612A CN201910794757.5A CN201910794757A CN110859612A CN 110859612 A CN110859612 A CN 110859612A CN 201910794757 A CN201910794757 A CN 201910794757A CN 110859612 A CN110859612 A CN 110859612A
- Authority
- CN
- China
- Prior art keywords
- chamber
- locations
- processor
- map
- local
- 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.)
- Pending
Links
- 230000002159 abnormal effect Effects 0.000 title abstract description 6
- 230000004044 response Effects 0.000 claims abstract description 20
- 230000002107 myocardial effect Effects 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 26
- 239000003086 colorant Substances 0.000 claims description 13
- 230000004913 activation Effects 0.000 claims description 11
- 238000004040 coloring Methods 0.000 claims description 2
- 230000000875 corresponding effect Effects 0.000 description 12
- 239000000523 sample Substances 0.000 description 12
- 238000001994 activation Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000002547 anomalous effect Effects 0.000 description 3
- 210000003748 coronary sinus Anatomy 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 241000204031 Mycoplasma Species 0.000 description 2
- 238000010009 beating Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 229930186657 Lat Natural products 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
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/339—Displays specially adapted therefor
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/06—Devices, other than using radiation, for detecting or locating foreign bodies ; determining position of probes within or on the body of the patient
- A61B5/065—Determining position of the probe employing exclusively positioning means located on or in the probe, e.g. using position sensors arranged on the probe
- A61B5/066—Superposing sensor position on an image of the patient, e.g. obtained by ultrasound or x-ray imaging
-
- 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/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
-
- 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]
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
- A61B5/743—Displaying an image simultaneously with additional graphical information, e.g. symbols, charts, function plots
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/06—Devices, other than using radiation, for detecting or locating foreign bodies ; determining position of probes within or on the body of the patient
- A61B5/061—Determining position of a probe within the body employing means separate from the probe, e.g. sensing internal probe position employing impedance electrodes on the surface of the body
- A61B5/062—Determining position of a probe within the body employing means separate from the probe, e.g. sensing internal probe position employing impedance electrodes on the surface of the body using magnetic field
Abstract
本发明题为“自动识别异常LAT差异”。本发明提供了一种设备,所述设备由导管组成,所述导管被配置为插入心脏的腔室中,并且具有一个或多个电极,所述一个或多个电极被配置为接触所述腔室中多个位置处的心肌组织。一个或多个电极接收响应于行进穿过所述组织的传导波的电信号。设备包括显示器和处理器,所述处理器被配置为接收来自所述导管的所述电信号,并且响应于所述电信号向所述显示器呈现所述腔室的标测图,所述腔室的标测图包括所述多个位置处的所述传导波的发生的本地时间的指示。处理器还被配置为响应于所述发生的本地时间来计算所述位置之间的所述传导波的速度,并且在所述标测图上标记其中所述速度低于预设阈值的、所述腔室的一个或多个区域。
Description
相关申请的交叉引用
本专利申请要求2018年8月27日提交的美国临时专利申请62/723,218的权益,该临时专利申请以引用方式并入本文。
技术领域
本发明整体涉及电生理电位,具体地讲,涉及具有异常电生理电势的组织的识别。
背景技术
一种用于心脏研究的标准方法,包括产生跳动心脏的腔室的本地激活时间(LAT)的标测图。在窦性节律中操作的“正常”心脏产生针对不同腔室的良好限定的LAT标测图。与这些标测图的偏差允许医师检查标测图以识别可能存在问题的腔室区。然而,从LAT标测图本身来看,此类识别可能不是显而易见的。
发明内容
本发明的一个实施方案提供了由导管组成的设备,该导管被配置为插入心脏的腔室中,并且具有一个或多个电极,所述一个或多个电极被配置为接触腔室中多个位置处的心肌组织并且接收响应于行进穿过组织的传导波的电信号。
该设备还具有显示器和处理器,该处理器被配置为接收来自导管的电信号并且响应于该电信号向显示器呈现腔室的标测图,该腔室的标测图包括所述多个位置处的传导波的发生的本地时间的指示。处理器还响应于发生的本地时间来计算位置之间的传导波的速度,并且在标测图上标记其中速度低于预设阈值的、腔室的一个或多个区域。
在一个公开的实施方案中,处理器被配置为计算在多个位置之间的相应位置处的内插的发生的本地时间,并且在相应位置处将内插的发生的本地时间的指示标记在标测图上的。通常,处理器被配置为不在所述一个或多个区域处将内插的发生的本地时间的指示标记在标测图上。
在另一个公开的实施方案中,处理器被配置为将腔室分成邻接表面,每个表面围绕所述多个位置中的相应一个位置,并且所述一个或多个区域位于两个邻接表面之间的边界处。
处理器可被配置为通过将腔室的表面形成为体素组而将腔室分成邻接表面,并且,对于所述多个位置中的每个给定位置,迭代地并且同时地将紧邻体素添加到先前体素,该先前体素最初包括对应于给定位置的本地原体素,以便形成给定的邻接表面。
处理器还可被配置为通过找到边界上的区与两个邻接表面的相应位置的距离的总和来计算区处的传导波的速度,并且将该总和除以在相应位置处的发生的本地时间之间的时间差。
在另一个公开的实施方案中,电信号包括心电图(ECG)信号,并且发生的本地时间对应于从ECG信号推出的本地激活时间(LAT)。通常,处理器被配置为将来自具有与LAT的值一一对应的相应颜色的色标的颜色结合到标测图中。处理器可被配置为以与色标的颜色不同的颜色来对所述一个或多个区域着色。
在一个另选实施方案中,预设阈值介于0.01mm/ms和1mm/ms之间。
在另一个另选实施方案中,预设阈值为0.1mm。
在又一个另选实施方案中,导管包括位置传感器,该位置传感器被配置为向处理器提供指示所述多个位置的信号。通常,位置传感器包括一个或多个线圈,该一个或多个线圈响应于穿越线圈的磁场来提供信号。
本发明的另一个实施方案提供方法,包括:
将导管插入心脏的腔室中,该导管具有一个或多个电极,所述一个或多个电极被配置为接触腔室中多个位置处的心肌组织并且接收响应于行进穿过组织的传导波的电信号;
提供显示器;
接收来自导管的电信号;
响应于电信号来向显示器呈现腔室的标测图,该腔室的标测图包括所述多个位置处的传导波的发生的本地时间的指示;
响应于发生的本地时间来计算所述位置之间的传导波的速度;以及在标测图上标记其中速度低于预设阈值的、腔室的一个或多个区域。
结合附图,通过以下对本公开的实施方案的详细描述,将更全面地理解本公开,其中:
附图说明
图1为根据本发明的实施方案的使用设备的侵入式医疗规程的示意图;
图2为根据本发明的实施方案的用于设备中的探针的远侧端部的示意图;
图3为在根据本发明的实施方案的设备的操作期间所执行的步骤的流程图;以及
图4A、图4B和图5示出了根据本发明的实施方案的流程图的一些步骤。
具体实施方式
概述
通常,在产生心脏的腔室的本地激活时间(LAT)标测图时,对腔室表面上的多个点进行取样,并且从取样的数据计算这些点中的每一个点的LAT值。在取样点中间的表面上的点的LAT值通常通过从取样点的插值来发现。LAT值通常以不同的颜色结合到腔室的3D标测图中,然后可将包括取样和内插LAT值的着色标测图呈现给医师进行分析。
然而,着色标测图本身并不一定提供对腔室的异常区的指示,尤其是因为在异常区的情况下,上述插值可给出错误的结果并因此隐藏任何异常的LAT值。例如,对着色标测图的检查可不给出任何指示行进穿过腔室的传导波的速度减小的指示,并且此类减小的速度指示异常区。
本发明的实施方案提供了一种使用对取样数据的进一步分析来通过估计传导波的速度来指示异常区的方法。该指示在“常规”LAT着色标测图上提供,使得不需要产生或查看示出传导波速度的不同标测图。
首先,将探针插入所研究的腔室中,并且将探针用于采集位置数据信号,该位置数据信号被分析以给出腔室表面的三维(3D)位置。探针上或另一探针上的一个或多个电极也用于在表面上的选定位置处采集定时数据信号,通常为心电图(ECG)信号。可分析ECG信号以计算选定位置处的LAT值,并且可生成3D腔室位置的标测图,其中LAT值作为颜色结合到标测图中。
除了生成标测图之外,还分析选定位置及其相应LAT值,以计算在选定位置之间的表面区处的传导波的速度。当其中一个区的一部分的速度低于预设阈值时,在该区的该部分的标测图上放置视觉指示。
具体实施方式
图1为使用设备12的侵入式医疗规程的示意图,并且图2为根据本发明的实施方案的用于该设备中的导管20的远侧端部22的示意图。导管20在本文中也称为探针20。该规程由医疗专业人员14执行,并且在下文的描述中,假定该规程包括人类患者18的心脏16的腔室15的三维(3D)表面17的电极电位(EP)研究。
为了进行研究,专业人员14将探针20插入到已预先定位在患者的内腔中的护套21中。护套21被定位成以便于探针的远侧端部22进入患者的心脏。远侧端部22包括位置传感器24,该位置传感器能够跟踪远侧端部的位置和取向。远侧端部22还包括用于从表面17采集电极电位的一个或多个电极30,如下所述。为清楚和简单起见,以下描述假定仅存在一个电极30。
设备12由系统处理器46控制,该系统处理器位于设备的操作控制台48中。控制台48包括由专业人员14用来与处理器通信的控件49。可将用于处理器46的软件通过例如网络以电子形式下载到处理器。另选地或除此之外,软件可通过非暂时性有形介质诸如光学、磁或电子存储介质提供。
从由电极30采集的EP信号推出的电极电位(3D)数据呈现在形成于显示器62上的患者18的心脏的三维(3D)表示60上。EP信号响应于行进穿过心脏的心肌组织的传导波而生成。
在本发明的实施方案中,使用色标64来呈现3D表示,其中以一对一方式对不同的颜色分配从EP数据推出的参数的相应不同值。本文所用的参数是表面17的区的本地激活时间(LAT),LAT对应于传导波的发生的本地时间。在一个实施方案中,LAT从约-150ms的低值变化至约90ms的高值,如相对于基准参考信号所测量的,该基准参考信号通常从定位在心脏16的冠状窦中的探针推出。在这种情况下,对于低LAT值,色标64被预设为红色,对于高LAT值为紫色,并且对于中LAT值的预设颜色为诸如橙色、黄色、绿色、蓝色的色调。
系统处理器46包括通常被配置为现场可编程门阵列(FPGA)的实时降噪电路45,之后是模数(A/D)信号转换集成电路47。处理器可将信号从A/D电路47传递到另一个处理器,和/或可被编程以执行本文所公开的至少一个算法,该算法包括下文所公开的步骤,并且被存储在与处理器46通信的存储器72中。处理器使用电路45和电路47,以及下文详述的模块的特征,以便执行算法。
为了操作设备12,处理器46的算法与模块库50通信,该模块库具有被处理器用来操作设备的多个模块。因此,库50包括采集和分析来自电极30的信号的心电图(ECG)模块56、以及接收和分析来自位置传感器24的信号并利用信号分析来生成远侧端部22的位置和取向的追踪模块58。在一些实施方案中,传感器24包括一个或多个线圈,所述一个或多个线圈响应于穿越线圈的磁场来提供传感器信号。在这些实施方案中,除了接收和分析来自传感器24的信号之外,追踪模块58还控制辐射器32、34和36,所述辐射器辐射穿越传感器24的磁场。这些辐射器被定位成接近心脏16,并且被配置成将交变磁场辐射到接近心脏的区中。由Biosense Webster(33Technology Drive,Irvine,CA 92618USA)制造的系统使用此类磁跟踪系统。
图3为根据本发明的实施方案的在设备12的操作期间由处理器46和专业人员14执行的算法的步骤的流程图,图4A、图4B和图5示出了根据本发明的实施方案的一些步骤。流程图描述了标测腔室15的3D表面17的步骤,并且还描述了从表面采集信号的步骤。虽然在一些实施方案中,这两个功能由独立的探针执行,但为简单起见,在以下描述中,假定由探针20执行这两个功能。针对用于不同功能的不同探针的情况的实施方案的具体实施也被认为是在本发明的范围内。
在初始步骤100中,产生了患者18的心脏的腔室15的表面17的3D标测图。通常通过将远侧端部22插入腔室中,并将远侧端部移动到腔室内来进行生产。当移动远侧端部时,跟踪模块58从位置传感器24采集邻接数据位置信号,并且该模块与处理器46一起分析每个信号,以产生用于远侧端部的相应3D位置。
处理器46分析有效地为点云的一组3D位置,以产生包封点云的网片,并且处理器从网片生成腔室的表面17作为一组体素。用于生成包封网片和针对点云的对应表面的方法是本领域所熟知的。
在信号采集步骤102中,远侧端部22的电极30在表面17的不同点上被触摸,并且在每个点处,被触摸的点的位置P从从传感器24采集的信号推出。
同时,在被触摸的点处生成的信号由ECG模块56从电极30采集,该模块也记录信号。对于任何给定的点“a”,采集的信号为一组电压时间有序对:{(Vt,t)}a,其中Vt为时间t处的电极电位。在一些实施方案中,可将{(Vt,t)}a呈现给专业人员14作为显示器62上的电压对时间的曲线图。
ECG模块56和处理器46分析所记录的信号中的每一个,{(Vt,t)},以导出电极30所触摸的每个点的相应本地激活时间(LAT)。由公式(1)表示该分析:
LATa=f[{(Vt,t)}a] (1)
其中a为表面17上的点的标识符,f表示应用于{(Vt,t)}以在点a,LATa处导出本地激活时间的函数。
如上所述,LATa的值通常是相对于预定参考时间信号测量的时间,诸如由插入心脏16的冠状窦中的参考探针所提供的时间。
因此,对于给定的点a,步骤102找到有序对(位置、本地激活时间):(Pa,LATa)。
由于心脏16的跳动,在步骤100中标测的表面不是静止表面,使得由电极30所触摸并在步骤102中导出的点的位置通常可不与在步骤100中生成的表面重合。投影步骤104解决了这个点。
在投影步骤104中,处理器46将在步骤102中找到的每个位置投影到在步骤100中导出的体素表面17上。图4A为投影的示意图,并且以举例的方式示出了接近3D表面17的三个点a、b、c的位置Pa、Pb和Pc。投影对应于表面的相应体素Va、Vb和Vc。在步骤104中,处理器存储通过投影在步骤102中确定的所有点而发现的体素组。从投影推出的体素在本文中被称为本地原体素。
本地原体素根据标度64使用被投影位置的相应LAT值来着色。表面17上的中间体素的LAT值(即,非本地原体素的体素)通常通过内插来确定,并且中间体素也根据标度64来着色。在一个实施方案中,中间体素的插值发现最近的本地原体素的LAT值的平均值,该平均值根据与本地原体素的距离来加权。
在传播步骤106中,处理器46限定3D表面17中的相应区,即,围绕每个本地原体素的相应体素组。为了定义每个区,处理器为该区中的每个体素限定共同标识符Vid,其中id是该区的数字标识符,并且也是该区的本地原体素的数字标识符。如下所述,处理器还测量区中每个体素与该区中的本地原体素在3D表面17上方的距离D。因此,给定区中的每个体素具有与之相关联的有序对(体素标识符,与本地原体素的距离),本文以(Vid,D)表示。区的每个本地原体素的有序对为(Vid,0)。
在限定每个本地原体素的有序对(V1,0)、(V2,0)、(V3,0),…之后,处理器分配表面17的其他体素以具有临时有序对,从而使得处理器能够区分非本地原体素与本地原体素。以举例的方式,在以下描述中,对所有非本地原体素分配有序对(0,0);即,对非本地原体素的标识符分配值0。然而,假定将非本地原体素与本地原体素区分开的其他方法被认为包括在本发明的范围内。
一旦表面17的所有体素已被分配如上所述的有序对,处理器就继续“传播”本地原体素中的每一个本地原体素的标识符。为了实现传播,对于标识符为Vg的给定本地原体素,处理器确定紧邻体素,并检查有序对的第一项(体素的标识符)是否为0。
如果为0,则处理器对该体素分配标识符Vg,并且还计算与本地原体素的距离D,该距离被分配给该体素的有序对的第二项。
如果该体素的标识符不为0,则处理器对另一紧邻体素进行检查,并且继续直至所有紧邻体素已被检查。
处理器继续检查邻近已检查的紧邻体素的另外的非本地原体素,在每种情况下均检查体素标识符是否为0。对于标识符为0每个体素,处理器对该体素分配标识符Vg。处理器还通过将要比较的两个体素之间的距离与具有已经存在的标识符Vg的体素的距离相加来计算到本地原体素在3D表面17上方的距离D。应当理解,分配给任何特定体素的距离等于所有紧邻体素之间的累积距离。
上述过程迭代,以便体素的邻接区围绕本地原体素生长。生长区的体素中的每一个体素均具有与本地原体素的标识符Vg共同的标识符,并且距离D等于从本地原体素计算出的在3D表面17上方的距离。
处理器46将上述过程同时应用于所有本地原体素,使得在由步骤106实施的过程中存在多个生长区,每个区源自本地原体素。应当理解,当处理器不再发现标识符为0的另外的相邻体素时,并且当相邻体素中的至少一个具有不同于本地原体素的标识符时,每个区终止生长。
条件108检查是否已对表面17的所有体素分配标识符和到本地原体素的距离。如果条件返回负值,则处理器继续传播步骤106。
如果条件返回正值,则控制继续进行如下所述的距离计算步骤110。应当理解,在完成上文针对步骤106所述的方法时,表面17被分成一组邻接区,该组邻接区是接触的,即,区之间没有间隙。每个区包括一组邻接体素,该组邻接体素具有与该区的本地原体素的标识符相对应的共同体素标识符。如上所述,给定区的每个体素也与在3D表面17上方测量的与该给定区的本地原体素的距离相关联。
图4B示意性地示出了步骤106的结论,以举例的方式示出了分别与表面17上的本地原体素Va、Vb和Vc相关联的三个区150、154和158。每个区由至少一个其他区邻接,并且在所示的示例中,每个区由两个区邻接。因此,区150和154之间存在边界160,区154和158之间存在边界164,并且区158和150之间存在边界168。
在图4B中,为清楚起见,区边界160、164和168显示为实线;然而,应当理解,在实施过程中,本地原体素周围的区的边界在呈现在显示器62上的表面17的表示上不可见。此外,应当理解,由于在步骤106中生成的区边界与边界体素的LAT值不相关,因此边界体素的内插颜色(如步骤104中所示)与边界体素的位置之间通常没有关系。
区边界包括一组n体素,其中n为正整数,每个边界体素具有与本地原体素的标识符相同的标识符。此外,每个边界体素具有至少一个最近的相邻体素,其具有不同于本地原体素的标识符。因此,如果区150的本地原体素Va具有标识符A,并且如果区154的本地原体素Vb具有标识符B,则边界160处的区150中的n个边界体素中的每一个均具有标识符A,并且这些边界体素中的每一个均具有标识符为B的最近的边界。在边界160处的区154中存在对应的体素组,这些体素中的每一个均具有标识符B以及标识符为A的最近的边界。
当条件108返回正值时,处理器记录每个边界的所有n个体素的位置。
在距离计算步骤110中,处理器分析与边界体素相关联的参数,以便计算在3D表面17上方测量的相邻本地原体素之间的距离。在边界160中,区150的n个边界体素具有相应的有序对:(A,DVa1)、(A,DVa2)、(A,DVa3)、…、(A,DVan),其中DVan为边界160的第n个边界体素与本地原体素Va的距离。边界160的n个边界体素中的每一个均具有标识符为B的至少一个最近的相邻体素。因此区154中还存在具有相应有序对(B,DVb1)、(B,DVb2)、(B,DVb3)、…、(B,DVbn)的一组n个边界体素,其中DVbn为边界160的第n个边界体素与本地原体素Vb的距离。
对于边界160,处理器46选择邻近边界体素对。对于一对给定的对,其中一个体素位于区150中,另一个体素位于区154中。每一对具有参数(A,DVak)、(B,Dvbk)Vbk),其中k为1和n之间的整数。从本地原体素Va到本地原体素Vb的距离经由具有参数(A,DVak)、(B,DVbk)的边界对体素由下式给出:
Dabk=Dvak+Dvbk (2)
当处理器使用公式(2)经由相应边界体素对已计算出在步骤106中生成的所有边界的相邻本地原体素之间的所有距离时,步骤110结束。对于两个通用的相邻区p、q,这些距离在本文中一般被称为Dpq。
在定时步骤112中,处理器计算邻近本地原体素之间的本地激活时间的差值。为了计算该差值,处理器首先估计心脏16的中值循环长度CL,即连续心室激活之间的中值时间段。处理器可通过任何便利的方式估计CL,诸如分析从上述冠状窦参比探针推出的时间,或者分析来自体表电极的时间。
为了考虑LAT值的循环性质,处理器然后根据以下公式来计算邻近本地原体素之间的LAT差值:
ΔLATpq=min(|LATp-LATq|,[CL-|LATp-LATq|]) (3)
其中LATp、LATq为相邻区p、q的相应本地原体素的LAT,并且
ΔLATpq为LAT的差值。
在步骤112中,处理器计算表面上的邻近区p、q的所有对的ΔLATpq值。
处理器然后检查所计算的速度是否低于预设阈值速度V阈值,即,处理器检查表达式(4)是否有效。V阈值的典型值为约0.1mm/ms,而V阈值的典型范围为约0.01mm/ms至约1mm/ms。
如果对表达式(4)的比较返回正值,即,计算的速度低于预设阈值速度,则在区分色步骤116中,处理器在考虑到与标度64的颜色不同的区分颜色的情况下对边界体素对进行着色。例如,如果褐色不在色标64中,则处理器可将表达式(4)为正值,棕色的一对边界体素进行着色。
图5示出了其中对表达(4)的比较返回正值的情况。图5大致类似于图4B,但与图4B相反,边界160、164、168未以实线示出;然而,边界仍然存在。
在图5中,线180被假定为在边界160中包括边界体素对,其中表达式(4)返回正值。因此,从步骤116开始,线180的体素对被着色为不同于标度64的颜色的区分颜色。
呈现不同着色的区,诸如对于线180的区,有助于专业人员14判定传导波穿过心脏的行进方向。该线有效地示出了对可能从区150行进至区154的波的阻断。因此,传导波可能沿着类似于实线箭头188的路径行进,而不是沿着作为可能的波路径的虚线箭头184行进。
返回图3的流程图,如果对表达式(4)的比较返回负值,即,计算的速度等于或高于预设阈值速度,则在标准着色步骤118中,处理器根据标度64来对边界体素对进行着色,即,通过使用它们的内插LAT值来从标度中选择对应的颜色。
如本文所用,针对任何数值或范围的术语“约”指示允许部件或元件的集合实现如本文所述的其预期要达到的目的的合适的尺寸公差。更具体地,“约”可指列举值的值±10%的范围,例如“约90%”可指81%至99%的值范围。
应当理解,在步骤116中施加区分颜色只是指示一组边界体素具有低于预设阈值的速度的一种方法,并且假定其他区分方法(诸如将虚线覆盖在体素上)被包括在本发明的范围内。
因此应当理解,上面描述的实施方案以举例的方式被引用,并且本发明不限于上文特定示出和描述的内容。相反,本发明的范围包括上述各种特征的组合和子组合以及它们的变型和修改,本领域的技术人员在阅读上述说明时应当想到所述变型和修改,并且所述变型和修改并未在现有技术中公开。
Claims (26)
1.一种设备,包括:
导管,所述导管被配置为插入心脏的腔室中并且包括一个或多个电极,所述一个或多个电极被配置为接触所述腔室中多个位置处的心肌组织并且接收响应于行进穿过所述组织的传导波的电信号;
显示器;和
处理器,所述处理器被配置为接收来自所述导管的所述电信号,并且响应于所述电信号向所述显示器呈现所述腔室的标测图,所述腔室的标测图包括所述多个位置处的所述传导波的发生的本地时间的指示,并且响应于所述发生的本地时间来计算所述位置之间的所述传导波的速度,并且在所述标测图上标记其中所述速度低于预设阈值的、所述腔室的一个或多个区域。
2.根据权利要求1所述的设备,其中所述处理器被配置为计算在所述多个位置之间的相应位置处的内插的发生的本地时间,并且在所述相应位置处将所述内插的发生的本地时间的指示标记在所述标测图上。
3.根据权利要求2所述的设备,其中所述处理器被配置为不在所述一个或多个区域处将所述内插的发生的本地时间的所述指示标记在所述标测图上。
4.根据权利要求1所述的设备,其中所述处理器被配置为将所述腔室分成邻接表面,每个表面围绕所述多个位置中的相应一个位置,并且其中所述一个或多个区域位于两个邻接表面之间的边界处。
5.根据权利要求4所述的设备,其中所述处理器被配置为通过将所述腔室的表面形成为体素组而将所述腔室分成所述邻接表面,并且,对于所述多个位置中的每个给定位置,迭代地并且同时地将紧邻体素添加到先前体素,所述先前体素最初包括对应于所述给定位置的本地原体素,以便形成给定的邻接表面。
6.根据权利要求4所述的设备,其中所述处理器被配置为通过找到所述边界上的区与所述两个邻接表面的所述相应位置的距离的总和来计算所述区处的所述传导波的所述速度,并且将所述总和除以在所述相应位置处的所述发生的本地时间之间的时间差。
7.根据权利要求1所述的设备,其中所述电信号包括心电图(ECG)信号,并且其中所述发生的本地时间对应于从所述ECG信号推出的本地激活时间(LAT)。
8.根据权利要求7所述的设备,其中所述处理器被配置为将来自具有与所述LAT的值一一对应的相应颜色的色标的颜色结合到所述标测图中。
9.根据权利要求8所述的设备,其中所述处理器被配置为以与所述色标的颜色不同的颜色来对所述一个或多个区域着色。
10.根据权利要求1所述的设备,其中所述预设阈值介于0.01mm/ms和1mm/ms之间。
11.根据权利要求1所述的设备,其中所述预设阈值为0.1mm/ms。
12.根据权利要求1所述的设备,其中所述导管包括位置传感器,所述位置传感器被配置为向所述处理器提供指示所述多个位置的信号。
13.根据权利要求12所述的设备,其中所述位置传感器包括一个或多个线圈,所述一个或多个线圈响应于穿越所述线圈的磁场来提供所述信号。
14.一种方法,包括:
将导管插入心脏的腔室中,所述导管包括一个或多个电极,所述一个或多个电极被配置为接触所述腔室中多个位置处的心肌组织并且接收响应于行进穿过所述组织的传导波的电信号;
提供显示器;
接收来自所述导管的所述电信号;
响应于所述电信号来向所述显示器呈现所述腔室的标测图,所述腔室的标测图包括所述多个位置处的所述传导波的发生的本地时间的指示;
响应于所述发生的本地时间来计算所述位置之间的所述传导波的速度;以及
在所述标测图上标记其中所述速度低于预设阈值的、所述腔室的一个或多个区域。
15.根据权利要求14所述的方法,还包括计算在所述多个位置之间的相应位置处的内插的发生的本地时间,并且在所述相应位置处将所述内插的发生的本地时间的指示标记在所述标测图上。
16.根据权利要求15所述的方法,其中还包括不在所述一个或多个区域处将所述内插的发生的本地时间的所述指示标记在所述标测图上。
17.根据权利要求14所述的方法,还包括将所述腔室分成邻接表面,每个表面围绕所述多个位置中的相应一个位置,并且其中所述一个或多个区域位于两个邻接表面之间的边界处。
18.根据权利要求17所述的方法,还包括将所述腔室的表面形成为体素组,并且,对于所述多个位置中的每个给定位置,迭代地并且同时地将紧邻体素添加到先前体素,所述先前体素最初包括对应于所述给定位置的本地原体素,以便形成给定的邻接表面。
19.根据权利要求17所述的方法,还包括通过找到所述边界上的区与所述两个邻接表面的所述相应位置的距离的总和来计算所述区处的所述传导波的所述速度,并且将所述总和除以在所述相应位置处的所述发生的本地时间之间的时间差。
20.根据权利要求14所述的方法,其中所述电信号包括心电图(ECG)信号,并且其中所述发生的本地时间对应于从所述ECG信号推出的本地激活时间(LAT)。
21.根据权利要求20所述的方法,还包括将来自具有与所述LAT的值一一对应的相应颜色的色标的颜色结合到所述标测图中。
22.根据权利要求21所述的方法,还包括以与所述色标的颜色不同的颜色来对所述一个或多个区域着色。
23.根据权利要求14所述的方法,其中所述预设阈值介于0.01mm/ms和1mm/ms之间。
24.根据权利要求14所述的方法,其中所述预设阈值为0.1mm/ms。
25.根据权利要求1所述的方法,其中所述导管包括位置传感器,所述位置传感器被配置为提供指示所述多个位置的信号。
26.根据权利要求25所述的方法,其中所述位置传感器包括一个或多个线圈,所述一个或多个线圈响应于穿越所述线圈的磁场来提供所述信号。
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201862723218P | 2018-08-27 | 2018-08-27 | |
US62/723218 | 2018-08-27 | ||
US16/512,561 US11357437B2 (en) | 2018-08-27 | 2019-07-16 | Automatic identification of abnormal LAT differences |
US16/512561 | 2019-07-16 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110859612A true CN110859612A (zh) | 2020-03-06 |
Family
ID=67875216
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910794757.5A Pending CN110859612A (zh) | 2018-08-27 | 2019-08-27 | 自动识别异常lat差异 |
Country Status (7)
Country | Link |
---|---|
US (1) | US11357437B2 (zh) |
EP (1) | EP3616610B1 (zh) |
JP (1) | JP7475826B2 (zh) |
CN (1) | CN110859612A (zh) |
AU (1) | AU2019219785A1 (zh) |
CA (1) | CA3053045A1 (zh) |
IL (1) | IL268832B1 (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6301496B1 (en) * | 1998-07-24 | 2001-10-09 | Biosense, Inc. | Vector mapping of three-dimensionally reconstructed intrabody organs and method of display |
US20090262992A1 (en) * | 2008-04-18 | 2009-10-22 | Markowitz H Toby | Method And Apparatus For Mapping A Structure |
US20170172508A1 (en) * | 2015-12-20 | 2017-06-22 | Boston Scientific Scimed Inc. | Automatic mapping using velocity information |
US20170311833A1 (en) * | 2010-12-30 | 2017-11-02 | St. Jude Medical, Atrial Fibrillation Division, Inc. | System and method for diagnosing arrhythmias and directing catheter therapies |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2011237661B2 (en) * | 2010-04-08 | 2015-05-21 | The Regents Of The University Of California | Methods, system and apparatus for the detection, diagnosis and treatment of biological rhythm disorders |
CN107072574B (zh) * | 2014-10-15 | 2020-06-12 | 圣犹达医疗用品心脏病学部门有限公司 | 用于标测局部传导速度的方法和系统 |
-
2019
- 2019-07-16 US US16/512,561 patent/US11357437B2/en active Active
- 2019-08-21 IL IL268832A patent/IL268832B1/en unknown
- 2019-08-21 AU AU2019219785A patent/AU2019219785A1/en not_active Abandoned
- 2019-08-23 EP EP19193440.5A patent/EP3616610B1/en active Active
- 2019-08-26 JP JP2019153497A patent/JP7475826B2/ja active Active
- 2019-08-26 CA CA3053045A patent/CA3053045A1/en not_active Abandoned
- 2019-08-27 CN CN201910794757.5A patent/CN110859612A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6301496B1 (en) * | 1998-07-24 | 2001-10-09 | Biosense, Inc. | Vector mapping of three-dimensionally reconstructed intrabody organs and method of display |
US20090262992A1 (en) * | 2008-04-18 | 2009-10-22 | Markowitz H Toby | Method And Apparatus For Mapping A Structure |
US20170311833A1 (en) * | 2010-12-30 | 2017-11-02 | St. Jude Medical, Atrial Fibrillation Division, Inc. | System and method for diagnosing arrhythmias and directing catheter therapies |
US20170172508A1 (en) * | 2015-12-20 | 2017-06-22 | Boston Scientific Scimed Inc. | Automatic mapping using velocity information |
Also Published As
Publication number | Publication date |
---|---|
AU2019219785A1 (en) | 2020-03-12 |
EP3616610A1 (en) | 2020-03-04 |
US20200060568A1 (en) | 2020-02-27 |
EP3616610C0 (en) | 2023-09-27 |
IL268832A (en) | 2020-02-27 |
JP7475826B2 (ja) | 2024-04-30 |
US11357437B2 (en) | 2022-06-14 |
JP2020032185A (ja) | 2020-03-05 |
CA3053045A1 (en) | 2020-02-27 |
EP3616610B1 (en) | 2023-09-27 |
IL268832B1 (en) | 2024-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11284830B2 (en) | Detecting conduction timing | |
CN107485389B (zh) | 多通道记录中心脏激活序列的识别和可视化 | |
JP7341637B2 (ja) | 電極信号に従った電極画像の強調表示 | |
EP2749213B1 (en) | Removal of artifacts from map data | |
CN106073780B (zh) | 一种用于提供关于患者心脏的电解剖表示的信息的系统 | |
CN106901718B (zh) | 在电解剖标测图上显示多激活区域 | |
CN114650780A (zh) | 用于跟踪消融设备和生成病变线的方法、系统和装置 | |
EP3818938A1 (en) | Using statistical characteristics of multiple grouped ecg signals to detect inconsistent signals | |
RU2764590C1 (ru) | Способ изображения внутрисердечной электрограммы | |
CN112998721A (zh) | 电生理参数的双选图形呈现 | |
CN104902808B (zh) | 在解剖映射系统中估计恢复曲线 | |
CN112401901A (zh) | 电解剖标测图重新标注 | |
Punshchykova et al. | Influence of torso model complexity on the noninvasive localization of ectopic ventricular activity | |
CN110859612A (zh) | 自动识别异常lat差异 | |
CN112842353A (zh) | 标测窦性和非窦性心动周期的局部激活时间 | |
JP2020081890A (ja) | 不整脈の起点の発見 | |
JP7471914B2 (ja) | Lat体積マップ | |
JP7135071B2 (ja) | 電気循環の少なくとも1つの特異な区画の検出のために心筋壁のモデルを生成するための3次元画像のセグメンテーション方法 | |
US10806359B2 (en) | Connectivity analysis for arrhythmia drivers | |
EP4079217A1 (en) | Improved electrophysiological (ep) map coloration by considering outliers | |
Unger et al. | Estimating refractory periods during atrial fibrillation based on electrogram cycle lengths in a heterogeneous simulation setup | |
US20230190171A1 (en) | Electrophysiological (ep) map points adjustments based on user clinical interpretation | |
CN115105087A (zh) | 可视化叠加在解剖标测图上的多个参数 | |
WO2023214234A1 (en) | Detecting potential slow-conduction cardiac tissue areas in stable arrhythmias |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |