US20060189876A1 - Method and system for analyzing an electrocardiographic signal - Google Patents

Method and system for analyzing an electrocardiographic signal Download PDF

Info

Publication number
US20060189876A1
US20060189876A1 US11/407,902 US40790206A US2006189876A1 US 20060189876 A1 US20060189876 A1 US 20060189876A1 US 40790206 A US40790206 A US 40790206A US 2006189876 A1 US2006189876 A1 US 2006189876A1
Authority
US
United States
Prior art keywords
electrocardiogram signal
set forth
reach
time
beat
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.)
Abandoned
Application number
US11/407,902
Inventor
Jean-Philippe Couderc
Arthur Moss
Wojciech Zareba
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to US11/407,902 priority Critical patent/US20060189876A1/en
Publication of US20060189876A1 publication Critical patent/US20060189876A1/en
Abandoned legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4836Diagnosis combined with treatment in closed-loop systems or methods
    • A61B5/4839Diagnosis combined with treatment in closed-loop systems or methods combined with drug delivery
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/352Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/366Detecting abnormal QRS complex, e.g. widening
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals

Definitions

  • the present invention is directed to a method and system for the quantification of repolarization changes when measured in a dynamic electrocardiogram (“ECG”) signal.
  • ECG electrocardiogram
  • An ECG signal represents changes in electrical potential produced by contractions of the heart recorded from the surface of the body.
  • An example of an ECG signal is illustrated in FIG. 1 .
  • Each P-QRS-T complex reflects the electrical depolarization and repolarization components of the heart beat.
  • the cardiac beat comprises three major waves identified as: the P wave representing the depolarization of the auricles, the QRS complex generated by the depolarization of both right and left ventricles, and the T-wave which represents the repolarization of the ventricles.
  • the T wave also includes an apex, T apex , and an endpoint, T offset .
  • the repolarization of the auricles is hidden in the QRS complex.
  • the RR intervals are defined by the interval between two consecutive R peaks in successive beats.
  • the RR interval is a direct measurement of the heart rate.
  • an ECG signal a modification of the T wave morphology can be observed when certain pharmaceutical agents are taken and also with certain types of cardiac disorders.
  • prolongation of the QT interval in an ECG signal is clearly associated with an increased risk for ventricular arrhythmias and sudden cardiac death.
  • prior systems have been developed to monitor the QT interval in ECG signals to detect potential heart problems.
  • a method and a computer readable medium with programmed instructions for analyzing an ECG signal in accordance with an embodiment of the present invention includes obtaining a measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal and detecting an altered ventricular repolarization based on the obtained measured area based repolarization interval.
  • a system for analyzing an ECG signal in accordance with an embodiment of the present invention includes a measurement system and a detection system.
  • the measurement system obtains a measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal.
  • the detection system detects an altered ventricular repolarization based on the obtained measured area based repolarization interval.
  • a method and a computer readable medium with programmed instructions for analyzing an effect of a pharmacological agent on an electrocardiogram signal in accordance with an embodiment of the present invention includes obtaining a first measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal.
  • the pharmacological agent is administered.
  • a second measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal is obtained from at least one of a first period during which the pharmacological agent is in effect and a second period after the pharmacological agent is no longer in effect.
  • An altered ventricular repolarization is detected based on the first and second measurements.
  • a system for analyzing an effect of a pharmacological agent on an electrocardiogram signal in accordance with an embodiment of the present invention includes a measurement system and a detection system.
  • the measurement system obtains a first measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal before administering the pharmacological agent.
  • the measurement system also obtains a second measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal from at least one of a first period during which the pharmacological agent is in effect and a second period after the pharmacological agent is no longer in effect.
  • the detection system detects an altered ventricular repolarization based on the first and second measurements.
  • the present invention provides an effective non-invasive method for dynamic quantification of ventricular repolarization.
  • the present invention's comprehensive evaluation of ventricular repolarization provides enhanced utility: in the clinical diagnosis of acquired and inherited repolarization disorders; in the evaluation of new chemical entities for pharmaceutical companies during drug development and during post-marketing surveillance of approved drugs that can adversely affect ventricular repolarization; and in the clinical identification of patients with potentially life-threatening ventricular repolarization disorders, such as those after myocardial infarction.
  • FIG. 1 is a diagram of a beat in an ECG signal where two vertical lines mark the apex and the offset of the T wave used to measure QT intervals;
  • FIG. 2 is a block diagram of a quantification system in accordance with one embodiment of the present invention.
  • FIG. 3A is a diagram of another ECG signal
  • FIG. 3B is a graph of a normalized absolute integral of the repolarization segment of the ECG signal shown in FIG. 3A ;
  • FIG. 4 is a flow chart of a method for analyzing an ECG signal in accordance with another embodiment of the present invention.
  • FIGS. 5A-5C are graphs of the magnitude, phase, and zeros and poles of a filter for the quantification system
  • FIG. 6A is a graph of six area based repolarization intervals of median beats which are superimposed
  • FIG. 6B is a graph of cumulative area of the six area based repolarization intervals shown in FIG. 6A ;
  • FIG. 7 is a table of RTx % and RR interval bins
  • FIG. 8A is a graph of a heart rate of one patient showing beats per minute (BPM) against time;
  • FIG. 8B is a graph of a RT intervals (Apex and offset of T waves using median value across leads) of the one patient;
  • FIGS. 9A-9D are graphs of trends for T-wave amplitude for the median across leads, the X leads, the Y leads, and the Z leads;
  • FIGS. 10A-10D are graphs of the trends for area-based measurements for 25% of the area for the median, , the X leads, the Y leads, and the Z leads;
  • FIGS. 11A-11D are graphs of the trends for area-based measurements for 50% of the area for the median, , the X leads, the Y leads, and the Z leads;
  • FIGS. 12A-12D are graphs of the trends for area-based measurements for 90% of the area for the median, , the X leads, the Y leads, and the Z leads;
  • FIGS. 13A-13D are graphs of the trends for area-based measurements for 97% of the area for the median,, the X leads, the Y leads, and the Z leads;
  • FIG. 14 is a table of averages for 24 hours, during the day and at night for: beats per minute (BPM); X, Y, and Z amplitudes; X RT apex, Y RT apex, and Z RT apex; X RT offset, Y RT offset, and Z RT offset; XT 25% area, YT 25% area and ZT 25% area; XT 50% area, YT 50% area and ZT 50%; area XT 90% area, YT 90% area and ZT 90%; and area XT 97% area, YT 97% area and ZT 97% area;
  • BPM beats per minute
  • FIG. 15 is an hourly table of heart rate (BPM), T wave ampl X, T wave ampl Y, and T wave ampl Z;
  • FIG. 16 is an hourly table of heart rate (BPM), T wave apex X, T wave apex Y, and T wave apex Z;
  • FIG. 17 is an hourly table of heart rate (BPM), T wave offset lead X, T wave offset lead Y, and T wave offset lead Z;
  • FIG. 18 is an hourly table of heart rate (BPM), ABRI RT25 X-lead, ABRI RT25 Y-lead, and ABRI RT25 Z-lead;
  • FIG. 19 is an hourly table of heart rate (BPM), ABRI RT50 X-lead, ABRI RT50 Y-lead, and ABRI RT50 Z-lead;
  • FIG. 20 is an hourly table of heart rate (BPM), ABRI RT90 X-lead, ABRI RT90 Y-lead, and ABRI RT90 Z-lead;
  • FIG. 21 is an hourly table of heart rate (BPM), ABRI RT97X-lead, ABRI RT97 Y-lead, and ABRI RT97 Z-lead;
  • FIGS. 2 and 4 A system 10 and method for analyzing an ECG signal in accordance with one embodiment of the present invention are illustrated in FIGS. 2 and 4 .
  • the system 10 includes a quantification system 12 coupled to an ECG device 24 and a storage device 26 .
  • the method includes filtering and then extracting a portion of the filtered ECG signal, substantially removing any baseline from the extracted portion of the ECG signal, obtaining a measurement of at least one area based repolarization interval in the extracted portion, and then detecting any altered ventricular repolarization based on the measured area based repolarization interval.
  • the present invention provides a number of advantages including providing an effective non-invasive method for dynamic quantification of ventricular repolarization.
  • the quantification system 12 has a central processing unit (“CPU”) or processor 14 , a memory 16 , an input/output interface 22 , a display 18 , and a user input device 20 which are coupled together by a bus system 23 or other link, although the quantification system 12 may comprise other types of components, other numbers of the components, and other combinations of the components.
  • the processor 14 executes a program of stored instructions for at least a portion of the method for analyzing an ECG signal in accordance with one embodiment of the present invention as described herein.
  • the preprocessing and processing steps are carried out by quantification system 12 , although other types and combinations of systems and/or devices can be used to execute the preprocessing and processing steps.
  • Memory 16 comprises a random access memory (RAM) and a read only memory (ROM), although other types and combinations of memory storage devices can be used, such as a floppy disk, hard disk, CD ROM, or other computer readable medium which is read from and/or written to by a magnetic, optical, or other reading and/or writing system that is coupled to the processor 14 .
  • RAM random access memory
  • ROM read only memory
  • other types and combinations of memory storage devices can be used, such as a floppy disk, hard disk, CD ROM, or other computer readable medium which is read from and/or written to by a magnetic, optical, or other reading and/or writing system that is coupled to the processor 14 .
  • the method in accordance with one embodiment of the present invention is stored as programmed instructions in memory 16 in the quantification system 12 for execution by the processor 14 , some or all of the programmed instructions could be stored elsewhere.
  • the display 18 enables an operator to observe information, such as the ECG reading for a patient.
  • the display 18 is a cathode ray tube device, although other types of displays can be used, such as a printer.
  • the user input device 20 enables an operator to generate and transmit signals or commands to the processor, such as a request to obtain or retrieve ECG signals for a particular patient for processing.
  • the user input device 20 is a keyboard, although other types of user input devices can be used, such as a computer mouse.
  • the input/output interface 22 is used to operatively couple the quantification system 12 to the ECG device 24 and the storage device 26 , although quantification system 12 can be coupled via input/output interface to other systems and/or devices.
  • the ECG device 24 records one or more ECG signals from one or more patients.
  • the ECG device can capture the electrical activity from the heart in a variety of different manners, such as Holter recordings, exercise ECG testing, bedside ECG monitoring, event monitoring, implantable ECG recorders using one or more leads and with a duration that may vary from few minutes to twenty-four hours, even few days.
  • a variety of different types of ECG devices can be used for ECG device 24 .
  • the ECG information can be transmitted to the quantification system 12 from the ECG device 24 for processing.
  • the storage device 26 comprises a RAM, although other types and combinations of memory storage devices can be used, such as a ROM or a floppy disk, hard disk, CD ROM, or other computer readable medium which is read from and/or written to by a magnetic, optical, or other reading and/or writing system.
  • the storage device 26 can store ECG readings for processing by the quantification system 12 and can store the results of any processing.
  • FIG. 3A an example of another ECG signal is illustrated.
  • the Q peak is reached at time Q 0
  • the T wave begins at time S 0
  • the T apex is reached at T m
  • the endpoint of the T wave is reached at T o .
  • FIG. 3B a graph of the measured area under the repolarization interval shown in FIG. 3A with the Y axis showing the percentage of the total area under the repolarization interval and the X-axis showing time is illustrated.
  • the graph begins at time S 0 which again corresponds to the time when the T wave begins and U 0 represents the time when substantially 100% of the area under the repolarization area is measured.
  • the times to accumulate 25%, 50% and 75% of the absolute T-wave area are identified as LA25, LA50, LA75, respectively.
  • the present invention does not rely upon accurately determining either the T apex, T m , or the endpoint of the T wave, T o .
  • the present system identify the end of the T wave when the slope of the repolarization signal is below a given threshold.
  • the present invention quantifies changes in ventricular repolarization based on measurements of the area under the repolarization interval obtained starting at a time S 0 , which is a first period of time after the R peak is detected.
  • the morphology of the repolarization T wave in an ECG signal can be influenced by any factor that alters ion-channel kinetics involved in ventricular repolarization, including inherited diseases, acquired cardiac diseases, metabolic disturbances, and drugs. Since the normal morphology of the ventricular repolarization interval is well known to those of ordinary skill in the art, an altered or abnormal morphology ventricular repolarization interval can be both visually and quantitatively assessed if it is accurately measured. A process for analyzing ventricular repolarization in accordance with the present invention is described below.
  • the quantification system 12 obtains an ECG signal or signals.
  • the ECG signals can be obtained from a variety of different sources, such as directly from the ECG device 24 or from the storage device 26 which stores previously recorded ECG signals.
  • step 30 the ECG signal undergoes low pass-filtering in quantification system 12 .
  • the filtering technique used should preserve the initial shape of the T wave.
  • the ECG signal is filtered in the forward direction, the filtered ECG signal is then reversed and run back through the filter again.
  • the result of this filtering operation has precisely zero phase distortion and magnitude modified by the square of the filter's magnitude response.
  • a Butterworth filter and a bidirectional filtering technique are used in order to obtain a high quality filtering free of distorted components, although other types of filters and filtering techniques can be used.
  • a filter of order two with a cutoff frequency at ⁇ 3 db equal to 18 Hz is used.
  • Graphs of the magnitude, phase (radians) and zeros and poles of this filter used for removing noise components of the ECG signal and reducing morphological changes of area based repolarization interval are illustrated in FIGS. 5A-5C .
  • a portion or limited length sequence of the filtered ECG signal is extracted by the quantification system 12 . More specifically, the ECG signal is cut into portions or sequences of consecutive continuous beats by quantification system 12 , excluding non-normal beats. In each of these portions, R(k) is the R peak location in the k th beat in each portion and ECG(n) is the n th portion of ECG signal where n represents the number of portions of the ECG signal.
  • Each ECG portion which is analyzed is ECG(n) for n ⁇ [R(k), R(k+NB)] where NB is the number of beats in each analyzed interval.
  • a repolarization segment is extracted for each of seventy consecutive beats in the filtered ECG signal.
  • the baseline B(k) of the extracted portion of the ECG signal ECG(n) is substantially removed by the quantification system 12 .
  • the baseline removal must be accurate and must avoid adding significant components to the repolarization segment in the ECG signal ECG(n).
  • the baseline removal can be obtained based using a Cubic-Spline interpolation of the baseline of the ECG signal ECG(n), although other techniques for removing the baseline can be used.
  • the interpolation used in this particular example is based on the average amplitude value of the samples from a limited squared time window W(n) of 80 ms length located just before the beginning of the QRS complex.
  • the signal B(k) is interpolated at a sampling period equal to the sampling period of the ECG.
  • the Cubic Spline interpolation of B(k) is used to insure continuous smoothness up to the second derivative and is less likely to generate erratic oscillations than polynomial interpolation.
  • the interpolated curve, corresponding to the very low components of the ECG signal ECG(n) is subtracted from the ECG signal ECG(n) in order to reduce the effect of respiratory components of the signal.
  • two beats are removed from the beginning and two beats from the end of the analyzed intervals in the ECG signal ECG(n).
  • two beats are removed from the extremities of the ECG signal with seventy beats resulting in sixty-six beats which are analyzed. This insures elimination of undesirable side effects of the Cubic Spline interpolation that may occur at the beginning and the end of the analyzed portion of the ECG signal ECG(n) which may alter the morphology of T wave;
  • step 34 the portion of the ECG signal ECG(n) with the baseline removed is now subjected to beat filtering to remove non-normal beats from the analysis by quantification system 12 . No interpolation of the missing beats is necessary since continuity of the measurement on a beat to beat basis is not needed.
  • step 36 repolarization segments for each of the beats portion of the ECG signal ECG(n) are extracted by quantification system 12 after the non-normal beats are removed in step 34 . More specifically, in one embodiment ST k (n) segments are extracted from the ECG signal ECG(n) based on two-level criterion: the first step identifies the end of the repolarization wave where the slope (first derivative) of the terminal portion of the repolarization equal or inferior to 0.1 mV. The scanning process is going forward relatively to time.
  • ST k (n) ECG(m) for: m ⁇ [R(k)+D, ⁇ (R(k+1) ⁇ R(k))] (1)
  • D is equal to sixteen in order to have a repolarization segment beginning 100 msec after the R peak of the complex k.
  • is equal to 4/7 for normal range of heart rate (50-100 bpm), but can be manually adjusted for heart rate outside of this range.
  • the repolarization interval is defined either by a flat segment at the end of the repolarization segment or as a portion of the RR interval of the currently analyzed beat. All ST k (n) intervals beginning D/SF (sixteen samples at 200 Hz for 80 msec) after the R peak.
  • ST k (n) segments are extracted from the ECG signal ECG(n) from a starting time S 0 for each segment until an endpoint of the T wave is reached at T o for each segment.
  • the starting time S 0 for each segment is established as a first period of time, such as 100 msec after the R peak, although this time can vary as needed.
  • the particular point for starting time S 0 can be modified in order to avoid including QRS complex components in the area based measurements when a patient has electrical abnormalities leading to an increased QRS complex duration, such as patient with bundle branch block.
  • the starting time S 0 is set by an operator after analyzing the ECG signal for a particular subject, although the starting time So could be established in other manners, such as by an analysis of a characteristic or characteristics of the ECG signal by the quantification system 12 .
  • the endpoint of the T wave T o for each segment is established when the slope of the T wave goes below a threshold slope which is set in the quantification system 12 , although the endpoint of the T wave can be established in other manners. This threshold slope can be changed by the operator.
  • the slope of the T wave portion of the segment is taken by the quantification system 12 by taking a derivative of the T wave portion of that segment, although the quantification system can determine the slope in other manners.
  • the index r identifies the set of NB beats, with the values of r ranging from three to NB-2.
  • the initial portion of the ECG(n) can include seventy beats with two beats removed at each of the extremities of the portion of the ECG signal resulting in sixty-six beats analyzed with each analyzed set being eleven beats, although the size of the analyzed set can vary as needed.
  • an area based repolarization interval for each of the extracted repolarization segments ST k (n) for the portion of the ECG signal ECG(n) are obtained or measured by the quantification system 12 .
  • the area is obtained by the mathematical integration of each of the extracted repolarization segments ST k (n) by the quantification system 12 from the starting time So to various percentages of the ST segment ST k (n). The time needed to reach different percentages of the total area from the starting time S 0 (shown in FIG. 3B ) is measured and recorded.
  • the times needed to reach 25%, 50%, 90% and 97% of the area based repolarization interval for each of the extracted repolarization segment ST k (n) are recorded, although other percentages can be used. If median signal ST r (n) are computed, then the area based repolarization intervals are based on median signal ST r (n), instead of on the individual segments ST k (n). By way of example, six area based repolarization intervals are superimposed as shown in FIG. 6A and their cumulative repolarization curves are shown in FIG. 6B .
  • step 40 the quantification system 12 determines if there are any more ECG signals to process. If there are more ECG signals to process, then the Yes branch is taken back to step 30 . If there are no more ECG signals to process, then the No branch is taken back to step 42 .
  • the quantified area based repolarization interval or intervals are calculated.
  • the calculated interval or intervals may be output or displayed in a variety of different formats, such as on an hourly basis, a day/night basis, or on a twenty-four hour basis.
  • the quantified area based repolarization interval or intervals may also be analyzed by quantification system 12 to detect any alterations in ventricular repolarization.
  • the morphology of the ventricular repolarization segment is analyzed by looking at the length of the interval needed to reach a certain percentage of A(n). For example, if the time need to reach a certain percentage deviates from a designated standard by more than a standard deviation, then the ventricular repolarization is designated as altered.
  • the quantified area based repolarization intervals allow for the identification of slow trends in the area based repolarization interval. For example, in the study of pharmacological agents, the morphology of the area based repolarization interval before, during and/or after exposure to the pharmacological agent can be examined to identify trends and problems.
  • Quantification system 12 may perform other types of analysis, such as repolarization duration vs. cycle length relationship—paired RTx % and RR intervals are determined, the RTx %/RR) slope and its 95% confidence interval are obtained and the scattergram of the raw RTx % vs. RR relationship is graphically displayed.
  • the slope of the regression line is a measure of QT dynamicity.
  • Quantification system 12 may compare RTx % by RR interval bins to determine if there has been a change in the RTx % interval with an intervention that also may have a small effect on the heart rate, the RTx % interval is grouped into selected RR interval ranges, thus avoiding the need to correct for heart rate.
  • RTx %b refers to the RT interval during the baseline ECG recording
  • RTx %d refers to the RT value after drug administration.
  • a statistical comparison of the difference between the RTx %b and RTx %d for values in the RR interval bin, such as the 601 to 700 msec RR interval range is represented by the P-value . Similar comparisons are provided for other RR interval bins, such as 701-800 msec and 801-900 msec as shown in FIG. 7 , with RTx % values represented by X 1 to X 4 .
  • RR interval bins can be arbitrarily selected by the operator of the quantification system 12 .
  • step 44 the output form the analysis of one patient with the present invention is illustrated. Referring back to FIG. 4 , this embodiment of the method ends in step 44 .
  • the present invention described above uses the principle of signal integration of the ventricular repolarization segment to accurately identify ventricular repolarization changes induced by physiologic conditions, such as exercise, cardiac disorders, such as heart disease, or pharmacological agents, such as the evaluation of the effect of a drug on the ventricular repolarization of the heart.
  • the present invention can be applied to ECG signals acquired in variety of dynamic settings including: short- and long-term Holter monitoring, exercise ECG testing, bedside telemetry ECG monitoring, ECG recording using event recorders, implantable devices, or intracardiac electrodes.

Abstract

A method for analyzing an ECG signal includes obtaining a measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal and detecting an altered ventricular repolarization based on the obtained measured area based repolarization interval.

Description

  • This application claims the benefit of U.S. Provisional Patent Application Ser. No. 60/311,921 filed Aug. 13, 2001 which is hereby incorporated by reference in its entirety.
  • FIELD OF THE INVENTION
  • The present invention is directed to a method and system for the quantification of repolarization changes when measured in a dynamic electrocardiogram (“ECG”) signal.
  • BACKGROUND OF THE INVENTION
  • An ECG signal represents changes in electrical potential produced by contractions of the heart recorded from the surface of the body. An example of an ECG signal is illustrated in FIG. 1. Each P-QRS-T complex reflects the electrical depolarization and repolarization components of the heart beat. The cardiac beat comprises three major waves identified as: the P wave representing the depolarization of the auricles, the QRS complex generated by the depolarization of both right and left ventricles, and the T-wave which represents the repolarization of the ventricles. The T wave also includes an apex, Tapex, and an endpoint, Toffset. The repolarization of the auricles is hidden in the QRS complex. The RR intervals are defined by the interval between two consecutive R peaks in successive beats. The RR interval is a direct measurement of the heart rate.
  • In an ECG signal, a modification of the T wave morphology can be observed when certain pharmaceutical agents are taken and also with certain types of cardiac disorders. In particular, it is known that prolongation of the QT interval in an ECG signal is clearly associated with an increased risk for ventricular arrhythmias and sudden cardiac death. Accordingly, prior systems have been developed to monitor the QT interval in ECG signals to detect potential heart problems.
  • Although these prior systems work, they have problems in accurately identifying the endpoint of the T wave when the endpoint of the T wave gradually approaches the baseline, when a U wave is present, or when the shape of T wave is biphasic. If the apex or the endpoint of the T wave are not accurately identified, the accuracy and robustness of repolarization analysis, i.e. analysis of the QT interval, especially with dynamic ECG signals, such as those obtained by exercise ECG testing, is compromised.
  • SUMMARY OF THE INVENTION
  • A method and a computer readable medium with programmed instructions for analyzing an ECG signal in accordance with an embodiment of the present invention includes obtaining a measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal and detecting an altered ventricular repolarization based on the obtained measured area based repolarization interval.
  • A system for analyzing an ECG signal in accordance with an embodiment of the present invention includes a measurement system and a detection system. The measurement system obtains a measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal. The detection system detects an altered ventricular repolarization based on the obtained measured area based repolarization interval.
  • A method and a computer readable medium with programmed instructions for analyzing an effect of a pharmacological agent on an electrocardiogram signal in accordance with an embodiment of the present invention includes obtaining a first measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal. The pharmacological agent is administered. A second measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal is obtained from at least one of a first period during which the pharmacological agent is in effect and a second period after the pharmacological agent is no longer in effect. An altered ventricular repolarization is detected based on the first and second measurements.
  • A system for analyzing an effect of a pharmacological agent on an electrocardiogram signal in accordance with an embodiment of the present invention includes a measurement system and a detection system. The measurement system obtains a first measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal before administering the pharmacological agent. The measurement system also obtains a second measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal from at least one of a first period during which the pharmacological agent is in effect and a second period after the pharmacological agent is no longer in effect. The detection system detects an altered ventricular repolarization based on the first and second measurements.
  • The present invention provides an effective non-invasive method for dynamic quantification of ventricular repolarization. The present invention's comprehensive evaluation of ventricular repolarization provides enhanced utility: in the clinical diagnosis of acquired and inherited repolarization disorders; in the evaluation of new chemical entities for pharmaceutical companies during drug development and during post-marketing surveillance of approved drugs that can adversely affect ventricular repolarization; and in the clinical identification of patients with potentially life-threatening ventricular repolarization disorders, such as those after myocardial infarction.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • FIG. 1 is a diagram of a beat in an ECG signal where two vertical lines mark the apex and the offset of the T wave used to measure QT intervals;
  • FIG. 2 is a block diagram of a quantification system in accordance with one embodiment of the present invention;
  • FIG. 3A is a diagram of another ECG signal;
  • FIG. 3B is a graph of a normalized absolute integral of the repolarization segment of the ECG signal shown in FIG. 3A;
  • FIG. 4 is a flow chart of a method for analyzing an ECG signal in accordance with another embodiment of the present invention;
  • FIGS. 5A-5C are graphs of the magnitude, phase, and zeros and poles of a filter for the quantification system;
  • FIG. 6A is a graph of six area based repolarization intervals of median beats which are superimposed;
  • FIG. 6B is a graph of cumulative area of the six area based repolarization intervals shown in FIG. 6A;
  • FIG. 7 is a table of RTx % and RR interval bins;
  • FIG. 8A is a graph of a heart rate of one patient showing beats per minute (BPM) against time;
  • FIG. 8B is a graph of a RT intervals (Apex and offset of T waves using median value across leads) of the one patient;
  • FIGS. 9A-9D are graphs of trends for T-wave amplitude for the median across leads, the X leads, the Y leads, and the Z leads;
  • FIGS. 10A-10D are graphs of the trends for area-based measurements for 25% of the area for the median, , the X leads, the Y leads, and the Z leads;
  • FIGS. 11A-11D are graphs of the trends for area-based measurements for 50% of the area for the median, , the X leads, the Y leads, and the Z leads;
  • FIGS. 12A-12D are graphs of the trends for area-based measurements for 90% of the area for the median, , the X leads, the Y leads, and the Z leads;
  • FIGS. 13A-13D are graphs of the trends for area-based measurements for 97% of the area for the median,, the X leads, the Y leads, and the Z leads;
  • FIG. 14 is a table of averages for 24 hours, during the day and at night for: beats per minute (BPM); X, Y, and Z amplitudes; X RT apex, Y RT apex, and Z RT apex; X RT offset, Y RT offset, and Z RT offset; XT 25% area, YT 25% area and ZT 25% area; XT 50% area, YT 50% area and ZT 50%; area XT 90% area, YT 90% area and ZT 90%; and area XT 97% area, YT 97% area and ZT 97% area;
  • FIG. 15 is an hourly table of heart rate (BPM), T wave ampl X, T wave ampl Y, and T wave ampl Z;
  • FIG. 16 is an hourly table of heart rate (BPM), T wave apex X, T wave apex Y, and T wave apex Z;
  • FIG. 17 is an hourly table of heart rate (BPM), T wave offset lead X, T wave offset lead Y, and T wave offset lead Z;
  • FIG. 18 is an hourly table of heart rate (BPM), ABRI RT25 X-lead, ABRI RT25 Y-lead, and ABRI RT25 Z-lead;
  • FIG. 19 is an hourly table of heart rate (BPM), ABRI RT50 X-lead, ABRI RT50 Y-lead, and ABRI RT50 Z-lead;
  • FIG. 20 is an hourly table of heart rate (BPM), ABRI RT90 X-lead, ABRI RT90 Y-lead, and ABRI RT90 Z-lead;
  • FIG. 21 is an hourly table of heart rate (BPM), ABRI RT97X-lead, ABRI RT97 Y-lead, and ABRI RT97 Z-lead;
  • FIGS. 22A-22B are graphs of PointCarre plots over 24-hours of RT offset=f(RR) (median across leads) and RT apex=f(RR) (median across leads);
  • FIGS. 23A-23B are graphs of PointCarre plots over 24-hours of RT area 25%=f(RR) (median across leads) and RT area 50%=f(RR) (median across leads); and
  • FIGS. 24A-24B are graphs of PointCarre plots over 24-hours of RT area 90%=f(RR) (median across leads) and RT area 97%=f(RR) (median across leads).
  • DETAILED DESCRIPTION
  • A system 10 and method for analyzing an ECG signal in accordance with one embodiment of the present invention are illustrated in FIGS. 2 and 4. In this particular embodiment, the system 10 includes a quantification system 12 coupled to an ECG device 24 and a storage device 26. The method includes filtering and then extracting a portion of the filtered ECG signal, substantially removing any baseline from the extracted portion of the ECG signal, obtaining a measurement of at least one area based repolarization interval in the extracted portion, and then detecting any altered ventricular repolarization based on the measured area based repolarization interval. The present invention provides a number of advantages including providing an effective non-invasive method for dynamic quantification of ventricular repolarization.
  • Referring to FIG. 2, the quantification system 12 has a central processing unit (“CPU”) or processor 14, a memory 16, an input/output interface 22, a display 18, and a user input device 20 which are coupled together by a bus system 23 or other link, although the quantification system 12 may comprise other types of components, other numbers of the components, and other combinations of the components. The processor 14 executes a program of stored instructions for at least a portion of the method for analyzing an ECG signal in accordance with one embodiment of the present invention as described herein. In this particular embodiment, the preprocessing and processing steps are carried out by quantification system 12, although other types and combinations of systems and/or devices can be used to execute the preprocessing and processing steps.
  • Memory 16 comprises a random access memory (RAM) and a read only memory (ROM), although other types and combinations of memory storage devices can be used, such as a floppy disk, hard disk, CD ROM, or other computer readable medium which is read from and/or written to by a magnetic, optical, or other reading and/or writing system that is coupled to the processor 14. Although in his particular embodiment, the method in accordance with one embodiment of the present invention is stored as programmed instructions in memory 16 in the quantification system 12 for execution by the processor 14, some or all of the programmed instructions could be stored elsewhere.
  • The display 18 enables an operator to observe information, such as the ECG reading for a patient. In this particular embodiment, the display 18 is a cathode ray tube device, although other types of displays can be used, such as a printer.
  • The user input device 20 enables an operator to generate and transmit signals or commands to the processor, such as a request to obtain or retrieve ECG signals for a particular patient for processing. In this particular embodiment, the user input device 20 is a keyboard, although other types of user input devices can be used, such as a computer mouse.
  • The input/output interface 22 is used to operatively couple the quantification system 12 to the ECG device 24 and the storage device 26, although quantification system 12 can be coupled via input/output interface to other systems and/or devices.
  • The ECG device 24 records one or more ECG signals from one or more patients. The ECG device can capture the electrical activity from the heart in a variety of different manners, such as Holter recordings, exercise ECG testing, bedside ECG monitoring, event monitoring, implantable ECG recorders using one or more leads and with a duration that may vary from few minutes to twenty-four hours, even few days. A variety of different types of ECG devices can be used for ECG device 24. The ECG information can be transmitted to the quantification system 12 from the ECG device 24 for processing.
  • The storage device 26 comprises a RAM, although other types and combinations of memory storage devices can be used, such as a ROM or a floppy disk, hard disk, CD ROM, or other computer readable medium which is read from and/or written to by a magnetic, optical, or other reading and/or writing system. The storage device 26 can store ECG readings for processing by the quantification system 12 and can store the results of any processing.
  • Referring to FIG. 3A, an example of another ECG signal is illustrated. In this particular example, the Q peak is reached at time Q0, the T wave begins at time S0, the T apex is reached at Tm, and the endpoint of the T wave is reached at To.
  • Referring to FIG. 3B, a graph of the measured area under the repolarization interval shown in FIG. 3A with the Y axis showing the percentage of the total area under the repolarization interval and the X-axis showing time is illustrated. The graph begins at time S0 which again corresponds to the time when the T wave begins and U0 represents the time when substantially 100% of the area under the repolarization area is measured. The times to accumulate 25%, 50% and 75% of the absolute T-wave area are identified as LA25, LA50, LA75, respectively. Unlike prior systems, the present invention does not rely upon accurately determining either the T apex, Tm, or the endpoint of the T wave, To. Instead, the present system identify the end of the T wave when the slope of the repolarization signal is below a given threshold. The present invention quantifies changes in ventricular repolarization based on measurements of the area under the repolarization interval obtained starting at a time S0, which is a first period of time after the R peak is detected.
  • The morphology of the repolarization T wave in an ECG signal, such as the one shown in FIG. 3A, can be influenced by any factor that alters ion-channel kinetics involved in ventricular repolarization, including inherited diseases, acquired cardiac diseases, metabolic disturbances, and drugs. Since the normal morphology of the ventricular repolarization interval is well known to those of ordinary skill in the art, an altered or abnormal morphology ventricular repolarization interval can be both visually and quantitatively assessed if it is accurately measured. A process for analyzing ventricular repolarization in accordance with the present invention is described below.
  • Referring to FIGS. 2, 3A, 3B, and 4, a method for analyzing an ECG signal in accordance with one embodiment is described below. Beginning in step 28, the quantification system 12 obtains an ECG signal or signals. The ECG signals can be obtained from a variety of different sources, such as directly from the ECG device 24 or from the storage device 26 which stores previously recorded ECG signals.
  • In step 30, the ECG signal undergoes low pass-filtering in quantification system 12. The filtering technique used should preserve the initial shape of the T wave. In this particular embodiment, the ECG signal is filtered in the forward direction, the filtered ECG signal is then reversed and run back through the filter again. The result of this filtering operation has precisely zero phase distortion and magnitude modified by the square of the filter's magnitude response.
  • By way of example only, in this particular embodiment a Butterworth filter and a bidirectional filtering technique are used in order to obtain a high quality filtering free of distorted components, although other types of filters and filtering techniques can be used. In this example, a filter of order two with a cutoff frequency at −3 db equal to 18 Hz is used. Graphs of the magnitude, phase (radians) and zeros and poles of this filter used for removing noise components of the ECG signal and reducing morphological changes of area based repolarization interval are illustrated in FIGS. 5A-5C.
  • Referring back to FIG. 4, in step 32 a portion or limited length sequence of the filtered ECG signal is extracted by the quantification system 12. More specifically, the ECG signal is cut into portions or sequences of consecutive continuous beats by quantification system 12, excluding non-normal beats. In each of these portions, R(k) is the R peak location in the kth beat in each portion and ECG(n) is the nth portion of ECG signal where n represents the number of portions of the ECG signal. Each ECG portion which is analyzed is ECG(n) for n ε[R(k), R(k+NB)] where NB is the number of beats in each analyzed interval. In this example, a repolarization segment is extracted for each of seventy consecutive beats in the filtered ECG signal.
  • In step 33, the baseline B(k) of the extracted portion of the ECG signal ECG(n) is substantially removed by the quantification system 12. The baseline removal must be accurate and must avoid adding significant components to the repolarization segment in the ECG signal ECG(n). By way of example only, in this particular embodiment the baseline removal can be obtained based using a Cubic-Spline interpolation of the baseline of the ECG signal ECG(n), although other techniques for removing the baseline can be used. The interpolation used in this particular example is based on the average amplitude value of the samples from a limited squared time window W(n) of 80 ms length located just before the beginning of the QRS complex. B(k) is defined such as: B(k)=ECG(R(k)−16) for a signal at 200 Hz sampling frequency. The signal B(k) is interpolated at a sampling period equal to the sampling period of the ECG. The Cubic Spline interpolation of B(k) is used to insure continuous smoothness up to the second derivative and is less likely to generate erratic oscillations than polynomial interpolation. The interpolated curve, corresponding to the very low components of the ECG signal ECG(n), is subtracted from the ECG signal ECG(n) in order to reduce the effect of respiratory components of the signal. After subtracting the estimated baseline, two beats are removed from the beginning and two beats from the end of the analyzed intervals in the ECG signal ECG(n). In this particular example, two beats are removed from the extremities of the ECG signal with seventy beats resulting in sixty-six beats which are analyzed. This insures elimination of undesirable side effects of the Cubic Spline interpolation that may occur at the beginning and the end of the analyzed portion of the ECG signal ECG(n) which may alter the morphology of T wave;
  • In step 34, the portion of the ECG signal ECG(n) with the baseline removed is now subjected to beat filtering to remove non-normal beats from the analysis by quantification system 12. No interpolation of the missing beats is necessary since continuity of the measurement on a beat to beat basis is not needed.
  • In step 36, repolarization segments for each of the beats portion of the ECG signal ECG(n) are extracted by quantification system 12 after the non-normal beats are removed in step 34. More specifically, in one embodiment STk(n) segments are extracted from the ECG signal ECG(n) based on two-level criterion: the first step identifies the end of the repolarization wave where the slope (first derivative) of the terminal portion of the repolarization equal or inferior to 0.1 mV. The scanning process is going forward relatively to time. If such criterion is not met because the signal does has such characteristic then: STk(n)=ECG(m) for:
    mε[R(k)+D, β(R(k+1)−R(k))]  (1)
    By way of example only, at a 200 Hz sampling frequency (SF), D is equal to sixteen in order to have a repolarization segment beginning 100 msec after the R peak of the complex k. β is equal to 4/7 for normal range of heart rate (50-100 bpm), but can be manually adjusted for heart rate outside of this range. As previously described, the repolarization interval is defined either by a flat segment at the end of the repolarization segment or as a portion of the RR interval of the currently analyzed beat. All STk(n) intervals beginning D/SF (sixteen samples at 200 Hz for 80 msec) after the R peak.
  • In another embodiment, STk(n) segments are extracted from the ECG signal ECG(n) from a starting time S0 for each segment until an endpoint of the T wave is reached at To for each segment. The starting time S0 for each segment is established as a first period of time, such as 100 msec after the R peak, although this time can vary as needed. For example, the particular point for starting time S0 can be modified in order to avoid including QRS complex components in the area based measurements when a patient has electrical abnormalities leading to an increased QRS complex duration, such as patient with bundle branch block. In this embodiment, the starting time S0 is set by an operator after analyzing the ECG signal for a particular subject, although the starting time So could be established in other manners, such as by an analysis of a characteristic or characteristics of the ECG signal by the quantification system 12. In this embodiment, the endpoint of the T wave To for each segment is established when the slope of the T wave goes below a threshold slope which is set in the quantification system 12, although the endpoint of the T wave can be established in other manners. This threshold slope can be changed by the operator. The slope of the T wave portion of the segment is taken by the quantification system 12 by taking a derivative of the T wave portion of that segment, although the quantification system can determine the slope in other manners.
  • In this embodiment, a median signal STr(n) is computed from a set of segments STk(n) such as: STr(n)=median([ST3(n), . . . , STNB-2(n)]). The index r identifies the set of NB beats, with the values of r ranging from three to NB-2. For instance, as described earlier the initial portion of the ECG(n) can include seventy beats with two beats removed at each of the extremities of the portion of the ECG signal resulting in sixty-six beats analyzed with each analyzed set being eleven beats, although the size of the analyzed set can vary as needed.
  • In step 38, in one embodiment an area based repolarization interval for each of the extracted repolarization segments STk(n) for the portion of the ECG signal ECG(n) are obtained or measured by the quantification system 12. The area is obtained by the mathematical integration of each of the extracted repolarization segments STk(n) by the quantification system 12 from the starting time So to various percentages of the ST segment STk(n). The time needed to reach different percentages of the total area from the starting time S0 (shown in FIG. 3B) is measured and recorded. In this particular embodiment, the times needed to reach 25%, 50%, 90% and 97% of the area based repolarization interval for each of the extracted repolarization segment STk(n) are recorded, although other percentages can be used. If median signal STr(n) are computed, then the area based repolarization intervals are based on median signal STr(n), instead of on the individual segments STk(n). By way of example, six area based repolarization intervals are superimposed as shown in FIG. 6A and their cumulative repolarization curves are shown in FIG. 6B.
  • In step 40, the quantification system 12 determines if there are any more ECG signals to process. If there are more ECG signals to process, then the Yes branch is taken back to step 30. If there are no more ECG signals to process, then the No branch is taken back to step 42.
  • In step 42, the quantified area based repolarization interval or intervals are calculated. The calculated interval or intervals may be output or displayed in a variety of different formats, such as on an hourly basis, a day/night basis, or on a twenty-four hour basis.
  • The quantified area based repolarization interval or intervals may also be analyzed by quantification system 12 to detect any alterations in ventricular repolarization. In this particular embodiment, the morphology of the ventricular repolarization segment is analyzed by looking at the length of the interval needed to reach a certain percentage of A(n). For example, if the time need to reach a certain percentage deviates from a designated standard by more than a standard deviation, then the ventricular repolarization is designated as altered.
  • The quantified area based repolarization intervals allow for the identification of slow trends in the area based repolarization interval. For example, in the study of pharmacological agents, the morphology of the area based repolarization interval before, during and/or after exposure to the pharmacological agent can be examined to identify trends and problems.
  • Quantification system 12 may perform other types of analysis, such as repolarization duration vs. cycle length relationship—paired RTx % and RR intervals are determined, the RTx %/RR) slope and its 95% confidence interval are obtained and the scattergram of the raw RTx % vs. RR relationship is graphically displayed. The slope of the regression line is a measure of QT dynamicity.
  • Quantification system 12 may compare RTx % by RR interval bins to determine if there has been a change in the RTx % interval with an intervention that also may have a small effect on the heart rate, the RTx % interval is grouped into selected RR interval ranges, thus avoiding the need to correct for heart rate. One example of this approach is shown in the table in FIG. 7. In this table in FIG. 7, RTx %b refers to the RT interval during the baseline ECG recording and RTx %d refers to the RT value after drug administration. A statistical comparison of the difference between the RTx %b and RTx %d for values in the RR interval bin, such as the 601 to 700 msec RR interval range is represented by the P-value . Similar comparisons are provided for other RR interval bins, such as 701-800 msec and 801-900 msec as shown in FIG. 7, with RTx % values represented by X1 to X4. RR interval bins can be arbitrarily selected by the operator of the quantification system 12.
  • Referring to FIGS. 8A-24B, the output form the analysis of one patient with the present invention is illustrated. Referring back to FIG. 4, this embodiment of the method ends in step 44.
  • The present invention described above uses the principle of signal integration of the ventricular repolarization segment to accurately identify ventricular repolarization changes induced by physiologic conditions, such as exercise, cardiac disorders, such as heart disease, or pharmacological agents, such as the evaluation of the effect of a drug on the ventricular repolarization of the heart. The present invention can be applied to ECG signals acquired in variety of dynamic settings including: short- and long-term Holter monitoring, exercise ECG testing, bedside telemetry ECG monitoring, ECG recording using event recorders, implantable devices, or intracardiac electrodes.
  • Having thus described the basic concept of the invention, it will be rather apparent to those skilled in the art that the foregoing detailed disclosure is intended to be presented by way of example only, and is not limiting. Various alterations, improvements, and modifications will occur and are intended to those skilled in the art, though not expressly stated herein. These alterations, improvements, and modifications are intended to be suggested hereby, and are within the spirit and scope of the invention. Additionally, the recited order of processing elements or sequences, or the use of numbers, letters, or other designations therefor, is not intended to limit the claimed processes to any order except as may be specified in the claims. Accordingly, the invention is limited only by the following claims and equivalents thereto.

Claims (46)

1. A method for analyzing an electrocardiogram signal, the method comprising:
obtaining a measurement of an area based repoloarization interval from at least one beat in the electrocardiogram signal; and
detecting an altered ventricular repolarization based on the obtained measured area based repolarization interval.
2. The method as set forth in claim 1 further comprising preprocessing the electrocardiogram signal, the preprocessing including at least one of:
filtering the electrocardiogram signal;
extracting a portion of the electrocardiogram signal with at least the beat; and
substantially removing any baseline from the portion electrocardiogram signal.
3. The method as set forth in claim 2 wherein the preprocessing comprises the filtering, the extracting and the substantially removing.
4. The method as set forth in claim 2 wherein the filtering further comprises:
filtering the electrocardiogram signal with a low-pass filter;
reversing the filtered electrocardiogram signal; and
filtering the reversed filtered electrocardiogram signal with the low-pass filter.
5. The method as set forth in claim 2 wherein substantially removing any baseline comprises:
estimating the baseline using a cubic spline interpolation technique; and
substantially removing the estimated baseline.
6. The method as set forth in claim 1 wherein the obtaining further comprises:
determining a starting point of a T wave in the beat in the electrocardiogram signal; and
determining an end point of the T wave in the beat in the electrocardiogram signal, wherein the portion of the electrocardiogram signal which is measured is from the starting point to the end point of the T wave.
7. The method as set forth in claim 6 wherein the starting point of the T wave is a first time after an R peak in the beat in the electrocardiogram signal.
8. The method as set forth in claim 6 wherein the end point of the T wave is when a slope of the T wave in the beat in the electrocardiogram signal goes below a first threshold.
9. The method as set forth in claim 1 wherein the obtaining further comprises:
obtaining a plurality of measurements, wherein each of the measurements is of a different beat in the electrocardiogram signal;
dividing the obtained plurality of measurements into one or more sets; and
determining a median measurement for each of the sets, wherein the detecting an altered ventricular repolarization is based on the median measurement for one or more of the sets.
10. The method as set forth in claim 9 wherein the detecting an altered ventricular repolarization further comprises:
determining a time need to reach at least a first set percentage less than one hundred percent of the median measurement for at least one of the sets; and
comparing the determined time need to reach the first set percentage against another time to reach the first set percentage, wherein an altered ventricular repolarization is detected if the determined time deviates from the another time by more than a first set amount.
11. The method according to claim 10 wherein the another time comprises at least one of a standard time and another determined time to reach the first set percentage.
12. The method according to claim 9 further comprising removing one or more of the beats which are not normal from the electrocardiogram signal.
13. The method as set forth in claim 9 further comprising removing one or more of the beats from at least one end of the electrocardiogram signal.
14. The method as set forth in claim 1 wherein the detecting an altered ventricular repolarization further comprises:
determining a time need to reach at least a first set percentage less than one hundred percent of the obtained measured area based repolarization interval; and
comparing the determined time need to reach the first set percentage against another time to reach the first set percentage, wherein an altered ventricular repolarization is detected if the determined time deviates from the standard time by more than a first set amount.
15. The method according to claim 14 wherein the another time comprises at least one of a standard time and another determined time to reach the first set percentage.
16. A computer readable medium having stored thereon instructions for analyzing an electrocardiogram signal which when executed by a processor, causes the processor to perform the steps of:
obtaining a measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal; and
detecting an altered ventricular repolarization based on the obtained measured area based repolarization interval.
17. The medium as set forth in claim 15 further comprising preprocessing the electrocardiogram signal, the preprocessing including at least one of:
filtering the electrocardiogram signal;
extracting a portion of the electrocardiogram signal with at least the beat; and
substantially removing any baseline from the portion electrocardiogram signal.
18. The medium as set forth in claim 17 wherein the preprocessing comprises the filtering, the extracting and the substantially removing.
19. The medium as set forth in claim 17 wherein the filtering further comprises:
filtering the electrocardiogram signal with a low-pass filter;
reversing the filtered electrocardiogram signal; and
filtering the reversed filtered electrocardiogram signal with the low-pass filter.
20. The medium as set forth in claim 17 wherein substantially removing any baseline comprises:
estimating the baseline using a cubic spline interpolation technique; and
substantially removing the estimated baseline.
21. The medium as set forth in claim 16 wherein the obtaining further comprises:
determining a starting point of a T wave in the beat in the electrocardiogram signal; and
determining an end point of the T wave in the beat in the electrocardiogram signal, wherein the portion of the electrocardiogram signal which is measured is from the starting point to the end point of the T wave.
22. The medium as set forth in claim 21 wherein the starting point of the T wave is a first time after an R peak in the beat in the electrocardiogram signal.
23. The medium as set forth in claim 21 wherein the end point of the T wave is when a slope of the T wave in the beat in the electrocardiogram signal goes below a first threshold.
24. The medium as set forth in claim 16 wherein the obtaining further comprises:
obtaining a plurality of measurements, wherein each of the measurements is of a different beat in the electrocardiogram signal;
dividing the obtained plurality of measurements into one or more sets; and
determining a median measurement for each of the sets, wherein the detecting an altered ventricular repolarization is based on the median measurement for one or more of the sets.
25. The medium as set forth in claim 24 wherein the detecting an altered ventricular repolarization further comprises:
determining a time need to reach at least a first set percentage less than one hundred percent of the median measurement for at least one of the sets; and
comparing the determined time need to reach the first set percentage against another time to reach the first set percentage, wherein an altered ventricular repolarization is detected if the determined time deviates from the standard time by more than a first set amount.
26. The medium according to claim 25 wherein the another time comprises at least one of a standard time and another determined time to reach the first set percentage.
27. The medium according to claim 24 further comprising removing one or more of the beats which are not normal from the electrocardiogram signal.
28. The medium as set forth in claim 24 further comprising removing one or more of the beasts from at least one end of the electrocardiogram signal.
29. The medium as set forth in claim 16 wherein the detecting an altered ventricular repolarization further comprises:
determining a time need to reach at least a first set percentage less than one hundred percent of the obtained measured area based repolarization interval; and
comparing the determined time need to reach the first set percentage against another time to reach the first set percentage, wherein an altered ventricular repolarization is detected if the determined time deviates from the standard time by more than a first set amount.
30. The medium according to claim 29 wherein another time comprises at least one of a standard time and another determined time to reach the first set percentage.
31. A system for analyzing an electrocardiogram signal, the system comprising:
a measurement system that obtains a measurement of an area based repolarization interval from at least one beat in the electrocardiogram signal; and
a detecting system that detects an altered ventricular repolarization based on the obtained measured area based repolarization interval.
32. The system as set forth in claim 31 further comprising a preprocessing system that includes at least one of:
a filtered system that filters the electrocardiogram signal;
an extraction system that extracts a portion of the electrocardiogram signal with at least the beat; and
a baseline removal system that substantially removes any baseline from the portion electrocardiogram signal.
33. The system as set forth in claim 32 wherein the preprocessing system comprises the filter system, the extraction system and the baseline removal system.
34. The system as set forth in claim 32 wherein the filter system further comprises a bidirectional, low pass filter that filters the electrocardiogram signal.
35. The system as set forth in claim 32 wherein baseline removal system further comprises an estimation system that estimates the baseline using a cubic spline interpolation technique, the baseline removal system substantially removing the estimated baseline.
36. The system as set forth in claim 32 wherein the measurement system further comprises a determining system that determines a starting point of a T wave in the beat in the electrocardiogram signal and an end point of the T wave in the beat in the electrocardiogram signal, wherein the portion of the electrocardiogram signal which is measured is from the starting point to the end point of the T wave.
37. The system as set forth in claim 36 wherein the starting point of T wave is a first time after an R peak in the beat in the electrocardiogram signal.
38. The system as set forth in claim 36 wherein the end point of the T wave is when a slope of the T wave in the beat in the electrocardiogram signal goes below a first threshold.
39. The system as set forth in claim 31 wherein the measurement system further comprises a median measurement determining system that obtains a plurality of measurements, wherein each of the measurements is of a different beat in the electrocardiogram signal, divides the obtained plurality of measurements into one or more sets, and determines a median measurement for each of the sets, wherein the detecting an altered ventricular repolarization is based on the median measurement for at least one of the sets.
40. The system as set forth in claim 39 wherein the detection system further comprises:
a determining system that determines a time need to reach at least a first set percentage less than one hundred percent of the median measurement for one or more of the sets; and
a comparison system that compares the determined time need to reach the first set percentage against another time to reach the first set percentage, wherein an altered ventricular repolarization is detected by the detection system if the determined tine deviates from the another time by more than a first set amount.
41. The system according to claim 40 wherein the another time comprises at least one of a standard time and another determined time to reach the first set percentage.
42. The system according to claim 39 further comprising a removal system that removes one or more of the beats which are not normal from the electrocardiogram signal.
43. The system sa set forth in claim 39 further comprising a removing system that removes one or more of the beats from at least one end of the electrocardiogram signal.
44. The system as set forth in claim 31 wherein the detection system further comprises:
a determining system that determines a time need to reach at least a first set percentage less than one hundred percent of the obtained measured area based repolorization interval; and
a comparison system that compares the determined time need to reach the first set percentage against another time to reach the first set percentage, wherein an altered ventricular repolarization is detected by the detection system if the determined time deviates from the another time by more than a first set amount.
45. The system according to claim 44 wherein the another time comprises at least one of a standard time and another determined time to reach the first set percentage.
46-59. (canceled)
US11/407,902 2001-08-13 2006-04-21 Method and system for analyzing an electrocardiographic signal Abandoned US20060189876A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US11/407,902 US20060189876A1 (en) 2001-08-13 2006-04-21 Method and system for analyzing an electrocardiographic signal

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US31192101P 2001-08-13 2001-08-13
US10/217,883 US7463921B2 (en) 2001-08-13 2002-08-13 Method and system for analyzing an electrocardiographic signal
US11/407,902 US20060189876A1 (en) 2001-08-13 2006-04-21 Method and system for analyzing an electrocardiographic signal

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US10/217,883 Division US7463921B2 (en) 2001-08-13 2002-08-13 Method and system for analyzing an electrocardiographic signal

Publications (1)

Publication Number Publication Date
US20060189876A1 true US20060189876A1 (en) 2006-08-24

Family

ID=23209073

Family Applications (2)

Application Number Title Priority Date Filing Date
US10/217,883 Active 2024-06-16 US7463921B2 (en) 2001-08-13 2002-08-13 Method and system for analyzing an electrocardiographic signal
US11/407,902 Abandoned US20060189876A1 (en) 2001-08-13 2006-04-21 Method and system for analyzing an electrocardiographic signal

Family Applications Before (1)

Application Number Title Priority Date Filing Date
US10/217,883 Active 2024-06-16 US7463921B2 (en) 2001-08-13 2002-08-13 Method and system for analyzing an electrocardiographic signal

Country Status (3)

Country Link
US (2) US7463921B2 (en)
AU (1) AU2002324691A1 (en)
WO (1) WO2003015607A2 (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060264768A1 (en) * 2005-05-13 2006-11-23 Cardiocore Lab, Inc. Method and apparatus for sequenced extraction from electrocardiogramic waveforms
US7676264B1 (en) * 2007-04-13 2010-03-09 Pacesetter, Inc. Systems and methods for use by an implantable medical device for evaluating ventricular dyssynchrony based on T-wave morphology
CN104634906A (en) * 2013-11-07 2015-05-20 苏州普源精电科技有限公司 Baseline adjusting method of resolved peaks, and chromatographic workstation with baseline adjusting function
US20150265175A1 (en) * 2014-03-18 2015-09-24 Jami SHAKIBI GILANI Non-Invasive Evaluation of Cardiac Repolarisation Instability for Risk Stratification of Sudden Cardiac Death
CN109381181A (en) * 2017-08-14 2019-02-26 深圳大学 The end-point detecting method of electrocardiosignal signature waveform
DE102007057885B4 (en) 2006-11-30 2021-12-09 General Electric Co. Method for analyzing physiological data and checking the quality of the measurement using the analysis of the main components

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7412282B2 (en) * 2005-01-26 2008-08-12 Medtronic, Inc. Algorithms for detecting cardiac arrhythmia and methods and apparatuses utilizing the algorithms
US8391964B2 (en) * 2009-05-11 2013-03-05 Medtronic, Inc. Detecting electrical conduction abnormalities in a heart
US9572505B2 (en) 2012-10-11 2017-02-21 Medtronic, Inc. Determining onsets and offsets of cardiac depolarization and repolarization waves
CN104173043B (en) * 2014-09-04 2017-02-15 东莞理工学院 Electrocardiogram (ECG) data analysis method suitable for mobile platform
US10610115B1 (en) 2017-05-24 2020-04-07 Eresearch Technology, Inc. Methods and systems for the delivery of accurate and precise measurements from the body-surface electrocardiogram
CN109645975A (en) * 2017-10-10 2019-04-19 深圳市理邦精密仪器股份有限公司 The processing method and processing device of fetal heart rate curve data
CN108236461B (en) * 2017-12-22 2021-02-19 天津天堰科技股份有限公司 Processing method of electrocardiosignal capable of being edited

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5560370A (en) * 1991-02-20 1996-10-01 Georgetown University Method and apparatus for prediction of cardiac electrical instability by simultaneous assessment of T-wave alternans and QT interval dispersion
US5827195A (en) * 1997-05-09 1998-10-27 Cambridge Heart, Inc. Electrocardiogram noise reduction using multi-dimensional filtering
US20020151806A1 (en) * 2000-06-26 2002-10-17 Starobin Joseph M. Method and system for evaluating cardiac ischemia with an abrupt stop exercise protocol
US20030013978A1 (en) * 2001-07-12 2003-01-16 Schlegel Todd T. Real-time, high frequency QRS electrocardiograph
US6847840B2 (en) * 2000-09-25 2005-01-25 Pfizer, Inc. System and method for statistical analysis of QT interval as a function of changes in RR interval
US6983183B2 (en) * 2001-07-13 2006-01-03 Cardiac Science, Inc. Method and apparatus for monitoring cardiac patients for T-wave alternans

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2669523B1 (en) * 1990-11-23 1997-06-06 Fred Zacouto HEART FAILURE PREVENTION DEVICE

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5560370A (en) * 1991-02-20 1996-10-01 Georgetown University Method and apparatus for prediction of cardiac electrical instability by simultaneous assessment of T-wave alternans and QT interval dispersion
US5827195A (en) * 1997-05-09 1998-10-27 Cambridge Heart, Inc. Electrocardiogram noise reduction using multi-dimensional filtering
US20020151806A1 (en) * 2000-06-26 2002-10-17 Starobin Joseph M. Method and system for evaluating cardiac ischemia with an abrupt stop exercise protocol
US6847840B2 (en) * 2000-09-25 2005-01-25 Pfizer, Inc. System and method for statistical analysis of QT interval as a function of changes in RR interval
US20030013978A1 (en) * 2001-07-12 2003-01-16 Schlegel Todd T. Real-time, high frequency QRS electrocardiograph
US6983183B2 (en) * 2001-07-13 2006-01-03 Cardiac Science, Inc. Method and apparatus for monitoring cardiac patients for T-wave alternans

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060264768A1 (en) * 2005-05-13 2006-11-23 Cardiocore Lab, Inc. Method and apparatus for sequenced extraction from electrocardiogramic waveforms
US8055331B2 (en) * 2005-05-13 2011-11-08 Cardiocore Lab, Inc. Method and apparatus for sequenced extraction from electrocardiogramic waveforms
DE102007057885B4 (en) 2006-11-30 2021-12-09 General Electric Co. Method for analyzing physiological data and checking the quality of the measurement using the analysis of the main components
US7676264B1 (en) * 2007-04-13 2010-03-09 Pacesetter, Inc. Systems and methods for use by an implantable medical device for evaluating ventricular dyssynchrony based on T-wave morphology
CN104634906A (en) * 2013-11-07 2015-05-20 苏州普源精电科技有限公司 Baseline adjusting method of resolved peaks, and chromatographic workstation with baseline adjusting function
US20150265175A1 (en) * 2014-03-18 2015-09-24 Jami SHAKIBI GILANI Non-Invasive Evaluation of Cardiac Repolarisation Instability for Risk Stratification of Sudden Cardiac Death
US20150265174A1 (en) * 2014-03-18 2015-09-24 Jami SHAKIBI GILANI Non-Invasive Evaluation of Cardiac Repolarisation Instability for Risk Stratification of Sudden Cardiac Death
US9474460B2 (en) * 2014-03-18 2016-10-25 Jami SHAKIBI GILANI Non-invasive evaluation of cardiac repolarisation instability for risk stratification of sudden cardiac death
CN109381181A (en) * 2017-08-14 2019-02-26 深圳大学 The end-point detecting method of electrocardiosignal signature waveform

Also Published As

Publication number Publication date
US7463921B2 (en) 2008-12-09
AU2002324691A1 (en) 2003-03-03
WO2003015607A3 (en) 2003-10-16
WO2003015607A2 (en) 2003-02-27
US20030050565A1 (en) 2003-03-13

Similar Documents

Publication Publication Date Title
US20060189876A1 (en) Method and system for analyzing an electrocardiographic signal
US7412282B2 (en) Algorithms for detecting cardiac arrhythmia and methods and apparatuses utilizing the algorithms
US6169919B1 (en) System and method for quantifying alternation in an electrocardiogram signal
JP6381444B2 (en) Method and system for automatically measuring ST level of electrocardiogram in real-time ECG
JP5539199B2 (en) Automatic identification of the responsible coronary artery
US8560054B2 (en) Method and apparatus for extracting optimum holter ECG reading
US7412283B2 (en) Method and system for processing electrocardial signals
JPS6134331B2 (en)
US8437839B2 (en) Electrocardiographic assessment of arrhythmia risk
EP1219236A2 (en) System and method for detecting new left branch bundle block for accelerating treatment of acute myocardial infarction
US6718197B1 (en) LV ectopic density trending
US20070010752A1 (en) System and method for analysing the p-wave of an ecg-signal
US20140330148A1 (en) Cycle length iteration for the detection of atrial activations from electrogram recordings of atrial fibrillation
US6246903B1 (en) Statistical mapping of the physiological state of the heart of a mammal
WO2010000009A1 (en) Improved detection of cardiac dysfunction
Gencel et al. Correlation of signal‐averaged P wave with electrophysiological testing for atrial vulnerability in strokes of unexplained etiology
Laguna et al. Automatic QT interval analysis in postmyocardial infarction patients
WO2020100129A1 (en) A system and a method for non-invasive detection of cad
Schluter et al. A bedside cardiac arrhythmia monitor
JP7245921B2 (en) Apparatus and method for QT correction
RU2791006C1 (en) System and method for automated electrocardiogram analysis and interpretation
US20230414153A1 (en) Methods and Apparatuses for Determining a QT Interval of an ECG Signal
Manis et al. The effect of white noise and false peak Detection on HRV Analysis
Surat-E-Mostafa et al. Estimating spectral heart rate variability (HRV) features with missing RR-interval data
Tirtawinata P-waves Detection in Normal and Pathologic Condition Using SciPy Library with Python Programming Language

Legal Events

Date Code Title Description
STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION