Method and apparatus for estimating physiological parameters using model-based adaptive filtering

• US 6,083,172 A
• Filed: 08/20/1998
• Issued: 07/04/2000
• Est. Priority Date: 08/07/1995
• Status: Expired due to Term

• Current Assignee: Nellcor Puritan Bennett Incorporated (Medtronic Plc)
• Original Assignee: Nellcor Puritan Bennett Incorporated (Medtronic Plc)
• Inventors: Yorkey, Thomas J., Baker, Clark R. Jr.
• Primary Examiner(s): Nasser, Robert L.
• Application Number:
  US09/137,479
• Time in Patent Office: 684 Days
• Field of Search: 600/300, 600/485, 600/500-504, 600/561, 600/490, 600/493-6
• US Class Current: 600/500
• CPC Class Codes: A61B 5/0075 by spectroscopy, i.e. measuring spectra, e.g. Raman spectroscopy, infrared absorption spectroscopy A61B5/0071 takes precedence
• CPC Class Codes: A61B 5/024 Detecting, measuring or recording pulse rate or heart rate A61B5/0205, A61B5/021 take precedence
• 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/14552 Details of sensors specially adapted therefor
• 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/7239 using differentiation including higher order derivatives
• CPC Class Codes: A61B 5/725 using specific filters therefor, e.g. Kalman or adaptive filters A61B5/316 take precedence
• CPC Class Codes: H03H 21/003 KALMAN filters

Abstract

A method and apparatus for reducing the effects of noise on a system for measuring physiological parameters, such as, for example, a pulse oximeter. The method and apparatus of the invention take into account the physical limitations on various physiological parameters being monitored when weighting and averaging a series of measurements. Varying weights are assigned different measurements, measurements are rejected, and the averaging period is adjusted according to the reliability of the measurements. Similarly, calculated values derived from analyzing the measurements are also assigned varying weights and averaged over adjustable periods. More specifically, a general class of filters such as, for example, Kalman filters, is employed in processing the measurements and calculated values. The filters use mathematical models which describe how the physiological parameters change in time, and how these parameters relate to measurement in a noisy environment. The filters adaptively modify a set of averaging weights to optimally estimate the physiological parameters.

Claims

1. A method for generating a pulse rate of a patient using data corresponding to at least one wavelength of electromagnetic energy transmitted through tissue of the patient, the method comprising:

defining a comb filter to isolate signal energy in the data corresponding to a fundamental frequency and related higher frequency components;

determining a particular frequency that optimized energy at an output of the comb filter; and

generating a filtered pulse rate corresponding to the particular frequency.

2. The method of claim 1, wherein the determining minimizes signal energy at the output of the comb filter.

3. The method of claim 1, wherein the determining maximizes noise energy at the output of the comb filter.

4. The method of claim 1, wherein the determining comprises:

calculating squared noise for the data;

calculating a second derivative of the squared noise with respect to the fundamental frequency; and

performing a Newton-Raphson search to determine the particular frequency.

5. The method of claim 1, further comprising:

comparing selected samples of the pulse rate to an expected measurement characteristic;

assigning a variable weight to each selected sample based on the comparing, thereby generating a plurality of weighted samples of the pulse rate;

averaging a first number of the weighted samples to obtain a filtered sample of the pulse rate; and

generating the filtered pulse rate from a plurality of filtered samples.

6. The method of claim 5, wherein the first number varies in response to the variable weights assigned to the selected samples.

7. The method of claim 1, wherein, preceding the defining, the method further comprises:

computing a logarithm of a signal representative of the at least one wavelength of electromagnetic energy to generate a first signal;

bandpass filtering the first signal to generate a second signal;

normalizing the second signal to generate a third signal; and

computing a derivative of the third signal to generate the data.

8. The method of claim 1, wherein the pulse rate is associated with a confidence level, the confidence level being based on at least one quality metric selected from the group consisting of pulse signal shape, signal-to-noise ratio, and correlation between at least two wavelengths of electromagnetic energy.

9. The method of claim 1, further comprising:

determining an alternate pulse rate using an alternate rate finder; and

selectively setting the pulse rate to the alternate pulse rate.

10. The method of claim 9, wherein the determining the alternate pulse rate comprises:

comparing the data to a predetermined waveform template;

identifying a sequence of waveform characteristics indicative of a waveform period;

averaging a number of successive waveform periods to determine an average waveform period; and

determining the alternate pulse rate from the average waveform period.

11. The method of claim 9, wherein the determining the alternate pulse rate comprises:

evaluating a power spectrum corresponding to the data to determine which of a plurality of peaks in the power spectrum corresponds to the fundamental frequency;

generating the alternate pulse rate corresponding to the fundamental frequency; and

calculating a confidence level for the alternate pulse rate.

12. The method of claim 9, further comprising:

generating an initial pulse rate, the initial pulse rate being set to the alternate pulse rate in the selectively setting.

13. The method of claim 9, wherein the determining the alternate pulse rate occurs periodically, and wherein the selectively setting occurs only when at least one condition is met.

14. The method of claim 13, wherein the at lease one condition is selected from the group comprising:

a) a confidence level corresponding to the filtered pulse rate being less than a first value;

b) the confidence level corresponding to the alternate pulse rate being greater than or equal to the first value;

c) a signal-to-noise ratio corresponding to the alternate pulse rate being greater than or equal to a second value;

d) the alternate pulse rate being a local maxima on an autocorrelation curve corresponding to the data; and

e) the alternate pulse rate being within a third value of a previously determined alternate pulse rate, the previously determined alternate pulse rate having been determined within a first interval.

15. The method of claim 14, wherein conditions a) through e) are met before the pulse rate is set to the alternate pulse rate.

16. An apparatus for generating a pulse rate of a patient using data corresponding to at least one wavelength of electromagnetic energy transmitted through tissue of the patient comprising:

means for defining a comb filter to isolate signal energy in the data corresponding to a fundamental frequency and related higher frequency components;

means for determining a particular frequency that optimizes energy at an output of the comb filter; and

means for generating a filtered pulse rate corresponding to the particular frequency.

17. The apparatus of claim 16, wherein the means for determining comprises:

means for calculating squared noise for the data;

means for calculating a second derivative of the squared noise with respect to the fundamental frequency; and

means for performing a Newton-Raphson search to determine the particular frequency.

18. The apparatus of claim 16, further comprising:

means for comparing selected samples of the pulse rate to an expected measurement characteristic;

means for assigning a variable weight to each selected sample based on the comparing, thereby generating a plurality of weighted samples of the pulse rate;

means for averaging a first number of the weighted samples to obtain a filtered sample of the pulse rate; and

means for generating a filtered pulse rate from a plurality of filtered samples.

19. The apparatus of claim 16, further comprising:

means for computing a logarithm of a signal representative of the at least one wavelength of electromagnetic energy to generate a first signal;

means for bandpass filtering the first signal to generate a second signal;

means for normalizing the second signal to generate a third signal; and

means for computing a derivative of the third signal to generate the data.

20. The apparatus of claim 16, further comprising:

means for determining an alternate pulse rate using an alternate rate finder; and

means for selectively setting the pulse rate to the alternate pulse rate.

21. The apparatus of claim 20, wherein the means for determining the alternate pulse rate comprises:

means for comparing the data to a predetermined waveform template;

means for identifying a sequence of waveform characteristics indicative of a waveform period;

means for averaging a number of successive waveform periods to determine an average waveform period; and

means for determining the alternate pulse rate from the average waveform period.

22. The apparatus of claim 20, wherein the means for determining the alternate pulse rate comprises:

means for evaluating a power spectrum corresponding to the data to determine which of a plurality of peaks in the power spectrum corresponds to the fundamental frequency;

means for generating the alternate pulse rate corresponding to the fundamental frequency; and

means for calculating a confidence level for the alternate pulse rate.

23. A circuit for generating a pulse rate of a patient using data corresponding to at least one wavelength of electromagnetic energy transmitted through tissue of the patient comprising:

a comb filter operative to isolate signal energy in the data corresponding to a fundamental frequency and related higher frequency components;

a processor coupled to the comb filter and operative to determine a particular frequency that optimizes energy at an output of the comb filter; and

an averaging filter coupled to the processor and operative to generate a filtered pulse rate corresponding to the particular frequency.

24. The circuit of claim 23, further comprising:

an alternate rate finder operative to determine an alternate pulse rate of the patient; and

a selector coupled to the alternative rate finder and operative to selectively set the pulse rate to the alternate pulse rate.

25. The circuit of claim 24, wherein the alternate rate finder is further operative to

compare the data to a predetermined waveform template,

identify a sequence of waveform characteristics indicative of a waveform period,

average a number of successive waveform periods to determine an average waveform period, and

determine the alternate pulse rate from the average waveform period.

26. The circuit of claim 24, wherein the alternate rate finder is further operative to

evaluate a power spectrum corresponding to the data to determine which of a plurality of peaks in the power spectrum corresponds to the fundamental frequency,

generate the alternate pulse rate corresponding to the fundamental frequency, and

calculate a confidence level for the alternate pulse rate.

27. The circuit of claim 24, wherein the alternate rate finder is further operative to generate an initial pulse rate.

28. The circuit of claim 24, wherein the alternate rate finder is operative to periodically determine the alternate pulse rate, and wherein the selector is operative to set the pulse rate to the alternate pulse rate only when at least one condition is met.

29. The circuit of claim 28, wherein the at lease one condition is selected from the group comprising:

a) a confidence level corresponding to the filtered pulse rate being less than a first value,

b) the confidence level corresponding to the alternate pulse rate being greater than or equal to the first value,

c) a signal-to-noise ratio corresponding to the alternate pulse rate being greater than or equal to a second value,

d) the alternate pulse rate being a local maxima on an autocorrelation curve corresponding to the data, and

e) the alternate pulse rate being within a third value of a previously determined alternate pulse rate, the previously determined alternate pulse rate having been determined within a first interval.

30. The circuit of claim 29, wherein the selector is operative to set the pulse rate to the alternate pulse rate when conditions a) through e) are met.

31. The circuit of claim 23, wherein the comb filter is operative to minimize signal energy at the output of the comb filter.

32. The circuit of claim 23, wherein the comb filter is operative to maximize noise energy at the output of the comb filter.

33. The circuit of claim 23, wherein the processor is operative to

calculate squared noise for the data,

calculate a second derivative of the squared noise with respect to the fundamental frequency, and

perform a Newton-Raphson search to determine the particular frequency.

34. The circuit of claim 23, wherein the processor is operative to

compare selected samples of the pulse rate to an expected measurement characteristic,

assign a variable weight to each selected sample based on the compare, thereby generating a plurality of weighted samples of the pulse rate,

average a first number of the weighted samples to obtain a filtered sample of the pulse rate, and

generate a filtered pulse rate from a plurality of filtered samples.

35. The circuit of claim 34, wherein the first number varies in response to the variable weights assigned to the selected samples.

36. The circuit of claim 23, wherein the processor is operative to

compute a logarithm of a signal representative of the at least one wavelength of electromagnetic energy to generate a first signal,

bandpass filter the first signal to generate a second signal,

normalize the second signal to generate a third signal, and

compute a derivative of the third signal to generate the data.

37. The circuit of claim 23, wherein the pulse rate is associated with a confidence level, the confidence level being based on at least one quality metric selected from the group consisting of pulse signal shape, signal-to-noise ratio, and correlation between at least two wavelengths of electromagnetic energy.