Improved method and apparatus for detecting optical pulses

Abstract: A method and apparatus for improving the calculation of oxygen saturation and other blood constituents by non-invasive pulse oximeters. The method and apparatus permit more accurate determina­ tion of blood flow by collecting time-measures of the absorption signal at two or more wavelengths and processing the collected time-measure to obtain composite pulsatile flow data from which artifacts have been filtered. The processing occurs in the frequency domain. The time-measure is Fourier trans­ formed into its spectral components to form the composite information.

Description: This invention relates to non-invasive pulse oximetry and specifically to an improvement on the method and apparatus for photoelectric determination of blood constituents disclosed in WO 86/05674. Non-invasive photoelectric pulse oximetry has been previously described in US-A-4,407,290, US-A-4,266,554, US-A-4,086,915, US-A-3,998,550, US-A-3,704,706, EP-A-102,816 EP-A-104,772 and EP-A-104,771. Pulse oximeters are commercially available from Nellcor Incorporated, Hayward, California, U.S.A., and are known as, for example, Pulse Oximeter Model N-100 (herein "N-100 oximeter"). Pulse oximeters typically measure and display various blood flow characteristics including but not limited to blood oxygen saturation of hemoglobin in arterial blood, volume of individual blood pulsations supplying the flesh, and the rate of blood pulsations corresponding to each heartbeat of the patient. The oximeters pass light through human or animal body tissue where blood perfuses the tissue such as a finger, an ear, the nasal septum or the scalp, and photoelectrically sense the absorption of light in the tissue. The amount of light absorbed is then used to calculate the amount of blood constituent being measured. The light passed through the tissue is selected to be of one or more wavelengths that is absorbed by the blood in an amount representative of the amount of the blood constituent present in the blood. The amount of transmitted light passed through the tissue will vary in accordance with the changing amount of blood constituent in the tissue and the related light absorption. For example, the N-100 oximeter is a microprocessor controlled device that measures oxygen saturation of hemoglobin using light from two light emitting diodes ("LED'S"), one having a discrete frequency of about 660 nanometers in the red light range and the other having a discrete frequency of about 925 nanometers in the infrared range. The N-100 oximeter microprocessor uses a four-state clock to provide a bipolar drive current for the two LED'S so that a positive current pulse drives the infrared LED and a negative current pulse drives the red LED to illuminate alternately the two LED'S so that the incident light will pass through, e.g., a fingertip, and the detected or transmitted light will be detected by a single photodetector. The clock uses a high strobing rate, e.g., one thousand five hundred cycles per second, to be easily distinguished from other light sources. The photodetector current changes in response to the red and infrared light transmitted in sequence and is converted to a voltage signal, amplified, and separated by a two-channel synchronous detector -- one channel for processing the red light waveform and the other channel for processing the infrared light waveform. The separated signals are filtered to remove the strobing frequency, electrical noise, and ambient noise and then digitized by an analog to digital converter ("ADC"). As used herein, incident light and transmitted light refers to light generated by the LED or other light source, as distinguished from ambient or environmental light. The light source intensity may be adjusted to accomodate variations among patients' skin color, flesh thickness, hair, blood, and other variants. The light transmitted is thus modulated by the absorption of light in the variants, particularly the arterial blood pulse or pulsatile component, and is referred to as the plethysmograph waveform, or the optical signal. The digital representation of the optical signal is referred to as the digital optical signal. The portion of the digital optical signal that refers to the pulsatile component is labeled the optical pulse. The detected digital optical signal is processed by the microprocessor of the N-100 oximeter to analyze and identify arterial pulses and to develop a history as to pulse periodicity, pulse shape, and determined oxygen saturation. The N-100 oximeter microprocessor decides whether or not to accept a detected pulse as corresponding to an arterial pulse by comparing the detected pulse against the pulse history. To be accepted, a detected pulse must meet certain predetermined criteria, for example, the expected size of the pulse, when the pulse is expected to occur, and the expected ratio of the red light to infrared light of the detected optical pulse in accordance with a desired degree of confidence. Identified individual optical pulses accepted for processing are used to compute the oxygen saturation from the ratio of maximum and minimum pulse levels as seen by the red wavelength compared to the maximum and minimum pulse levels as seen by the infrared wavelength. Several alternate methods of processing and interpreting optical signal data have been disclosed in the patents and references cited above. A problem with non-invasive pulse oximeters is that the plethysmograph signal and the optically derived pulse rate may be subject to irregular variants in the blood flow, including but not limited to motion artifact, that interfere with the detection of the blood flow characteristics. Motion artifact is caused by the patient's muscle movement proximate to the oximeter sensor, for example, the patient's finger, ear or other body part to which the oximeter sensor is attached, and may cause spurious pulses that are similar to pulses caused by arterial blood flow. These spurious pulses, in turn, may cause the oximeter to process the artifact waveform and provide erroneous data. This problem is particularly significant with infants, fetuses, or patients that do not remain still during monitoring. A second problem exists in circumstances where the patient is in poor condition and the pulse strength is very weak. In continuously processing the optical data, it can be difficult to separate the true pulsatile component from artifact pulses and noise because of a low signal to noise ratio. Inability to reliably detect the pulsatile component in the optical signal may result in a lack of the information needed to calculate blood constituents. It is well known that electrical heart activity occurs simultaneously with the heartbeat and can be monitored externally and characterized by the electrocardiogram ("ECG") waveform. The ECG waveform, as is known to one skilled in the art, comprises a complex waveform having several components that correspond to electrical heart activity. The QRS component relates to ventricular heart contraction. The R wave portion of the QRS component is typically the steepest wave therein, having the largest amplitude and slope, and may be used for indicating the onset of cardiovascular activity. The arterial blood pulse flows mechanically and its appearance in any part of the body typically follows the R wave of the electrical heart activity by a determinable period of time that remains essentially constant for a given patient. See, e.g., Goodlin et al., "Systolic Time Intervals in the Fetus and Neonate", Obstetrics and Gynecology , Vol. 39, No. 2, February 1972, where it is shown that the scalp pulse of fetuses lag behind the ECG "R" wave by 0.03-0.04 second, and U.S. Patent 3,734,086. In WO 86/05674, there is disclosed an invention for measuring the patient's heart activity and correlating it with the patient's detected blood flow signals to calculate more accurately the patient's oxygen saturation and pulse rate. The correlation includes auto- and cross correlation techniques to enhance the periodic information contained in each individual waveform as well as determine the time relationship of one waveform to another. Correlating the occurrence of cardiovascular activity with the detection of arterial pulses occurs by measuring an ECG signal, detecting the occurrence of the R-wave portion of the ECG signal, determining the time delay by which an optical pulse in the detected optical signal follows the R-wave, and using the determined time delay between an R-wave and the following optical pulse so as to evaluate arterial blood flow only when it is likely to present a true blood pulse for waveform analysis. The measured time delay is used to determine a time window when, following the occurrence of an R-wave, the probability of finding an optical pulse corresponding to a true arterial pulse is high. The time window provides an additional criterion to be used in accepting or rejecting a detected pulse as an optical pulse. Any spurious pulses caused by motion artifact or noise occurring outside of that time window are typically rejected and are not used to calculate the amount of blood constituent. Correlating the ECG with the detected optical pulses thereby provides for more reliable measurement of oxygen saturation. That publication refers to a modified N-100 oximeter (the "enhanced N-100 oximeter") whereby the device is provided with an additional heart activity parameter in the form of a detected R-wave from the patient's ECG waveform, in addition to the N-100 pulse oximeter functions, and the microprocessor is modified to include software and memory for controlling and processing the optical signal and heart activity information. The additional heart activity parameter is independent of the detection of peripheral arterial pulses, e.g., ECG signals, ultrasound, ballistocardiogram, and maybe, accelerometers, nuclear magnetic resonators, electrical impedance techniques, and the like, and provides an identifiable and detectable signal in response to each heartbeat for use by the signal processing of the oximeter. It is an object of this invention to provide for improved processing of the detected optical signal containing periodic information corresponding to arterial pulsatile blood flow and aperiodic information corresponding to noise, spurious signals, and motion artifact unrelated to the beating heart and arterial pulsatile blood flow, to improve further the reliability and accuracy of the determination of blood constitutent, particularly oxygen saturation of hemoglobin by a non-invasive oximeter device. It is another object of this invention to provide an improved method and apparatus for collecting successive portions of detected optical signals encompassing periodic information for more than one heartbeat and processing the collected portions to attenuate and filter therefrom aperiodic signal waveforms to provide enhanced periodic information from which the patient's blood constituent can be accurately determined. It is another object to maintain the enhanced periodic information updated by continuing to add new portions of detected optical signals as they are obtained. It is another object of this invention to evaluate the collected periodic information for a predetermined number of successive portions of the detected optical signal corresponding to a predetermined number of heartbeats in the frequency domain to obtain enhanced periodic information. It is another object of this invention to Fourier transform a time-measure of detected optical signals including periodic information for N heartbeats to determine the relative maxima at the fundamental frequency N and the minima at the zero frequency for use in determining the light modulation ratio for the amount of blood constituents. It is another object of this invention to correlate the Fourier Transform of the time-measure of detected optical signals with the Fourier Transform of a time-measure of the ECG signal, and more particularly the R-wave events of the ECG signal, to determine the maxima at the fundamental heart frequency. It is another object of this invention to correlate the periodic information in a time-measure of the detected optical signal with a time-measure of the detected heart activity, preferably in the form of the ECG signal and more preferably in the form of the R-wave of the ECG signal, to define a predetermined number of samples in a data set and use frequency domain analysis techniques to evaluate the collected predetermined number of sample data sets to determine the relative maxima at the fundamental frequency. This invention provides enhanced periodic information with improved rejection of noise, spurious pulses, motion artifact, and other undesired aperiodic waveforms and thereby improves the ability of oximeters to accurately determine amounts of blood constituents. The present invention provides methods and apparatus for collecting a time-measure of the detected optical signal waveform containing a plurality of periodic information corresponding to arterial pulses caused by the patient's heartbeat and aperiodic information unrelated to pulsatile flow, and processing the collected time-measure of information to obtain enhanced periodic information that is closely related to the most recent arterial pulsatile blood flow. The time-measure may comprise a continuous portion of detected optical signals including a plurality of periodic information from successive heartbeats, or a plurality of discrete portions of detected optical signals including a corresponding plurality of periodic information. By updating the time-measure of information to include the most recently detected periodic information, and processing the updated measure collectively, an updated enhanced periodic information is obtained (including the new and historical data) from which aperiodic information (including any new aperiodic information) is attenuated. Applicants have discovered that by collectively processing a time-measure including successive periodic information to obtain the enhanced periodic information, and using the enhanced periodic information as the basis for making oxygen saturation calculations, the accuracy and reliability of oxygen-saturation determinations can be significantly increased. Applicants also have discovered that a time-measure of detected optical signals containing a plurality of periodic information corresponding to successive heartbeats can be collectively processed and analyzed using frequency domain techniques. These frequency domain techniques utilize the synchronous nature of the heartbeat and the asynchronous characteristics of noise, spurious signals, and motion artifacts. The amount of a blood constituent, for example, oxygen saturation, can be then determined from this enhanced periodic information (also referred to as composite signal information) by determining the relative maxima and minima in the enhanced periodic information for the respective wavelengths for use in determining the modulation ratios of the known Lambert-Beers equations. In the preferred embodiment, the detected optical signals are conventionally obtained by passing red (660 namometers) and infrared (910 namometers) light through a patient's blood perfused tissue, detecting the transmitted light which is modulated by the blood flow, and providing red and infrared detected optical signals that are preferably separately processed and optionally converted from analog to digital signals, for example, as described above for the Nellcor N-100 oximeter. Portions of the corresponding red and infrared digital signals are then collectively processed in accordance with the present invention and the light modulation ratios are determined based on the resulting enhanced periodic information and used to calculate oxygen saturation. In the frequency domain, the optical signals for a given wavelength corresponding to the pulsatile arterial blood flow have spectral components including a zero frequency at the background intensity level, a fundamental frequency at the frequency of the beating heart, and additional harmonic frequencies at multiples of the fundamental frequency. Noise, spurious signals, and motion artifact that appear in the detected optical signal have frequencies that spread across the spectrum. Transient changes in the average background intensity level have frequencies that appear spread out between the zero frequency and the fundamental frequency. The present invention provides a method and apparatus for collecting a time-measure of detected optical signals including a predetermined number of optical pulses, converting the collected detected optical signals into the frequency domain, and analyzing the spectral components of the frequency spectrum thereby to determine the red and infrared relative maxima intensity at the fundamental frequency, and relative minima at the background intensity zero frequency, for use as maxima and minima in the percentage modulation ratio for calculating oxygen saturation. Applicants have discovered that if the digitized values of the time domain detected optical signals are stored in memory for a period of N heartbeats, and the stored data set is transformed into the frequency domain using Fourier Transforms, the amplitude of the fundamental heartrate is summed for the N heartbeats and appears in the frequency spectrum at a location of N cycles. In contrast, the amplitude of asynchronous signals is 1/m, where m is the number of data points in the digitized stored data set, and appear spread across the frequency domain spectrum. The average intensity of the detected optical signal background intensity appears at the spectral line corresponding to zero cycles and corresponds to the average background intensity for that wavelength. If the detected optical signal for the red and infrared signals is considered as a single complex data set, i.e., having real and imaginary components, only a single Fourier transform is required to analyze the spectral contents of the collective time-measure of the two signals. If F(s) represents the Fourier Transform of the complex data set f(t) = Red(t) + jIR(t) (for Red(t) being the red detected optical signal and IR(t) being the infrared detected optical signal), the Fourier Transform of the real component of f(t) is found by F{Re[f(t)]} = 1/2{F(s) + F*(-s)}. Similarly, the Fourier transform of the imaginary component of f(t) is found by F{Im[f(t)]} = 1/2{F(s) - F*(-s)}. F*(-s) is the complex conjugate of F(s) with the index s reversed. The relative amplitudes of the red and infrared fundamentals at the heartrate may be found by searching the frequency spectrum in the region of expected heart rates for a relative maximum and insuring that this is the fundamental by determining the existence of another relative maximum at twice this rate. This provides a technique for obtaining relative modulation data to calculate arterial oxygen saturation without the need to identify the heart rate independently, e.g., by detecting the ECG. Alternately, the amplitude data at the fundamental may be found by the use of independent heart rate determining mechanism such as ECG or phonoplethysmography or the like to determine a heart rate. However, unlike the time domain techniques, the precise time of occurrence of each heartbeat need not be determined and the optical signal and a heart rate parameter need not be correlated to obtain accurate saturation values. Rather, it is sufficient to obtain an approximate indicator of heart rate, which will facilitate identification of the fundamental frequency and improve saturation reliability. The number of spectral lines computed is preferably optimized to include the expected range of clinically applicable heartbeats (from 20-250 beats per minute), while the length of the data set is selected by the allowable equivalent delay in displaying measured arterial oxygen saturaton. A time-measure of data of, for example, 9-10 seconds represent delays of only 4-5 seconds in the display of computed saturations, and, depending upon the computational speed of the oximeter microprocessor, the time-measure can be updated in timely fashion every 1 to 2 seconds. In the preferred embodiment, the optical signal is digitized at 57 samples per second for each red and infrared signal. When 512 data points are accumulated, the data is Fourier transformed, and the red and infrared fundamental maxima are located. The percentage modulation ratio (red/infrared) is computed by dividing the energy at each maxima by the zero cycle background intensity for that wavelength, then dividing the red modulation by the infrared modulation. The resultant ratio, R, is then used in the manner set forth in the Lambert-Beers equations for calculating arterial saturation of hemoglobin. The collective data can be updated so that new data points replace the oldest data points by using a push down stack memory or equivalent so that the transform, evaluation and saturation calculation could be made after each new data set was obtained. An alternative embodiment of the frequency domain analysis technique includes sampling the real time ECG waveform and the real time detected optical signal at high rates, e.g. 1000 samples per second. By examining the ECG wave, the time of occurrence for each heartbeat and the appropriate sample rate to obtain m samples during that heartbeat could be determined. Thus, the data set for each heartbeat can be selected to contain the same number of m samples, where each sample is a fraction of the heartbeat period, and N heartbeats contains mxN samples. Taking the Fourier transform of this mxN data set and processing the spectral components of the transform in the same manner as described previously, results in a spectral analysis having several additional advantages. First, the fundamental maximum would always occur at the spectral line for N cycles in "heartbeat" space. Second, any signal present in the data set which did not remain synchronous with the heart, including noise, artifact and transient background intensity changes, would be spread over the heartbeat spectrum. Third, the enhancement in signal-to-noise would be the same for all heart rates. Fourth, because only two spectral lines are of interest, the zero spectral line corresponding to the zero frequency background intensity and the N spectral line corresponding to the number of heartbeats for the data set, the Fourier Transform need only be made at the two frequency components and not of the entire spectrum, and the computation efforts required by the microprocessor are significantly diminished. The apparatus of the present invention includes inputs for the plethysmographic detected optical signals and, optionally, ECG signals of a patient, an analog to digital converter for converting the analog plethysmographic signal to the digital optical signals and for converting the analog ECG signals into digital ECG signals (unless the plethysmographic or ECG signals are provided in digital form), and a digital signal processing section for receiving the digital signals and processing the digital detected optical signal in accordance with one of the foregoing analysis techniques of the present invention, including a microprocessor, memory devices, buffers, software for controlling the microprocessor, and display devices. In its context, the apparatus of the present invention is a part of an oximeter device which has the capability to detect the red and infrared light absorption, and receive an ECG signal from the patient. In the preferred embodiment, the apparatus of this invention is a part of the Nellcor N-200 Pulse Oximeter (herein the "N-200 oximeter"), a commercially available noninvasive pulse oximeter device manufactured and sold by Nellcor, Incorporated, Hayward, California U.S.A. The N-200 oximeter is an improved version of the enhanced N-100 oximeter described above and in International Publication No. WO86/05674. The N-200 includes circuits that perform many of the same functions as in the N-100 device, but includes some changes, for example, to expand the dynamic range of the device over the N-100 device and to include a 16 bit microprocessor manufactured by Intel Corporation, Model No. 8088. The N-100 oximeter uses an 8 bit microprocessor manufactured by Intel Corporation, Model 8085. The N-200 oximeter includes software for controlling the microprocessor to perform the operations of the conventional oximeter functions, and has some structure and processing methods that are unrelated to the present invention, and therefore are not discussed herein. The software could be modified to perform the frequency domain analysis techniques of the present invention. The invention is described in detail in connection with the drawings in which Figs. 1A and 1B are a block diagram of the apparatus of this invention and the apparatus associated with the present invention. Fig. 2A is a detailed circuit schematic of the saturation preamplifier in the patient module of Fig. 1. Fig. 2B is a detailed circuit schematic of the ECG preamplifier and input protection circuit in the patient module of Fig. 1. Figs. 3A and 3B are a detailed circuit schematic of the saturation analog front end circuit of Fig. 1. Fig. 4 is a detailed circuit schematic of the LED drive circuit of Fig. 1. Fig. 5 is a detailed circuit schematic of the ECG analog front end circuit of Fig. 1. Figs. 6A and 6B are a detailed circuit schematic of the analog to digital converter section of Fig. 1. Figs. 7A, 7B and 7C are a detailed circuit schematic of the digital signal processing section of Fig. 1. Fig. 8 is a detailed circuit schematic of the external ECG circuitry of Fig. 1. Fig. 9 is a flow chart for the frequency domain optical pulse processing of this invention. Figs. 10A, 10B, 10C, 10D and 10E are a series of waveforms corresponding to the flow chart of Fig. 10. Referring to Figs. 1A and 1B, the preferred embodiment of the present invention relates to the apparatus for processing the detected analog optical signal and the analog ECG signal and comprises portions of analog to digital conversion section ("ADC converter") 1000 and digital signal processing section ("DSP") 2000, including the software for driving microprocessor 2040, which processes the digitized optical signals and ECG signals to determine the oxygen saturation of hemoglobin in arterial blood. Associated with the invention, but not forming a part of the invention, is the apparatus for obtaining the detected analog optical signals and the analog ECG signals from the patient that is part of or is associated with the commercially available Nellcor N-200 Pulse Oximeter. Such apparatus include plethysmograph sensor 100 for detecting optical signals including periodic optical pulses, patient module 200 for interfacing plethysmograph sensor 100 and the conventional ECG electrodes with saturation analog front end circuit 300 and ECG analog front end circuit 400 respectively, saturation analog circuit 300 for processing the detected optical signals into separate red and infrared channels that can be digitized, and ECG analog front end circuit 400 for processing the ECG signal so that it can be digitized. The N-200 oximeter also includes external ECG input circuit 500 for receiving an external ECG signal and processing the signal so that it is compatible with the N-200 processing techniques (as explained below), LED drive circuit 600 for strobing the red and infrared LEDs in plethysmograph sensor 100 at the proper intensity to obtain a detected optical signal that is acceptable for processing, and various regulated power supplies (not shown) for driving or biasing the associated circuits, as well as ADC 1000 and DSP 2000, from line current or storage batteries. The associated elements are straightforward circuits providing specified functions which are within the skill of the ordinary engineer to design and build. The associated elements are briefly described here, and reference is made to the corresponding detailed schematics in the Figures and circuit element tables set forth below, to place the apparatus for using the present invention in its operating context in the preferred embodiment. In the preferred embodiment, the invention requires two input signals, the two plethysmograph or detected optical signals (e.g., red and infrared) and, optionally, a third signal, the ECG signal of the patient. If analog signals are provided, they must be within or be adjusted by, for example, offset amplifiers, to be within the voltage input range for the ADC. In circumstances where the signals have been digitized already, they must be bit compatible with the digital signal processing devices, DSP. The plethysmograph signal is obtained in a conventional manner for a non-invasive oximeter, typically by illuminating the patients tissue with red and infrared light in an alternating fashion, in the manner described above for the N-100 oximeter. Referring to Figs. 1A and 1B, sensor circuit 100 has red LED 110 and infrared LED 120 connected in parallel, anode to cathode, so that the LED drive current alternately illuminates one LED and then the other LED. Circuit 100 also includes photodetector 130, preferably a photodiode, which detects the level of light transmitted through the patient's tissue, e.g., finger 140, as a single, analog optical signal containing both the red and infrared light plethysmographic, detected optical signal waveforms. Referring to Figs. 1A, 1B, 2A, and 2B, patient module 200 includes preamplifier 210 for preamplifying the analog detected optical signal of photodetector 130, ECG preamplifer 220 for preamplifying the analog ECG signal detected from the ECG electrodes that would be attached to the patient in a conventional manner, and protection circuitry 250 interposed between instrumentation amplifier 220 and inverter 230 and the three ECG signal leads, to prevent high voltage transients from damaging the ECG preamplifier electronics. Preamplifier 210 may be an operational amplifier configured as a current to voltage converter, biased by a positive voltage to extend the dynamic range of the system, thereby converting the photocurrent of photodiode 130 into a usable voltage signal. ECG preamplifier 220 is preferably a high quality instrumentation amplifier which amplifies the differential signal present on the two ECG signal electrodes. The common-mode signal present on the two signal electrodes is inverted by inverter 230 and returned to the patient by the third ECG lead, effectively nulling the common-mode signals. A biasing network on the two ECG signal leads is provided to aid in the detection of when an ECG electrode lead becomes disconnected from patient module 200 or the patient. Patient module 200 also includes leads for passing the LED drive voltages to LEDs 110 and 120. Referring to Figs. 1A, 1B, 3A and 3B, saturation analog front end circuit 300 receives the analog optical signal from patient module 200 and filters and processes the detected signal to provide separate red and infrared analog voltage signals corresponding to the detected red and infrared optical pulses. The voltage signal is passed through low pass filter 310 to remove unwanted high frequency components above, for example, 100 khz, AC coupled through capacitor 325 to remove the DC component, passed through high pass filter 320 to remove any unwanted low frequencies below, for example, 20 hertz, and passed through programmable gain stage 330 to amplify and optimize the signal level presented to synchronous detector 340. Synchronous detector 340 removes any common mode signals present and splits the time multiplexed optical signal into two channels, one representing the red voltage signals and the other representing the infrared voltage signals. Each signal is then passed through respective filter chains having two 2-pole 20 hertz low pass filters 350 and 360, and offset amplifier 370 and 380. The filtered voltage signals now contain the signal information corresponding to the red and infrared detected optical signals. Additionally, circuits for use in preventing overdriving the amplifiers in saturation circuit 300 may be applied, for example, level-sensing circuits 312 and 314 (located after low pass filter 310) for indicating unacceptable LED drive current, and level sensing circuit 315 (located after programmable gain amplifier 330) for indicating unacceptable input amplifier gain setting. Referring to Figs. 1A, 1B, and 5, ECG analog front end circuit 400 receives the preamplified ECG signal from patient module 200 and processes it for use with the present invention. The analog ECG signal is passed through 2-pole 40 hertz low pass filter 410 for removing unwanted frequencies above 40 hertz, and programmable notch filter 420 for removing unwanted line frequency components. Optionally, circuitry may be provided to measure the line frequency and to select an appropriate clock frequency for the notch filter. The ECG signal is then passed through low pass filter 430, preferably configured to remove further unwanted components above about 40 hertz, and in particular any frequency components that may have been generated by notch filter 420. Thereafter, the ECG signal is passed through 2-pole 0.5 hertz high pass filter 440 to remove any low-frequency baseline shifts present in the original ECG signal, and then passed through offset amplifier 450 to add an offset voltage so that the voltage is within the input signal specifications of the analog to digital converter device and the complete waveform will be properly digitized. It also is desirable to pass the signal output from low pass filter 410 into a circuit that detects whether or not the ECG signal is being detected to identify a "leads-off" condition. The signal voltage is passed through absolute value circuit 480 to take the absolute value of the low pass filter output voltage and sends the value to comparator 490. Comparator 490 compares the absolute value voltage to a reference threshold or range and, when the absolute value voltage is not within the acceptable range, comparator 490 changes state which change is input to latch 495, to indicate this condition to, for example, the microprocessor. Referring to Figs. 1A, 1B, and 8, the Nellcor N-200 device also is equipped with external ECG circuit 500 for receiving the ECG signal of a stand alone ECG detector device and processing the ECG signal so that it can be used with the N-200 oximeter and the present invention. Circuit 500 receives the external analog ECG signal, passes it across capacitor 510 to remove any DC offset voltage and then passes the signal through peak detection circuit 530. A portion of the AC coupled signal also is passed through buffer amplifier 520 and input to comparator 570. The held peak voltage is used as the reference threshold voltage that is fed to the other input of comparator 570 so that subsequent QRS complexes in the ECG signal that rise above the threshold generate a trigger signal that is transferred to DSP 2000 by an electrically isolated optical serial communication link comprising serial driving opto-isolator 580, electrically isolated optical link 590, and corresponding serial driving opto-isolator 2590 in DSP 2000. Referring to Figs. 1A, 1B, 6A, and 6B, ADC 1000 provides the analog to digital conversions required by the N-200 oximeter. The aforementioned three voltage signals, the red detected optical signal, the infrared detected optical signal, and the ECG signal (preferably the ECG signal from patient module 200), are input to ADC 1000. These three signals are conventionally multiplexed and digitized by an expanded range 12-bit analog to digital conversion technique, yielding 16-bit resolution. The input signals are passed through multiplexor 1010 and buffer amplifier 1020. The converter stage includes offset amplifier 1030, programmable gain circuitry 1040 which allows a portion of the signal to be removed and the remainder to be further amplified for greater resolution, sample and hold circuit 1050, comparator 1060, and 12-bit digital to analog converter 1080. The buffered signal is passed through offset amplifier 1030 to add a DC bias to the signal wherein a portion of the signal is removed and the balance is amplified by being passed through programmable gain circuitry 1040 to improve the resolution. The amplified signal is then passed through sample and hold circuit 1050, the output of which is fed to one input of comparator 1060. The other input of comparator 1060 is the output of digital to analog ("DAC") converter 1080 so that when the inputs to comparator 1060 are the same, the analog voltage at the sample and hold circuit is given the corresponding digital word in DAC converter 1080 which is then stored in an appropriate memory device as the digitized data for the sample, and the next sample is sent to sample and hold circuit 105 to be digitized. Referring to Figs. 1A, 1B, 4, 6A, 6B, 7A, 7B, and 7C, DAC 1080 also generates the sensor LED drive voltages, under the control of microprocessor 2040, using analog multiplexor 610, which separates the incoming analog signal into one of two channels for respectively driving the red and infrared LEDs, having respective sample and hold circuits 620 and 630, and LED driver circuit 640 for converting the respective analog voltage signals into the respective positive and negative bipolar current signals for driving LEDs 110 and 120. Alternate techniques of converting the analog signals to digital signals could be used, for example, a 16-bit analog to digital converter. Referring to Figs. 1A, 1B, 7A, 7B, and 7C, DSP 2000 controls all aspects of the signal processing operation including the signal input and output and intermediate processing. The apparatus includes 16-bit microprocessor 2040 and its associated support circuitry including data bus 10, random access memory (RAM) 2020, read only memory (ROM) 2030, a conventional LED display device 2010 (not shown in detail), system timing circuit 2050 for providing the necessary clock synchronizing and notch filter frequency signals. In the preferred embodiment, microprocessor 2040 is a model 8088 microprocessor manufactured by Intel Corporation, Santa Clara, California. Alternate microprocessors may be used, such as any of model nos. 8086, 80186, and 80286, also made by Intel Corporation. Referring to Figs. 9, 10A, 10B, 10C, 10D, 10E, the flow chart for the software operation of the preferred embodiment of the present invention are shown. A corresponding software written is the Asyst computer language, which is a commercially available language, is published in EP-A-335 357. The routine begins at 4000 with the acquisition of 512 data points for each of the digitized red and infrared optical signals, which are shown graphically at Fig. 10A. At 4010, the complex data set, f(t) = Red(t) + jIR(t) , is formed. At 4020, the "D.C.", component is formed by summing all of the data points, and the "D.C." component is then removed from the complex data set by subtraction at 4030, which is graphically shown at Fig. 10B. The resulting data is then decimated in time to 64 samples at 4040, which is illustrated in Fig. 10C, and the time decimated data is then processed by the Hamming Window function at 4050, which result is illustrated in Fig. 10D. Thereafter, the Fourier Transform is taken at 4060. The spectral components of the transform are shown in Fig. 10E. The Fourier Transforms of the red and infrared components are then calculated at 4070 in accordance with the aforementioned equations, and at 4080 the maximum value at the fundamental heart rate and the minimum value at the zero frequency are determined for each of the red and infrared transforms. The saturation ratio R is calculated as: R = peak at heartrate for red Re("D.C.") peak at heartrate for infrared Im("D.C.") The minimum values for the red and infrared waveforms are taken from the respective real and imaginary components of the "D.C." component. Thereafter, the pulse data is declared ready and saturation is calculated in accordance with the known saturation formula. With each occurrence of the heartbeat, new data is acquired, the 512 data point set is updated and the routine operates to determine the saturation ratio R. In the preferred embodiment, the blood constituent measured is the oxygen saturation of the blood of a patient. The calculation of the oxygen saturation is made based on the ratio of the pulse seen by the red light compared to the pulse seen by the infrared light in accordance with the following equation: Saturation = 100% x BR2 - R(BR1) R(BO1 - BR1) + BR2 - BO2 wherein BO1 is the extinction coefficient for oxygenated hemoglobin at light wavelength 1 (Infrared) BO2 is the extinction coefficient for oxygenated hemoglobin at light wavelength 2 (red) BR1 is the extinction coefficient for reduced hemoglobin at light wavelength 1 BR2 is the extinction coefficient for reduced hemoglobin at light wavelength 2 light wavelength 1 is infrared light light wavelength 2 is red light and R is the saturation ratio as defined above. The various extinction coefficients are determinable by empirical study as is well known to those of skill in the art. For convenience of calculation, the natural log of the ratios may be calculated by use of the Taylor expansion series for the natural log.

Claims: A method for calculating the concentration of a blood constituent from the blood flow characteristics of a patient by detecting an absorption signal corresponding to the absorption of light measured at two or more wavelengths in the patient's tissue including periodic changes in amplitude caused by periodic arterial pulses in the blood flow characteristics related to the patient's heartbeat and aperiodic changes in amplitude unrelated to the patient's heartbeat, characterized by, for each of the measured wavelengths: obtaining a time-measure of the absorption signal including periodic information and aperiodic information; processing the time-measure collectively to determine a composite waveform having relative maximum and minimum amplitudes corresponding to the periodic information in the time-measure, the composite waveform corresponding to a frequency domain transform of a complex data set derived from the time-measure; and thereafter calculating the concentration of the blood constituent from the relative maximum and minimum amplitudes of the composite periodic waveforms of the detected wavelengths. The method of claim 1 characterized in that obtaining a time-measure of the absorption signal is further characterized by updating the time-measure with the occurrence of new heartbeats so that the time-measure and the determined composite waveform include periodic information related to the new heartbeats. The method of claim 2 characterized in that obtaining a time-measure of the absorption signal is further characterized by obtaining a predetermined time-measure including a predetermined number of periodic information corresponding to the predetermined number of heartbeats. The method of claim 3 characterized in that processing the collected time-measure is further characterized by Fourier transforming the collected time-measure into the frequency domain having spectral components corresponding to the frequency components of the collected time-measure, whereby the relative maximum amplitude corresponds to the amplitude at the spectral line for the predetermined number of heartbeats in the collected time-measure, and the relative minimum amplitude and the average background intensity amplitude correspond to the zero frequency spectral component. The method of claim 4 characterized in that processing the collected time-measure is further characterized by identifying the amplitude at the spectral component corresponding to the predetermined number of heartbeats by detecting a significant spectral component amplitude at the frequency corresponding to an integral multiple of the predetermined number of heartbeats. The method of claim 4 characterized in that processing the collected time-measure is further characterized by identifying the amplitude at the spectral component corresponding to the predetermined number of heartbeats by detecting the patient's heartrate and correlating the detected heartrate for the time-measure to the spectral components of the transformed time-measure. The method of claims 2 to 6, characterized in that the absorption signal includes two wavelengths, further characterized in that collecting the time-measure is further characterized by converting the absorption signals for each of the wavelengths into digital data; collecting a first predetermined number of digitized data points for each of the wavelengths of the absorption signal; forming a complex data set wherein one of the wavelengths data correspond to the real component and the other of the wavelengths data correspond to the imaginary component; determining the background absorption signal corresponding to the zero frequency component from the complex data set and subtracting the determined background absorption component from the complex data set, thereby forming a modified data set; decimating the modified data set in time into a second predetermined number of samples; processing the second predetermined number of samples using a function selected from among the group comprising Hamming windows and similar artifact reduction window functions, thereby forming a processed data set; Fourier transforming the processed data set into the frequency domain; determining the spectral components of the first and second wavelengths at the fundamental frequency for the first pretermined number of heartbeats in the time sample from the transform; and determining the relative maximum and minimum amplitudes for the first and second wavelengths using the amplitude of the zero frequency and fundamental frequency spectral components. The method of claims 3 to 7 further characterized by determining the time of the occurrence of a heartbeat from the patient's ECG signal, and sampling the time-measure of the absorption signal to obtain a second predetermined number of samples per heartbeat, based on the determined time of occurrence of the heartbeats for each of the predetermined number of heartbeats in the time-measure. An apparatus for calculating the concentration of a blood constituent from a photoelectrically detected absorption signal corresponding to the absorption of light measured at two or more wavelengths in a patient's tissue, including periodic changes in amplitude caused by periodic arterial pulses in the blood flow characteristics that are related to the patient's heartbeat, and including aperiodic changes in amplitude unrelated to the patient's heartbeat, and including a means (130) for receiving the photoelectrically detected absorption signals of each of the measured wavelengths characterized by: means (200, 300, 400) for obtaining a time-measure of the detected absorption signals including periodic information and aperiodic information; means (1000 to 1080) for processing the obtained time-measures collectively to determine a composite waveform having relative maximum and minimum amplitudes corresponding to the periodic information in the time-measure, the composite waveform corresponding to a frequency domain transform of a complex data set derived from the time-measure, and means (2000) for calculating the concentration of the blood constituent from the relative maximum and minimum amplitudes of the composite periodic wave forms of the detected wavelengths. The apparatus of claim 9 characterized in that the means for obtaining a time-measure of the detected absorption signal is further characterized by means (220, 230, 250) for updating the time-measure with the occurrence of new heartbeats so that the time-measure and the determined composite waveform includes periodic information related to the new heartbeats. The apparatus of claim 10 characterized in that the means for obtaining a time-measure of the detected absorption signal obtains a time-measure including the periodic information corresponding to a predetermined number of heartbeats. The apparatus of claim 11 characterized in that the means for processing the collected time-measure is further characterized by means (Fig. 9) for Fourier transforming the collected time-measure into a frequency domain waveform having spectral components corresponding to the frequency components of the collected time-measure, whereby the relative maximum amplitude corresponds to the amplitude at the spectral line for the predetermined number of heartbeats in the collected time-measure, and the relative minimum amplitude and the average background intensity amplitude corresponds to the amplitude at the zero frequency spectral component. The apparatus of claim 12 characterized in that the means for processing the collected time-measure detects the amplitude at the spectral component corresponding to the predetermined number of heartbeats by detecting a significant spectral component amplitude at a frequency that is about an integral multiple of the predetermined number of heartbeats. The apparatus of claim 12 further characterized by means for detecting the heartrate of the patient, characterized in that the means for processing the collected time-measure detects the amplitude at the spectral component corresponding to the predetermined number of heartbeats and the detected heartrate. The apparatus of claims 10 to 14 characterized in that the absorption signal has two wavelengths, and the means for collecting the time-measure is further characterized by: means (1060, 1080) for converting the analog absorption signals for each of the wavelengths into digital data; means for collecting (1070) a first predetermined number of digitized data points for each of the wavelengths of the absorption signal; first processor means (2000) for forming a complex data set wherein one of the wavelengths data correspond to the real component and the other of the wavelengths data correspond to the imaginary component; second processor means for determining the background absorption signal corresponding to the zero frequency component from the complex data set and subtracting the determined background absorption component from the complex data set, thereby forming a modified data set; means for decimating the modified data set in time into a second predetermined number of samples; third processing means for processing the second predetermined number of samples using Hamming Window functions, thereby forming a processed data set; means (Fig. 9) for Fourier transforming the processed data set into the frequency domain; means for determining the spectral components of the first and second wavelengths at the fundamental frequency for the first predetermined number of heartbeats in the time sample from the transform; and means for determining the relative maximum and minimum amplitudes for the first and second wavelengths using the amplitude of the zero frequency and fundamental frequency spectral components.

Inventors: BOROSS ANDRAS, BRIGGS DEBORAH A, CORENMAN JAMES E, GOODMAN DAVID E, STONE ROBERT T, BOROSS, ANDRAS, BRIGGS, DEBORAH A., CORENMAN, JAMES E., GOODMAN, DAVID E., STONE, ROBERT T.
Applicants: NELLCOR INC, NELLCOR INCORPORATED
Application year: 1989
Main classification: A61B5
IPC codes: A61B5