Method and apparatus for reducing motion artifacts in ECG signals

• US 9,826,937 B2
• Filed: 12/13/2013
• Issued: 11/28/2017
• Est. Priority Date: 12/31/2012
• Status: Active Grant

• Current Assignee: Koninklijke Philips N.V.
• Original Assignee: Koninklijke Philips N.V.
• Inventors: Wang, Jin, Zhao, Dan, Shi, Cheng, Chen, Shiyang
• Primary Examiner(s): Carey, Michael
• Assistant Examiner(s): Patel, Natasha
• Application Number:
  US14/654,570
• Publication Number:
  US 20150342534A1
• Time in Patent Office: 1,446 Days
• CPC Class Codes: A61B 5/316 Modalities, i.e. specific diagnostic methods
• CPC Class Codes: A61B 5/349 Detecting specific parameters of the electrocardiograph cycle
• CPC Class Codes: A61B 5/7207 of noise induced by motion artifacts

Abstract

The present invention provides a method and apparatus for reducing motion artifacts in ECG signals. According to an aspect of the present invention, there is proposed a method of reducing motion artifacts in ECG signals, comprising: acquiring a current beat from a continuously measured ECG signal of a patient; calculating a correlation coefficient between a previous mean value beat and the current beat in the ECG signal; determining the weights to be assigned to the previous mean value beat and the current beat based on the correlation coefficient; and calculating a current mean value beat based on the previous mean value beat, the current beat, and the weights thereof. Accordingly, the novel method of deriving the current mean value beat may reduce ECG artifacts due to patient movement in such a manner that the SNR of the ECG signal can be improved substantially.

Claims

1. A method of reducing motion artifacts in ECG signals, comprising:

acquiring, via an acquiring unit, a current beat from a continuously measured ECG signal of a patient;

calculating, via a first calculating unit, a correlation coefficient between a previous mean value beat and the current beat in the ECG signal;

determining, via a determining unit, weights to be assigned to the previous mean value beat and the current beat based on the correlation coefficient; and

calculating, via a second calculating unit, a current mean value beat based on (i) the previous mean value beat, (ii) the current beat, and (iii) the determined weights respectively assigned to the previous mean value beat and the current beat, wherein the current mean value beat comprises a motion artifact reduced beat having reduced motion artifacts over the current beat in the ECG signal.

2. The method according to claim 1, further comprising a step of comparing, via a comparator, the correlation coefficient with a predetermined coefficient.

3. The method according to claim 1, further comprising:

subtracting, via a subtracting unit, the current mean value beat from the current beat to obtain a residual signal;

performing, via a filter unit, piecewise filtering on the residual signal based on signal characteristics of the ECG signal to obtain a filtered signal; and

adding, via an adding unit, the filtered signal back to the current mean value beat to obtain a modified current beat.

4. The method according to claim 3, further comprising:

extracting, via an extracting unit, from the ECG signal, the signal characteristics, wherein the signal characteristics include (i) a noise envelope, (ii) correlation information, (iii) peak positions of P wave, QRS complex wave, and T wave, and (iv) segment intervals, and

selecting, via a selection unit, (i) one of P wave, QRS complex wave, and T wave, or (ii) any combination thereof, as the current beat to be processed.

5. The method according to claim 4, wherein, responsive to selecting the QRS complex wave as the current beat to be processed, the step of performing piecewise filtering further comprises the sub-steps of:

identifying, via an identification unit, a QRS segment and T and P segments based on the signal characteristics; and

performing filtering on the QRS, T and P segments, respectively, via different filters.

6. The method according to claim 1, further comprising

deriving, via a derivation unit, an initial mean value beat from a predetermined number of consecutive beats in the ECG signal.

7. The method according to claim 6, wherein the predetermined number of consecutive beats comprises (i) first predetermined number of consecutive beats and (ii) a second predetermined number of consecutive beats following the first predetermined number of consecutive beats, and wherein the step of deriving an initial mean value beat further comprises the sub-steps of:

deriving an initial mean value beat from the first predetermined number of consecutive beats; and

updating the initial mean value beat only in response to the correlation coefficient between a previous initial mean value beat and a current beat in the second predetermined number of consecutive beats being larger than a predetermined coefficient.

8. The method according to claim 6, further comprising:

reactivating, in response to a user operation of a reactivation button, the step of deriving an initial mean value beat from a predetermined number of consecutive beats in the ECG signal.

9. The method according to claim 1, wherein the steps of acquiring via the acquiring unit, calculating via the first calculating unit, determining via the determining unit, and calculating via the second calculating unit are applied to single channel ECG signals.

10. An apparatus for reducing motion artifacts in ECG signals, comprising:

an acquiring unit for acquiring a current beat from a continuously measured ECG signal of a patient;

a first calculating unit for calculating a correlation coefficient between a previous mean value beat and the current beat in the ECG signal;

a determining unit for determining the weights to be assigned to the previous mean value beat and the current beat based on the correlation coefficient; and

a second calculating unit for calculating a current mean value beat based on (i) the previous mean value beat, (ii) the current beat, and (iii) the determined weights respectively assigned to the previous mean value beat and the current beat, wherein the current mean value beat comprises a motion artifact reduced beat having reduced motion artifacts over the current beat in the ECG signal.

11. The apparatus according to claim 10, further comprising a comparator for comparing the correlation coefficient with a predetermined coefficient.

12. The apparatus according to claim 10, further comprising:

a subtracting unit for subtracting the current mean value beat from the current beat to obtain a residual signal;

a filter unit for performing piecewise filtering on the residual signal based on signal characteristics of the ECG signal to obtain a filtered signal; and

an adding unit for adding the filtered signal back to the current mean value beat to obtain a modified current beat.

13. The apparatus according to claim 12, further comprising:

an extracting unit for extracting, from the ECG signal, the signal characteristics, wherein the signal characteristics include (i) a noise envelope, (ii) correlation information, (iii) peak positions of P wave, QRS complex wave, and T wave, and (iv) segment intervals; and

a selection unit for selecting (i) one of P wave, QRS complex wave, and T wave, or (ii) any combination thereof, as the current beat to be processed.

14. The apparatus according to claim 13, wherein responsive to the selection unit selecting the QRS complex wave as the current beat to be processed, the apparatus further comprising:

an identifying unit for identifying a QRS segment and T and P segments based on the signal characteristics, and wherein the filtering unit comprises a first group of low-pass filters and high-pass filters for the QRS segment and a second group of low-pass filters and high-pass filters for T and P segments, respectively.

15. The apparatus according to claim 14, wherein the first group of low-pass filters and high-pass filters for the QRS segment and the second group of low-pass filters and high-pass filters are different in cutoff frequency.