Blind source separation of pulse oximetry signals

• US 20060111623A1
• Filed: 12/27/2005
• Published: 05/25/2006
• Est. Priority Date: 11/02/2001
• Status: Active Grant

• Current Assignee: Covidien LP (Medtronic Plc)
• Original Assignee: Nellcor Puritan Bennett Incorporated (Medtronic Plc)
• Inventors: Stetson, Paul F.
• Granted Patent: US 7,343,187 B2
• Time in Patent Office: Days
• US Class Current: 600/336
• 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 computer useable medium having computer readable code, the computer readable code comprising:

code configured to process a plurality of signals to obtain a plurality of principal components, wherein the plurality of signals corresponds to sensed optical energies from a plurality of wavelengths and wherein each of the plurality of signals comprises a source component corresponding to a physiological parameter and an interference component;

code configured to process the plurality of principle components to obtain a plurality of independent components, wherein a matrix of the plurality of signals corresponds to a matrix product of a matrix of the plurality of independent components and a matrix of mixing coefficients, and wherein the code configured to process the plurality of principal components comprises

code configured to separate the source component from the interference component; and

code configured to obtain a ratio of mixing coefficients from the matrix of mixing coefficients, wherein the ratio corresponds to a ratio of modulation ratios of red to infrared wavelengths of the plurality of wavelengths; and

code configured to extract a first measure of the physiological parameter corresponding to the source component from one of the plurality of independent components.

2. The medium of claim 1 wherein the physiological parameter comprises oxygen saturation of a person'"'"'s blood.

3. The medium of claim 1 wherein the plurality of signals corresponds to the time derivative of the sensed optical energies from the plurality of wavelengths.

4. The medium of claim 1 wherein the interference component comprises signal components caused by motion, respiratory artifact, ambient light, optical scattering, or other interference between a tissue location being sensed and a sensor, or any combination thereof.

5. The medium of claim 1 wherein the code configured to process the plurality of signals comprises:

code configured to decorrelate the plurality of signals by minimizing a cross-correlation of the plurality of signals to obtain a plurality of decorrelated signals; and

code configured to normalize the plurality of decorrelated signals to obtain the plurality of principal components.

6. The medium of claim 1 wherein the code configured to process the plurality of signals comprises code configured to decorrelate the plurality of signals by singular-value decomposition of the plurality of signals to obtain the plurality of principal components.

7. The medium of claim 1 wherein the code configured to process the plurality of signals comprises code configured to decorrelate the plurality of signals by multiplying the plurality of signals by the inverse square root of the covariance matrix of the plurality of signals to obtain the plurality of principal components.

8. The medium of claim 1 wherein the code configured to process the plurality of principal components comprises code configured to perform successive transformations to simultaneously minimize higher-order correlations among the outputs of the transformations.

9. The medium of claim 1 wherein the code configured to process the plurality of principal components comprises code configured to perform successive rotations to minimize estimated mutual information among outputs of the rotations.

10. The medium of claim 1 comprising code configured to extract a second measure of the physiological parameter from the ratio, wherein the second measure of the physiological parameter corresponds to an oxygen saturation.

11. The medium of claim 1 wherein the first measure of a physiological parameter corresponds to a pulse rate.

12. A pulse oximeter, comprising:

a processor configured to;

process a plurality of signals to obtain a plurality of principal components, wherein each of the plurality of signals comprises a source component corresponding to a physiological parameter and an interference component;

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

extract a first measure of the physiological parameter corresponding to the source component from one of the plurality of independent components, wherein the physiological parameter comprises a pulse rate.

13. A pulse oximeter, comprising:

a processor configured to;

process the plurality of signals to obtain a plurality of principal components, wherein each of the plurality of signals comprises a source component corresponding to a physiological parameter and an interference component;

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

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

14. The pulse oximeter of claim 13 wherein the physiological parameter comprises oxygen saturation.

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

16. The pulse oximeter of claim 13 wherein the interference component comprises signal components caused by motion, respiratory artifact, ambient light, optical scattering, or other interference between a tissue location being sense and a sensor, or any combination thereof.

17. The pulse oximeter of claim 13 wherein the processor is configured to:

decorrelate the plurality of signals by minimizing a cross-correlation of the plurality of signals to obtain a plurality of decorrelated signals; and

normalize the plurality of decorrelated signals to obtain the plurality of principal components.

18. The pulse oximeter of claim 13 wherein the processor is configured to decorrelate the plurality of signals by singular-value decomposition of the plurality of signals to obtain the plurality of principal components.

19. The pulse oximeter of claim 13 wherein the processor is configured to decorrelate the plurality of signals by multiplying the plurality of signals by the inverse square root of the covariance matrix of the plurality of signals to obtain the plurality of principal components.

20. The pulse oximeter of claim 13 wherein the processor is configured to perform a higher-order decorrelation of the plurality of principal components.

21. The pulse oximeter of claim 13 wherein the processor is configured to maximize a function of the higher-order cumulants of a mixture of the plurality of signals, thus separating the source component from the interference component.

22. The pulse oximeter of claim 21, wherein the higher-order cumulant is a cumulant having an order greater than two.

23. The pulse oximeter of claim 21, wherein the higher-order cumulant is a third-order cumulant of the plurality of signals

24. The pulse oximeter of claim 21 wherein the higher-order cumulant is a fourth-order cumulant of the plurality of signals.

25. The pulse oximeter of claim 13 wherein the processor is configured to perform successive transformations to simultaneously minimize secondhand higher-order correlations among the outputs of the transformations.

26. The pulse oximeter of claim 13 wherein the processor is configured to perform successive rotations to minimize estimated mutual information among outputs.

27. The pulse oximeter of claim 13 wherein the processor is configured to obtain a ratio of mixing coefficients from the matrix of mixing coefficients, wherein the ratio corresponds to a ratio of modulation ratios of red to infrared signals.

28. The pulse oximeter of claim 27 wherein the processor is configured to extract a second measure of the physiological parameter from the ratio, wherein the second measure of the physiological parameter corresponds to an oxygen saturation.

29. The pulse oximeter of claim 13 wherein the first measure of a physiological parameter corresponds to a pulse rate.

30. The pulse oximeter of claim 13 wherein the processor is configured to extract the interference component from another one of the plurality of independent components.

31. A pulse oximeter, comprising:

a processor configured to;

process a plurality of signals to obtain a plurality of principal components, wherein each of the plurality of signals comprises a source component corresponding to a physiological parameter and an interference component;

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

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