US20140371554A1 - Method for improved oxygen saturation estimation in the presence of noise - Google Patents

Method for improved oxygen saturation estimation in the presence of noise Download PDF

Info

Publication number
US20140371554A1
US20140371554A1 US14/472,038 US201414472038A US2014371554A1 US 20140371554 A1 US20140371554 A1 US 20140371554A1 US 201414472038 A US201414472038 A US 201414472038A US 2014371554 A1 US2014371554 A1 US 2014371554A1
Authority
US
United States
Prior art keywords
oxygen saturation
estimate
processor
neural network
modified
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Abandoned
Application number
US14/472,038
Inventor
Clark R. Bake, JR.
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Covidien LP
Original Assignee
Covidien LP
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Covidien LP filed Critical Covidien LP
Priority to US14/472,038 priority Critical patent/US20140371554A1/en
Assigned to COVIDIEN LP reassignment COVIDIEN LP ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: BAKER, CLARK R., JR
Publication of US20140371554A1 publication Critical patent/US20140371554A1/en
Abandoned legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A61B5/14551Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7282Event detection, e.g. detecting unique waveforms indicative of a medical condition
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/7405Details of notification to user or communication with user or patient ; user input means using sound
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/742Details of notification to user or communication with user or patient ; user input means using visual displays
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/7455Details of notification to user or communication with user or patient ; user input means characterised by tactile indication, e.g. vibration or electrical stimulation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems

Definitions

  • the present disclosure relates, in some embodiments, to devices, systems, and methods for estimating oxygen saturation (SpO2).
  • Pulse oximeters typically measure and display various blood flow characteristics including but not limited to the oxygen saturation of hemoglobin in arterial blood. Oximeters pass light through blood perfused tissue such as a finger or an ear, and photoelectrically sense the absorption of light in the tissue. The amount of light absorbed is then used to calculate the amount of the blood constituent (e.g., oxyhemoglobin) being measured.
  • the light passed through the tissue is selected to be of one or more wavelengths that are absorbed by the blood in an amount representative of the amount of the blood constituent present in the blood.
  • the amount of light passed through the tissue varies in accordance with the changing amount of blood constituent in the tissue and the related light absorption.
  • the optical signal through the tissue may be degraded by noise or other artifacts.
  • One source of noise is ambient light which reaches the light detector.
  • Another source of noise is electromagnetic coupling from other electronic instruments.
  • Motion of the patient may also introduce noise onto the signal.
  • the contact between the detector and the skin, or the emitter and the skin may be temporarily disrupted when motion causes either to move away from the skin.
  • blood since blood is a fluid, it may respond differently than the surrounding tissue to inertial effects, thus resulting in momentary changes in volume at the point to which the oximeter probe is attached. Noise may degrade a pulse oximetry signal relied upon by a physician, without the physician's awareness.
  • An oximetry algorithm may include noise metrics to allow it to quantify artifacts that degrade SpO2 accuracy. These metrics may be used to reduce the impact of the artifacts on accuracy. For example, these metrics may be used to adapt multiple internal filters and to select from multiple internal SpO2 estimates. In cases with motion artifact, the noise metrics generally have a significant negative correlation to the SpO2 level for typical sensor designs. That is, increased motion artifact usually results in a negative SpO2 bias for transmission-type sensor designs, rather than simply generating zero-mean SpO2 errors.
  • the present disclosure relates, according to some embodiments, to devices, systems, and methods for estimating a physiological parameter in the presence of noise.
  • the disclosure relates, in some embodiments, to devices, systems, and methods for assessing (e.g., estimating, measuring, calculating) oxygen saturation (SpO2).
  • Methods of assessing SpO2 may include assessing a noise metric associated with motion artifact.
  • a percentage e.g., an empirically determined percentage
  • a noise metric may be simply added to the SpO2 estimate to produce a corrected SpO2 estimate. Indeed, by adding the percentage of the noise metric, the negative SpO2 bias may be corrected.
  • An oximetry algorithm may include, according to some embodiments, combining multiple internal SpO2 estimates and associated noise and/or signal quality metrics (e.g., using a radial basis neural network) to produce a modified (e.g., corrected) SpO2 estimate (e.g., rather than merely selecting the estimate from a finite number of candidates).
  • a modified SpO2 estimate may include little or no movement-based error.
  • the device may include a processor in communication with a radial basis neural network, wherein the processor is configured to receive optical oximetry data, process the optical oximetry data to produce at least one oxygen saturation estimate, communicate the at least one oxygen saturation estimate to the radial basis neural network, and receive a modified oxygen saturation estimate comprising less noise bias than the at least one oxygen saturation estimate from the radial basis neural network, wherein the radial basis neural network is configured to define the modified oxygen saturation estimate based at least in part on combining the at least one oxygen saturation estimate with corresponding signal quality metrics.
  • the device may also include a display in communication with the processor, wherein the display is operable to display a representation of the modified oxygen saturation estimate.
  • One embodiment includes a method for estimating oxygen saturation in the presence of noise.
  • the method may include determining a change in estimated oxygen saturation comprising calculating a difference between a first estimated oxygen saturation value at a first time and a second estimated oxygen saturation value at a second time. Additionally, the method may include determining a change in a saturation noise estimate between a first saturation noise estimate and a second saturation noise estimate, wherein the change in the saturation noise estimate attends the change in estimated oxygen saturation. Further, the method may include displaying either the first estimated oxygen saturation value or the second oxygen saturation value based at least in part on a comparison between the change in the estimated oxygen saturation and the change in the saturation noise estimate.
  • One embodiment includes a method for reducing affects of noise in an oxygen saturation estimate.
  • the method may include solving a computer-implemented radial basis neural network using a training set of oxygen saturation estimates, wherein the training set of oxygen saturation estimates include a target saturation estimate from noise-free data, a noise-perturbed saturation estimate, and metrics associated with calculating the noise-perturbed saturation estimate. Further, the method may include processing at least one test oxygen saturation estimate that is biased by noise using the computer-implemented radial basis neural network to produce a modified test oxygen saturation estimate comprising less noise bias than the at least one test oxygen saturation estimate.
  • One embodiment includes a method of detecting hypoxia in a subject.
  • the method may include collecting optical oximetry data from a subject, processing the optical oximetry data to produce at least one oxygen saturation estimate, processing the at least one oxygen saturation estimate using a solved, computer-implemented radial basis neural network to produce at least one modified oxygen saturation estimate, and comparing the at least one modified oxygen saturation estimate with a hypoxia oxygen saturation threshold.
  • FIG. 1 illustrates a transfer function F(x) according to a specific example embodiment of the disclosure
  • FIG. 2 illustrates oxygen saturation in a non-moving finger as a function of time using an example algorithm
  • FIG. 3 illustrates oxygen saturation from moving fingers as a function of time using a prior art algorithm
  • FIG. 4 illustrates oxygen saturation (SpO2) from moving fingers as a function of time using a prior art algorithm
  • FIG. 5 illustrates oxygen saturation from moving fingers as a function of time using an algorithm according to a specific example embodiment of the disclosure
  • FIG. 6 illustrates oxygen saturation from moving fingers as a function of time using an algorithm according to a specific example embodiment of the disclosure.
  • FIG. 7 is a block diagram of a device or system for estimating oxygen saturation in the presence of motion artifact.
  • the present disclosure relates, in some embodiments, to devices, systems, and methods for reducing the effects of motion artifact and/or noise on a physiological parameter measurement.
  • methods according to some embodiments of the disclosure may take into account physical limitations on various physiological parameters being monitored when weighting and combining a series of samples and/or measurements. Varying weights may be assigned different measurements. Optionally, measurements may be rejected if unduly corrupt. The averaging period may be adjusted according to the reliability of the measurements in some embodiments. More specifically, a general class of filters may be employed in processing the measurements. These filters may use mathematical models which (a) describe how physiological parameters change in time and/or (b) how these parameters relate to measurement in a noisy environment. The filters may adaptively modify a set of averaging weights and/or averaging times to optimally estimate the physiological parameters.
  • a method for estimating SpO2 may comprise (a) collecting optical oximetry data from a subject, and (b) processing the data to produce at least one oxygen saturation estimate.
  • a system in some embodiments, may take the natural logarithm of the optical oximetry data and then bandpass filter the data to get absorption-like data.
  • a bandpass filter may strongly attenuate data below about 0.5 Hz and/or above about 10 Hz in an attempt to remove out-of-band noise. In other words, the bandpass filter may be designed to filter out frequencies that are not due to a cardiac pulse.
  • This filtered data may be processed through two algorithms: a rate calculator and a saturation calculator.
  • an oxygen saturation estimate may be further processed to produce a modified saturation estimate (e.g., corrected to improve accuracy).
  • This post-processing may comprise modifying the saturation estimate using a radial basis neural network.
  • the present disclosure relates to methods for reducing noise effects in a system for measuring a physiological parameter.
  • a plurality of measurements may be generated corresponding to at least one wavelength of electromagnetic energy transmitted through living tissue. Selected measurements may be compared with at least one expected measurement characteristic.
  • a variable weight may be assigned to each of the selected measurements based on the comparison, thereby generating a plurality of differently weighted measurements for each wavelength.
  • a first number of weighted measurements may be averaged to obtain a filtered measurement, the first number varying according to the manner in which weights are assigned to a plurality of successive weighted measurements.
  • a plurality of filtered measurements may be generated for each wavelength.
  • the filtered measurements for each wavelength may be then combined and calculations resulting therefrom may be adaptively filtered using variable weights based on comparing the calculations to an expected calculation.
  • a second number of the weighted calculations may be averaged to obtain a filtered calculation, the second number varying according to the manner in which weights are assigned to a plurality of successive weighted calculations.
  • a method for estimating SpO2 may comprise post-processing and/or alarm handling strategies that include using a saturation noise estimate.
  • the saturation noise estimate may be used to correct the saturation estimate.
  • Empirical data includes a number of events where the saturation estimate and saturation noise estimate are mirror images of each other. There are entire cases where the saturation estimate and saturation noise estimate are fairly well correlated. Thus, where this is true across a population, a corrected saturation estimate may be calculated as follows:
  • a saturation noise estimate may be used to determine whether a change in the saturation estimate is significant. For example, if the saturation noise estimate is high, small changes in the saturation estimate (e.g., changes that are within the range of the saturation noise estimate) may not be significant. On the other hand, if the saturation noise estimate is small, large changes in the saturation estimate (e.g., changes that exceed the saturation noise estimate) may be significant. For example, it may be desirable to show only those changes in the saturation estimate that have a 95% probability of being real.
  • points may refer to percent saturation, as oxygen saturation is usually displayed in increments of 1%, from 0% to 100%. In other words, a single point may correspond to 1% oxygen saturation. It may not be desirable to continue to hold a saturation estimate if the saturation noise estimate drops below what it had been when the initial decision was made to hold it.
  • a saturation noise estimate may be used, according to some embodiments, to determine the likelihood that the saturation has crossed the saturation threshold and exceeded some sat-time integral.
  • the sat-time integral may be defined as the difference between the oxygen saturation value and its associated alarm threshold, integrated over time, such as is described in U.S. Pat. No. 5,865,736, which is incorporated herein by reference.
  • an alarm may only be raised when the likelihood that the alarm condition has occurred is high.
  • a saturation of 80 may mean that there is a 50% probability that the saturation is above 80 and a 50% probability that it is below 80. In the absence of motion, with a 2 point accuracy spec, the probability that that saturation estimate of 80 is below an alarm threshold of 85 is virtually 100%.
  • the saturation noise estimate is 5 points, then that saturation is only one standard deviation below the threshold of 85. There is about an 84 percent probability that the saturation is below the threshold and a 16 percent probability that it is above the threshold. So, a sat-time integral may be incremented only about 68 percent (i.e., the difference between 84% and 16%) as much as in the no motion case. The noisier the saturation estimate becomes, and the closer the saturation estimate gets to the alarm threshold, the slower this alarm integral may be incremented. The alarm integral may be decremented when the saturation estimate goes above the threshold, as this may indicate a greater probability that the physiological saturation is above rather than below the alarm threshold.
  • the alarm integral may be set to zero when the saturation value indicates a very high probability that the physiological saturation is above the alarm threshold (i.e. the difference between the saturation value and the alarm threshold is more than one or two times the estimate of saturation accuracy derived from one or more signal quality or noise metrics indicative of saturation accuracy).
  • a radial basis neural network may be used for saturation estimation, according to some embodiments.
  • a radial basis network may be a computer-implemented, single-hidden-layer network with n unique inputs, m nodes in the hidden layer, and a single linear node in the output layer.
  • a unique feature of a radial basis network in some embodiments, may be the equation used by its hidden nodes. Each node (a) may receive n inputs, (b) may have n associated weights, and (c) may have a single bias.
  • the equation for a node's output may be:
  • a transfer function F(x) may only have a high output when its input is near zero. This is shown in FIG. 1 . Therefore, the output of each radial basis neuron may be a function of how close x 1 . . . x n is to the ideal inputs w 1 . . . w n , and the size of the region within which the neuron will output a high value (>0.5) may be proportional to b.
  • the output of a radial basis network therefore, may be a function of how close x 1 . . . x n is to the weights for each hidden neuron.
  • Matrix-based methods for solving (determining) the weights for radial basis networks are known to those skilled in the art of signal processing, and are implemented in commonly used signal processing software, such as MATLAB®.
  • a radial-basis neural network with the previously determined weights may then be implemented in system or apparatus embodying the invention.
  • one skilled in the art of signal processing may define the value of b, thus controlling how much overlap there will be between the areas in the input space that each neuron is sensitive to. More overlap may be chosen to maximize coverage of the input space with a small network. Less overlap may be chosen if the neural net needed to respond very differently within narrow regions of the input space.
  • the ability to solve the weights for a radial basis neural network may pose an advantage over designs that apply some other types of neural networks to the same problem (i.e. accurate estimation of oxygen saturation). For example, training a backpropagation network may undesirably require that its initial weights be initialized to random values, followed by iterative training. Consequently, the weights may converge to very different values in different training sessions on the same data, and these different values may significantly impact the performance of the neural net. In contrast a radial basis neural network presented with the same training data will always have the same weights determined, making its behavior repeatable.
  • An engineer skilled in the art of signal processing may observe the weights of each hidden neuron to determine the region of data to which a neuron is sensitive.
  • a backpropagation network in some embodiments, may require a monotonic transfer function that causes a neuron's output to be high over one half of the input space, low in the other half, and changing rapidly over limited region near the boundary.
  • a radial basis network may be desirable for input spaces containing many local peaks and a backpropagation network may be desirable for inputs containing many edges.
  • a radial basis neural network may receive an input vector that is well outside of the network's training set. If this circumstance is not otherwise addressed, the output of all the hidden nodes may be virtually zero, and the output layer may output whatever its default value (bias) is. However, where an outlying data set is identified (e.g., by a skilled engineer), an additional node centered on this data may be specified and added to the hidden layer. This outlying data may then be added to the original training set. The output layer, which is always linear, may then be re-solved for the expanded training set, without having to redefine the rest of the hidden layer.
  • a method may comprise initiating or ending an alarm when an oxygen saturation estimate is below a threshold (e.g., below about 80%) or above a threshold (e.g., above about 95%).
  • a method may comprise initiating or ending an alarm when an oxygen saturation estimate is within a selected and/or desired range.
  • An alarm may include any type of alert that may be perceived by a human and/or a machine including, for example, a visible alert, an audible alert, and/or a tactile alert.
  • Non-limiting examples of a visible alert include the appearance or change in appearance of text or an icon on a display and/or a blinking light (e.g., a red light).
  • Non-limiting examples of an audible alert include an intermittent and/or continuous tone and/or playback of a recorded warning message.
  • Non-limiting examples of a tactile alert include vibration or other movement of an object.
  • Embodiments of the present invention may facilitate calculation of corrected saturation values when only a single saturation value is available along with associated signal quality metrics. Further, when multiple saturation estimates and associated signal quality metrics are combined, a corrected saturation value calculated in accordance with present embodiments may not have a value between the multiple saturation estimates. Indeed, present embodiments may result in extrapolation as well as interpolation. It should be noted that, in accordance with present embodiments, the combining of saturation values and associated metrics to reduce bias may be specific to a sensor model or design. For example, combination for a reflectance sensor may differ from combination for a transmission sensor.
  • Desaturation (hypoxia) was induced by reducing the amount of oxygen in the breathing gas delivered to the subject.
  • Normoxia resaturation
  • Oxygen saturation was assessed on a moving and a non-moving hand.
  • FIG. 3 shows the non-pulse-synchronous ensemble averaging saturation computed using the data concatenated from the four moving fingers. Saturations below 80 are calculated (detected) for only 11 of the 20 desaturations.
  • FIG. 4 shows the pulse-synchronous ensemble averaging saturation computed using the same motion data. Saturations below 80 are calculated for only 12 of the 20 desaturations. Different desaturations were missed by the two saturation estimates.
  • each “error” estimate quantifies the degree to which processed two-wavelength photoplethysmographs used to calculate a given saturation estimate are not correlated, and each “saturation age” estimate quantifies the cumulative filtering (i.e. averaging occurring in each of multiple filtering steps) used in the calculation of a given saturation estimate.
  • One of the “saturation age” estimates is highly correlated with its corresponding error estimates, but the other is not. It has been observed that these age or error estimates are sometimes mirror images of the saturation trace. Therefore, the prior art algorithm computes five semi-independent inputs that are correlated to saturations: 2 saturations, 2 error estimates, and 2 rates.
  • a Radial Basis network with 30 hidden neurons was solved for this training set of 5 inputs X 20 desaturation events, wherein the solving comprised determining the set of weights that minimizes the squared difference between the output of the Radial Basis neural net and the non-motion saturation estimate.
  • the output is shown in FIG. 5 .
  • FIG. 5 shows the output of the neural net after further processing by a 20-second median filter to mitigate this noise.
  • FIGS. 3 and 4 show two estimates of oxygen saturation calculated with different signal processing, each of which demonstrate unique perturbations due to motion artifact.
  • each trace of oxygen saturations is associated with metrics unique to that trace.
  • FIG. 5 demonstrates that by combining multiple estimates of oxygen saturation, such as those shown in FIGS. 3 and 4 , and signal quality metrics uniquely associated with each of the multiple estimates, that embodiments of the present invention may produce a corrected saturation estimate that is significantly more accurate than either of the input saturations by itself This is apparent when the trace in FIG. 2 is compared with the traces represented in FIGS. 3-5 .
  • FIG. 7 is a block diagram of a device or system for estimating oxygen saturation in the presence of motion artifact.
  • the system is generally indicated by reference number 600 .
  • the system 600 e.g., a computer
  • the system 600 includes a processor 602 (e.g., a microprocessor) and a radial basis neural network 604 in communication with one another.
  • the system 600 also includes a display 606 in communication with the processor 602 .
  • the radial basis neural network 602 may include a single-hidden-layer network and may be computer-implemented. Indeed, the radial basis neural network 604 may be stored within a memory component of the processor 602 .
  • the processor 602 may be configured to receive optical oximetry data from a patient or from another source (e.g., a database). Further, the processor 602 may be configured to process the optical oximetry data to produce at least one oxygen saturation estimate. Once the processor 602 has produced the at least one oxygen saturation estimate, the processor 602 may be configured to communicate the at least one oxygen saturation estimate to the radial basis neural network 604 .
  • the radial basis neural network 604 may be configured to modify the oxygen saturation estimate received from the processor 602 to reduce bias caused by motion artifact.
  • the radial basis neural network 604 may be configured to define a modified oxygen saturation estimate based on combining the at least one oxygen saturation estimate received from the processor 602 with one or more additional oxygen saturation estimates and corresponding signal quality metrics to reduce motion artifact bias.
  • the radial basis neural network 604 may be configured to communicate the modified oxygen saturation estimate to the processor 602 , which may be configured to present a representation of the modified oxygen saturation estimate on the display 606 , which may be communicatively coupled with the processor 602 .
  • a device and/or system may include a single computer-implemented network per subject (e.g., patient) or a single computer-implemented network that serves two or more subjects.
  • the size of a device and/or system may be scaled up or down to suit the needs and/or desires of a practitioner.
  • the disclosed endpoints may be treated as exact and/or approximations as desired or demanded by the particular embodiment.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Artificial Intelligence (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Optics & Photonics (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Physics (AREA)
  • Fuzzy Systems (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

The present disclosure relates, according to some embodiments, to devices, systems, and methods for estimating a physiological parameter in the presence of noise. For example, the disclosure relates, in some embodiments, to devices, systems, and methods for assessing (e.g., estimating, measuring, calculating) oxygen saturation (SpO2). Methods of assessing SpO2 may include assessing a noise metric associated with motion artifact. In some embodiments, a percentage (e.g., an empirically determined percentage) of a noise metric may be simply added to the SpO2 estimate to produce a corrected SpO2 estimate. An oximetry algorithm may include, according to some embodiments, combining multiple internal SpO2 estimates and associated noise and/or signal quality metrics (e.g., using a radial basis neural network) to produce a modified (e.g., corrected) SpO2 estimate (e.g., rather than merely selecting the estimate from a finite number of candidates). A modified SpO2 estimate may include little or no movement-based error.

Description

  • This application is a divisional of U.S. patent application Ser. No. 12/165,256, entitled “Method for Improved Oxygen Saturation Estimation in the Presence of Noise,” filed Jul. 3, 2008, which is incorporated herein by reference in its entirety.
  • FIELD OF THE DISCLOSURE
  • The present disclosure relates, in some embodiments, to devices, systems, and methods for estimating oxygen saturation (SpO2).
  • BACKGROUND OF THE DISCLOSURE
  • Pulse oximeters typically measure and display various blood flow characteristics including but not limited to the oxygen saturation of hemoglobin in arterial blood. Oximeters pass light through blood perfused tissue such as a finger or an ear, and photoelectrically sense the absorption of light in the tissue. The amount of light absorbed is then used to calculate the amount of the blood constituent (e.g., oxyhemoglobin) being measured.
  • The light passed through the tissue is selected to be of one or more wavelengths that are absorbed by the blood in an amount representative of the amount of the blood constituent present in the blood. The amount of light passed through the tissue varies in accordance with the changing amount of blood constituent in the tissue and the related light absorption.
  • The optical signal through the tissue may be degraded by noise or other artifacts. One source of noise is ambient light which reaches the light detector. Another source of noise is electromagnetic coupling from other electronic instruments. Motion of the patient may also introduce noise onto the signal. For example, the contact between the detector and the skin, or the emitter and the skin, may be temporarily disrupted when motion causes either to move away from the skin. In addition, since blood is a fluid, it may respond differently than the surrounding tissue to inertial effects, thus resulting in momentary changes in volume at the point to which the oximeter probe is attached. Noise may degrade a pulse oximetry signal relied upon by a physician, without the physician's awareness.
  • An oximetry algorithm may include noise metrics to allow it to quantify artifacts that degrade SpO2 accuracy. These metrics may be used to reduce the impact of the artifacts on accuracy. For example, these metrics may be used to adapt multiple internal filters and to select from multiple internal SpO2 estimates. In cases with motion artifact, the noise metrics generally have a significant negative correlation to the SpO2 level for typical sensor designs. That is, increased motion artifact usually results in a negative SpO2 bias for transmission-type sensor designs, rather than simply generating zero-mean SpO2 errors.
  • SUMMARY
  • Accordingly, a need has arisen for improved noise metrics to correct negative SpO2 bias associated with artifacts (e.g., movement artifacts).
  • The present disclosure relates, according to some embodiments, to devices, systems, and methods for estimating a physiological parameter in the presence of noise. For example, the disclosure relates, in some embodiments, to devices, systems, and methods for assessing (e.g., estimating, measuring, calculating) oxygen saturation (SpO2). Methods of assessing SpO2 may include assessing a noise metric associated with motion artifact. In some embodiments, a percentage (e.g., an empirically determined percentage) of a noise metric may be simply added to the SpO2 estimate to produce a corrected SpO2 estimate. Indeed, by adding the percentage of the noise metric, the negative SpO2 bias may be corrected. An oximetry algorithm may include, according to some embodiments, combining multiple internal SpO2 estimates and associated noise and/or signal quality metrics (e.g., using a radial basis neural network) to produce a modified (e.g., corrected) SpO2 estimate (e.g., rather than merely selecting the estimate from a finite number of candidates). A modified SpO2 estimate may include little or no movement-based error.
  • One embodiment includes a device for estimating oxygen saturation in the presence of noise. The device may include a processor in communication with a radial basis neural network, wherein the processor is configured to receive optical oximetry data, process the optical oximetry data to produce at least one oxygen saturation estimate, communicate the at least one oxygen saturation estimate to the radial basis neural network, and receive a modified oxygen saturation estimate comprising less noise bias than the at least one oxygen saturation estimate from the radial basis neural network, wherein the radial basis neural network is configured to define the modified oxygen saturation estimate based at least in part on combining the at least one oxygen saturation estimate with corresponding signal quality metrics. The device may also include a display in communication with the processor, wherein the display is operable to display a representation of the modified oxygen saturation estimate.
  • One embodiment includes a method for estimating oxygen saturation in the presence of noise. The method may include determining a change in estimated oxygen saturation comprising calculating a difference between a first estimated oxygen saturation value at a first time and a second estimated oxygen saturation value at a second time. Additionally, the method may include determining a change in a saturation noise estimate between a first saturation noise estimate and a second saturation noise estimate, wherein the change in the saturation noise estimate attends the change in estimated oxygen saturation. Further, the method may include displaying either the first estimated oxygen saturation value or the second oxygen saturation value based at least in part on a comparison between the change in the estimated oxygen saturation and the change in the saturation noise estimate.
  • One embodiment includes a method for reducing affects of noise in an oxygen saturation estimate. The method may include solving a computer-implemented radial basis neural network using a training set of oxygen saturation estimates, wherein the training set of oxygen saturation estimates include a target saturation estimate from noise-free data, a noise-perturbed saturation estimate, and metrics associated with calculating the noise-perturbed saturation estimate. Further, the method may include processing at least one test oxygen saturation estimate that is biased by noise using the computer-implemented radial basis neural network to produce a modified test oxygen saturation estimate comprising less noise bias than the at least one test oxygen saturation estimate.
  • One embodiment includes a method of detecting hypoxia in a subject. The method may include collecting optical oximetry data from a subject, processing the optical oximetry data to produce at least one oxygen saturation estimate, processing the at least one oxygen saturation estimate using a solved, computer-implemented radial basis neural network to produce at least one modified oxygen saturation estimate, and comparing the at least one modified oxygen saturation estimate with a hypoxia oxygen saturation threshold.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • Some embodiments of the disclosure may be understood by referring, in part, to the present disclosure and the accompanying drawings, wherein:
  • FIG. 1 illustrates a transfer function F(x) according to a specific example embodiment of the disclosure;
  • FIG. 2 illustrates oxygen saturation in a non-moving finger as a function of time using an example algorithm;
  • FIG. 3 illustrates oxygen saturation from moving fingers as a function of time using a prior art algorithm;
  • FIG. 4 illustrates oxygen saturation (SpO2) from moving fingers as a function of time using a prior art algorithm;
  • FIG. 5 illustrates oxygen saturation from moving fingers as a function of time using an algorithm according to a specific example embodiment of the disclosure;
  • FIG. 6 illustrates oxygen saturation from moving fingers as a function of time using an algorithm according to a specific example embodiment of the disclosure; and
  • FIG. 7 is a block diagram of a device or system for estimating oxygen saturation in the presence of motion artifact.
  • DETAILED DESCRIPTION
  • The present disclosure relates, in some embodiments, to devices, systems, and methods for reducing the effects of motion artifact and/or noise on a physiological parameter measurement. For example, methods according to some embodiments of the disclosure may take into account physical limitations on various physiological parameters being monitored when weighting and combining a series of samples and/or measurements. Varying weights may be assigned different measurements. Optionally, measurements may be rejected if unduly corrupt. The averaging period may be adjusted according to the reliability of the measurements in some embodiments. More specifically, a general class of filters may be employed in processing the measurements. These filters may use mathematical models which (a) describe how physiological parameters change in time and/or (b) how these parameters relate to measurement in a noisy environment. The filters may adaptively modify a set of averaging weights and/or averaging times to optimally estimate the physiological parameters.
  • The present disclosure relates, in some embodiments, to devices, systems, and methods for estimating oxygen saturation (SpO2). For example, a method for estimating SpO2 may comprise (a) collecting optical oximetry data from a subject, and (b) processing the data to produce at least one oxygen saturation estimate. A system, in some embodiments, may take the natural logarithm of the optical oximetry data and then bandpass filter the data to get absorption-like data. A bandpass filter may strongly attenuate data below about 0.5 Hz and/or above about 10 Hz in an attempt to remove out-of-band noise. In other words, the bandpass filter may be designed to filter out frequencies that are not due to a cardiac pulse. This filtered data may be processed through two algorithms: a rate calculator and a saturation calculator. In some embodiments, an oxygen saturation estimate may be further processed to produce a modified saturation estimate (e.g., corrected to improve accuracy). This post-processing may comprise modifying the saturation estimate using a radial basis neural network.
  • According to some embodiments, the present disclosure relates to methods for reducing noise effects in a system for measuring a physiological parameter. A plurality of measurements may be generated corresponding to at least one wavelength of electromagnetic energy transmitted through living tissue. Selected measurements may be compared with at least one expected measurement characteristic. A variable weight may be assigned to each of the selected measurements based on the comparison, thereby generating a plurality of differently weighted measurements for each wavelength. A first number of weighted measurements may be averaged to obtain a filtered measurement, the first number varying according to the manner in which weights are assigned to a plurality of successive weighted measurements. A plurality of filtered measurements may be generated for each wavelength. The filtered measurements for each wavelength may be then combined and calculations resulting therefrom may be adaptively filtered using variable weights based on comparing the calculations to an expected calculation. A second number of the weighted calculations may be averaged to obtain a filtered calculation, the second number varying according to the manner in which weights are assigned to a plurality of successive weighted calculations.
  • According to some embodiments, a method for estimating SpO2 may comprise post-processing and/or alarm handling strategies that include using a saturation noise estimate. For example, the saturation noise estimate may be used to correct the saturation estimate. Empirical data includes a number of events where the saturation estimate and saturation noise estimate are mirror images of each other. There are entire cases where the saturation estimate and saturation noise estimate are fairly well correlated. Thus, where this is true across a population, a corrected saturation estimate may be calculated as follows:

  • Sat′=Sat+k(max(noise−m,0)),
  • where k and m are determined empirically so as to minimize the saturation error.
  • In some embodiments, a saturation noise estimate may be used to determine whether a change in the saturation estimate is significant. For example, if the saturation noise estimate is high, small changes in the saturation estimate (e.g., changes that are within the range of the saturation noise estimate) may not be significant. On the other hand, if the saturation noise estimate is small, large changes in the saturation estimate (e.g., changes that exceed the saturation noise estimate) may be significant. For example, it may be desirable to show only those changes in the saturation estimate that have a 95% probability of being real. If a change in the saturation estimate is attended by a change in the saturation noise estimate from 2 points to 12 points, then the saturation estimate may be held unless the change in the saturation estimate changed by twice that delta (i.e., 2×(12 points−2 points)=20 points). It should be noted that the term “points” may refer to percent saturation, as oxygen saturation is usually displayed in increments of 1%, from 0% to 100%. In other words, a single point may correspond to 1% oxygen saturation. It may not be desirable to continue to hold a saturation estimate if the saturation noise estimate drops below what it had been when the initial decision was made to hold it. For example, if a saturation estimate changes by 15 points, but the saturation noise estimate increases from 2 points to 12 points, the saturation estimate may be held (i.e., 15<20). However, if the saturation noise estimate subsequently drops to 7 points (now, 2×delta=10 points), then the 15-point change in the saturation estimate is more likely to be real. Thus, it may be desirable to discontinue holding the saturation estimate. Alternatively, an assumption may be made that, on average, a held saturation estimate becomes progressively less accurate at a rate of perhaps 0.2 saturation points per second. In such a case, if the saturation noise estimate goes from 2 points to 12 points, the saturation noise estimate may be held for up to 50 seconds (i.e., 10 points/0.2 points/second=50 seconds).
  • A saturation noise estimate may be used, according to some embodiments, to determine the likelihood that the saturation has crossed the saturation threshold and exceeded some sat-time integral. The sat-time integral may be defined as the difference between the oxygen saturation value and its associated alarm threshold, integrated over time, such as is described in U.S. Pat. No. 5,865,736, which is incorporated herein by reference. During motion, an alarm may only be raised when the likelihood that the alarm condition has occurred is high. For example, a saturation of 80 may mean that there is a 50% probability that the saturation is above 80 and a 50% probability that it is below 80. In the absence of motion, with a 2 point accuracy spec, the probability that that saturation estimate of 80 is below an alarm threshold of 85 is virtually 100%. But if the saturation noise estimate is 5 points, then that saturation is only one standard deviation below the threshold of 85. There is about an 84 percent probability that the saturation is below the threshold and a 16 percent probability that it is above the threshold. So, a sat-time integral may be incremented only about 68 percent (i.e., the difference between 84% and 16%) as much as in the no motion case. The noisier the saturation estimate becomes, and the closer the saturation estimate gets to the alarm threshold, the slower this alarm integral may be incremented. The alarm integral may be decremented when the saturation estimate goes above the threshold, as this may indicate a greater probability that the physiological saturation is above rather than below the alarm threshold. Furthermore, the alarm integral may be set to zero when the saturation value indicates a very high probability that the physiological saturation is above the alarm threshold (i.e. the difference between the saturation value and the alarm threshold is more than one or two times the estimate of saturation accuracy derived from one or more signal quality or noise metrics indicative of saturation accuracy).
  • A radial basis neural network may be used for saturation estimation, according to some embodiments. A radial basis network may be a computer-implemented, single-hidden-layer network with n unique inputs, m nodes in the hidden layer, and a single linear node in the output layer. A unique feature of a radial basis network, in some embodiments, may be the equation used by its hidden nodes. Each node (a) may receive n inputs, (b) may have n associated weights, and (c) may have a single bias. The equation for a node's output may be:

  • y=F(b*distance(x 1 . . . x n , w 1 . . . w n)),
    • where x1 . . . xn denotes the inputs,
    • w1 . . . wn denotes the weights,
    • b is the bias,
    • distance ( ) is the distance between the x and w vectors in an n-dimensional space, and

  • F(x)=e −x 2 .
  • In some embodiments, a transfer function F(x) may only have a high output when its input is near zero. This is shown in FIG. 1. Therefore, the output of each radial basis neuron may be a function of how close x1 . . . xn is to the ideal inputs w1 . . . wn, and the size of the region within which the neuron will output a high value (>0.5) may be proportional to b. The output of a radial basis network, therefore, may be a function of how close x1 . . . xn is to the weights for each hidden neuron. Matrix-based methods for solving (determining) the weights for radial basis networks are known to those skilled in the art of signal processing, and are implemented in commonly used signal processing software, such as MATLAB®. A radial-basis neural network with the previously determined weights may then be implemented in system or apparatus embodying the invention.
  • Depending on the desired behavior, one skilled in the art of signal processing may define the value of b, thus controlling how much overlap there will be between the areas in the input space that each neuron is sensitive to. More overlap may be chosen to maximize coverage of the input space with a small network. Less overlap may be chosen if the neural net needed to respond very differently within narrow regions of the input space.
  • The ability to solve the weights for a radial basis neural network may pose an advantage over designs that apply some other types of neural networks to the same problem (i.e. accurate estimation of oxygen saturation). For example, training a backpropagation network may undesirably require that its initial weights be initialized to random values, followed by iterative training. Consequently, the weights may converge to very different values in different training sessions on the same data, and these different values may significantly impact the performance of the neural net. In contrast a radial basis neural network presented with the same training data will always have the same weights determined, making its behavior repeatable.
  • An engineer skilled in the art of signal processing, in some embodiments, may observe the weights of each hidden neuron to determine the region of data to which a neuron is sensitive. A backpropagation network, in some embodiments, may require a monotonic transfer function that causes a neuron's output to be high over one half of the input space, low in the other half, and changing rapidly over limited region near the boundary. Thus, in some embodiments, a radial basis network may be desirable for input spaces containing many local peaks and a backpropagation network may be desirable for inputs containing many edges.
  • In some embodiments, a radial basis neural network may receive an input vector that is well outside of the network's training set. If this circumstance is not otherwise addressed, the output of all the hidden nodes may be virtually zero, and the output layer may output whatever its default value (bias) is. However, where an outlying data set is identified (e.g., by a skilled engineer), an additional node centered on this data may be specified and added to the hidden layer. This outlying data may then be added to the original training set. The output layer, which is always linear, may then be re-solved for the expanded training set, without having to redefine the rest of the hidden layer.
  • The present disclosure relates, according to some embodiments, to devices, systems, and methods for detecting oxygen saturation that crosses a threshold and/or falls within a range. For example, a method may comprise initiating or ending an alarm when an oxygen saturation estimate is below a threshold (e.g., below about 80%) or above a threshold (e.g., above about 95%). A method may comprise initiating or ending an alarm when an oxygen saturation estimate is within a selected and/or desired range. An alarm may include any type of alert that may be perceived by a human and/or a machine including, for example, a visible alert, an audible alert, and/or a tactile alert. Non-limiting examples of a visible alert include the appearance or change in appearance of text or an icon on a display and/or a blinking light (e.g., a red light). Non-limiting examples of an audible alert include an intermittent and/or continuous tone and/or playback of a recorded warning message. Non-limiting examples of a tactile alert include vibration or other movement of an object.
  • Embodiments of the present invention may facilitate calculation of corrected saturation values when only a single saturation value is available along with associated signal quality metrics. Further, when multiple saturation estimates and associated signal quality metrics are combined, a corrected saturation value calculated in accordance with present embodiments may not have a value between the multiple saturation estimates. Indeed, present embodiments may result in extrapolation as well as interpolation. It should be noted that, in accordance with present embodiments, the combining of saturation values and associated metrics to reduce bias may be specific to a sensor model or design. For example, combination for a reflectance sensor may differ from combination for a transmission sensor.
  • Some specific example embodiments of the disclosure may be illustrated by one or more of the examples provided herein.
  • EXAMPLE 1 Breathdown Procedure
  • Desaturation (hypoxia) was induced by reducing the amount of oxygen in the breathing gas delivered to the subject. Normoxia (resaturation) was restored by reverting to the original breathing gas mixture. Oxygen saturation was assessed on a moving and a non-moving hand.
  • Five desaturations and resaturations were caused in quick succession. Simultaneous data was collected from the four fingers of the moving hand and from one finger of the non-moving hand. The non-motion saturations were concatenated from each finger of the non-moving hand to form a target vector containing twenty desaturations (FIG. 2), each of which has a minima below 80 percent. The saturations computed from the moving hand by a prior-art algorithm or an algorithm embodying the present invention may likewise be concatenated (FIGS. 3-6) to form a test set. Thus, an ideal algorithm should report twenty desaturations below 80 percent on both the non-motion and motion data.
  • EXAMPLE 2 Prior Art Algorithm with Pulse-Synchronous Ensemble Averaging
  • A prior art algorithm that computes saturation with and without pulse-synchronous ensemble averaging was tested. FIG. 3 shows the non-pulse-synchronous ensemble averaging saturation computed using the data concatenated from the four moving fingers. Saturations below 80 are calculated (detected) for only 11 of the 20 desaturations. FIG. 4 shows the pulse-synchronous ensemble averaging saturation computed using the same motion data. Saturations below 80 are calculated for only 12 of the 20 desaturations. Different desaturations were missed by the two saturation estimates.
  • In addition to these two saturation estimates, the same prior art algorithm also computes two associated error and “saturation age” estimates, wherein each “error” estimate quantifies the degree to which processed two-wavelength photoplethysmographs used to calculate a given saturation estimate are not correlated, and each “saturation age” estimate quantifies the cumulative filtering (i.e. averaging occurring in each of multiple filtering steps) used in the calculation of a given saturation estimate. One of the “saturation age” estimates is highly correlated with its corresponding error estimates, but the other is not. It has been observed that these age or error estimates are sometimes mirror images of the saturation trace. Therefore, the prior art algorithm computes five semi-independent inputs that are correlated to saturations: 2 saturations, 2 error estimates, and 2 rates.
  • EXAMPLE 3 A Radial Basis Network
  • A Radial Basis network with 30 hidden neurons was solved for this training set of 5 inputs X 20 desaturation events, wherein the solving comprised determining the set of weights that minimizes the squared difference between the output of the Radial Basis neural net and the non-motion saturation estimate. The output is shown in FIG. 5.
  • Although the tracing in FIG. 5 contains some short-term noise, the network successfully estimated all 20 desaturations below 80. FIG. 6 shows the output of the neural net after further processing by a 20-second median filter to mitigate this noise.
  • FIGS. 3 and 4 show two estimates of oxygen saturation calculated with different signal processing, each of which demonstrate unique perturbations due to motion artifact. As the text describes, each trace of oxygen saturations is associated with metrics unique to that trace. FIG. 5 demonstrates that by combining multiple estimates of oxygen saturation, such as those shown in FIGS. 3 and 4, and signal quality metrics uniquely associated with each of the multiple estimates, that embodiments of the present invention may produce a corrected saturation estimate that is significantly more accurate than either of the input saturations by itself This is apparent when the trace in FIG. 2 is compared with the traces represented in FIGS. 3-5.
  • FIG. 7 is a block diagram of a device or system for estimating oxygen saturation in the presence of motion artifact. The system is generally indicated by reference number 600. The system 600 (e.g., a computer) includes a processor 602 (e.g., a microprocessor) and a radial basis neural network 604 in communication with one another. The system 600 also includes a display 606 in communication with the processor 602. The radial basis neural network 602 may include a single-hidden-layer network and may be computer-implemented. Indeed, the radial basis neural network 604 may be stored within a memory component of the processor 602.
  • The processor 602 may be configured to receive optical oximetry data from a patient or from another source (e.g., a database). Further, the processor 602 may be configured to process the optical oximetry data to produce at least one oxygen saturation estimate. Once the processor 602 has produced the at least one oxygen saturation estimate, the processor 602 may be configured to communicate the at least one oxygen saturation estimate to the radial basis neural network 604. The radial basis neural network 604 may be configured to modify the oxygen saturation estimate received from the processor 602 to reduce bias caused by motion artifact. For example, the radial basis neural network 604 may be configured to define a modified oxygen saturation estimate based on combining the at least one oxygen saturation estimate received from the processor 602 with one or more additional oxygen saturation estimates and corresponding signal quality metrics to reduce motion artifact bias. Once the modified oxygen saturation estimate has been established by the radial basis neural network 604, the radial basis neural network 604 may be configured to communicate the modified oxygen saturation estimate to the processor 602, which may be configured to present a representation of the modified oxygen saturation estimate on the display 606, which may be communicatively coupled with the processor 602.
  • As will be understood by those skilled in the art who have the benefit of the instant disclosure, other equivalent or alternative compositions, devices, methods, and systems for estimating oxygen saturation can be envisioned without departing from the description contained herein. Accordingly, the manner of carrying out the disclosure as shown and described is to be construed as illustrative only.
  • Persons skilled in the art may make various changes in the shape, size, number, and/or arrangement of parts without departing from the scope of the instant disclosure. For example, a device and/or system may include a single computer-implemented network per subject (e.g., patient) or a single computer-implemented network that serves two or more subjects. In addition, the size of a device and/or system may be scaled up or down to suit the needs and/or desires of a practitioner. Also, where ranges have been provided, the disclosed endpoints may be treated as exact and/or approximations as desired or demanded by the particular embodiment. In addition, it may be desirable in some embodiments to mix and match range endpoints. All or a portion of an oxygen saturation device and/or system may be configured and arranged to be disposable, serviceable, interchangeable, and/or replaceable. These equivalents and alternatives along with obvious changes and modifications are intended to be included within the scope of the present disclosure. Accordingly, the foregoing disclosure is intended to be illustrative, but not limiting, of the scope of the disclosure as illustrated by the following claims.

Claims (20)

What is claimed is:
1. A device for estimating oxygen saturation in the presence of noise, the device comprising:
a processor in communication with a neural network, wherein the processor is configured to receive oximetry data, process the oximetry data to produce at least one oxygen saturation estimate, communicate the at least one oxygen saturation estimate to the neural network, and receive a modified oxygen saturation estimate from the neural network comprising less noise bias than the at least one oxygen saturation estimate, wherein the neural network is configured to define the modified oxygen saturation estimate based at least in part on combining the at least one oxygen saturation estimate with corresponding signal quality metrics; and
a display in communication with the processor, wherein the display is configured to display a representation of the modified oxygen saturation estimate.
2. The device of claim 1, wherein the neural network comprises a radial basis neural network.
3. The device of claim 2, wherein the radial basis neural network comprises a single-hidden-layer network.
4. The device of claim 3, wherein the single-hidden-layer network has a plurality of nodes, and each node of the plurality of nodes includes one or more associated weights and a single bias.
5. The device of claim 4, wherein the single bias of at least one node of the plurality of nodes is in the range of about 5% oxygen saturation to about 10% oxygen saturation.
6. The device of claim 4, wherein the single bias of at least one node in the plurality of nodes is in the range of about 1% oxygen saturation to about 5% oxygen saturation.
7. The device of claim 1, wherein the neural network is configured to define the modified oxygen saturation estimate based at least in part on combining the at least one oxygen saturation estimate with one or more additional oxygen saturation estimates.
8. The device of claim 7, wherein the modified oxygen saturation estimate has a value that is not between the values of the at least one oxygen saturation estimate and the one or more additional oxygen saturation estimates.
9. The device of claim 1, wherein the processor is configured to compare the modified oxygen saturation estimate with a hypoxia oxygen saturation threshold and to initiate an alarm if the modified oxygen saturation estimate is less than or equal to the hypoxia oxygen saturation threshold.
10. The device of claim 9, wherein the alarm comprises one or more of a visible alert, an audible alert, a tactile alert, or any combination thereof
11. A method for estimating oxygen saturation in the presence of noise, the method comprising:
generating at least one oxygen saturation estimate using the processor;
communicating the at least one oxygen saturation estimate to a neural network using the processor;
receiving, at the processor, a modified oxygen saturation estimate from the neural network comprising less noise bias than the at least one oxygen saturation estimate, wherein the neural network defines the modified oxygen saturation estimate based at least in part on combining the at least one oxygen saturation estimate with corresponding signal quality metrics; and
presenting, via the processor, a representation of the modified oxygen saturation estimate on a display that is in communication with the processor.
12. The method of claim 11, wherein the neural network comprises a radial basis neural network.
13. The method of the claim 12, wherein the radial basis neural network comprises a single-hidden-layer network.
14. The method of claim 13, wherein the single-hidden-layer network has a plurality of nodes, and each node of the plurality of nodes includes one or more associated weights and a single bias.
15. The method of claim 11, wherein the neural network defines the modified oxygen saturation estimate based at least in part on combining the at least one oxygen saturation estimate with one or more additional oxygen saturation estimates.
16. The method of claim 15, wherein the modified oxygen saturation estimate has a value that is not between the values of the at least one oxygen saturation estimate and the one or more additional oxygen saturation estimates.
17. The method of claim 11, comprising comparing the modified oxygen saturation estimate with a hypoxia oxygen saturation threshold using the processor and initiating an alarm, via the processor, if the modified oxygen saturation estimate is less than or equal to the hypoxia oxygen saturation threshold.
18. The method of claim 11, comprising receiving oximetry data from a pulse oximeter at the processor and generating the at least one oxygen saturation estimate based on the received oximetry data, using the processor.
19. A system for estimating oxygen saturation in the presence of noise, the system comprising:
a processor in communication with a radial basis neural network, wherein the processor is configured to process oximetry data to produce at least one oxygen saturation estimate, the radial basis neural network is configured to receive the at least one oxygen saturation estimate from the processor and to modify the at least one oxygen saturation estimate based on corresponding signal quality metrics indicative of an accuracy of the at least one oxygen saturation estimate to define a modified oxygen saturation estimate having less noise bias than the at least one oxygen saturation, and the processor is configured to provide an output indicative of the modified oxygen saturation estimate.
20. The system of claim 19, wherein the processor is configured to compare the modified oxygen saturation estimate with a hypoxia oxygen saturation threshold and to initiate an alarm if the modified oxygen saturation estimate is less than or equal to the hypoxia oxygen saturation threshold.
US14/472,038 2008-06-30 2014-08-28 Method for improved oxygen saturation estimation in the presence of noise Abandoned US20140371554A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US14/472,038 US20140371554A1 (en) 2008-06-30 2014-08-28 Method for improved oxygen saturation estimation in the presence of noise

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US12/165,256 US8862194B2 (en) 2008-06-30 2008-06-30 Method for improved oxygen saturation estimation in the presence of noise
US14/472,038 US20140371554A1 (en) 2008-06-30 2014-08-28 Method for improved oxygen saturation estimation in the presence of noise

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US12/165,256 Continuation US8862194B2 (en) 2008-06-30 2008-06-30 Method for improved oxygen saturation estimation in the presence of noise

Publications (1)

Publication Number Publication Date
US20140371554A1 true US20140371554A1 (en) 2014-12-18

Family

ID=41448292

Family Applications (2)

Application Number Title Priority Date Filing Date
US12/165,256 Active 2033-07-31 US8862194B2 (en) 2008-06-30 2008-06-30 Method for improved oxygen saturation estimation in the presence of noise
US14/472,038 Abandoned US20140371554A1 (en) 2008-06-30 2014-08-28 Method for improved oxygen saturation estimation in the presence of noise

Family Applications Before (1)

Application Number Title Priority Date Filing Date
US12/165,256 Active 2033-07-31 US8862194B2 (en) 2008-06-30 2008-06-30 Method for improved oxygen saturation estimation in the presence of noise

Country Status (1)

Country Link
US (2) US8862194B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120123231A1 (en) * 2010-11-11 2012-05-17 O'reilly Michael Monitoring cardiac output and vessel fluid volume
US10342466B2 (en) 2015-03-24 2019-07-09 Covidien Lp Regional saturation system with ensemble averaging
US10634081B2 (en) * 2018-02-05 2020-04-28 Toyota Jidosha Kabushiki Kaisha Control device of internal combustion engine
CN118098474A (en) * 2024-04-19 2024-05-28 自贡市第一人民医院 Data optimization processing method for medical nursing clinical monitoring equipment

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5339818A (en) * 1989-09-20 1994-08-23 University Of Utah Research Foundation Method for determining blood pressure utilizing a neural network
US6371114B1 (en) * 1998-07-24 2002-04-16 Minnesota Innovative Technologies & Instruments Corporation Control device for supplying supplemental respiratory oxygen
US20050049470A1 (en) * 2003-08-27 2005-03-03 Terry Alvin Mark Multi-domain motion estimation and plethysmographic recognition using fuzzy neural-nets

Family Cites Families (102)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3638640A (en) 1967-11-01 1972-02-01 Robert F Shaw Oximeter and method for in vivo determination of oxygen saturation in blood using three or more different wavelengths
US4063824A (en) 1975-08-05 1977-12-20 E. I. Du Pont De Nemours And Company Chemical dosimeter having a constant flow air sampling pump
US4123932A (en) 1977-05-25 1978-11-07 E. I. Du Pont De Nemours And Company Dosimeter for monitoring working areas
US4714341A (en) 1984-02-23 1987-12-22 Minolta Camera Kabushiki Kaisha Multi-wavelength oximeter having a means for disregarding a poor signal
US4911167A (en) 1985-06-07 1990-03-27 Nellcor Incorporated Method and apparatus for detecting optical pulses
US4936679A (en) 1985-11-12 1990-06-26 Becton, Dickinson And Company Optical fiber transducer driving and measuring circuit and method for using same
US4805623A (en) 1987-09-04 1989-02-21 Vander Corporation Spectrophotometric method for quantitatively determining the concentration of a dilute component in a light- or other radiation-scattering environment
CA1331483C (en) 1988-11-02 1994-08-16 Britton Chance User-wearable hemoglobinometer for measuring the metabolic condition of a subject
US4972331A (en) 1989-02-06 1990-11-20 Nim, Inc. Phase modulated spectrophotometry
US5122974A (en) 1989-02-06 1992-06-16 Nim, Inc. Phase modulated spectrophotometry
US5564417A (en) 1991-01-24 1996-10-15 Non-Invasive Technology, Inc. Pathlength corrected oximeter and the like
US5119815A (en) 1988-12-21 1992-06-09 Nim, Incorporated Apparatus for determining the concentration of a tissue pigment of known absorbance, in vivo, using the decay characteristics of scintered electromagnetic radiation
US6708048B1 (en) 1989-02-06 2004-03-16 Non-Invasive Technology, Inc. Phase modulation spectrophotometric apparatus
US6183414B1 (en) 1999-04-26 2001-02-06 Michael S. Wysor Technique for restoring plasticity to tissues of a male or female organ
US6681128B2 (en) 1990-10-06 2004-01-20 Hema Metrics, Inc. System for noninvasive hematocrit monitoring
US6266546B1 (en) 1990-10-06 2001-07-24 In-Line Diagnostics Corporation System for noninvasive hematocrit monitoring
US5372136A (en) 1990-10-06 1994-12-13 Noninvasive Medical Technology Corporation System and method for noninvasive hematocrit monitoring
US6246894B1 (en) 1993-02-01 2001-06-12 In-Line Diagnostics Corporation System and method for measuring blood urea nitrogen, blood osmolarity, plasma free hemoglobin and tissue water content
CA2105682C (en) 1991-03-07 2003-09-02 Mohamed K. Diab Signal processing apparatus and method
MX9702434A (en) * 1991-03-07 1998-05-31 Masimo Corp Signal processing apparatus.
JP3720046B2 (en) 1991-05-16 2005-11-24 ノン−インヴェイシヴ テクノロジイ,インク. Instruments such as hemoglobin meters for measuring the metabolic state of subjects
US6549795B1 (en) 1991-05-16 2003-04-15 Non-Invasive Technology, Inc. Spectrophotometer for tissue examination
US5385143A (en) 1992-02-06 1995-01-31 Nihon Kohden Corporation Apparatus for measuring predetermined data of living tissue
JP2608828B2 (en) 1992-02-06 1997-05-14 日本光電工業株式会社 Non-invasive oximeter
US5297548A (en) 1992-02-07 1994-03-29 Ohmeda Inc. Arterial blood monitoring probe
US5337745A (en) 1992-03-10 1994-08-16 Benaron David A Device and method for in vivo qualitative or quantative measurement of blood chromophore concentration using blood pulse spectrophotometry
US6785568B2 (en) 1992-05-18 2004-08-31 Non-Invasive Technology Inc. Transcranial examination of the brain
US5355880A (en) 1992-07-06 1994-10-18 Sandia Corporation Reliable noninvasive measurement of blood gases
GB9216431D0 (en) 1992-08-01 1992-09-16 Univ Swansea Optical monitoring or measuring artefact suppression
US7376453B1 (en) * 1993-10-06 2008-05-20 Masimo Corporation Signal processing apparatus
US5485847A (en) 1993-10-08 1996-01-23 Nellcor Puritan Bennett Incorporated Pulse oximeter using a virtual trigger for heart rate synchronization
US5758644A (en) 1995-06-07 1998-06-02 Masimo Corporation Manual and automatic probe calibration
US5853364A (en) 1995-08-07 1998-12-29 Nellcor Puritan Bennett, Inc. Method and apparatus for estimating physiological parameters using model-based adaptive filtering
WO1997049330A1 (en) 1996-06-27 1997-12-31 Falcon Medical, Inc. Motion artifact resistant oximeter using three wavelengths
US6163715A (en) 1996-07-17 2000-12-19 Criticare Systems, Inc. Direct to digital oximeter and method for calculating oxygenation levels
US6120460A (en) 1996-09-04 2000-09-19 Abreu; Marcio Marc Method and apparatus for signal acquisition, processing and transmission for evaluation of bodily functions
US6544193B2 (en) 1996-09-04 2003-04-08 Marcio Marc Abreu Noninvasive measurement of chemical substances
US6847490B1 (en) 1997-01-13 2005-01-25 Medispectra, Inc. Optical probe accessory device for use in vivo diagnostic procedures
US6487439B1 (en) 1997-03-17 2002-11-26 Victor N. Skladnev Glove-mounted hybrid probe for tissue type recognition
EP1006863B1 (en) 1997-03-21 2003-10-15 Nellcor Puritan Bennett Incorporated Method and apparatus for adaptively averaging data signals
US20080287756A1 (en) * 1997-07-14 2008-11-20 Lynn Lawrence A Pulse oximetry relational alarm system for early recognition of instability and catastrophic occurrences
JP2002501803A (en) 1998-02-05 2002-01-22 イン−ラインダイアグノスティックスコーポレイション Non-invasive blood component monitoring method and apparatus
US6662030B2 (en) 1998-05-18 2003-12-09 Abbott Laboratories Non-invasive sensor having controllable temperature feature
ATE521277T1 (en) * 1998-06-03 2011-09-15 Masimo Corp STEREO PULSE OXIMETER
US6671526B1 (en) 1998-07-17 2003-12-30 Nihon Kohden Corporation Probe and apparatus for determining concentration of light-absorbing materials in living tissue
JP2000083933A (en) 1998-07-17 2000-03-28 Nippon Koden Corp Instrument for measuring concentration of light absorptive material in vital tissue
EP1104254A2 (en) 1998-08-13 2001-06-06 Whitland Research Limited Optical device
US6949081B1 (en) 1998-08-26 2005-09-27 Non-Invasive Technology, Inc. Sensing and interactive drug delivery
US6064898A (en) 1998-09-21 2000-05-16 Essential Medical Devices Non-invasive blood component analyzer
US6519486B1 (en) * 1998-10-15 2003-02-11 Ntc Technology Inc. Method, apparatus and system for removing motion artifacts from measurements of bodily parameters
US6438399B1 (en) 1999-02-16 2002-08-20 The Children's Hospital Of Philadelphia Multi-wavelength frequency domain near-infrared cerebral oximeter
US6675031B1 (en) 1999-04-14 2004-01-06 Mallinckrodt Inc. Method and circuit for indicating quality and accuracy of physiological measurements
US6675029B2 (en) 1999-07-22 2004-01-06 Sensys Medical, Inc. Apparatus and method for quantification of tissue hydration using diffuse reflectance spectroscopy
US7904139B2 (en) 1999-08-26 2011-03-08 Non-Invasive Technology Inc. Optical examination of biological tissue using non-contact irradiation and detection
US20050107676A1 (en) 2003-03-07 2005-05-19 Acosta George M. Method and apparatus for noninvasive glucose concentration estimation through near-infrared spectroscopy
JP2001149349A (en) 1999-11-26 2001-06-05 Nippon Koden Corp Sensor for living body
US6622095B2 (en) 1999-11-30 2003-09-16 Nihon Kohden Corporation Apparatus for determining concentrations of hemoglobins
US6415236B2 (en) 1999-11-30 2002-07-02 Nihon Kohden Corporation Apparatus for determining concentrations of hemoglobins
EP1251772B1 (en) 1999-12-22 2012-11-14 Orsense Ltd. A method of optical measurements for determining various parameters of the patient's blood
US6594513B1 (en) 2000-01-12 2003-07-15 Paul D. Jobsis Method and apparatus for determining oxygen saturation of blood in body organs
EP1274343B1 (en) 2000-04-17 2012-08-15 Nellcor Puritan Bennett LLC Pulse oximeter sensor with piece-wise function
US6430525B1 (en) * 2000-06-05 2002-08-06 Masimo Corporation Variable mode averager
IL138683A0 (en) 2000-09-25 2001-10-31 Vital Medical Ltd Apparatus and method for monitoring tissue vitality parameters
IL138884A (en) 2000-10-05 2006-07-05 Conmed Corp Pulse oximeter and a method of its operation
US6594512B2 (en) * 2000-11-21 2003-07-15 Siemens Medical Solutions Usa, Inc. Method and apparatus for estimating a physiological parameter from a physiological signal
US6501974B2 (en) 2001-01-22 2002-12-31 Datex-Ohmeda, Inc. Compensation of human variability in pulse oximetry
US6606509B2 (en) 2001-03-16 2003-08-12 Nellcor Puritan Bennett Incorporated Method and apparatus for improving the accuracy of noninvasive hematocrit measurements
US6591122B2 (en) 2001-03-16 2003-07-08 Nellcor Puritan Bennett Incorporated Device and method for monitoring body fluid and electrolyte disorders
US7239902B2 (en) 2001-03-16 2007-07-03 Nellor Puritan Bennett Incorporated Device and method for monitoring body fluid and electrolyte disorders
US7657292B2 (en) 2001-03-16 2010-02-02 Nellcor Puritan Bennett Llc Method for evaluating extracellular water concentration in tissue
SG126677A1 (en) 2001-06-26 2006-11-29 Meng Ting Choon Method and device for measuring blood sugar level
IL145445A (en) 2001-09-13 2006-12-31 Conmed Corp Signal processing method and device for signal-to-noise improvement
US6701170B2 (en) * 2001-11-02 2004-03-02 Nellcor Puritan Bennett Incorporated Blind source separation of pulse oximetry signals
CN1298284C (en) 2002-02-14 2007-02-07 加藤俊德 Apparatus for evaluating biological function
WO2003071939A1 (en) 2002-02-22 2003-09-04 Masimo Corporation Active pulse spectraphotometry
DE10213692B4 (en) 2002-03-27 2013-05-23 Weinmann Diagnostics Gmbh & Co. Kg Method for controlling a device and device for measuring ingredients in the blood
US6690958B1 (en) 2002-05-07 2004-02-10 Nostix Llc Ultrasound-guided near infrared spectrophotometer
US6711425B1 (en) 2002-05-28 2004-03-23 Ob Scientific, Inc. Pulse oximeter with calibration stabilization
US6909912B2 (en) * 2002-06-20 2005-06-21 University Of Florida Non-invasive perfusion monitor and system, specially configured oximeter probes, methods of using same, and covers for probes
US20040087836A1 (en) 2002-10-31 2004-05-06 Green Michael R. Computer system and method for closed-loop support of patient self-testing
WO2004047631A2 (en) 2002-11-22 2004-06-10 Masimo Laboratories, Inc. Blood parameter measurement system
US7006856B2 (en) 2003-01-10 2006-02-28 Nellcor Puritan Bennett Incorporated Signal quality metrics design for qualifying data for a physiological monitor
JP4284674B2 (en) 2003-01-31 2009-06-24 日本光電工業株式会社 Absorbent concentration measuring device in blood
US20050105554A1 (en) 2003-11-18 2005-05-19 Michael Kagan Method and switch system for optimizing the use of a given bandwidth in different network connections
US20050267346A1 (en) 2004-01-30 2005-12-01 3Wave Optics, Llc Non-invasive blood component measurement system
CA2555807A1 (en) 2004-02-12 2005-08-25 Biopeak Corporation Non-invasive method and apparatus for determining a physiological parameter
US7162288B2 (en) 2004-02-25 2007-01-09 Nellcor Purtain Bennett Incorporated Techniques for detecting heart pulses and reducing power consumption in sensors
US8611977B2 (en) 2004-03-08 2013-12-17 Covidien Lp Method and apparatus for optical detection of mixed venous and arterial blood pulsation in tissue
US7194293B2 (en) 2004-03-08 2007-03-20 Nellcor Puritan Bennett Incorporated Selection of ensemble averaging weights for a pulse oximeter based on signal quality metrics
US7534212B2 (en) 2004-03-08 2009-05-19 Nellcor Puritan Bennett Llc Pulse oximeter with alternate heart-rate determination
US7039538B2 (en) * 2004-03-08 2006-05-02 Nellcor Puritant Bennett Incorporated Pulse oximeter with separate ensemble averaging for oxygen saturation and heart rate
US7277741B2 (en) 2004-03-09 2007-10-02 Nellcor Puritan Bennett Incorporated Pulse oximetry motion artifact rejection using near infrared absorption by water
US7551950B2 (en) 2004-06-29 2009-06-23 O2 Medtech, Inc,. Optical apparatus and method of use for non-invasive tomographic scan of biological tissues
US7343186B2 (en) 2004-07-07 2008-03-11 Masimo Laboratories, Inc. Multi-wavelength physiological monitor
US7392075B2 (en) 2005-03-03 2008-06-24 Nellcor Puritan Bennett Incorporated Method for enhancing pulse oximetry calculations in the presence of correlated artifacts
GB2424270A (en) 2005-03-14 2006-09-20 Spectrum Medical Llp Optical monitoring of predetermined substances in blood
US8092379B2 (en) 2005-09-29 2012-01-10 Nellcor Puritan Bennett Llc Method and system for determining when to reposition a physiological sensor
US7869850B2 (en) 2005-09-29 2011-01-11 Nellcor Puritan Bennett Llc Medical sensor for reducing motion artifacts and technique for using the same
US20070077200A1 (en) 2005-09-30 2007-04-05 Baker Clark R Method and system for controlled maintenance of hypoxia for therapeutic or diagnostic purposes
US20070100220A1 (en) 2005-10-28 2007-05-03 Baker Clark R Jr Adjusting parameters used in pulse oximetry analysis
US7706852B2 (en) 2006-01-30 2010-04-27 Nellcor Puritan Bennett Llc System and method for detection of unstable oxygen saturation
US7848891B2 (en) * 2006-09-29 2010-12-07 Nellcor Puritan Bennett Llc Modulation ratio determination with accommodation of uncertainty

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5339818A (en) * 1989-09-20 1994-08-23 University Of Utah Research Foundation Method for determining blood pressure utilizing a neural network
US6371114B1 (en) * 1998-07-24 2002-04-16 Minnesota Innovative Technologies & Instruments Corporation Control device for supplying supplemental respiratory oxygen
US20050049470A1 (en) * 2003-08-27 2005-03-03 Terry Alvin Mark Multi-domain motion estimation and plethysmographic recognition using fuzzy neural-nets

Also Published As

Publication number Publication date
US8862194B2 (en) 2014-10-14
US20090326348A1 (en) 2009-12-31

Similar Documents

Publication Publication Date Title
CA2283858C (en) Method and apparatus for adaptively averaging data signals
US9220440B2 (en) Determining a characteristic respiration rate
EP2408354B1 (en) Systems and methods for monitoring heart rate and blood pressure correlation
CA2728437C (en) Consistent signal selection by signal segment selection techniques
US20170143272A1 (en) Rejecting Noise in a Signal
US9198582B2 (en) Determining a characteristic physiological parameter
CA2735121C (en) Detecting a signal quality decrease in a measurement system
Yang et al. Estimation and validation of arterial blood pressure using photoplethysmogram morphology features in conjunction with pulse arrival time in large open databases
WO2010128275A1 (en) Selection of signal regions for parameter extraction
JP2001517128A (en) Method and apparatus for arbitrating and removing harmonics to obtain the best evaluation of blood component values
WO2012072849A2 (en) System and method for the simultaneous, non-invasive estimation of blood glucose, glucocorticoid level and blood pressure
WO1999025241A1 (en) Determining the state of an oximetry sensor
US20110098933A1 (en) Systems And Methods For Processing Oximetry Signals Using Least Median Squares Techniques
WO2014100084A1 (en) Systems and methods for distinguishing between central apnea and obstructive apnea
US20140371554A1 (en) Method for improved oxygen saturation estimation in the presence of noise
US10898087B2 (en) Motion detection and cancellation using ambient light
US20230148961A1 (en) Systems and methods for computationally efficient non-invasive blood quality measurement
CA2283856C (en) Method and apparatus for harmonically filtering data
WO2017100519A1 (en) Systems and methods for adaptable presentation of sensor data
CN110811582A (en) Heart rate detection method, device and equipment
JP4571220B2 (en) Data signal adaptive averaging method and apparatus
WO2020012807A1 (en) Biological information processing device and biological information processing method

Legal Events

Date Code Title Description
AS Assignment

Owner name: COVIDIEN LP, MASSACHUSETTS

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:BAKER, CLARK R., JR;REEL/FRAME:033805/0266

Effective date: 20080507

STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION