Method and system for determining respiration rate

• US 10,420,516 B2
• Filed: 09/15/2014
• Issued: 09/24/2019
• Est. Priority Date: 11/20/2013
• Status: Active Grant

• Current Assignee: General Electric Company
• Original Assignee: General Electric Company
• Inventors: Swamy, Gokul, Genc, Sahika, Rao, Ashish Anil, Patil, Abhijit Vishwas, Anandkumar, Amod Jai Ganesh
• Primary Examiner(s): Luan, Scott
• Application Number:
  US15/035,632
• Publication Number:
  US 20160262705A1
• Time in Patent Office: 1,835 Days
• CPC Class Codes: A61B 5/0205 Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
• CPC Class Codes: A61B 5/021 Measuring pressure in heart or blood vessels A61B5/0205 takes precedence
• CPC Class Codes: A61B 5/0535 Impedance plethysmography for measuring blood flow A61B5/0295
• CPC Class Codes: A61B 5/0816 Measuring devices for examining respiratory frequency measuring frequency of electric signals G01R23/00
• CPC Class Codes: A61B 5/318 Heart-related electrical modalities, e.g. electrocardiography [ECG]
• CPC Class Codes: A61B 5/7207 of noise induced by motion artifacts
• CPC Class Codes: A61B 5/7225 Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation A61B5/30, A61B5/316 take precedence; transformers for specific applications H01F38/00; amplifiers H03F
• CPC Class Codes: A61B 5/725 using specific filters therefor, e.g. Kalman or adaptive filters A61B5/316 take precedence
• CPC Class Codes: A61B 5/7278 Artificial waveform generation or derivation, e.g. synthesising signals from measured signals A61B5/327 takes precedence

Abstract

A method for determining a respiration rate of a subject, includes receiving a first signal and a second signal, each signal being representative of a physiological parameter of the subject. The method includes removing a cardiac artifact signal from the first signal and the second signal to generate a first processed signal and a second processed signal respectively. The method includes removing a motion artifact signal from the first processed signal and the second processed signal to generate a first periodic signal and the second processed signal respectively. The method further includes removing a residual noise signal from the first periodic signal and the second periodic signal to generate a first noise free signal and the second noise free signal respectively. The method includes generating a combined value from a first value and a second value based on the first noise free signal and the second noise free signal respectively.

Claims

1. A method, comprising:

receiving a plurality of signals comprising a first signal and a second signal, each signal being representative of physiological parameters of a subject;

(i) processing the first signal to remove a cardiac artifact signal from the first signal and generate a first processed signal and processing the second signal to remove the cardiac artifact signal from the second signal and generate a second processed signal;

(ii) processing the first processed signal using a polynomial fitting technique to remove a motion artifact signal from the first processed signal and generate a first periodic signal and processing the second processed signal using the polynomial fitting technique to remove the motion artifact signal from the second processed signal and generate a second periodic signal;

(iii) processing the first periodic signal using a spectral decomposition technique to remove a residual noise signal from the first periodic signal and generate a first noise free signal and processing the second periodic signal using the spectral decomposition technique to remove the residual noise signal from the second periodic signal and generate a second noise free signal;

(iv) processing the first noise free signal using a subspace-based MUSIC technique or ESPRIT technique to generate a first value and processing the second noise free signal using the subspace-based technique to generate a second value; and

determining a combined value based on the first value and the second value, wherein the combined value is representative of a respiration rate of the subject.

2. The method of claim 1, wherein the plurality of signals further comprises a third signal, wherein the first signal is an impedance plethysmograph signal, the second signal is an electro cardiograph signal, and the third signal is an arterial blood pressure signal.

3. The method of claim 2, further comprising:

processing the third signal to remove the cardiac artifact signal from the third signal and generate a third processed signal;

processing the third processed signal using the polynomial fitting technique to remove the motion artifact signal from the third processed signal and generate a third periodic signal;

processing the third periodic signal using the spectral decomposition technique to remove the residual noise signal from the third periodic signal and generate a third noise free signal;

processing the third noise free signal using the subspace-based technique to generate a third value; and

determining the combined value based on the first value, the second value, and the third value.

4. The method of claim 1, wherein determining the combined value comprises:

determining a first weighted value for the first value based on a first signal-to-noise ratio of the first signal;

determining a second weighted value for the second value based on a second signal-to-noise ratio of the second signal; and

determining a weighted average value based on the first weighted value, the first value, the second weighted value, and the second value.

5. The method of claim 1, wherein the polynomial fitting technique comprises a linear least squares method.

6. The method of claim 1, wherein the spectral decomposition technique comprises determining a spectral component corresponding to at least one of the residual noise signal and the first noise free signal, in the first periodic signal.

7. The method of claim 6, wherein determining the spectral component comprises:

constructing a time series matrix based on a plurality of sample values of the first periodic signal; and

performing a singular value decomposition of the time series matrix.

8. The method of claim 6, wherein determining the spectral component comprises determining a principal mode of the first periodic signal, using an empirical mode decomposition (EMD) technique.

9. A system comprising:

a signal acquisition system configured to receive a plurality of signals comprising a first signal and a second signal, each signal being representative of physiological parameters of a subject; and

at least one processor coupled to the signal acquisition system, the at least one processor encoded with instructions to perform signal processing of the first and second signals, the instructions including;

a cardiac-artifact removal module to process each signal among the first signal and the second signal to remove a cardiac artifact signal and generate a first processed signal and a second processed signal respectively;

a motion-artifact removal module to process each signal among the first processed signal and the second processed signal, using a polynomial fitting technique, to remove a motion artifact signal and generate a first periodic signal and a second periodic signal respectively;

a noise removal module to process each signal among the first periodic signal and the second periodic signal, using a spectral decomposition technique, to remove a residual noise signal and generate a first noise free signal and a second noise free signal respectively;

an estimation module to process each signal among the first noise free signal and the second noise free signal, using a subspace-based MUSIC technique or ESPRIT technique, and generate a first value and a second value respectively; and

a combining module to determine a combined value based on the first value and the second value, wherein the combined value is representative of a respiration rate of the subject.

12. The system of claim 9, wherein the combining module determines the combined value by:

determining a first weighted value for the first value based on a first signal-to-noise ratio of the first signal;

determining a second weighted value for the second value based on a second signal-to-noise ratio of the second signal; and

determining a weighted average value based on the first weighted value, the first value, the second weighted value, and the second value.

13. The system of claim 9, wherein the cardiac-artifact removal module comprises a notch filter.

14. The system of claim 9, wherein the motion-artifact removal module comprises a Savitzky Golay filter.

15. The system of claim 9, wherein the noise removal module uses the spectral decomposition technique comprising an empirical mode decomposition (EMD) technique.

16. The system of claim 9, wherein the noise removal module uses the spectral decomposition technique comprising a singular spectrum analysis (SSA) technique.

10. The system of 9, wherein the signal acquisition system is configured to receive a third signal, wherein the first signal is an impedance plethysmograph signal, the second signal is an electro cardiograph signal, and the third signal is an arterial blood pressure signal.

11. The system of claim 10, wherein the cardiac-artifact removal module, the motion-artifact removal module, the noise removal module, and the estimation module generates a third processed signal, a third periodic signal, a third noise free signal, and a third value respectively, and determines the combined value based on the first value, the second value, and the third value.

17. A non-transitory computer readable medium encoded with a program to instruct at least one processor to:

receive a plurality of signals comprising a first signal and a second signal, each signal being representative of physiological parameters of a subject;

(i) process the first signal to remove a cardiac artifact signal from the first signal and generate a first processed signal and process the second signal to remove the cardiac artifact signal from the second signal and generate a second processed signal;

(ii) process the first processed signal using a polynomial fitting technique to remove a motion artifact signal from the first processed signal and generate a first periodic signal and process the second processed signal using the polynomial fitting technique to remove the motion artifact signal from the second processed signal and generate a second periodic signal;

(iii) process the first periodic signal using a spectral decomposition technique to remove a residual noise signal from the first periodic signal and generate a first noise free signal and process the second periodic signal using the spectral decomposition technique to remove the residual noise signal from the second periodic signal and generate a second noise free signal;

(iv) process the first noise free signal using a subspace-based MUSIC technique or ESPRIT technique to generate a first value and process the second noise free signal using the subspace-based technique to generate a second value; and

determine a combined value based on the first value and the second value, wherein the combined value is representative of a respiration rate of the subject.