Wave form filter pulse detector and method for modulated signal

Information

  • Patent Grant
  • 4867571
  • Patent Number
    4,867,571
  • Date Filed
    Tuesday, February 28, 1989
    35 years ago
  • Date Issued
    Tuesday, September 19, 1989
    35 years ago
Abstract
A device and method for determining pulses in a modulated signal is disclosed. The modulated signal is processed and converted into a quantized analog type of digital data stream. The data stream is evaluated over time by considering preceeding and subsequent values in the data stream to generate a filtered wave form. By using extreme values in the filtered wave form pulse detection is accurately determined regardless of whether the modulated signal has complicating features, such as dicrotic notch, or high noise levels, or both.
Description

BACKGROUND AND SUMMARY OF THE INVENTION
The present invention relates to a device and method for detecting the pulse in a modulated signal. Specifically, the invention relates to the accurate detection of a pulse even if the modulated signal has complicating features such as dicrotic notch, or high noise levels, or both.
A modulated signal is generated, for example, by transmitting light through living tissue and monitoring changes in intensity of the transmitted light. Circuit means are provided for signal processing and a microprocessor for mathematical evaluation of the changes in the transmitted light. The processing includes signal separation, noise reduction, amplification and an analog to digital conversion. The processed signal is mathematically evaluated over time by considering preceeding and subsequent values in the data stream to generate a filtered wave form. By using an extreme value in the filtered wave form pulse detection can be accurately determined.





BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a data stream showing the A.C. component, the D.C. component, and the D.C. offset.
FIG. 2 is the data stream of FIG. 1 following removal of the D.C. offset and following amplification.
FIG. 3 is a noisy data stream with a complicating feature known as baseline drift.
FIG. 4 is a data stream following removal of the D.C. offset and following amplification with a complicating feature known as a dicrotic notch.
FIG. 5a is a data stream similar to FIG. 4, but with substantial noise.
FIG. 5b is the filtered data stream of FIG. 5a.
FIG. 6 is a block diagram of the hardware comprising an oximeter and pulse detector.
FIG. 7 is a process flow block diagram of the microprocessor unit of an oximeter and pulse detector.





DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT
The invention can be best understood by first examining typical analog modulated signals, or values of the data points in the data stream following the analog to digital conversion. In each instance, FIGS. 1-5b, the y-axis represents the electrical signal, e.g., voltage, and the x-axisrepresents time. FIG. 1 depicts a relatively noise free data stream withouta baseline drift. The data stream comprises a DC component and an AC component. The DC component further comprises a DC offset and DC remainder. The AC component to the DC component, is small. To simplify theevaluation of the AC signal, necessary for the pulse determination, the DC offset is removed. The remaining signal is thereafter amplified. The data stream of FIG. 1 after removal of the D.C. offset and amplification is shown in FIG. 2.
The signal of FIGS. 1 and 2 is a relatively clean data stream. However, many data streams are substantially more complicated. FIG. 3 depicts a complex data stream, that is both noisy and has a substantial baseline drift. Another data stream, shown in FIG. 4, following removal of the DC offset and following amplification has a complicating feature known as a dicrotic notch. The data stream of FIG. 5a is similar to that shown in FIG. 4. However, the FIG. 5a data stream has the additional complicating factor of a high noise level. FIG. 5b shows the data stream of FIG. 5a after filtering.
The operating principles of the pulse detector are described first. Associated with every given data point in the data stream and equidistant from that data point are multiple pairs of data points. That is, each datapoint has multiple pairs of associated data points. For each pair, the first associated data point occurs some time prior to the given data pointand the second associated data point occurs an equal amount of time after the given data point. A difference is found by subtracting the detected values for the associated data points one from another. The difference foreach pair of associated data points is then summed to form the output of a wave form filter for the given point. The total time spanned by these associated points is called the wave form filter length. This can be expressed in the following fashion: ##EQU1##Where F.sub.i is the wave form filter output for a given time t.sub.i having a data stream value P.sub.i, P.sub.i-j is the data stream value at time t.sub.i -t.sub.j ; P.sub.i+j is the data stream value at time t.sub.i +t.sub.j ; and the wave form filter length is 2*L+1. Using this approach to finding F.sub.i requires L subtractions and L-1 additions.
A simpler calculation of F.sub.i+l is possible if F.sub.i has already been computed. That is:
F.sub.i+l =F.sub.i +P.sub.i+L+l +P.sub.i-L -P.sub.i -P.sub.i+l
This calculation requires only two additions and two subtractions regardless of the length of the wave form filter. Present microprocessors are able to make this calculation in real time if the discrete points in the data stream occur at, for example, 15 millisecond intervals. For this calculation, memory of 2L+1 values is required.
The accuracy of the output of the wave form filter for pulse detection is best when the wave form filter length and the pulse length of the signal are the same. When there is a large mismatch in these two quantities, the accuracy of the filter is diminished. Two methods have been found to overcome this problem. The first is to use two or more filters and examineeach of them separately to determine which most closely matches the pulse length. The other is to combine two filters such that their combined output will work on any signal of interest. The second method requires four additions and four subtractions for each point.
The modulated signal from which the pulse is determined may be one of the two signals used in determining oxygen saturation by means of an oximeter.The functioning of the pulse detector and, to some extent, an associated oximeter is now described. In FIG. 6, there is depicted a block diagram ofthe hardware used in an oximeter. A photoelectric transducer or photodiode 10 receives the light transmitted through a measuring point in the human body such as an ear lobe or finger. For the determination of oxygen saturation two light components are transmitted through the measuring point. Light component A is transmitted from LED A and light component B is transmitted from LED B. If only the pulse is determined, only one light, either A or B, need be transmitted, processed and evaluated.
The data streams detected by photodiode 10 are amplified by pre-amplifer 12and passed through the synchronized demodulator 14 to separate the data streams for each of the two light components. For each of the two data streams, the data stream is further separated (16 and 18) into a DC offsetand the DC remainder plus the AC component. The values of the DC offset aresampled and held in the microprocessor 40 for further processing. Alternatively, the DC offset can be present at a fixed value. Once the DC offset is removed, data streams A and B are passed through operational amplifiers 20 and 22. The signal streams are amplified by fixed gains relative to the signal strengths of channel A and channel B. If the A channel processes the red signal, the fixed gain may be approximately a multiple of 200-250 of the preexisting data stream while, the B channel, if processing an infrared signal, the fixed gain may be approximately a multiple of 40-60 of the preexisting data stream.
The data streams of both channels A and B are passed through the filters 24and 26 to reduce gross extraneous noise. The signal streams are then passedthrough variable attenuators 30 and 32, the control of which is performed by an evaluation of the signal strength made by the microprocessor 40. Thesignal streams are then inputted to multiplexor 34 where they are sampled and held until the analog to digital convertor 36 has converted each incoming analog signal into an outgoing digital also referred to as a quantized analogue signal signal.
The data streams are thereafter processed as shown in FIG. 7. As each data point in the data stream is inputted into the microprocessor 40, they are stored in buffers 52 and 54. As inputs are received the values are sequentially stored in the buffers replacing previous values which are shifted through and eventually out of the buffers. When the process is first started or after data is lost, the microprocessor 40 holds until thebuffer is full before commencing a calculation.
Once the buffers fill, the wave form filter or filter outputs, if more thanone wave form filter is used, are calculated (90, 91). The initial wave form filter output is determined by subtracting each successive trailing data point in the data streamm from each successive leading data point in the data stream and then summing these values. Each successive wave form filter output is determined by summing (i) the most recent wave form filter output, (ii) the trailing data point in the data stream, and (iii) the leading data point in the data stream for the most recent wave form filter output calculation and subtracting the sum of (i) the data point inthe data stream halfway through the buffer and (ii) the data point in the data stream one data point beyond the halfway mark in the buffer. The foregoing mathematical calculation is performed for each wave form filter used.
Once the wave form filter outputs are determined, a detector 92 is used to determine an extreme value, such as a maximum or minimum. The rate of the extreme values are compared (93), by using quality criteria, with an expected range of values. If the extreme value is within the range of expected values, the pulse rate is determined (94). The pulse rate is alsocompared (95), by using quality criteria, with an expected range of pulse rates. If the pulse rate is within the range of expected values, the pulserate is outputted to the pulse display driver (96).
While the above embodiments have been disclosed as the best mode presently contemplated by the inventor, it should be realized that these examples should not be interpreted as limiting, because an artisan skilled in this field, once given the present teachings, can vary from these specific embodiments. Accordingly, the scope of the present invention should be determined solely from the following claims.
Claims
  • 1. An apparatus for the determination of a pulse in an amplitude modulated quantized analog signal having a signal amplitude complicating feature, comprising:
  • means for transmitting light through a specimen,
  • means for monitoring the transmitted light to produce a data stream comprising values representative of the attenuation of light by the specimen over time,
  • means for storing the series of values whereby each of the values is a discrete moment in time and the stored values are continually updated over time,
  • wave form filter for evaluating the stored values over time whereby a series of filtered values are determined the filtered values including a plurality of extreme values representing the pulse, and
  • means for determining a frequency of the extreme values in the series of filtered values.
  • 2. A method for determining pulses in an amplitude modulated quantized analog signal having a signal complicating feature, comprising:
  • transmitting light through a specimen,
  • monitoring the transmitted light to produce a data stream comprising a series of values representative of the attenuation of light by the specimen over time,
  • storing the series of values whereby each of the stored values is a discrete moment in time,
  • filtering the stored values whereby a series of filtered values are determined the filtered values including a plurality of extreme values representing the pulse, and
  • evaluating the series of filtered values to determine a frequency of the extreme values.
  • 3. An apparatus for the determination of a pulse in an amplitude modulated quantized analog signal having a signal amplitude complicating feature, comprising:
  • means for monitoring a data stream comprising a series of time spaced values,
  • means for storing at least a portion of the series of time spaced values whereby the stored values are continually updated over time,
  • wave form filter for evaluating the stored values over time whereby a series of wave form filter outputs are determined, the filter outputs including a plurality of extreme values representing the pulse, by solving the expression: ##EQU2## where F.sub.i is a wave form filter output for a given time t.sub.i having a data stream value P.sub.i, P.sub.i,j is the data stream value at time t.sub.i -t.sub.j, and p.sub.i+j is the data stream value at time t.sub.i +t.sub.j, and data stream value at time t.sub.i +t.sub.j, and
  • means for determining a frequency of the extreme values in the series of wave form filter outputs.
  • 4. A wave form filter for the determination of pulses in an amplitude modulated quantized analog signal having a signal amplitude complicating feature comprising a means for generating the digital signal, the quantized analog signal comprising a series of time spaced values, a means for storing and continually updating at least three of the time spaced values, a filtering means to determine a filtered wave form by using at least two of the time spaced values the filtered wave form including a plurality of extreme values representing the pulses, and a means for determining the extreme values in the filtered wave form.
  • 5. A method for determining pulses in an amplitude modulated quantized analog signal having a signal amplitude complicating feature comprising generating a modulated signal comprised of a series of time space data values, storing some of the data values which are continually updated over time, filtering the stored data values over time to determine a series of filtered values the filtered values including a plurality of extreme value in the signal representing the pulses, and evaluating the filtered values to determine the extreme values.
  • 6. The apparatus of claim 1 wherein the amplitude complicating feature is a high noise level.
  • 7. The apparatus of claim 1 wherein the amplitude complicating feature is a dicrotic notch.
  • 8. The apparatus of claim 1 wherein the amplitude complicating feature is a baseline drift.
  • 9. The method of claim 2 wherein the amplitude complicating feature is a high noise level.
  • 10. The method of claim 2 wherein the amplitude complicating feature is a dicrotic notch.
  • 11. The method of claim 2 wherein the amplitude complicating feature is a baseline drift.
  • 12. The wave form filter of claim 4 wherein the amplitude complicating feature is a high noise level.
  • 13. The wave form filter of claim 4 wherein the amplitude complicating feature is a dicrotic notch.
  • 14. The wave form filter of claim 4 wherein the amplitude complicating feature is a baseline drift.
  • 15. The method of claim 5 wherein the amplitude complicating feature is a high noise level.
  • 16. The method of claim 5 wherein the amplitude complicating feature is a dicrotic notch.
  • 17. The method of claim 5 wherein the amplitude complicating feature is a baseline drift.
Parent Case Info

This is a continuation of co-pending application Ser. No. 911,956 filed on Sept. 26, 1986, now abandoned.

US Referenced Citations (18)
Number Name Date Kind
2706927 Wood Apr 1955
3522139 Coor et al. Aug 1970
3638640 Shaw Feb 1972
3647299 Lavallee Mar 1972
3799672 Vurek Mar 1974
3842355 Fleagle Oct 1974
3938052 Glasson et al. Feb 1976
3992672 Fasching Nov 1976
3998550 Konishi et al. Dec 1976
4063551 Sweeny Dec 1977
4086915 Kofsky et al. May 1978
4093988 Scott Jun 1978
4167331 Nielson Jul 1980
4214213 Ferrie Jul 1980
4260951 Lewyn Apr 1981
4266554 Hamaguri May 1981
4402290 Wilber Oct 1983
4523279 Sperinde et al. Jun 1985
Foreign Referenced Citations (2)
Number Date Country
0102816 Mar 1984 EPX
128387 Oct 1975 JPX
Non-Patent Literature Citations (5)
Entry
Grover, "Improved Extracorporeal Reflectance Oximeter", Conference Proceedings of the 26th Annual Conference on Engineering in Medicine and Biology, Sep. 30-Oct. 4, 1973, p. 275.
Kramer et al., "Influence of Oxygen Saturation, Erythocyte Concentration and Optical Depth Upon the Red and Near-Infrared Light Transmittance of Whole Blood", American Journal of Physiology, 165:229-246 (1951).
Schibili et al., "An Electronic Circuit for Red/Infrared Oximeters", IEEE Transactions on Bio-Medical Engineering, vol. BME-25, No. 1, Jan., 1978.
Tait et al., "A Theoretical Analysis of Some Errors in Oximetry", IEE Transactions on Bio-Medical Engineering, vol. BME-13, No. 4, Oct., 1966.
Wood and Garci, "Photoelectric Determination of Arterial Oxygen Saturation in Man", J. Lab. Clin. Med. 34, pp. 387-401 (1949).
Continuations (1)
Number Date Country
Parent 911956 Sep 1986