Separating motion from cardiac signals using second order derivative of the photo-plethysmogram and fast fourier transforms

• US 8,073,516 B2
• Filed: 01/30/2006
• Issued: 12/06/2011
• Est. Priority Date: 01/31/2002
• Status: Active Grant

• Current Assignee: Spacelabs Healthcare LLC (OSI Systems Incorporated)
• Original Assignee: Dolphin Medical, Inc. (OSI Systems Incorporated)
• Inventors: Scharf, John E., Shah, Bhavin B.
• Primary Examiner(s): WINAKUR, ERIC FRANK
• Application Number:
  US11/343,914
• Publication Number:
  US 20060211930A1
• Time in Patent Office: 2,136 Days
• Field of Search: 600/310, 600/322, 600/323, 600/330, 600/336, 600/481, 600/500
• US Class Current: 600/336
• CPC Class Codes: A61B 5/14551 for measuring blood gases
• CPC Class Codes: A61B 5/7207 of noise induced by motion artifacts
• CPC Class Codes: A61B 5/7239 using differentiation including higher order derivatives
• CPC Class Codes: A61B 5/7253 characterised by using transforms
• CPC Class Codes: A61B 5/7257 using Fourier transforms
• CPC Class Codes: A61B 5/726 using Wavelet transforms

Abstract

The present invention is directed toward a pulse oximetry system for the determination of a physiological parameter capable of removing motion artifacts from physiological signals comprises a hardware subsystem and a software subsystem. The software subsystem is used in conjunction with the hardware subsystem to perform a method for removing a plurality of motion artifacts from the photo-plethysmographic data and for obtaining a measure of at least one physiological parameter from the data. The method comprises acquiring the raw photo-plethysmographic data, transforming the data into the frequency domain, analyzing the transformed data to locate a series of candidate cardiac spectral peaks (primary plus harmonics), reconstructing a photo-plethysmographic signal in the time domain with only the candidate cardiac spectral peaks (primary plus harmonics), computing the second order derivative of the reconstructed photo-plethysmographic signal, analyzing the candidate second order derivative photo-plethysmographic signal to determine the absence or presence of cardiac physiologic signal characteristics, and finally selecting the best physiologic candidate from the series of potential cardiac spectral peaks (primary plus harmonics) based upon a second derivative scoring system. This scoring system is preferentially based upon second derivative processing analysis, but can be equally applied using the first, third, fourth or other similar derivative processing analysis.

Claims

1. A pulse oximetry system for the determination of a physiological parameter wherein said pulse oximetry system is configured to remove motion artifacts from physiological signals, comprising:

an optical sensor for providing said pulse oximetry system with photo-plethysmographic data, wherein said sensor comprises a light emitting diode for emitting light and a photo-detector for detecting said light, wherein said detected light contains information on the photo-plethysmographic data;

a processor for processing the photo-plethysmographic data; and

a memory unit, wherein said memory unit is operably coupled to said processor, and wherein said processor is configured to;

(a) acquire said photo-plethysmographic data;

(b) transform the photo-plethysmographic data into frequency domain data;

(c) analyze the frequency domain data;

(d) reconstruct a time domain data photo-plethysmographic signal from said frequency domain data;

(e) analyze a second order derivative of the time domain photo-plethysmographic signal to identify parameters;

(f) use said parameters to separate a physiological signal from a motion artifact in the time domain photo-plethysmographic signal; and

(g) compute a physiological parameter from said motion artifact independent physiological signal.

2. The system of claim 1 wherein the processor is configured to condition the photo-plethysmographic data prior to transforming the data into the frequency domain.

3. The system of claim 2 wherein conditioning said photo-plethysmographic data reduces spectral leakage from subsequent signal processing and frequency analysis.

4. The system of claim 1 wherein the processor is configured to analyze the frequency domain data to locate candidate spectral peaks before transforming said frequency domain data into time domain photo-plethysmographic signal data, wherein analyzing said transformed data to locate candidate spectral peaks comprises the steps of assigning a largest power amplitude from said data as a primary candidate spectral peak and assigning a next largest power amplitude that is not a harmonic of said primary candidate spectral peak as a secondary candidate spectral peak.

5. The system of claim 1 wherein the processor is configured to analyze the frequency domain data to locate candidate spectral peaks before transforming said frequency domain data into time domain photo-plethysmographic signal data, wherein said analyzing transformed data to locate candidate spectral peaks is performed using a first derivative criteria, wherein said candidate spectral peak is a frequency bin having a slope value equal to approximately zero.

6. The system of claim 1 wherein the processor is configured to compute a second order derivative of the time domain photo-plethysmographic signal after reconstructing the time domain data photo-plethysmographic signal.

7. The system of claim 6 wherein the second order derivative photo-plethysmographic signal is analyzed to determine the absence or presence of cardiac physiologic signal characteristics.

8. The system of claim 7 wherein analyzing said second order derivative photo-plethysmographic signal comprises the step of performing a morphological analysis of the second order derivative photo-plethysmographic signal to determine a unique cardiac morphology, wherein said cardiac morphology consists of a plurality of systolic and a plurality of diastolic waves, and wherein said analysis accounts for a plurality of wave transitions.

9. The system of claim 8 wherein analyzing said second order derivative photo-plethysmographic signal further comprises the step of performing a morphological analysis of the second order derivative photo-plethysmographic signal to determine at least one height of said wave transitions.

10. The system of claim 6 wherein analyzing said second order derivative photo-plethysmographic signal comprises the steps of:

determining a plurality of wave transitions, calculating a plurality of heights of said wave transitions, and applying a criterion for detecting a predefined type of wave transition.

11. The system of claim 10 wherein the criterion for detecting a predefined type of wave transition among a first wave, a second wave, a third wave, a fourth wave, and a fifth wave comprises identifying the largest transition between a first wave and a second wave and determining whether a set of first, second, third, fourth, and fifth waves adhere to a formula.

12. The system of claim 1 wherein the time domain photo-plethysmographic data is reconstructed by including only a plurality of primary candidate spectral peaks and a plurality of corresponding harmonic spectral peaks.

13. The system of claim 1 wherein using said parameters to separate a physiological signal from a motion artifact comprises removing motion artifacts from the photo-plethysmographic signal by classifying a first primary-harmonic set with a first highest score as a likely cardiac geometry and classifying a second primary-harmonic set with a highest score as a likely motion geometry.

14. A pulse oximetry system for the determination of a physiological parameter wherein said pulse oximetry system is configured to remove motion artifacts from physiological signals, comprising:

an optical sensor for providing said pulse oximetry system with photo-plethysmographic data, wherein said sensor comprises a light emitting diode for emitting light and a photo-detector for detecting said light, wherein said detected light contains information on the photo-plethysmographic data;

a processor for processing the photo-plethysmographic data; and

a memory unit, wherein said memory unit is operably coupled to said processor, and wherein said processor is configured to;

(a) analyze frequency domain data derived from said photo-plethysmographic data, (b) transform the frequency domain data into a time domain photo-plethysmographic signal;

(c) analyze a second order derivative of the time domain photo-plethysmographic signal to identify parameters; and

(d) use the identified parameters to separate a physiological signal from a motion artifact in the time domain photo-plethysmographic signal.