`US 6,393,311 B1
`(10) Patent No.:
`May21, 2002
`(45) Date of Patent:
`Edgar, Jr. et al.
`
`US006393311B1
`
`(54)
`
`(75)
`
`METHOD, APPARATUS AND SYSTEM FOR
`REMOVING MOTION ARTIFACTS FROM
`MEASUREMENTS OF BODILY
`PARAMETERS
`
`Inventors: Reuben W. Edgar, Jr.; August J. Allo,
`Jr., both of San Antonio, TX (US);
`Jesus D. Martin, Wallingford, CT
`(US); John R. DelFavero, East
`Hampton, CT (US); Michael B. Jaffe,
`Cheshire, CT (US)
`
`(73)
`
`Assignee:
`
`NTC Technology Inc., Wilmington, DE
`(US)
`
`Notice:
`
`Subject to any disclaimer, the term of this
`patent is extended or adjusted under 35
`US.C. 154(b) by 0 days.
`
`(21)
`
`(22)
`
`(60)
`
`(61)
`(52)
`
`(58)
`
`(56)
`
`Appl. No.:
`Filed:
`
`09/410,991
`
`Oct. 1, 1999
`
`Related U.S. Application Data
`Provisional application No. 60/104,422, filed on Oct. 15,
`1998.
`
`Tint, C07 ieee cscs ces ce eeesteseeeseeseeeess A61B 5/00
`U.S. Ch. veces 600/323; 600/336; 600/310;
`600/324.
`Field of Search ................ 600/309-311, 322-326,
`600/330-331, 336, 473, 476; 356/39-41;
`369/60.01
`
`References Cited
`US. PATENT DOCUMENTS
`
`4,800,495 A *
`4,942,877 A *
`4,955,379 A *
`5,025,791 A
`eerie A
`;266,
`5,299,120 A *
`5,349,952 A *
`5,351,685 A
`
`1/1989 Smith wc 600/322
`
`.. 600/323
`7/1990 Sakaietal.
`se eeeesceseecnsceescescees 600/366
`9/1990 Hall
`6/1991 Niwa
`5hoo, en et a \
`wedlow et al.
`3/1994 Kaestle veccceccscecssseeeseees 600/310
`9/1994 McCarthy et al.
`.......... 600/473
`10/1994 Potratz
`
`5,368,026 A
`5,368,224 A
`5,398,680 A
`
`11/1994 Swedlowetal.
`11/1994 Richardsonetal.
`3/1995 Polsonetal.
`
`(List continued on next page.)
`
`OTHER PUBLICATIONS
`
`Dowla, et al., Neural Networks and Wavelet Analysis in the
`Computer Interpretation of Pulse Oximetry Data, Neural
`Networks for Signal Processing VI—Proc.
`IEEE, 1996
`IEEE
`Signal
`Process,
`Soc.
`TEEE Workshop,
`0-7803-3550-3 (1996).
`
`Primary Examiner—Enic F. Winakur
`Assistant Examiner—Matthew Kremer
`(74) Attorney, Agent, or Firm—TraskBritt
`
`(57)
`
`ABSTRACT
`
`A method for removing motion artifacts from devices for
`sensing bodily parameters and apparatus and system for
`effecting same. The method includes analyzing segments of
`measured data representing bodily parameters and possibly
`noise from motionartifacts. Each segment of measured data
`may correspond to a single light signal
`transmitted and
`detected after transmission or reflection through bodily
`tissue. Each data segment is frequency analyzed to deter-
`mine dominant frequency components. The frequency com-
`ponent which represents at least one bodily parameter of
`interest is selected for further processing. The segment of
`data is subdivided into subsegments, each subsegment rep-
`resenting one heartbeat. The subsegments are used to cal-
`culate a modified average pulse as a candidate output pulse.
`The candidate output pulse is analyzed to determine whether
`it is a valid bodily parameter and,if yes, it is output for use
`in calculating the at least one bodily parameter of interest
`without any substantial noise degradation. The above
`method may be applied to red and infrared pulse oximetry
`signals prior to calculating pulsatile blood oxygen concen-
`tration. Apparatus and systems disclosed incorporate meth-
`.
`.
`.
`.
`dS disclosed according to the invention.
`
`34 Claims, 11 Drawing Sheets
`
`Start
`
`vy
`of raw data
`100 | nequire a segment
`
`
`y
`
`ToAesa
`
`
`;
`y
`120° “a
`
`¥v
`130 ~™ ‘Ouiput average
`pulse signal
`
`y
`
`
`
`
`
`
`140ene>
`NO
`
`End
`
`1
`
`APPLE 1005
`
`APPLE 1005
`
`1
`
`
`
`US 6,393,311 B1
`
`Page 2
`
`9/1998
`9/1998
`10/1998
`* 12/1998
`12/1998
`3/1999
`7/1999
`8/1999
`5/2000
`9/2000
`9/2000
`
`*
`*
`*
`
`Kaestle
`Potratz
`Polson etal.
`Chen et al. wee 375/344
`Baker, Jr. et al.
`Richardsonetal.
`Diab
`Mortz
`Diab et al. vo. eee 600/310
`Hermanskyet al.
`........ 704/226
`Kaestle et al. we. 600/322
`
`AAAAAAAAAAA
`
`5,800,348
`5,803,910
`5,820,550
`5,852,638
`5,853,364
`5,885,213
`5,919,134
`5,934,277
`6,067,462
`6,098,038
`6,122,535
`
`* cited by examiner
`
`5,431,170
`5,448,991
`5,482,036
`5,490,505
`5,555,882
`5,588,427
`5,632,272
`5,645,060
`5,685,299
`5,713,355
`5,743,263
`5,769,785
`
`PPPrrrrrrrrpy
`
`U.S. PATENT DOCUMENTS
`
`7/1995
`9/1995
`1/1996
`2/1996
`9/1996
`12/1996
`5/1997
`7/1997
`11/1997
`2/1998
`4/1998
`6/1998
`
`Matthews
`Polsonetal.
`Diab etal.
`Diab etal.
`Richardsonetal.
`Tien
`Diab etal.
`Yorkey
`Diab etal.
`Richardsonetal.
`Baker, Jr.
`Diab etal.
`
`2
`
`
`
`U.S. Patent
`
`May 21, 2002
`
`Sheet 1 of 11
`
`US 6,393,311 B1
`
`100
`
`110
`
`120
`
`130
`
`Acquire a segment
`of raw data
`
`Analyze frequency
`componenis
`
`Output average
`pulse signal
`
`
`
`
`
`
`
`
`
`Determine best
`frequency
`component
`
`New segment
`140 oN of raw data?
`
`FIG.
`
`|
`
`3
`
`
`
`U.S. Patent
`
`May 21, 2002
`
`Sheet 2 of 11
`
`US 6,393,311 B1
`
`
`
`€>
`
`
`
`
`
`RawInfrared Data Segment
`@ 3.60E+05
`3.50E+05
`QO23.40E+05
`@ 3.30E+05
`a
`
`(
`
`@ 3.20E+05
`Cc
`
`180 200
`16.0
`140
`Time (Seconds)
`
`3.10E+Oo 9 120
`
`+
`
`>“
`
`4
`
`
`
`U.S. Patent
`
`May 21, 2002
`
`Sheet 3 of 11
`
`US 6,393,311 B1
`
`
`
`Power Spectrum of Data from FIG. 1
`50000
`
`
`
`~
`
`40000
`
`
`
`
`
`
`5 30000
`
`5a
`
`. 20000
`
`10000
`
`0
`
`0.5
`
`2.5
`1.5
`Frequency (Hertz)
`
`3.5
`
`5
`
`
`
`U.S. Patent
`
`May 21, 2002
`
`Sheet 4 of 11
`
`US 6,393,311 B1
`
`
`
`FIG. 4
`BandpassFiltered and Raw IR Data
`2 3.6E+05
`
`
`
`TUN)
`
`
`
`8evo ae iy
`B32660 AH
`is
`SAT NING
`AWN
`= 3,2E+05 +
`531E+05
`D
`41.0
`
`!
`
`19.0
`
`13.0
`
`17.0
`15.0
`Time (Seconds)
`
`—RawIRData — Filtered IR Data
`
`
`6
`
`
`
`U.S. Patent
`
`May 21, 2002
`
`Sheet 5 of 11
`
`US 6,393,311 B1
`
`12 Superimposed Pulse Subsegments
`1.5E+04
`
` FIG. 5
`
`
`iLg
`po AWEWBA
`Beoraoa JNpeaY!Gut
`5 -1.0E+04
`Us
`
`
` é 1.0E+04 LN,a»
`
`Q
`
`-1.5E+04
`
`10 20 30 40 50 60 70 80 90 100
`Data Points
`
`7
`
`
`
`U.S. Patent
`
`May 21, 2002
`
`Sheet 6 of 11
`
`US 6,393,311 B1
`
`
`
`FIG. 6
`Modified Average Pulse
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`0 010203040506 07 08 0.9
`Time (Seconds)
`Po
`
`,
`
`_ 6000 +
`co)
`® 4000
`—!
`
`®@ 2000 -
`
`0
`
`-2000
`
`-4000
`6000
`
`Da
`
`o©
`
`c©5
`
`8
`
`
`
`
`
`FIG. 7
`
`9
`
`
`
`U.S. Patent
`
`May 21, 2002
`
`Sheet 8 of 11
`
`US 6,393,311 B1
`
`
`
`
`ie)
`
`
`Processor
`Device
`
`
`
`
` Output
`Input
`
`Device
`
`
`Device
`
`Motion Artifact Circuit Card
`
`10
`
`
`
`U.S. Patent
`
`May 21, 2002
`
`Sheet 9 of 11
`
`US 6,393,311 B1
`
`Input
`Device
`
`Device
`
`Processor
`
`storage
`
`11
`
`
`
`U.S. Patent
`
`May 21, 2002
`
`Sheet 10 of 11
`
`US 6,393,311 B1
`
`
`
`200 Segmentwith FIR(B)
`
`Bandpass Filter Data
`
`Filter
`
`210
`
`220
`
`230
`
`240
`
`Subsegment Data
`
`“| 120
`
`on Newest Pulse
`
`Measure Dispersion
`
`of Saturation Based
`
`
` YES
`
`Is the
`
`
`Dispersion Low?
`
`
`NO
`
`Pulse
`
`
`
`
`
`Use Newest Pulse
`Instead of
`Computing Average
`
`Compute Average
`Pulse
`
`250
`
`FIG. 11
`
`12
`
`
`
`U.S. Patent
`
`May 21, 2002
`
`Sheet 11 of 11
`
`US 6,393,311 B1
`
`Lia
`
`260
`
`2/0
`
`280
`
`290
`
`
`
`
`Evaluate Pulse
`
`Shape Against Rules
`to Obtain Quality
`
`Measure
`
`
`
`
`
`
`
`
`Minimize the
`
`
`Dispersion of
`Saturation Values
`
`
`
`
`
`
`Last Candidate
`Frequency?
`
`
`
`
`
`
`Select Candidate
`
`
`Frequencywith Best
`Quality Measure
`
`FIG. 12
`
`13
`
`
`
`US 6,393,311 B1
`
`1
`
`METHOD, APPARATUS AND SYSTEM FOR
`REMOVING MOTION ARTIFACTS FROM
`MEASUREMENTSOF BODILY
`PARAMETERS
`
`CROSS-REFERENCE TO RELATED
`APPLICATION
`
`This utility patent application claims the benefit of U.S.
`provisional patent application Ser. No. 60/104,422, filed
`Oct. 15, 1998.
`
`BACKGROUND OF THE INVENTION
`
`1. Technical Field
`
`This invention relates to the field of signal processing.
`Moreparticularly, this invention relates to processing mea-
`sured signals to remove unwanted signal components caused
`by noise and especially noise caused by motionartifacts.
`2. State of the Art
`
`The measurement of physiological signals can often be
`difficult because the underlying physiological processes may
`generate very low level signals. Furthermore, interfering
`noise is inherent in the body and the interface between the
`body and sensors of physiological processes, Examples of
`physiological measurements include: measurementof elec-
`trocardiogram (ECG) signals based on the electrical depo-
`larization of the heart muscle, blood pressure, blood oxygen
`saturation, partial pressure of CO,, heart rate, respiration
`rate, and depth of anesthesia. ECG signals, for example, are
`typically detected by surface electrodes mounted on the
`chest of a patient. ECG signals are weakat the signal source
`(i.e., the heart) and are even weaker at the surface of the
`chest. Furthermore, electrical interference from the activity
`of other muscles (e.g., noise caused by patient breathing,
`general movement, etc.) causes additional interference with
`physiological signals such as an ECG. Thus, considerable
`care must be taken in the design and use of physiological
`processors to enhance the quality of the true signal and
`reduce the effects of interfering noise signals.
`It is convenient to characterize a measured signal as being
`a composite signal composedof a true signal component and
`a noise signal component. The terms “measured signal” and
`“composite signal” will be used interchangeably hereinafter.
`Signal processors are frequently used to removenoise signal
`components from a composite measured signal in order to
`obtain a signal which closely, if not identically, represents
`the true signal. Conventional filtering techniques such as
`low pass, band pass, and high passfiltering can be used to
`remove noise signal components from the measured com-
`posite signal where the noise signal component occupies a
`frequency range outside the true signal component. More
`sophisticated techniques for conventional noise filtering
`include multiple notch filters, which are suitable for use
`where the noise signal componentexists at multiple, distinct
`frequencies, all outside the true signal frequency band.
`However, it is often the case that the frequency spectrum
`of the true and noise signal components overlap and that the
`statistical properties of both signal components change with
`time. More importantly, there are many cases wherelittle is
`known about the noise signal component. In such cases,
`conventional filtering techniques are ineffective in extract-
`ing the true signal.
`The measurement of oxygen saturation in the blood of a
`patient is a common physiological measurement the accu-
`racy of which may be compromised by the presence of
`noise. Knowledge of blood oxygen saturation can becritical
`
`10
`
`15
`
`20
`
`25
`
`30
`
`40
`
`45
`
`50
`
`55
`
`60
`
`65
`
`2
`during surgery. There are means of obtaining blood oxygen
`saturation by invasive techniques, such as extracting and
`testing blood removed from a patient using a co-oximeter.
`But, such invasive means are typically time consuming,
`expensive, and uncomfortable for the patient. Fortunately,
`non-invasive measurements of blood oxygen saturation can
`be made using known properties of energy attenuation as a
`selected form of energy passes through a bodily medium.
`Such non-invasive measurements are performed routinely
`with a pulse oximeter.
`The basic idea behind such energy attenuation measure-
`ments is as follows. Radiant energy is directed toward a
`bodily medium, where the medium is derived from or
`contained within a patient, and the amplitude of the energy
`transmitted through or reflected from the medium is then
`measured. The amountof attenuation of the incident energy
`caused by the medium is strongly dependent on the thick-
`ness and composition of the medium through which the
`energy must pass, as well as the specific form of energy
`selected. Information about a physiological system can be
`derived from data taken from the attenuated signal of the
`incident energy transmitted or reflected. However, the accu-
`racy of such information is reduced where the measured
`signal includes noise. Furthermore, non-invasive measure-
`ments often do not afford the opportunity to selectively
`observe the interference causing the noise signal component,
`making it difficult to remove.
`A pulse oximeter is one example of a physiological
`monitoring system which is based upon the measurementof
`energy attenuated by biological
`tissues and substances.
`Morespecifically, a pulse oximeter measures the variable
`absorption caused by arterial blood volume changes. Pulse
`oximeters transmit electromagnetic energy at two different
`wavelengths, typically at 660 nm (red) and 940 nm (infrared,
`hereinafter IR) into the tissue and measure the attenuation of
`the energy as a function of time. The output signal of a pulse
`oximeter is sensitive to the pulsatile portion of the arterial
`blood flow and contains a component which is a waveform
`representative of the patient’s arterial pulse. This type of
`signal, which contains a componentrelated to the patient’s
`pulse, is called a plethysmographic waveform or plethys-
`mogram.
`Pulse oximetry measurements typically use a digit, such
`as a finger, or an ear lobe or other elementof the body, where
`blood flows close to the skin as the medium through which
`light energy is transmitted. The finger, for example,
`is
`composed of various tissues and substances including skin,
`fat, bone, muscle, blood, etc. The extent to which each of
`these biological tissues and substances attenuate incident
`electromagnetic energy is generally known. However, the
`effect of motion can cause changesin the optical coupling of
`the sensor (or probe)
`to the finger,
`the underlying
`physiology, the local vasculature, optical properties of tis-
`sues due to changing optical path length as well as combi-
`nations andinteractions of the all of the above. Thus, patient
`motion may cause erratic energy attenuation.
`A typical pulse oximeter includes a sensor, cabling from
`the sensor to a computer for signal processing and visual
`display,
`the computer and visual display typically being
`included in a patient monitor. The sensor typically includes
`two light emitting diodes (LEDs) placed across a finger tip
`and a photodetector on the side opposite the LEDs. Each
`LED emits a light signal at different frequencies. The
`detector measures both transmitted light signals once they
`have passed through the finger. The signals are routed to a
`computer for analysis and. display of the various parameters
`measured.
`
`14
`
`14
`
`
`
`US 6,393,311 B1
`
`3
`The underlying physical basis of a pulse oximeter is
`Beer’s law (also referred to as Beer-Lambert’s or Bouguer’s
`law) which described attenuation of monochromatic light
`traveling through a uniform medium which absorbs light
`with the equation:
`
`Transmitted=!incident®
`
`—doa()
`
`@)
`
`where I,,,nsmitred IS the intensity of the light transmitted
`through the uniform medium,I,,,.;aen, 1S the intensity of
`incident light, d is the distance light is transmitted through
`the uniform medium,c is the concentration of the absorbing
`substance in the uniform medium, expressed in units of
`mmol L-1, and a() is the extinction or absorption coeffi-
`cient of the absorbing substance at wavelength i, expressed
`in units of L(mmol cm). The properties of Beer’s law are
`valid even if more than one substance absorbs light in the
`medium. Each light absorbing substance contributesits part
`to the total absorbance.
`
`Each LED emits light at different wavelengths, typically
`red (centered at about 660 nm) and IR (centered at about 940
`nm) frequency bands. The intensity of light
`transmitted
`throughtissue, I,nsmirea1S different for each wavelength of
`light emitted by the LEDs. Oxyhemoglobin (oxygenated
`blood) tends to absorb IR light, whereas deoxyhemoglobin
`(deoxygenated blood) tends to absorb red light. Thus, the
`absorption of IR lightrelative to the red light increases with
`oxyhemoglobin. The ratio of the absorption coefficients can
`be used to determine the oxygen saturation of the blood.
`To estimate blood oxygen saturation, SpO., a two-solute
`concentration is assumed. A measure of functional blood
`
`oxygen saturation level, SpO., can be defined as:
`
`Co
`SpO, =
`pe Cr +Co
`
`;
`
`(2)
`
`where c, represents oxyhemoglobin solute concentration,
`and C, represents reduced or deoxyhemoglobin solute con-
`centration.
`Noise signal components in a measured pulse oximetry
`light signal can originate from both AC and DC sources. DC
`noise signal components may be caused by transmission of
`electromagnetic energy through tissues of relatively con-
`stant thickness within the body, e.g., bone, muscle, skin,
`blood, etc. Such DC noise signal components maybe easily
`removed with conventional filtering techniques. AC noise
`signal components may occur whentissues being measured
`are perturbed and, thus, change in thickness while a mea-
`surement is being made. Such AC noise signal components
`are difficult
`to remove with conventional filtering tech-
`niques. Since most materials in and derived from the body
`are easily compressed, the thickness of such matter changes
`if the patient moves during a non-invasive physiological
`measurement. Thus, patient movement can cause the prop-
`erties of energy attenuation to vary erratically. The erratic or
`unpredictable nature of motion artifacts induced by noise
`signal components is a major obstacle in removing them.
`Various approaches to removing motion artifacts from
`measured physiological signals, and particularly for use in
`pulse oximeters, have been proposed. U.S. Pat. Nos. 5,482,
`036, 5,490,505, 5,632,272, 5,685,299, and 5,769,785, all to
`Diab et al., and U.S. Pat. No. 5,919,134 to Diab, disclose
`methods and apparatuses for removing motion artifacts
`using adaptive noise cancellation techniques. The basic
`proposition behind the Diab et al. approach is to first
`generate a noise reference signal from the two measured
`signals, and then use the noise reference signal as an input
`
`10
`
`15
`
`20
`
`25
`
`30
`
`35
`
`40
`
`45
`
`50
`
`55
`
`60
`
`65
`
`4
`to an adaptive noise canceler along with either or both of the
`measured signals to remove the reference noise signal from
`the measured signals, thus approximating the actual para-
`metric signals of interest. The Diab et al. approach appears
`to require the use of both measured input signals to generate
`a noise reference signal.
`Another approachto noise artifact elimination is disclosed
`in USS. Pat. No. 5,588,427 to Tien. Tien uses fractal dimen-
`sion analysis to determine the complexity of waveforms in
`order to determine the proper value of the ratio of true
`intensities based on signal complexity. The Tien approach
`employs a fractal analyzer to determine values for two
`ratios, « and B, based on the measured time varying intensity
`of the transmitted red and IR light signals including noise.
`a. is defined as the ratio of the time varying true intensity of
`light transmitted from the red LED andthe time varying true
`intensity of the light transmitted from the IR LED. n is a
`similar ratio relating the noise introduced during the mea-
`surement of the light transmitted by the red LED and the
`noise introduced during the measurementof the light trans-
`mitted by the IR LED. According to Tien, a fractal analyzer
`then determines values for a and B and provides (a8)pairs
`to a statistical analyzer. The statistical analyzer performs
`analysis of one or more (0,8) pairs to determine the best
`value for a, which is then provided to a look-up table. The
`look-up table provides a value corresponding to the arterial
`oxygen saturation in the patient. While the Tien approach
`appears to be an innovative use of fractal analysis, it also
`appears to be computationally complex.
`Yet another approach to noise artifact elimination is
`disclosed in U.S. Pat. Nos. 5,885,213, 5,713,355, 5,555,882
`and 5,368,224, all to Richardsonet al. The basic proposition
`behind the Richardsonet al. approach is to switch operative
`frequencies periodically based on evaluating the noise level
`associated with various possible frequencies of operation in
`order to select the frequency of operation that has the lowest
`associated noise level. It would appear that data measured at
`a noisy frequency, using the Richardsonet al. approach may
`be invalid or useless for calculating arterial oxygen satura-
`tion. Furthermore, Richardson et al. requires a computa-
`tional overhead to constantly monitor which frequency of
`operation provides the least noise.
`Another approachto noise artifact elimination is disclosed
`in U.S. Pat. No. 5,853,364 to Baker, Jr et al. The Baker, Jr.
`et al. approach first calculates the heart rate of the patient
`using an adaptive comb filter, power spectrum and pattern
`matching. Once the heart rate is determined, the oximetty
`data is adaptively combfiltered so that only energyat integer
`multiples of the heart rate are processed. The combfiltered
`data and the raw oximetry data are filtered using a K nan
`filter to adaptively modify averaging weights and averaging
`times to attenuate motion artifact noise. The adaptivefilter-
`ing of the Baker,Jr. et al. approach appears to add significant
`computational complexity to solve the problem of motion
`artifact rejection.
`Still another approach to noise artifact elimination is
`disclosed in U.S. Pat. No. 5,431,170 to Mathews. Mathews
`couples a conventional pulse oximeter light transmitter and
`receiver with a transducer responsive to movementor vibra-
`tion of the body. The transducer provides an electrical signal
`varying according to the body movements or vibrations,
`which is relatively independent of the blood or other fluid
`flow pulsations. Mathews then provides means for compar-
`ing the light signals measured with the transducer output and
`performing adaptive noise cancellation. An apparent disad-
`vantage of the Mathews approachis the need for a secondary
`sensor to detect motion.
`
`15
`
`15
`
`
`
`5
`Thus, a need in the art exists for a method, apparatus and
`system to eliminate motion-induced noise artifacts from
`light signals, that is relatively simple computationally, and
`that does not require more than one sensor.
`SUMMARYOF THE INVENTION
`
`The present invention includes methods, apparatuses and
`systems for removing noise in physiological measurements
`caused by motion or other similar artifacts. The methods,
`apparatuses and systems of the present invention eliminate
`noise from light signals using a single conventional sensor
`and are relatively simple computationally.
`In a method embodiment, a segment of data from a
`measured pulse oximetry signal is conventionally filtered,
`and frequency analyzed for major frequency components.
`The frequency components with the largest power spectral
`density are selected for subdividing into subsegments, each
`comprising an individual heartbeat. The subsegments are
`averaged and then analyzed to determine if the averaged
`subsegmentis a valid pulse oximetry signal. Additionally,
`various quality or confidence measures may be used to
`evaluate the validity of such signal. Valid averaged subseg-
`ments become outputs for further processing to calculate
`physiological parameters such as blood oxygen saturation
`levels.
`
`includes a processor and
`A circuit card embodiment
`memory for storing a computer program capable of execut-
`ing instructions embodying the above method.
`Asystem embodimentincludes an input device, an output
`device, a memory device and a motion artifact rejection
`circuit card capable of executing instructions stored in the
`memory device implementing the methods described herein.
`Finally, a system embodiment includes an input device,
`and output device, a memory device and a processor, which
`may be a digital signal processor, capable of executing
`instructions stored in the memory device implementing the
`methods described herein.
`
`BRIEF DESCRIPTION OF THE DRAWING
`FIGURES
`
`is currently
`In the drawings, which illustrate what
`regarded as the best mode for carrying out the invention and
`in which like reference numerals refer to like parts in
`different views or embodiments:
`
`15
`
`20
`
`25
`
`30
`
`35
`
`40
`
`FIG. 1 is a high-level flowchart of a method embodiment
`of the invention.
`
`45
`
`FIG. 2 is a graph of a segment of data representing
`transmitted infrared (IR) light data signal received by a pulse
`oximetry sensor suitable for signal processing in accordance
`with the invention.
`
`FIG. 3 is a graph of the power spectrum of the IR data
`segment in FIG. 2 in accordance with the invention.
`FIG. 4 illustrates a bandpass filtered IR data segment
`superimposed on the original IR data segment of FIG. 2 in
`accordance with the invention.
`
`FIG. 5 is a graph of each pulse subsegment detrended and
`superimposed upon one another in accordance with the
`second method of computing an average pulse.
`FIG. 6 shows the modified average pulse as computed
`from the detrended pulse subsegments shown in FIG. 5.
`FIG. 7 illustrates the acquisition of a new pulse segment
`of raw IR data and the removalof the oldest segment of raw
`IR data in accordance with the invention.
`
`FIG. 8 is a block diagram of a motion artifact rejection
`circuit card configured to removenoise artifacts from signals
`representing bodily parameters in accordance with the
`invention.
`
`50
`
`55
`
`60
`
`65
`
`US 6,393,311 B1
`
`6
`FIG. 9 is a block diagram of a pulse oximetry system
`including a motion artifact rejection circuit card capable of
`removing noise from pulse oximetry data in accordance with
`the invention.
`
`FIG. 10 is a block diagram of a pulse oximetry system
`including a processor device programmed to remove noise
`from pulse oximetry data in accordance with the invention.
`FIGS. 11 and 12 are a detailed flowchart of the method of
`
`10
`
`determining the best frequency component, block 120, of
`FIG. 1.
`
`DETAILED DESCRIPTION OF THE
`INVENTION
`
`The following detailed description discloses methods,
`apparatuses and systems for removing motionartifacts from
`measured plethysmographic waveforms, particularly, but
`without limitation, those used in pulse oximetry. A system
`embodimentof the invention includes pulse oximetry hard-
`ware and associated software to perform the motion artifact
`suppression. A method embodiment of this invention
`includes a series of steps which exploits certain character-
`istics of plethysmographic waveforms. The methods, appa-
`ratuses and systems described below are suitable for light
`transmitted or reflected through bodily tissues and sub-
`stances. For convenience, the following detailed description
`will assume measurementof light which has been transmit-
`ted through a finger of a human. The terms “signal” and
`“waveform” are used interchangeably herein.
`FIG. 1 is a high-level flowchart of a method embodiment
`of the invention. The method steps include acquiring a
`segment of raw data 100, either red or IR, analyzing the data
`segment for dominant frequency components 110, determin-
`ing the frequency component which represents a valid
`plethysmographic pulse 120, computing an average pulse
`based on the correct frequency component 130 and repeating
`for new raw data segments 140. In order to calculate blood
`oxygen concentration, SpO., the method embodimentof the
`invention may beapplied to both red and IR data signals to
`eliminate or reduce noise from the data signals prior to
`calculating SpO,.
`The method of this invention begins with acquiring a
`segment of data (e.g., five or more pulses or approximately
`ten seconds) measured from a single light source transmitted
`through a finger and detected with a sensor on the opposite
`side of the finger. Acquiring a data segmentis illustrated by
`block 100 of FIG. 1. For convenience, a 10.24 second
`segment of data will be used to illustrate the method. This
`corresponds to 1024 data points with a sampling rate of 100
`data points per second. It should be readily apparent that the
`method of the invention is not limited to data segments of
`this size. FIG. 2 is an example of such a data segmentfor an
`IR light source. The signal processing steps described herein
`may be performed on both red and IR data segments
`independently and simultaneously. Thus, while the steps of
`the method may be illustrated with data from an IR light
`signal, the same steps are applicable to data from a red light
`signal and vice versa. The terms “data segment’, “input
`waveform”, “data signal” and “signal” are used interchange-
`ably herein.
`A segment of data may be received from a sensor that
`converts transmitted or reflected light signals into electrical
`signals. Once a segment of data from a single electrical
`signal has been acquired, it may befiltered to reduce spectral
`leakage resulting from frequency analysis. There are several
`windowfilters which may be suitable for such purposes. For
`example, and not by way of limitation, a Hanning window
`16
`
`16
`
`
`
`US 6,393,311 B1
`
`10
`
`20
`
`8
`7
`heartbeat pulse within the data segment in order to subseg-
`may be used to reduce spectral leakage. It will be readily
`apparentto one of skill in the art of digital signal processing
`ment the data into its consecutive heartbeat pulses.
`that other window filters and methods of filtering data to
`The period of rhythmic contraction of the heart by which
`reduce spectral leakage may be selected for reducing spec-
`blood is driven through the aorta and pulmonary artery is
`tral leakage. As methods of filtermg and variousfilters are
`knownas systole. Maximum light absorbance occurs during
`knownto oneof skill in the art of signal processing, they will
`the systole of a cardiac cycle and is indicated on a plethys-
`not be further detailed herein. The filtered data is then
`mogram bya low point or systolic valley. Conversely, the
`frequency analyzed to determine the dominant frequency
`period of rhythmic relaxation and dilation of the heart
`components, see block 110 of FIG. 1. FIG. 3 illustrates the
`cavities occurs during diastole when blood is drawn into the
`power spectrum of the IR data segment of FIG. 2 after
`heart cavities. Minimum light absorbance occurs during the
`filtering.
`diastole of a cardiac cycle and is indicated on a plethysmo-
`Signal processing as described herein is generally per-
`gram by a high point or diastolic peak. While it is theoreti-
`formed in the frequency domain. The segment of data is
`cally possible to useafirst derivative of the data segmentto
`converted into the frequency domain by, for example, per-
`identify transitions between diastole and systole and vice
`15
`forming the well-known Fast Fourier Transform (FFT) on
`versa, in practice, the data may be so noisy that the first
`the data. Other common techniques of converting time-
`derivative provides useless information.
`domain data to the frequency domain mayalso be used,e.g.,
`Subsegmenting each pulse within a data segment accord-
`classical methods using the FFT such as the periodogram or
`ing to the invention begins with a narrow bandpassfilter at
`correlogram, autoregressive methods, Prony’s method,
`(or near) F, with a frequency spread (bandpass window
`minimum variance methods, maximum likelihood methods.
`width) of approximately +0.25 Hertz about the center fre-
`Additionally,
`time domain data may be converted to the
`quency of F,, see block 200 of FIG. 11. To improve
`frequency domain using transforms such as discrete cosine
`discrimination, especially with closely spaced peaks, the
`transform, wavelet transform, discrete Hartley transform,
`and Gabortransform.
`bandpassfilter coefficients may be generated and adjusted as
`neededso that the center frequency is nearly identical to the
`candidate frequency, F,. The resulting bandpassfiltered data
`will resemble a sinusoidal waveform andis used to identify
`the point of diastolic peak for each pulse. Diastolic peaks in
`the data segment will occurat or near peaks in the sinusoidal
`waveform of the bandpassfiltered data. FIG. 4 illustrates the
`bandpassfiltered data superimposed on the original IR data
`segment of FIG. 2.
`Each pulse subsegmentis defined as starting one quarter
`pulse width before a diastolic peak and ending approxi-
`mately one quarter pulse width after the next diastolic peak,
`see block 210 of FIG. 11. This ensures that each heartbeat
`
`Both transient and periodic noise artifacts can induce
`peaks in the frequency domain that may be larger than the
`peak caused by the patient’s heart rate. The frequency peak
`which actually represents the patient’s heart rate must then
`be determined, see block 120 of FIG. 1. One approach to
`determining the correct frequencyis to order the frequencies
`by peak amplitude from largest to smallest, F, to F,,, where
`F, through F,, are not harmonics of each other, and analyze
`them one by one to find the correct frequency,
`i.e., the
`patient’s heart rate. For purposes ofillustration, only the
`frequencies associated with the two largest power spectrum
`amplitude (peaks), F, and F,, will be used to explain the
`signal processing in accordance with the invention. It will be
`readily apparent that the signal processing described herein
`may be extended from 2 to n candidate frequencies. For
`convenience of notation, F, is the frequency of the largest
`amplitude peak, and F, is the next largest peak, which is not
`a harmonic of F,.
`Wherethis is not the first analysis cycle, an additional
`check is made to determine if one of the two potential
`frequencies, F, and F,,
`is similar to a known valid
`frequency, F,, of the patient’s heart rate as determined
`during the previous analysis cycle. Otherwise, the signal
`processing proceeds as described in the next paragraph.
`Historical trends in heart rate may be used to select the
`proper frequency peak.
`In the