ECG APPARATUS

Abstract: Apparatus for classifying QRS complexes as dual chamber paced, atrially paced or ventricularly paced by assigning (10) complexes to groups in accordance with their shape, deriving the intervals (20) between pace pulses and their respective QRS complexes, determining (22) if the intervals for a given number of recent complexes from the group to which the current QRS complex is assigned are clustered on a histogram plot, and classifying (26) the current QRS complex by the correlation between the intervals of its pace pulses and the intervals of said clusters.

Description: Background of the Invention The analysis of ECG waveforms using computerized arrhythmia monitoring systems consists of several steps. First each QRS complex must be detected and classified into normal and abnormal categories. Subsequent analysis on the timing and pattern of these beats will yield alarms with several priorities to be annunciated to the user for treatment of the patient. The patient population can be broken into two groups, patients with pacemakers and patients without pacemakers. Non-paced patients generate abnormal QRS complexes i.e. beats that are quite different in shape from their normal counterparts. However, paced patients generate paced beats initiated by the pacemaker that are treated as normal beats. These paced beats are substantially similar in wave shape to their abnormal counterparts. Thus computerized arrhythmia monitoring systems designed for non-paced patients generate a high level of false alarms. Some systems employ a simplistic approach to reduce false alarms by classifying any beat that is preceded by a pace pulse, as a paced beat. However this does not significantly reduce false alarms. Pace pulses can appear to be generated by noise in the ECG signals or be ineffective pace pulses that have not initiated any beat, as to indicate an abnormal pacemaker operation. Such ineffectual pace pulses and false pace pulse indications may preceed abnormal beats which are misclassified as paced beats . A group of shaped QRS complexes, that are pacemaker initiated, have pace pulses at a constant distance to the beats in this group. The interval between the QRS complex and the preceding pace pulse is larger for atrially paced beats than ventricularly paced beats because of the delay in normal conduction of the pacing stimulus from the atrium to the ventricles. This relationship between atrial and ventricular pace pulse intervals to the QRS complex remains the same for dually paced beats. In the prior art, EP - A0080821 discloses a heart analysing system for an ambulatory patient, which uses QRS sequence classification and also describes pacemaker/cardiac activity classifications. However, this document is not concerned with the criteria for distinguishing between different pacing modes, nor does it describe ways of achieving this. Document WO-A 86/05698 describes an analysis system in a pacemaker, so the pacing modes are inherent in the system, and the problem behind the present invention is not addressed. Brief Summary of the Invention If the times of occurrence of pace pulses and QRS complexes could be determined without error, and if the intervals between them remained constant, it would not be difficult to classify the complexes as being dual chamber paced, atrially paced, ventricularly paced, normal or other and therefore not paced, but means for detecting the time of occurrence of a pace pulse can give false positive or false negative indications because of the presence of noise, artifacts and body movement as well as physiological factors that can change the shape of a QRS complex so as to make it difficult to say precisely when it occurs. Furthermore, the pacemakers themselves may provide erroneous signals or may not produce the desired signals so that classification on a beat by beat basis would be in error. In fact one of the reasons for classifying the QRS complexes in an ECG wave is to monitor the operation of the pacer. In accordance with this invention, there is provided apparatus for classifying the type of pacing of QRS complexes contained in an ECG wave signal comprising, first means for detecting the time of occurrence of said QRS complexes in an ECG wave signal; second means for detecting the times of occurrence of pace pulses in an ECG wave; third means coupled to the second means for determining the time intervals between each pace pulse and the next QRS complex; and grouping means for assigning said QRS complexes to groups in accordance with their shape; characterised by clustering means coupled to the third and the grouping means for detecting the existence of at least one cluster of said intervals occurring during a given number of the most recent QRS complexes; means for deriving an average interval value for each said cluster; and classifying means for classifying the current QRS complex by the correlation between said intervals associated with it and the averages for the clusters of the group to which it has been assigned. In practice, the time at which a QRS complex is received is provided and the complex is stored with a group of complexes having the same shape in accordance with existing techniques. A pace pulse detector provides the times of occurrence of pace pulses, and means are provided for determining the interval between each pace pulse and the next QRS complex. These intervals are stored with that complex. Means are provided for determining whether the intervals for the last complexes of the group to which the current complex has been assigned are clustered about different values, and other means derive the average value of the intervals occurring in each cluster. A small average value may be designated as ds, and a large average value may be designated as dl. These values are stored in a template buffer. If no clusters are found, there are no average values and nothing to be stored in the template buffer. The occurrence of a value ds is an indication that the QRS's in the group to which the current complex has been assigned are ventricularly paced; the occurrence of a value dl is an indication that they are atrially paced; this occurrence of ds and dl is an indication that they are dual chamber paced, and the absence of both ds and dl is an indication that they are derived from normal beats in which the heart is not reacting to a pace pulse even though it may be present. These indications are stored in the template buffer. Whereas the probable classification of the current QRS complex can be derived from these indications, a more reliable classification is obtained from certain known facts and the correlation between these indications and the intervals for pace pulses associated with the current QRS complex. By way of example, if two clusters at ds and dl are found for the group of QRS complexes to which the current complex has been assigned and the complex has a pace pulse at an interval that is very close to dl, it will be classified as dual-chamber paced. Other correlations will be described in the detailed description of the invention. Brief Description of the Drawings Figure 1 illustrates a series of QRS complexes of two different groups or shapes A and B and their accompanying pace pulses, Figure 1A shows the complexes of the group having the shape A in vertical alignment so that the distances of the single pace pulses from their respective QRS complexes may be more readily compared, Figure 1A′ illustrates the clustering of the distances for the group A, Figure 1B shows the complexes of the group B in vertical alignment so that the distances of the dual pace pulses from their respective QRS complexes may be more readily compared, Figure 1B′ shows the clustering of the distances for the group of complexes shown in Figure 1B, Figure 2 shows an ECG signal and the classifications of the QRS complexes that would be made by monitoring apparatus incorporating this invention, Figure 3 is a block diagram of a monitoring apparatus incorporating this invention, Figure 4 illustrates what is stored in a template storage buffer of figure 3, Figure 5 is a flow chart for identifying Pace Pulses associated with QRS complexes. Figure 6A is a flow chart for determining certain values related to pace pulse clustering, Figure 6B is a flow chart for determining other values related to pace pulse clustering, Figure 7 is a flow chart for utilizing the values provided in accordance with the flow charts of Figures 6A and 6B for the purpose of updating a template classification for each group, and Figure 8 is a flow chart for determining the classification of a current QRS complex. Detailed Description of the Invention In Figure 1 six QRS complexes Q1 through Q6 are shown. The complexes Q1, Q3 and Q6 form a group having shapes that while not identical are sufficiently similar to be classified as a shape A, and the complexes Q2, Q4 and Q5 form a different group having shapes that while not identical are sufficiently similar to be classified as a shape B. Single pace pulses are shown ahead of each of the complexes having the shape A. These complexes are vertically aligned in Figure 1A so that it can be seen that the intervals between the single pace pulses and their respective QRS complexes are very close to the same values. Figure 1A′ is a histogram showing at CA the clustering or distribution of these intervals. A pair of pace pulses are shown ahead of each of the complexes having the shape B. These complexes are vertically aligned in Figure 1B so that it can be seen that the intervals between corresponding ones of each pair of pace pulses and their respective QRS complex are very close to the same value. Figure 1B′ is a histogram showing at CB the clustering or distribution of short intervals for pace pulses closer to their respective complexes and at CB′ the clustering or distribution of long intervals for the pace pulses farther from their respective complexes. Figure 2 shows a series of QRS complexes having an N over them when they would be classified as normal by the apparatus of this invention because the distance between the pace pulses and their respective QRS complexes varies. This is a proper classification, but those that are underlined would be erroneously classified as paced by an algorithm that asserts such a classification merely because the pace pulse accompanies the QRS complex. Figure 3 is a block diagram of a QRS classification apparatus that operates in accordance with this invention. ECG signals are applied to an input terminal Ti. After passing through a band pass filter 2 that may have a lower cut-off frequency of 0.5 Hz for eliminating base line wander, the ECG signals are applied to an A/D converter 4. If its sampling rate is 500 Hz, the provision of an upper cut-off frequency of 125 Hz for the band pass filter 2 has been found satisfactory. The samples emerging from the A/D converter 4 are applied via a shift register 6 to a QRS detector 8 that, as is well known by those skilled in the art, supplies a signal at its output indicating the time of the peak of the R wave of a QRS complex. A QRS complex grouper 10 that is coupled to the shift register 6 and the output of the QRS detector 8 updates beat shapes it previously stored in a template storage buffer 12, representing different shapes of QRS complexes. The function is well understood by those skilled in the art so as to require no further description. The contents of the template buffer 12 illustrated in Figure 4 contain updated templates of the various shapes of QRS complexes encountered, the number of them that have occured, and the classification of each shape. A pace pulse detector 14 is connected between the input terminal Ti and an A/D converter 16, and its output is applied to a shift register 18. A means 20 is provided for establishing the differences in time between the current QRS complex detected by the QRS detector 8 and the pace pulses detected by the detector 14. The pace pulses are those that occurred within a given interval between the current QRS complex and the previous one. The output of the QRS detector 8 and the shift register 18 are coupled to the means 20. The difference in time is stored in the template storage buffer 12 in association with its QRS complex. A means 22 is provided for deriving the clustering information for each of the last QRS complexes of the shape group to which the current QRS belongs. In particular, it provides a signal OPs=1 if at least six of the last 8 of the QRS complexes of the shape group including the current QRS have pace pulse intervals within a narrow cluster centered at a short time prior to their respective QRS complexes, in which event the average ds of all these intervals is derived. If the condition does not exist set OPs=0. Similarly, the means 22 provides a signal OPl=1 if at least six of the last eight of the QRS complexes have pace pulse intervals within a narrow cluster centered at a long time prior to their respective QRS complexes, in which event the average dl of all these times is derived. If the condition does not exist, set OPl=0. For each shape group, a means 24 derives a paced template classification determined from the signals OPs, OPl, ds and dl whether the classification for the group should be dual-chamber paced, atrial chamber paced, ventricular chamber paced or not paced and stores it in the template buffer 12. Finally, a means 26 is provided for classifying the current QRS complex. If the last eight QRS complexes for a given shape including the current one causes the paced template classification means 24 to indicate a non paced condition, i.e. if OPs=0 and OPl=0, the classification must be made by some other means. FLOW CHARTS Pace Pulse to Beat Association Means 20 The purpose of this procedure is to store in the template storage buffer 12 the intervals between each pace pulse and the peak of the succeeding R wave that falls within a time window. In this example the time window extends from 48 to 600 msec before the the R wave peak. Every time a pace pulse is detected its time of occurrence is stored in the shift register 18, and every time a QRS is detected its time of occurrence is known. The procedure starts with block 28 that sets the index or pointer for the pace pulse shift register 18 at the end where the time of the latest pace pulse is stored. A block 32 then computes the difference PPR between the time of the current QRS and the time at time index in the shift register 18. It also computes the difference PPRP between the pace pulse time at time index and the time of the previous QRS complex. A decision block D1 then checks to see if PPR ≧ E1 where E1 is the shortest interval between a pace pulse and a QRS complex that will be used. In this embodiment E1 is 48 msec. If PPR is less than E1, D1 gives a false indication and the procedure goes to a block 34 which increases the time index by one. A decision block D2 then indicates whether the procedure has reached the end of the shift register 18. If not, the procedure is repeated until the output of D1 is true. When this occurs, D3 checks to see if PPRP > 0, then a decision block D4 checks to see if PPR ≦ E2, which in this embodiment is 600 msec. If the time of a pace pulse fetched from the shift register 18 is between 48 and 600 ms earlier than the time of the current QRS a block 36 stores the value of PPR in the template buffer 12 with the current QRS complex. The block 34 then increases the time index and the process is repeated until D2, D3 or D4 indicates a false condition. Pace Pulse Clustering Analyses Means 22 The purpose of this means is to examine the PPR values that have been stored in the template storage buffer 12 by the block 36 of Figure 5 for the last n QRS complexes having the same shape classification as the current complex to see if there are m PPR intervals differing by less than E3 seconds from each other. In one embodiment of this invention n=8, m=6 and E3=40 msec. If so, a cluster exists. A cluster of short intervals can be found by examining the PPR intervals in ascending order of magnitude as in Figure 6A and a cluster of long intervals can be found by examining the PPR intervals in descending order of magnitude as in Figure 6B. If a cluster of smaller intervals is found, a signal OPs=1 is produced and the average ds of these intervals is provided. If no such cluster is found, a signal OPs=0 is produced. If a cluster of larger intervals is found, a signal OPl=1 is produced and the average dl of these intervals is provided. If no such cluster is found, a signal OPl=0 is produced. The procedure of Figure 6A is as follows. A block 38 fetches the pace pulse PPR intervals stored in the buffer 12 for the last 8 QRS complexes having the same shape classification as the current QRS complex, and a block 40 arrange these intervals in ascending order in the buffer 12, and provides their total number. A block 42 initializes the procedure by setting an index1=0, and a block 44 increments index1 by one and sets an index2 to the same value. Decision block D5 checks to see if index1 exceeds the total number of PPR intervals. If it does, a block 45 sets OPS=0, but if not, a block 46 sets a counter equal to 1 and a summer equal to the first PPR interval in the buffer, which is at index1. A block 48 increments index2 by one so that it is at the next PPR interval beyond index1 and D6 checks to see if index2 exceeds the number of stored PPR intervals. If not, D7 checks to see if the difference between the PPR interval at index1 and the PPR interval at index2 is less than a given value E3, which in this embodiment is 40 msec. If so, a block 50 increments the counter by one and increments the summer by the value of the PPR interval at index2. This process can continue until D6 indicates that there are no more PPR intervals in the buffer, or if the difference between the PPR values at index2 and index1 is greater than E3. Decision block D8 checks to see if the counter is greater than or equal to 6. If not, the procedure reverts to block 44 that increments index1 by one and sets index2 equal to index1. Whenever D8 indicates a count of 6 or more, a block 52 set OPs=1 and derives ds by dividing the sum in the summer by the count in the counter. In Figure 6B components corresponding to Figure 6A have the same designations primed. The only difference is that the block 40′ arranges the PPR values in descending order. It, of course, provides the signals OPl=0, OPl=1 and dl. Paced Template Classification 24 The purpose of this procedure, Figure 7, is to utilize the information about clusters and average values thereof provided by the procedures of Figures 6A and 6B to store information in the template storage buffer 12 as to whether QRS complexes having the shape of the current QRS complex are paced or not and if so, whether they are dually paced, atrially paced or ventricularly paced. If D9 indicates that OPs=1 and that OPl=1, it means that there are two clusters so as to indicate the probability of dual chamber pacing, but a further check is made by D10 to see if the difference dl-ds between the average values of the clusters is less than a given amount herein indicated as being 120 ms. If the difference is not less than this amount, block 54 provides an indication to the template storage buffer 12 that the shape of the QRS complex just examined should be classified as dual chamber paced. On the other hand if D10 indicates that the time dl-ds is less than a given amount, the clusters are too close together for dual chamber pacing. D11 determines whether atrial pacing or ventricular pacing is involved by checking to see if dl is greater than a given time, herein indicated as being 152 msec. If it is, a block 56 stores in the template buffer 12 information, that atrial pacing is involved, but if it is not, a block 58 stores information in the buffer indicating that ventricular pacing is involved. Should D9 give a false indication meaning that both clusters do not exist at the same time D12,checks to see if a cluster of long intervals exists. If OPl=1, it does and the procedure goes to D11 again. But if OPl=0, D13 checks to see if a cluster of short intervals exists by checking OPs. If OPs=1 a cluster of short intervals does exist and D14 checks ds. If ds>152 msec, the single cluster indicates atrial pacing so that the procedure goes to the block 56, but if ds is less than or equal to 152 msec, a block 58 indicates that ventricular pacing is involved. If D13 indicates that OPs=0, no clusters exist so that pacing is not involved, and a block 60 stores this fact in the template buffer 12. QRS Classification 26 The purpose of this procedure, Figure 8, is to make a final classification of the current QRS complex. When this new complex arrived and was found by the QRS grouping means 10 to have a shape of one of the existing groups, the last eight complexes of that group with the current complex being the most recent one, are subjected to pace pulse clustering analysis 22, and paced template classification 24. D15 checks to see if the template classification is paced. If not, another classification is made, but if it is, D16 checks to see if the classification was dual-chamber paced. If so, D17 checks to see if the current QRS complex has a pace pulse within 48 msec of dl. If so the current QRS complex is classified as dual chamber paced by a block 62. If D17 indicates that the current QRS complex does not have a pace pulse within 48 msec of dl, D18 checks to see if it has a pace pulse within 48 msec of ds. If so, a block 64 indicates that the current complex is ventricularly paced. If not, D19 checks to see if the QRS complex is premature i.e. early compared to paced beat intervals. If so, a block 66 classifies the complex as ectopic. In the event that D16 indicates that the template classification is not dual chamber paced, D20 checks to see if the template classification was ventricularly paced. If it is the procedure goes to D18. If D20 indicates that the template classification is not ventricularly paced, D21 checks to see if there is a pace pulse within 48 msec from dl. If so, block 68 classifies the complex as being atrially paced. If not, block 70 classifies the complex as normal. If only one cluster is present, it can be detected by the procedure of Figure 6A and the procedure of Figure 6B so that the procedure shown in Figure 7 for providing a group classification must take this into account. D9 of Figure 7 will give a true signal and D10 will find that dl-ds is zero so that it will provide a true signal, but if two clusters have been detected that are separated by more than 120 msec, D10 will give a false signal and the block 54 will classify the group as being dual chamber paced. Figure 7 would not have to perform this function if the procedures for detecting clusters took into account the expected different intervals for atrial pace pulses and ventricular pace pulses.

Claims: Apparatus for classifying the type of pacing of QRS complexes contained in an ECG wave signal, comprising, first means (8) for detecting the time of occurrence of said QRS complexes in an ECG wave signal; second means (14,18) for detecting the times of occurrence of pace pulses in an ECG wave; third means (20) coupled to the second means for determining the intervals between each pace pulse and the next QRS complex, and grouping means (10) for assigning said QRS complexes to groups in accordance with their shape; characterised by clustering means (22) coupled to the third and the grouping means (20,10) for detecting the existence of at least one cluster of said intervals occurring during a given number of the most recent QRS complexes; means for deriving an average interval value for each said cluster, and classifying means (26) for classifying the current QRS complex by the correlation between said intervals associated with it and the averages for the clusters of the group to which it has been assigned. Apparatus as set forth in claim 1 wherein said classifying means (26) classifies a current QRS complex as dual chamber paced when two of said averages are separated by at least a given amount and the interval of a pace pulse associated with the current QRS complex is within a given time of the cluster having the larger interval. Apparatus as set forth in claim 1 wherein said classifying means (26) classifies a current QRS complex as as ventricularly paced when said clustering means (22) derives an average of only one cluster that is less than a given amount and the current QRS complex has a pace pulse associated with it that is within a given time of said average. Apparatus as set forth in Claim 1 wherein said classifying means (26) classifies a current QRS complex as atrially paced when said clustering means derives an average of only one cluster that is greater than said given amount and the interval of a pace pulse associated with the current QRS complex is within a given time of the cluster having the larger interval. Apparatus as set forth in claim 1 wherein said classifying means classifies a current QRS complex as normal when said clustering means has detected only one cluster for the group to which the current QRS complex has been assigned that has an average interval greater than a given amount of time and any interval of the current QRS complex differs from the average by at least a different given amount of time. Apparatus as set forth in claim 1 wherein said classifying means classifies a current QRS complex as ectopic when said clustering means has detected two clusters for the group to which the current QRS complex has been assigned that have respective averages differing by at least a given amount of time and any intervals of the current QRS complex differ from both averages by at least a different given amount of time and the current QRS complex is premature in time.

Inventors: SHAYA, MOUSA NOORY, WANG, JYH-YUN
Applicants: HEWLETT-PACKARD COMPANY
Application year: 1988
Main classification: A61B5
IPC codes: A61B5, A61N1