US11071491B2 - Methods and systems for mapping cardiac repolarization - Google Patents
Methods and systems for mapping cardiac repolarization Download PDFInfo
- Publication number
- US11071491B2 US11071491B2 US15/766,283 US201615766283A US11071491B2 US 11071491 B2 US11071491 B2 US 11071491B2 US 201615766283 A US201615766283 A US 201615766283A US 11071491 B2 US11071491 B2 US 11071491B2
- Authority
- US
- United States
- Prior art keywords
- time
- tick
- repolarization
- tick time
- slm
- 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, expires
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
- 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
-
- 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
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2505/00—Evaluating, monitoring or diagnosing in the context of a particular type of medical care
- A61B2505/05—Surgical care
Definitions
- the instant disclosure relates to electrophysiological mapping, such as may be performed in cardiac diagnostic and therapeutic procedures.
- the instant disclosure relates to systems, apparatuses, and methods for mapping cardiac repolarization activity, including action potential duration (“APD”) and activation recovery interval (“ARI”).
- APD action potential duration
- ARI activation recovery interval
- a method of mapping cardiac repolarization activity including the following steps: receiving an electrogam signal S(t) at a signal processor; and using the signal processor: computing a second derivative S′′(t) of the electrogram signal S(t) for a preset window about a reference time T ref ; identifying a depolarization tick time T D within the preset window using S′′(t); identifying a repolarization tick time T R using S(t), wherein the repolarization tick time T R is identified as a point in time when S(t) passes below a preset threshold; and computing an action potential duration for the preset window using the depolarization tick time T D and the repolarization tick time T R .
- the electrogram signal S(t) can be a bipolar electrogram signal from a monophasic action potential (“MAP”) catheter.
- MAP monophasic action potential
- the reference time T ref can be detected using a reference cardiac signal.
- the step of computing a second derivative S′′(t) of the electrogram S(t) for a preset window about a reference time T ref includes computing a first derivative S′(t) of the electrogram S(t) for the preset window about the reference time T ref according to a formula
- the step of identifying a depolarization tick time T D within the preset window using S′′(t) can include: identifying a time point within the preset window when S′′(t) reaches its maximum; and designating the time point when S′′(t) reaches its maximum as the depolarization tick time T D .
- the preset threshold can be computed as a function of an amplitude of S(t), such as 10% of the amplitude of S(t).
- the preset threshold can be computed according to the following steps: identifying a maximum tick time T max and a minimum tick time T min within a preset refractory window surrounding the depolarization tick time T D using S(t); identifying an isopotential tick time T iso occurring before T min using S′′(t); identifying a plateau tick time T plateau occurring after T max using a first derivative S′(t) of S(t); defining the amplitude A of S(t) as S(T plateau ) ⁇ S(T iso ); and defining the preset threshold as S(T plateau ) ⁇ x*A, 0 ⁇ x ⁇ 1.
- the method optionally also includes computing a diastolic interval for the preset window using the depolarization tick time T D within the preset window and an immediately previous repolarization tick time T R-1 .
- a method of mapping cardiac repolarization activity includes: receiving an electrogram signal S(t) at a signal processor; and using the signal processor: identifying a local activation time (“LAT”) tick time within a preset window about a reference time T ref ; defining the LAT tick time as a depolarization tick time T D for the preset window; computing a first derivative S′(t) of the electrogram signal S(t) for the preset window; identifying a repolarization tick time T R using local maxima and local minima of S′(t); and computing an activation recovery interval for the preset window using the depolarization tick time T D and the repolarization tick time T R .
- LAT local activation time
- the electrogram signal S(t) can include a unipolar electrogram signal.
- the step of identifying a local activation time (“LAT”) tick time within a preset window about a reference time T ref can include identifying the LAT tick time within the preset window about the reference time T ref using a bipolar electrogram signal.
- LAT local activation time
- the step of computing a first derivative S′(t) of the electrogram signal S(t) for the preset window includes computing the first derivative S′(t) according to a formula
- the step of identifying a repolarization tick time T R using local maxima and local minima of S′(t) includes: identifying a maximum tick time T max and a minimum tick time T min within a preset refractory window ⁇ RW surrounding the depolarization tick time T D using S(t); computing a LAT amplitude L 1 using S(T max ) and S(T min ); identifying a next cycle tick T next occurring after T max using S(t); determining a search period following the depolarization tick time T D for the repolarization tick time T R ; identifying a time T SLM during the search period at which S′(t) reaches a local maximum; identifying a time T slm during the search period at which S′(t) reaches a local minimum; and defining the repolarization tick time T R using one or more of T SLM and T slm .
- the step of identifying a next cycle tick T next occurring after T max using S(t) can include: defining a local minimum threshold higher than S(T min ); defining a local maximum threshold higher than S(T max ); analyzing S(t) after T max ; and identifying T next as a first to occur of: a time at which S(t) drops below the local minimum threshold; a time at which S(t) exceeds the local maximum threshold; and a time at which S(t) ends.
- the local minimum threshold can be computed as S(T min )+0.3*L 1 ; and the local maximum threshold can be computed as S(T max )+0.3*L 1 .
- the step of determining a search period following the depolarization tick time T D for the repolarization tick time T R includes: computing a search interval L 2 as T next ⁇ T D ⁇ 2( ⁇ RW); defining a start of the search period at T D + ⁇ RW; and defining an end of the search period at T D + ⁇ RW+0.75*L 2 .
- the step of defining the repolarization tick time T R using one or more of T SLM and T slm includes: defining T R as
- T SLM T slm + T SLM + T slm 2 for a biphasic T-wave; defining T R as T SLM for a negative T-wave; and defining T R as T slm for a positive T-wave.
- the method optionally includes computing a diastolic interval for the preset window using the depolarization tick time T D for the preset window and an immediately previous repolarization tick time T R-1 .
- the system includes: a cardiac repolarization detection processor configured: to receive an electrogram S(t); to identify a depolarization tick time T D using the electrogram S(t); to identify a repolarization tick time T R using either a point in time when S(t) passes below a preset threshold or local maxima and local minima of a first derivative S′(t) of the electrogram S(t); and to compute a cardiac repolarization time period using the depolarization tick time T D and the repolarization tick time T R .
- the cardiac repolarization detection processor can also be configured to compute a diastolic interval using the depolarization tick time T D and a an immediately previous repolarization tick time T R-1 .
- the repolarization time period can be an action potential duration and/or an activation recovery interval.
- FIG. 1 is a schematic of an electrophysiology system, such as may be used in an electrophysiology study including mapping cardiac repolarization activity.
- FIG. 2 depicts an exemplary multi-electrode catheter used in an electrophysiology study.
- FIG. 3 is a flowchart of representative steps that can be followed to map cardiac repolarization activity using a monophasic action potential (“MAP”) catheter to determine action potential duration (“APD”).
- MAP monophasic action potential
- APD action potential duration
- FIG. 4 is a representation of cardiac electrical activity with various tick times and other values referenced in the flowchart of FIG. 3 annotated thereon.
- FIG. 5 is a flowchart of representative steps that can be followed to map cardiac repolarization activity using unipolar electrograms (“EGMs”) to determine activation recovery interval (“ARI”).
- EMMs unipolar electrograms
- ARI activation recovery interval
- FIG. 6 is a representation of cardiac electrical activity with various tick times and other values referenced in the flowchart of FIG. 5 annotated thereon.
- the present disclosure provides methods, apparatuses, and systems for the creation of electrophysiology maps (e.g., electrocardiographic maps) that provide information regarding cardiac repolarization activity.
- electrophysiology maps e.g., electrocardiographic maps
- Certain embodiments of the disclosure will be explained with reference to the use of bipolar electrograms, such as measured by a MAP catheter, to determine APD.
- Other embodiments of the disclosure will be explained with reference to the use of unipolar electrograms to determine ARI. It should be understood, however, that the teachings herein can be applied to good advantage in other contexts, such as mapping local electrogram QRS activity width, which can be representative of the pattern of wave conduction.
- FIG. 1 shows a schematic diagram of an electrophysiology system 8 for conducting cardiac electrophysiology studies by navigating a cardiac catheter and measuring electrical activity occurring in a heart 10 of a patient 11 and three-dimensionally mapping the electrical activity and/or information related to or representative of the electrical activity so measured.
- System 8 can be used, for example, to create an anatomical model of the patient's heart 10 using one or more electrodes.
- System 8 can also be used to measure electrophysiology data, including, but not limited to, electrical activation data (e.g., local activation time (“LAT”)), at a plurality of points along a cardiac surface and store the measured data in association with location information for each measurement point at which the electrophysiology data was measured, for example to create an electrophysiology map of the patient's heart 10 (or a portion thereof).
- electrical activation data e.g., local activation time (“LAT”)
- LAT local activation time
- system 8 can determine the location, and in some aspects the orientation, of objects, typically within a three-dimensional space, and express those locations as position information determined relative to at least one reference.
- the patient 11 is depicted schematically as an oval.
- three sets of surface electrodes e.g., patch electrodes
- a surface of the patient 11 defining three generally orthogonal axes, referred to herein as an x-axis, a y-axis, and a z-axis.
- the electrodes could be positioned in other arrangements, for example multiple electrodes on a particular body surface.
- the electrodes do not need to be on the body surface, but could be positioned internally to the body or on an external frame.
- the x-axis surface electrodes 12 , 14 are applied to the patient along a first axis, such as on the lateral sides of the thorax region of the patient (e.g., applied to the patient's skin underneath each arm) and may be referred to as the Left and Right electrodes.
- the y-axis electrodes 18 , 19 are applied to the patient along a second axis generally orthogonal to the x-axis, such as along the inner thigh and neck regions of the patient, and may be referred to as the Left Leg and Neck electrodes.
- the z-axis electrodes 16 , 22 are applied along a third axis generally orthogonal to both the x-axis and the y-axis, such as along the sternum and spine of the patient in the thorax region, and may be referred to as the Chest and Back electrodes.
- the heart 10 lies between these pairs of surface electrodes 12 / 14 , 18 / 19 , and 16 / 22 .
- An additional surface reference electrode (e.g., a “belly patch”) 21 provides a reference and/or ground electrode for the system 8 .
- the belly patch electrode 21 may be an alternative to a fixed intra-cardiac electrode 31 , described in further detail below.
- the patient 11 may have most or all of the conventional electrocardiogram (“ECG” or “EKG”) system leads in place.
- ECG electrocardiogram
- a standard set of 12 ECG leads may be utilized for sensing electrocardiograms on the patient's heart 10 .
- This ECG information is available to the system 8 (e.g., it can be provided as input to computer system 20 ). Insofar as ECG leads are well understood, and for the sake of clarity in the figures, only one lead 6 and its connection to computer system 20 is illustrated in FIG. 1 .
- a representative catheter 13 having at least one electrode 17 is also depicted in schematic fashion in FIG. 1 .
- This representative catheter electrode 17 can be referred to as a “measurement electrode” or a “roving electrode.”
- a measurement electrode or a “roving electrode.”
- multiple electrodes on catheter 13 or on multiple such catheters, will be used.
- system 8 may utilize sixty-four electrodes on twelve catheters disposed within the heart and/or vasculature of the patient.
- system 8 may utilize a single catheter that includes multiple (e.g., eight) splines, each of which in turn includes multiple (e.g., eight) electrodes.
- multiple splines each of which in turn includes multiple (e.g., eight) electrodes.
- these embodiments are merely exemplary, and any number of electrodes and catheters may be used.
- a high density mapping catheter such as the EnSiteTM ArrayTM non-contact mapping catheter of St. Jude Medical, Inc., can be utilized.
- catheter 13 (or multiple such catheters) are typically introduced into the heart and/or vasculature of the patient via one or more introducers and using familiar procedures.
- a segment of an exemplary multi-electrode catheter 13 is shown in FIG. 2 .
- catheter 13 extends into the left ventricle 50 of the patient's heart 10 through a transseptal sheath 35 .
- transseptal sheath 35 The use of a transseptal approach to the left ventricle is well known and will be familiar to those of ordinary skill in the art, and need not be further described herein.
- catheter 13 can also be introduced into the heart 10 in any other suitable manner.
- Catheter 13 includes electrode 17 on its distal tip, as well as a plurality of additional measurement electrodes 52 , 54 , 56 spaced along its length in the illustrated embodiment. Typically, the spacing between adjacent electrodes will be known, though it should be understood that the electrodes may not be evenly spaced along catheter 13 or of equal size to each other. Since each of these electrodes 17 , 52 , 54 , 56 lies within the patient, location data may be collected simultaneously for each of the electrodes by system 8 .
- each of electrodes 17 , 52 , 54 , and 56 can be used to gather electrophysiological data from the cardiac surface.
- catheter 13 is capable of sensing MAPs (e.g., from bipolar electrograms).
- MAPs e.g., from bipolar electrograms.
- electrophysiology data points including, for example, both contact and non-contact electrophysiological mapping and the collection of both unipolar and bipolar electrograms
- various techniques familiar in the art can be used to generate a graphical representation from the plurality of electrophysiology data points. Insofar as the ordinarily skilled artisan will appreciate how to create electrophysiology maps from electrophysiology data points, the aspects thereof will only be described herein to the extent necessary to understand the maps disclosed herein.
- a fixed reference electrode 31 (e.g., attached to a wall of the heart 10 ) is shown on a second catheter 29 .
- this electrode 31 may be stationary (e.g., attached to or near the wall of the heart) or disposed in a fixed spatial relationship with the roving electrodes (e.g., electrodes 17 , 52 , 54 , 56 ), and thus may be referred to as a “navigational reference” or “local reference,”
- the fixed reference electrode 31 may be used in addition or alternatively to the surface reference electrode 21 described above.
- a coronary sinus electrode or other fixed electrode in the heart 10 can be used as a reference for measuring voltages and displacements; that is, as described below, fixed reference electrode 31 may define the origin of a coordinate system.
- Each surface electrode is coupled to a multiplex switch 24 , and the pairs of surface electrodes are selected by software running on a computer 20 , which couples the surface electrodes to a signal generator 25 .
- switch 24 may be eliminated and multiple (e.g., three) instances of signal generator 25 may be provided, one for each measurement axis (that is, each surface electrode pairing).
- the computer 20 may comprise a conventional general-purpose computer, a special-purpose computer, a distributed computer, or any other type of computer.
- the computer 20 may comprise one or more processors 28 , such as a single central processing unit (CPU), or a plurality of processing units, commonly referred to as a parallel processing environment, which may execute instructions to practice the various aspects disclosed herein.
- processors 28 such as a single central processing unit (CPU), or a plurality of processing units, commonly referred to as a parallel processing environment, which may execute instructions to practice the various aspects disclosed herein.
- three nominally orthogonal electric fields are generated by a series of driven and sensed electric dipoles (e.g., surface electrode pairs 12 / 14 , 18 / 19 , and 16 / 22 ) in order to realize catheter navigation in a biological conductor.
- these orthogonal fields can be decomposed and any pairs of surface electrodes can be driven as dipoles to provide effective electrode triangulation.
- the electrodes 12 , 14 , 18 , 19 , 16 , and 22 (or any other number of electrodes) could be positioned in any other effective arrangement for driving a current to or sensing a current from an electrode in the heart.
- multiple electrodes could be placed on the back, sides, and/or belly of patient 11 .
- the potentials measured across the roving electrodes resulting from a predetermined set of drive (source-sink) configurations may be combined algebraically to yield the same effective potential as would be obtained by simply driving a uniform current along the orthogonal axes.
- any two of the surface electrodes 12 , 14 , 16 , 18 , 19 , 22 may be selected as a dipole source and drain with respect to a ground reference, such as belly patch 21 , while the unexcited electrodes measure voltage with respect to the ground reference.
- the roving electrodes 17 , 52 , 54 , 56 placed in the heart 10 are exposed to the field from a current pulse and are measured with respect to ground, such as belly patch 21 .
- the catheters within the heart 10 may contain more or fewer electrodes than the four shown, and each electrode potential may be measured.
- At least one electrode may be fixed to the interior surface of the heart to form a fixed reference electrode 31 , which is also measured with respect to ground, such as belly patch 21 , and which may be defined as the origin of the coordinate system relative to which localization system 8 measures positions. Data sets from each of the surface electrodes, the internal electrodes, and the virtual electrodes may all be used to determine the location of the roving electrodes 17 , 52 , 54 , 56 within heart 10 .
- the measured voltages may be used by system 8 to determine the location in three-dimensional space of the electrodes inside the heart, such as roving electrodes 17 , 52 , 54 , 56 , relative to a reference location, such as reference electrode 31 . That is, the voltages measured at reference electrode 31 may be used to define the origin of a coordinate system, while the voltages measured at roving electrodes 17 , 52 , 54 , 56 may be used to express the location of roving electrodes 17 , 52 , 54 , 56 relative to the origin.
- the coordinate system is a three-dimensional (x, y, z) Cartesian coordinate system, although other coordinate systems, such as polar, spherical, and cylindrical coordinate systems, are contemplated.
- the data used to determine the location of the electrode(s) within the heart is measured while the surface electrode pairs impress an electric field on the heart.
- the electrode data may also be used to create a respiration compensation value used to improve the raw location data for the electrode locations as described in U.S. Pat. No. 7,263,397, which is hereby incorporated herein by reference in its entirety.
- the electrode data may also be used to compensate for changes in the impedance of the body of the patient as described, for example, in U.S. Pat. No. 7,885,707, which is also incorporated herein by reference in its entirety.
- the system 8 first selects a set of surface electrodes and then drives them with current pulses. While the current pulses are being delivered, electrical activity, such as the voltages measured with at least one of the remaining surface electrodes and in vivo electrodes, is measured and stored. Compensation for artifacts, such as respiration and/or impedance shifting, may be performed as indicated above.
- system 8 is the EnSiteTM VelocityTM cardiac mapping and visualization system of St. Jude Medical, Inc., which generates electrical fields as described above, or another such system that relies upon electrical fields.
- Other systems may be used in connection with the present teachings, including for example, the CARTO navigation and location system of Biosense Webster, Inc., the AURORA® system of Northern Digital Inc., or Sterotaxis' NIOBE® Magnetic Navigation System, all of Which utilize magnetic fields rather than electrical fields.
- the localization and mapping systems described in the following patents can also be used with the present invention: U.S. Pat. Nos. 6,990,370; 6,978,168; 6,947,785; 6,939,309; 6,728,562; 6,640,119; 5,983,126; and 5,697,377.
- flowchart 300 may represent several exemplary steps that can be carried out by the computer 20 of FIG. 1 (e.g., by one or more processors 28 ) to generate a map of cardiac repolarization activity as described herein.
- processors 28 e.g., by one or more processors 28
- the representative steps described below can be either hardware- or software-implemented.
- the term “signal processor” is used herein to describe both hardware- and software-based implementations of the teachings herein.
- flowchart 300 in FIG. 3 illustrates a series of representative steps that may be carried out to map cardiac repolarization activity using APD.
- an electrogram signal denoted S(t) (and illustrated as lower trace 402 in FIG. 4 )
- a signal processor e.g., by one or more processors 28 within computer 20 .
- the electrogram signal 402 be a bipolar electrogram signal, such as from a MAP catheter (see, e.g., U.S. Pat. No. 5,398,683, which is hereby incorporated by reference as though fully set forth herein).
- the first and second derivatives of S(t), denoted S′(t) and S′′(t), respectively, are computed.
- the first derivative Sit) is calculated using the formula
- T ref 406 is set to correspond to QRS activity detected using a user-defined reference cardiac signal, such as the signal from an EKG lead or the signal from an in vivo reference electrode.
- the upper trace 400 in FIG. 4 is a representative reference cardiac signal, such as may be used to determine T ref 406 .
- RAI 404 can be user-defined. According to aspects of the disclosure, RAI 404 is between about 100 ms and about 300 ms wide.
- a depolarization tick time 408 is identified within the RAI using S′′(t).
- T D 408 can be designated as the time point within RAI 404 where S′′(t) reaches its maximum.
- T R a repolarization tick time 410
- the current cardiac trigger that is, the current beat.
- the repolarization tick time 410 for the current cardiac trigger need not occur within RAI 404 . Rather, according to some aspects of the disclosure, the repolarization tick time 410 can be identified within a user-defined minimum cycle length threshold (“CL min ”) of T plateau 412 (described further below). In any event, however, T R 410 is typically identified as a point in time where S(t) 402 passes below a preset threshold, as further described below.
- CL min minimum cycle length threshold
- a maximum tick time 414 (denoted T max ) and a minimum tick time 416 (denoted T min ) are identified within a preset refractory window 418 (denoted ⁇ RW) of T D 408 .
- T max 414 is the time point, within the refractory window 418 , at which S(t) 402 achieves its local maximum; conversely, T min 416 is the time point, within the refractory window 418 , at which S(t) 402 achieves its local minimum.
- the preset refractory window 418 is 12 ms to either side of the depolarization tick time 408 (i.e., T D ⁇ 12 ms).
- T D depolarization tick time 408
- teachings herein are not limited, however, to this specific refractory window, and other refractory windows, including refractory windows that are not centered on the depolarization tick time 408 , can be employed without departing from the scope of the teachings herein.
- An isopotential tick time 420 is detected in block 312 .
- the isopotential tick time 420 is detected as a time point before the minimum tick time 416 at which the absolute value of S′′(t) drops below a preset threshold, such as 0.01. If no such point exists, then flowchart 300 can terminate for the current beat.
- a plateau tick time 412 (T plateau ) is detected.
- the plateau tick time 412 is detected as a time point after the maximum tick time 414 at which the absolute value of S′(t) goes below a preset threshold, such as 0.5. If no such point exists, then flowchart 300 can terminate for the current beat.
- a MAP amplitude 420 (represented by the letter A in this disclosure and shown as “MAP amplitude” in FIG. 4 ) is computed as the difference between S(T plateau ) 422 and S(T iso ) 424 .
- the MAP amplitude 420 is used to determine the preset threshold that will be used to determine the repolarization tick time 410 (block 318 ), for example by setting the preset threshold as a fraction of the MAP amplitude 420 .
- the preset threshold is set at 10% of the MAP amplitude 420 (that is, the repolarization tick time is detected when S(t) drops below S(T plateau ) ⁇ 0.9(A)). If no such point exists within the minimum cycle length threshold of the plateau tick time 412 (that is, if the amplitude of S(t) 402 does not diminish sufficiently within a preset time of when it plateaus), then flowchart 300 can terminate for the current beat.
- APD can be computed as the difference between the depolarization and repolarization tick times 408 , 410 respectively (block 320 ).
- a diastolic interval can be computed as the difference between the depolarization tick time and an immediately previous repolarization tick time (i.e., T R-1 ).
- flowchart 500 may represent several exemplary steps that can be carried out by the computer 20 of FIG. 1 (e.g. by one or more processors 28 ) to generate a map of cardiac repolarization activity as described herein.
- processors 28 e.g. by one or more processors 28
- FIG. 5 Another basic methodology of mapping cardiac repolarization activity will be explained with reference to the flowchart 500 of representative steps presented as FIG. 5 . That is, flowchart 500 may represent several exemplary steps that can be carried out by the computer 20 of FIG. 1 (e.g. by one or more processors 28 ) to generate a map of cardiac repolarization activity as described herein.
- the representative steps described below can be either hardware- or software-implemented, and the term “signal processor” will again be used to describe both hardware- and software-based implementations of the teachings herein.
- flowchart 500 in FIG. 5 illustrates a series of representative steps that may be carried out to map cardiac repolarization activity using ARI.
- an electrogram signal again denoted S(t) and illustrated as lower trace 602 in FIG. 6 , is received at a signal processor (e.g., by one or more processors 28 within computer 20 ).
- a signal processor e.g., by one or more processors 28 within computer 20 .
- the electrogram signal be a unipolar electrogram signal.
- a local activation time (“LAT”) tick time 604 is identified within a preset window, again referred to as a “Roving Activation Interval” (“RAI”) 606 , about a reference time point T ref 608 .
- the upper trace 600 in FIG. 6 is a representative reference cardiac signal, such as may be used to identify T ref 608 .
- the LAT tick time 604 is identified using a bipolar electrogram signal 602 . Once identified, the LAT tick time 604 is defined as the depolarization tick time T D for the RAI 606 .
- the first derivative S′(t) of SW 602 is computed.
- S′(t) can be computed as
- S ′ ⁇ ( t ) S ⁇ ( t + ⁇ ⁇ ⁇ t ) - S ⁇ ( t ) ⁇ ⁇ ⁇ t , where ⁇ t is a preset time interval, such as 3 ms.
- the next several steps illustrated in FIG. 5 are directed at identifying a repolarization tick time T R for the current cardiac trigger (that is, the current beat).
- T R a repolarization tick time for the current cardiac trigger
- the repolarization tick time for the current cardiac trigger need not occur within the RAI 606 .
- a search period following the depolarization tick time 604 is defined, and the repolarization tick time is detected within this search period.
- T R is typically identified using local maxima and local minima of S′(t), as described in further detail below.
- a maximum tick time 610 (denoted T max ) and a minimum tick time 612 (denoted T min ) are identified within a preset refractory window 614 (denoted ⁇ RW) of T D 604 .
- T max 610 is the time point, within the refractory window 614 , at Which S(t) 602 achieves its local maximum; conversely.
- T min 612 is the time point, within the refractory window 614 , at which S(t) 602 achieves its local minimum.
- the preset refractory window 614 is 50 ms to either side of the depolarization tick time 604 (i.e., T D ⁇ 50 ms).
- T D depolarization tick time
- the teachings herein are not limited, however, to this specific refractory window, and other refractory windows, including refractory windows that are not centered on the depolarization tick time 604 , can be employed without departing from the spirit of the teachings herein.
- a LAT amplitude 616 (L 1 ) is computed as the difference between S(T max ) 618 and S(T min ) 620 in block 512 .
- T next 622 is identified in block 514 . It is contemplated that T next 622 can be identified with reference to a local minimum threshold, which is higher than S(T min ) 620 , and a local maximum threshold, which is higher than S(T max ) 618 .
- the LAT amplitude 616 can be used to compute the local minimum and local maximum thresholds.
- the local minimum threshold can Be computed as S(T min )+0.3*L 1
- the local maximum threshold can be computed as S(T max )+0.3*L 1 .
- T next 622 is designated as the time point at which the first of the following occurs after T max 610 : (1) S(t) 602 drops below the local minimum threshold; (2) S(t) 602 exceeds the local maximum threshold; or (3) S(t) 602 ends.
- a search period, following the depolarization tick time 604 , for the repolarization tick time is determined in block 516 .
- the search period begins at T D + ⁇ RW and ends at TD+ ⁇ RW+0.75(L 2 ), where L 2 ( 624 ) is a search interval defined as T next ⁇ T D ⁇ 2( ⁇ RW).
- T SLM is the time point during the search period at which S′(t) reaches a local maximum
- T slm is the time point during the search period at which S′(t) reaches a local minimum.
- the repolarization tick time is defined using one or more of T SLM and T slm .
- the definition of the repolarization tick time depends on whether the T-wave is biphasic, negative, or positive.
- the repolarization tick time can be defined as T SLM 626 for a negative T-wave and as T slm 628 for a positive T-wave.
- the repolarization tick time 630 can be defined using the formula
- ART 632 can be computed (block 522 ) as the difference between the depolarization and repolarization tick times.
- a diastolic interval can be computed (block 524 ) as the difference between the depolarization tick time and an immediately previous repolarization tick time T R-1 ).
- All directional references e.g., upper, lower, upward, downward, left, right, leftward, rightward, top, bottom, above, below, vertical, horizontal, clockwise, and counterclockwise
- Joinder references e.g., attached, coupled, connected, and the like
- Joinder references are to be construed broadly and may include intermediate members between a connection of elements and relative movement between elements. As such, joinder references do not necessarily infer that two elements are directly connected and in fixed relation to each other.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Cardiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Description
wherein Δt comprises a preset time interval, such as about 3 ms.
wherein Δt comprises a preset time interval, such as about 3 ms.
for a biphasic T-wave; defining TR as TSLM for a negative T-wave; and defining TR as Tslm for a positive T-wave.
where Δt is a preset time interval, such as 3 ms. Similarly, the second derivative S″(t) can be computed according to the formula
where Δt is a preset time interval, such as 3 ms.
Claims (9)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/766,283 US11071491B2 (en) | 2015-10-07 | 2016-09-29 | Methods and systems for mapping cardiac repolarization |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201562238323P | 2015-10-07 | 2015-10-07 | |
US15/766,283 US11071491B2 (en) | 2015-10-07 | 2016-09-29 | Methods and systems for mapping cardiac repolarization |
PCT/US2016/054280 WO2017062247A1 (en) | 2015-10-07 | 2016-09-29 | Methods and systems for mapping cardiac repolarization |
Publications (2)
Publication Number | Publication Date |
---|---|
US20180289276A1 US20180289276A1 (en) | 2018-10-11 |
US11071491B2 true US11071491B2 (en) | 2021-07-27 |
Family
ID=57113796
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/766,283 Active 2037-11-23 US11071491B2 (en) | 2015-10-07 | 2016-09-29 | Methods and systems for mapping cardiac repolarization |
Country Status (5)
Country | Link |
---|---|
US (1) | US11071491B2 (en) |
EP (1) | EP3344136B1 (en) |
JP (2) | JP2018534035A (en) |
CN (1) | CN108289631B (en) |
WO (1) | WO2017062247A1 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019152986A1 (en) * | 2018-02-05 | 2019-08-08 | Mayo Foundation For Medical Education And Research | Systems and methods for mapping and modulating repolarization |
CN113795196A (en) * | 2019-05-24 | 2021-12-14 | 圣犹达医疗用品心脏病学部门有限公司 | Systems and methods for cardiac mapping |
US11439339B1 (en) * | 2020-10-16 | 2022-09-13 | Neucures Inc. | Method and system for measuring cardiac tissue health based on DV/DTmin of a depolarization wave within a cardiac electrogram |
Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5423878A (en) | 1984-03-06 | 1995-06-13 | Ep Technologies, Inc. | Catheter and associated system for pacing the heart |
US5697377A (en) | 1995-11-22 | 1997-12-16 | Medtronic, Inc. | Catheter mapping system and method |
WO2001013992A1 (en) | 1999-08-20 | 2001-03-01 | Impulse Dynamics N.V. | Apparatus and method for controlling the delivery of contractility modulating non-excitatory signals to the heart |
US6640119B1 (en) | 1992-09-23 | 2003-10-28 | Endocardial Solutions, Inc. | Method for orienting an electrode array |
WO2004096354A1 (en) | 2003-04-29 | 2004-11-11 | Medtronic, Inc. | Method and apparatus for detecting myocardial electrical recovery and controlling extra-systolic stimulation |
US6939309B1 (en) | 1993-09-23 | 2005-09-06 | Endocardial Solutions, Inc. | Electrophysiology mapping system |
US20050261599A1 (en) * | 2004-05-21 | 2005-11-24 | Beth Israel Deaconess Medical Center | Differentiating ischemic from non-ischemic T-wave inversion |
US6978168B2 (en) | 1992-09-23 | 2005-12-20 | Endocardial Solutions, Inc. | Software for mapping potential distribution of a heart chamber |
US7263397B2 (en) | 1998-06-30 | 2007-08-28 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Method and apparatus for catheter navigation and location and mapping in the heart |
WO2008008361A2 (en) | 2006-07-11 | 2008-01-17 | Covance, Inc. | Rms electrocardiography system and method |
US20080125665A1 (en) | 2006-11-27 | 2008-05-29 | Nigam Indra B | Implantable Medical Device |
US20090299424A1 (en) | 2008-05-13 | 2009-12-03 | Narayan Sanjiv M | Methods and systems for treating heart instability |
US20100094274A1 (en) * | 2008-10-09 | 2010-04-15 | Sanjiv Narayan | Methods, system and appartus for the detection, diagnosis and treatment of biological rhythm disorders |
WO2010054409A1 (en) | 2008-11-10 | 2010-05-14 | Cardioinsight Technologies, Inc. | Visualization of electrophysiology data |
US7885707B2 (en) | 2005-09-15 | 2011-02-08 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Method of scaling navigation signals to account for impedance drift in tissue |
US8265756B1 (en) * | 2006-12-07 | 2012-09-11 | Pacesetter, Inc. | System and method for adjusting automatic sensitivity control parameters based on intracardiac electrogram signals |
US20140187989A1 (en) | 2013-01-02 | 2014-07-03 | Boston Scientific Scimed, Inc. | Estimating restitution curves in an anatomical mapping system |
US20150342488A1 (en) * | 2014-05-09 | 2015-12-03 | Boston Scientific Scimed, Inc. | Medical devices for mapping cardiac tissue |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE112007002930T5 (en) * | 2006-11-30 | 2010-04-29 | Aalborg Universitet | System and method for analyzing the complex curvature of ECG waveforms |
US10398326B2 (en) * | 2013-03-15 | 2019-09-03 | The Regents Of The University Of California | System and method of identifying sources associated with biological rhythm disorders |
US10039502B2 (en) * | 2011-04-12 | 2018-08-07 | Medtronic Ablation Frontiers Llc | Electrophysiological signal processing and utilization |
CN103207168A (en) * | 2012-12-26 | 2013-07-17 | 西安交通大学 | Optical reconstruction method for potential functional information of myocardial cell membranes |
US9132274B2 (en) * | 2013-07-26 | 2015-09-15 | Medtronic, Inc. | Determining onsets and offsets of cardiac depolarization and repolarization waves |
-
2016
- 2016-09-29 EP EP16778637.5A patent/EP3344136B1/en active Active
- 2016-09-29 WO PCT/US2016/054280 patent/WO2017062247A1/en active Application Filing
- 2016-09-29 JP JP2018517839A patent/JP2018534035A/en not_active Withdrawn
- 2016-09-29 US US15/766,283 patent/US11071491B2/en active Active
- 2016-09-29 CN CN201680066767.1A patent/CN108289631B/en active Active
-
2020
- 2020-03-09 JP JP2020039891A patent/JP6914383B2/en active Active
Patent Citations (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5423878A (en) | 1984-03-06 | 1995-06-13 | Ep Technologies, Inc. | Catheter and associated system for pacing the heart |
US6978168B2 (en) | 1992-09-23 | 2005-12-20 | Endocardial Solutions, Inc. | Software for mapping potential distribution of a heart chamber |
US6728562B1 (en) | 1992-09-23 | 2004-04-27 | Endocardial Solutions, Inc. | Method for creating a virtual electrogram |
US6990370B1 (en) | 1992-09-23 | 2006-01-24 | Endocardial Solutions, Inc. | Method for mapping heart electrophysiology |
US6640119B1 (en) | 1992-09-23 | 2003-10-28 | Endocardial Solutions, Inc. | Method for orienting an electrode array |
US6939309B1 (en) | 1993-09-23 | 2005-09-06 | Endocardial Solutions, Inc. | Electrophysiology mapping system |
US6947785B1 (en) | 1993-09-23 | 2005-09-20 | Endocardial Solutions, Inc. | Interface system for endocardial mapping catheter |
US5983126A (en) | 1995-11-22 | 1999-11-09 | Medtronic, Inc. | Catheter location system and method |
US5697377A (en) | 1995-11-22 | 1997-12-16 | Medtronic, Inc. | Catheter mapping system and method |
US7263397B2 (en) | 1998-06-30 | 2007-08-28 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Method and apparatus for catheter navigation and location and mapping in the heart |
WO2001013992A1 (en) | 1999-08-20 | 2001-03-01 | Impulse Dynamics N.V. | Apparatus and method for controlling the delivery of contractility modulating non-excitatory signals to the heart |
WO2004096354A1 (en) | 2003-04-29 | 2004-11-11 | Medtronic, Inc. | Method and apparatus for detecting myocardial electrical recovery and controlling extra-systolic stimulation |
US20050261599A1 (en) * | 2004-05-21 | 2005-11-24 | Beth Israel Deaconess Medical Center | Differentiating ischemic from non-ischemic T-wave inversion |
US7885707B2 (en) | 2005-09-15 | 2011-02-08 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Method of scaling navigation signals to account for impedance drift in tissue |
WO2008008361A2 (en) | 2006-07-11 | 2008-01-17 | Covance, Inc. | Rms electrocardiography system and method |
US20080125665A1 (en) | 2006-11-27 | 2008-05-29 | Nigam Indra B | Implantable Medical Device |
US8265756B1 (en) * | 2006-12-07 | 2012-09-11 | Pacesetter, Inc. | System and method for adjusting automatic sensitivity control parameters based on intracardiac electrogram signals |
US20090299424A1 (en) | 2008-05-13 | 2009-12-03 | Narayan Sanjiv M | Methods and systems for treating heart instability |
US20100094274A1 (en) * | 2008-10-09 | 2010-04-15 | Sanjiv Narayan | Methods, system and appartus for the detection, diagnosis and treatment of biological rhythm disorders |
WO2010054409A1 (en) | 2008-11-10 | 2010-05-14 | Cardioinsight Technologies, Inc. | Visualization of electrophysiology data |
US20140187989A1 (en) | 2013-01-02 | 2014-07-03 | Boston Scientific Scimed, Inc. | Estimating restitution curves in an anatomical mapping system |
US20150342488A1 (en) * | 2014-05-09 | 2015-12-03 | Boston Scientific Scimed, Inc. | Medical devices for mapping cardiac tissue |
Non-Patent Citations (2)
Title |
---|
International Search Report and Written Opinion for PCT/US2016/054280, dated Jan. 3, 2017. |
P. Colli Franzone et al. "Monophasic action potentials generated by bidomain modeling as a tool for detecting cardiac repolarization times," American Journal of Physiology: Heart and Circulatory Physiology, vol. 293, No. 5., dated Jan. 1, 2007. |
Also Published As
Publication number | Publication date |
---|---|
JP2018534035A (en) | 2018-11-22 |
CN108289631A (en) | 2018-07-17 |
JP2020103940A (en) | 2020-07-09 |
JP6914383B2 (en) | 2021-08-04 |
EP3344136B1 (en) | 2020-07-01 |
WO2017062247A1 (en) | 2017-04-13 |
EP3344136A1 (en) | 2018-07-11 |
CN108289631B (en) | 2021-03-02 |
US20180289276A1 (en) | 2018-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2598202B1 (en) | System for identifying cardiac arrhythmia driver sites | |
US9901271B2 (en) | System and method for analyzing biological signals and generating electrophysiology maps | |
US20160317094A1 (en) | Coaxial electrode catheters for extracting electrophysiologic parameters | |
US10799188B2 (en) | System and method for real-time electrophysiological mapping | |
US11931158B2 (en) | Methods and systems for mapping cardiac activity | |
US9220435B2 (en) | System and method for generating electrophysiology maps | |
EP3344135B1 (en) | Methods and systems for mapping cardiac restitution | |
JP6914383B2 (en) | Methods and systems for mapping cardiac repolarization | |
US11103177B2 (en) | System and method for mapping cardiac activity | |
US11071486B2 (en) | System and method for generating activation timing maps | |
EP3609400A1 (en) | System and method for electrophysiological mapping |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
AS | Assignment |
Owner name: ST JUDE MEDICAL, CARDIOLOGY DIVISION, INC., MINNESOTA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:RELAN, JATIN S.;AFONSO, VALTINO X.;SIGNING DATES FROM 20151030 TO 20151201;REEL/FRAME:047970/0691 Owner name: ST JUDE MEDICAL, CARDIOLOGY DIVISION, INC., MINNES Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:RELAN, JATIN S.;AFONSO, VALTINO X.;SIGNING DATES FROM 20151030 TO 20151201;REEL/FRAME:047970/0691 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: FINAL REJECTION MAILED |
|
STCV | Information on status: appeal procedure |
Free format text: NOTICE OF APPEAL FILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: AWAITING TC RESP., ISSUE FEE NOT PAID |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT RECEIVED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |