CN103549948A - 心脏激动时间检测 - Google Patents

心脏激动时间检测 Download PDF

Info

Publication number
CN103549948A
CN103549948A CN201310142932.5A CN201310142932A CN103549948A CN 103549948 A CN103549948 A CN 103549948A CN 201310142932 A CN201310142932 A CN 201310142932A CN 103549948 A CN103549948 A CN 103549948A
Authority
CN
China
Prior art keywords
signal
bipolar
time
complex wave
unipolar
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.)
Granted
Application number
CN201310142932.5A
Other languages
English (en)
Other versions
CN103549948B (zh
Inventor
M.A.埃尔哈达
G.哈亚姆
M.巴－塔尔
R.P.M.霍本
M.F.杜伊特谢弗
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.)
Biosense Webster Israel Ltd
Original Assignee
Biosense Webster Israel Ltd
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 Biosense Webster Israel Ltd filed Critical Biosense Webster Israel Ltd
Publication of CN103549948A publication Critical patent/CN103549948A/zh
Application granted granted Critical
Publication of CN103549948B publication Critical patent/CN103549948B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • 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
    • 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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7239—Details of waveform analysis using differentiation including higher order derivatives

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Physiology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明提供了一种用于表征心电图的方法，包括从心脏的第一位置接收第一单极信号以及从所述心脏的第二位置接收第二单极信号。所述方法还包括由所述第一单极信号和第二单极信号生成双极信号，以及分析所述双极信号以描绘所述第一位置和第二位置生成双极复合波的时间周期。所述方法还包括分析所述时间周期内的所述第一单极信号以确定所述第一位置的激动时间。

Description

心脏激动时间检测
技术领域
本发明一般涉及信号分析，并且具体地讲涉及对心脏搏动所生成的信号的分析。
背景技术
用于表征心脏活动的方法中的一种依赖于对心脏搏动时由心脏生成的电信号的分析。这些信号通常具有相对低的毫伏级的电平，使得对信号的精确分析可能比较困难。虽然存在困难，但是精确分析可实现对心脏活动的改善表征，包括确定可能有缺陷的心脏区域。
以引用方式并入本专利申请的文献将视为本专利申请的整体部分，但是如果这些并入的文献中定义任何术语与本说明书中明确或隐含地给出的定义相冲突，则应只考虑本说明书中的定义。
发明内容
本发明的一个实施例提供了用于表征心电图的方法，包括：
从心脏的第一位置接收第一单极信号以及从心脏的第二位置接收第二单极信号；
由第一单极信号和第二单极信号生成双极信号；
分析双极信号以描绘第一位置和第二位置生成双极复合波的时间周期；以及
分析时间周期内的第一单极信号以确定第一位置的激动时间。
通常，分析双极信号包括确定将被应用于双极信号的搜索窗口界限。分析第一单极信号可包括将搜索窗口界限应用于第一单极信号。
在一个本发明所公开的实施例中，描绘时间周期包括将双极信号的数据馈送到双态状态机，以便确定时间周期的界限。
在又一个本发明所公开的实施例中，分析双极信号包括对双极信号的数据分类以确定双极复合波的阈电平。
在又一个本发明所公开的实施例中，分析双极信号包括区分然后修正双极信号的数据，以便生成差异化数据。描绘时间周期可包括将差异化数据馈送到四态状态机，以便确定时间周期的界限。确定激动时间可包括形成第一单极信号的一阶导数，并且将单极起始激动时间指定为其中一阶导数为最小值的时间点。
在一个可供选择的实施例中，激动时间包括第一激动时间，并且所述方法还包括分析时间周期内的第二单极信号以确定第二位置的第二激动时间。
在又一个可供选择的实施例中，双极复合波包括第一双极复合波和第二双极复合波，并且时间周期包括生成第一双极复合波的第一时间周期以及生成第二双极复合波的第二时间周期，并且分析第一单极信号包括分别确定第一时间周期和第二时间周期内的第一激动时间和第二激动时间。
根据本发明的一个实施例，还提供了用于表征心电图的设备，包括：
探针，其被配置成从心脏的第一位置接收第一单极信号以及从心脏的第二位置接收第二单极信号；和
处理器，其被配置成：
由第一单极信号和第二单极信号生成双极信号，
分析双极信号以描绘第一位置和第二位置生成双极复合波的时间周期，以及
分析时间周期内的第一单极信号以确定第一位置的激动时间。
根据本发明的一个实施例，还提供了用于表征心电图的计算机软件产品，包括其中存储有计算机程序指令的有形计算机可读介质，该指令在被计算机读取时，使计算机：
从心脏的第一位置接收第一单极信号以及从心脏的第二位置接收第二单极信号；
由第一单极信号和第二单极信号生成双极信号；
分析双极信号以描绘第一位置和第二位置生成双极复合波的时间周期；以及
分析时间周期内的第一单极信号以确定第一位置的激动时间。
结合附图，通过以下对本发明实施例的详细说明，将更全面地理解本发明。
附图说明
图1是根据本发明的实施例的激动时间检测系统的示意图；
图2是根据本发明的实施例示出了操作系统的整体过程的示意性方框图；
图3是根据本发明的实施例示出了搜索窗口模块的示意性方框图；
图4是根据本发明的实施例示出了用于搜索窗口参数块中的参数之间的关系的时间线；
图5A是根据本发明的实施例示出了在第一相块中执行的第一组动作的示意性方框图，并且图5B和5C是根据本发明的实施例的动作前后数据的示意性电压与时间曲线图；
图6根据本发明的实施例示出了由过滤块输出的窗口化平滑数据；
图7A是根据本发明的实施例示出了在第一相块中执行的第二组动作的示意性方框图，并且图7B是根据本发明的实施例的由动作产生的数据的示意性曲线图；
图8是根据本发明的实施例的四态状态机的示意图；
图9根据本发明的实施例示出了状态机的操作；
图10根据本发明的实施例示出了在时间线上绘制的时间实例的值；
图11A和11B是根据本发明的实施例的示意性双极和单极曲线图；
图12A、12B和12C是根据本发明的实施例由一次心脏搏动内发生的多个双极复合波衍生的信号曲线图；并且
图13是根据本发明的实施例确定激动时间所遵循的步骤的流程图。
具体实施方式
综述
本发明的一个实施例通过在两个阶段处理心电图数据提供了用于表征心电图的方法。数据的形式为来自心脏中两个不同位置的两个单极信号，并且表征能够确定提供数据的心脏中的位置的激动时间。
在方法的第一阶段中，将数据作为双极信号分析，以确定描绘信号内双极复合波的信号的时间实例。在方法的第二阶段中，将时间实例用作可在其中单独分析单极信号每一个的界限。
为了确定不同位置的激动时间，对单极信号每一个的一阶导数进行评估。假设一阶导数最小处的时间为起始激动时间，即生成单极信号的组织开始激活时的时间。该方法可用于找到两个不同位置中的每一个的起始激动时间。
该方法可用于分析每次心脏搏动具有一个双极复合波的信号，还可用于分析每次心脏搏动具有多于一个双极复合波的信号。
发明人已实时操作了本方法，并在临床上证实了本方法可提供精确的结果。
对系统的描述
现在参考图1，其为根据本发明的实施例的激动时间检测系统10的示意图。系统10分析心电图信号，以便尤其测量在给定信号的时间内的起始点。为了简洁和清楚起见，除非另外指明，否则以下说明均假设其中系统10使用探针14对心脏12(本文假设包括人类心脏)进行测量的研究程序。
通常，探针14包括导管，该导管在研究程序的过程中插入受检者16体内。探针的远侧末端18包括从心脏12中的相应位置23和25接收心电图(ECG)信号的第一电极20和第二电极21。所述位置通常在心脏的组织27内。如本文所述，来自两个电极的信号形成由系统10分析的双极信号。研究程序由系统10的使用者22执行，并且在本文的描述中，以举例的方式假设使用者22为医疗专业人员。
在程序的过程中使用一个或多个其他电极29。其他电极可附接到探针14、附接到类似于探针14并位于心脏内的另一探针、和/或附接到受检者16的皮肤。该其他电极用作参考电极，以提供来自电极20和21的信号的参考接地，在此情况下，相应电极的两个信号为单极信号。
系统10通常由系统处理器24控制，该系统处理器可被实现为通用计算机。系统处理器包括与存储器28通信的处理单元26。处理器24可安装在控制台30中，控制台包括操作控制器32，该操作控制器通常包括专业人员22用于与处理器互动的小键盘和定位装置，例如鼠标或轨迹球。处理器24所执行的操作的结果在屏幕34上提供给专业人员，该屏幕可显示系统所执行的分析的结果的图36。作为另外一种选择或除此之外，系统10将结果用于向专业人员22展示其他参数，例如心脏12的局部激动时间(LAT)的标测图。在系统10的操作中，专业人员22能够使用控制器32来输入处理器24所使用的参数的值。
处理器24使用存储器28中存储的软件来操作系统10。软件可以电子形式通过网络下载至处理器24，例如作为另外一种选择或除此之外，软件可以提供于和/或存储在非临时性有形计算机可读介质(例如，磁存储器、光学存储器或电子存储器)上。
系统10可实现为购自Biosense Webster，Inc.(3333 Diamond CanyonRoad，Diamond Bar，CA 91765)的CARTO XP EP导航和消融系统，并经适当的修改以执行本文所述的程序。
在一些情况下，电极20和/或21可提供ECG和其他信号，或者电极可用于其他目的。例如，上文提及的CARTO系统使用检测ECG信号的电极、测量电极的阻抗以进行跟踪，以及使用电极提供射频消融。
图2是根据本发明的实施例示出了处理器24在操作系统10时所遵循的总体过程的示意性方框图70。在双极阶段72中，处理器从电极20和21接收原始未过滤的信号作为电压电平，并对它们进行操作以形成双极信号数据。处理器对双极数据进行分析，以确定限定双极复合波的时间周期或窗口。为了简洁和清楚起见，在以下描述中，除非另外指明，否则均假设每次心脏搏动中存在一个双极复合波。
双极复合波由初始时间实例T起始和最终时间实例T终止界定。处理器使用双极复合波的时间界限来限定将在其中执行单极分析的窗口。
在单极阶段74中，处理器将电极20和21信号的每一个单独地视为单极电压与时间信号，并分析在双极阶段中发现的时间窗口内的单极信号。分析使得处理器能够确定与电极20和21接触的区域激活时的相应单极激动时间。激动时间通常包括单极信号的导数具有最大负值时的时间。
双极阶段72由三个模块形成：搜索窗口模块76和两个后续模块，即第一相模块78和第二相模块80。处理器针对每个模块所执行的操作在下文描述。在该描述中，假设来自电极20和21的信号在大约2.5s的周期内以大约1kHz的频率采样，从而得到大约2,500个将由系统10分析的样本。然而，系统10可通过任何合适的样本周期和采样率进行操作。
图3是根据本发明的实施例更详细地示出了搜索窗口模块76的示意性方框图。在R波检测块90中，处理器24分析一组来样的值，以识别样本中R波出现时的时间。通常，对于一组在2.5s内采集的样本，存在大约二至四个R波，但是患有心动过速的受检者在2.5s的时间周期内可具有五个或更多个R波。该识别过程通常通过发现样本峰值出现时的时间而执行。
在RR间隔块92中，处理器发现在块92中识别的峰值之间的平均时间周期
在搜索窗口参数块94中，处理器计算搜索窗口的开始时间(SW开始)和结束时间(SW结束)，以用于输入数据的进一步分析。在上文提及的CARTO系统中，专业人员22能够对感兴趣窗口(WOI)的中心时间(WOI中心)和宽度(WOI宽度)进行编程。为了在CARTO系统中执行计算，块94使用由专业人员22提供的参数WOI中心、WOI宽度以及附加时间周期WOIDELTA(本文也使用符号Δ提及)的值。WOI中心通常由专业人员专门设定，以逼近平均时间周期
Figure BDA00003089057200062
的预期中间时间点，但也可将WOI中心设定为任何其他合适的时间点。WOI宽度通常也由专业人员专门设定，以逼近预期平均时间周期
Figure BDA00003089057200063
但也可将其设定为任何合适的时间周期。使用WOI中心、WOI宽度和WOIΔ的值，块94计算搜索窗口的SW开始、SW结束的值。
图4是根据本发明的实施例示出了用于搜索窗口参数块94中的参数之间的关系的时间线。如时间线所示，由块94描绘的搜索窗口具有
Figure BDA00003089057200064
的总宽度，其在SW开始时间处开始并在SW结束时间处结束。
应当理解，虽然已结合CARTO系统阐述了由块94生成的搜索窗口的开始和结束时间的计算，但是专业人员22可以使用本领域已知的任何合适的方法描绘合适的搜索窗口。
Δ的典型值为大约20ms。
Figure BDA00003089057200071
的典型值取决于受检者16。对于心动过速的受检者，
Figure BDA00003089057200072
可为大约240ms，在此情况下，Δ具有20ms的值，搜索窗口为大约220ms宽。
图5A是根据本发明的实施例示出了在第一相块78中由处理器24执行的第一组动作的示意性方框图，并且图5B和5C是根据本发明的实施例的动作前后数据的示意性电压与时间曲线图。(为了简洁起见，未示出曲线图的电压和时间轴。)在修正和过滤块120中，首先对来自电极20和21并在图5B中示出的双极原始数据进行修正，然后低通过滤以移除数据中的高频分量并产生平滑数据。在一个实施例中，本发明人使用截止频率为大约20Hz的二阶巴特沃斯滤波器。
然后使用块94(图3)中的搜索窗口时间SW开始和SW结束使过滤后的平滑数据窗口化，以生成一组样本数据{X(n)}，其中n为数据的指数，X为数据值。这组平滑数据在图5C中示意性地示出。假设上文针对心动过速的受检者给出的示例搜索窗口宽度和大约1000Hz的采样率，在窗口化数据中存在大约220个平滑样本，使得在此情况下n为介于1和大约220之间的正整数。
在分类块122中，按值对平滑样本分类，并将其布置成频率分布。从频率分布中，提取出将被应用于数据分析的阈电压电平THR。对电平THR进行选择，以接近但高于平滑基线数据的电平。在一个实施例中，将电平选择为对应于频率分布的第5个百分率的基础值，加上一个平滑信号幅度5％的因子。作为另外一种选择，可通过任何其他合适的方法选择电平THR，以限定接近但高于平滑基线数据的电平。
此外，分类块122确定平滑数据的峰值样本X(np1)。
处理器向双态状态机124提供电平THR和采样的平滑值X(n)。状态机两状态A和B之间的过渡的条件在图5A的方括号[]内指定；在过渡期间执行的动作在大括号{}内指定。从峰值样本X(np1)开始，将数据X(n)及时按顺序后馈，直到第一过渡在指数underTHR开始处发生。此外，从峰值样本X(np1)开始将数据及时前馈，直到第二过渡在指数underTHR结束处发生。参数cnt对状态机所操作的样本的数量计数。由使用者设定的变量CNTMAX(表示在underTHR开始和underTHR结束两过渡之间的样本的可接受数量)通常设定为大约100，但也可以设定为任何其他合适的数量。
图6根据本发明的实施例示出了由过滤块120(也如图5C中示意性地示出)输出的窗口化平滑数据。曲线图130代表由过滤块输出的窗口化平滑样本X(n)。状态机124将样本分成三个部分：低于阈THR的两个基线部分132和134，以及双极复合波部分136。双极复合波由状态机生成的两个过渡指数underTHR开始和underTHR结束界定。
图7A是根据本发明的实施例示出了处理器24在第一相块78(图2)中执行的第二组动作的示意性方框图，并且图7B是根据本发明的实施例由动作所产生的数据的示意曲线图。在过滤块150中，将得自电极20和21的双极原始数据进行低通过滤，以移除高频分量并产生平滑数据。在一个实施例中，发明人使用截止频率为大约35Hz的二阶巴特沃斯滤波器。在区分块152中，对平滑数据进行区分，然后在修正块154中修正，以产生经修正的差异化数据。
使用得自块94(图3)的搜索窗口时间SW开始和SW结束，将得自块154的数据在窗口块156中进行窗口化。窗口化生成一组差异化平滑数据{D(n)}，其中D为数据值。图7B是块154的数据输出的图示，并在图9中更详细地示出。
在双极阶段的第二相80(图2)中，这组差异化平滑数据传输到分类块158以及四态状态机160。在分类块158中，由双态状态机124确定并在图6中示出的指数underTHR开始和underTHR结束用于将{D(n)}分成差异化二元复合波部分和两个噪声部分。处理器24将两个噪声部分中的值分类成频率分布，并从该分布中提取出将被应用于分析差异化平滑数据的差异化噪声电平噪声。将电平噪声选择为接近但高于两个噪声部分的电平，并在图7B中示意性地示出。在一个实施例中，该电平基于频率分布的第95个百分率。
分类块158还确定差异化二元复合波的峰值D(np2)和指数np2，并将D(np2)传输到四态状态机。
图8是根据本发明的实施例的四态状态机160的示意图。状态机包括四个状态A、B、C和D，以及两个退出状态E和F。状态之间过渡的条件在图8的方括号[]中指定；在过渡期间执行的动作在大括号{}内指定。从峰值样本D(np2)开始，并且当状态机处于状态A时，将样本数据D(n)及时后馈，直到达到退出状态F。达到状态F时的时间(即指数值)为双极复合波的起始时间T起始。此外，通过将样本数据D(n)及时前馈直到达到退出状态E，找到双极复合波的终止时间T终止。
在状态机中，参数cnt和gcnt对状态机所操作的样本的数量计数。变量CNT状态2、CNT状态3和CNT状态4可由专业人员22设定，它们代表通过差异化噪声电平噪声出现过渡时的状态机的状态之间的样本的可接受数量。CNT状态2、CNT状态3和CNT状态4的典型值分别为8、18和4，但是这些值可由专业人员22设定为任何合适的值。
图9根据本发明的实施例示出了状态机160的操作。曲线图170(与图7B相似)代表从窗口块156传输到状态机的平滑数据D(n)。从分类块158传输的噪声电平噪声和峰D(np2)的值也显示在曲线图170上。
曲线图172显示了在确定T起始值时的状态机的状态以及状态之间的过渡。如曲线图中所示，处理器24(图1)开始在状态A中在样本np2处从峰值D(np2)操作状态机。随着后续的及时向后样本馈送到状态机，机器最初在状态A和B之间交替后继而转到状态C、D、A、B和C。在最后一个状态C，机器转到退出状态F(图8)。对于及时向前馈送的样本，一组相似的过渡也从峰值D(np2)开始，在状态D和退出状态E结束，并确定T终 止的值。
图10根据本发明的实施例示出了在时间线上绘制的T起始和T终止的值。时间线示出了T起始和T终止值之间的典型关系，以及用于研究双极复合波和上文结合图4所述的时间值。
通过T起始和T终止的值，系统10能够根据公式(1)评估双极复合波的信噪比(SNR)：
SNR = 20 · log ( S - N N ) - - - ( 1 )
其中S是位于T起始和T终止之间的未过滤双极数据的均方根(RMS)值，并且
N是T起始之前和T终止之后的未过滤双极数据的RMS值。
专业人员22能够使用SNR的值以便确立T起始和T终止的评估值的置信度。
回到图2，处理器24将T起始和T终止的值传输到单极阶段74。在阶段74中，处理器形成由T起始和T终止界定的时间窗口，并分析窗口内得自电极20和21每一个的平滑单极电压(V)与时间(t)信号。在窗口内，处理器计算每个单极信号的斜率的值，即一阶导数
Figure BDA00003089057200101
的值。对于每个信号，处理器选择一阶导数
Figure BDA00003089057200102
具有其最负值(即，其最小值)时的时间，并将该时间假设为生成信号的组织开始激活时的时间。
图11A和11B是根据本发明的实施例的示意性双极和单极曲线图。曲线图180是双极信号的电压与时间曲线图，而曲线图182和184是形成双极信号的相应单极信号的电压与时间曲线图。两套曲线图均具有如上确定的标记在曲线图上的时间T起始和T终止。就曲线图182和184而言，示出了相应的激动时间186和188，它们是由T起始和T终止限定的窗口内的相应单极信号的最负导数的时间。激动时间186和188是生成单极信号的组织开始激活的时间，在本文也称为单极起始激动时间。
为清楚起见，上文的描述考虑的是评估信号参数的系统10的实施例，其中每次心脏搏动存在一个双极复合波。系统10不限于此类评估，并可用于识别其中每次心脏搏动出现多个双极复合波的信号，此外还可用于评估多个双极复合波的信号参数。识别多个双极信号的出现通常可以通过测量相邻复合波之间的间隔，因为与每次心脏搏动具有一个双极复合波的信号相比，该间隔将发生变化。
本领域的普通技术人员将能够对上文的说明加以必要的变更，以评估每次心脏搏动出现多个双极复合波的单极信号的参数。此类参数包括但不限于评估给定心脏搏动中每个双极复合波的相应单极起始激动时间。
图12A、12B和12C是根据本发明的实施例由一次心脏搏动内发生的多个双极复合波衍生的信号曲线图。曲线图190(图12A)是表现出心房双极复合波194以及心室双极复合波192和196的双极信号。每个双极复合波可通过首先限定给定复合波的搜索窗口进行分析。限定每个复合波的搜索窗口的方法基本如上文结合图3所述，加以必要的变更以允许区分双极信号内的RR间隔。
曲线图200(图12B)是具体心室双极复合波192的放大曲线图。复合波的起始时间202和终止时间204已标记在曲线图上。通过将来源于复合波的平滑数据馈送到状态机160，基本如上文结合图8所述对时间进行评估。
曲线图210(图12C)示出了对应于图12B的双极复合波192的单极信号212和214。如上所述，各信号的相应单极起始激动时间216和218在其中各信号的一阶导数(在起始时间202和终止时间204之间测量)具有其最负值(即，最小值)时的时间发生。
系统10还可用于评估与在一次心脏搏动内发生的具有多个双极复合波的信号相关的其他参数，这对本领域的普通技术人员将显而易见。此类参数包括但不限于度量复合波之间的平均RR间隔的第一心房和第二心房双极复合波之间的持续时间。假设所有此类参数均包括在本发明的范围内。
图13是根据本发明的实施例由处理器24在操作系统10中用以确定激动时间时所遵循的步骤的流程图250。为了简洁和清楚起见，除非另外指明，否则流程图步骤的描述假设接收到的信号在每次心脏搏动具有一个双极复合波。本领域的普通技术人员将能够针对每次心脏搏动具有多个双极复合波的情况对这些描述进行调整。
步骤252至260是在双极阶段72中执行的动作，步骤262在单极阶段74(图2)。
在初始步骤252中，处理器从电极20和21接收作为采样数据的信号。处理器对信号进行分析，以识别R波、值和搜索窗口的界限，如上文结合图3和4所述。
在第一过滤步骤254中，对采样数据进行修正、过滤和窗口化，并将所得的平滑数据馈送到双态状态机124。在划分步骤256中，双态状态机将其接收到的数据分成基线部分和双极复合波部分。步骤254和256在上文结合图5A-5C和图6进行了描述。
在第二过滤步骤258中，对双极复合波的采样数据进行过滤、区分和窗口化，以导出第二平滑信号，如上文结合图7A和7B所述。
在双极复合波分析步骤260中，处理器通过将第二平滑信号数据馈送到四态状态机160而评估复合波的起始和终止时间，如结合图8和9所述。
在激动时间步骤262中，通过分析在由步骤260的双极起始和终止时间限定的窗口内得自各电极的单极信号，确定与电极20和21接触的组织的激动时间。在步骤262中由处理器执行的动作结合图11A和11B以及(对于在一次心脏搏动中有多个双极复合波的情形)结合图12A-12C进行描述。
分析对窗口内的单极信号进行区分，并找到一阶导数最负(即，最小)时的相应时间。这些时间对应于与电极20接触的组织的起始激动时间以及与电极21接触的组织的起始激动时间。
应当理解，上述实施例仅以举例的方式引用，并且本发明并不限于上面具体示出和描述的内容。相反，本发明的范围包括上述各种特征的组合和子组合以及它们的变型和修改形式，本领域的技术人员在阅读上述说明时将会想到所述变型和修改形式，并且所述变型和修改形式并未在现有技术中公开。

Claims (20)

1.一种用于表征心电图的方法，包括：
从心脏的第一位置接收第一单极信号以及从所述心脏的第二位置接收第二单极信号；
由所述第一单极信号和第二单极信号生成双极信号；
分析所述双极信号以描绘所述第一位置和第二位置生成双极复合波的时间周期；以及
分析所述时间周期内的所述第一单极信号以确定所述第一位置的激动时间。
2.根据权利要求1所述的方法，其中分析所述双极信号包括确定将被应用于所述双极信号的搜索窗口界限。
3.根据权利要求2所述的方法，其中分析所述第一单极信号包括将所述搜索窗口界限应用于所述第一单极信号。
4.根据权利要求1所述的方法，其中描绘所述时间周期包括将所述双极信号的数据馈送到双态状态机，以便确定所述时间周期的界限。
5.根据权利要求1所述的方法，其中分析所述双极信号包括对所述双极信号的数据分类，以确定所述双极复合波的阈电平。
6.根据权利要求1所述的方法，其中分析所述双极信号包括区分然后修正所述双极信号的数据，以便生成差异化数据。
7.根据权利要求6所述的方法，其中描绘所述时间周期包括将所述差异化数据馈送到四态状态机，以便确定所述时间周期的界限。
8.根据权利要求7所述的方法，其中确定所述激动时间包括形成所述第一单极信号的一阶导数，并且将单极起始激动时间指定为其中所述一阶导数为最小值的时间点。
9.根据权利要求1所述的方法，其中所述激动时间包括第一激动时间，所述方法还包括分析所述时间周期内的所述第二单极信号以确定所述第二位置的第二激动时间。
10.根据权利要求1所述的方法，其中所述双极复合波包括第一双极复合波和第二双极复合波，并且其中所述时间周期包括生成所述第一双极复合波的第一时间周期以及生成所述第二双极复合波的第二时间周期，并且其中分析所述第一单极信号包括分别确定所述第一时间周期和第二时间周期内的第一激动时间和第二激动时间。
11.一种用于表征心电图的设备，包括：
探针，所述探针被配置成从心脏的第一位置接收第一单极信号以及从所述心脏的第二位置接收第二单极信号；和
处理器，所述处理器被配置成：
由所述第一单极信号和第二单极信号生成双极信号，
分析所述双极信号以描绘所述第一位置和第二位置生成双极复合波的时间周期，以及
分析所述时间周期内的所述第一单极信号以确定所述第一位置的激动时间。
12.根据权利要求11所述的设备，其中分析所述双极信号包括确定将被应用于所述双极信号的搜索窗口界限。
13.根据权利要求11所述的设备，其中描绘所述时间周期包括将所述双极信号的数据馈送到双态状态机，以便确定所述时间周期的界限。
14.根据权利要求11所述的设备，其中分析所述双极信号包括对所述双极信号的数据分类，以确定所述双极复合波的阈电平。
15.根据权利要求11所述的设备，其中分析所述双极信号包括区分然后修正所述双极信号的数据，以便生成差异化数据。
16.根据权利要求15所述的设备，其中描绘所述时间周期包括将所述差异化数据馈送到四态状态机，以便确定所述时间周期的界限。
17.根据权利要求16所述的设备，其中确定所述激动时间包括形成所述第一单极信号的一阶导数，并且将单极起始激动时间指定为其中所述一阶导数为最小值的时间点。
18.根据权利要求16所述的设备，其中所述激动时间包括第一激动时间，并且还包括分析所述时间周期内的所述第二单极信号以确定所述第二位置的第二激动时间。
19.根据权利要求11所述的设备，其中所述双极复合波包括第一双极复合波和第二双极复合波，并且其中所述时间周期包括生成所述第一双极复合波的第一时间周期和生成所述第二双极复合波的第二时间周期，并且其中分析所述第一单极信号包括分别确定所述第一时间周期和第二时间周期内的第一激动时间和第二激动时间。
20.一种用于表征心电图的计算机软件产品，包括有形计算机可读介质，在所述有形计算机可读介质中存储有计算机程序指令，所述指令在被计算机读取时，使所述计算机：
从心脏的第一位置接收第一单极信号以及从所述心脏的第二位置接收第二单极信号；
由所述第一单极信号和第二单极信号生成双极信号；
分析所述双极信号以描绘所述第一位置和第二位置生成双极复合波的时间周期；以及
分析所述时间周期内的所述第一单极信号以确定所述第一位置的激动时间。
CN201310142932.5A 2012-04-23 2013-04-23 心脏激动时间检测 Active CN103549948B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US13/453,249 US9833157B2 (en) 2012-04-23 2012-04-23 Cardiac activation time detection
US13/453249 2012-04-23
US13/453,249 2012-04-23

Publications (2)

Publication Number Publication Date
CN103549948A true CN103549948A (zh) 2014-02-05
CN103549948B CN103549948B (zh) 2018-09-04

Family

ID=48190721

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310142932.5A Active CN103549948B (zh) 2012-04-23 2013-04-23 心脏激动时间检测

Country Status (7)

Country Link
US (1) US9833157B2 (zh)
EP (2) EP2656784B1 (zh)
JP (1) JP6309205B2 (zh)
CN (1) CN103549948B (zh)
AU (1) AU2013204274B2 (zh)
CA (1) CA2813586A1 (zh)
IL (1) IL225773B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107714031A (zh) * 2016-08-11 2018-02-23 韦伯斯特生物官能(以色列)有限公司 波前的注释
CN107714025A (zh) * 2016-08-11 2018-02-23 韦伯斯特生物官能(以色列)有限公司 分类ecg信号

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10758141B2 (en) * 2012-04-23 2020-09-01 Biosense Webster (Israel) Ltd. Cardiac activation time detection
WO2015066113A1 (en) 2013-10-30 2015-05-07 St. Jude Medical, Cardiology Division, Inc. Cardiac mapping system and method for bi-directional activation detection of electrograms
JP6273355B2 (ja) 2013-10-30 2018-01-31 セント・ジュード・メディカル，カーディオロジー・ディヴィジョン，インコーポレイテッド 電位図の電圧ベースの評価のための心臓マッピングのシステム及び方法
US9554718B2 (en) * 2014-01-29 2017-01-31 Biosense Webster (Israel) Ltd. Double bipolar configuration for atrial fibrillation annotation
US9380953B2 (en) * 2014-01-29 2016-07-05 Biosense Webster (Israel) Ltd. Hybrid bipolar/unipolar detection of activation wavefront
US9833161B2 (en) 2015-02-09 2017-12-05 Biosense Webster (Israel) Ltd. Basket catheter with far-field electrode
CN107837083B (zh) * 2017-10-31 2019-05-10 太原理工大学 基于最小二乘支持向量机的j波自动检测方法
CN112004463B (zh) 2018-04-26 2024-03-22 圣犹达医疗用品心脏病学部门有限公司 用于标测心律失常驱动器位点的系统和方法
US10932686B2 (en) 2018-05-22 2021-03-02 Biosense Webster (Israel) Ltd. Identifying activations in an atrial fibrillation electrogram
JP7285958B2 (ja) * 2019-05-24 2023-06-02 セント・ジュード・メディカル，カーディオロジー・ディヴィジョン，インコーポレイテッド 心臓マッピングするためのシステム及び方法
US11311226B2 (en) 2019-12-12 2022-04-26 Biosense Webster (Israel) Ltd. Detection of ventricular activity using unipolar and bipolar signals

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1647301A1 (en) * 2004-10-12 2006-04-19 Pacesetter, Inc. Mode switching heart stimulation apparatus and method
US7218960B1 (en) * 2003-06-24 2007-05-15 Pacesetter, Inc. System and method for detecting cardiac ischemia based on T-waves using an implantable medical device
US20100274150A1 (en) * 2009-04-23 2010-10-28 Rhythmia Medical, Inc. Multi-Electrode Mapping System

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6399840A (ja) 1986-10-17 1988-05-02 テルモ株式会社 生体信号計測装置
US5778881A (en) * 1996-12-04 1998-07-14 Medtronic, Inc. Method and apparatus for discriminating P and R waves
US7225015B1 (en) 2003-06-24 2007-05-29 Pacesetter, Inc. System and method for detecting cardiac ischemia based on T-waves using an implantable medical device
US20050059897A1 (en) 2003-09-17 2005-03-17 Snell Jeffery D. Statistical analysis for implantable cardiac devices
US7194300B2 (en) 2004-01-21 2007-03-20 Cardionet, Inc. Cardiac monitoring
JP2008523929A (ja) 2004-12-21 2008-07-10 シドニー ウエスト エリア ヘルス サービス 電気生理学データの自動処理
US7672722B1 (en) 2005-12-21 2010-03-02 Pacesetter, Inc. Hardware-based state machine for use in discriminating near field signals from far field signals for use in an implantable cardiac stimulation device
US7912535B2 (en) 2006-03-02 2011-03-22 University Of Rochester Method and system for assessing repolarization abnormalities
JP2008028830A (ja) * 2006-07-24 2008-02-07 General Res Of Electronics Inc 位相比較信号処理回路
GB0624081D0 (en) 2006-12-01 2007-01-10 Oxford Biosignals Ltd Biomedical signal analysis method
US8676303B2 (en) 2008-05-13 2014-03-18 The Regents Of The University Of California Methods and systems for treating heart instability
CA2739838C (en) 2008-10-09 2016-11-01 The Regents Of The University Of California Methods, system and apparatus for the detection, diagnosis and treatment of biological rhythm disorders
US8326419B2 (en) 2009-04-07 2012-12-04 Pacesetter, Inc. Therapy optimization via multi-dimensional mapping
EP2453975B1 (en) * 2009-07-15 2016-11-02 Cardiac Pacemakers, Inc. Remote sensing in an implantable medical device
US8521275B2 (en) * 2010-04-28 2013-08-27 Medtronic, Inc. Method of dual EGM sensing and heart rate estimation in implanted cardiac devices
US8700136B2 (en) * 2011-11-11 2014-04-15 Biosense Webster (Israel), Ltd. Accurate time annotation of intracardiac ECG signals

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7218960B1 (en) * 2003-06-24 2007-05-15 Pacesetter, Inc. System and method for detecting cardiac ischemia based on T-waves using an implantable medical device
EP1647301A1 (en) * 2004-10-12 2006-04-19 Pacesetter, Inc. Mode switching heart stimulation apparatus and method
US20100274150A1 (en) * 2009-04-23 2010-10-28 Rhythmia Medical, Inc. Multi-Electrode Mapping System

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
GJIN NDREPEPA M D D PH ET AL: "Activation Time Determination by High-Resolution Unipolar and Bipolar Extracellular Electrograms in the Canine Heart Activation Time Determination by High一Resolution Electrograms.Introduction", 《JOURNAL OF CARDIOVASCULAR ELECTROPHYSIOLOGY》 *
周益锋等: "单极和双极心内电图的比较", 《实用医学杂志》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107714031A (zh) * 2016-08-11 2018-02-23 韦伯斯特生物官能(以色列)有限公司 波前的注释
CN107714025A (zh) * 2016-08-11 2018-02-23 韦伯斯特生物官能(以色列)有限公司 分类ecg信号
CN107714031B (zh) * 2016-08-11 2021-12-24 韦伯斯特生物官能(以色列)有限公司 波前的注释
CN107714025B (zh) * 2016-08-11 2022-06-07 韦伯斯特生物官能(以色列)有限公司 分类ecg信号

Also Published As

Publication number Publication date
EP2656784B1 (en) 2022-06-22
CA2813586A1 (en) 2013-10-23
EP3925535A1 (en) 2021-12-22
US9833157B2 (en) 2017-12-05
EP2656784A1 (en) 2013-10-30
CN103549948B (zh) 2018-09-04
IL225773B (en) 2018-11-29
AU2013204274A1 (en) 2013-11-07
JP6309205B2 (ja) 2018-04-11
JP2013223730A (ja) 2013-10-31
AU2013204274B2 (en) 2014-10-09
IL225773A0 (en) 2013-09-30
US20130281870A1 (en) 2013-10-24

Similar Documents

Publication Publication Date Title
CN103549948A (zh) 心脏激动时间检测
JP6675467B2 (ja) 活性化波面のバイポーラ／ユニポーラハイブリッド検出のための装置
AU2012247057B2 (en) Accurate time annotation of intracardiac ecg signals
CN107837082A (zh) 基于人工智能自学习的心电图自动分析方法和装置
CN101828918B (zh) 基于波形特征匹配的心电信号r波峰检测方法
CN107714025B (zh) 分类ecg信号
JP7002881B2 (ja) 波面のアノテーション
JP2013223730A5 (zh)
EP3281579B1 (en) Identifying ecg signals having the same morphology
CN103300848B (zh) 控制用于检测生理信号的顶峰的阈值的设备和方法
CN111345814B (zh) 一种心电信号中心拍的分析方法、装置、设备和存储介质
US10758141B2 (en) Cardiac activation time detection
CN111345815A (zh) 一种心电信号中qrs波的检测方法、装置、设备和存储介质
CN111685754B (zh) 针对可穿戴ecg采集设备的心率计算方法及系统
Gallón et al. Development of a low-cost sEMG Acquisition and processing platform
Silva Filho et al. A new parameter in the nonlinear dynamics of the heart: the higher reconstruction step
CN113842127A (zh) 脉率检测装置、系统及存储介质

Legal Events

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