Physiological metric estimation rise and fall limiting

• US 10,390,762 B2
• Filed: 12/24/2012
• Issued: 08/27/2019
• Est. Priority Date: 01/16/2012
• Status: Active Grant

• Current Assignee: Yukka Magic LLC
• Original Assignee: Valencell, Inc.
• Inventors: Romesburg, Eric Douglas
• Primary Examiner(s): Roy, Baisakhi
• Assistant Examiner(s): Mehl, Patrick M
• Application Number:
  US14/370,658
• Publication Number:
  US 20150011898A1
• Time in Patent Office: 2,437 Days
• CPC Class Codes: A61B 2560/0223 of calibration, e.g. protocols for calibrating sensors
• CPC Class Codes: A61B 5/02416 using photoplethysmograph signals, e.g. generated by infrared radiation A61B5/14552 takes precedence
• CPC Class Codes: A61B 5/02438 with portable devices, e.g. worn by the patient
• CPC Class Codes: A61B 5/14551 for measuring blood gases
• CPC Class Codes: A61B 5/6801 specially adapted to be attached to or worn on the body surface
• CPC Class Codes: A61B 5/6815 Ear
• CPC Class Codes: A61B 5/7203 for noise prevention, reduction or removal
• CPC Class Codes: A61B 5/7235 Details of waveform analysis A61B5/349 takes precedence
• 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/7257 using Fourier transforms
• CPC Class Codes: A61B 5/742 using visual displays A61B5/339 takes precedence

Abstract

Methods and apparatus disclosed herein use a filtering technique to improve the accuracy of the results achieved when processing data provided by a physiological sensor. The disclosed filtering technique corrects many of the accuracy problems associated with physiological sensors, particularly PPG sensors. Broadly, the filtering technique adjusts a current filtered estimate of a physiological metric as a function of a rate limit based on a comparison between an instantaneous estimate of the physiological metric and the current filtered estimate.

Claims

1. A physiological processor configured to process a time domain photoplethysmography waveform provided by a photoplethysmography sensor, the processor comprising:

a spectral transformer configured to;

convert the photoplethysmography waveform from the time domain to a frequency domain to generate a spectrum; and

determine, responsive to the spectrum, an instantaneous estimate of a heart rate; and

a filter configured to;

compare the instantaneous estimate to a current filtered estimate of the heart rate; and

output, to an output interface, a revised filtered estimate of the heart rate computed based on the current filtered estimate and a rate limit and responsive to the comparison between the instantaneous estimate and the current filtered estimate, said revised filtered estimate comprising a sum of the current filtered estimate and an extremum value, said extremum value comprising the extremum of the rate limit and a difference between the instantaneous estimate and the current filtered estimate.

2. The physiological processor of claim 1 wherein the physiological processor is further configured to process information provided by the photoplethysmography sensor to determine at least one additional dynamically changing vital sign different from the heart rate.

3. A method of processing data provided by a physiological sensor, the method comprising:

receiving a time domain physiological waveform from the physiological sensor;

converting the received physiological waveform from the time domain to a frequency domain to generate a physiological spectrum;

determining, using the physiological spectrum, an instantaneous estimate of a physiological metric;

comparing the instantaneous estimate to a current filtered estimate of the physiological metric;

computing a revised filtered estimate of the physiological metric based on the current filtered estimate and a rate limit and responsive to the comparison between the instantaneous estimate and the current filtered estimate, said revised filtered estimate comprising a sum of the current filtered estimate and an extremum value, said extremum value comprising the extremum of the rate limit and a difference between the instantaneous estimate and the current filtered estimate; and

outputting the revised filtered estimate to an output interface;

wherein the rate limit comprises one of a first rate limit for an increasing physiological metric and a second rate limit for a decreasing physiological metric, where a magnitude of the second rate limit relative to the current filtered estimate is different from a magnitude of the first rate limit relative to the current filtered estimate; and

wherein the revised filtered estimate comprises;

the sum of the current filtered estimate and the extremum value wherein the extremum value comprises the minimum of the first rate limit and the difference between the instantaneous estimate and the current filtered estimate when the instantaneous estimate is greater than the current filtered estimate; and

the sum of the current filtered estimate and the extremum value wherein the extremum value comprises the maximum of the second rate limit and the difference between the instantaneous estimate and the current filtered estimate when the instantaneous estimate is less than the current filtered estimate.

4. The method of claim 3 wherein outputting the revised filter estimate comprises outputting the revised filter estimate to at least one of a user interface, a transceiver, a database, a processor, and a processing function.

5. The method of claim 3 wherein the physiological metric comprises a heart rate.

6. The method of claim 3 wherein receiving the physiological waveform comprises receiving the physiological waveform from a photoplethysmography (PPG) sensor.

7. A physiological processor configured to process a time domain physiological waveform provided by a physiological sensor, the processor comprising:

a spectral transformer configured to;

convert the physiological waveform from the time domain to a frequency domain to generate a physiological spectrum; and

determine, using the physiological spectrum, an instantaneous estimate of a physiological metric; and

a filter configured to;

compare the instantaneous estimate to a current filtered estimate of the physiological metric; and

output, to an output interface, a revised filtered estimate of the physiological metric based on the current filtered estimate and a rate limit and responsive to the comparison between the instantaneous estimate and the current filtered estimate, said revised filtered estimate comprising a sum of the current filtered estimate and an extremum value, said extremum value comprising the extremum of the rate limit and a difference between the instantaneous estimate and the current filtered estimate;

wherein the rate limit comprises one of a first rate limit for an increasing physiological metric and a second rate limit for a decreasing physiological metric, where a magnitude of the second rate limit relative to the current filtered estimate is different from a magnitude of the first rate limit relative to the current filtered estimate; and

wherein the revised filtered estimate comprises;

the sum of the current filtered estimate and the extremum value wherein the extremum value comprises the minimum of the first rate limit and the difference between the instantaneous estimate and the current filtered estimate when the instantaneous estimate is greater than the current filtered estimate; and

the sum of the current filtered estimate and the extremum value wherein the extremum value comprises the maximum of the second rate limit and the difference between the instantaneous estimate and the current filtered estimate when the instantaneous estimate is less than the current filtered estimate.

8. The physiological processor of claim 7 wherein the filter comprises:

an adjustment processor configured to compute an adjustment parameter as a function of the current filtered estimate, the instantaneous estimate, and the rate limit, said function derived responsive to the comparison between the instantaneous estimate and the current filtered estimate; and

a combiner configured to combine the current filtered estimate with the adjustment parameter to compute the revised filter estimate of the physiological metric.

9. The physiological processor of claim 7 further comprising an output element configured to receive the revised filter estimate from the filter, said output element comprising at least one of a user interface, transceiver, a database, a processor, and a processing function.

10. The physiological processor of claim 7 wherein the physiological metric comprises a heart rate.

11. The physiological processor of claim 7 wherein the physiological sensor comprises a photoplethysmography (PPG) sensor.

12. A method of processing data provided by a physiological sensor, the method comprising:

receiving a time domain physiological waveform from the physiological sensor;

converting the received physiological waveform from the time domain to a frequency domain to generate a physiological spectrum;

determining, using the physiological spectrum, an instantaneous estimate of a physiological metric;

comparing the instantaneous estimate to a current filtered estimate of the physiological metric;

determining a rate limit relative to the current filtered estimate responsive to one or more properties of the physiological waveform;

computing a revised filtered estimate of the physiological metric based on the current filtered estimate and the rate limit and responsive to the comparison between the instantaneous estimate and the current filtered estimate, said revised filtered estimate comprising a sum of the current filtered estimate and an extremum value, said extremum value comprising the extremum of the rate limit and a difference between the instantaneous estimate and the current filtered estimate; and

outputting the revised filtered estimate to an output interface.

13. The method of claim 12 further comprising computing an adjustment parameter as a function of the current filtered estimate, the instantaneous estimate, and the rate limit, said function derived responsive to the comparison between the instantaneous estimate and the current filtered estimate, wherein computing the revised filter estimate comprises combining the current filtered estimate with the adjustment parameter.

14. The method of claim 12 wherein determining the rate limit relative to the current filtered estimate as a function of one or more properties of the physiological waveform comprises determining the rate limit relative to the current filtered estimate as a function of—

a sum of two or more of the properties of the physiological waveform.

15. The method of claim 12 wherein determining the rate limit relative to the current filtered estimate as a function of one or more properties of the physiological waveform comprises determining the rate limit relative to the current filtered estimate as a function of—

a product of two or more of the properties of the physiological waveform.

16. The method of claim 12 further comprising computing one or more modifiers as a function of the one or more properties of the physiological waveform, wherein determining the rate limit relative to the current filtered estimate as a function of one or more properties of the physiological waveform comprises determining the rate limit relative to the current filtered estimate as a function of the one or more modifiers and an initial rate limit.

17. The method of claim 16 wherein computing the one or more modifiers comprises computing a first modifier as a function of a boundary value associated with the physiological metric and the one or more properties of the physiological waveform, and wherein determining the rate limit relative to the current filtered estimate as a function of the one or more modifiers and the initial rate limit comprises applying the first modifier to the initial rate limit.

18. The method of claim 16 wherein at least one of the one or more properties of the physiological waveform comprise spectral characteristics of the physiological spectrum, wherein computing the one or more modifiers comprises computing a first modifier as a function of the spectral characteristics of the physiological spectrum, wherein determining the rate limit relative to the current filtered estimate as a function of the one or more modifiers and the initial rate limit comprises applying the first modifier to the initial rate limit.

19. The method of claim 18 wherein computing the one or more modifiers further comprises computing a second modifier as a function of a boundary value associated with the physiological metric and the one or more properties of the physiological waveform, and wherein determining the rate limit relative to the current filtered estimate as a function of the one or more modifiers and the initial rate limit comprises applying the first and second modifiers to the initial rate limit.

20. The method of claim 18 wherein determining the instantaneous estimate comprises identifying one or more spectral peaks in the physiological spectrum, and determining the instantaneous estimate using one or more characteristics of a first spectral peak, said first spectral peak having a larger amplitude than the remaining spectral peaks, and wherein computing the first modifier comprises computing the first modifier as a function of a ratio between a second spectral peak and the first spectral peak, said second spectral peak having a smaller amplitude than the first spectral peak and a larger amplitude than the remaining spectral peaks.

21. The method of claim 12 wherein the rate limit comprises at least one of a first rate limit for an increasing physiological metric and a second rate limit for a decreasing physiological metric, where a magnitude of the first rate limit relative to the current filtered estimate is different from a magnitude of the second rate limit relative to the current filtered estimate, the method further comprising selecting at least one of the first and second rate limits responsive to the comparison between the instantaneous estimate and the current filtered estimate.

22. A method of processing data provided by a physiological sensor, the method comprising:

receiving a time domain physiological waveform from the physiological sensor;

converting the received physiological waveform from the time domain to a frequency domain to generate a physiological spectrum;

determining, using the physiological spectrum, an instantaneous estimate of a physiological metric by processing the physiological spectrum responsive to a second waveform from a second sensor to reduce noise in the instantaneous estimate;

comparing the instantaneous estimate to a current filtered estimate of the physiological metric;

computing a revised filtered estimate of the physiological metric based on the current filtered estimate and a rate limit and responsive to the comparison between the instantaneous estimate and the current filtered estimate, said revised filtered estimate comprising a sum of the current filtered estimate and an extremum value, said extremum value comprising the extremum of the rate limit and a difference between the instantaneous estimate and the current filtered estimate; and

outputting the revised filtered estimate to an output interface.

23. The method of claim 22 wherein the second sensor comprises an inertial sensor and the second waveform comprises a time domain inertial waveform, the method further comprising converting the inertial waveform from the time domain to the frequency domain to generate an inertial spectrum, and wherein determining the instantaneous estimate comprises processing the physiological spectrum as a function of the inertial spectrum to reduce the noise in the instantaneous estimate.

24. A physiological processor configured to process a time domain physiological waveform provided by a physiological sensor, the processor comprising:

a spectral transformer configured to;

convert the physiological waveform from the time domain to a frequency domain to generate a physiological spectrum; and

determine, using the physiological spectrum, an instantaneous estimate of a physiological metric; and

a filter configured to;

compare the instantaneous estimate to a current filtered estimate of the physiological metric; and

output, to an output interface, a revised filtered estimate of the physiological metric based on the current filtered estimate and a rate limit and responsive to the comparison between the instantaneous estimate and the current filtered estimate, said revised filtered estimate comprising a sum of the current filtered estimate and an extremum value, said extremum value comprising the extremum of the rate limit and a difference between the instantaneous estimate and the current filtered estimate;

wherein the filter comprises an adjustment processor configured to determine the rate limit relative to the current filtered estimate responsive to one or more properties of the physiological waveform.

25. The physiological processor of claim 24 wherein the adjustment processor is configured to determine the rate limit relative to the current filtered estimate as a function of a sum of two or more of the properties of the physiological waveform.

26. The physiological processor of claim 24 wherein the adjustment processor is configured to determine the rate limit relative to the current filtered estimate as a function of a product of two or more of the properties of the physiological waveform.

27. The physiological processor of claim 24 wherein the adjustment processor comprises a modifier processor configured to calculate one or more modifiers as a function of the one or more properties of the physiological waveform, and wherein the adjustment processor determines the rate limit relative to the current filtered estimate by determining the rate limit relative to the current filtered estimate as a function of the one or more modifiers and an initial rate limit.

28. The physiological processor of claim 27 wherein the modifier processor comprises:

a calculator configured to compute a first modifier as a function of a boundary value associated with the physiological metric and the one or more properties of the physiological waveform; and

a modifier applicator configured to apply the first modifier to the initial rate limit to determine the rate limit relative to the current filtered estimate.

29. The physiological processor of claim 27 wherein at least one of the one or more properties of the physiological waveform comprise spectral characteristics of the physiological spectrum, and wherein the modifier processor comprises:

a calculator configured to compute a first modifier as a function of the spectral characteristics of the physiological spectrum; and

a modifier applicator configured to apply the first modifier to the initial rate limit to determine the rate limit relative to the current filtered estimate.

30. The physiological processor of claim 29 wherein the calculator is further configured to compute a second modifier as a function of a boundary value associated with the physiological metric and the one or more properties of the physiological waveform, and wherein the modifier applicator is configured to apply the first and second modifiers to the initial rate limit to determine the rate limit relative to the current filtered estimate.

31. The physiological processor of claim 29 wherein the spectral transformer is further configured to identify one or more spectral peaks of the physiological spectrum, and to determine the instantaneous estimate using one or more characteristics of a first spectral peak of the one or more spectral peaks, said first spectral peak having a larger amplitude than the remaining spectral peaks of the one or more spectral peaks, and wherein the calculator is configured to compute the first modifier as a function of a ratio between a second spectral peak of the one or more spectral peaks and the first spectral peak, said second spectral peak having a smaller amplitude than the first spectral peak and a larger amplitude than the remaining spectral peaks.

32. The physiological processor of claim 24 wherein the rate limit comprises at least one of a first rate limit for an increasing physiological metric and a second rate limit for a decreasing physiological metric, where a magnitude of the first rate limit relative to the current filtered estimate is different from a magnitude of the second rate limit relative to the current filtered estimate, and wherein the filter comprises a rate processor configured to select at least one of the first rate limit and the second rate limit responsive to the comparison between the instantaneous estimate and the current filtered estimate.

33. A physiological processor configured to process a time domain physiological waveform provided by a physiological sensor, the processor comprising:

a spectral transformer configured to;

convert the physiological waveform from the time domain to a frequency domain to determine a physiological spectrum; and

determine, based on the physiological spectrum, an instantaneous estimate of a physiological metric;

a filter configured to;

compare the instantaneous estimate to a current filtered estimate of the physiological metric; and

output, to an output interface, a revised filtered estimate of the physiological metric based on the current filtered estimate and a rate limit and responsive to the comparison between the instantaneous estimate and the current filtered estimate, said revised filtered estimate comprising a sum of the current filtered estimate and an extremum value, said extremum value comprising the extremum of the rate limit and a difference between the instantaneous estimate and the current filtered estimate; and

a second sensor configured to receive a time domain second waveform;

wherein the spectral transformer is further configured to convert the second waveform from the time domain to the frequency domain to determine a second spectrum, and to determine the instantaneous estimate responsive to the second spectrum to reduce the noise in the instantaneous estimate.

34. The physiological processor of claim 33 wherein the second sensor comprises an inertial sensor and wherein the time domain second waveform comprises a time domain inertial waveform and the second spectrum comprises an inertial spectrum, the processor further comprising a noise filter configured to filter the physiological spectrum responsive to the inertial spectrum to generate a noise-reduced waveform, wherein the spectral transformer is configured to determine the instantaneous estimate using the noise-reduced waveform to reduce the noise in the instantaneous estimate.