Method and apparatus for signal averaging and analyzing high resolution P wave signals from an electrocardiogram

• US 5,840,038 A
• Filed: 05/29/1997
• Issued: 11/24/1998
• Est. Priority Date: 05/29/1997
• Status: Expired due to Term

• Current Assignee: Marquette Medical Systems, Inc.
• Original Assignee: Marquette Medical Systems, Inc.
• Inventors: Xue, Qiuzhen, Reddy, Shankara
• Primary Examiner(s): Getzow, Scott M.
• Application Number:
  US08/865,155
• Time in Patent Office: 544 Days
• Field of Search: 600/508, 600/509, 600/510, 600/512, 600/515, 600/516, 600/517, 600/521
• US Class Current: 600/512
• CPC Class Codes: A61B 5/349 Detecting specific parameters of the electrocardiograph cycle
• CPC Class Codes: A61B 5/353 Detecting P-waves

Abstract

A method and apparatus for signal averaging and analyzing high resolution P wave signals from an electrocardiogram. The method and apparatus for averaging the P wave portion of ECG signals includes sensing human ECG signals having at least a QRS portion and a P wave portion, detecting the QRS portion of the ECG signals, suppressing the QRS portion of the ECG signals, generating a P wave detection function from the wave portion of the ECG signals, cross correlating the P wave portions of the ECG signals and averaging the cross correlated P wave signals. The method and apparatus for analyzing high resolution P wave signals includes averaging at least a portion of the ECG signals, filtering the averaged signals, deriving the vector magnitude (VM) of the filtered averaged signals, determining the base noise level, and setting the P onset and P offset of the VM signals above the noise level.

Claims

1. A method of averaging the P wave portion of ECG signals comprising the steps of,

sensing human ECG signals having at least a QRS portion and a P wave portion,

detecting the QRS portion of the ECG signals,

suppressing the QRS portion of the ECG signals,

generating a P wave detection function from the P wave portion of the ECG signals,

cross correlating the detection functions of the P wave portions of the ECG signals, and

averaging the cross correlated P wave detection function signals.

2. The method set forth in claim 1 and including the step of generating a P wave detection function template and cross the correlating the detection function of the P wave portions of the ECG signals with the template.

3. The method set forth in claim 2 wherein the template is generated by sampling a plurality of successive ECG wave forms, detecting the QRS portions of the ECG signals, suppressing the QRS portions of the ECG signals, generating a P detection function from the P wave portions of the ECG signals, selecting a seed beat from the P wave detection functions of the ECG signals, cross correlating the detection functions of the P wave portions of the ECG signals with the seed beat, averaging the cross correlated signals to provide the template.

4. The method set forth in claim 1 wherein the step of detecting the QRS portion of the ECG signals comprised defining the QRS portion and the step of suppressing the QRS portion of the ECG signals comprises subtracting the QRS portion of the signal from the ECG signal.

5. The method set forth in claim 1 wherein the ECG signals sensed include X, Y and Z ECG lead signals.

6. The method set forth in claim 1 wherein the step of generating the P detection function comprises filtering the ECG signal after suppression of the QRS portion, generating the vector magnitude of the P portions of X, Y and Z ECG lead signals, and low pass filtering the resultant signals.

7. The method set forth in claim 1 and including the steps of defining a P wave search zone, locating the signal peaks in the search zone, selecting the maximum peak, defining a P wave threshold and determining P wave on and P wave off markers within the search zone.

8. The method set forth in claim 1 wherein the step of cross correlating the P waves comprises the steps of computing the P-Q interval of the ECG signals, setting flags at the maximum and minimum levels of acceptable signals, and finding the maximum cross correlation co-efficient of the signals.

9. The method set forth in claim 1 wherein the ECG signals sensed include X, Y and Z lead signals, and including the step of generating a template of the P wave detection functions and cross correlating the detection functions of the P wave portions of the ECG signals with the template.

10. The method set forth in claim 9 and including the step of defining a P wave search zone, selecting the maximum peak in the search zone, and defining a P wave threshold in determining P wave on and off markers within the search zone.

11. The method set forth in claim 9 wherein the step of cross correlating the P waves comprises the steps of computing the P-Q interval of the ECG signals, setting flags at the maximum minimum levels of acceptable signals, and finding maximum cross correlation coefficient of the signals.

12. The method set forth in claim 11 wherein the step of generating the P detection function comprises filtering the ECG signals after suppression of the QRS portion, generating the vector magnitude of the P portion of the X, Y and Z ECG lead signals, and low pass filtering the resultant signals.

13. The method set forth in claim 12 wherein the step of detecting the QRS portion of the ECG signals comprises defining the QRS portion and the step of suppressing the QRS portion of ECG signals comprises subtracting the QRS portion of the signals from the ECG signals.

14. An apparatus for averaging the P wave portions of ECG signals and including, acquisition means for acquiring human ECG signals taken along the X, Y and Z axes and having at least a QRS portion and a P wave portion, means for detecting the QRS portions of the ECG signals, means for suppressing the QRS portion of the ECG signals, means for generating a P detection function of the ECG signals after suppression of the QRS portion, means for cross correlating the detection functions of the P wave portion of the ECG signals with a template and means for averaging the cross correlated detection functions of the P wave portions.

15. The apparatus set forth in claim 14 including means for defining a P wave search zone, means for locating the maximum signal peak in the search zone and means for defining a P wave threshold and for determining P wave on and P wave off markers within the search zone.

16. An apparatus for averaging the P wave portions of ECG signals and including,

acquisition means for acquiring human ECG signals taken along the X, Y and Z axes and having at least a QRS portion and a P wave portion,

processor means programmed for detecting the QRS portions of the ECG signals, for suppressing the QRS portions of the ECG signals, for generating a P detection function, for cross correlating the detection functions of the P waves with a template, and for averaging the cross correlated detection functions of the P waves.

17. A method of determining the filtered onset and offset of the P wave portion of ECG signals, comprising the steps of averaging at least a portion of the ECG signals, filtering the averaged signals, deriving the vector magnitude (VM) of the filtered averaged signals, determining the base noise level, and setting the P onset and P offset of the VM signals above the noise level.

18. The method set forth in claim 17 including the steps of setting a P offset threshold and a P onset threshold, setting the P offset by locating a first sample point of the vector magnitude signal at which a predetermined number of successive sample points from the first sample point are less than the initial P offset, and setting the P onset threshold by locating a second sample point of the vector magnitude signal which is a predetermined number of successive sample points forwardly from the second sample point are less than the P onset threshold.

19. The method set forth in claim 18 and including the step of low pass filtering the ECG signal portion prior to determining the vector magnitude of said signals.

20. The method set forth in claim 19 wherein said ECG signals are taken along the X, Y and Z axes.

21. The method set forth in claim 18 and including the step of band pass filtering the ECG signals prior to determining the vector magnitude of said signals.

22. The method set forth in claim 21 wherein said ECG signals are taken along the X, Y and Z axes.

23. An apparatus for determining the onset and offset of the P wave portion of ECG signals, including means for averaging at least a portion of the ECG signals, means for filtering the average signals, means for deriving the vector magnitude of the filter average signals, means for determining the base noise level, and means for setting the P onset and P offset of the vector magnitude signals above the noise level.

24. An apparatus for averaging the P wave portions of ECG signals and including,

acquisition means for acquiring human ECG signals,

processor means programmed to detect the QRS portions of the ECG signals, to suppress the QRS portions of the ECG signals, to generate a P detection function, to derive a template of the P wave detection functions, to cross correlate the P wave with the template, and to average the cross correlated P waves.

25. The apparatus set forth in claim 24 wherein said processor is further programmed to generate a P wave detection function template and to cross the correlating the detection function of the P wave portions of the ECG signals with the template.

26. The apparatus set forth in claim 24 wherein said processor is further programmed to detect the QRS portion of the ECG signals by defining the QRS portion and to suppress the QRS portion of the ECG signals by subtracting the QRS portion of the signal from the ECG signal.

27. The apparatus set forth in claim 24 wherein said processor is further programmed to generate the P detection function by filtering the ECG signal after suppression of the QRS portion, to generate the vector magnitude of the P wave detection functions portions of X, Y and Z ECG lead signals, and to low pass filter the resultant signals.

28. The apparatus set forth in claim 24 wherein said processor is further programmed to define a P wave search zone, to locate the signal peaks in the search zone, to select the maximum peak, to define a P wave threshold and to determine P wave on and P wave off markers within the search zone.

29. The apparatus set forth in claim 24 wherein said processor is further programmed to cross correlate the P waves by computing the P-Q interval of the ECG signals, to set flags at the maximum and minimum levels of acceptable signals, and to find the maximum cross correlation co-efficient of the signals.