Blind source separation of pulse oximetry signals

• US 6,701,170 B2
• Filed: 11/02/2001
• Issued: 03/02/2004
• Est. Priority Date: 11/02/2001
• Status: Expired due to Term

• Current Assignee: Covidien LP (Medtronic Plc)
• Original Assignee: Nellcor Puritan Bennett Incorporated (Medtronic Plc)
• Inventors: Stetson, Paul F.
• Primary Examiner(s): Winakur, Eric F.
• Assistant Examiner(s): Kremer, Matthew
• Application Number:
  US10/033,703
• Publication Number:
  US 20030088164A1
• Time in Patent Office: 851 Days
• Field of Search: 600/309-310, 600/300, 600/322-326, 600/709, 600/481, 600/500-502
• US Class Current: 600/323
• CPC Class Codes: A61B 5/02416 using photoplethysmograph signals, e.g. generated by infrared radiation A61B5/14552 takes precedence
• 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

Abstract

A method and apparatus for the application of Blind Source Separation (BSS), specifically independent Component Analysis (ICA) to mixture signals obtained by a pulse oximeter sensor. In pulse oximetry, the signals measured at different wavelengths represent the mixture signals, while the plethysmographic signal, motion artifact, respiratory artifact and instrumental noise represent the source components. The BSS is carried out by a two-step method including an ICA. In the first step, the method uses Principal Component Analysis (PCA) as a preprocessing step, and the Principal Components are then used to derive sat and the Independent Components, where the Independent Components are determined in a second step. In one embodiment, the independent components are obtained by high-order decorrelation of the principal components, achieved by maximizing the sum of the squares of the higher-order cumulants of the plurality of mixture signals.

Claims

1. A method for measuring a physiological parameter, comprising:

measuring a plurality of signals, wherein each of said signals comprises a source component corresponding to said physiological parameter and an interference component;

processing said plurality of signals to obtain a plurality of principal components;

processing said plurality of principal components to obtain a plurality of independent components, wherein a matrix of said plurality of signals corresponds to a matrix product of a matrix of said plurality of independent components and a matrix of mixing coefficients; and

extracting a first measure of said physiological parameter corresponding to said source component from one of said plurality of independent components, wherein said plurality of signals corresponds to sensed optical energies from a plurality of wavelengths.

2. The method of claim 1 wherein said physiological parameter is a function of an oxygen saturation.

3. The method of claim 1 wherein said processing said plurality of signals further comprises

obtaining a time derivative of the sensed optical energies from a plurality of wavelengths.

4. The method of claim 1 wherein said interference component comprises signal components caused by motion, respiratory artifact, ambient light, optical scattering and other interference between a tissue location being sensed and a sensor.

5. The method of claim 1 wherein said processing said plurality of signals further comprises decorrelating said plurality of signals by minimizing a cross-correlation of said plurality of signals, to obtain a plurality of decorrelated signals;

and

normalizing said plurality of decorrelated signals to obtain the plurality of principal components.

6. The method of claim 1 wherein said processing said plurality of signals comprises decorrelating said plurality of signals by singular-value decomposition of said plurality of signals, to obtain the plurality of principal components.

7. The method of claim 1 wherein said processing said plurality of signals comprises decorrelating said plurality of signals by multiplying said plurality of signals the inverse square root of the covariance matrix of said plurality of signals to obtain the plurality of principal components.

8. The method of claim 1 wherein said processing of said plurality of principal components comprises higher-order decorrelation of said plurality of principal components.

9.The method of claim 1 wherein said processing said plurality of principal components comprises maximizing a function of the higher-order cumulants of a mixture of said plurality of signals, thus separating said source component from said interference component.

12. The method of claim 1 further comprising obtaining a ratio of mixing coefficients from said matrix of mixing coefficients, wherein said ratio corresponds to a ratio of modulation ratios of red to infrared signals, wherein said plurality of signals comprise modulated optical signals in the red and infrared ranges.

14. The method of claim 1 further comprising extracting said interference component from another one of said plurality of independent components.

9. The method of claim 9 wherein said higher-order cumulant is cumulant having order greater than two.

10. The method of claim 9 wherein said higher-order cumulant is a third-order cumulant of said plurality of signals.

11. The method of claim 9 wherein said higher-order cumulant is a fourth-order cumulant of said plurality of signals.

13. The method of claim 13 further comprising extracting a second measure of said physiological parameter from said ratio, wherein said second measure of said physiological parameter corresponds to an oxygen saturation.

15. A method for measuring a physiological parameter, comprising:

measuring a plurality of signals, wherein each of said signals comprises a source component corresponding to said physiological parameter and an interference component;

processing said plurality of signals to obtain a plurality of principal components;

processing said plurality of principal components to obtain a plurality of independent components, wherein a matrix of said plurality of signals corresponds to a matrix product of a matrix of said plurality of independent components and a matrix of mixing coefficients; and

extracting a first measure of said physiological parameter corresponding to said source component from one of said plurality of independent components, wherein said physiological parameter is a function of a pulse rate.

16. A method for measuring a physiological parameter, comprising:

measuring a plurality of signals, wherein each of said signals comprises a source component corresponding to said physiological parameter and an interference component;

processing said plurality of signals to obtain a plurality of principal components;

processing said plurality of principal components to obtain a plurality of independent components, wherein a matrix of said plurality of signals corresponds to a matrix product of a matrix of said plurality of independent components and a matrix of mixing coefficients; and

extracting a first measure of said physiological parameter corresponding to said source component from one of said plurality of independent components, wherein said plurality of signals corresponds to sensed optical energies from a plurality of wavelengths from different times.

17. A method for measuring a physiological parameter, comprising:

measuring a plurality of signals, wherein each of said signals comprises a source component corresponding to said physiological parameter and an interference component;

processing said plurality of signals to obtain a plurality of principal components;

processing said plurality of principal components to obtain a plurality of independent components, wherein a matrix of said plurality of signals corresponds to a matrix product of a matrix of said plurality of independent components and a matrix of mixing coefficients; and

extracting a first measure of said physiological parameter corresponding to said source component from one of said plurality of independent components, wherein said first measure of a physiological parameter corresponds to a pulse rate.

18. A pulse oximeter, comprising:

a sensor configured for measuring a plurality of signals, wherein each of said signals comprises a source component corresponding to said physiological parameter and an interference component;

a computer useable medium having computer readable code embodied therein for measuring a physiological parameter, said computer readable code configured to execute functions comprising;

processing said plurality of signals to obtain a plurality of principal components;

processing said plurality of principle components to obtain a plurality of independent components, wherein a matrix of said plurality of signals corresponds to a matrix product of a matrix of said plurality of independent components and a matrix of mixing coefficients; and

extracting a first measure of said physiological parameter corresponding to said source component from one of said plurality of independent components, wherein said physiological parameter is a pulse rate.

19. A pulse oximeter, comprising:

a sensor configured for measuring a plurality of signals, wherein each of said signals comprises a source component corresponding to said physiological parameter and an interference component;

a computer useable medium having computer readable code embodied therein for measuring a physiological parameter, said computer readable code configured to execute functions comprising;

processing said plurality of signals to obtain a plurality of principal components;

processing said plurality of principle components to obtain a plurality of independent components, wherein a matrix of said plurality of signals corresponds to a matrix product of a matrix of said plurality of independent components and a matrix of mixing coefficients; and

extracting a first measure of said physiological parameter corresponding to said source component from one of said plurality of independent components, wherein said plurality of signals corresponds to sensed optical energies from a plurality of wavelengths.

20. The pulse oximeter of claim 20 wherein said physiological parameter is an oxygen saturation.

21. The pulse oximeter of claim 20 wherein said plurality of signals corresponds to the time derivative of the sensed optical energies from a plurality of wavelengths.

22. The pulse oximeter of claim 20 wherein said interference component comprises signal components caused by motion, respiratory artifact, ambient light, optical scattering and other interference between a tissue location being sense and a sensor.

23. The pulse oximeter of claim 20 wherein said processing said plurality of signals comprises decorrelating said plurality of signals by minimizing a cross-correlation of said plurality of signals, to obtain a plurality of decorrelated signals;

and

normalizing said plurality of decorrelated signals to obtain the plurality of principal components.

24. The pulse oximeter of claim 20 wherein said processing said plurality of signals comprises decorrelating said plurality of signals by singular-value decomposition of said plurality of signals, to obtain the plurality of principal components.

25. The pulse oximeter of claim 20 wherein said processing said plurality of signals comprises decorrelating said plurality of signals by multiplying said plurality of signals by the inverse square root of the covariance matrix of said plurality of signals to obtain the plurality of principal components.

26. The pulse oximeter of claim 20 wherein said processing of said plurality of principal components comprises higher-order decorrelation of said plurality of principal components.

27. The pulse oximeter of claim 20 wherein said processing said plurality of principal components comprises maximizing a function of the higher-order cumulants of a mixture of said plurality of signals, thus separating said source component from said interference component.

31. The pulse oximeter of claim 20 wherein said processing said plurality of principal components comprises successive transformations to simultaneously minimize second- and higher-order correlations among the outputs of the transformations.

32. The pulse oximeter of claim 20 wherein said processing said plurality of principal components comprises successive rotations to minimize estimated mutual information among outputs of the invention.

33. The pulse oximeter of claim 20 further comprising obtaining a ratio of mixing coefficients from said matrix of mixing coefficients, wherein said ratio corresponds to a ratio of modulation ratios of red to infrared signals.

35. The pulse oximeter of claim 20 wherein said first measure of a physiological parameter corresponds a pulse rate.

36. The pulse oximeter of claim 20 further comprising extracting said interference component from another one of said plurality of independent components.

28. The pulse oximeter of claim 28, wherein said higher-order cumulant is cumulant having order greater than two.

29. The pulse oximeter of claim 28, wherein said higher-order order cumulant is a third-order cumulant of said plurality of signals.

30. The pulse oximeter of claim 28 wherein said higher-order order cumulant is a fourth-order cumulant of said plurality of signals.

34. The pulse oximeter of claim 34 further comprising extracting a second measure of said physiological parameter from said ratio, wherein said second measure of said physiological parameter corresponds to an oxygen saturation.

37. A pulse oximeter, comprising:

a sensor configured for measuring a plurality of signals, wherein each of said signals comprises a source component corresponding to said physiological parameter and an interference component;

a computer useable medium having computer readable code embodied therein for measuring a physiological parameter, said computer readable code configured to execute functions comprising;

processing said plurality of signals to obtain a plurality of principal components;

processing said plurality of principle components to obtain a plurality of independent components, wherein a matrix of said plurality of signals corresponds to a matrix product of a matrix of said plurality of independent components and a matrix of mixing coefficients; and

extracting a first measure of said physiological parameter corresponding to said source component from one of said plurality of independent components, wherein said plurality of signals corresponds to sensed optical energies from a plurality of wavelengths from different times.