Source: https://patents.google.com/patent/US9278219B2/en
Timestamp: 2019-04-19 03:42:14
Document Index: 22030782

Matched Legal Cases: ['art 112', 'art 112', 'art 112', 'art 112', 'art 112', 'art 112', 'art 112', 'art 112', 'art 112', 'art 112', 'art 112', 'art 112', 'art 112', 'art 300', 'art 500']

US9278219B2 - Closed loop optimization of control parameters during cardiac pacing - Google Patents
Closed loop optimization of control parameters during cardiac pacing Download PDF
US9278219B2
US9278219B2 US13/833,882 US201313833882A US9278219B2 US 9278219 B2 US9278219 B2 US 9278219B2 US 201313833882 A US201313833882 A US 201313833882A US 9278219 B2 US9278219 B2 US 9278219B2
US13/833,882
US20140277233A1 (en
2013-03-15 Priority to US13/833,882 priority Critical patent/US9278219B2/en
2013-03-15 Assigned to MEDTRONIC, INC. reassignment MEDTRONIC, INC. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: GHOSH, Subham
2014-09-18 Publication of US20140277233A1 publication Critical patent/US20140277233A1/en
2016-03-08 Publication of US9278219B2 publication Critical patent/US9278219B2/en
A system and method control a pacing parameter in a closed-loop manner by determining a value of an EGM-based index corresponding an optimal electrical activation condition of a patient's heart and adjusting a pacing therapy to maintain the EGM-based index value. The closed loop control method performed by the system may establish a relationship between an EGM-based index and multiple settings of a pacing control parameter. Values of the EGM-based index are stored with corresponding setting shifts relative to a previously established optimal setting. A processor of an implantable medical device monitors the EGM-based index during cardiac pacing. Responsive to detecting an EGM-based index value corresponding to a non-optimal setting of the control parameter, the processor determines an adjustment of the control parameter from the stored index values and corresponding setting shifts.
The disclosure relates generally to optimizing control parameters during cardiac pacing therapies and more particularly to optimizing pacing therapy control parameters based on relationships between an index determined from intracardiac electrogram (EGM) signals and ventricular activation.
Cardiac resynchronization therapy (CRT) is a treatment for heart failure patients in which one or more heart chambers are electrically stimulated (paced) to restore or improve heart chamber synchrony. Improved heart chamber synchrony is expected to improve hemodynamic performance of the heart, such as assessed by ventricular pressure and the rate of change in ventricular pressure or other hemodynamic parameters, thereby alleviating symptoms of heart failure. Achieving a positive clinical benefit from CRT is dependent on several therapy control parameters, such as the atrio-ventricular (AV) delay and the inter-ventricular (VV) delay. The AV delay controls the timing of ventricular pacing pulses relative to an atrial depolarization, intrinsic or paced. The VV delay controls the timing of a pacing pulse in one ventricle relative to a paced or intrinsic sensed event in the other ventricle.
Numerous methods for selecting optimal AV and VV delays for use in controlling CRT pacing pulses have been proposed. For example, clinicians may select an optimal AV or VV delay using Doppler echocardiography or other imaging modalities to optimize a hemodynamic variable. Such clinical techniques are time-consuming, require an expert technician to perform, and are performed at a discreet time, such as at device implantation or during a clinical visit, without ongoing adjustment to the CRT therapy parameters as the need may arise. Other methods may be based on a hemodynamic sensor signal or a sensor of mechanical heart function. Such methods may require additional sensors and therefore add cost to the overall system. A need remains, therefore, for a device and method that enables closed loop optimization of CRT therapy parameters for restoring ventricular synchrony and alleviating the symptoms of heart failure associated with ventricular dyssynchrony.
In general, the disclosure is directed towards techniques for controlling a pacing control parameter during a pacing therapy. In particular, a control parameter such as a timing parameter, e.g., an AV delay or a VV delay, is controlled by an implantable medical device processor using a relationship established between an index derived from an EGM signal and pacing control parameter settings corresponding to an optimal electrical activation condition and corresponding to non-optimal electrical activation conditions. An implantable medical device processor may be configured to determine a value of the EGM-based index during an unknown electrical activation condition of a patient's heart and adjust a pacing therapy to maintain the EGM-based index at a value previously established as corresponding to optimal electrical activation synchrony.
In some examples, an optimal control parameter setting is established by a processor determining a ventricular dyssynchrony index from surface ECG electrode signals. Activation times are determined from surface ECG signals that are acquired from multiple surface electrodes. A ventricular dyssynchrony index is determined from the activation times. An optimal setting of the control parameter is determined based on the ECG-based ventricular dyssynchrony index, establishing an optimal electrical activation condition for which the EGM-based index corresponding to the optimal activation condition is determined.
The relationship between an EGM-based index and multiple settings of the control parameter is established by delivering cardiac pacing according to each of a plurality of control parameter settings, including a previously established optimal setting and at least one non-optimal setting different than the optimal setting. The EGM-based index is determined from a EGM signal acquired during each of the settings. Values of the determined index are stored with corresponding parameter setting differences or shifts relative to the optimal setting for each of the multiple settings. The EGM-based index is monitored during cardiac pacing. Responsive to detecting an EGM-based index value corresponding to a non-optimal setting of the control parameter, an adjustment of the control parameter is determined from the stored relationship between the EGM-based index and corresponding setting shifts.
In one example, a method includes establishing a relationship between an EGM-based index and a plurality of settings of a pacing control parameter. Establishing this relationship includes delivering cardiac pacing according to each of a plurality of settings comprising a previously established optimal setting and at least one non-optimal setting different than the optimal setting by a setting shift, determining the EGM-based index from an EGM signal acquired for each of the plurality of settings, and storing values of the determined index with corresponding setting shifts relative to the optimal setting for each of the plurality of settings. The method further includes monitoring the EGM-based index during cardiac pacing, and, responsive to detecting an EGM-based index value corresponding to a non-optimal setting of the control parameter, determining an adjustment of the control parameter from the stored index values and corresponding setting shifts. The control parameter is adjusted by the determined adjustment.
In another example, a system includes a plurality of cardiac electrodes, a cardiac signal generator coupled to the plurality of cardiac electrodes, a sensing module coupled to the plurality of electrodes, and a processor. The cardiac signal generator delivers cardiac pacing according to each of a plurality of settings of a pacing control parameter, the plurality of settings comprising a previously established optimal setting and at least one non-optimal setting different than the optimal setting by a setting shift. The sensing module senses an EGM signal during the cardiac pacing. The processor is configured to establish a relationship between an EGM-based index and the plurality of settings by determining the EGM-based index from the EGM signal acquired for each of the plurality of settings and storing values of the determined index with corresponding setting shifts of each of the plurality of settings relative to the optimal setting. The processor is further configured to monitor the EGM-based index during cardiac pacing; and, responsive to detecting an EGM-based index value corresponding to a non-optimal setting of the control parameter, determine an adjustment of the control parameter from the stored index values and corresponding setting shifts. The processor adjusts the control parameter by the determined adjustment
In another example a non-transitory, computer-readable storage medium includes instructions that, when executed, cause a processor to establish a relationship between an EGM-based index and a plurality of settings of a pacing control parameter, monitor the EGM-based index during cardiac pacing, responsive to detecting an EGM-based index value corresponding to a non-optimal setting of the control parameter, determine an adjustment of the control parameter from the stored index values and corresponding setting shifts, and adjust the control parameter by the determined adjustment.
FIG. 1 is a schematic diagram of one embodiment of an implantable medical device (IMD) system in which techniques disclosed herein may be implemented.
FIG. 3 is a flow chart of a method for generating EGM-based index values according to one embodiment.
FIG. 4 depicts look-up table data that may be generated using the method shown in FIG. 3 according to one embodiment.
FIG. 5 is a flow chart of a method for controlling CRT control parameters in a closed-loop manner using an EGM-based index.
FIGS. 6A through 6D are a timing diagrams illustrating closed-loop control of an AV delay based on an EGM-based index.
Fundamentally, CRT alters electrical activation of the ventricles, improving spatial synchronization of electrical conduction in hearts with electrical conduction disorders such as left bundle branch block, right bundle branch block or other disorders. Optimal electrical activation of the heart may therefore be important for CRT efficacy. Optimal electrical activation can depend on a number of factors including the location of the pacing electrodes and pacing timing parameters such as AV delay and VV delay. Techniques disclosed herein enable an IMD to perform closed loop optimization of electrical activation of the heart.
IMD 10 senses RA EGM signals and delivers RA pacing pulses using RA lead 122, carrying tip electrode 148 and ring electrode 150. RA lead 122 is shown to be carrying coil electrode 166 which may be positioned along the superior vena cava (SVC) for use in delivering cardioversion/defibrillation shocks. In other embodiments, RV lead 118 carries both the RV coil electrode 162 and the SVC coil electrode 166. IMD 10 may detect tachyarrhythmias of heart 112, such as fibrillation of ventricles 128 and 132, and deliver high voltage cardioversion or defibrillation therapy to heart 112 in the form of electrical shock pulses. Pacing and sensing of the cardiac chambers is typically achieved using the pace/sense electrodes 140, 142, 144 148 and 150, however in some embodiments coil electrodes 162 and/or 166 may be used in sensing and/or pacing electrode vectors.
IMD 10 is configured for delivering CRT by delivering pacing pulses in one or both ventricles 128 and 132 for controlling and improving ventricular synchrony. LV pacing may be delivered using a selected pacing vector that utilizes at least one electrode 144 on multipolar LV lead 120. RV pacing is delivered using RV tip electrode 140 and ring electrode 142. CRT may be delivered by pacing in a single ventricular chamber (LV or RV) or both chambers (biventricular pacing) depending on patient need. The methods described herein are implemented in a dual or multi-chamber pacemaker or ICD delivering pacing pulses to the right and/or left ventricles using programmable pacing pulse timing parameters and selected pacing vectors.
In some embodiments, IMD 10 is configured to provide “adaptive CRT” which automatically switches between biventricular pacing and LV-only pacing in response to changes in the patient's intrinsic AV conduction. When AV conduction is impaired or blocked, or more generally when AV conduction time is slowed, biventricular pacing is delivered. When normal AV conduction returns, LV-only pacing is delivered. In this way, RV pacing is delivered only when needed based on the patient's own AV conduction status, which may fluctuate over time.
While a multi-chamber ICD is shown in FIG. 1, it is recognized that techniques disclosed herein may be implemented in a single chamber, dual chamber or multi-chamber pacemaker, with or without anti-arrhythmia therapies such as cardioversion and defibrillation shock capabilities. For example, techniques disclosed herein for closed-loop optimization of a CRT control parameter may be used to optimize the AV delay applied between an atrial event, sensed or paced, and a ventricular pacing pulse delivered in one ventricle (RV or LV only) or biventricular pacing pulses (RV and LV).
Programmer 170 includes a display 172, a processor 174, a user interface 176, and a communication module 178 including wireless telemetry circuitry for communication with IMD 10. In some examples, programmer 170 may be a handheld device or a microprocessor-based home monitor or bedside programming device. A user, such as a physician, technician, nurse or other clinician, may interact with programmer 170 to communicate with IMD 10. For example, the user may interact with programmer 170 via user interface 176 to retrieve currently programmed operating parameters, physiological data collected by IMD 10, or device-related diagnostic information from IMD 10. A user may also interact with programmer 170 to program IMD 10, e.g., select values for operating parameters of the IMD. A user interacting with programmer 170 can initiate a CRT optimization procedure performed by IMD 10 automatically or semi-automatically, to establish data for closed-loop optimization of CRT control parameters. As will be described herein, an EGM-based index value, corresponding to varying states of ventricular activation synchrony as assessed from surface ECG signals, is determined from an EGM signal sensed using selected cardiac electrodes 140, 142, 144, 162, 166 or housing electrode 158. The EGM-based index values are determined during CRT as it is delivered using varying parameter settings to establish EGM-based index data corresponding to differing states of ventricular dyssynchrony. The EGM-based index data is used for optimizing the control parameter in a closed-loop manner by IMD 10. This EGM-based index, derived from the EGM signal, is also referred to herein as an “EGM-based dyssynchrony index” because its relationship to CRT control parameter settings associated with optimal ventricular electrical synchrony and with non-optimal ventricular electrical synchrony is established. Knowledge of a patient-specific optimal electrical synchrony condition is established using surface ECG electrode signals, and the EGM-based index is recorded for an optimal CRT parameter setting, such as AV delay, as well as multiple non-optimal settings.
Programmer 170 includes a communication module 178 to enable wireless communication with IMD 10. Examples of communication techniques used by system 100 include low frequency or radiofrequency (RF) telemetry, which may be an RF link established via Bluetooth, WiFi, or MICS, for example as described in U.S. Pat. No. 5,683,432 (Goedeke, et al). In some examples, programmer 170 may include a programming head that is placed proximate to the patient's body near the IMD 10 implant site, and in other examples programmer 170 and IMD 10 may be configured to communicate using a distance telemetry algorithm and circuitry that does not require the use of a programming head and does not require user intervention to maintain a communication link.
It is contemplated that programmer 170 may be coupled to a communications network via communications module 178 for transferring data to a remote database or computer to allow remote monitoring and management of patient 114 using the techniques described herein. Remote patient management systems may be configured to utilize the presently disclosed techniques to enable a clinician to review CRT therapy parameters and authorize programming of IMD 10. Reference is made to commonly-assigned U.S. Pat. No. 6,599,250 (Webb et al.), U.S. Pat. No. 6,442,433 (Linberg et al.), U.S. Pat. No. 6,418,346 (Nelson et al.), and U.S. Pat. No. 6,480,745 (Nelson et al.) for general descriptions and examples of network communication systems for use with implantable medical devices for remote patient monitoring and device programming, all of which patents are hereby incorporated herein by reference in their entirety.
System 100 further includes an array of surface electrodes 182, which may be carried by a belt or strap 180 adapted to be wrapped around the torso of patient 114 to position electrodes 182 in the vicinity of heart 112. Strap 180 is shown inferior to heart 112 in FIG. 1, but it is understood that belt 180 may be positioned in a relatively more superior position to surround heart 112 such that electrodes 180 are positioned nearer to heart 112. Electrodes 182 are used to acquire surface signals from heart 112 during a CRT optimization session. A CRT control parameter, such as the AV delay and/or VV delay may be optimized by adjusting the parameter until the surface ECG-based determinations of ventricular activation indicate optimally synchronized ventricular activation. EGM data is then generated by IMD 10 at the optimized parameter setting and multiple increments/decrements from the optimal delay setting to establish EGM-based dyssynchrony index data and its relationship to increments or decrements from optimal control parameter settings, specific to the patient. The EGM-based data is used by IMD 10 to adjust the control parameter in a closed loop to maintain optimized ventricular activation.
In one example illustrated in FIG. 1, strap 180 is wrapped around the torso of patient 114 such that the electrodes 182 surround heart 112. Electrodes 182 may be positioned around the circumference of patient 114, including the posterior, lateral, and anterior surfaces of the torso of patient 114. In other examples, electrodes 182 may be positioned on any one or more of the posterior, lateral, and anterior surfaces of the torso. Electrodes 182 may be electrically connected to an ECG processing unit 184 via a wired connection 186. Some configurations may use a wireless connection to transmit the signals sensed by electrodes 182 to ECG processing unit 184, e.g., as channels of data.
Although in the example of FIG. 1, strap 180 is shown carrying surface electrodes 182, in other examples any of a variety of mechanisms, e.g., tape or adhesives, may be employed to aid in the spacing and placement of electrodes 182. In some examples, strap 180 may include an elastic band, strip of tape, or cloth. In some examples, electrodes 182 may be placed individually on the torso of patient 114.
Electrodes 182 may surround heart 112 of patient 114 and record the electrical signals associated with the depolarization and repolarization of heart 112. Each of electrodes 182 may be used in a unipolar configuration to sense the surface potentials that reflect the cardiac signals. ECG processing unit 184 may also be coupled to a return or indifferent electrode (not shown) which may be used in combination with each of electrodes 182 for unipolar sensing. In some examples, there may be 12 to 16 electrodes 182 spatially distributed around the torso of patient 114. Other configurations may have more or fewer electrodes 182.
ECG processing unit 184 may record and analyze the surface potential signals, referred to generally herein as “ECG” signals, sensed by electrodes 182. Processing unit 184 may be configured to provide an output to a user indicating electrical dyssynchrony of heart 112. The user may make a diagnosis, prescribe CRT, position therapy devices, e.g., leads, or adjust or select treatment parameters based on the indicated electrical dyssynchrony.
ECG processing unit 184 may compute activation times directly from sensed surface potential signals. An activation time for each electrode location (of electrodes 182) may be determined as a time period between two events, such as between the QRS complex onset and the minimum derivative during the QRS signal (i.e., the steepest negative slope of the sensed potential signal) at the respective electrode. Values of one or more indices indicative of the temporal and/or spatial distribution of the activation times may be determined as measures or indicators of electrical dyssynchrony. These indicators of electrical dyssynchrony may be used to evaluate different CRT control parameters and identify an optimal CRT control parameter.
Examples of indices of cardiac electrical dyssynchrony that may be calculated from surface potential signals sensed by electrodes 182 include a standard deviation of the determined activation times, a range of activation times, and a percentage of late activations. All or a subset of the surface electrodes (e.g., only electrodes located on the left anterior, left lateral and left posterior regions of the torso) may be used for calculation or computation of the activation times. The range of activation times may be computed as the difference between the maximum and the minimum cardiac activation times determined from all or a subset of electrodes 182. The percentage of late activations estimates the percentage of electrodes 182 whose associated activation times are greater than a certain percentile, for example the 70th percentile, of the QRS complex duration or the determined activation times for electrodes 182. Techniques for determining indices of electrical dyssynchrony based on surface activation times are generally disclosed in commonly-assigned pre-grant U.S. Patent Publication No. 2012/0283587 A1 (Ghosh, et al.) hereby incorporated herein by reference in its entirety. Indices of electrical dyssynchrony derived from external surface ECG leads are generally described. CRT optimization based on such indices derived from surface ECG leads can be performed at implant or at patient follow-up visits. Techniques disclosed herein, however, enable tuning of pacing timing parameters in an ongoing closed-loop manner to maintain optimal electrical activation of the ventricles in a patient-specific manner.
One or more indices of ventricular dyssynchrony based on the surface potential signals sensed by electrodes 182 is used to identify an optimal CRT parameter setting, such as an AV or VV delay. A user may program the AV delay into IMD 10 using programmer 170. In some embodiments, ECG processing unit 184 and programmer 170 are in wired or wireless communication or integrated in a common device that enables system 100 to automatically step through multiple CRT parameter settings, record and analyze surface potential signals to obtain one or more ECG-based indices of ventricular dyssynchrony, and identify and program an optimal setting for the CRT parameter based on analysis of ventricular electrical activations determined from surface ECG signals and the synchrony or dyssynchrony thereof.
The strap 180 carrying electrodes 182 is one illustrative embodiment of an apparatus that is useful in recording surface ECG signals from which ventricular activation times can be determined. Other surface cardiac signal recording apparatus may be used for acquiring cardiac signal data from which ventricular activation times can be computed and used in computing an index of ventricular dyssynchrony for establishing an optimal setting of one or more CRT control parameters. Other signal recording apparatus and techniques may include 12-lead ECG electrodes, a vest carrying an array of electrodes, and vectorcardiography.
Once an optimal CRT parameter is established based on optimal synchronization of electrical activation signals of the ventricles derived from surface ECG signals, CRT is delivered by IMD 10 using the optimal parameter setting and multiple non-optimal settings increased or decreased, i.e. shifted, from the optimal setting. IMD 10 acquires EGM signals for the optimal setting and multiple non-optimal settings to establish EGM parameter data for differing states of ventricular activation, i.e. different states of optimal electrical synchrony and non-optimal electrical synchrony corresponding to different increments and/or decrements from the optimal setting of the control parameters. This EGM parameter data is stored by IMD 10 and used in closed-loop CRT control parameter optimization as will be described herein.
FIG. 2 is a block diagram illustrating one example configuration of IMD 10. In the example illustrated by FIG. 2, IMD 10 includes a processor and control unit 80, also referred to herein as “processor” 80, memory 82, signal generator 84, sensing module 86, and telemetry module 88. IMD 10 further includes cardiac signal analyzer 90.
Processor and control unit 80 includes a therapy control unit that controls signal generator 84 to deliver electrical stimulation therapy, e.g., cardiac pacing or CRT, to heart 112 according to a selected one or more therapy programs, which may be stored in memory 82. Signal generator 84 is electrically coupled to electrodes 140, 142, 144A-144D (collectively 144), 148, 150, 158, 162, and 166 (all of which are shown in FIG. 1), e.g., via conductors of the respective leads 118, 120, 122, or, in the case of housing electrode 158, via an electrical conductor disposed within housing 160 of IMD 10. Signal generator 84 is configured to generate and deliver electrical stimulation therapy to heart 112 via selected combinations of electrodes 140, 142, 144, 148, 150, 158, 162, and 166. Signal generator 84 delivers cardiac pacing pulses according to AV and/or VV delays during CRT. These delays are set based on an analysis of cardiac signals by analyzer 90 as will be described herein.
The occurrence of R-waves in the ventricles, e.g. in the RV, may be used in monitoring intrinsic AV conduction time. In particular, prolongation of the AV conduction time or the detection of AV block based on R-wave sensing during no ventricular pacing (or pacing at an extended AV delay that allows intrinsic conduction to take place) is used to control adaptive CRT in some embodiments. When AV conduction is impaired, signal generator 84 is controlled by processor 80 to deliver biventricular pacing, i.e. pacing pulses are delivered in the RV and the LV using a selected AV delay and a selected VV delay. When AV conduction is intact, signal generator 84 is controlled by processor 80 to deliver LV-only pacing at a selected AV delay to optimally improve ventricular synchrony according to an EGM-based parameter whose relationship to ventricular electrical activation synchrony has been previously established.
As described herein, the AV delay may be optimized uniquely for different heart rhythm states such as rhythm states involving atrial sensing, atrial pacing, LV-only pacing, or biventricular pacing. For example, four distinct atrioventricular rhythm states may be evaluated: 1) atrial-sensed, biventricular paced 2) atrial-paced, biventricular paced 3) atrial-sensed, LV-only paced and 4) atrial-pace, LV-only paced. EGM-based ventricular dyssynchrony data may be established for different rhythm states and used to adjust the AV delay according to the EGM-based data and the current atrial sensing or pacing rhythm state and/or LV-only or biventricular pacing state.
Memory 82 stores intervals, counters, or other data used by processor 80 to control the delivery of pacing pulses by signal generator 84. Such data may include intervals and counters used by processor 80 to control the delivery of pacing pulses to one or both of the left and right ventricles for CRT. The intervals and/or counters are, in some examples, used by processor 80 to control the timing of delivery of pacing pulses relative to an intrinsic or paced event in another chamber. Memory 82 stores look-up tables and/or equations established for adjusting CRT control parameters such as AV and VV delays as will be described herein. Equations may be stored in the form of coefficient and intercept values defining a relationship between an EGM-based ventricular dyssynchrony parameter and different settings of a control parameter.
FIG. 3 is a flow chart 300 of a method for generating EGM-based ventricular dyssynchrony data for use in closed-loop adjustments of a CRT control parameter according to one embodiment. At block 302, ventricular dyssynchrony is assessed by analyzing activation times determined at multiple surface electrodes as described above and in the incorporated ‘587 published application. One or more indices of ventricular dyssynchrony are determined for each CRT parameter setting to be tested. Based on the analysis of activation times determined from surface electrodes, an optimal CRT parameter setting is identified at block 304. For example, a CRT control parameter resulting in a minimized activation time range, a minimum activation time standard deviation, and/or a minimized percentage of late activation times may be identified as an optimal CRT parameter.
One or more CRT parameters may be evaluated. Parameters that may be evaluated include AV delay and VV delay. Furthermore, each of these timing parameters may be evaluated during different rhythm conditions or states, e.g. during atrial sensing and during atrial pacing (with ventricular pacing occurring in one or both ventricles). In one embodiment, AV delay is set to multiple settings during an atrial pacing rhythm and during an atrial sensing rhythm to identify an optimal AV delay during atrial sensing (SAV delay) and during atrial pacing (PAV delay).
Additionally, AV delay may be set to multiple settings during pacing in one ventricle, RV-only and/or LV-only, and during biventricular pacing. During biventricular pacing, multiple settings of VV delay may be tested. For example, once an optimal AV delay is identified using a nominal VV delay, multiple VV delay settings may be applied to determine the optimal VV delay according to an ECG-based index of ventricular dyssynchrony.
In some embodiments, multiple pacing vectors may be available. For example, as shown in FIG. 1, a multi-polar CS lead may include multiple electrodes available for pacing the LV. Accordingly, an optimal AV delay resulting in a minimized index of ventricular dyssynchrony based on surface activation time determinations may be identified for each available pacing vector. When a pacing vector is changed, for example due to a change in lead impedance or other condition, the AV or VV delay may be adjusted to an optimal setting identified for the new pacing vector.
Once an optimal setting is established for a CRT parameter block 304, CRT is delivered at the optimal parameter setting(s) at block 306. An EGM signal is acquired by sensing module 86 and provided to cardiac signal analyzer 90. Signal analyzer 90 determines a ventricular dyssynchrony parameter from the EGM signal at block 308.
One or more sensing vectors may be selected from the available cardiac electrodes to provide EGM signals for determining a parameter of the EGM signal during different control parameter settings. In one embodiment, an EGM sensing vector selected for computing an EGM-based parameter for assessing a ventricular electrical dyssynchrony state is from an LV cathode electrode 144 to RV coil electrode 162 (shown in FIG. 1). An index may be computed as the time interval between two fiducial or reference points of the EGM signal. In one embodiment, the EGM-based index is determined as the time interval between a maximum amplitude and the minimum amplitude of an EGM signal during a monitoring window, e.g. a 200 ms window, immediately following a pacing pulse (or after a blanking interval following the pacing pulse). Alternatively, the ventricular dyssynchrony index may be determined as the time interval from a pacing pulse to a maximum amplitude of the EGM signal sensed from the selected sensing vector. The index may be an averaged value of multiple indices obtained from multiple cardiac cycles.
After recording the EGM-based index for the optimal CRT parameter setting at block 308, the parameter setting is increased or decreased at block 312 to a non-optimal setting by a known increment/decrement. The EGM-based index is determined for the new setting by returning to block 308. This process of increasing or decreasing the CRT parameter by known increments/decrements is repeated until an EGM-based index has been recorded for all desired CRT parameter test settings as determined at block 310.
The EGM-based index is stored for each test setting, along with the difference between the test (non-optimal) setting and the optimal setting for the CRT parameter at block 314. In this way, an EGM-based index is characterized for a known optimal ventricular activation condition, i.e. the optimal CRT parameter setting identified in response to surface ECG analysis, and for multiple non-optimal settings. In other words, a relationship is established between the EGM-based index and multiple parameter settings, including the optimal setting and one or more non-optimal settings different than, i.e. shifted from, the optimal setting.
Knowing the value of the EGM-based index during optimized electrical activation of the ventricles, adjustments to the CRT control parameter may be made to return the EGM-based index toward the value associated with electrical activation synchrony in a closed-loop control method. The optimal control parameter setting, such as AV delay, may change with changes in heart rate, activity or other conditions. This variation in an optimal setting occurs when intrinsic AV conduction timing changes. To maintain optimal ventricular activation under changing conditions, the optimal control parameter setting, like the A-V delay, needs to be adjusted so that the relationship of the timing of a CRT ventricular pacing pulse and the timing of intrinsic ventricular conduction remains consistent. However, to determine the intrinsic AV conduction and its changes directly, ventricular pacing needs to be inhibited temporarily, suspending CRT therapy. Even short disruptions in CRT therapy may be undesirable in some patients. By monitoring the EGM-based index and its changes during CRT pacing, it is possible to detect a need to adjust a timing control parameter without temporary suspension of CRT therapy.
Adjustments to control parameters, like the A-V delay, may be made based on the stored patient-specific relationship of the EGM-based index to increments and decrements of the control parameters. In this way, the EGM-based index can be restored to the value associated with optimal electrical activation synchrony and maintained at this value regardless of heart rate, intrinsic conduction changes or other changing conditions.
FIG. 4 depicts a look-up table 400 that may be generated using the method shown in FIG. 3 according to one embodiment. If AV delay is the CRT parameter being optimized, after setting an optimal AV delay based on a ventricular dyssynchrony index computed from ECG surface potential signals, the IMD 10 may be triggered to establish EGM-based index data, for example by a command received from programmer 170.
The signal generator 84 is controlled to deliver a ventricular pacing pulse to one or both ventricles using the previously established optimal AV delay. In the illustrative example, CRT is delivered during four different rhythm states 402, including single ventricular pacing (LV only) and biventricular pacing (LV and RV) during both atrial sensing and atrial pacing. Ventricular pacing pulses may be delivered simultaneously to the RV and LV during atrial sensing (SAV-BV) 402 a and during atrial pacing (PAV-BV) 402 c or at a nominal VV interval. For establishing the index data during biventricular pacing, a nominal V-V delay, e.g. 0 ms, or other previously established value may be used. An LV-only pulse is delivered during atrial sensing (SAV-LV) 402 b and during atrial pacing (PAV-LV) 402 d.
During pacing, the ventricular pacing pulses are delivered at varying AV delays. The cardiac signal analyzer 90 receives a selected EGM signal, e.g. from an LV electrode 144 to RV coil 162 sensing vector or any other available sensing vector, and determines an index from the EGM signal. In one embodiment the EGM-based index is the time interval between a maximum and minimum amplitude in a 200 ms window following the ventricular pacing pulse(s).
The EGM-based index 406 is computed during pacing using the optimal AV delay for each of the rhythm states and stored in the look-up table 400 with an associated 0 ms shift from the optimal AV delay setting. The AV delay shift 404 is the difference between a non-optimal AV delay and the optimal A-V delay. The EGM-based index is computed and stored for multiple increments and decrements relative to the optimal AV delay. For each AV delay setting, the EGM-based index 406 is stored for each corresponding AV delay shift 404 relative to the optimal AV delay. For example, the AV delay may be decreased by decrements of 10 ms to obtain EGM-based index values at AV delay shifts of −10, −20, −30, −40 and −50 ms relative to the optimal AV delay. The AV delay may be additionally increased by increments of 10 ms to obtain EGM-based indexes for each AV delay shift relative to the optimal AV delay of +10, +20, +30, +40 and +50 ms. The size of the AV delay shifts relative to the optimal AV delay and the number of different settings greater and less than the optimal AV delay may vary based on what the optimal delay is, what the intrinsic AV timing interval is for a particular patient, and what range and resolution of settings are available in the IMD.
The shift or difference between the optimal delay and a non-optimal delay is stored in the look up table 400 instead of the actual delay values since the optimal delay resulting in an optimized EGM-based index value will vary with changing heart rate, intrinsic conduction, or other factors. It is the relationship between the EGM-based index and relative shifts (increments or decrements) of the control parameter setting from an optimal setting that is established by collecting the data shown in look up table 400.
The EGM-based index 406 may be recorded in the look up table 400 in IMD memory 82 for each rhythm type 402 for multiple AV delay shifts 404 relative to the optimal AV delay. It is contemplated that EGM-based indexes may be determined for establishing look up table data relating the EGM-based index values to other CRT control parameters. For example, using an optimal AV delay setting, EGM-based index values can be obtained during pacing using an optimal VV delay (identified from surface ECG signal analysis) and multiple shifts from the optimal VV delay during atrial pacing and/or atrial sensing.
A look-up table of EGM-based index values for different AV and/or VV delay shifts may be established for multiple pacing vectors. For example different LV pacing vectors may be selected when a multipolar coronary sinus lead is used. A look-up table 400 may be stored in memory for each available pacing vector or a preferred subset of available pacing vectors.
FIG. 5 is a flow chart 500 of a method for controlling CRT control parameters in a closed-loop manner using EGM-based index data. After establishing one or more look up tables in IMD memory 82, processor 80 uses the look up table data for controlling signal generator 84 to adjust a CRT parameter in response to determinations of the EGM-based ventricular dyssynchrony index made during CRT delivery. In the illustrative example, the CRT parameter being adjusted is AV delay, however, VV delay could be controlled in a closed loop manner based on look up table data and EGM-based index determinations made during CRT in a similar manner.
At block 502, CRT is delivered using an optimal parameter setting, e.g. AV delay for the current rhythm state (atrial sensing or pacing and LV-only or biventricular pacing) and currently selected pacing vector. The EGM-based index of ventricular dyssynchrony is monitored during CRT delivery at block 504. For example, sensing module 86 may provide cardiac signal analyzer 90 EGM signal data on every ventricular-paced cardiac cycle to enable continuous monitoring of the ventricular dyssynchrony index. Alternatively, the EGM-based index may be determined every Nth paced beat, at regular time intervals or according to any desired monitoring schedule. In some embodiments, the EGM-based index is at least determined each time a change in a CRT pacing occurs, including changes in atrial rate. For example, if the rhythm changes from a sensed atrial rhythm to a paced atrial rhythm or vice versa, if the rhythm changes between biventricular pacing and single ventricle pacing, if the heart rate changes, or if a pacing vector changes, the ventricular synchrony condition or state is monitored by determining an index from the EGM signal to determine if an adjustment to the AV delay is needed.
At block 506, the value of the EGM-based index determined during CRT delivery is looked up in a table 400 for the current rhythm state (e.g. SAV-BV, PAV-LV, SAV-BV, or PAV-BV) and stored for the selected pacing vector. If the index is approximately equal to the optimal value stored for a 0 ms AV delay shift, the AV delay remains unchanged and monitoring continues at block 504. If the EGM-based index has changed from the optimal value, indicating a change in ventricular dyssynchrony, an AV delay adjustment is determined at block 508. The AV delay adjustment is determined from the look up table 400 by looking up the EGM-based index value in the relevant portion of table 400 for the current rhythm state 402 and selected pacing vector. The AV delay adjustment is determined as the same magnitude but opposite polarity or direction of the AV delay shift associated with the current value of the EGM-based dyssynchrony index.
For example, in reference to look up table 400 of FIG. 4, if the current rhythm is SAV-BV 402A, and the EGM-based index is determined as being 34 ms, the associated AV delay shift is +30 ms. Since this value of the EGM-based index occurred for an AV delay that was shifted +30 ms from the optimal AV delay setting, the AV delay adjustment required to improve ventricular synchrony is −30 ms. By adjusting the current AV delay setting by −30 ms, the EGM-based index is expected to return to approximately 50 ms, corresponding to optimally synchronized electrical activation of the ventricles.
In some cases, an EGM-based index value 406 may be associated with more than one AV delay shift 404 for a given rhythm 402. For example, in look up table 400, the EGM-based dyssynchrony index is 80 ms for an AV delay shift of −40 ms and −50 ms during an SAV-BV rhythm 402 a. In this case, the smaller adjustment may be selected first. For example, the AV delay may be adjusted by +40 ms in response to a determination of the EGM-based dyssynchrony index being 80 ms. Continued monitoring of the EGM-based dyssynchrony index will result in one or more additional adjustments of the AV delay if the first adjustment does not restore the EGM-based index to its optimal value, e.g. 50 ms in the illustrative example. The AV delay is adjusted at block 510 by the adjustment determined at block 508.
If a change in the EGM-based index is detected at block 512 in response to the adjusted parameter, CRT delivery may continue using the adjusted control parameter. The EGM-based index continues to be monitored at block 504 and readjustments to the control parameter are made as needed.
In some cases, the EGM-based index may not change in response to multiple adjustments to a control parameter. In such cases, intrinsic AV conduction block may be detected. In response to no change in the EGM-based index (block 512) following an adjustment of an AV delay, the processor and control 80 may control the signal generator 84 to alter the pacing therapy, e.g. by switching a pacing mode, because of the loss of intrinsic AV conduction. For example, AV block may be detected as evidenced by an EGM-based index value that remains substantially unchanged in response to a threshold number of adjustments to AV delay, e.g. three or more adjustments. If LV-only pacing is being delivered, as determined at block 514, the signal generator 84 may be switched to a biventricular pacing mode at block 516. Biventricular pacing may be delivered at a nominal AV delay setting since the EGM-based index is unresponsive to changes in AV delay. The nominal value of AV delay may range from 80 to 200 ms for sensed A-V delay and 120 to 250 ms for paced AV delay. In another embodiment, the IMD may deliver biventricular pacing at the optimal AV delay that was determined initially. In yet another embodiment the device may deliver pacing at a nominal AV delay without switching pacing mode.
The processor 80 may record and store the corresponding value of the EGM-based index at the onset of biventricular pacing (or nominal AV delay pacing without switching pacing mode) at block 518. Processor 80 continues to monitor the EGM-based index for any changes in the EGM index value.
If a change in the EGM-based index is detected during subsequent monitoring at block 520, (indicating a possible return of intrinsic AV conduction), the processor 80 may control the signal generator 84 to switch back to LV-only pacing at block 522 (if a pacing mode switch to biventricular pacing had been made previously). The EGM-based index is monitored again relative to its optimal value at blocks 504 and 506. Control parameter adjustments continue as needed to restore and maintain the EGM-based index at the established value corresponding to optimal electrical activation.
Using the disclosed techniques, a value for an EGM-based index that corresponds to a known state of optimal ventricular electrical synchrony, as initially determined by ECG-based surface potential analysis, is established for one or more rhythm states and/or one or more pacing vectors. Values of the index associated with non-optimized control parameter setting shifts, resulting in various states of ventricular dyssynchrony, are also established. This enables adjustments to a control parameter to be determined when deviations occur in the value of the EGM-based index from its established value at an optimal control parameter setting. An EGM-based index that deviates from the optimal value during a given rhythm state can be corrected and restored by adjusting a CRT timing control parameter based on an established relationship between the EGM-based index and control parameter shifts. This relationship may be stored in the form of a look up table as shown in FIG. 4 or in the form of a best fit equation to the data. An adjustment to a CRT control parameter is determined in response to a change in the EGM-based index by computing a shift from a best fit equation or looking up the associated shift in a look up table. The timing control parameter adjustment is the same magnitude but opposite polarity of the shift relative to the optimal setting.
FIGS. 6A through 6D are timing diagrams illustrating closed-loop control of the AV delay based on an EGM-based index. In FIG. 6A, timing diagram 600 shows the timing of an intrinsic atrial sense event 602 followed by the time of an expected intrinsic ventricular sense event 604. An optimal AV delay 606 for delivering biventricular pacing pulses 608 during atrial sensing (SAV) is identified as 100 ms based on analysis of activation times determined from surface ECG signals. This optimal SAV delay 606 may be established during a resting heart rate. For this optimal SAV delay 606, the EGM-based index P is determined to be 50 ms, which may be a time interval determined between a maximum and minimum peak amplitude of the EGM signal following a pacing pulse. This determination establishes the expected value of the EGM-based index when the ventricular activation is known to be in an optimized state of synchrony based on surface ECG activation times.
FIG. 6B shows a timing diagram 610 for one example of a non-optimal SAV delay 616 applied during the same resting heart rate and biventricular pacing conditions as timing diagram 600. The SAV delay has been increased to 140 ms, a shift of +40 ms relative to the optimal SAV delay. Biventricular pacing pulses 618 are delivered after the shifted SAV delay 616 following a sense atrial event 612. As can be seen, the biventricular pacing pulses are delivered just prior to an expected (but unknown) time of an intrinsic ventricular depolarization.
The EGM-based index (P) determined at this shifted SAV delay 616 is 30 ms. Thus, in a look up table, an EGM-based index value of 30 ms is stored for an SAV-BV rhythm and delay shift of +40 ms.
FIG. 6C is a timing diagram 620 during an increased heart rate, for example during exercise. The intrinsic ventricular depolarization 624 would occur at a shorter intrinsic AV interval after the atrial sense event 622. However, the timing of the intrinsic ventricular depolarization 624 is unknown since the ventricles are being paced by biventricular pacing pulses 628. Without having to detect a change in heart rate or a change in the intrinsic AV interval, the IMD processor 80 can adjust the SAV-BV delay 626 automatically by determining the EGM-based index (P) of ventricular dyssynchrony. For the higher heart rate condition shown in FIG. 7C, the SAV delay of 100 ms, which was optimal at a resting heart rate, results in an EGM-based index of ventricular dyssynchrony equal to 30 ms. This change in the EGM-based index from the optimal 50 ms to 30 ms indicates non-optimal ventricular electrical activation.
FIG. 6D is a timing diagram 630 illustrating a closed loop adjustment of the SAV delay. The EGM-based index value of 30 ms determined during the SAV-BV rhythm shown in FIG. 6C corresponds to an SAV delay shift of +40 ms based on the previously obtained data shown in FIG. 6B and stored in the form of a look up table in IMD memory 82. For this EGM-based index value, an AV delay adjustment 640 of −40 ms (same magnitude but opposite polarity of the SAV delay shift in FIG. 6B) is indicated. In order to restore the EGM-based index to its optimal value of 50 ms, or as close as possible to the optimal value, the SAV delay 636 is adjusted by −40 ms, to 60 ms. The EGM-based index may be re-determined to verify that it has returned to a value of approximately 50 ms. A relative increment or decrement of the AV delay that is needed to restore an optimal EGM-based index value is ascertained from the look up table and can be applied regardless of the actual AV delay setting at the time of the deviated, non-optimal EGM-based index determination.
The relationship between the EGM-based index and ventricular electrical dyssynchrony as a CRT control parameter changes is established by determining the EGM-based index for a known optimal ventricular activation condition and recording the index for multiple conditions shifted from the optimal ventricular activation condition. By establishing this relationship, determinations of the EGM-based index during CRT delivery enable closed loop adjustment of a CRT control parameter without additional sensors and without interrupting CRT therapy, for example to determine an intrinsic AV conduction time. In this way, CRT control parameters are adjusted to restore the EGM-based index to a value known to correspond to an optimal electrical activation state, i.e. optimal ventricular electrical synchrony.
Thus, various embodiments of an IMD system and method for closed loop adjustment of a CRT control parameter have been described. However, one of ordinary skill in the art will appreciate that various modifications may be made to the described embodiments without departing from the scope of the claims. For example, although specific examples of indices of ventricular dyssynchrony determined from surface ECG activation time and from an EGM signal, other indices or determinations of ventricular dyssynchrony may be determined from ECG and EGM signals according the techniques of this disclosure. Furthermore, although an EGM-based index has been described herein, it is recognized that an index may be determined from other cardiac signals and corresponding index data may be determined in the manner described herein for establishing a relationship between the index and an optimal ventricular synchrony condition and multiple non-optimal states. Adjustments to pacing therapy control parameters may be made to maintain the index at or near an established optimal value. These and other examples are within the scope of the following claims.
1. A method for optimizing a pacing control parameter, comprising:
determining a value of an EGM-based index corresponding to an optimal electrical activation condition of a patient's heart;
delivering a pacing therapy to the heart;
monitoring the EGM-based index during delivery of the pacing therapy to the heart; and
adjusting the pacing therapy to maintain the EGM-based index at the value corresponding to the optimal electrical activation condition in response to detection of a value of the EGM-based index corresponding to a non-optimal electrical activation condition during the monitoring of the EGM-based index.
establishing a relationship between the EGM-based index and a plurality of settings of a pacing control parameter, wherein establishing the relationship comprises:
delivering cardiac pacing according to each of the plurality of settings comprising a previously established optimal setting and at least one non-optimal setting different than the optimal setting by a setting shift;
determining a value of the EGM-based index from an EGM signal acquired during cardiac pacing at each of the plurality of settings; and
storing the values of the determined index with corresponding setting shifts of each of the plurality of settings relative to the optimal setting; and
enabling a processor of an implantable medical device to:
monitor the EGM-based index during cardiac pacing,
responsive to detecting an EGM-based index value corresponding to a non-optimal setting of the control parameter, determine an adjustment of the control parameter from the stored index values and corresponding setting shifts, and
adjust the control parameter by the determined adjustment.
3. The method of claim 2, further comprising establishing the optimal setting by sensing surface potentials from a plurality of surface electrodes positioned on a patient's body and computing an index of ventricular dyssynchrony for each of a plurality of pacing control parameter settings in response to the sensed surface potentials.
4. The method of claim 2, wherein acquiring the EGM signal comprises selecting a sensing vector comprising a left ventricular electrode and a right ventricular coil electrode.
5. The method of claim 2, further comprising determining a magnitude and polarity of a setting shift corresponding to the detected EGM-based index corresponding to a non-optimal setting; and
determining the adjustment as the magnitude and a polarity opposite to the setting shift corresponding to the detected EGM-based index.
6. The method of claim 2, wherein establishing the optimal pacing control parameter comprises establishing an optimal atrio-ventricular (AV) delay during cardiac resynchronization therapy.
7. The method of claim 2, further comprising establishing an optimal pacing control parameter and relationship between the EGM-based index and a plurality of settings of the control parameter for each of a plurality of different heart rhythm states.
8. The method of claim 2, further comprising establishing an optimal pacing control parameter and a relationship between the EGM-based index and a plurality of settings of the control parameter for each of a plurality of pacing electrode vectors.
9. The method of claim 2, further comprising storing the relationship as a look-up table of determined EGM-based index values and associated setting shifts in a memory of an implantable medical device accessed by the processor.
detecting an unchanged value of the EGM-based index in response to an adjustment of the control parameter;
and switching a pacing mode in response to the unchanged value.
11. The method of claim 1, wherein determining the EGM-based index comprises detecting a maximum peak amplitude of the EGM signal following a pacing pulse.
12. The method of claim 1, wherein determining the EGM-based index comprises determining a time interval between a maximum peak amplitude and a minimum peak amplitude of the EGM signal following a pacing pulse.
13. A non-transitory, computer-readable storage medium comprising instructions that, when executed, cause a processor included in a medical device system to:
determine a value of an EGM-based index corresponding an optimal electrical activation condition of a patient's heart;
control a signal generator to deliver a pacing therapy to the heart;
monitor the EGM-based index during delivery of the pacing therapy to the heart; and
adjust a pacing therapy control parameter to maintain the EGM-based index value corresponding to the optimal electrical activation condition in response to detection of a value of the EGM-based index corresponding to a non-optimal electrical activation condition during the monitoring of the EGM-based index.
14. The computer-readable storage medium of claim 13 further comprising instructions that, when executed, cause the processor to:
establish a relationship between an EGM-based index and a plurality of settings of a pacing control parameter, wherein establishing the relationship comprises:
delivering cardiac pacing pulses according to each of the plurality of settings comprising a previously established optimal setting and at least one non-optimal setting different than the optimal setting by a setting shift;
determining a value of the EGM-based index from the EGM signal acquired during cardiac pacing at each of the plurality of settings; and
storing the values of the determined index with corresponding setting shifts of each of the plurality of settings relative to the optimal setting;
monitor the EGM-based index during cardiac pacing;
responsive to detecting an EGM-based index value corresponding to a non-optimal setting of the control parameter, determine an adjustment of the control parameter from the stored index values and corresponding setting shifts; and
US13/833,882 2013-03-15 2013-03-15 Closed loop optimization of control parameters during cardiac pacing Active US9278219B2 (en)
US13/833,882 US9278219B2 (en) 2013-03-15 2013-03-15 Closed loop optimization of control parameters during cardiac pacing
US15/062,945 US9649497B2 (en) 2013-03-15 2016-03-07 Closed loop optimization of control parameters during cardiac pacing
US15/062,945 Continuation US9649497B2 (en) 2013-03-15 2016-03-07 Closed loop optimization of control parameters during cardiac pacing
US20140277233A1 US20140277233A1 (en) 2014-09-18
US9278219B2 true US9278219B2 (en) 2016-03-08
ID=51531184
US13/833,882 Active US9278219B2 (en) 2013-03-15 2013-03-15 Closed loop optimization of control parameters during cardiac pacing
US15/062,945 Active US9649497B2 (en) 2013-03-15 2016-03-07 Closed loop optimization of control parameters during cardiac pacing
US (2) US9278219B2 (en)
2013-03-15 US US13/833,882 patent/US9278219B2/en active Active
2016-03-07 US US15/062,945 patent/US9649497B2/en active Active
"Analysis of Ventricular Activation Using Surface Electrocardiography to Predict Left Ventricular Reverse Volumetric Remodeling During Cardiac Resynchronization Therapy", by Michael O. Sweeney, M.D., et al., Circulation Journal of the American Heart Association, Jan. 25, 2010, 121; pp. 626-634.
"Electrical and Mechanical Ventricular Activation During Left Bundle Branch Block and Resynchronization", by Marc Strik, et al., J. of Cardiovasc. Trans. Res. (2012) 5:117-126.
"Electrophysiological Substrate and Intraventricular LV dyssynchrony in Non-ischemic Heart Failure Patients Undergoing Cardiac Resynchronization Therapy", by Subham Ghosh et al., The Heart Rhythm Society, May 2011; 8 (5): 692-699.
"Mathematical Modeling and Simulation of Ventricular Activation Sequences: Implications for Cardiac Resynchronization Therapy", by Mark Potse, J. Cardiovasc Transl Res., 2012 Apr, 5(2): 146-158.
"State-of-the-Art of Current Applications, Guidelines, Ongoing Trials, and Areas of Controversy", Arrhythmia Management Devices, Cardiac Resynchronization Therapy, Frits W. Prinzen, PhD., et al., 2013; 128: 2407-2418.
Cuculich, P.S., et al., "The Electrophysiological Cardiac Ventricular Substrate in Patients After Myocardial Infection" J. Am. Coll. Cardiol. 2011; 58:1893-1902.
Czerwinska et al., "Method of Segmentation of Thorax Organs Images Applied to Modeling the Cardiac Electrical Field," Engineering in Medicine and Biology Society, Proceedings of the 22nd Annual International Conference of the IEEE, vol. 1, 23, Jul. 23, 2000; pp. 402-405.
Dawoud, F. et al., "Inverse Electrocardiographic Imaging to Assess Electrical Dyssynchrony in Cardiac Resynchronization Therapy Patients," Computing in Cardiology, 2012; 39:993-996.
Hansen, "Regularization Tools: A Matlab Package for Analysis and Solution of Discrete Ill-Posed Problems," Version 4.1 for Matlab 7.3; Mar., 2008; 128 pages. Retrieved from the Internet: Jun. 19, 2014 http://www.mathworks.com/matlabcentral/fileexchange/52-regtools.
US9649497B2 (en) 2017-05-16
US20140277233A1 (en) 2014-09-18
US20160184590A1 (en) 2016-06-30
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:GHOSH, SUBHAM;REEL/FRAME:030011/0698