`
`(19) United States
`(12) Patent Application Publication (10) Pub. No.: US 2009/0105556 A1
`Fricke et al.
`(43) Pub. Date: Apr. 23, 2009
`
`
`(54) MEASUREMENT OF PHYSIOLOGICAL
`SIGNALS
`
`(75)
`
`Inventors:
`
`John Robert Fricke, Lexington,
`MA (US); Matthew Corbin
`Wiggins, Concord, MA (US)
`
`Correspondence Address:
`OCCHIUTI ROHLICEK & TSAO, LLP
`10 FAWCETT STREET
`
`CAMBRIDGE, MA 02138 (US)
`
`(73) Assignee:
`
`Tiax LLC, Cambridge, MA (US)
`
`(21) App]. No.:
`
`12/240,651
`
`(22)
`
`Filed:
`
`Sep. 29, 2008
`
`Related US. Application Data
`
`(60)
`
`Provisional application No. 60/995,723, filed on Sep.
`28, 2007.
`
`Publication Classification
`
`(51)
`
`Int. Cl.
`(2006.01)
`A613 5/00
`(2006.01)
`A613 5/1455
`(52) U.S.Cl. ......................................... 600/301;600/310
`
`(57)
`
`ABSTRACT
`
`A system includes an optical sensor and a signal processing
`module. The optical sensor is configured to be positioned on
`an area of skin of a patient. The optical sensor includes a light
`source for illuminating a capillary bed in the area of skin and
`a photodetector. The photodetector is configured to receive an
`optical signal from the capillary bed resulting from the illu-
`mination and to convert the optical signal into an electrical
`signal, the optical signal characterizing a fluctuation in a level
`ofblood in the capillary bed. The signal processing module is
`configured to process the electric signal using a nonstationary
`frequency estimation method to obtain a processed signal
`related to at least one of a heart rate and a respiration rate of
`the patient. Another aspect relates to obtaining a quantity
`related to the blood pressure of the patient in addition to or
`instead of obtaining a processed signal related to at least one
`of the heart rate and the respiration rate of the patient.
`
`
`
`126
`
`124
`
`001
`
`APLl 104
`
`Apple v.
`IPR2017-00318
`
`Valencell
`
`APL1104
`Apple v. Valencell
`IPR2017-00318
`
`001
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 1 0f 20
`
`US 2009/0105556 A1
`
`124 Processing Detector 126
`
`
`102
`
`Fig.1
`
`002
`
`002
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 2 0f 20
`
`US 2009/0105556 A1
`
` VON
`
`suunfiJemqla
`
`003
`
`N.3
`
`003
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 3 0f 20
`
`US 2009/0105556 A1
`
`‘i
`
`‘,
`
`mom
`
`
`
`4.mi259:08*me0M0_m_5ammo.ucmm
`
`
`.‘mi259:0cwfiwwcmmwa3533N:m.m-m.oEsmoEmEmE
` mum.EmcFm.
`
`m.GEOwzomcmEmqu55n—wmmluucmmW:
`
`mcmEmm
`
`vmm«pm
`
`canmomcom
`
`
`
`
`
`mcmmm25259:0cofiemfi2me_m_“E
`
`A
`
`mi959:0wEmEmmw
`.:o£me28
`
`
`
`0:9“HIP:_m_
`
`mNmmpm
`
`
`
`EmuoEmEm3mm:3:
`
`A.
`
`:mi233:0
`
` fifloEEmm
`
`:092mo_momzomcflcmmeEEmmmmficmm
`
`mgcoEmmm58:8:N:3-5.0
`
`
`coafifiwom
`
`m.9“.
`
`004
`
`004
`
`
`
`
`
`
`
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 4 0f 20
`
`US 2009/0105556 A1
`
`fiat?
`
` main.
`
`005
`
`005
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 5 0f 20
`
`US 2009/0105556 A1
`
`mmm
`
`Em
`
`can
`
`R.we:
`
`
`
`cofiufimfli=5;
`
`"‘“T
`
`--_;__.,....-
`
`___J--q_A___.l—_
`
` f
`
`006
`
`006
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 6 0f 20
`
`US 2009/0105556 A1
`
`
`
`
`
`5533+uEmEEmmEman:nmfigomam!535...
`
`
`
`
`
`m.9”.
`
`007
`
`007
`
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 7 0f 20
`
`US 2009/0105556 A1
`
`
`
`_7
`
`
`
`B=Emtm>fibuonmHERE.53.735
`
`
`
`
`
`
`
`
`
`
`
`BE:>995wd
`
`md......_..ono.0o77
`
`New
`
`n.3
`
`
`
`AMI$0-228“=7|.|_111:I.V
`
`777777
`
`com
`
`ova
`(89) time
`
`mtg:35cmmam
`
`
`
`7N:30-3.95‘i‘‘
`
`008
`
`008
`
`
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 8 0f 20
`
`US 2009/0105556 A1
`
`
`
`haawfimuoEmEmEhEm3EcEfiLEflE
`
`__+__—_q___..
`
`-(
`‘1:
`
`F.
`
`A:
`__p.....
`nun-J;
`
`on“
`
`uzuxsuuaalq
`
`m.9“.
`
`009
`
`009
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 9 0f 20
`
`US 2009/0105556 A1
`
`:2323391E3.8.5:.
`
`
`
`E>=mcmEm>_mm_m>2.:mm
`
`mgawmm
`
`3m:
`
`En:2%:3.3
`
`
`
`En:9mm33¢
`
`
`
`a.mi
`
`
`
`A13mE:_o>tam:
`
`2359:086.50
`
`010
`
`
`
`
`
`Aha:$5.me22205sz3;85m
`
`010
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 10 0f 20
`
`US 2009/0105556 A1
`
`3.;ng
`
`EmEm8m3“EnSm.
` 5.3$3325:
`
`E2::
`
`sum Mangqlz
`
`E.5
`
`011
`
`QQQF
`
`weep
`
`. I
`
`IlII
`
`|
`«
`‘D-W; -
`
`—-—-.|__.._.I____J——-—
`
`IIIII
`
`SE3.2
`
`,—
`
`,_
`
`011
`
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 11 0f 20
`
`US 2009/0105556 A1
`
`
`
`:.mE
`
`012
`
`012
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 12 0f 20
`
`US 2009/0105556 A1
`
`:9.
`
`3L Om c (
`
`Dm
`
`(D
`Q.
`E].
`
`1202
`
`:-'=5.
`
`{2.1:.Y,‘
`PJ
`
`it?)
`
`3:1
`
`activeclampingmE‘Chanism
`
`
`
`f;
`a
`i
`
`l:
`i
`
`5
`
`(.53
`r
`“"-'2';ir;mJ{r_-)z‘r
`
`
`
`g
`
`3
`E
`g.55
`
`N
`Y.
`
`013
`
`013
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 13 0f 20
`
`US 2009/0105556 A1
`
`2.2.
`
`m r.
`
`9“.
`
`mam“
`
`can»
`
`.805_m:mt<
`
`
`
`:EmSmomofizcmcoaxm
`
`
`
`.9828562me
`
`8magamo
`
`E26E82m
`
`6%9652,
`
`_m_8_>_
`
`SJ9385
`
`
`6:235-Asimt
`im.aa
`
`-Lowcmw
`
`whmh
`
`85:5"85:5
`
`014
`
`014
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 14 0f 20
`
`US 2009/0105556 A1
`
`
`
`
`
`:omtugooacct.«Emma:—uoofiwas“523:.“Montoya.
`
`
`
`
`
`
`
`
`
`“Ems:.522.$3.".fiemoeifiuafioE
`
`
`
`
`
`Purina.
`
`
`
`wI.I..ficwEmgsmmetonmm0:335
`
`s.
`
`3an
`
`mmmbmhe
`
`
`8.3;,_SS..8cm8%8.:
`Ia.i..E95...
`
`
`
`35E333E3333mflm>E.3335EELSE
`
`
`
`
`
`
`
`3553:23E3.933$32....522,33m0%.!!!
`
`
`
`
`
`
`
`
`
`“EmE9523E33.3;mEmSE.390...w“in.u
`
`
`
`flaw5333Etau“mm2.3»>m,‘,__,m,u,..,.,m
`
`
`
`92m5333Etau”2:32.”._a_._mt<cams.
`
`
`
`
`
`v
`
`s1uauiainsealu “n3 pun dpma
`JO aJoas mnpuals
`
`“$9.35
`
`323:.
`
`.6”
`m5
`
`mmwma
`xGas.
`
`"Mg.83
`
`E.3
`
`015
`
`015
`
`
`
`
`
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 15 0f 20
`
`US 2009/0105556 A1
`
`moms
`
`vamp
`
`
`
`<n_n_umoEBoUEmu.m.m
`
`
`
`
`
`5.2%.o...cozmoEzfiEoo538051
`
`
`
`
`
`momtmgLomcwm
`
`
`
`
`
`“E:mothum—m.233.51
`
`momr
`
`
`
`_o._EoDEmum>w
`
`van
`
`
`
`mom»new»
`
`
`
` womtflEEm:‘0.:m3-ucmLcot-ODD“mmmhowmEmamEmmmuEn..mcmEi@293.Il-whomcmw0...
`
`
`
`
`
`EmfiinzwEgon.
`
`EmEmmmcmE530.".ENBotmm
`
`“mm:o._.E->235
`
`
`
`,Npmhohm“
`
`
`
`2.9“.
`
`016
`
`016
`
`
`
`
`
`
`
`Patent Application Publication
`
`m2
`
`e
`
`61
`
`US 2009/0105556 A1
`
`3mmtme
`
`ofigzmm
`
`
`
`B25%:W“NewmmEn:mao_m>:m_335
`
`
`
`2mmcozmtammmmsomcmEEmE
`
`E“BmEzmmEye55:59;
`
`o_§_mc<
`
`8833:.
`
`52Emmmuficmm
`
`35
`
`22.
`
`VS“
`
`NE:
`
`ammmh
`
`vwcp
`
`«New
`
`92mmacumtamommsomcmEEwE
`
`
`
`
`%Nu"BmEzmmEN:5532“.
`
`o=>_mc<
`
`Egewcmfi
`
`
`
`55;wmmaficmm
`
`ES
`
`
`
` mugEmanmamenm3mm:o:m.._nmmmAmmo55:69“.M9033535
`
`3H:
`
`dew
`
`25mg
`
`F2833...
`
`
`
`55“.mmmqficmm
`
`3:
`
`mach
`
`vowF
`
`New“
`
`wsomcficflwE
`
`EIV>ocmzum£
`
`23.9%
`
`E535;
`
`mmmuécmm
`
`BE“.
`
`1::
`
`mun—n.
`
`lace“
`
`E.E
`
`017
`
`017
`
`
`
`
`
`
`
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 17 0f 20
`
`US 2009/0105556 A1
`
`
`nmflmEzwm..2[29:002TE..3.A
`
`
`
`
`tmmI“El.w>zm>=mm29853:55m_m:<meEoO“Em
`2mmNINo-01NIm.m-m.oz
`
`TNx69“.$512005025
`82¢32l«a:
`
`
`
`€829...tmEE
`
`:o_EE_xo.aa<
`
`a
`
`“:0mun—n—
`
`
`
`
`
`
`
`-23mm.53:...13mm.632.m>m_mn_Emmi;‘N>m_mn_Emmi;
`+ITIII22mm.cm?Wm.-
`
`2.3.
`
`mobF35h
`
`EEG
`
`52:
`
`018
`
`t.5
`
`018
`
`
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 18 0f 20
`
`US 2009/0105556 A1
`
`who“2.2.
`
`8883mm.
`
`3mmcozfiawmm
`
`NImodéz
`
`
`
`
`mm:295waE00
`
`
`
`m>wm>zmo99353.5::2waano2m-mE_-__mm
`
`
`
`5.22mi.tmnzr
`
`“Em.
`
`NImdéodz
`
`E9mmtum...
`
`50m.metmmI
`
`So0.5.flmo08.P-0?IN.aacosz_xoaa<E9.:>28
`.£660562:.iflmnmeBE£560Emmi:N>m'_mn_._NwmmE_
`
`
`
`
`
`
`
`019
`
`ONE.
`
`8.9;
`
`019
`
`
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009
`
`Sheet 19 of 20
`
`US 2009/0105556 A1
`
`SOmun—n—
`
`
`
`.F:o_ng_xo.:E<
`
`SOwe:22m3.--T-
`
`mwoh_
`
`umwmémwm_._mE.on
`
`
`
`:ozmtamwman:u>IE>=wQ2285agc:S”MW—thooSNWMHMEMHM00in2meE002mac:Em.”Em
`
`3mma.4.fl\“EmmuzzcmmimeEooN:m.wrm.a2
`NIM05102NIm.o.mo.ol
`8‘“New
`IIfl:2ai.
`SHIINVxFEIED._Nmmmt:
`
`
`
`
`Eoficmfi55:1
`
`QNQFVwehNE:
`
`
`
`
`
`Eoucmc.tunzIcozmgxeag
`
`
`
`
`
`
`
`
`
`932mm53::9.2mmEmmi...REED._maEEm>m_mn_.5me1233?:92mm.53:...
`
`
`
`
`
`
`
`
`
`
`
`2.3
`
`020
`
`020
`
`
`
`
`
`
`Patent Application Publication
`
`Apr. 23, 2009 Sheet 20 0f 20
`
`US 2009/0105556 A1
`
`BINEE
`
`
`
`I\..
`
`82N2:
`
`
`
`EDP/«DuI32IN.«2:
`
`_oo._.<n_.._
`
`cozmtqmmmLn:m_mc<meEoo“Em
`
`
`
`
`
`3mmNIno.0-ozNImdlmodz
`
`
`
`
`m>zm>zmo$9853.3::oyx-mEEooBmmE_-.mmm
`N»E_.829.
`
`Nv_29200Bme
`
`
`
`Eoficmfitwnzz
`
`
`
`
`
`“3.0Gan.:—a9660[092138[$95928lwmmfi:N329“.NINmmBEm:_._acosz_ona<
`
`
`
`
`
`
`
`021
`
`03:
`
`om.9“.
`
`021
`
`
`
`
`
`
`US 2009/0105556 A1
`
`Apr. 23, 2009
`
`MEASUREMENT OF PHYSIOLOGICAL
`SIGNALS
`
`CROSS-REFERENCE TO RELATED
`APPLICATIONS
`
`[0001] This application claims priority to US. provisional
`application No. 60/995,723, filed Sep. 28, 2007, entitled
`“Method and Devices for Measurement of Multi-modal
`
`Physiological Signals,” which is incorporated herein by ref-
`erence.
`
`STATEMENT REGARDING FEDERALLY
`SPONSORED RESEARCH
`
`[0002] The subject matter described in this application was
`partially funded by the Government of the United States
`under Contract No. W9lZLK-04-P-0239 awarded by the
`US. Department of the Army. The government has certain
`rights in the invention.
`
`FIELD OF THE INVENTION
`
`[0003] The invention relates to measurement ofphysiologi-
`cal signals.
`
`BACKGROUND
`
`Physiological signals are important for monitoring a
`[0004]
`subject’s physical and cognitive state. Often, heart rate
`parameters are measured directly via electrocardiogram
`(ECG) measurements of a heart beat. Respiration rate data
`can be obtained from a respiration chest strap. Physiological
`signals can also be extracted from infrared (IR) photoplethys-
`mographs (PPG). The signals of interest include heart rate,
`respiration rate, continuous blood pressure, and intrathoracic
`pressure. With respect to blood pressure, there is technology
`related to collecting data at two locations on the body and
`using pulse transit time and other parameters as the basis of
`the pressure estimate.
`
`SUMMARY
`
`In a general aspect, a system includes an optical
`[0005]
`sensor and a signal processing module. The optical sensor is
`configured to be positioned on an area of skin of a patient. The
`optical sensor includes a light source for illuminating a cap-
`illary bed in the area of skin and a photodetector. The photo-
`detector is configured to receive an optical signal from the
`capillary bed resulting from the illumination and to convert
`the optical signal into an electrical signal, the optical signal
`characterizing a fluctuation in a level ofblood in the capillary
`bed. The signal processing module is configured to process
`the electric signal using a nonstationary frequency estimation
`method to obtain a processed signal related to at least one of
`a heart rate and a respiration rate of the patient.
`[0006] Embodiments may include one or more of the fol-
`lowing. The system includes an output for providing infor-
`mation determined from the processed signal. The nonsta-
`tionary frequency estimation method includes a Hilbert
`transform method or an instantaneous frequency estimation
`method. The processed signal includes at least one of instan-
`taneous heart rate, inter-beat interval, heart rate variability,
`high-low heart rate ratios, respiration rate, inter-breath inter-
`val, and respiration rate variability. The fluctuation in the level
`ofblood in the capillary bed relates to a change in at least one
`
`ofvolume and pressure ofthe thoracic cavity or to a change in
`at least one ofvolume and pressure of an organ in the thoracic
`cavity.
`[0007] The system includes an auxiliary sensor configured
`to detect an ambient signal. The auxiliary sensor includes at
`least one of accelerometer, a pressure sensor, an optical detec-
`tor, a temperature sensor, and a piezoelectric device. The
`signal processing module is configured to remove an effect of
`the ambient signal from the electrical signal. The optical
`signal is a reflectance or a transmittance of the capillary bed.
`[0008]
`In another general aspect, a method includes illumi-
`nating a capillary bed in an area of skin of a patient, receiving
`an optical signal from the capillary bed resulting from the
`illumination, converting the optical signal into an electrical
`signal, and processing the electrical signal using a nonstation-
`ary frequency estimation method to obtain a processed signal
`related to at least one of a heart rate and a respiration rate of
`the patient. The optical signal characterizes a fluctuation in a
`level of blood in the capillary bed.
`[0009] Embodiments may include one or more of the fol-
`lowing. The method includes outputting information deter-
`mined from the processed signal. Processing the electrical
`signal using the nonstationary frequency estimation method
`includes performing a Hilbert transform or processing the
`electrical signal using an instantaneous frequency estimation
`method. Processing the electrical signal using the instanta-
`neous frequency method includes band pass filtering the elec-
`trical signal, determining an instantaneous frequency of the
`electrical signal, and using the instantaneous frequency to
`obtain the processed signal.
`[0010] The method further includes processing the electri-
`cal signal using a model to obtain a blood pressure signal
`related to a blood pressure of the patient. The optical signal
`characterizes a capillary refill time in the capillary bed. Pro-
`cessing the electrical signal includes processing the electrical
`signal in real time.
`[0011]
`In another aspect, a method for monitoring blood
`pressure includes illuminating a capillary bed in an area of
`skin of a patient, receiving an optical signal from the capillary
`bed resulting from the illumination, converting the optical
`signal into an electrical signal, and processing the electrical
`signal using a model characterizing a relationship ofthe fluc-
`tuation in the level of blood and the blood pressure of the
`patient to obtain a quantity related to the blood pressure ofthe
`patient. The optical signal characterizes a fluctuation in a
`level of blood in the capillary bed of the patient.
`[0012] Embodiments may include one or more of the fol-
`lowing. The method includes outputting information deter-
`mined based on the quantity related to the blood pressure of
`the patient. The optical signal characterizes a capillary refill
`time. The method further includes engaging a device to
`restrict circulation in the capillary bed of the patient and
`disengaging the device prior to receiving the optical signal
`from the capillary bed. The disengaging of the device occurs
`gradually. The device is an active clamping device.
`[0013] The quantity related to the blood pressure of the
`patient is a quantity related to the continuous blood pressure
`ofthe patient. Applying the model includes applying a model
`including circuit elements or properties of the capillary bed.
`The method further includes calibrating the model on the
`basis of a blood pressure of the patient determined by using a
`blood pressure cuff.
`
`022
`
`022
`
`
`
`US 2009/0105556 A1
`
`Apr. 23, 2009
`
`[0014] Embodiments may include one or more of the fol-
`lowing advantages.
`[0015] A system or method as described above can be used
`for both military and civilian applications. Combat casualty
`care requires close monitoring ofvital signs from the moment
`that a medic first attends to a wounded soldier in the battle-
`
`field and thence through the many transfer stages to the point
`of full hospital care, generally removed from the combat
`scene. Similar needs are evident in the civilian community
`where critical care is administered by first responders at the
`scene of accidents, by emergency room staff, and by intensive
`care unit staff. It is often desirable to obtain maximum infor-
`
`mation using as little equipment as possible. The system and
`method described herein support this need. They reduce the
`burden of equipment logistics, the burden of extra wires and
`sensors on and around the patient, and the complexity and
`cost of using multiple devices.
`[0016]
`For both military and civilian applications, a dispos-
`able, wearable device in keeping with the system and method
`described herein can be adapted to stay with a patient and to
`report vital signs throughout the care and transport processes.
`Further, the system can be configured to provide medical
`personnel with real-time visibility of vital signs as well as
`recording ofthis information for playback by attending medi-
`cal staff at a later time. The disposability of the device allows
`it to be fabricated with low cost parts and eliminates the need
`for sanitization and asset tracking logistics in large scale
`clinical or military uses.
`[0017]
`Such system and methods additionally support
`applications in fitness monitoring, where their ease ofuse and
`robustness make them a compelling alternative to chest strap
`monitors for the monitoring ofcardiac and respiratory param-
`eters during exercise. An ear-wom device can also integrate a
`speaker unit for mobile electronic devices such as mobile
`phones or music players.
`[0018] An advantage of applying a nonstationary fre-
`quency estimation method (e. g., analysis involving monitor-
`ing the frequency changes of the signal over time, such as
`monitoring changes in the instantaneous principal frequency
`over time) is that it is possible to avoid a tradeoff inherent in
`many stationary estimation methods between frequency reso-
`lution and duration of data signals being analyzed. For
`example, ifthe signal is assumed to be stationary within each
`of a series of data windows, the frequency resolution is gen-
`erally inversely proportional to the duration of the window.
`As the window duration increases, the assumption of a sta-
`tionary signal is increasingly violated and/or nonstationary
`events (e.g., transients) are more difficult to detect. At least
`some nonstationary frequency analysis methods, which may
`be based, without limitation, on a Hilbert transform approach,
`tracking of a nonstationary model, nonstationary principal
`frequency analysis, or other time-frequency methods, miti-
`gate the effects of such a time-frequency tradeoff. Further-
`more, use of such nonstationary techniques, as opposed to use
`of time domain peak picking and/or threshold based tech-
`niques, can provide robustness of algorithm against artifacts,
`and provide sensitivity to periodicity without being burdened
`by a window that can reduce the time resolution.
`[0019] Other features and advantages are apparent from the
`following description and from the appended claims.
`
`BRIEF DESCRIPTION OF DRAWINGS
`
`FIG. 1 is a schematic diagram of a photoplethysmo-
`[0020]
`graph (PPG) sensor system.
`
`FIG. 2 is a graph of a PPG detector signal taken over
`[0021]
`a 25 second period by an earlobe PPG device.
`[0022]
`FIG. 3 is a flow diagram of signal processing ofa
`detector signal from a PPG device to obtain heart rate and
`respiration rate parameters.
`[0023]
`FIG. 4 is a graph of a result ofband-pass filtering the
`data shown in FIG. 2 between 0.5 Hz and 5.5 Hz to extract a
`
`cardiac signal.
`[0024]
`FIG. 5 is a graph of a result ofband-pass filtering the
`data shown in FIG. 2 between 0.17 Hz and 0.5 Hz to extract a
`
`respiration signal.
`[0025]
`FIG. 6 is a graph ofan inter-beat interval obtained by
`applying an instantaneous frequency method to the cardiac
`signal shown in FIG. 4.
`[0026]
`FIG. 7 is a graph of a spectral analysis of the inter-
`beat interval data shown in FIG. 6.
`
`FIG. 8 is a graph of the respiration rate obtained by
`[0027]
`applying an instantaneous frequency method to the respira-
`tion signal shown in FIG. 5.
`[0028]
`FIG. 9 is a diagram of PPG measurements related to
`physiological states used to determine intrathoracic pressure.
`[0029]
`FIG. 10 is a graph of the output of a matched filter-
`ing process using the PPG detector signal shown in FIG. 2 and
`a pulse pilot signal.
`[0030]
`FIG. 11 is a block diagram of a least mean squares
`(LMS) adaptive filter.
`[0031]
`FIG. 12 is a schematic diagram of an active clamp-
`ing mechanism used to stimulate capillary refill.
`[0032]
`FIG. 13 is a diagram of a system model relating a
`PPG signal to blood pressure.
`[0033]
`FIG. 14 is a graph oftrends in various physiological
`parameters before and during a stress event.
`[0034]
`FIG. 15 is a block diagram of a portable electronics
`unit.
`
`FIG. 16 is a flow diagram of methods to estimate a
`[0035]
`heart rate and a respiration rate.
`[0036]
`FIG. 17 is a flow diagram of a processing delay in
`the estimation of a heart rate.
`
`FIG. 18 is a flow diagram ofa processing delay in a
`[0037]
`first method for the estimation of a respiration rate.
`[0038]
`FIG. 19 is a flow diagram ofa processing delay in a
`second method for the estimation of a respiration rate.
`[0039]
`FIG. 20 is a flow diagram ofa processing delay in a
`third method for the estimation of a respiration rate.
`
`DETAILED DESCRIPTION
`
`[0040] Referring to FIG. 1, examples of an infrared pho-
`toplethysmograph (PPG) device 100 are used to obtain physi-
`ological signals related to one or more of heart rate, respira-
`tion rate, blood pressure, and intrathoracic pressure. Such
`signals may be relevant for monitoring a person’s state,
`including one or more of the person’s physical state, long-
`terrn health, psychological state, and/or cognitive state. More
`generally, the physiological signals may provide information
`about the activity of the person’s sympathetic and parasym-
`pathetic nervous system. The PPG device 100 illustrated in
`FIG. 1 is attached to an earlobe 102 of a person, for example,
`using a clamping or adhesive approach. However, in other
`embodiments, PPG device 100 is used on other areas of the
`skin of a person, including but not limited to a portion of a
`forehead, a neck, an arm, a forearm, a finger, a leg, a back, an
`abdomen, or a stomach. In general, a requirement for the
`positioning of PPG device 100 is that the PPG sensor be
`located such that it can obtain a measurement via the skin that
`
`023
`
`023
`
`
`
`US 2009/0105556 A1
`
`Apr. 23, 2009
`
`is related to blood flow or pressure, for example to measure a
`level of blood in a capillary bed 104, for example, a blood
`volume, a rate of blood flow, or a rate of change of blood
`volume. Note also that the approach is not limited to use of a
`single PPG device on an individual. In some embodiments,
`multiple PPG devices are used, for example, on the torso
`and/or at different extremities, and signals obtained at the
`different PPG devices are processed independently or in com-
`bination to determine underlying characteristics of the indi-
`vidual’s state.
`
`In some embodiments, such as that shown in FIG. 1,
`[0041]
`an infrared light source 106 illuminates the earlobe 102. The
`blood level in capillary bed 104 affects the amount of light
`108 that is backscattered or reflected by earlobe 102. Light
`108 backscattered by earlobe 102 is received by an optical
`transducer such as a photodetector 110 and converted into a
`detector signal 112. Since the blood flow in capillary bed 104
`is controlled by the heart beat ofthe person and thus the blood
`level in the capillary bed varies with time, the backscattered
`light 108 and hence the detector signal 112 are also time-
`varying. In another embodiment, the PPG sensor operates in
`transmission mode and the light transmitted through the cap-
`illary bed is received by the photodetector.
`[0042] The detector signal 112 is sent to a signal processing
`unit 114 which processes the detector signal, which contains
`information about the person’s pulse, to extract desired physi-
`ological data, in various embodiments including one or more
`of instantaneous heart rate, inter-beat interval, heart rate vari-
`ability, high-low heart rate ratio, respiration rate, inter-breath
`interval, respiration rate variability, blood pressure, and
`intrathoracic pressure. A single PPG device 100, referred to
`below as an Integrated Multi-Modal Physiological Sensor
`(IMMPS), is capable of producing multiple (or all) of such
`types of physiological data.
`[0043]
`In some embodiments, the PPG device 100 provides
`real-time visibility of physiological parameters and vital
`signs, which can be transmitted to other equipment for real-
`time processing or for playback or off-line processing at a
`later time. In some embodiments, the PPG device includes
`user output devices, such as a set of light emitting diodes
`(LEDs) (e.g., a red LED 116, a yellow LED 118, and a green
`LED 120) or an audio device for producing alert sounds,
`which provide on-device status on PPG device 100. As an
`example for use of such output devices, when a selected
`physiological parameter is in a normal range, green LED 120
`is turned on; when the physiological parameter is in a slightly
`abnormal range, yellow LED 118 is turned on; when the
`physiological parameter is in a dangerous range, red LED 116
`is turned on. In some embodiments, the audio output device is
`used to provide other audio output, such as the output for an
`electronic device such as a mobile phone or a music player. In
`some embodiments, a wireless link 122 to an external moni-
`toring system 124, such as a bedside system or a wearable
`system, provides sensor data to the external system enabling
`a numeric readout 126 of various physiological parameters.
`In some embodiments, the PPG device, or at least some
`wearable portion of the device, is disposable. In such dispos-
`able embodiments, the bedside system can be designed to be
`sterilized and reused; in another embodiment, the bedside
`system itself is also disposable. In some embodiments, the
`bedside system includes or communicates with a centralized
`monitoring system that monitors PPG devices of multiple
`patients.
`
`the photodetector based
`In some embodiments,
`[0044]
`detector signal is augmented with other signals, for example,
`accelerometer or pressure sensor signals. For example, aux-
`iliary sensors 130 are connected to signal processing unit 114
`via a wired connection 132. In other embodiments, auxiliary
`sensors 130 are connected to signal processing unit 114 via a
`wireless connection. Auxiliary sensors 130, such as tempera-
`ture sensors, accelerometers, pressure transducers, optical
`detectors, or piezoelectric films or matrices can provide aux-
`iliary signals 132 related to ambient sources ofnoise to signal
`processing unit 114. Signal processing unit 114 incorporates
`auxiliary signals 132 into the signal processing, for example,
`to increase the signal-to-noise ratio ofthe desired physiologi-
`cal data.
`
`Heart and Respiration Rate Signals
`
`[0045] Referring to FIG. 2, in an embodiment of the PPG
`device that uses a photodetector signal obtained from at an
`earlobe location, a detector signal 200 obtained from the PPG
`device 100 has a high-frequency pulse signal 202 whose local
`peaks in the time domain have a one-to-one correspondence
`with cardiac beats. Significant low-frequency amplitude vari-
`ability in the detector signal is due in large part a respiration
`signal 204, which modulates the baseline of the pulse signal.
`[0046] The time varying heart and respiration components
`in detector signal 200 can be modeled as
`3(1):AH(I)COS [OJH(Z)I+¢H(T)]+AR(I)COS [mR(l+¢R(T)]+
`N0),
`
`(1)
`
`where AH(t) is the amplitude modulation of pulse signal 202,
`mH(t) is the frequency modulation ofpulse signal 202, ¢H(t) is
`the phase modulation of pulse signal 202, AR(t) is the ampli-
`tude modulation of respiration signal 204, (nR(t) is the fre-
`quency modulation of respiration signal 204, ¢R(t) is the
`phase modulation of respiration signal 204, and N(t) is the
`time varying noise, which includes baseline drift and broad-
`band noise in the overall signal band.
`[0047] Given the measured detector signal 200 (s(t)), signal
`processing is performed to estimate the slowly varying com-
`ponents of the heart rate mH(t) and the respiration rate (nR(t).
`It is known that u)H(t)zu)H031 beat per second 1 Hz for heart
`rate and u)R(t)zu) 10312 breaths per minutez0.2 Hz for respi-
`ration rate.
`
`[0048] Amplitude, phase, and frequency modulation cause
`spectral spread that broadens the pure tones implied by these
`frequencies. Amplitude and phase modulation and rapid fluc-
`tuations of the frequency modulation are confounding com-
`ponents of detector signal 200. The slowly varying compo-
`nents of u)(t) are the desired components for obtaining heart
`and respiration rate information.
`[0049]
`For both heart and respiration rate, u)(t), is com-
`posed of three parts: constant frequency 000, which is the
`nominal heart or respiration rate; a zero-mean, slowly varying
`frequency component S25 having a time scale of minutes; and
`a zero -mean, rapidly varying frequency component thaving
`a time scale of seconds. In this case, the composite heart rate
`or respiration rate is written as
`
`m(z):030+£25(z)+£21(1).
`
`(2)
`
`The signal of interest is the combination of the constant and
`the slowly varying component (for both heart rate and respi-
`ration rate) and is written as
`
`03(1):030+£25(z).
`
`(3)
`
`024
`
`024
`
`
`
`US 2009/0105556 A1
`
`Apr. 23, 2009
`
`The phase modulation noted is Eq. (1), ([)(t), is assumed to be
`small, since large phase modulation can be represented as
`frequency modulation and this is already captured in £25 or 92f.
`[0050] Applying Eq. (3) to Eq. (1), a new equation for the
`measured raw signal 200 is obtained:
`
`a<z>:AH<z>cos [&H<z>z+¢.H<z>]+AR<z>cos [
`mR<z>z+¢.R<z>]+N<zx
`
`<4)
`
`where for both heart and respiration rate, small components
`¢E(t) are defined as
`
`¢E(l):lioo’9fi)dt+¢(l)
`
`(5)
`
`and by construction |¢E(t)l<<l and ¢E(t) is zero mean.
`[0051] Using the Law of Cosines and the fact that |¢E(t)
`|<<l, the following equation for raw signal 200 is obtained:
`
`s<z>:4H<z>{ cos [&H<z>u;¢.y<z>sin [03w>z]}+AR<z>{
`[mR<z>z]—¢.R<z>sin [wR<z>z]}+N<z>.
`
`<6)
`
`[0052] The formulation in Eq. (6) of PPG detector signal
`200 suggests a number of methods to estimate the desired
`slowly varying heart and respiration rate signals (hH(t) and
`(1)R(t), respectively. Such estimation techniques can include
`instantaneous frequency determination Via analytic signals,
`moving averages, band pass filtering, synchronous detection,
`correlation detection, narrowband processes (e.g., demodu-
`lation), matched filtering, wavelet filtering, short-time fre-
`quency analysis (e.g., short-time fast Fourier transform,
`Wigner-Ville transform), state estimation (e.g., Kalman fil-
`tering, unscented filtering), Doppler processing, or a combi-
`nation of the above methods. A number of these techniques
`can be implemented to account for the nonstationary nature of
`the detector signal, which relates to the time variation of the
`frequency modulation signals, u)(t).
`[0053] Referring to FIG. 3, one example of a procedure for
`obtaining heart rate parameters and/or respiration rate param-
`eters involves receiving the detector signal 200, such as that
`shown in FIG. 2, from a PPG photodetector (IR Plethysmo-
`graph 300). To emphasize the pulse signal in order to extract
`cardiac parameters from the detector signal, the respiratory
`modulation is removed by band-pass filtering (302) the detec-
`tor signal between approximately 0.5 Hz and approximately
`5.5 HZ. Band-pass filtering of the detector signal shown in
`FIG. 2 results in a cardiac signal shown in FIG. 4, which has
`significantly less amplitude variability due to respiratory con-
`tamination. The band-pass filtered signal pulse signal can be
`expressed as
`
`§<z>:AH<z>{ cos [&H<z>U—¢.H<z>sin [Mammal
`
`(7)
`
`where N(t) is the content of the noise, N(t), within the heart
`rate bandwidth. Likewise, referring again to FIG. 3, to extract
`parameters related to respiration from a PPG detector signal,
`the pulse signal is removed by band-pass filtering (304) the
`detector signal between approximately 0.17 Hz (equivalent to
`10 breaths per minute) and 0.5 Hz (30 breaths per minute).
`Band-pass filtering the detector signal of FIG. 2 produces a
`respiration signal shown in FIG. 5, which retains primarily
`low frequency respiration components of the original signal.
`[0054]
`Focusing now on the cardiac signal, one approach to
`detecting heart beats in the cardiac signal is by threshold-
`based peak picking, which can be used to determine the time
`of specific events such as heart beats. In some implementa-
`tions, peak picking can be sensitive to the threshold selected.
`For example, if the threshold is set too low, then false beats
`can be counted, and the inter-beat interval (IBI) is determined
`to be shorter than it really is. If the threshold is set too high,
`
`then true beat peaks may be missed. For general processing it
`is convenient to select a single threshold, but if the data are
`sufficiently variable then selecting a single threshold is not
`possible.
`[0055] Other approaches to processing the cardiac signal
`use nonstationary frequency estimation methods. For
`instance, instantaneous frequency computation component
`306, 318 are used to determine the time-varying principal
`frequencies present in bandpass filtered signals. Approaches
`to instantaneous frequency computation include Hilbert
`Transform methods, which are particularly effective because
`the band-pass filtered cardiac signal shown in FIG. 4 has a
`strong sinusoidal content and because systematic changes in
`the IBI over long periods of time are often of interest. Model-
`based approaches, such as Kalman filtering approaches
`described later in this document can also be used, with state
`variables (or alternatively time-varying system parameters)
`that are estimated corresponding to the instantaneous fre-
`quencies of interest. Yet other approaches to tracking instan-
`taneous frequency can be based on adaptive model