Search Images Maps Play YouTube News Gmail Drive More »
Sign in
Screen reader users: click this link for accessible mode. Accessible mode has the same essential features but works better with your reader.

Patents

  1. Advanced Patent Search
Publication numberUS7289845 B2
Publication typeGrant
Application numberUS 10/734,327
Publication date30 Oct 2007
Filing date12 Dec 2003
Priority date31 Oct 2000
Fee statusPaid
Also published asEP1330289A2, EP2269691A2, EP2269691A3, US6684100, US20040127806, WO2002040094A2, WO2002040094A3
Publication number10734327, 734327, US 7289845 B2, US 7289845B2, US-B2-7289845, US7289845 B2, US7289845B2
InventorsRobert J. Sweeney, Eric G. Lovett
Original AssigneeCardiac Pacemakers, Inc.
Export CitationBiBTeX, EndNote, RefMan
External Links: USPTO, USPTO Assignment, Espacenet
Curvature based method for selecting features from an electrophysiologic signal for purpose of complex identification and classification
US 7289845 B2
Abstract
A method for curvature based complex identification and classification comprises sensing a cardiac signal and computing curvatures at sample points on the sensed cardiac signal. Then to extract features from the computed curvatures, and compare the extracted features with a set of predetermined templates, and then to classify the sensed cardiac signal based on the outcome of the comparison.
Images(5)
Previous page
Next page
Claims(20)
1. A cardiac rhythm management system, comprising:
at least one electrode;
a signal sensing circuit coupled to the electrode to sense a cardiac signal, the cardiac signal having signal complexes;
a controller coupled to the sensing circuit, wherein the controller receives the sensed cardiac signal, and wherein the controller includes:
memory configured to store signal complexes and to store sets of predetermined templates, the predetermined templates including sets of template zones; and
an analyzing module configured to:
compute curvatures at sample points on the sensed cardiac signal,
extract features from the computed curvatures, and
classify the signal complexes by comparing the extracted features to the template zones of the predetermined templates.
2. The system of claim 1, wherein the template zones include a center amplitude, a time width and an amplitude width.
3. The system of claim 2, wherein the template zones include different centers and widths for different types of cardiac signals.
4. The system of claim 3, wherein the controller further includes a scoring module, wherein the scoring module assigns a score to an extracted feature based on how well the feature matches a set of template zones.
5. The system of claim 4, wherein the template zones are modified over time to improve matching with the extracted features.
6. The system of claim 1, wherein the analyzing module is configured to assign an unknown classification to a signal complex when the extracted features fail to match the sets of templates.
7. A method comprising:
sensing a voltage versus time cardiac signal using an implantable medical device;
converting the voltage versus time cardiac signal into a time versus time signal;
computing curvature at sample points on the time versus time signal;
extracting features from the computed curvature;
comparing the extracted features with a set of predetermined templates; and
classifying the cardiac signal based on an outcome of the comparison.
8. The method of claim 7, wherein extracting features from the computed curvatures includes separating a set of features associated with a first cardiac signal complex upon detecting a subsequent second cardiac signal complex.
9. The method of claim 8, wherein separating the set of features includes identifying a fiducial feature from the set of separated features based on comparing the time when the first signal complex occurs with the time when the separated set of features associated with the first signal complex occurs.
10. The method of claim 7, wherein comparing the extracted features with a set of predetermined templates includes comparing the extracted features with a set of predetermined time versus time templates.
11. A cardiac rhythm management system, comprising:
at least one electrode;
a signal sensing circuit coupled to the electrode to sense a cardiac signal, the cardiac signal having signal complexes;
a controller coupled to the sensing circuit, wherein the controller receives the sensed cardiac signal, and wherein the controller includes:
memory configured to store at least a first signal complex; and
an analyzing module configured to:
compute curvatures at sample points of the sensed cardiac signal,
continuously process the sensed cardiac signal to produce an absolute value of the first derivative of the sensed cardiac signal, and
compare the absolute value of the first derivative to a predetermined decaying threshold value to detect a second subsequent signal complex.
12. The system of claim 11, wherein the analyzing module separates a set of features associated with the first signal complex upon detecting the second signal complex.
13. The system of claim 11, wherein the second signal complex is detected when the absolute value of the first derivative of the cardiac signal exceeds the predetermined decaying threshold value.
14. The system of claim 12, wherein the analyzer module is configured to identify a fiducial feature from the set of separated features associated with the first signal complex.
15. The system of claim 14, wherein the fiducial feature is identified based on comparing a time when the first signal complex occurs with the time when the separated features occur.
16. The system of claim 12, wherein the memory stores the features associated with the first signal complex, and the stored features are replaced with features of the second signal complex when the second signal complex is detected.
17. A method comprising:
sensing a cardiac signal using an implantable medical device;
computing curvature at sample points on the sensed cardiac signal;
computing average curvature between two adjacent sample points by computing an integral of the curvature between the two adjacent sample points;
extracting features from the computed curvatures;
comparing the extracted features with a set of predetermined templates; and
classifying the sensed cardiac signal based on an outcome of the comparison.
18. The method of claim 17, wherein computing the curvature at sample points on the sensed cardiac signal includes computing a curvature versus time curve, and wherein computing an integral of the curvature includes finding the area under a curvature versus time curve.
19. The method of claim 18, wherein extracting features from the computed curvatures includes defining the features by a value of the area under the curvature versus time curve.
20. The method of claim 19, wherein extracting features from the computed curvatures includes comparing the value of the area for the features to an area threshold value and eliminating features with area values less than the threshold area value.
Description
CROSS-REFERENCE TO RELATED APPLICATION(S)

This application is a continuation of U.S. patent application Ser. No. 09/703,269, filed on Oct. 31, 2000, now issued as U.S. Pat. No. 6,684,100, the specification of which is incorporated herein by reference.

FIELD OF THE INVENTION

The present invention relates generally to cardiac rhythm management systems, and more particularly, it pertains to a system and method of classification and detection of cardiac signals.

BACKGROUND OF THE INVENTION

Analysis of cardiac signals, which is routinely performed in electrocardiography is generally based on visual inspection to quantify or qualify wave morphology for the purpose of identifying and classifying abnormal patterns. Certain morphological characteristics of commonly recorded signals have high diagnostic value. The shape and inter arrival times of R-waves recorded in the electrocardiogram generally provide a wealth of information about the state of the heart. Accordingly, automated approaches for identifying and classifying abnormalities in signals such as cardiac signals have sought to use a signal's significant morphologic characteristics.

However, given the wide diversity of possible shapes for cardiac signals, it is usually not possible for an automatic approach to identify significant characteristics that can be used for unambiguous classification. Rather, automated classification approaches generally compare the entire morphological shape of a signal with the shape of similar signals with known abnormalities but without particular regard to the specific characteristics that the signals contain. Alternatively, automated classification approaches restrict the automated examination only to those signals which are essentially normal and use detailed metrics (for example QRS width, QT interval or ST segment amplitude) of the essentially normal morphology for classifying abnormalities.

Despite its importance in the analysis of biologic signals, the automated and accurate identification and quantification of the significant morphological characteristics (for example turns, peaks, knees, inflection points, and the like) in any cardiac signal (both abnormal as well as normal) is still in a developing stage. Existing methods have used the concept of sharpness (for example to detect R-waves) but have had limited success. This is due in part to the overly simplistic mathematical treatment this concept has received, as reflected in the rudimentary algorithms used for these measurements. Most of the current detection methods rely on three point interpolations to measure sharpness. The simplest and most commonly used methods for measuring peaks of R-waves are based upon Taylor-series approximations to estimate the second derivative of the sensed signal. This formula utilizes highly local information (the point at the peak and its two close neighbors) ignoring nearby points which may contribute to signal peak. Other popular approaches utilize less local data, such as the peak and two adjacent extrema. All of these methods, which rely on three-point estimates of sharpness, may produce inaccurate estimates, if waveforms are complex or are contaminated with noise. Thus, a need exists for automated identification and classification of peaks, knees, inflection points, and the like in sensed cardiac signals that takes into account wave scale and complexity that can yield a more accurate estimate of peaks for identifying and classifying abnormalities in cardiac signals.

SUMMARY OF THE INVENTION

The present subject matter provides a curvature based method of selecting features from electrophysiologic signals for purpose of complex identification and classification. According to one aspect of the present subject matter, this is accomplished by sensing a cardiac signal (sensing the cardiac signal includes sensing complexes continuously on a real-time basis) and computing curvatures on a sample point-by-sample point basis (X1, X2, X3, . . . XI) on the sensed cardiac signal on a continuous basis. In one embodiment, the curvature at the sample points X1, X2, X3, . . . XI are computed by fitting a cubic least square error curve (using N number of sample points) to the sensed cardiac signal. In this embodiment, N is an odd number, and the sample point (where the curvature is computed) is at a mid point of the N number of sample points.

Also on a continuous basis, features of significant interest are extracted from the computed curvatures. In some embodiments, this is accomplished by comparing the computed curvatures at the sample points X1, X2, X3, . . . XI to a set of predetermined threshold values. In some embodiments, the features are extracted based on computing features such as a time when the feature occurs, an amplitude of the feature, and other similar features at each of the sample points and comparing them to a set of threshold values.

Also on a continuous basis a set of features associated with a first complex in the sensed cardiac signal are identified and separated from the continuously computed and extracted features upon detecting a second subsequent complex. The second subsequent complex is a complex that is adjacent to the first complex and occurs substantially immediately after the first complex. This process of identifying and separating extracted features repeats itself from one sensed complex to another subsequent sensed complex on a real time basis. One reason for identifying and separating the set of features associated with the first complex is to prevent the features associated with the first complex from mixing with the features associated the second subsequent complex. Separating the features associated with first complex aids in classifying the sensed first complex.

Next, the process includes identifying a fiducial feature from the separated set of extracted features associated with the first complex and aligning the separated set of features with respect to the identified fiducial feature. In one embodiment, fiducial feature is identified based on comparing the times when each of the separated set of features occur with a time when a complex associated with the separated set of features is detected on the sensed cardiac signal, and selecting a feature from the separated set of features that is closest in time to the time when the complex associated with the separated set of features was detected. One reason for using a time when a complex is detected (such as R wave) in identifying the fiducial feature, is because the detection of a complex in a sensed cardiac signal is generally more reliable and consistent.

Next, the process includes aligning the separated set of features around the identified fiducial feature. Aligning the separated set of features around the identified fiducial feature aids in normalizing each of the separated set of feature around a datum such as the associated detected complex, and further aids in comparing the separated set of features with a set of predetermined templates.

Next the process includes comparing the aligned set of features to a set of predetermined templates to classify the associated complex. In some embodiments, the predetermined templates are a set of identified complexes associated with known cardiac arrhythmias that would assist in comparing and classifying the extracted set of features. In some embodiments, a therapy is provided to the heart based on the outcome of the classification. The above described process repeats itself on a continuous basis for a real-time classification of complexes from the sensed cardiac signal.

Other aspects of the invention will be apparent on reading the following detailed description of the invention and viewing the drawings that form a part thereof.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a cardiac signal according to one aspect of the present invention.

FIG. 2 is a curvature figure derived from the cardiac signal of FIG. 1 according to one aspect of the present invention.

FIGS. 3A and 3B illustrate separating and aligning features derived from the curvature figure of FIG. 2 according to one aspect of the present invention.

FIG. 4 is a flow diagram illustrating generally one embodiment of operation of the present subject matter.

FIG. 5 is a schematic/block diagram illustrating generally one embodiment of portions of a cardiac rhythm management system of the present subject matter.

DETAILED DESCRIPTION

In the following detailed description, reference is made to the accompanying drawings which form a part hereof, and in which is shown by way of illustration specific embodiments in which the invention may be practiced. These embodiments are described in sufficient detail to enable those skilled in the art to practice the invention, and it is to be understood that the embodiments may be combined, or that other embodiments may be utilized and that structural, logical and electrical changes may be made without departing from the spirit and scope of the present invention. The following detailed description is, therefore, not to be taken in a limiting sense, and the scope of the present invention is defined by the appended claims and their equivalents.

Analysis of cardiac signals, which is routinely performed in electrocardiography is generally based on visual inspection to quantify or qualify wave morphology for the purpose of identifying and classifying abnormal patterns. Certain morphological characteristics of commonly recorded signals have high diagnostic value. The shape and inter arrival times of R-waves recorded in the electrocardiogram generally provide a wealth of information about the state of the heart. Accordingly, automated approaches for identifying and classifying abnormalities in signals such as cardiac signals have sought to use a signal's significant morphologic characteristics.

However, given the wide diversity of possible shapes for cardiac signals, it is usually not possible for an automatic approach to identify significant characteristics that can be used for unambiguous classification. Rather, automated classification approaches generally compare the entire morphological shape of a signal with the shape of similar signals with known abnormalities but without particular regard to the specific characteristics that the signals contain. Alternatively, automated classification approaches restrict the automated examination only to those signals which are essentially normal and use detailed metrics (for example QRS width, QT interval or ST segment amplitude) of the essentially normal morphology for classifying abnormalities.

Despite its importance in the analysis of biologic signals, the automated and accurate identification and quantification of the significant morphological characteristics (for example turns, peaks, knees, inflection points, and the like) in any cardiac signal (both abnormal as well as normal) is still in a developing stage. Existing methods have used the concept of sharpness (for example to detect R-waves) but have had limited success. This is due in part to the overly simplistic mathematical treatment this concept has received, as reflected in the rudimentary algorithms used for these measurements. Most of the current detection methods rely on three point interpolations to measure sharpness, may produce inaccurate estimates, if waveforms are complex or are contaminated with noise.

GENERAL SYSTEM OVERVIEW

The present subject matter provides, among other things, a system of classifying a cardiac signal by computing a curvature versus time signal from a sensed cardiac signal and comparing the computed curvature versus time signal with a set of predetermined templates. The present system also provides an automated identification and classification of peaks, knees, inflection points, and the like in cardiac signals that takes into account wave scale and complexity that can yield a more accurate classification of cardiac signals.

The first step in the process includes computing continuously a curvature versus time signal from a sensed cardiac signal on a real-time basis. The process also includes extracting features (features are generally significant points of interest on the sensed cardiac signal) continuously on a real-time basis from the computed curvature versus time signal. Next, the process includes continuously separating a set of features associated with a first complex upon detecting a second subsequent complex so that each individual complex can be identified and classified on a real-time basis. Then the process includes aligning and classifying the extracted and separated set of features by comparing the separated set of features with a set of predetermined templates. Other aspects of the invention will be apparent on reading the following detailed description of the invention and viewing the drawings that form a part thereof.

Referring now to FIG. 1, there is one embodiment of the first step in the process of sensing a cardiac signal 100 according to the present invention. Shown in FIG. 1 are detection of heart beats associated with first and second complexes 110 and 120 at times t1′ and t2′ associated with the sensed cardiac signal 100 on a real-time basis on a time line 130. Also, shown on the time line 130 are the various sample points X1, X2, X3, . . . XI used in computing a curvature versus time signal on the first complex 110. In the example embodiment shown in FIG. 1, the sample points can be some fixed interval of time, such as 2, 4, 6 . . . milliseconds, at which curvatures are computed on the sensed cardiac signal. Also shown in this embodiment, is the detection of the complexes (110, 120, and so on) on a continuous basis on the sensed cardiac signal 100.

Referring now to FIG. 2, there is shown one embodiment of generating a curvature versus time signal 200 on a real-time basis derived from the sensed cardiac signal 100 shown in FIG. 1 according to the present invention. FIG. 2 is also on the same time line 130 as in FIG. 1. Curvature FIG. 200 shown in FIG. 2 is a computed curvature (expressed in time) versus time signal generated continuously by computing curvatures at the sample points X1, X2, X3, . . . XI on the sensed cardiac signal 100 shown in FIG. 1. Curvature FIGS. 210 and 220 are computed curvature versus time signals associated with first and second complexes 110 and 120 shown in FIG. 1, respectively. In this example embodiment, the curvature versus time signal is computed by overcoming the dimensionality of the curvature by changing the dimensionality of the curvature to time. The reason for overcoming the dimensionality of the curvature is to convert the time versus voltage signal, which is time versus curvature, to a time versus time signal to facilitate easier computation of curvature (for example, when a curve exists in length-length space, the quantity has dimension of length and is traditionally called the radius of curvature, and if the curve exists in time-time space, the analogous quantity will be called the radius of time). In this time-time space the curvature has a dimension of 1/Time. Also, shown in FIG. 2 are some example sample points 230 at which curvatures are computed.

Referring now to FIGS. 3A and 3B, there is shown one example embodiment of extracting and aligning 300 sets of features associated with respective detected complexes. FIG. 3A, illustrates continuously extracting features 320 from the computed curvature versus time signal 200 shown in FIG. 2. Also shown in this example embodiment, are some sample times in milliseconds 330 when the extracted features occur on the time line 130. FIG. 3A is also on the same time as that of FIGS. 1 and 2. In the example embodiment shown in FIG. 3A, the sample times (where the curvatures are computed) have a fixed interval of 2 milliseconds. Also, in this example embodiment, the beat associated with first complex 110 is detected at 11 milliseconds (t1) and the beat associated with the second subsequent complex 120 is detected at 101 milliseconds (t2).

FIG. 3B shows separating sets of extracted features 310 associated with the complexes 110, 120 . . . on a continuous basis from the continuously extracted features 320 shown in FIG. 3A. In the example embodiment shown in FIG. 3B, a set of features 340 associated with first complex is separated upon detecting the second subsequent complex 120. Also shown are separating a set of features 350 associated with second complex upon detecting a third subsequent complex. This process of separating a set of features associated with a complex is done on a continuous basis. Features associated with first complex is separated from the features of subsequent second complex in order to (so that the features are not mingled with the extracted features of the second subsequent complex) classify the sensed first complex 110. The next step is to identify fiducial features 360, 370, and so on. In the example embodiment shown in FIG. 3B, fiducial feature 360 is identified based on comparing times when the separated sets of extracted features associated with the first complex 110 occur with a time when the heart beat associated with the first complex 110 is detected (t1′), and choosing a feature from the set of features 340 that is closest in time to the time when the heart beat associated with the first complex 110 is detected. In the example shown in FIG. 1, the first complex 110 is detected at 11 milliseconds. The times when the separated set of features 340 associated with the first complex 110 occur are 2, 4, 8, 10, 12, 14, 20, and 24 milliseconds as shown in FIG. 3A. Based on comparing the times when the separated set of features 340 occur with the time when first complex 110 was detected, the feature 360 that occurs at a time of 10 milliseconds is the closest in time to the detected first complex 110 which occurs at 11 milliseconds. Therefore, extracted feature 360 is selected as the fiducial feature from the separated set of features 340 associated with the first complex 110. Similarly, the extracted feature 370 is selected as the fiducial feature for the separated set of features 350 associated with subsequent second complex 120. This process repeats itself continuously on a real time basis. In some embodiments, the fiducial feature is identified from the set of separated features associated with the first complex based on a predetermined deviation value. The predetermined deviation value is based on a sample point having an amplitude farthest from a predetermined reference point.

The next step is to align the separated set of features 340 associated with the first complex 110, and to classify the separated set of features 340 by comparing the aligned set of features with a set of templates. Separated sets of features 340 and 350 are aligned to normalize the separated sets of feature 340 and 350 so that the aligned features aid in comparing the separated sets of features with a template to classify the sensed complexes 110 and 120. FIG. 3B, illustrates one example embodiment of aligning the set of extracted features 340 (associated with first complex 110) with respect to an identified fiducial feature 360. In this example embodiment, the time when the identified fiducial feature occurs (at 10 milliseconds) is subtracted from each of the times when the separated set of features 340 associated with the first complex 110 occurs (at 2, 4, 8, 10, . . . milliseconds) to arrive at aligned times −8, −6, −2, 0, 2, 4, 10, and 14 milliseconds for the separated set of features 340. Similarly, the aligned times −8, −2, 0, 2, 4, and 10 are computed for the separated set of features 350 associated with the second complex 120. This process of aligning separated set of features continues on a real time basis for the detected complexes on the sensed cardiac signal 100. It can be seen from the above illustrated computation, that the process aligning helps normalize the times associated with each of the separated sets of features, and how it aids in comparing the separated set of features 340 and 350 with a set of templates and in classifying the sensed complexes 110 and 120.

FIG. 4, illustrates a method 400 of identifying and classifying a sensed cardiac signal on a real-time basis. Method 400, as shown in FIG. 4, begins with step 410 by sensing a cardiac signal 410 from the one or more electrodes disposed in or around a heart 515 on a continuous basis.

The next step 420 in the process comprises computing curvatures at sample points X1, X2, X3, . . . XI on the sensed cardiac signal. In some embodiments, the sample points X1, X2, X3, . . . XI are fixed intervals of time at which the sensed cardiac signal is computed for its curvature. The sensed cardiac signal can be an atrial or ventricular signal of a heart. In one embodiment, the curvatures at sample point X1, X2, X3, . . . XI is computed by fitting a cubic least square error curve to the sensed cardiac signal by using N sample points. In this embodiment, N is an odd number greater than or equal to 5. Also in this embodiment, the cardiac rhythm management system computes the curvature at a mid point of the N sample points. In one embodiment, the cardiac rhythm management system computes curvature (K) at a sample point XI on the cardiac signal by using 5 sample points over which the cubic fit is (N being 5 in this example) made using
K=(2C I/(1+B I 2)3/2)
where
B I =T 1 Y(I−2)+T 2 Y(I−1)+T 3 Y(I)+.T 4 Y(I+1)+T 5 Y(I+2)
C I =S 1 Y(I−2)+S 2 Y(I−2)+S 3 Y(I)+S 4 Y(I+1)+S 5 Y(I+2)
where S1, S2, S3 . . . and T1, T2, T3 . . . are constant terms in the summation. In this embodiment, I corresponds to the sample point at which the curvature is being computed. For example, if the sample point is X3, then I will be equal to 3 and so on.

In another embodiment, the system computes average curvature between two adjacent sample points by using a linear interpolation of the BI and CI between the two adjacent sample points (this is done by doing the integral of the curvature between the two adjacent sample points). This integration helps capture any curvatures that become larger between the two adjacent sample points, so that the point-wise curvature values may more robustly represent the curvatures in the sensed cardiac signal. The system also integrates curvature versus time curve to find the area under the curve. This area under the curve represent the total change in angle between two points on the original cardiac signal. Thus, peaks, turns, inflections, and the like in the cardiac signals can be identified when the curvature versus time curve has peaks or nadirs. These peaks and nadirs denote points along the sensed cardiac signal where the turns are locally sharpest. Also, slow but steady turns on the sensed cardiac signal can be characterized by a small peak having a substantial area under the curvature curve. Thus, integrating curvatures between two adjacent sample points can help capture any curvatures that may become larger between any two adjacent sample points.

In one embodiment, the cardiac rhythm management system computes curvatures at sample points X1, X2, X3, . . . XI on the sensed cardiac signal on a continuous basis.

Then the next step 430 in the process includes extracting features from the computed curvatures by comparing the computed curvatures with a set of predetermined threshold values. FIG. 3A, illustrates extracting features from the computed curvatures at the sample points X1, X2, X3, . . . XI. In some embodiments, the set of predetermined threshold values are based on a previous curvature value, a first curvature value, and a curvature threshold limit. In one embodiment, the step of extracting features from the curvature versus time signal is done on a continuous basis on the sensed cardiac signal as shown in FIG. 3A. In some embodiments, the features are defined by one or more metrics. In this embodiment, the one or more metrics are defined as area under a computed curvature, a time of centroid of the area, and a value of original signal amplitude at a time of centroid of the area associated with the feature.

Then the next step 440 comprises separating a set of features associated with a first complex upon detecting a second subsequent complex from the sensed cardiac signal. Separating of extracted features associated with the detected complex is illustrated in detail FIGS. 1, 2, 3A and 3B. In some embodiments, the second subsequent complex is detected based on continuously processing the sensed cardiac signal to produce a signal that is the absolute value of the first derivative of a sensed cardiac signal versus time. Then comparing the produced absolute value versus time signal with a predetermined decaying threshold value to detect an occurrence of the second subsequent complex. In some other embodiments, the second subsequent complex is detected when the absolute value of the first derivative exceeds a predetermined decaying threshold value. In some embodiments, the features associated with the first complex are stored in a memory until the second subsequent complex is detected, and at that point, the features associated with the first complex are replaced by the features associated with the second subsequent complex in the memory. The step of separating can comprise separating the set of features based on selecting features associated with the first complex having a predetermined time earlier than the detected second subsequent complex. In some embodiments, each of the set of features are further compared to an area threshold value to eliminate features less than the area threshold value. This can result in a set of revised features that can represent the first complex.

Then the next step 450 in the process is to identify a fiducial feature from the set of separated features associated with the first complex based on a predetermined deviation value. In some embodiments, fiducial feature is identified based on comparing the time when the first complex occurs with the times when the separated set of features associated with the first complex occur. The process of identifying a fiducial feature is discussed in more detail with respect to FIGS. 3A and 3B.

Then the next step 460 in the process comprises aligning the set of separated features with respect to the identified fiducial feature. The process of aligning is also discussed in more detail in FIG. 3B. The process of aligning helps normalize the times associated with each of the separated set of features and further facilitates in comparing the separated set of features with a set of predetermined templates to classify the detected complex.

Then the next step 470 comprises comparing the aligned set of separated features with a set of predetermined templates. In one embodiment, the set of predetermined templates are one or more known types of heart beat signals. In another embodiment, the set of predetermined templates consists of sets of template boxes or zones for each of the known types of heart beat signals. Generally, each template box or zone can have center amplitude, a time width, and an amplitude width. Also generally, the template boxes or zones for each known type of heart beat can consist of one or more template boxes or zones each of which can have different centers and widths. In this embodiment, the comparing step comprises assigning a score based on how well each feature matches with each of the set of template boxes or zones. In this embodiment, if the separated set of features associated with the first complex fails to match with the set of templates then the complex is assigned with an unknown heart beat classification. In one embodiment, the set of template boxes or zones are modified over time to reflect any changes in sets of separated features used in matching with the template boxes or zones. This can permit known beat type templates to follow slow trends or long-term changes in the cardiac signal.

The next step 480 in the process comprises classifying the detected first complex based on the outcome of the comparison. Then the next step 490 comprises providing a therapy based on the outcome of the classification. In one embodiment, the next step includes guiding a therapy to a heart based on the outcome of the classification. In another embodiment, the next step in the process includes storing classifications for future diagnostic purposes.

FIG. 5, is a schematic/block diagram 500 illustrating one embodiment of portions of a cardiac rhythm management system 505, which is coupled to a heart 515. Cardiac rhythm management system 505 includes a power source 580, a controller 540, a sensing circuit 520, a therapy circuit 590, and a via node/bus 530. The sensing circuit 520 is coupled by a lead 510 to the heart 515 for receiving, sensing, and or detecting electrical heart signals. The sensing circuit 520 provides one or more sensed cardiac signals to the controller 540, via the node/bus 530. The controller 540 also controls the delivery of a therapy provided by the therapy circuit 590 and/or other circuits, as discussed below.

The controller 540 includes various modules, which are implemented either in hardware or as one or more sequences of steps carried out on a microprocessor or other controller. Such modules are illustrated separately for conceptual clarity; it is understood that the various modules of the controller 540 need not be separately embodied, but may be combined and/or otherwise implemented, such as in software/firmware.

In general terms, the sensing circuit 520 senses a cardiac signal on a continuous basis from a heart tissue in contact with a catheter lead 510 to which the sensing circuit 520 is coupled. In one embodiment, the sensed cardiac signal is an atrial signal. In another embodiment, the sensed cardiac signal is a ventricular signal. Sensed cardiac signal from the sensing circuit 520 is then received and processed by an analyzer 550 of the controller 540 to compute curvatures continuously on a real-time basis at sample point X1, X2, X3, . . . XI on the sensed cardiac signal. In some embodiments, the sample points are spaced at fixed intervals of time preset by the analyzer 550 on the sensed cardiac signal. In one embodiment, the analyzer 550 computes curvature at sample points X1, X2, X3, . . . XI by fitting a cubic least square error curve to the sensed cardiac signal using N sample points to fit the cubic least square error curve. In one embodiment, the analyzer 550 computes the computing curvature (K) at a sample point XI on the cardiac signal when using 5 sample points over which the cubic fit is (N being 5 in this example) made using
K=(2C I/(1+B I 2)3/2)
where
B I =T 1 Y(I−2)+T 2 Y(I−1)+T 3 Y(I)+T 4 Y(I+1)+T 5 Y(I+2)
C I =S 1 Y(I−2)+S 2 Y(I−2)+S 3 Y(I)+S 4 Y(I+1)+S 5 Y(I+2)
where S1, S2, S3 . . . and T1, T2, T3 . . . are constant terms in the summation. In this example embodiment, I is equal to the point at which the curvature is being computed. For example, if the curvature is being computed at sample point X3, then I will be equal to 3 and so on. In one embodiment, N is an odd number greater than or equal to 5. In this embodiment, the analyzer 550 computes curvature at a midpoint of the N number of sample points.

In another embodiment, the analyzer 550 computes curvature based on an average curvature between two adjacent sample points based on a linear interpolation of the BI and CI between two adjacent points (integrating the computed curvatures between the two adjacent sample points). In one embodiment, the analyzer 550 computes the curvatures continuously on a real-time basis from the sensed cardiac signal. Then the analyzer 550 extracts features by comparing the computed curvatures to a set of predetermined threshold values. In one embodiment, the predetermined threshold values are based on a previous curvature value, a curvature value, and a curvature threshold limit. In one embodiment, the analyzer 550 extracts features on a continuous basis. Then the analyzer 550 separates a set of features associated with a first complex upon detecting a second subsequent complex from the sensed cardiac signal. Then the analyzer 550 identifies a fiducial feature from the set of separated features associated with the first complex based on a predetermined deviation value. Then the analyzer 550 further aligns the set of separated features with respect to the identified fiducial feature. In one embodiment, the analyzer 550 separates the set of features associated with the first complex based on a predetermined time earlier than the detected second subsequent complex. In one embodiment, the predetermined deviation value is based on a sample point having an amplitude farthest from a predetermined reference point.

Then a comparator 570 coupled to the analyzer 550, compares the aligned set of features associated with the first complex with a predetermined set of templates, and classifies the first complex based on the outcome of the comparison. In one embodiment, the comparator 570 issues a command signal based on the outcome of the classification. In another embodiment, the therapy circuit 590 coupled to the comparator 570 delivers the electrical energy through the lead 510 to at least one electrode disposed in or around the heart 515 upon receiving the command signal from the comparator 570. The electrical energy delivered by the therapy circuit 590 can be a pacing pulse electrical energy. In one embodiment, the analyzer 550 includes a variable gain circuitry that can adapt according to the changes in the sensed cardiac signal. In one embodiment, the controller 540 includes a memory 560 coupled to the comparator 570 and the analyzer 550 to store the extracted features of the first complex. In another embodiment, the memory 560 stores the classified first complexes for diagnostic purposes.

It is to be understood that the above description is intended to be illustrative, and not restrictive. Many other embodiments will be apparent to those of skill in the art upon reviewing the above description. The scope of the invention should, therefore, be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled.

CONCLUSION

The above described method provides, among other things, a curvature based complex identification and classification. The process comprises sensing a cardiac signal and computing curvatures at sample points on the sensed cardiac signal, and then to extract features from the computed curvatures. Then to compare the extracted features with a set of predetermined templates, and then to classify the cardiac signal based on the outcome of the comparison.

Patent Citations
Cited PatentFiling datePublication dateApplicantTitle
US382876813 Jul 197213 Aug 1974Physiological Electronics CorpMethod and apparatus for detecting cardiac arrhythmias
US39398249 Jan 197524 Feb 1976General Electric CompanyPhysiological waveform detector
US399821419 May 197521 Dec 1976Brondy Laboratories, Inc.Premature ventricular contraction detector and method
US40004614 Oct 197328 Dec 1976Textronix, Inc.R-wave detector
US402356426 Jan 197617 May 1977Spacelabs, Inc.Arrhythmia detector
US433681030 Sep 198029 Jun 1982Del Mar AvionicsMethod and apparatus for arrhythmia analysis of ECG recordings
US453152723 Apr 198230 Jul 1985Survival Technology, Inc.Ambulatory monitoring system with real time analysis and telephone transmission
US458355315 Nov 198322 Apr 1986Medicomp, Inc.Ambulatory ECG analyzer and recorder
US458942013 Jul 198420 May 1986Spacelabs Inc.Method and apparatus for ECG rhythm analysis
US463740013 Aug 198520 Jan 1987Marcus Edward NFor characterizing a heart condition
US472111421 Feb 198626 Jan 1988Cardiac Pacemakers, Inc.Method of detecting P-waves in ECG recordings
US480249130 Jul 19867 Feb 1989Massachusetts Institute Of TechnologyMethod and apparatus for assessing myocardial electrical stability
US480969714 Oct 19877 Mar 1989Siemens-Pacesetter, Inc.Interactive programming and diagnostic system for use with implantable pacemaker
US482586928 Sep 19872 May 1989Telectronics N.V.System for automatically performing a clinical assessment of an implanted pacer based on information that is telemetrically received
US483203820 Jan 198723 May 1989The Board Of Trustees Of University Of IllinoisApparatus for monitoring cardiovascular regulation using heart rate power spectral analysis
US483827826 Feb 198713 Jun 1989Hewlett-Packard CompanyPaced QRS complex classifier
US488434512 May 19885 Dec 1989Siemens-Pacesetter, Inc.Adjustable template for pacemaker ECG analysis and method of use
US49248759 Oct 198715 May 1990Biometrak CorporationMethod of detecting heart disorders
US494591723 May 19887 Aug 1990Solange AkselrodMethod and apparatus for indicating repetition intervals of a specified component of a composite electrical signal, particularly useful for displaying fetal r-waves
US49478571 Feb 198914 Aug 1990Corazonix CorporationMethod and apparatus for analyzing and interpreting electrocardiograms using spectro-temporal mapping
US498961016 Nov 19875 Feb 1991Spacelabs, Inc.Method and system of ECG data review and analysis
US500018915 Nov 198919 Mar 1991Regents Of The University Of MichiganMethod and system for monitoring electrocardiographic signals and detecting a pathological cardiac arrhythmia such as ventricular tachycardia
US501428430 Jun 19897 May 1991Cardiac Telecom CorporationDiscrete slope delta modulation with recovery means
US50146982 Oct 198914 May 1991Leonard BloomMethod of and system for monitoring and treating a malfunctioning heart
US502054022 Aug 19894 Jun 1991Biometrak CorporationCardiac biopotential analysis system and method
US50465046 Apr 199010 Sep 1991Corazonix CorporationMethod and apparatus for analyzing and interpreting electrocardiograms using spectro-temporal mapping
US51078502 Nov 199028 Apr 1992Cardiac Pacemakers, Inc.Method and apparatus for classifying and treating cardiac arrhythmias based on atrial and ventricular activity
US510984224 Sep 19905 May 1992Siemens Pacesetter, Inc.Implantable tachyarrhythmia control system having a patch electrode with an integrated cardiac activity system
US513902826 Oct 199018 Aug 1992Telectronics Pacing Systems, Inc.Heart rejection monitoring apparatus and method
US515614813 May 199120 Oct 1992Leonard BloomIn an implantable system
US51846158 Mar 19919 Feb 1993Telectronics Pacing Systems, Inc.Apparatus and method for detecting abnormal cardiac rhythms using evoked potential measurements in an arrhythmia control system
US519355030 Nov 199016 Mar 1993Medtronic, Inc.Method and apparatus for discriminating among normal and pathological tachyarrhythmias
US520132111 Feb 199113 Apr 1993Fulton Keith WMethod and apparatus for diagnosing vulnerability to lethal cardiac arrhythmias
US521509812 Aug 19911 Jun 1993Telectronics Pacing Systems, Inc.Data compression of cardiac electrical signals using scanning correlation and temporal data compression
US521702130 Jul 19918 Jun 1993Telectronics Pacing Systems, Inc.Detection of cardiac arrhythmias using correlation of a cardiac electrical signals and temporal data compression
US524000925 Mar 199131 Aug 1993Ventritex, Inc.Medical device for monitoring intracardiac electrograms
US524702127 Aug 199221 Sep 1993Kanegafuchi Kagaku Kogyo Kabushiki KaishaProcess for preparation of a polymer having reactive terminal group
US52551866 Aug 199119 Oct 1993Telectronics Pacing Systems, Inc.Signal averaging of cardiac electrical signals using temporal data compression and scanning correlation
US526930119 Oct 199214 Dec 1993Leonard BloomMultimode system for monitoring and treating a malfunctioning heart
US527141110 Sep 199121 Dec 1993Colin Electronics Co., Ltd.Method and apparatus for ECG signal analysis and cardiac arrhythmia detection
US52730499 Apr 199228 Dec 1993Telectronics Pacing Systems, Inc.Detection of cardiac arrhythmias using template matching by signature analysis
US527562113 Apr 19924 Jan 1994Medtronic, Inc.Method and apparatus for terminating tachycardia
US528079218 Sep 199225 Jan 1994The University Of SydneyMethod and system for automatically classifying intracardiac electrograms
US52914009 Apr 19921 Mar 1994Spacelabs Medical, Inc.System for heart rate variability analysis
US52923412 Mar 19928 Mar 1994Siemens Pacesetter, Inc.Method and system for determining and automatically adjusting the sensor parameters of a rate-responsive pacemaker
US529234813 Apr 19928 Mar 1994Telectronics Pacing Systems, Inc.Implantable cardioverter/defibrillator and method employing cross-phase spectrum analysis for arrhythmia detection
US531187418 May 199217 May 1994Cardiac Pacemakers, Inc.Method for tachycardia discrimination
US53124453 Feb 199217 May 1994Telectronics Pacing Systems, Inc.Implantable cardiac stimulating apparatus and method employing detection of P-waves from signals sensed in the ventricle
US531395314 Jan 199224 May 1994Incontrol, Inc.For monitoring the physiology of the human heart
US53305043 May 199319 Jul 1994Telectronics Pacing Systems, Inc.Cardioverting defibrillating device with off-line ECG analysis
US533196616 Dec 199326 Jul 1994Medtronic, Inc.Subcutaneous multi-electrode sensing system, method and pacer
US534240229 Jan 199330 Aug 1994Medtronic, Inc.Method and apparatus for detection and treatment of tachycardia and fibrillation
US535040622 Dec 199227 Sep 1994Ela MedicalApparatus and method of control for an implanted anti-tachycardia device
US535169623 Apr 19934 Oct 1994Medtronic, Inc.Method and apparatus for intracardiac electrogram morphologic analysis
US53604363 Apr 19921 Nov 1994Intermedics, Inc.Cardiac pacing responsive to multiple activity types
US53664874 Nov 199322 Nov 1994Angeion CorporationFibrillation and tachycardia detection
US53885786 Apr 199314 Feb 1995Incontrol, Inc.Electrode system for use with an implantable cardiac patient monitor
US540079522 Oct 199328 Mar 1995Telectronics Pacing Systems, Inc.Method of classifying heart rhythms by analyzing several morphology defining metrics derived for a patient's QRS complex
US541103124 Nov 19932 May 1995Incontrol, Inc.Implantable cardiac patient monitor
US541152913 Sep 19932 May 1995Medtronic, Inc.Waveform discriminator for cardiac stimulation devices
US54172217 May 199323 May 1995Psytech, Inc.Method and apparatus for distinguishing electric signal waveforms
US542183027 Aug 19936 Jun 1995Pacesetter, Inc.Programming system having means for recording and analyzing a patient's cardiac signal
US543169310 Dec 199311 Jul 1995Intermedics, Inc.Method of verifying capture of the heart by a pacemaker
US544751919 Mar 19945 Sep 1995Medtronic, Inc.Method and apparatus for discrimination of monomorphic and polymorphic arrhythmias and for treatment thereof
US544752425 Jul 19945 Sep 1995Intermedics, Inc.Cardiac pacing method and apparatus responsive to multiple activity types
US545626116 Dec 199310 Oct 1995Marquette Electronics, Inc.Cardiac monitoring and diagnostic system
US545862014 Feb 199417 Oct 1995Angeion CorporationCardiac arrhythmia detection using interdependent multiple parameter settings
US54586234 Mar 199417 Oct 1995Telectronics Pacing Systems, Inc.Automatic atrial pacing threshold determination utilizing an external programmer and a surface electrogram
US546443314 Jun 19947 Nov 1995Incontrol, Inc.Atrial defibrillator and method providing dual reset of an interval timer
US547880723 Jun 199226 Dec 1995Genentech, Inc.Increases rate of stimulus provided by subsidiary pacemakers
US548041216 Nov 19932 Jan 1996Pacesetter, Inc.System and method for deriving hemodynamic signals from a cardiac wall motion sensor
US550315917 Mar 19952 Apr 1996Hewlett-Packard CompanyMethod for enhancement of late potentials measurements
US55099271 Feb 199523 Apr 1996Pacesetter, Inc.Programming system having means for recording and analyzing a patient's cardiac signal
US551155427 Jul 199430 Apr 1996Hewlett Packard CompanyReal-time artifact removal from waveforms using a dynamic filter having a fixed delay
US551364422 Sep 19947 May 1996Pacesetter, Inc.Cardiac arrhythmia detection system for an implantable stimulation device
US55201917 Oct 199428 May 1996Ortivus Medical AbMyocardial ischemia and infarction analysis and monitoring method and apparatus
US554243016 Sep 19946 Aug 1996Telectronics Pacing Systems, Inc.Apparatus and method for discriminating between cardiac rhythms on the basis of their morphology using a neural network
US554518630 Mar 199513 Aug 1996Medtronic, Inc.Prioritized rule based method and apparatus for diagnosis and treatment of arrhythmias
US554964131 Oct 199427 Aug 1996Incontrol, Inc.Atrial fibrillation type selective cardiovertor and method
US554965415 Apr 199427 Aug 1996Medtronic, Inc.Interactive interpretation of event markers in body-implantable medical device
US556036922 Sep 19941 Oct 1996Pacesetter, Inc.Cardiac arrhythmia detection system for an implantable stimulation device and method
US560515916 Feb 199625 Feb 1997Smith; Joseph M.System and method for determining spatial organization of atrial activation
US56091581 May 199511 Mar 1997Arrhythmia Research Technology, Inc.Apparatus and method for predicting cardiac arrhythmia by detection of micropotentials and analysis of all ECG segments and intervals
US56221784 May 199422 Apr 1997Spacelabs Medical, Inc.System and method for dynamically displaying cardiac interval data using scatter-plots
US562832629 Nov 199513 May 1997Hewlett-Packard CompanyCalculating a heart rate from an ECG waveform by discarding a percentage of R-R intervals prior to averaging
US56344682 Apr 19933 Jun 1997Micromedical Industries LimitedEcg monitoring system
US564507025 Sep 19958 Jul 1997Ventritex, Inc.Method and apparatus for determining the origins of cardiac arrhythmias morphology dynamics
US567615325 Oct 199614 Oct 1997Smith; Joseph M.System and method for determining spatial organization of atrial activation
US568290029 Nov 19954 Nov 1997Hewlett-Packard CompanyMethod and apparatus for obtaining heartbeat measurements from a ECG waveform
US568290216 Oct 19954 Nov 1997Hewlett-Packard CompanyECG pace pulse detection and processing
US568342512 Dec 19954 Nov 1997Pacesetter AbHeart pacemaker with improved detection of electrical signals
US568531529 Mar 199611 Nov 1997Pacesetter, Inc.Cardiac arrhythmia detection system for an implantable stimulation device
US568773731 Oct 199418 Nov 1997Washington UniversityComputerized three-dimensional cardiac mapping with interactive visual displays
US570436514 Nov 19956 Jan 1998Cambridge Heart, Inc.Using related signals to reduce ECG noise
US571280125 Sep 199527 Jan 1998Pacesetter, Inc.Method for characterizing dynamical systems
US571336616 Sep 19963 Feb 1998Sulzer Intermedics Inc.Method and apparatus for dual chamber cardiac analysis
US571336726 Jan 19953 Feb 1998Cambridge Heart, Inc.Measuring and assessing cardiac electrical stability
US571824226 Jun 199617 Feb 1998Pacesetter, Inc.Cardiac arrhythmia detection system for an implantable stimulation device and method
USRE3514820 Jan 199523 Jan 1996Riverside Research InstituteFrequency diversity for image enhancement
Non-Patent Citations
Reference
1"PCT International Search Report", from International Application No. PCT/US01/45763, dated Sep. 3, 2002,8.
2Afonso, V. X., et al., "ECG Beat Detection Using Filter Banks", IEEE Transactions on Biomedical Engineering, 46, (Feb. 1999),192-202.
3Afonso, V. X., et al., "Filter Bank-based ECG Beat Classification", Proceedings of the 19th International Conference-IEEE/EMBS, Chicago, IL,(Oct. 30-Nov. 2, 1997),97-100.
4Afonso, V. X., et al., "Filter Bank-based ECG Beat Detection", Proceedings of the 18th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Amsterdam, The Netherlands,(1996),1037-1038.
5Afonso, V. X., et al., "Filter Bank-based Processing of the Stress ECG", 1995 IEEE-EMBC and CMBEC-Theme 4: Signal Processing, (1995),887-888.
6Afonso, V. X., et al., "Multirate Processing of the ECG using Filter Banks", Computers in Cardiology, (1996),245-248.
7Duru, Firat , et al., "The Potential for Inappropriate Ventricular Tachycardia Confirmation Using the Intracardiac Electrogram (EGM) Width Criterion", PACE, vol. 22, (Jul. 1999),1039-1046.
8Grady, Thomas A., et al., "Prognostice Significance of Exercise-Induced Left Bundle-Branch Block", JAMA, vol. 279, No. 2, Jan. 14, 1998,153-156.
9IEEE, "Microcomputer-Based Telemetry System for ECG Monitoring", IEEE Proc. of the Ann. Int'l Conf. of the Engineering in Medicine and Biology Society, vol. Conf. 9, XP000015425, (1987),1492-193.
10Kinoshita, Shinji , et al., "Transient Disapperance of Complete Right Bundle Branch (BBB) During Exercise", Journal of Electrocardiology, vol. 29, No. 3, 1996,255-256.
11Sweeney, Robert J., "Signal Compression Based on Curvature Parameters", U.S. Appl. No. 10/607,820, filed Jun. 27, 2003, 50 pgs.
12Sweeney, Robert J., et al., "Tachyarrahythmia Detection and Discrimination Based on Curvature Parameters", U.S. Appl. No. 10/607,818, filed Jun. 27, 2003, 59 pgs.
13Vapnik, V. N., "Chapter 5 Constructing Learning Algorithms", The Nature of Statistical Learning Theory, Springer-Verlag New York, Inc.,(1995),119-166.
Referenced by
Citing PatentFiling datePublication dateApplicantTitle
US7761162 *13 Dec 200420 Jul 2010Cardiac Pacemakers, Inc.Capture verification with intrinsic response discrimination
US816566624 Aug 201124 Apr 2012Topera, Inc.System and method for reconstructing cardiac activation information
US852127611 Jan 201127 Aug 2013Cardiac Pacemakers, Inc.Use of significant point methodology to prevent inappropriate therapy
US85947773 Apr 201226 Nov 2013The Reagents Of The University Of CaliforniaSystem and method for reconstructing cardiac activation information
US874455624 Jan 20123 Jun 2014Cardiac Pacemakers, Inc.Noise detection in implantable medical devices
Classifications
U.S. Classification600/517, 600/515
International ClassificationA61B5/0452, A61B5/04, G06F17/00, A61B5/0472, G06K9/00, A61B5/0456, A61N1/362
Cooperative ClassificationG06K9/00523, A61B5/04525, A61B5/7264, A61N1/362
European ClassificationA61B5/72K12, A61N1/362, G06K9/00M2, A61B5/0452B
Legal Events
DateCodeEventDescription
30 Mar 2011FPAYFee payment
Year of fee payment: 4