METHOD AND SYSTEM FOR DETECTION AND REJECTION OF MOTION/NOISE ARTIFACTS IN PHYSIOLOGICAL MEASUREMENTS

• US 20130191035A1
• Filed: 10/12/2011
• Published: 07/25/2013
• Est. Priority Date: 10/12/2010
• Status: Abandoned Application

• Current Assignee: Ki H. Chon, Nandakumar Selvaraj
• Original Assignee: Ki H. Chon, Nandakumar Selvaraj
• Inventors: Chon, Ki H., Selvaraj, Nandakumar
• Application Number:
  US13/822,750
• Publication Number:
  US 20130191035A1
• Time in Patent Office: Days
• US Class Current: 702/19
• CPC Class Codes: A61B 5/02416 using photoplethysmograph signals, e.g. generated by infrared radiation A61B5/14552 takes precedence
• CPC Class Codes: A61B 5/7203 for noise prevention, reduction or removal
• CPC Class Codes: A61B 5/7207 of noise induced by motion artifacts
• CPC Class Codes: A61B 5/721 using a separate sensor to detect motion or using motion information derived from signals other than the physiological signal to be measured
• CPC Class Codes: A61B 5/7214 using signal cancellation, e.g. based on input of two identical physiological sensors spaced apart, or based on two signals derived from the same sensor, for different optical wavelengths
• CPC Class Codes: A61B 5/7253 characterised by using transforms
• CPC Class Codes: G06F 2218/20 by applying autoregressive analysis
• CPC Class Codes: G16H 40/63 for local operation

Abstract

Methods and systems for quantitatively detecting the presence of artifacts in physiological measurement data and for determining usable data among those that have been designated to be corrupted with artifacts are presented.

Claims

1. A method for detection and amelioration of the effects of motion/noise artifacts in physiological measurements, the method comprising the steps of:

a) preprocessing a segment of a signal from a physiological measurement;

b) obtaining a value of at least one indicator of volatility for the preprocessed segment;

c) including the segment in analyses of physiological measurements, if comparison of the value of the at least one indicator of volatility with a predetermined threshold indicates noise/motion artifacts are not present;

d) selecting another segment of the signal from the physiological measurement and proceeding to step (a), if the value of the at least one indicator of volatility is less than a predetermined threshold and another segment is available;

e) performing a time-frequency spectrum analysis for the preprocessed segment, if comparison of the value of the at least one indicator of volatility with the predetermined threshold indicates noise/motion artifacts are present;

f) comparing a predetermined measure of the time-frequency spectrum analysis to a predetermined measure'"'"'s threshold;

g) including the segment in analyses of physiological measurements, if the predetermined measure is within limits determined by the predetermined measure'"'"'s threshold;

h) discarding the segment, if the predetermined measure is not within the limits determined by the predetermined measure'"'"'s threshold, and

i) selecting another segment of the signal from the physiological measurement and proceeding to step (a), if another segment is available;

j) repeat steps a) through i), if another segment is available.

2. The method of claim 1 wherein said at least one measure of volatility comprises kurtosis.

3. The method of claim 1 wherein said at least one measure of volatility comprises Shannon entropy.

4. The method of claim 1 wherein said at least one measure of volatility comprises kurtosis and Shannon entropy.

5. The method of claim 1 wherein said physiological measurement is a pulse oximeter waveform, referred to as a Photoplethysmogram (PPG).

6. The method of claim 5 wherein said at least one measure of volatility comprises a quadratic phase coupling between a fundamental heart rate frequency and a first harmonic of the fundamental heart rate frequency.

7. The method of claim 5 wherein said at least one measure of volatility comprises kurtosis.

8. The method of claim 7 wherein the predetermined threshold is determined using receiver operator characteristic (ROC) analysis.

9. The method of claim 5 wherein said at least one measure of volatility comprises Shannon entropy.

10. The method of claim 9 wherein the predetermined threshold is determined using receiver operator characteristic (ROC) analysis.

11. The method of claim 5 wherein said at least one measure of volatility comprises kurtosis and Shannon entropy.

12. The method of claim 11 wherein a first decision indicator is one if the kurtosis is less than a predetermined kurtosis threshold;

wherein a second decision indicator for is one if the Shannon entropy is less than another predetermined kurtosis threshold; and

wherein a joint decision indicator is 1 if the sum of the first decision indicator and the second decision indicator is equal to 2;

the joint decision indicator having a value of 1 indicates that noise/motion artifacts are not present.

13. The method of claim 5 wherein said at least one measure of volatility comprises kurtosis and Shannon entropy and a quadratic phase coupling between a fundamental heart rate frequency and a first harmonic of the fundamental heart rate frequency.

14. The method of claim 5 wherein the time-frequency spectrum analysis is performed using a variable frequency complex demodulation method;

and wherein the predetermined measure is a largest instantaneous amplitude within a frequency band centered on the heart rate.

22. The method of claim 7 wherein the predetermined threshold is determined using receiver operator characteristic (ROC) analysis.

15. A system for detection and amelioration of the effects of motion/noise artifacts in physiological measurements, the system comprising:

at least one processor; and

computer usable media having computer readable code embodied therein, the computer readable code causing said at least one processor to;

a) preprocess a segment of a signal from a physiological measurement;

b) obtain a value of at least one indicator of volatility for the preprocessed segment;

c) include the segment in analyses of physiological measurements, if comparison of the value of the at least one indicator of volatility with a predetermined threshold indicates noise/motion artifacts are not present;

d) select another segment of the signal from the physiological measurement and proceeding to step (a), if the value of the at least one indicator of volatility is less than a predetermined threshold and another segment is available;

e) perform a time-frequency spectrum analysis for the preprocessed segment, if comparison of the value of the at least one indicator of volatility with the predetermined threshold indicates noise/motion artifacts are present;

f) compare a predetermined measure of the time-frequency spectrum analysis to a predetermined measure'"'"'s threshold;

g) include the segment in analyses of physiological measurements, if the predetermined measure is within limits determined by the predetermined measure'"'"'s threshold;

h) discard the segment, if the predetermined measure is not within the limits determined by the predetermined measure'"'"'s threshold, and

i) select another segment of the signal from the physiological measurement and proceeding to step (a), if another segment is available;

j) repeat steps a) to i), if another segment is available.

16. The system of claim 15 wherein said at least one measure of volatility comprises kurtosis.

17. The system of claim 15 wherein said at least one measure of volatility comprises Shannon entropy.

18. The system of claim 15 wherein said at least one measure of volatility comprises kurtosis and Shannon entropy.

19. The system of claim 15 wherein said physiological measurement is a pulse oximeter waveform, referred to as a Photoplethysmogram (PPG).

20. The system of claim 19 wherein said at least one measure of volatility comprises a quadratic phase coupling between a fundamental heart rate frequency and a first harmonic of the fundamental heart rate frequency.

21. The system of claim 19 wherein said at least one measure of volatility comprises kurtosis.

23. The system of claim 19 wherein said at least one measure of volatility comprises Shannon entropy.

24. The system of claim 23 wherein the predetermined threshold is determined using receiver operator characteristic (ROC) analysis.

25. The system of claim 19 wherein said at least one measure of volatility comprises kurtosis and Shannon entropy.

26. The system of claim 25 wherein a first decision indicator is one if the kurtosis is less than a predetermined kurtosis threshold;

wherein a second decision indicator for is one if the Shannon entropy is less than another predetermined kurtosis threshold; and

wherein a joint decision indicator is 1 if the sum of the first decision indicator and the second decision indicator is equal to 2;

the joint decision indicator having a value of 1 indicates that noise/motion artifacts are not present.

27. The system of claim 19 wherein said at least one measure of volatility comprises kurtosis and Shannon entropy and a quadratic phase coupling between a fundamental heart rate frequency and a first harmonic of the fundamental heart rate frequency.

28. The system of claim 19 wherein the time-frequency spectrum analysis is performed using a variable frequency complex demodulation method;

and wherein the predetermined measure is a largest instantaneous amplitude within a frequency band centered on the heart rate.

29. A computer program product comprising:

a non-transitory computer usable medium having computer readable code embodied therein for detection and amelioration of the effects of motion/noise artifacts in physiological measurements, the computer readable code causing at least one processor to;

a. preprocess a segment of a signal from a physiological measurement;

b. obtain a value of at least one indicator of volatility for the preprocessed segment;

c. include the segment in analyses of physiological measurements, if comparison of the value of the at least one indicator of volatility with a predetermined threshold indicates noise/motion artifacts are not present;

d. select another segment of the signal from the physiological measurement and proceeding to step (a), if the value of the at least one indicator of volatility is less than a predetermined threshold and another segment is available;

e. perform a time-frequency spectrum analysis for the preprocessed segment, if comparison of the value of the at least one indicator of volatility with the predetermined threshold indicates noise/motion artifacts are present;

f. compare a predetermined measure of the time-frequency spectrum analysis to a predetermined measure'"'"'s threshold;

g. include the segment in analyses of physiological measurements, if the predetermined measure is within limits determined by the predetermined measure'"'"'s threshold;

h. discard the segment, if the predetermined measure is not within the limits determined by the predetermined measure'"'"'s threshold, and

i. select another segment of the signal from the physiological measurement and proceeding to step (a), if another segment is available;

j. repeat steps a) to i), if another segment is available.

30. The computer program product of claim 29 wherein said at least one measure of volatility comprises kurtosis.

31. The computer program product of claim 29 wherein said at least one measure of volatility comprises Shannon entropy.

32. The computer program product of claim 29 wherein said at least one measure of volatility comprises kurtosis and Shannon entropy.

33. The computer program product of claim 29 wherein said physiological measurement is a pulse oximeter waveform, referred to as a Photoplethysmogram (PPG).

34. The computer program product of claim 33 wherein said at least one measure of volatility comprises a quadratic phase coupling between a fundamental heart rate frequency and a first harmonic of the fundamental heart rate frequency.

35. The computer program product of claim 33 wherein the time-frequency spectrum analysis is performed using a variable frequency complex demodulation method;

and wherein the predetermined measure is a largest instantaneous amplitude within a frequency band centered on the heart rate.