US20190343406A1 - System and Method for Noninvasively Measuring Ventricular Stroke Volume and Cardiac Output - Google Patents
System and Method for Noninvasively Measuring Ventricular Stroke Volume and Cardiac Output Download PDFInfo
- Publication number
- US20190343406A1 US20190343406A1 US16/522,487 US201916522487A US2019343406A1 US 20190343406 A1 US20190343406 A1 US 20190343406A1 US 201916522487 A US201916522487 A US 201916522487A US 2019343406 A1 US2019343406 A1 US 2019343406A1
- Authority
- US
- United States
- Prior art keywords
- state
- dssm
- model
- time
- computer implemented
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 98
- 230000000747 cardiac effect Effects 0.000 title claims abstract description 9
- 230000002861 ventricular Effects 0.000 title claims description 4
- 238000013178 mathematical model Methods 0.000 claims abstract description 14
- 210000000748 cardiovascular system Anatomy 0.000 claims abstract description 7
- 230000008569 process Effects 0.000 claims description 41
- 239000008280 blood Substances 0.000 claims description 26
- 210000004369 blood Anatomy 0.000 claims description 25
- 238000005259 measurement Methods 0.000 claims description 25
- 238000012545 processing Methods 0.000 claims description 20
- 239000013598 vector Substances 0.000 claims description 16
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 14
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 13
- 230000009977 dual effect Effects 0.000 claims description 13
- 229910052760 oxygen Inorganic materials 0.000 claims description 13
- 239000001301 oxygen Substances 0.000 claims description 13
- 238000005315 distribution function Methods 0.000 claims description 9
- 230000036772 blood pressure Effects 0.000 claims description 7
- 230000036581 peripheral resistance Effects 0.000 claims description 6
- 230000035790 physiological processes and functions Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000036387 respiratory rate Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 2
- 210000000653 nervous system Anatomy 0.000 claims description 2
- 230000001457 vasomotor Effects 0.000 claims description 2
- 230000000007 visual effect Effects 0.000 claims description 2
- 230000036962 time dependent Effects 0.000 claims 2
- 210000005259 peripheral blood Anatomy 0.000 claims 1
- 239000011886 peripheral blood Substances 0.000 claims 1
- 238000009826 distribution Methods 0.000 description 16
- 238000013186 photoplethysmography Methods 0.000 description 16
- 238000002565 electrocardiography Methods 0.000 description 14
- 210000001519 tissue Anatomy 0.000 description 12
- 238000002106 pulse oximetry Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 10
- 238000004422 calculation algorithm Methods 0.000 description 7
- 206010019280 Heart failures Diseases 0.000 description 6
- 239000003814 drug Substances 0.000 description 6
- 229940079593 drug Drugs 0.000 description 6
- 239000002245 particle Substances 0.000 description 6
- 206010007559 Cardiac failure congestive Diseases 0.000 description 5
- 210000001147 pulmonary artery Anatomy 0.000 description 5
- 102000001554 Hemoglobins Human genes 0.000 description 4
- 108010054147 Hemoglobins Proteins 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 230000035479 physiological effects, processes and functions Effects 0.000 description 4
- 238000003672 processing method Methods 0.000 description 4
- 210000000779 thoracic wall Anatomy 0.000 description 4
- 239000000654 additive Substances 0.000 description 3
- 230000000996 additive effect Effects 0.000 description 3
- 210000002376 aorta thoracic Anatomy 0.000 description 3
- 230000004087 circulation Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000010801 machine learning Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 230000009885 systemic effect Effects 0.000 description 3
- 238000012952 Resampling Methods 0.000 description 2
- 238000002835 absorbance Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000013477 bayesian statistics method Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 230000008081 blood perfusion Effects 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 210000003743 erythrocyte Anatomy 0.000 description 2
- 210000003238 esophagus Anatomy 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000004089 microcirculation Effects 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 230000036284 oxygen consumption Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000002685 pulmonary effect Effects 0.000 description 2
- 230000004043 responsiveness Effects 0.000 description 2
- 238000013175 transesophageal echocardiography Methods 0.000 description 2
- 241001535291 Analges Species 0.000 description 1
- 206010007687 Carotid artery stenosis Diseases 0.000 description 1
- 208000028399 Critical Illness Diseases 0.000 description 1
- 208000032843 Hemorrhage Diseases 0.000 description 1
- 208000001953 Hypotension Diseases 0.000 description 1
- 206010021137 Hypovolaemia Diseases 0.000 description 1
- 206010021143 Hypoxia Diseases 0.000 description 1
- 208000004880 Polyuria Diseases 0.000 description 1
- 208000004301 Sinus Arrhythmia Diseases 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000000137 annealing Methods 0.000 description 1
- 210000000709 aorta Anatomy 0.000 description 1
- 210000001765 aortic valve Anatomy 0.000 description 1
- 230000004872 arterial blood pressure Effects 0.000 description 1
- 210000002565 arteriole Anatomy 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 210000003403 autonomic nervous system Anatomy 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000002902 bimodal effect Effects 0.000 description 1
- 238000001574 biopsy Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000001715 carotid artery Anatomy 0.000 description 1
- 208000006170 carotid stenosis Diseases 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013075 data extraction Methods 0.000 description 1
- 238000013501 data transformation Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000035619 diuresis Effects 0.000 description 1
- 238000002651 drug therapy Methods 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001839 endoscopy Methods 0.000 description 1
- 230000012953 feeding on blood of other organism Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 235000012631 food intake Nutrition 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000004190 glucose uptake Effects 0.000 description 1
- 230000003862 health status Effects 0.000 description 1
- 238000005534 hematocrit Methods 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 208000021822 hypotensive Diseases 0.000 description 1
- 230000001077 hypotensive effect Effects 0.000 description 1
- 230000007954 hypoxia Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000000543 intermediate Substances 0.000 description 1
- 210000005240 left ventricle Anatomy 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005399 mechanical ventilation Methods 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 230000003278 mimic effect Effects 0.000 description 1
- 238000002324 minimally invasive surgery Methods 0.000 description 1
- 210000004115 mitral valve Anatomy 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 230000004220 muscle function Effects 0.000 description 1
- 210000004165 myocardium Anatomy 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000010412 perfusion Effects 0.000 description 1
- 230000036544 posture Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000004088 pulmonary circulation Effects 0.000 description 1
- 230000000541 pulsatile effect Effects 0.000 description 1
- 230000035485 pulse pressure Effects 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 210000003491 skin Anatomy 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000005654 stationary process Effects 0.000 description 1
- 238000002948 stochastic simulation Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 210000000115 thoracic cavity Anatomy 0.000 description 1
- 238000002627 tracheal intubation Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 210000005166 vasculature Anatomy 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0295—Measuring blood flow using plethysmography, i.e. measuring the variations in the volume of a body part as modified by the circulation of blood therethrough, e.g. impedance plethysmography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring 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/1455—Measuring 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/14551—Measuring 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6813—Specially adapted to be attached to a specific body part
- A61B5/6825—Hand
- A61B5/6826—Finger
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
-
- G06K9/00523—
-
- G06K9/6226—
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/021—Measuring pressure in heart or blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/0816—Measuring devices for examining respiratory frequency
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/7405—Details of notification to user or communication with user or patient ; user input means using sound
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/7455—Details of notification to user or communication with user or patient ; user input means characterised by tactile indication, e.g. vibration or electrical stimulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
Definitions
- the present invention relates generally to apparatus, systems, and methods for noninvasively measuring cardiac output (CO) and left ventricular stroke volume (SV). More specifically, the invention includes systems and methods that measure CO and/or SV using plethysmographic waveform data collected by a pulse oximeter.
- CO cardiac output
- SV left ventricular stroke volume
- SV left ventricular stroke volume
- CO cardiac output
- Patients being treated for heart failure are normally medicated with drugs that regulate diuresis and heart muscle function.
- An important aim of drug therapy is to maintain a CO that is sufficient to perfuse tissues with oxygenated blood. It is advantageous to use the lowest possible doses of drugs to manage CHF because the drugs used produce unwanted side effects.
- To optimally manage the dosage and selection of drugs in the treatment of CHF one must monitor CO to assess the efficacy of the drugs and dosages being administered and/or to monitor patient compliance.
- CO is measured using invasive techniques such as the Fick method, the Thermodilution method, and implantable microelectromechanical devices (MEMs), also called CardioMEMs.
- the Fick method involves the measurement of oxygen consumption and computing the arteriovenous difference using samples of arterial blood and mixed venous blood from the pulmonary artery.
- the thermodilution method measures the rate at which cold saline solution is diluted in the blood. Both of these methods are performed in a hospital setting because they require the placement of a catheter in the pulmonary artery. Additionally, the use of a pulmonary artery catheter use may also increase morbidity in critically ill patients.
- CardioMEMs are surgically implanted into patients in a hospital setting and, once implanted, provide measurements of SV and CO that can be used to monitor patients.
- the implantation of the cardioMEMs device into the pulmonary artery is expensive and is performed in a hospital setting and involves the risks associated with heart catheterization.
- Esophageal Doppler measures blood flow velocity in the descending thoracic aorta using a flexible ultrasound probe that is inserted into the esophagus. The blood flow velocity is combined with an estimate of the cross-sectional area of the aorta estimated from the patient's age, height, and weight to calculate SV and CO. This technique requires someone with technical skill to insert an esophageal Doppler monitor, which must be properly aligned with respect to the thoracic aorta to provide accurate measurements.
- Transesophageal echocardiography involves measuring SV using flow velocity calculated from the area under the measured Doppler velocity waveform at the pulmonary artery, the mitral valve, or the aortic valve. This technique requires a highly trained operator to position and place the esophageal Doppler monitor. These procedures are not truly noninvasive because accessing the esophagus is perceived by patients as invasively uncomfortable.
- a pulse oximeter is a device that obtains photoplethysmography (PPG) data, which measures changes in blood volume within a tissue caused by the pulse of blood pressure through the vasculature in the tissue.
- PPG photoplethysmography
- the blood volume change is detected by measuring the amount of light transmitted or reflected to a sensor from a light source used to illuminate the skin.
- the shape of the PPG waveform varies with the location and manner in which the pulse oximeter is contacted with the body.
- a pulse oximeter measures peripheral oxygen saturation (SpO2). Most often, the device operates in a transmission mode in which two wavelengths of light are passed through a body part to a photodetector. Changes in absorbance at each of the wavelengths are measured, which allow the determination of the absorbance due to pulsing arterial blood, excluding venous blood, skin, bone, muscle, and other tissues. Alternatively, reflectance pulse oximetry can be used.
- a typical pulse oximeter comprises a data processor and a pair of light-emitting diodes (LEDs) facing a photodiode. One LED produces red light having a wavelength of 660 nm, the other infrared light having a wavelength of 940 nm.
- Oxygenated hemoglobin absorbs more infrared light and less red light than hemoglobin.
- the transmission signals fluctuate over time because of changes in the amount of arterial blood present in the tissue caused by the blood pulse associated with each cycle of the heart.
- the ratio of red light measured to infrared light measured is calculated by the processor and is converted by the processor to SpO2 using a lookup table based on the Beer-Lambert law.
- Awad et al. (J. Clinical Monitoring and Computing, 2006, 20:175-184) reports that researchers have been attempting to understand the relationship between central cardiac hemodynamics and the resulting measured peripheral waveforms.
- Awad et al. studied ear pulse oximeter waveforms in order to understand the underlying physiology reflected in these waveforms and to extract information about cardiac performance.
- Multi-linear regression analysis of ear plethysmographic waveform components were used to estimate CO from the ear plethysmograph and it was determined that ear plethysmographic width correlates with CO.
- Awad et al. does not suggest that this correlation allows pulse oximetry or plethysmographic data to be used to calculate SV or CO.
- US 2013/0310669 discloses a method for determining mixed venous oxygen saturation (SvO2) using a photoplethysmography pulse oximeter (PPG PO) device that measures changes in pulmonary circulation using a light source and a light detector applied to the thoracic wall of a patient.
- the light source and the detector are separated by at least 15-20 mm so that the region of illumination overlaps a portion of the pulmonary microcirculation beneath the PPG device.
- the contribution of circulation in the thoracic wall to the PPG signal must be assessed using an additional detector and/or light source that is attached to the thoracic wall less than 8 mm apart.
- the SvO2 value can be used to calculate CO using the Fick method from the values of total oxygen consumption, arterial oxygen content and venous oxygen content.
- One drawback of this method is that the pulmonary microcirculation is surrounded by bone, muscle, and other tissues that make the contribution of circulation in the thoracic wall to the PPG signal difficult to measure reliably.
- Another drawback is that the measured values depend on the accurate placement of the light emitters and sensors, which may be difficult to reproduce for each subsequent measurement.
- U.S. Pat. No. 9,289,133 discloses a method and apparatus for monitoring proportional changes in CO from a blood pressure signal measurement obtained by fingertip PPG.
- a time constant of the arterial tree is defined as the product of the total peripheral resistance (TPR) and a constant arterial compliance and is determined by analyzing long time scale variations of more than one cardiac cycle.
- a value proportional to CO is determined from the ratio of the blood pressure signal to the estimated time constant using Ohm's law.
- An invasive, absolute CO calibrating measurement is required to derive absolute CO values from the proportional CO change values obtained using PPG. This method and apparatus cannot measure values for SV or CO without an invasive CO measurement and therefore requires a hospital setting and skilled medical personnel.
- WO 2010/0274102 A1 discloses a data processing method and associated apparatus and systems that measures SV and CO from pulse oximetry data.
- the pulse oximeter system comprises a data processor configured to perform a method that combines a probabilistic processor and a physiological model of the cardiovascular system in a Dynamic State-Space Model (DSSM) that can remove contaminating noise and artifacts from the pulse oximeter sensor output and measure blood oxygen saturation, HR, SV, aortic pressure and systemic pressures.
- DSSM Dynamic State-Space Model
- This pulse oximeter and associated method provides truly noninvasive measurement of CO and SV suitable for monitoring patients with CHF.
- the DSSM comprises a mathematical model of the cardiovascular system that models the physiological processes which produce the pulses measured by the pulse oximeter.
- the model comprises parameters including aortic pressure, radial pressure, peripheral resistance, aortic impedance, and blood density.
- the present invention fills a need in the art for truly noninvasive apparatus and methods for measuring SV and CO by providing a system and method for measuring and reporting SV and/or CO using plethysmographic (PG) waveform data from a photoplethysomgraphic pulse oximeter (PPG PO).
- PG plethysmographic
- PPG PO photoplethysomgraphic pulse oximeter
- FIG. 1 is a flow chart showing the path of information flow from a biomedical sensor to a data processor and on to an output display according to one embodiment of the invention.
- FIG. 2 is a flow chart showing inputs, outputs, and conceptual division of model parts for a dynamic state-space model (DSSM).
- DSSM dynamic state-space model
- FIG. 3 is a block diagram showing mathematical representations of inputs, output and conceptual divisions of the DSSM shown in FIG. 2 .
- FIG. 4 is a mathematical representation of the process of dual estimation.
- FIG. 5 is a schematic diagram showing the process steps involved in a dual estimation process.
- FIG. 6 is a mathematical representation of the process of joint estimation.
- FIG. 7 is a schematic diagram showing the process steps involved in a joint estimation process.
- FIG. 8 is a flow chart showing the components of a DSSM used for pulse oximetry data processing.
- FIG. 9 is a flow chart showing examples of parameter inputs and outputs for a DSSM used for pulse oximetry data processing.
- FIG. 10 is a flow chart showing the components of a DSSM used for electrocardiography data processing.
- FIG. 11 is a chart showing input sensor data from a fingertip pulse oximeter and processed output data from a data processor configured to process pulse oximetry data.
- FIG. 12 is a chart showing input sensor data from a fingertip pulse oximeter and processed output data from a data processor configured to process pulse oximetry data under a low blood perfusion condition.
- FIG. 13 is a chart showing noisy non-stationary ECG sensor data input and processed heart rate and ECG output for a data processor configured to process ECG sensor data.
- FIG. 14 is a chart showing input ECG sensor data and comparing output data from a data processor according to the present invention with output data generating using a Savitzky-Golay FIR data processing algorithm.
- noninvasively and noninvasive are used to indicate that a method, system, or apparatus does not involve an invasive or minimally invasive procedure such as surgery or catheterization. Additionally, “noninvasively” and “noninvasive” are used to indicate that a method, system, or apparatus does not involve endoscopy, phlebotomy, biopsy, or a procedure that requires intubation or artificial respiration.
- FIG. 1 is a top-level schematic for data processing according to the present invention.
- a biomedical sensor such as pulse oximeter, normally produces a raw analog output signal that is converted to a raw digital output signal.
- the analog to digital conversion may also be accompanied by signal filtering or conditioning.
- Digital signals are received by a data processor configured to process the digital data and produce a processed (or clean) signal comprising an estimated true value for the physiological parameter being measured.
- the processed signal is then displayed, for example, in the form of an electronic, hard copy, audible, visual, and/or tactile output.
- the output may be used, for example, by a user to monitor a patient, by a user for self monitoring, or by a user as biofeedback process.
- the data processor shown in FIG. 1 is configured to receive, as input data, digital signals from one or more biomedical sensors, and enter the data into a dynamic state-space model (DSSM) integrated with a processor engine.
- the integrated DSSM/processor engine produces transformed output data that may correspond to a physiological parameter measured by the biomedical sensor(s), in the form of an estimated true value for the physiological parameter.
- the processor engine may operate in a dual estimation mode (a dual estimation engine) or in a joint estimation mode (a join estimation engine).
- the output may include additional outputs corresponding to physiological parameters not measured by the physiological sensor(s), diagnostic information, and a confidence interval representing the probability that the output estimated value(s) for the physiological parameter(s) is accurate.
- the data transformation process performed by the data processor may be used to remove artifacts from the input data to produce output data having higher accuracy than the input data and/or to extract information from the input data to generate output data estimating the values for physiological parameters that are not otherwise measured or reported using data from the sensor(s).
- the physiological parameter(s) that are not otherwise measured may include SV, and/or CO.
- a mathematical model or a computer model, as used herein, involves the use of state parameters and model parameters in a mathematical representation of physiological processes that give rise to a physiological parameter being measured and processes through which sensor data is detected.
- a mathematical model may include model and/or state parameters that correspond directly to physiological parameters including vital signs such as oxygen saturation of blood (SpO 2 ), heart rate (HR), respiratory rate (RR), and blood pressure (BP) that can be directly measured; physiological parameters not directly measured such as total blood volume (TBV), left-ventricular stroke volume (SV), vasomotor tone (VT), autonomous nervous system (ANS) tone, and cardiac output (CO); and hemoglobin-bound complexes, concentrations of metabolic intermediates, and concentrations of drugs present in one or more tissues or organs.
- vital signs such as oxygen saturation of blood (SpO 2 ), heart rate (HR), respiratory rate (RR), and blood pressure (BP) that can be directly measured
- physiological parameters not directly measured such as total blood volume (TBV), left-ventricular stroke volume (SV), vasomotor tone (VT), autonomous nervous system (ANS) tone, and cardiac output (CO); and hemoglobin-bound complexes, concentrations of metabolic intermediates, and concentrations of drugs present in one or more tissues or organs.
- the fusion of two or more biomedical signals follows the same principle. For instance, if the intent is to measure blood pressure waves and electrocardiogram signals simultaneously, the use of a heart model describing both the electrical and mechanical aspects of the organ should be used. Initially, the model may also accept manual data input as a complement to data from sensors. Nonlimiting examples of manually entered data include food consumption over time vital signs, gender, age, weight, and height.
- Non-physiological models may be included in and/or coupled to the DSSM in cases where non-biomedical signals are measured. For instance, one may use non-biomedical measurements to enhance or complement biomedical measurements.
- a non-limiting example is the use of accelerometer data to enhance motion artifact rejection in biomedical measurements.
- the physiological model is extended to describe both measurements, which may include, in this example, cardiovascular circulation at rest, at different body postures (standing, supine, etc), and in motion.
- FIG. 2 and FIG. 3 show schematics of one embodiment of a dynamic state-space model (DSSM) used in the processing of data according to the present invention.
- the DSSM comprises a process model F that mathematically represents physiological processes involved in generating one or more physiological parameters measured by a biomedical sensor and describes the state of the subject over time in terms of state parameters.
- This mathematical model optimally includes mathematical representations accounting for process noise such as physiologically caused artifacts that may cause the sensor to produce a digital output that does not produce an accurate measurement for the physiological parameter being sensed.
- the DSSM also comprises an observational model H that mathematically represents processes involved in collecting sensor data measured by the biomedical sensor. This mathematical model optimally includes mathematical representations accounting for observation noise produced by the sensor apparatus that may cause the sensor to produce a digital output that does not produce an accurate measurement for a physiological parameter being sensed. Noise terms in the mathematical models are not required to be additive.
- process and observational mathematical models F and H may be conceptualized as separate models, they are normally integrated into a single mathematical model that describes processes that produce a physiological parameter and processes involved in sensing the physiological parameter. That model, in turn, is integrated with a processing engine within an executable program stored in a data processor that is configured to receive digital data from one or more sensors and to output data to a display or other output formats.
- FIG. 3 provides mathematical descriptions of the inputs and outputs corresponding to FIG. 2 .
- values for state parameters preferably in the form of a state vector x k
- Wk input model parameters
- Process noise v k and observation noise n k are also received by the DSSM, which updates the state parameter vector and model parameter vector and produces an output observation vector y k .
- the updated state vector X k+1 , updated model parameters W k+! , and time-specific sensor data are used as input for each calculation for subsequent iterations, or time steps.
- the DSSM is integrated in a dual estimation processing engine or a joint estimation processing engine.
- the most favored embodiment makes use of a DSSM built into a Sigma point Kalman filter (SPKF) or Sequential Monte Carlo (SMC) processing engine.
- Sigma point Kalman filter (SPKF) refers to the collective name used for derivativeless Kalman filters that employ the deterministic sampling based sigma point approach to calculate approximations of the optimal terms of the Gaussian approximate linear Bayesian update rule, including unscented, central difference, square-root unscented, and square-root central difference Kalman filters.
- SMC and SPKF processing engines operate on a general nonlinear DSSM having the form:
- a hidden system state, x k propagates over time index, k, according to the system model, f.
- the process noise is v k ⁇ 1
- W is the vector of model parameters.
- Observations, y k about the hidden state are given by the observation model h and n k is the measurement noise.
- W is fixed, only state estimation is required and either SMC or SPKF can be used to estimate the hidden states.
- Unsupervised machine learning sometimes referred to as system identification or parameter estimation, involves determining the nonlinear mapping:
- the nonlinear map g( ⁇ ) is parameterized by the model parameter vector W.
- the nonlinear map may be a feed-forward neural network, recurrent neural network, expectation maximization algorithm, or enhanced Kalman filter algorithm. Learning corresponds to estimating W in some optimal fashion.
- SPKF or SMC is used for updating parameter estimates. One way to accomplish this is to write a new state-space representation
- w k correspond to a stationary process with identity state transition matrix, driven by process noise r k .
- the desired output d k corresponds to a nonlinear observation on w k .
- the state and parameter estimation steps may be coupled in an iterative dual-estimation mode as shown in FIGS. 4 and 5 .
- This formulation for a state estimator operates on an adaptive DSSM.
- states x k and parameters W are estimated sequentially inside a loop.
- states x k and parameters W are estimated sequentially inside a loop.
- Parameter estimates are passed from the previous iteration to state estimation for the current iteration.
- SMC and SPKF methods including the sigma-point, Gaussian-sum, and square-root forms. The particular choice may be influenced by the application.
- the current estimate of the parameters W k is used in the state estimator as a given (known) input, and likewise the current estimate of the state x k is used in the parameter estimator. This results in a step-wise stochastic optimization within the combined state-parameter space.
- the flow chart shown FIG. 5 provides a summary of the steps involved in dual estimation process.
- Initial probability distributions for state and model parameters are provided to the DSSM to produce an initial probability distribution function (first PDF or prior PDF) representing the initial state.
- first PDF or prior PDF data for a time t 1 from a sensor (new measurement) and the initial PDF are combined using a Bayesian statistical process to generate a second, posterior PDF that represents the state at the time of the measurement for the first sensor data.
- Expectation values for the second PDF are calculated, which may represent the most likely true value.
- Expectation values may also represent, for instance, the confidence interval or any statistical measure of uncertainty associated with the value.
- updated state parameters for time t 1 are combined with sensor data for time t 1 to update the model parameters for time t 2 in the DSSM by the process shown in FIG. 4 .
- the expectation values are also fed into the DSSM as the state, in the form of a vector of state parameters (new PDF) as shown in FIG. 4 .
- new PDF new PDF
- the process is repeated with timed data for time t 2 to produce updated parameters for time t 2 and so forth.
- the time interval between time steps is usually constant such that time points may be described as t, t+n, t+2n, etc. If the time interval is not constant, then the time may be described using two or more time intervals as t, t+n, t+n+m, etc.
- the state and parameter estimation steps may also be performed in a simultaneous joint-estimation mode as shown in FIGS. 6 and 7 .
- the calculated variables for the state parameters and model parameters of the physiological model are concatenated into a single higher-dimensional joint state vector:
- x k are the state parameters and w k the model parameters.
- the joint state space is used to produce simultaneous estimates of the states and parameters.
- the flow chart shown FIG. 7 provides a summary of the steps involved in dual estimation process.
- the process is similar to that shown for dual estimation in FIG. 5 , with the exception that model and state parameters are not separated into two separate vectors, but are represented together in a single vector.
- the process is initiated by entering a vector representing initial state and model parameter value distributions into the DSSM and producing an initial first PDF.
- the first PDF is combined with sensor data (new measurement) for time t 1 in a Bayesian statistical process to generate a second, posterior PDF that represents the state and model parameters at time t 1 .
- Expectation values for the second PDF are calculated, which may represent the most likely true value.
- Expectation values may also represent, for instance, the confidence interval or any statistical measure of uncertainty associated with the value.
- updated state and model parameters for time t 1 are entered into the DSSM by the process shown in FIG. 6 .
- the process is repeated with timed data for time t 2 to produce updated parameters for time t 2 and so forth.
- both state and model parameters are concatenated into a single vector that is transformed by the DSSM.
- no machine learning step is necessary in order to update model parameters.
- Joint estimation may be performed using a sequential Monte Carlo method or sigma-point Kalman method. These may take the form of unscented, central difference, square-root unscented, and square-root central difference forms. The optimal method will depend on the particular application.
- SMC methods estimate the probability distributions of all the model unknowns by propagating a large number of samples called probability particles in accordance with the system models (typically nonlinear, non-Gaussian, non-stationary) and the rules of probability. Artifacts are equivalent to noise with short-lived probability distributions, also called non-stationary distributions. The number of simulated particles scales linearly with computational power, with ⁇ 100 particles being reasonable for real time processing with presently available processors.
- the system model describes pertinent physiology and the processor engine uses the system model as a “template” from which to calculate, using Bayesian statistics, posterior probability distribution functions (processed data). From this, the expectation values (e.g. the mean) and confidence intervals can be estimated FIG. 7 .
- the combination of SMC with Bayesian Statistics to calculate posterior probability distribution functions is often referred to as a Particle Filter.
- SMC process nonlinear and non-Gaussian problems by discretizing the posterior into weighted samples, or probability particles, and evolving them using Monte Carlo simulation.
- Monte Carlo simulation uses weighted particles to map integrals to discrete sums:
- the optimal Bayesian solution can be outlined by the following recursive algorithm. Suppose the required PDF p(x k ⁇ 1
- x k ⁇ 1 ) is defined by the system equation (1) and the known statistics of v k ⁇ 1 .
- a measurement y k becomes available, and this may be used to update the prior (updated stage) via Bayes' rule:
- Equation 5 The known distribution is introduced into Equation 5 to yield:
- Importance sampling is made sequential by reiterating the Markov 1 st order assumption, resulting in the assumption that the current state is not dependent on future observations:
- SIS Sequential Importance Sampling
- SIR Sampling-Importance Resampling
- a pseudo-code for a generic SMC (also called bootstrap filter or condensation algorithm) can be written as:
- the posterior statistics are calculated using a procedure resembling the linear Kalman filter.
- the time-update equations are:
- X h ⁇ h - 1 x f ⁇ ( X h - 1 x , X h - 1 v , u h - 1 ) ( 25 )
- the mathematical structure for sequential Monte Carlo and SPKF represent two examples of a family of probabilistic inference methods exploiting Monte Carlo simulation and the sigma point transform, respectively, in conjunction with a Bayesian statistical process.
- SPKF are generally inferior to SMC but are computationally cheaper. Like SMC, SPKF evolve the state using the full nonlinear DSSM, but represent probability distributions using a sigma-point set. This is a deterministic step that replaces the stochastic Monte Carlo step in the SMC. As a result, SPKF lose accuracy when posterior distributions depart heavily from the Gaussian form, such as with bimodal or heavily-tailed distributions, or with strong nonstationary distributions such as those caused by motion artifacts in pulse oximeters. For these cases SMC are more suitable.
- SPKF yields higher-order accuracy than the extended Kalman filter (EKF) and its related variants with equal algorithm complexity, O (L 2 ).
- EKF returns 2 nd order accuracy for nonlinear and non-Gaussian problems, and 3 rd order for Gaussian problems.
- EKF has only 1 st order accuracy for nonlinear problems.
- Both EKF and SPKF approximate state distributions with Gaussian random variables (GRV). However, the EKF propagates the GRV using a single measure (usually the mean) and the 1 st order Taylor expansion of the nonlinear system.
- the SPKF decomposes the GRV into distribution moments (sigma points) and propagates those using the unmodified nonlinear system. SPKF implementation is simpler than EKF since it is derivativeless. That is, it uses the unmodified DSSM form, and therefore does not require lengthy Jacobian derivations.
- the data processing method is also capable of prediction because the method can operate faster than real time measurements.
- the measurement PDF obtained either from SPKF or SMC, embodies all available statistical information up to that point in time. It is therefore possible to march the system model forwards in time, for instance, using the same sequential Monte Carlo method, to obtain deterministic or stochastic simulations of future signal trajectories. In this way, the future health status (physiological state) of a patient can be predicted with attached probabilities indicating the confidence of each prediction.
- the data processing method may benefit from a noise adaptation method if timed sensor data contains noise and/or artifact that changes its spectral qualities over time. That is, if timed sensor data has a non-stationary probability distribution function.
- a known algorithm such as the Robins-Monro or Annealing methods may be added to the data processing method in order to adapt the probability distribution functions of noise terms (stochastic terms) in the DSSM according to changing noise and artifact present in sensor data.
- the output may include estimates of the true measured signals (i.e. processed data), and estimates of values for one or more physiological parameters measured by one or more sensors from which data was received, and estimates of values for one or more physiological parameters not measured by the sensors from which data was received (data extraction).
- a state parameter estimate is the processed data from the physiological sensor. Both noise and artifacts can be attenuated or rejected even though they may have very distinct probability distribution functions and may mimic the real signal.
- a model parameter estimate may be also used to produce a physiological parameter. For example, an estimate of total blood volume may be used to diagnose hemorrhage or hypovolemia. An estimate of tissue oxygen saturation may indicate poor tissue perfusion and/or hypoxia. Estimates of glucose uptake in several tissues may differentiate between diabetes mellitus types and severities; and estimates of carotid artery radius may be indicative of carotid artery stenosis.
- FIG. 8 shows the components of a DSSM suitable for processing data from a pulse oximeter, including components required to describe processes occurring in a subject.
- FIG. 9 illustrates the DSSM broken down into process and observation models, and including input and output variables.
- Heart rate (HR), stroke volume (SV) and whole-blood oxygen saturation (SpO 2 ) are estimated from input noisy red and infrared intensity ratios (R).
- Radial (Pw) and aortic (Pao) pressures are also available as state estimates.
- the DSSM comprises the following function to represent CO:
- the CO function pumps blood into a Windkessel 3-element model of the vascular system including two state parameters: aortic pressure, Pao, and radial (Windkessel) pressure, P W :
- Rp and Zo are the peripheral resistance and characteristic aortic impedance, respectively. The sum of these two terms is the total peripheral resistance due to viscous (Poiseuille-like) dissipation:
- the elastic component due to vessel compliance is a nonlinear function including thoracic aortic cross-sectional area, A:
- a ⁇ ( P CO ) A max ⁇ [ 1 2 + 1 ⁇ ⁇ arctan ⁇ ( P CO - P 0 P 1 ) ] ( 35 )
- the time-varying Windkessel compliance, Cw, and the aortic compliance per unit length, Cl, are:
- the peripheral resistance is defined as the ratio of average pressure to average flow.
- R p P set ( HR ⁇ SV ) / 60 ( 40 )
- Va is the arterial blood volume, which changes as the cross-sectional area of illuminated blood vessels, ⁇ A w , changes as:
- the optical absorption cross sections for red blood cells containing totally oxygenated (HbO 2 ) and totally deoxygenated (Hb) hemoglobin are ⁇ a 100% and ⁇ a 0% , respectively.
- the function K( ⁇ ,d,r) contains, along with the scattering coefficient, the wavelength, sensor geometry and oxygen saturation dependencies that alter the effective optical pathlengths:
- the attenuation coefficient ⁇ is:
- ⁇ a and ⁇ s are whole-tissue absorption and scattering coefficients, respectively, which are calculated from Mie Theory.
- Red and infrared K values as a function of SpO 2 may be represented by two linear fits:
- R r,k+1 c ⁇ ′ z,r K r ⁇ a,r art ⁇ A w +R r,k +v r (49)
- R ir,k+1 c ⁇ ′ z,ir K ir ⁇ a,ir art ⁇ A w +R ir,k +v ir (50)
- v are Gaussian-distributed process noises intended to capture the baseline wander of the two channels.
- the constant c subsumes all factors common to both wavelengths and is treated as a calibration constant.
- the observation model adds Gaussian-distributed noises, n, to R r and R ir :
- FIG. 11 plots estimates for a 15 s stretch of data.
- Photoplethysmographic (PPG) waveforms (A) were used to extract heart rate (B), SV (C), CO (D), blood oxygen saturation (E), and aortic and systemic (radial) pressure waveforms (F). Results of processing pulse oximetry at low blood perfusion are shown in FIG. 12 .
- Low signal-to-noise PPG waveforms (A) were used to extract heart rate (B), SV (C), blood oxygen saturation (D), and aortic and systemic (radial) pressure waveforms (E).
- FIG. 10 is a schematic of a DSSM suitable for processing electrocardiograph data, including components required to describe the processes occurring in a subject.
- the combination of SPKF or SMC in state, joint or dual estimation modes can be used to filter electrocardiography (ECG) data.
- ECG electrocardiography
- Any physiology model adequately describing the ECG signal can be used, as well as any model of noise and artifact sources interfering or contaminating the signal.
- One non-limiting example of such a model is the ECG signal generator proposed by McSharry (IEEE Transactions on Biomedical Engineering, 2003. 50(3):289-294).
- this model uses a sum of Gaussians with amplitude, center and standard deviation, respectively, for each wave (P, Q, R, S, T ⁇ , T+15) in an ECG.
- the observation model comprises the state plus additive Gaussian noise, but more realistic pink noise or any other noise distributions can be used.
- FIG. 13 shows the results of processing a noisy non-stationary ECG signal.
- Heart rate oscillations representative of normal respiratory sinus arrhythmia are present in the ECG.
- the processor accomplishes accurate, simultaneous estimation of the true ECG signal and heart rate that follows closely the true values.
- the performance of the processor in a noise and artifact-corrupted signal is shown in FIG. 14 .
Abstract
A method for non-invasively measuring cardiac output, stroke volume, or both comprises collecting plethysmographic waveform data of a patient, providing the plethysmographic waveform to a data processor, and calculating measured values for SV, CO, or both. Software of the data processor comprises a mathematical model of the cardiovascular system integrated in a dynamic state space model (DSSM).
Description
- This application is a divisional application of U.S. Ser. No. 15/225,803 filed Aug. 1, 2016, which is is a continuation-in-part of U.S. Ser. No. 12/640,278 filed Dec. 17, 2009, which is non-provisional of U.S. 61/171,802 filed Apr. 22, 2009. This application is additionally related to the following US patent applications:
- Ser. No. 12/796,512 filed Jun. 8, 2010, issued as U.S. Pat. No. 9,060,722;
Ser. No. 13/096,845 filed Apr. 28, 2011, issued as U.S. Pat. No. 9,173,574;
Ser. No. 13/096,876 filed Apr. 28, 2011, issued as U.S. Pat. No. 9,275,171;
Ser. No. 13/096,904 filed Apr. 28, 2011, issued as U.S. Pat. No. 9,649,036
Ser. No. 13/181,027 filed Jul. 12, 2011, issued as U.S. Pat. No. 8,494,829;
Ser. No. 13/181,140 filed Jul. 12, 2011, issued as U.S. Pat. No. 9,451,886; and Ser. No.
Ser. No. 13/181,247 filed Jul. 12, 2011. - The U.S. Government may have certain rights to this invention pursuant to Contract Number IIP-0839734 awarded by the National Science Foundation.
- The present invention relates generally to apparatus, systems, and methods for noninvasively measuring cardiac output (CO) and left ventricular stroke volume (SV). More specifically, the invention includes systems and methods that measure CO and/or SV using plethysmographic waveform data collected by a pulse oximeter.
- There has been a long felt need in the medical arts for a noninvasive way to measure left ventricular stroke volume (SV) and cardiac output (CO). SV is the volume of blood pumped by the left ventricle of the heart with a single heart cycle. CO is the product of SV and hear rate (HR). SV and CO are important physiological parameters for a number of medical conditions, including congestive heart failure (CHF).
- Patients being treated for heart failure are normally medicated with drugs that regulate diuresis and heart muscle function. An important aim of drug therapy is to maintain a CO that is sufficient to perfuse tissues with oxygenated blood. It is advantageous to use the lowest possible doses of drugs to manage CHF because the drugs used produce unwanted side effects. To optimally manage the dosage and selection of drugs in the treatment of CHF, one must monitor CO to assess the efficacy of the drugs and dosages being administered and/or to monitor patient compliance.
- Currently, CO is measured using invasive techniques such as the Fick method, the Thermodilution method, and implantable microelectromechanical devices (MEMs), also called CardioMEMs. The Fick method involves the measurement of oxygen consumption and computing the arteriovenous difference using samples of arterial blood and mixed venous blood from the pulmonary artery. The thermodilution method measures the rate at which cold saline solution is diluted in the blood. Both of these methods are performed in a hospital setting because they require the placement of a catheter in the pulmonary artery. Additionally, the use of a pulmonary artery catheter use may also increase morbidity in critically ill patients. CardioMEMs are surgically implanted into patients in a hospital setting and, once implanted, provide measurements of SV and CO that can be used to monitor patients. The implantation of the cardioMEMs device into the pulmonary artery, however, is expensive and is performed in a hospital setting and involves the risks associated with heart catheterization.
- Less invasive techniques for measuring SV and CO include esophageal Doppler and transesophageal echocardiography. Esophageal Doppler measures blood flow velocity in the descending thoracic aorta using a flexible ultrasound probe that is inserted into the esophagus. The blood flow velocity is combined with an estimate of the cross-sectional area of the aorta estimated from the patient's age, height, and weight to calculate SV and CO. This technique requires someone with technical skill to insert an esophageal Doppler monitor, which must be properly aligned with respect to the thoracic aorta to provide accurate measurements. Transesophageal echocardiography involves measuring SV using flow velocity calculated from the area under the measured Doppler velocity waveform at the pulmonary artery, the mitral valve, or the aortic valve. This technique requires a highly trained operator to position and place the esophageal Doppler monitor. These procedures are not truly noninvasive because accessing the esophagus is perceived by patients as invasively uncomfortable.
- Truly noninvasive methods, apparatus, and systems are needed that can measure SV and CO, preferably in the homes of patients without the need for skilled caregivers. Pulse oximetry has been investigated for decades as a possible tool for the noninvasive measurement of SV and CO. A pulse oximeter (PO) is a device that obtains photoplethysmography (PPG) data, which measures changes in blood volume within a tissue caused by the pulse of blood pressure through the vasculature in the tissue. The blood volume change is detected by measuring the amount of light transmitted or reflected to a sensor from a light source used to illuminate the skin. The shape of the PPG waveform varies with the location and manner in which the pulse oximeter is contacted with the body. In addition to PPG data, a pulse oximeter measures peripheral oxygen saturation (SpO2). Most often, the device operates in a transmission mode in which two wavelengths of light are passed through a body part to a photodetector. Changes in absorbance at each of the wavelengths are measured, which allow the determination of the absorbance due to pulsing arterial blood, excluding venous blood, skin, bone, muscle, and other tissues. Alternatively, reflectance pulse oximetry can be used. A typical pulse oximeter comprises a data processor and a pair of light-emitting diodes (LEDs) facing a photodiode. One LED produces red light having a wavelength of 660 nm, the other infrared light having a wavelength of 940 nm. Oxygenated hemoglobin absorbs more infrared light and less red light than hemoglobin. The transmission signals fluctuate over time because of changes in the amount of arterial blood present in the tissue caused by the blood pulse associated with each cycle of the heart. The ratio of red light measured to infrared light measured is calculated by the processor and is converted by the processor to SpO2 using a lookup table based on the Beer-Lambert law.
- Awad et al. (J. Clinical Monitoring and Computing, 2006, 20:175-184) reports that researchers have been attempting to understand the relationship between central cardiac hemodynamics and the resulting measured peripheral waveforms. Awad et al. studied ear pulse oximeter waveforms in order to understand the underlying physiology reflected in these waveforms and to extract information about cardiac performance. Multi-linear regression analysis of ear plethysmographic waveform components were used to estimate CO from the ear plethysmograph and it was determined that ear plethysmographic width correlates with CO. Awad et al. does not suggest that this correlation allows pulse oximetry or plethysmographic data to be used to calculate SV or CO.
- Natalini et al. (Anesth. Analg. 2006, 103:1478-1484) reports the use of pulse oximetry to predict which hypotensive patients are likely to respond positively to increasing blood volume. Arterial blood pressure changes during mechanical ventilation are reported as accurately predicting fluid responsiveness. Photoplethysmographic (PPG) waveform variations measured by pulse oximetry showed a correlation with measured pulse pressure variation values associated with fluid responsiveness but neither SV nor CO were calculated. Natilini does not indicate that SV or CO can be calculated using PPG data.
- US 2013/0310669 discloses a method for determining mixed venous oxygen saturation (SvO2) using a photoplethysmography pulse oximeter (PPG PO) device that measures changes in pulmonary circulation using a light source and a light detector applied to the thoracic wall of a patient. The light source and the detector are separated by at least 15-20 mm so that the region of illumination overlaps a portion of the pulmonary microcirculation beneath the PPG device. The contribution of circulation in the thoracic wall to the PPG signal must be assessed using an additional detector and/or light source that is attached to the thoracic wall less than 8 mm apart. The SvO2 value can be used to calculate CO using the Fick method from the values of total oxygen consumption, arterial oxygen content and venous oxygen content. One drawback of this method is that the pulmonary microcirculation is surrounded by bone, muscle, and other tissues that make the contribution of circulation in the thoracic wall to the PPG signal difficult to measure reliably. Another drawback is that the measured values depend on the accurate placement of the light emitters and sensors, which may be difficult to reproduce for each subsequent measurement.
- U.S. Pat. No. 9,289,133 discloses a method and apparatus for monitoring proportional changes in CO from a blood pressure signal measurement obtained by fingertip PPG. A time constant of the arterial tree is defined as the product of the total peripheral resistance (TPR) and a constant arterial compliance and is determined by analyzing long time scale variations of more than one cardiac cycle. A value proportional to CO is determined from the ratio of the blood pressure signal to the estimated time constant using Ohm's law. An invasive, absolute CO calibrating measurement is required to derive absolute CO values from the proportional CO change values obtained using PPG. This method and apparatus cannot measure values for SV or CO without an invasive CO measurement and therefore requires a hospital setting and skilled medical personnel.
- WO 2010/0274102 A1 discloses a data processing method and associated apparatus and systems that measures SV and CO from pulse oximetry data. The pulse oximeter system comprises a data processor configured to perform a method that combines a probabilistic processor and a physiological model of the cardiovascular system in a Dynamic State-Space Model (DSSM) that can remove contaminating noise and artifacts from the pulse oximeter sensor output and measure blood oxygen saturation, HR, SV, aortic pressure and systemic pressures. This pulse oximeter and associated method provides truly noninvasive measurement of CO and SV suitable for monitoring patients with CHF. The DSSM comprises a mathematical model of the cardiovascular system that models the physiological processes which produce the pulses measured by the pulse oximeter. In one embodiment, the model comprises parameters including aortic pressure, radial pressure, peripheral resistance, aortic impedance, and blood density.
- The present invention fills a need in the art for truly noninvasive apparatus and methods for measuring SV and CO by providing a system and method for measuring and reporting SV and/or CO using plethysmographic (PG) waveform data from a photoplethysomgraphic pulse oximeter (PPG PO).
-
FIG. 1 is a flow chart showing the path of information flow from a biomedical sensor to a data processor and on to an output display according to one embodiment of the invention. -
FIG. 2 is a flow chart showing inputs, outputs, and conceptual division of model parts for a dynamic state-space model (DSSM). -
FIG. 3 is a block diagram showing mathematical representations of inputs, output and conceptual divisions of the DSSM shown inFIG. 2 . -
FIG. 4 is a mathematical representation of the process of dual estimation. -
FIG. 5 is a schematic diagram showing the process steps involved in a dual estimation process. -
FIG. 6 is a mathematical representation of the process of joint estimation. -
FIG. 7 is a schematic diagram showing the process steps involved in a joint estimation process. -
FIG. 8 is a flow chart showing the components of a DSSM used for pulse oximetry data processing. -
FIG. 9 is a flow chart showing examples of parameter inputs and outputs for a DSSM used for pulse oximetry data processing. -
FIG. 10 is a flow chart showing the components of a DSSM used for electrocardiography data processing. -
FIG. 11 is a chart showing input sensor data from a fingertip pulse oximeter and processed output data from a data processor configured to process pulse oximetry data. -
FIG. 12 is a chart showing input sensor data from a fingertip pulse oximeter and processed output data from a data processor configured to process pulse oximetry data under a low blood perfusion condition. -
FIG. 13 is a chart showing noisy non-stationary ECG sensor data input and processed heart rate and ECG output for a data processor configured to process ECG sensor data. -
FIG. 14 is a chart showing input ECG sensor data and comparing output data from a data processor according to the present invention with output data generating using a Savitzky-Golay FIR data processing algorithm. - As used herein, “noninvasively” and “noninvasive” are used to indicate that a method, system, or apparatus does not involve an invasive or minimally invasive procedure such as surgery or catheterization. Additionally, “noninvasively” and “noninvasive” are used to indicate that a method, system, or apparatus does not involve endoscopy, phlebotomy, biopsy, or a procedure that requires intubation or artificial respiration.
-
FIG. 1 is a top-level schematic for data processing according to the present invention. A biomedical sensor, such as pulse oximeter, normally produces a raw analog output signal that is converted to a raw digital output signal. The analog to digital conversion may also be accompanied by signal filtering or conditioning. Digital signals are received by a data processor configured to process the digital data and produce a processed (or clean) signal comprising an estimated true value for the physiological parameter being measured. The processed signal is then displayed, for example, in the form of an electronic, hard copy, audible, visual, and/or tactile output. The output may be used, for example, by a user to monitor a patient, by a user for self monitoring, or by a user as biofeedback process. - The data processor shown in
FIG. 1 is configured to receive, as input data, digital signals from one or more biomedical sensors, and enter the data into a dynamic state-space model (DSSM) integrated with a processor engine. The integrated DSSM/processor engine produces transformed output data that may correspond to a physiological parameter measured by the biomedical sensor(s), in the form of an estimated true value for the physiological parameter. The processor engine may operate in a dual estimation mode (a dual estimation engine) or in a joint estimation mode (a join estimation engine). The output may include additional outputs corresponding to physiological parameters not measured by the physiological sensor(s), diagnostic information, and a confidence interval representing the probability that the output estimated value(s) for the physiological parameter(s) is accurate. The data transformation process performed by the data processor may be used to remove artifacts from the input data to produce output data having higher accuracy than the input data and/or to extract information from the input data to generate output data estimating the values for physiological parameters that are not otherwise measured or reported using data from the sensor(s). In the case of a pulse oximeter, the physiological parameter(s) that are not otherwise measured may include SV, and/or CO. - A mathematical model or a computer model, as used herein, involves the use of state parameters and model parameters in a mathematical representation of physiological processes that give rise to a physiological parameter being measured and processes through which sensor data is detected.
- A mathematical model may include model and/or state parameters that correspond directly to physiological parameters including vital signs such as oxygen saturation of blood (SpO2), heart rate (HR), respiratory rate (RR), and blood pressure (BP) that can be directly measured; physiological parameters not directly measured such as total blood volume (TBV), left-ventricular stroke volume (SV), vasomotor tone (VT), autonomous nervous system (ANS) tone, and cardiac output (CO); and hemoglobin-bound complexes, concentrations of metabolic intermediates, and concentrations of drugs present in one or more tissues or organs.
- While the scope of a mathematical model used in the context of the present invention cannot possibly encompass every single process of human physiology, it should have the capacity to interpret the measured observable(s). For instance, if the intent is to process electrocardiography (ECG) signals, a model describing the generation and propagation of electrical impulses in the heart should be included.
- The fusion of two or more biomedical signals follows the same principle. For instance, if the intent is to measure blood pressure waves and electrocardiogram signals simultaneously, the use of a heart model describing both the electrical and mechanical aspects of the organ should be used. Initially, the model may also accept manual data input as a complement to data from sensors. Nonlimiting examples of manually entered data include food consumption over time vital signs, gender, age, weight, and height.
- Non-physiological models may be included in and/or coupled to the DSSM in cases where non-biomedical signals are measured. For instance, one may use non-biomedical measurements to enhance or complement biomedical measurements. A non-limiting example is the use of accelerometer data to enhance motion artifact rejection in biomedical measurements. In order to accomplish this, the physiological model is extended to describe both measurements, which may include, in this example, cardiovascular circulation at rest, at different body postures (standing, supine, etc), and in motion.
-
FIG. 2 andFIG. 3 show schematics of one embodiment of a dynamic state-space model (DSSM) used in the processing of data according to the present invention. The DSSM comprises a process model F that mathematically represents physiological processes involved in generating one or more physiological parameters measured by a biomedical sensor and describes the state of the subject over time in terms of state parameters. This mathematical model optimally includes mathematical representations accounting for process noise such as physiologically caused artifacts that may cause the sensor to produce a digital output that does not produce an accurate measurement for the physiological parameter being sensed. The DSSM also comprises an observational model H that mathematically represents processes involved in collecting sensor data measured by the biomedical sensor. This mathematical model optimally includes mathematical representations accounting for observation noise produced by the sensor apparatus that may cause the sensor to produce a digital output that does not produce an accurate measurement for a physiological parameter being sensed. Noise terms in the mathematical models are not required to be additive. - While the process and observational mathematical models F and H may be conceptualized as separate models, they are normally integrated into a single mathematical model that describes processes that produce a physiological parameter and processes involved in sensing the physiological parameter. That model, in turn, is integrated with a processing engine within an executable program stored in a data processor that is configured to receive digital data from one or more sensors and to output data to a display or other output formats.
-
FIG. 3 provides mathematical descriptions of the inputs and outputs corresponding toFIG. 2 . Initially, values for state parameters, preferably in the form of a state vector xk, are received by the DSSM together with input model parameters Wk. Process noise vk and observation noise nk are also received by the DSSM, which updates the state parameter vector and model parameter vector and produces an output observation vector yk. Once the model is initialized, the updated state vector Xk+1, updated model parameters Wk+!, and time-specific sensor data are used as input for each calculation for subsequent iterations, or time steps. - The DSSM is integrated in a dual estimation processing engine or a joint estimation processing engine. The most favored embodiment makes use of a DSSM built into a Sigma point Kalman filter (SPKF) or Sequential Monte Carlo (SMC) processing engine. Sigma point Kalman filter (SPKF), as used herein, refers to the collective name used for derivativeless Kalman filters that employ the deterministic sampling based sigma point approach to calculate approximations of the optimal terms of the Gaussian approximate linear Bayesian update rule, including unscented, central difference, square-root unscented, and square-root central difference Kalman filters. SMC and SPKF processing engines operate on a general nonlinear DSSM having the form:
-
x k =f(x k−1 ,v k−1 ;W) (1) -
y k =h(x k ,n k ;W) (2) - A hidden system state, xk, propagates over time index, k, according to the system model, f. The process noise is vk−1, and W is the vector of model parameters. Observations, yk, about the hidden state are given by the observation model h and nk is the measurement noise. When W is fixed, only state estimation is required and either SMC or SPKF can be used to estimate the hidden states.
- Unsupervised machine learning, sometimes referred to as system identification or parameter estimation, involves determining the nonlinear mapping:
-
y k =g(x k ;w k) (3) - where xk is the input, yk is the output, and the nonlinear map g(·) is parameterized by the model parameter vector W. The nonlinear map, for example, may be a feed-forward neural network, recurrent neural network, expectation maximization algorithm, or enhanced Kalman filter algorithm. Learning corresponds to estimating W in some optimal fashion. In the preferred embodiment, SPKF or SMC is used for updating parameter estimates. One way to accomplish this is to write a new state-space representation
-
w k+1 =w k +r k (4) -
d k =g(x k ;w k)+e k (5) - where wk correspond to a stationary process with identity state transition matrix, driven by process noise rk. The desired output dk corresponds to a nonlinear observation on wk.
- The state and parameter estimation steps may be coupled in an iterative dual-estimation mode as shown in
FIGS. 4 and 5 . This formulation for a state estimator operates on an adaptive DSSM. In the dual estimation process, states xk and parameters W are estimated sequentially inside a loop. When used in a data processor for a pulse oximeter, the current state xk from pulse oximeter sensor input yk. States xk, and parameters W are estimated sequentially inside a loop. Parameter estimates are passed from the previous iteration to state estimation for the current iteration. Several different implementations or variants of the SMC and SPKF methods exist, including the sigma-point, Gaussian-sum, and square-root forms. The particular choice may be influenced by the application. - The current estimate of the parameters Wk is used in the state estimator as a given (known) input, and likewise the current estimate of the state xk is used in the parameter estimator. This results in a step-wise stochastic optimization within the combined state-parameter space.
- The flow chart shown
FIG. 5 provides a summary of the steps involved in dual estimation process. Initial probability distributions for state and model parameters are provided to the DSSM to produce an initial probability distribution function (first PDF or prior PDF) representing the initial state. Data for a time t1 from a sensor (new measurement) and the initial PDF are combined using a Bayesian statistical process to generate a second, posterior PDF that represents the state at the time of the measurement for the first sensor data. Expectation values for the second PDF are calculated, which may represent the most likely true value. Expectation values may also represent, for instance, the confidence interval or any statistical measure of uncertainty associated with the value. Based upon the expectation values, usually but not necessarily the values for state parameters having the highest probability of being correct, updated state parameters for time t1 are combined with sensor data for time t1 to update the model parameters for time t2 in the DSSM by the process shown inFIG. 4 . The expectation values are also fed into the DSSM as the state, in the form of a vector of state parameters (new PDF) as shown inFIG. 4 . Once the state parameters and model parameters for the DSSM are updated to time t1, the process is repeated with timed data for time t2 to produce updated parameters for time t2 and so forth. The time interval between time steps is usually constant such that time points may be described as t, t+n, t+2n, etc. If the time interval is not constant, then the time may be described using two or more time intervals as t, t+n, t+n+m, etc. - The state and parameter estimation steps may also be performed in a simultaneous joint-estimation mode as shown in
FIGS. 6 and 7 . The calculated variables for the state parameters and model parameters of the physiological model are concatenated into a single higher-dimensional joint state vector: -
X=[x k T w k T]T (6) - where xk are the state parameters and wk the model parameters. The joint state space is used to produce simultaneous estimates of the states and parameters.
- The flow chart shown
FIG. 7 provides a summary of the steps involved in dual estimation process. The process is similar to that shown for dual estimation inFIG. 5 , with the exception that model and state parameters are not separated into two separate vectors, but are represented together in a single vector. The process is initiated by entering a vector representing initial state and model parameter value distributions into the DSSM and producing an initial first PDF. The first PDF is combined with sensor data (new measurement) for time t1 in a Bayesian statistical process to generate a second, posterior PDF that represents the state and model parameters at time t1. Expectation values for the second PDF are calculated, which may represent the most likely true value. Expectation values may also represent, for instance, the confidence interval or any statistical measure of uncertainty associated with the value. Based upon the expectation values, updated state and model parameters for time t1 are entered into the DSSM by the process shown inFIG. 6 . Once the state parameters and model parameters for the DSSM are updated to time t1, the process is repeated with timed data for time t2 to produce updated parameters for time t2 and so forth. - Compared to dual estimation, both state and model parameters are concatenated into a single vector that is transformed by the DSSM. Hence, no machine learning step is necessary in order to update model parameters. Joint estimation may be performed using a sequential Monte Carlo method or sigma-point Kalman method. These may take the form of unscented, central difference, square-root unscented, and square-root central difference forms. The optimal method will depend on the particular application.
- SMC methods estimate the probability distributions of all the model unknowns by propagating a large number of samples called probability particles in accordance with the system models (typically nonlinear, non-Gaussian, non-stationary) and the rules of probability. Artifacts are equivalent to noise with short-lived probability distributions, also called non-stationary distributions. The number of simulated particles scales linearly with computational power, with ≥100 particles being reasonable for real time processing with presently available processors. The system model describes pertinent physiology and the processor engine uses the system model as a “template” from which to calculate, using Bayesian statistics, posterior probability distribution functions (processed data). From this, the expectation values (e.g. the mean) and confidence intervals can be estimated
FIG. 7 . The combination of SMC with Bayesian Statistics to calculate posterior probability distribution functions is often referred to as a Particle Filter. - SMC process nonlinear and non-Gaussian problems by discretizing the posterior into weighted samples, or probability particles, and evolving them using Monte Carlo simulation. For discretization, Monte Carlo simulation uses weighted particles to map integrals to discrete sums:
-
- where the random samples {x(i); i=1, 2, . . . , N}, are drawn from p(xk|y1:k) and δ(·) is the Dirac delta function. Expectations of the form
-
E[g(x k)]=∫g(x k)p(x k |y 1:k)dx k (8) - can be approximated by the estimate:
-
- if the distribution has finite support. As N approaches infinity, the estimate converges to the true expectation.
- The optimal Bayesian solution can be outlined by the following recursive algorithm. Suppose the required PDF p(xk−1|y1:k−1) at time k−1 is available. In the prediction stage, the prior PDF at time k is obtained using the DSSM via the Chapman-Kolmogorov equation:
-
p(x k |y 1:k−1)=∫p(x k |y k−1)p(x k−1 |y 1:k−1)dx k−1 (10) - The DSSM model describing the state evolution p(xk|xk−1) is defined by the system equation (1) and the known statistics of vk−1. At time step k a measurement yk becomes available, and this may be used to update the prior (updated stage) via Bayes' rule:
-
- where the normalizing constant)
-
p(x k |y 1:k−1)=∫p(y k |x k)p(x k |y 1:k−1)dx k (12) - depends on the likelihood function p(yk|xk) defined by the measurement model (equation 2) and the known statistics of nk.
- It is not possible to sample directly from the posterior density function so importance sampling from a known proposal distribution π(x0:k|y1:k) is used. One may use sigma-point Kalman filters, for example, to generate the proposal.
- The known distribution is introduced into Equation 5 to yield:
-
- where the variables wk(x0:k) are unnormalized importance weights, which are written as wk(x0:k)=wk:
-
- resulting in a weighted expectation:
-
- where {tilde over (w)}k (i) are normalized importance weights:
-
- Importance sampling is made sequential by reiterating the
Markov 1st order assumption, resulting in the assumption that the current state is not dependent on future observations: -
π(x 0:k |y 1:k)=π(x 0:k−1 |y 1:k−1)π(x k |x 0:k−1 ,y 1:k) (17) - and that observations are conditionally independent given the states:
-
- A recursive estimate for the importance weights is:
-
- which is called Sequential Importance Sampling (SIS). SIS suffers from degeneracy so that, over a few iterations, all but one of the importance weights will be zero, effectively removing a large number of samples. To remedy this, samples with low importance weights may be eliminated while high importance samples may be multiplied. One way to accomplish this is Sampling-Importance Resampling (SIR), which involves mapping the Dirac random measure
-
{x k (i) ,{tilde over (w)} k (i) ;i=1, . . . ,N} (21) - into a measure with equal weights, 1/N:
-
- A pseudo-code for a generic SMC (also called bootstrap filter or condensation algorithm) can be written as:
- 1. Importance sampling step. For i=1, . . . , N, do:
-
- i) sample xk (i)˜p(xk|xk−1 (i))
- ii) evaluate wk (i)˜wk−1 (i)p(yk|xk (i))
- iii) normalize {tilde over (w)}k (i)=wk (i)/Σj=1 Nwk (j)
2. Importance resampling step - i) eliminate or multiply samples according to weights to obtain N random samples approximately distributed according to
- ii) For i=1, . . . , N, set wk (i)={tilde over (w)}k (i)=N−1
-
-
- i) any expectation, for instance:
-
- SPKF may be used to approximate probability distributions. Assuming that x has a mean
x covariance Px, and dimension L, a set of 2L+1 weighted sigma-points, Si={wi,Xi}, is chosen according to: -
- where h is a scaling parameter. Each sigma-point is propagated through the DSSM to yield the posterior sigma-point set, Yi:
-
Y i =h(f(X i)), i=0, . . . ,2L (24) - From this, the posterior statistics are calculated using a procedure resembling the linear Kalman filter. For instance, for the unscented Kalman filter case, a SPKF variant, the time-update equations are:
-
- and the measurement-update equations are:
-
- where x, v and n superscripts denote the state, process noise and measurement noise dimensions, respectively.
- The mathematical structure for sequential Monte Carlo and SPKF represent two examples of a family of probabilistic inference methods exploiting Monte Carlo simulation and the sigma point transform, respectively, in conjunction with a Bayesian statistical process.
- SPKF are generally inferior to SMC but are computationally cheaper. Like SMC, SPKF evolve the state using the full nonlinear DSSM, but represent probability distributions using a sigma-point set. This is a deterministic step that replaces the stochastic Monte Carlo step in the SMC. As a result, SPKF lose accuracy when posterior distributions depart heavily from the Gaussian form, such as with bimodal or heavily-tailed distributions, or with strong nonstationary distributions such as those caused by motion artifacts in pulse oximeters. For these cases SMC are more suitable.
- SPKF yields higher-order accuracy than the extended Kalman filter (EKF) and its related variants with equal algorithm complexity, O (L2). SPKF returns 2nd order accuracy for nonlinear and non-Gaussian problems, and 3rd order for Gaussian problems. EKF has only 1st order accuracy for nonlinear problems. Both EKF and SPKF approximate state distributions with Gaussian random variables (GRV). However, the EKF propagates the GRV using a single measure (usually the mean) and the 1st order Taylor expansion of the nonlinear system. The SPKF, on the other hand, decomposes the GRV into distribution moments (sigma points) and propagates those using the unmodified nonlinear system. SPKF implementation is simpler than EKF since it is derivativeless. That is, it uses the unmodified DSSM form, and therefore does not require lengthy Jacobian derivations.
- The data processing method is also capable of prediction because the method can operate faster than real time measurements. At any given time during data processing, the measurement PDF, obtained either from SPKF or SMC, embodies all available statistical information up to that point in time. It is therefore possible to march the system model forwards in time, for instance, using the same sequential Monte Carlo method, to obtain deterministic or stochastic simulations of future signal trajectories. In this way, the future health status (physiological state) of a patient can be predicted with attached probabilities indicating the confidence of each prediction.
- The data processing method may benefit from a noise adaptation method if timed sensor data contains noise and/or artifact that changes its spectral qualities over time. That is, if timed sensor data has a non-stationary probability distribution function. Here, a known algorithm such as the Robins-Monro or Annealing methods may be added to the data processing method in order to adapt the probability distribution functions of noise terms (stochastic terms) in the DSSM according to changing noise and artifact present in sensor data.
- In general, the output may include estimates of the true measured signals (i.e. processed data), and estimates of values for one or more physiological parameters measured by one or more sensors from which data was received, and estimates of values for one or more physiological parameters not measured by the sensors from which data was received (data extraction). A state parameter estimate is the processed data from the physiological sensor. Both noise and artifacts can be attenuated or rejected even though they may have very distinct probability distribution functions and may mimic the real signal. A model parameter estimate may be also used to produce a physiological parameter. For example, an estimate of total blood volume may be used to diagnose hemorrhage or hypovolemia. An estimate of tissue oxygen saturation may indicate poor tissue perfusion and/or hypoxia. Estimates of glucose uptake in several tissues may differentiate between diabetes mellitus types and severities; and estimates of carotid artery radius may be indicative of carotid artery stenosis.
- Pulse Oximeter with Probabilistic Data Processing
-
FIG. 8 shows the components of a DSSM suitable for processing data from a pulse oximeter, including components required to describe processes occurring in a subject.FIG. 9 illustrates the DSSM broken down into process and observation models, and including input and output variables. Heart rate (HR), stroke volume (SV) and whole-blood oxygen saturation (SpO2) are estimated from input noisy red and infrared intensity ratios (R). Radial (Pw) and aortic (Pao) pressures are also available as state estimates. - In this example, the DSSM comprises the following function to represent CO:
-
- wherein CO QCO(t), is expressed as a function of HR and SV and where QCO=(HR×SV)/60. The CO function pumps blood into a Windkessel 3-element model of the vascular system including two state parameters: aortic pressure, Pao, and radial (Windkessel) pressure, PW:
-
P ao,k+1 =P w,k+1 +Z 0 Q CO (33) - Rp and Zo are the peripheral resistance and characteristic aortic impedance, respectively. The sum of these two terms is the total peripheral resistance due to viscous (Poiseuille-like) dissipation:
-
Z 0=√{square root over (η/AC 1)} (34) - where ρ is blood density. The elastic component due to vessel compliance is a nonlinear function including thoracic aortic cross-sectional area, A:
-
- where Amax, P0 and P1 are fitting constants correlated with age and gender:
-
A max=(5.62−1.5(gender))·cm2 (36) -
P 0=(76−4(gender)−0.89(age))·mmHg (37) -
P 1(57−0.44(age))·mmHg (38) - The time-varying Windkessel compliance, Cw, and the aortic compliance per unit length, Cl, are:
-
- where l is the aortic effective length. The peripheral resistance is defined as the ratio of average pressure to average flow. A set-point pressure, Pset, and the instantaneous flow:
-
- are used to provide compensation autonomic nervous system responses. The value for Pset is adjusted manually to obtain 120 over 75 mmHg for a healthy individual at rest. The compliance of blood vessels changes the interactions between light and tissues with pulse. This is accounted for using a homogenous photon diffusion theory for a reflectance or transmittance pulse oximeter configuration. For the reflectance case:
-
- for each wavelength. In this example, the red and infrared bands are centered at ˜660 nm and ˜880 nm. I denotes the detected intensities: total reflected (no subscript), and the pulsating (ac) and background (dc) components. Va is the arterial blood volume, which changes as the cross-sectional area of illuminated blood vessels, ΔAw, changes as:
-
ΔV a ≈r·ΔA w (42) - where r is the source-detector distance. The tissue scattering coefficient, Σs′, is assumed constant but the arterial absorption coefficient, Σaart, depends on blood oxygen saturation, SpO2:
-
- which is the Beer-Lambert absorption coefficient, with hematocrit, H, and red blood cell volume, vi. The optical absorption cross sections for red blood cells containing totally oxygenated (HbO2) and totally deoxygenated (Hb) hemoglobin are σa 100% and σa 0%, respectively.
- The function K(α,d,r) contains, along with the scattering coefficient, the wavelength, sensor geometry and oxygen saturation dependencies that alter the effective optical pathlengths:
-
- The attenuation coefficient α is:
-
- where Σa and Σs are whole-tissue absorption and scattering coefficients, respectively, which are calculated from Mie Theory.
- Red and infrared K values as a function of SpO2 may be represented by two linear fits:
-
K r ≈−4.03·SpO 2−1.17 (46) -
K ir ≈0.102·SpO 2−0.753 (47) - in mm2. The overbar denotes the linear fit of the original function. The pulsatile behavior of ΔAw, which couples optical detection with the cardiovascular system model, is:
-
- with Pw,0=(1/3)P0 and Pw,1=(1/3)P1 to account for the poorer compliance of arterioles and capillaries relative to the thoracic aorta. Third and fourth state variables, the red and infrared reflected intensity ratios, R=lac/ldc, are:
-
R r,k+1 =cΣ′ z,rK r Σa,r art ΔA w +R r,k +v r (49) -
R ir,k+1 =cΣ′ z,irK ir Σa,ir art ΔA w +R ir,k +v ir (50) - Here, v are Gaussian-distributed process noises intended to capture the baseline wander of the two channels. The constant c subsumes all factors common to both wavelengths and is treated as a calibration constant. The observation model adds Gaussian-distributed noises, n, to Rr and Rir:
-
- A calibration constant c was used to match the variance of the real lac/ldc signal with the variance of the DSSM-generated signal for each wavelength. After calibration, the age and gender of the patient was entered. Estimates for the means and covariances of both state and parameter PDFs were entered.
FIG. 11 plots estimates for a 15 s stretch of data. Photoplethysmographic (PPG) waveforms (A) were used to extract heart rate (B), SV (C), CO (D), blood oxygen saturation (E), and aortic and systemic (radial) pressure waveforms (F). Results of processing pulse oximetry at low blood perfusion are shown inFIG. 12 . Low signal-to-noise PPG waveforms (A) were used to extract heart rate (B), SV (C), blood oxygen saturation (D), and aortic and systemic (radial) pressure waveforms (E). - Electrocardiograph with Probabilistic Data Processing
-
FIG. 10 is a schematic of a DSSM suitable for processing electrocardiograph data, including components required to describe the processes occurring in a subject. The combination of SPKF or SMC in state, joint or dual estimation modes can be used to filter electrocardiography (ECG) data. Any physiology model adequately describing the ECG signal can be used, as well as any model of noise and artifact sources interfering or contaminating the signal. One non-limiting example of such a model is the ECG signal generator proposed by McSharry (IEEE Transactions on Biomedical Engineering, 2003. 50(3):289-294). Briefly, this model uses a sum of Gaussians with amplitude, center and standard deviation, respectively, for each wave (P, Q, R, S, T−, T+15) in an ECG. The observation model comprises the state plus additive Gaussian noise, but more realistic pink noise or any other noise distributions can be used. -
FIG. 13 shows the results of processing a noisy non-stationary ECG signal. Heart rate oscillations representative of normal respiratory sinus arrhythmia are present in the ECG. The processor accomplishes accurate, simultaneous estimation of the true ECG signal and heart rate that follows closely the true values. The performance of the processor in a noise and artifact-corrupted signal is shown inFIG. 14 . A clean ECG signal representing one heart beat (truth) was contaminated with additive noise and an artifact in the form of a plateau at R and S peaks (beginning at time=10 s). Estimates by the processor remain close to the true signal despite the noise and artifact. - While a specific DSSMs and input and output parameters are provided for the purpose of describing the present method, the present invention is not limited to the DSSMs, sensors, biological monitoring devices, inputs, outputs, except as defined by the following claims.
Claims (15)
1. A computer implemented method for non-invasively measuring left ventricular stoke volume (SV) and/or cardiac output (CO) of a subject using a pulse oximeter, said method comprising:
receiving measured plethysmographic waveform data of a subject into a data processor comprising software, said software comprising a mathematical model of a cardiovascular system integrated in a dynamic state space model (DSSM);
calculating a measured value for CO and/or SV using the plethysmographic waveform data; and
reporting the measured SV and/or CO values.
2. The computer implemented method of claim 1 , further comprising collecting plethysmographic waveform data from the subject.
3. The computer implemented method of claim 2 , further comprising transmitting the plethysmographic waveform data of the subject to the data processor in real time.
4. The computer implemented method of claim 1 , wherein said calculating and reporting are performed in real time.
5. The computer implemented method of claim 1 , wherein the mathematical model of the cardiovascular system comprises aortic pressure, radial pressure, peripheral resistance, aortic impedance, heart rate, stroke volume, and blood density as state or model parameters.
6. The computer implemented method of claim 1 , wherein the DSSM is integrated with a dual estimation processor engine or a joint estimation processor engine.
7. The computer implemented method of claim 1 , wherein said reporting the measured SV and/or CO values to a user comprises producing an electronic visual display, a hard copy display, an audible sound, or a tactile output.
8. The computer implemented method of claim 1 , wherein the mathematical model of the cardiovascular system comprises model and/or state parameters that correspond directly to one or more of: peripheral blood oxygen saturation, heart rate, respiratory rate, and blood pressure.
9. The computer implemented method of claim 8 , further comprising receiving into the data processor, data of a direct measurement for one or more of SpO2, HR, RR, and BP.
10. The computer implemented method of claim 1 , wherein said software comprises code directing the data processor to:
a) receive system and model parameters for a time t into the DSSM to produce a first probability distribution function (PDF) vector comprising state and model parameters for time t+n;
b) use the first PDF vector and timed data obtained for time t+n from the plethysmographic waveform data in a Bayesian statistical process to produce a second PDF vector for state and model parameters for time t+n;
c) calculate probabilistic expectation values for the state and model parameters for time t+n from the second probability distribution function; and
d) determine a measured value for SV and/or CO for time t+n from probabilistic expectation values for the state and/or model parameters for time t+n
wherein:
the state and model parameters for a time t received into the DSSM in step a) are in the form of a probability distribution function produced from a sampling of expectation values calculated in step c) for an immediately preceding time t−n′; and n and n′ are time intervals that may be the same of different;
the DSSM mathematically represents physiological processes responsible for the measured plethysmographic waveform data and stroke volume to produce a time dependent state representing a time dependent physiological state of the subject;
the DSSM comprises at least one model parameter and/or state parameter representing at least one of total blood volume (TBV), stroke volume (SV), vasomotor tone (VT), and autonomous nervous system (ANS) tone; and
the software comprises code that determines a measured value for at least one of SV, CO, TBV, VT, and ANS tone.
11. The computer implemented method of claim 10 , wherein said measured value determined in step d) is equal to the value of a model parameter or a state parameter of said DSSM, or is calculated from the value of a model parameter and/or a state parameter of said DSSM.
12. The computer implemented method of claim 10 , wherein the DSSM is integrated in a joint estimation processing engine.
13. The computer implemented method of claim 10 , wherein the first PDF is produced using a Sequential Monte Carlo or Sigma Point Kalman Filter method.
14. The computer implemented method of claim 13 , wherein the Sigma Point Kalman Filter method is an unscented Kalman Filter, a central difference Kalman Filter, a squareroot unscented Kalman Filter, a square-root central difference Kalman Filter, or a combination thereof; and the Sequential Monte Carlo method is an unscented Monte Carlo method, a central difference Monte Carlo method, a square-root unscented Monte Carlo method, a square-root central difference Monte Carlo method, Gaussian Sum Monte Carlo method, Bayes Monte Carlo method, a Gaussian Mixture Sigma Point Monte Carlo method, or any combination thereof.
15. A non-transitory computer-readable storage medium storing a program that causes a computer to execute the method of claim 1 .
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US16/522,487 US20190343406A1 (en) | 2009-04-22 | 2019-07-25 | System and Method for Noninvasively Measuring Ventricular Stroke Volume and Cardiac Output |
Applications Claiming Priority (11)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US17180209P | 2009-04-22 | 2009-04-22 | |
US12/640,278 US20100274102A1 (en) | 2009-04-22 | 2009-12-17 | Processing Physiological Sensor Data Using a Physiological Model Combined with a Probabilistic Processor |
US12/796,512 US9060722B2 (en) | 2009-04-22 | 2010-06-08 | Apparatus for processing physiological sensor data using a physiological model and method of operation therefor |
US13/096,904 US9649036B2 (en) | 2009-04-22 | 2011-04-28 | Biomedical parameter probabilistic estimation method and apparatus |
US13/096,876 US9375171B2 (en) | 2009-04-22 | 2011-04-28 | Probabilistic biomedical parameter estimation apparatus and method of operation therefor |
US13/096,845 US9173574B2 (en) | 2009-04-22 | 2011-04-28 | Mechanical health monitor apparatus and method of operation therefor |
US13/181,027 US8494829B2 (en) | 2010-07-21 | 2011-07-12 | Sensor fusion and probabilistic parameter estimation method and apparatus |
US13/181,140 US9451886B2 (en) | 2009-04-22 | 2011-07-12 | Probabilistic parameter estimation using fused data apparatus and method of use thereof |
US13/181,247 US20120022336A1 (en) | 2010-07-21 | 2011-07-12 | Iterative probabilistic parameter estimation apparatus and method of use therefor |
US15/225,803 US10405762B2 (en) | 2009-04-22 | 2016-08-01 | System and method for noninvasively measuring ventricular stroke volume and cardiac output |
US16/522,487 US20190343406A1 (en) | 2009-04-22 | 2019-07-25 | System and Method for Noninvasively Measuring Ventricular Stroke Volume and Cardiac Output |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/225,803 Division US10405762B2 (en) | 2009-04-22 | 2016-08-01 | System and method for noninvasively measuring ventricular stroke volume and cardiac output |
Publications (1)
Publication Number | Publication Date |
---|---|
US20190343406A1 true US20190343406A1 (en) | 2019-11-14 |
Family
ID=58638049
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/225,803 Active 2030-12-24 US10405762B2 (en) | 2009-04-22 | 2016-08-01 | System and method for noninvasively measuring ventricular stroke volume and cardiac output |
US16/522,487 Pending US20190343406A1 (en) | 2009-04-22 | 2019-07-25 | System and Method for Noninvasively Measuring Ventricular Stroke Volume and Cardiac Output |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/225,803 Active 2030-12-24 US10405762B2 (en) | 2009-04-22 | 2016-08-01 | System and method for noninvasively measuring ventricular stroke volume and cardiac output |
Country Status (1)
Country | Link |
---|---|
US (2) | US10405762B2 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2875673A1 (en) * | 2020-05-07 | 2021-11-10 | Rgb Medical Devices Sa | Method, device and system to determine the trend of cardiac output (Machine-translation by Google Translate, not legally binding) |
EP3928689A1 (en) * | 2020-06-26 | 2021-12-29 | Ectosense NV | Apparatus and method for compensating assessment of peripheral arterial tone |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA3125985A1 (en) | 2019-01-16 | 2020-07-23 | Abiomed, Inc. | Left ventricular volume and cardiac output estimation using machine learning model |
CN109620198B (en) * | 2019-02-21 | 2022-04-12 | 天津惊帆科技有限公司 | Cardiovascular index detection and model training method and device |
JPWO2021066028A1 (en) * | 2019-09-30 | 2021-04-08 | ||
CN111493855B (en) * | 2020-04-21 | 2023-01-06 | 重庆理工大学 | System and method for non-invasive measurement of individualized cardiac output |
CN112580612B (en) * | 2021-02-22 | 2021-06-08 | 中国科学院自动化研究所 | Physiological signal prediction method |
CN113057617B (en) * | 2021-04-30 | 2022-08-26 | 重庆理工大学 | Non-invasive monitoring system for cardiac output |
CN114271805B (en) * | 2021-12-31 | 2023-07-25 | 四川大学 | Cardiac output measurement method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5626140A (en) * | 1995-11-01 | 1997-05-06 | Spacelabs Medical, Inc. | System and method of multi-sensor fusion of physiological measurements |
US20080183232A1 (en) * | 2007-01-30 | 2008-07-31 | Voss Gregory I | Method and system for determining cardiac function |
US7407486B2 (en) * | 2002-10-14 | 2008-08-05 | Ge Healthcare Finland Oy | Method and an apparatus for pulse plethysmograph based detection of nociception during anesthesia or sedation |
EP2047794A1 (en) * | 2007-10-11 | 2009-04-15 | Lidco Group Plc | Hemodynamic monitor |
-
2016
- 2016-08-01 US US15/225,803 patent/US10405762B2/en active Active
-
2019
- 2019-07-25 US US16/522,487 patent/US20190343406A1/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5626140A (en) * | 1995-11-01 | 1997-05-06 | Spacelabs Medical, Inc. | System and method of multi-sensor fusion of physiological measurements |
US7407486B2 (en) * | 2002-10-14 | 2008-08-05 | Ge Healthcare Finland Oy | Method and an apparatus for pulse plethysmograph based detection of nociception during anesthesia or sedation |
US20080183232A1 (en) * | 2007-01-30 | 2008-07-31 | Voss Gregory I | Method and system for determining cardiac function |
EP2047794A1 (en) * | 2007-10-11 | 2009-04-15 | Lidco Group Plc | Hemodynamic monitor |
Non-Patent Citations (2)
Title |
---|
Awad et al. (J. Clinical Monitoring and Computing, 2006, 20, 175-184 * |
Jason Thomas Ambadan and Youmin Tang. Sigma-Point Kalman Filter Data Assimilation Methods for Strongly Nonlinear Systems. American Meteorological Society 2009 (Year: 2009) * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2875673A1 (en) * | 2020-05-07 | 2021-11-10 | Rgb Medical Devices Sa | Method, device and system to determine the trend of cardiac output (Machine-translation by Google Translate, not legally binding) |
EP3928689A1 (en) * | 2020-06-26 | 2021-12-29 | Ectosense NV | Apparatus and method for compensating assessment of peripheral arterial tone |
WO2021260190A1 (en) * | 2020-06-26 | 2021-12-30 | Ectosense NV | Apparatus and method for compensating assessment of peripheral arterial tone |
Also Published As
Publication number | Publication date |
---|---|
US20170119261A1 (en) | 2017-05-04 |
US10405762B2 (en) | 2019-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10405762B2 (en) | System and method for noninvasively measuring ventricular stroke volume and cardiac output | |
US9060722B2 (en) | Apparatus for processing physiological sensor data using a physiological model and method of operation therefor | |
EP2512325B1 (en) | Processing physiological sensor data using a physiological model combined with a probabilistic processor | |
US9375171B2 (en) | Probabilistic biomedical parameter estimation apparatus and method of operation therefor | |
US9220440B2 (en) | Determining a characteristic respiration rate | |
US9649036B2 (en) | Biomedical parameter probabilistic estimation method and apparatus | |
US9451886B2 (en) | Probabilistic parameter estimation using fused data apparatus and method of use thereof | |
US10699206B2 (en) | Iterative probabilistic parameter estimation apparatus and method of use therefor | |
US9198582B2 (en) | Determining a characteristic physiological parameter | |
US20140275886A1 (en) | Sensor fusion and probabilistic parameter estimation method and apparatus | |
US20080183232A1 (en) | Method and system for determining cardiac function | |
US20130012823A1 (en) | Methods and Systems for Non-Invasive Measurement of Blood Pressure | |
US10460843B2 (en) | Probabilistic parameter estimation using fused data apparatus and method of use thereof | |
WO2008055173A2 (en) | Methods and systems for non-invasive, internal hemorrhage detection | |
JP2003511101A (en) | Apparatus and method for continuous non-invasive determination of physiological properties | |
US9949696B2 (en) | Apparatus and methods for computing cardiac output of a living subject via applanation tonometry | |
WO2010043728A1 (en) | System and apparatus for the non-invasive measurement of blood pressure | |
US11363994B2 (en) | Cardiovascular state determination apparatus and method of use thereof | |
US20210177275A1 (en) | Non-cerebral organ autoregulation status determination | |
Pielmus et al. | Surrogate based continuous noninvasive blood pressure measurement | |
WO2011110491A1 (en) | A non-invasive system and method for diagnosing and eliminating white coat hypertention and white coat effect in a patient | |
US20220400960A1 (en) | Autoregulation monitoring using deep learning | |
Alqudah et al. | Multiple time and spectral analysis techniques for comparing the PhotoPlethysmography to PiezoelectricPlethysmography with electrocardiography | |
US20230000445A1 (en) | Cardiovascular state determination apparatus and method of use thereof | |
Arnold et al. | Pulse oximeter plethysmograph waveform and automated oscillometric sphygmomanometer for ankle-brachial index measurement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STPP | Information on status: patent application and granting procedure in general |
Free format text: APPLICATION DISPATCHED FROM PREEXAM, NOT YET DOCKETED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: FINAL REJECTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |