US20060009710A1 - Apparatus and method for determination of stroke volume using the brachial artery - Google Patents

Apparatus and method for determination of stroke volume using the brachial artery Download PDF

Info

Publication number
US20060009710A1
US20060009710A1 US11/158,521 US15852105A US2006009710A1 US 20060009710 A1 US20060009710 A1 US 20060009710A1 US 15852105 A US15852105 A US 15852105A US 2006009710 A1 US2006009710 A1 US 2006009710A1
Authority
US
United States
Prior art keywords
time
patient
ttp
value
derivative
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.)
Granted
Application number
US11/158,521
Other versions
US7806830B2 (en
Inventor
Donald Bernstein
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Cordeus Inc
Original Assignee
Cordeus Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from US10/870,281 external-priority patent/US7261697B2/en
Application filed by Cordeus Inc filed Critical Cordeus Inc
Priority to US11/158,521 priority Critical patent/US7806830B2/en
Assigned to CORDEUS, INC. reassignment CORDEUS, INC. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: BERNSTEIN, DONALD P.
Priority to PCT/US2005/044671 priority patent/WO2006063255A2/en
Publication of US20060009710A1 publication Critical patent/US20060009710A1/en
Priority to US12/489,335 priority patent/US7740590B2/en
Application granted granted Critical
Publication of US7806830B2 publication Critical patent/US7806830B2/en
Active - Reinstated legal-status Critical Current
Adjusted expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/026Measuring blood flow
    • A61B5/029Measuring or recording blood output from the heart, e.g. minute volume
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/02028Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/026Measuring blood flow
    • A61B5/0295Measuring 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/053Measuring electrical impedance or conductance of a portion of the body
    • A61B5/0535Impedance plethysmography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7239Details of waveform analysis using differentiation including higher order derivatives
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/352Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval

Definitions

  • This present invention relates to the determination of the volumetric output of the left ventricle of a person's heart per beat, known as stroke volume (SV) (mL), and, the volumetric output of a person's heart per minute, otherwise known as the cardiac output (CO) (L/Min). More particularly, this invention relates to the determination of SV and CO by transbrachial electrical bioimpedance.
  • SV stroke volume
  • CO cardiac output
  • transthoracic method also known as transthoracic or thoracic electrical bioimpedance plethysmography (or cardiography), or by total body (whole body) electrical bioimpedance plethysmography, also known as whole body electrical bioimpedance cardiography (Moshkovitz Y, et al. Curr Opin Cardiol 2004; 19:229-237).
  • transthoracic method also known as transthoracic or thoracic electrical bioimpedance plethysmography (or cardiography)
  • total body electrical bioimpedance plethysmography also known as whole body electrical bioimpedance cardiography
  • the present invention is an apparatus for determining stroke volume by bioimpedance from a patient, including two or more spaced apart alternating current flow electrodes positionable on a patient, two or more spaced apart voltage sensing electrodes positionable on the patient and between the alternating current flow electrodes, an alternating current source electrically connected to the alternating current flow electrodes, a voltmeter electrically connected to the voltage sensing electrodes, and a processing unit in communication with the voltage sensing electrodes.
  • the processing unit is capable of using a voltage sensed by the voltage sensing electrodes to calculate a mean value of a second time-derivative of a cardiogenically induced impedance variation of the patient and to determine therefrom a stroke volume of the patient.
  • a method of determining stroke volume by bioimpedance from a patient includes positioning two or more spaced apart alternating current flow electrodes on a patient, positioning two or more spaced apart voltage sensing electrodes on the patient and between the alternating current flow electrodes, providing an alternating current flow (I(t)) through the electrically conductive electrodes creating a current field, measuring a voltage (U(t)) between the voltage sensing electrodes within the current field, calculating a mean value of a second time-derivative of a cardiogenically induced impedance variation of the patient using the measured voltage (U(t)), and calculating the stroke volume (SV) of the patient using the calculated second time-derivative mean value.
  • FIG. 1 shows placement of electrodes on a patient.
  • A.C.(I) is injected through a segment of the upper arm, otherwise known as the brachium, the boundaries of which are the deltoid muscles of the shoulder and axilla proximally, and the elbow and antecubital fossa distally.
  • the upper arm, including the connective tissue, bone, nervous tissue, veins, and the brachial artery comprise an aggregate impedance (Z) to current flow.
  • A.C across the brachium generates a quasi-static voltage, (U 0 ), and, concordant with every pressure pulse of the brachial artery, a time-dependent drop in voltage, ( ⁇ U(t)), this pressure pulse following every onset of left ventricular ejection with a short time delay.
  • FIGS. 3 a and 3 b show the relationship between the dZ/dt curve and the dP/dt or d(SpO 2 )/dt curve.
  • FIG. 3 a further shows an example where points B and X are apparent on the dZ/dt curve and
  • FIG. 3 b shows an example where point B is not detectable, but point X is detectable on the dZ/dt curve.
  • FIG. 4 shows the primary waveforms of ⁇ SpO 2 (t) and/or ⁇ P(t), aligned in time with the dZ/dt waveform.
  • FIG. 5 shows an example where points B and X are distinguishable and that point B corresponds with aortic valve opening (AVO) on the first time-derivatives of either the ⁇ SpO 2 (t)or ⁇ P(t) waveforms and point X corresponds with aortic valve closing (AVC) of either derivative.
  • AVO aortic valve opening
  • AVC aortic valve closing
  • FIG. 6 shows a dZ/dt waveform where points B and X are not distinguishable, and where point C (dZ/dt max ) is aligned in time with dP/dt max(radial) .
  • FIG. 7 shows placement of the electrodes on a patient implementing the transthoracic approach, which, as described herein, is required for calibration of the transbrachial approach if auto-calibration is not employed.
  • A.C. (I) is injected through a segment of the thorax (chest) between the base of the neck (laterally) 210 and lower thorax (laterally) 212 at the level of the xiphoid process (inferior portion of the sternum, or breast bone) in the mid-axillary line.
  • an A.C. field is applied to the thoracic volume between points 214 and 216 , forcing an A.C.
  • the A.C. causes, in the direction of the electrical field, and between the current injecting electrodes, a measured voltage, U(t).
  • U(t) is further comprised of a static D.C. component, U 0 , and a dynamic A.C. component, ⁇ U(t).
  • the voltage, U 0 , and voltage drop, ⁇ U(t) are sensed by electrodes proximate the current injecting electrodes, and within the current field.
  • An A.C. generator 218 and voltmeter 220 are shown.
  • FIG. 8 shows waveform examples of the following: ECG, ⁇ Z(t), dZ/dt, and d 2 Z/dt 2 ( ⁇ /s 3 ) from the transthoracic approach.
  • Fiducial landmarks noted on the dZ/dt waveform are point B, denoting aortic valve opening, point C, denoting dZ/dt max , and point X, denoting aortic valve closure.
  • Fiducial landmarks noted on the d 2 Z/dt 2 waveform are point B-1 indicating d 2 Z/dt 2 max and corresponding proximate in time to aortic valve opening; point C-1 corresponding to the first zero crossing and thus dZ/dt max ; and point X-1, corresponding to the second zero crossing and aortic valve closure (point X on the dZ/dt waveform).
  • the magnitude, d 2 Z/dt 2 max is noted.
  • the terms “comprises,” “comprising,” “includes,” “including,” “has,” “having” or any other variation thereof, are intended to cover a non-exclusive inclusion.
  • a process, method, article, or apparatus that comprises a list of elements is not necessarily limited to only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus.
  • “or” refers to an inclusive or and not to an exclusive or. For example, a condition A or B is satisfied by any one of the following: A is true (or present) and B is false (or not present), A is false (or not present) and B is true (or present), and both A and B are true (or present).
  • the present invention discloses a method and apparatus for the determination of stroke volume (SV) and cardiac output (CO) by transbrachial electrical bioimpedance, wherein the signal source is the brachial artery.
  • SV and CO while not sensitive indices of the overall intrinsic force generation capacity, or contractility of the heart muscle, are the best indicators of the overall performance of the heart considered as a muscular pump.
  • the apparatus and method disclosed involve the application of a constant magnitude alternating current of high frequency and small amplitude across a segment of a person's upper extremity, and more specifically, the upper arm, otherwise known as the brachium.
  • the present invention may also provide for calibrating the transbrachial method and apparatus by determining SV/CO from the transthoracic approach.
  • the present invention relates to the acquisition and signal processing of the cardiogenically-induced, pulsatile transbrachial bioimpedance signal for the purpose of SV/CO determination.
  • Baseline transbrachial quasi-static impedance, Z 0 is not affected by pulmonary (lung) ventilation, thereby obviating the necessity for sophisticated stabilizing adaptive filtering techniques to obtain a steady baseline for measurement of the cardiac-induced transbrachial impedance change, ⁇ Z(t), and the magnitudes and fiducial landmarks on its first time-derivative, transbrachial dZ/dt and on its second time-derivative, d 2 Z/dt 2 .
  • the present invention relates to the measurement of stroke volume (SV) and cardiac output (CO) from the transbrachial method, using the brachial artery as the cardiogenically-induced signal source.
  • the transbrachial method is similar to the transthoracic technique for determining SV.
  • signal acquisition is effected over a segment of thorax by placement of voltage sensing electrodes 214 at the base of the neck, bilaterally, and voltage sensing electrodes 216 at the lower thorax at the xiphoid level, bilaterally (see FIG. 7 ).
  • the transbrachial technique uses a segment of brachium between voltage sensing electrodes 112 / 114 positioned proximate the axilla (arm pit) and junction of the upper and lower arm at the level of the olecranon process of the elbow. (see FIG. 1 ).
  • FIG. 1 schematically shows one apparatus embodiment according to the present invention, and its electrical interface with a subject 100 .
  • Signal acquisition from the upper arm 102 (brachium) requires application of a constant magnitude alternating current (A.C.) 104 of high frequency and small amplitude to electrodes 106 , 108 that are spaced apart, with one or more electrodes affixed to the skin of the axilla, as well as one or more electrodes placed medially at the level of the antecubital fossa creating a current field.
  • the electrodes are applied to the subject's left arm. In other embodiments, the electrodes may be positioned on the right arm.
  • the potential difference between the current injecting electrodes or alternating current flow electrodes 106 , 108 is measured by a voltmeter 110 connected to the voltage sensing electrodes 112 , 114 placed within the current field (see FIG. 1 ).
  • a baseline impedance between the voltage sensing electrodes 112 , 114 , as well as a change in impedance, ⁇ Z(t) can be measured transbrachially.
  • dZ/dt (brachium) results, its peak systolic magnitude being dZ/dt max(brachium) .
  • d 2 Z/dt 2 and d 2 Z/dt 2 max are equivalent to the rate of change and peak negative rate of change of dZ/dt and dZ/dt min , respectively.
  • dZ/dt max being equivalent to dZ/dt min
  • Many different methods of applying the electrodes or electrode arrays to the arm are envisioned, such as spot electrodes, arm band(s) both circumferential and non-circumferential, adhesive strips or other attachment means known in the art.
  • an 8 spot-electrode array can be implemented.
  • a 4 spot-electrode array placed on the inner, or medial aspect of the upper arm, proximate the brachial artery, can be implemented.
  • the voltages measured by the Voltmeter 110 not only contains a signal caused by the AC applied, but may also contain a signal component from which an electrocardiogram (ECG) can be derived.
  • ECG electrocardiogram
  • the application of filters separates the AC related and ECG related signal components.
  • EKG 116 may also be measured by placing EKG electrodes 118 on the patient 100 .
  • a 3-lead EKG is shown and EKG is measured by known means.
  • the magnitude of the alternating current (A.C.) 104 and voltmeter 110 may be components of an apparatus 120 .
  • the apparatus 120 may also include an input device and a processor.
  • the input device may be any suitable device that provides information to the apparatus, such as a keyboard.
  • the input device may also receive information from external sources, such as the EKG 116 .
  • the processor is in communication with the data input device, the alternating current source 104 and electrodes 106 , 108 , and the voltmeter 110 and electrodes 112 , 114 .
  • the processor is capable of receiving the information and calculating the stroke volume and cardiac output of the patient 100 .
  • the stroke volume and cardiac output of the patient may be displayed on a screen or be sent to other devices via a data output device of the apparatus.
  • dZ/dt waveform to the right shows dZ/dt max remaining constant throughout the ejection interval, T LVE .
  • Rationale for use of the brachium as an appropriate anatomic site for SV measurement by the bioimpedance technique is as follows.
  • A.C. (I) is injected through a segment of upper arm, otherwise known as the brachium, the boundaries of which are the deltoid muscles of the shoulder and axilla, proximally, and the elbow and antecubital fossa, distally, a quasi-static voltage, U 0 , and voltage change, ⁇ U(t), can be measured between the current injecting electrodes.
  • the brachial artery is a large artery, continuous with both the subclavian and axillary arteries, and, whereas the left subclavian artery is a major branch of the arch of the thoracic aorta, the right subclavian artery is a branch of the brachiocephalic artery.
  • the contents of the upper arm, including connective tissue, bone, nervous tissue, veins, and the brachial artery, comprise an impedance (Z) to current flow.
  • Z impedance
  • dZ/dt max(brachial) represents the ohmic analog of peak blood acceleration in the brachial artery.
  • brachial artery blood velocity is affected by downstream peripheral vasoactivity (vasodilation, vasoconstriction)
  • the magnitude of brachial artery blood acceleration is modulated only by beta ( ⁇ ) adrenergic stimulation or depression of the cardiac adrenoceptors (Chemla D, et al. Am J Cardiol 1990; 65:494-500).
  • dZ/dt max /Z 0 (1/s 2 ) undergoes square root transformation
  • ohmic mean velocity results, ⁇ Z v (t) max /Z 0 (1/s).
  • dZ/dt max represents ohmic mean acceleration
  • dZ/dt max /Z 0 is herein referred to as the acceleration index (ACI).
  • ACI acceleration index
  • SV TB V C ⁇ [ ( d 2 ⁇ Z / d t mean 2 Z 0 ) ⁇ 10 - 2 ] 3 ⁇ T LVE equation ⁇ ⁇ 1
  • SV TB transbrachial SV (mL); V c(brachium) equals the volume conductor (mL); Z 0 equals the quasi-static transbrachial base impedance (Ohm, ⁇ ); T lve equals left ventricular ejection time (s), and d 2 Z/dt 2 mean equals the mean value of the second time-derivative of the cardiogenically induced transbrachial impedance variation ( ⁇ /s 3 ).
  • dZ(t)/dt max the present technique assumes dZ(t)/dt max to represent the ohmic equivalent of the peak acceleration of red blood cells.
  • dv/dt max the peak red blood cell acceleration
  • dV/dt max the peak rate of change of volume
  • the mean value of the second time-derivative of ⁇ Z(t) has certain advantages over processing the first time-derivative (e.g. time to peak values calculated more accurately because of using cube root calculations over square root calculations), which are provided in another preferred embodiment using the transthoracic approach.
  • the cube root transformation is thus implemented in the preferred embodiment.
  • the signal processing technique comprising part of the invention, implies that the proper designation for the transbrachial approach is correctly stated as Transbrachial Electrical Bioimpedance Cardiovelocimetry or, simply, Transbrachial bioimpedance velocimetry.
  • FIGS. 3 a and 3 b show the relationship between the dZ/dt curve and either the dP/dt or d(SpO 2 )/dt curves.
  • FIG. 3 a further shows an example where points B and X are apparent on the dZ/dt curve and
  • FIG. 3 b shows an example where point B is not detectable, but point X is detectable on the dZ/dt curve.
  • the said alternative/obligatory means for determining T lve are those obtained from means such as from the waveform corresponding to the photoplethysmographic pulse oximetry waveform, ⁇ SpO 2 (t), or its first time-derivative, d(SpO 2 )/dt, and/or by the waveform obtained from a non-invasive applanated radial arterial pressure pulse waveform, ⁇ P(t)(radial), or its first time-derivative, dP/dt (radial) . (see FIG. 4 )
  • the said means for determining point B on the transbrachial dZ/dt curve are those methods used for determining T lve when point X on the transbrachial dZ/dt curve, or its first time derivative (d 2 Z/dt 2 ), are identifiable by those skilled in the art of bioimpedance curve analysis.
  • point X is not identifiable on the transbrachial dZ/dt curve, or its first time-derivative, d 2 Z/dt 2 , then alternative means for point B detection are necessary.
  • said means for point B detection include use of the first time-derivative of the applanated radial pressure waveform tracing, dP/dt (radial) . It will be clear to those skilled in the art of curve analysis, why the aforementioned said means are superior to those disclosed by others, and most recently by Baura et al. (U.S. Pat. No. 6,561,986 B2).
  • dZ/dt waveforms such as those shown in FIG. 3 a
  • point B and point X are readily distinguishable by one skilled in the art of curve analysis.
  • these fiducial landmarks are frequently distorted by motion and ventilation artifacts (especially using the transthoracic approach), as well as by certain disease processes.
  • LVET may be more accurately measured by curve analysis of the pulse oximetry and applanation tonometry waveforms ( FIG. 3 b ), or their first time-derivatives. In one embodiment of the invention, either or both methods may be implemented. Of these techniques, applanation tonometry is most likely to demonstrate a dicrotic notch, and, therefore, is considered the preferred technique.
  • the points coinciding with the beginning and end of ejection can be readily identified from the first time-derivative curves of both the oximetry and applanation tonometry waveforms; namely, d(SpO 2 )/dt and dP/dt (radial) .
  • the best method constitutes computer analysis of the first time-derivatives.
  • regression equations for heart rate versus LVET may be implemented.
  • Point B on the transbrachial dZ/dt waveform is known to coincide with aortic valve opening, albeit with a time delay.
  • Exemplary dZ/dt waveforms demonstrate a distinct change in slope at, or not uncommonly above the zero baseline impedance, followed by a steep, positive linear segment ending at point C, or dZ/dt max .
  • dZ/dt max When a distinct change in slope leading to point C is detected at or above the baseline, one skilled in the art of curve analysis can readily identify point B.
  • Debski T T et al. Biol Psychol 1993; 36:63-74
  • using the transthoracic method despite using fiducial landmarks on the time-derivatives of dZ/dt (i.e.
  • the method disclosed herein provides a new and innovative solution for point B detection.
  • the new method employs one, or a combination of methods disclosed under determination of LVET; namely, ⁇ SpO 2 (t) and/or ⁇ P(t) (radial) (as shown in FIG. 4 ), or, respectively, their time derivatives, d(SpO 2 )/dt and/or dP/dt (radial) (as shown in FIG.
  • the technique of point B detection involves computerized curve fitting and alignment in time of temporal landmark X on the transbrachial dZ/dt curve with the dicrotic notch equivalent of one or both of the measured aforementioned oximetry and pressure curves, and/or preferably with one or both of their first time-derivatives.
  • one or both of the first time-derivative curves can be aligned in time with the transbrachial dZ/dt curve, such that the temporal point of the termination of flow, or aortic valve closure (AVC) equivalent on the first derivative oximetry or pressure curves, can be aligned in time with point X of the transbrachial dZ/dt curve.
  • AVC aortic valve closure
  • Point B coinciding with aortic valve opening (AVO), and the beginning of flow, albeit with a time delay, is identified by determining the temporal point on the transbrachial dZ/dt curve, intersecting, and coinciding in time with the point of onset of flow/pressure on the ⁇ SpO 2 (t)/AP(t) curves, and/or on their first time-derivatives.
  • This temporal point is identified as a discreet point at the baseline occurring before the first positive maximum upslope measured from foot of the respective baselines of the ⁇ SpO 2 (t) and/or ⁇ P(t) (radial) curves, and/or from their first time-derivatives ( FIGS. 3 a , 3 b , 4 , 5 , 6 ).
  • said means requires alignment in time of the earliest maximum positive peak of dP/dt (dP/dt max ) with point C of the transbrachial dZ/dt curve.
  • point B can be identified by applying a perpendicular through, and coinciding in time with the onset at baseline of the first positive deflection of dP/dt, where said perpendicular line must intersect the dZ/dt curve at or above baseline impedance.
  • the point of intersection of the perpendicular with the transbrachial dZ/dt curve is designated point B.
  • transbrachial d 2 Z/dt 2 max is the measured peak positive deflection of the d 2 Z/dt 2 (d 2 Z/dt 2 max ) curve usually occurring temporally proximate point B on the transbrachial dZ/dt curve.
  • ⁇ ⁇ V c ⁇ ( thorax ) C 1 ⁇ ( thorax ) ⁇ [ W ⁇ ( kg ) ⁇ C 2 ] equation ⁇ ⁇ 4
  • 0.10 ⁇ C 1(thorax) ⁇ 0.75, wherein C 1 0.25 in the preferred embodiment,
  • SV transbrachial V C ⁇ ( cal ) ⁇ [ ( d 2 ⁇ Z / d t mean 2 Z 0 ) transbrachial ⁇ 10 - 2 ] 3 ⁇ T LVE equation ⁇ ⁇ 5
  • Equation 7 determining V c(cal) from equation 3, the mean value of the constant, C 1(brachium) , can be found for the general population as follows;
  • C 1(brachium) (mean) [( C 1-1 +C 1-2 +C 1-4 + . . . C 1-n )/ n] equation 8
  • SV transbrachial [ C 1 ⁇ ( mean , brachial ) ⁇ W ⁇ C 2 ] ⁇ [ ( d 2 ⁇ Z / d t mean 2 Z 0 ) transbrachial ⁇ 10 - 2 ] 3 ⁇ T LVE equation ⁇ ⁇ 9
  • the SV of the left ventricle can be determined using the transthoracic approach and the mean value of the second time-derivative of the cardiogenically induced transthoracic impedance variation ( ⁇ /s 3 ).
  • ⁇ /s 3 the mean value of the second time-derivative of the cardiogenically induced transthoracic impedance variation
  • a tetrapolar spot electrode array can be applied to a person's body. The description of signal acquisition and processing is precisely that described above with respect to FIG. 2 .
  • d 2 Z/dt 2 mean is determined identically as described above with regard to equations 1 and 1a-1g.
  • T LVE left ventricular ejection time (seconds, sec, s).

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physics & Mathematics (AREA)
  • Physiology (AREA)
  • Hematology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

Provided herein are methods and apparatus for stroke volume determination by bioimpedance from a patient's upper arm, or brachium, or a patient's thorax, utilizing pulsations of the arteries contained therein. The apparatus includes two or more spaced apart alternating current flow electrodes positioned on the patient's arm or thorax and two or more spaced apart voltage sensing electrodes positioned on the patient's arm or thorax and in-between alternating current flow electrodes. The system and method utilizes the mean value of the second time-derivative of the cardiogenically induced impedance variation of the patient using the measured voltage from the voltage sensors in calculating the stroke volume of the patient.

Description

  • This application is a continuation-in-part of U.S. application Ser. No. 10/870,281, filed Jun. 16, 2004, and claims priority from U.S. Provisional Patent Application No. 60/634,616, filed Dec. 8, 2004.
  • BACKGROUND OF THE INVENTION
  • 1. Field of the Invention
  • This present invention relates to the determination of the volumetric output of the left ventricle of a person's heart per beat, known as stroke volume (SV) (mL), and, the volumetric output of a person's heart per minute, otherwise known as the cardiac output (CO) (L/Min). More particularly, this invention relates to the determination of SV and CO by transbrachial electrical bioimpedance.
  • 2. Background Information
  • All methods, apparatus and inventions related to the measurement of SV/CO by the electrical bioimpedance method have heretofore been implemented either by the transthoracic method, also known as transthoracic or thoracic electrical bioimpedance plethysmography (or cardiography), or by total body (whole body) electrical bioimpedance plethysmography, also known as whole body electrical bioimpedance cardiography (Moshkovitz Y, et al. Curr Opin Cardiol 2004; 19:229-237). Apart from a velocimetric method and apparatus described by Bernstein et al. (U.S. Pat. No. 6,511,438 B2), all prior art assumes a plethysmographic origin for the measured impedance change with respect to time (ΔZ(t)), and its peak rate of change (dZ/dtmax), coinciding with each beat of the heart (Moshkovitz Y, et al. Curr Opin Cardiol 2004; 19:229-237). The plethysmograghic-based transthoracic SV equations used clinically basically comprise two methods; they are described in U.S. Pat. No. 6,511,438 B2, and are known as the Nyboer-Kubicek equation (Kubicek equation) and the Sramek-Bernstein equation. The deficiencies of the method and apparatus invented by Bernstein et al., disclosed in U.S. Pat. No. 6,511,438 B2, include the following:
      • 1. A volume conductor, Vc, which underestimates the intrathoracic blood volume (ITBV) by approximately 15-20%
      • 2. The implementation of a square root function for heart rate (H.R.) frequency (i.e. √{square root over (f)}0=1/(TRR)0.5=(H.R./60)0.5) which is superfluous and unnecessary.
      • 3. A best method in the preferred embodiment for determining left ventricular ejection time, Tlve, is not disclosed.
      • 4. A best method in the preferred embodiment for determining point B is not disclosed
      • 5. A best method in the preferred embodiment for determining dZ/dtmax, based on the accurate determination of point B, is not disclosed
  • There are numerous drawbacks to the current methods and apparatus used for measurement of the transthoracic electrical bioimpedance stroke volume parameters. What is needed is an alternative approach to the transthoracic electrical bioimpedance determination of stroke volume; specifically, an alternative site for signal acquisition, and better methods to measure the independent variables comprising the stroke volume equation.
  • SUMMARY OF THE INVENTION
  • The present invention is an apparatus for determining stroke volume by bioimpedance from a patient, including two or more spaced apart alternating current flow electrodes positionable on a patient, two or more spaced apart voltage sensing electrodes positionable on the patient and between the alternating current flow electrodes, an alternating current source electrically connected to the alternating current flow electrodes, a voltmeter electrically connected to the voltage sensing electrodes, and a processing unit in communication with the voltage sensing electrodes. The processing unit is capable of using a voltage sensed by the voltage sensing electrodes to calculate a mean value of a second time-derivative of a cardiogenically induced impedance variation of the patient and to determine therefrom a stroke volume of the patient.
  • In another aspect of the present invention, a method of determining stroke volume by bioimpedance from a patient includes positioning two or more spaced apart alternating current flow electrodes on a patient, positioning two or more spaced apart voltage sensing electrodes on the patient and between the alternating current flow electrodes, providing an alternating current flow (I(t)) through the electrically conductive electrodes creating a current field, measuring a voltage (U(t)) between the voltage sensing electrodes within the current field, calculating a mean value of a second time-derivative of a cardiogenically induced impedance variation of the patient using the measured voltage (U(t)), and calculating the stroke volume (SV) of the patient using the calculated second time-derivative mean value.
  • Other objects and features of the present invention will become apparent by a review of the specification, claims and appended figures.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • The drawings accompanying and forming part of this specification are included to depict certain aspects of the invention. A clearer conception of the invention, and of the components and operation of systems provided with the invention, will become more readily apparent by referring to the exemplary, and therefore nonlimiting, embodiments illustrated in the drawings, wherein identical reference numerals designate the same components. The invention may be better understood by reference to one or more of these drawings in combination with the description presented herein. It should be noted that the features illustrated in the drawings are not necessarily drawn to scale.
  • FIG. 1 shows placement of electrodes on a patient. A.C.(I) is injected through a segment of the upper arm, otherwise known as the brachium, the boundaries of which are the deltoid muscles of the shoulder and axilla proximally, and the elbow and antecubital fossa distally. Embedded within the brachial musculature and connective tissue, and anatomically situated medial to the brachial bone, otherwise known as the humerus, is the brachial artery. The upper arm, including the connective tissue, bone, nervous tissue, veins, and the brachial artery comprise an aggregate impedance (Z) to current flow. The passage of A.C across the brachium generates a quasi-static voltage, (U0), and, concordant with every pressure pulse of the brachial artery, a time-dependent drop in voltage, (ΔU(t)), this pressure pulse following every onset of left ventricular ejection with a short time delay.
  • FIG. 2 shows an example of ECG, ΔZ(t) (Ω) and dZ/dt (Ω/s2) waveforms obtained transthoracically from a human subject, where TRR=the R-R interval, or the time for one cardiac cycle (seconds, s); Q=onset of ventricular depolarization; ---------=maximum systolic upslope extrapolation of ΔZ(t); B=aortic valve opening; C=peak rate of change of the thoracic cardiogenic impedance variation, dZ/dtmax (Ω/s2); X=aortic valve closure; Y=pulmonic valve closure; O=rapid ventricular filling wave; Q-B interval=pre-ejection period, TPE (seconds, s); B-C interval=time-to-peak dZ/dt, TTP (seconds, s); B-X interval=left ventricular ejection period, TLVE (seconds, s). dZ/dt waveform to the right shows dZ/dtmax remaining constant throughout the ejection interval, TLVE, which represents outflow compensation.
  • FIGS. 3 a and 3 b show the relationship between the dZ/dt curve and the dP/dt or d(SpO2)/dt curve. FIG. 3 a further shows an example where points B and X are apparent on the dZ/dt curve and FIG. 3 b shows an example where point B is not detectable, but point X is detectable on the dZ/dt curve.
  • FIG. 4 shows the primary waveforms of ΔSpO2(t) and/or ΔP(t), aligned in time with the dZ/dt waveform.
  • FIG. 5 shows an example where points B and X are distinguishable and that point B corresponds with aortic valve opening (AVO) on the first time-derivatives of either the ΔSpO2(t)or ΔP(t) waveforms and point X corresponds with aortic valve closing (AVC) of either derivative.
  • FIG. 6 shows a dZ/dt waveform where points B and X are not distinguishable, and where point C (dZ/dtmax) is aligned in time with dP/dtmax(radial).
  • FIG. 7 shows placement of the electrodes on a patient implementing the transthoracic approach, which, as described herein, is required for calibration of the transbrachial approach if auto-calibration is not employed. As shown, A.C. (I) is injected through a segment of the thorax (chest) between the base of the neck (laterally) 210 and lower thorax (laterally) 212 at the level of the xiphoid process (inferior portion of the sternum, or breast bone) in the mid-axillary line. As operationally implemented, an A.C. field is applied to the thoracic volume between points 214 and 216, forcing an A.C. of high frequency (50-100 kH) and low magnitude (1.0-4.0 mA (rms)) to flow longitudinally between the neck and lower thorax. The A.C. causes, in the direction of the electrical field, and between the current injecting electrodes, a measured voltage, U(t). U(t) is further comprised of a static D.C. component, U0, and a dynamic A.C. component, ΔU(t). The voltage, U0, and voltage drop, ΔU(t), are sensed by electrodes proximate the current injecting electrodes, and within the current field. An A.C. generator 218 and voltmeter 220 are shown.
  • FIG. 8 shows waveform examples of the following: ECG, ΔZ(t), dZ/dt, and d2Z/dt2 (Ω/s3) from the transthoracic approach. Fiducial landmarks noted on the dZ/dt waveform are point B, denoting aortic valve opening, point C, denoting dZ/dtmax, and point X, denoting aortic valve closure. Fiducial landmarks noted on the d2Z/dt2 waveform are point B-1 indicating d2Z/dt2 max and corresponding proximate in time to aortic valve opening; point C-1 corresponding to the first zero crossing and thus dZ/dtmax; and point X-1, corresponding to the second zero crossing and aortic valve closure (point X on the dZ/dt waveform). The magnitude, d2Z/dt2 max, is noted.
  • DESCRIPTION OF PREFERRED EMBODIMENTS
  • The invention and the various features and advantageous details thereof are explained more fully with reference to the nonlimiting embodiments that are illustrated in the accompanying drawings and detailed in the following description. Descriptions of well known starting materials, processing techniques, components, and equipment are omitted so as not to unnecessarily obscure the invention in detail. It should be understood, however, that the detailed description and the specific examples, while indicating preferred embodiments of the invention, are given by way of illustration only and not by way of limitation. Various substitutions, modifications, additions and/or rearrangements within the spirit and/or scope of the underlying inventive concept will become apparent to those skilled in the art from this disclosure.
  • As used herein, the terms “comprises,” “comprising,” “includes,” “including,” “has,” “having” or any other variation thereof, are intended to cover a non-exclusive inclusion. For example, a process, method, article, or apparatus that comprises a list of elements is not necessarily limited to only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Further, unless expressly stated to the contrary, “or” refers to an inclusive or and not to an exclusive or. For example, a condition A or B is satisfied by any one of the following: A is true (or present) and B is false (or not present), A is false (or not present) and B is true (or present), and both A and B are true (or present).
  • Also, use of the “a” or “an” are employed to describe elements and components of the invention. This is done merely for convenience and to give a general sense of the invention. This description should be read to include one or at least one and the singular also includes the plural unless it is obvious that it is meant otherwise.
  • Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Although methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present invention, suitable methods and materials are described below. All publications, patent applications, patents, and other references mentioned herein are incorporated by reference in their entirety. In case of conflict, the present specification, including definitions, will control. In addition, the materials, methods, and examples are illustrative only and not intended to be limiting.
  • The present invention discloses a method and apparatus for the determination of stroke volume (SV) and cardiac output (CO) by transbrachial electrical bioimpedance, wherein the signal source is the brachial artery. SV and CO, while not sensitive indices of the overall intrinsic force generation capacity, or contractility of the heart muscle, are the best indicators of the overall performance of the heart considered as a muscular pump. The apparatus and method disclosed involve the application of a constant magnitude alternating current of high frequency and small amplitude across a segment of a person's upper extremity, and more specifically, the upper arm, otherwise known as the brachium. The present invention may also provide for calibrating the transbrachial method and apparatus by determining SV/CO from the transthoracic approach. Thus, in contradistinction to the generally accepted transthoracic bioimpedance method for SV/CO determination, the present invention relates to the acquisition and signal processing of the cardiogenically-induced, pulsatile transbrachial bioimpedance signal for the purpose of SV/CO determination.
  • Advantages of the transbrachial method include:
      • 1. Stroke volume (SV) and cardiac output (CO) values are not affected by excess, extra-vascular, intrathoracic liquids; namely, pulmonary edema fluid.
  • 2. Baseline transbrachial quasi-static impedance, Z0, is not affected by pulmonary (lung) ventilation, thereby obviating the necessity for sophisticated stabilizing adaptive filtering techniques to obtain a steady baseline for measurement of the cardiac-induced transbrachial impedance change, ΔZ(t), and the magnitudes and fiducial landmarks on its first time-derivative, transbrachial dZ/dt and on its second time-derivative, d2Z/dt2.
      • 3. The cumbersome and user-unfriendly transthoracic band, or tetrapolar spot-electrode array, is replaced by a circumferential or non-circumferential arm band or bands, an adhesive strip or other appropriate means for positioning the electrodes near the brachial artery containing a bipolar, or alternatively, a tetrapolar spot (or band) electrode array positioned on the medial aspect of the brachium between the axilla (arm pit) and a point distal on the brachium at the level of the olecranon process (elbow).
      • 4. With the arm at rest, motion artifacts are minimized as compared to the transthoracic approach, and thus, adaptive filtering techniques are less critical.
      • 5. Long-term monitoring of SV/CO in the surgical operating room, or intensive care unit, is facilitated by application of the apparatus to the arm, containing the bipolar, or, alternatively, the tetrapolar montage.
      • 6. The bioimpedance signal obtained from the brachium is unaffected by the presence of chest thoracostomy tubes, external pacemaker wires, surgical bandages or appliances, and percutaneously placed central venous access catheters located in the neck or upper chest.
      • 7. Without the perturbing influence of pulmonary ventilation, and pulmonary artery and other intrathoracic large vessel venous pulsations, the signal to noise ratio (S/N) relating to those portions the transbrachial dZ/dt and d2Z/dt2 signals pertaining only to left ventricular ejection are substantially higher than that of the transthoracic approach.
  • As disclosed above, the present invention relates to the measurement of stroke volume (SV) and cardiac output (CO) from the transbrachial method, using the brachial artery as the cardiogenically-induced signal source. Methodologically, the transbrachial method is similar to the transthoracic technique for determining SV. However, in the transthoracic technique, signal acquisition is effected over a segment of thorax by placement of voltage sensing electrodes 214 at the base of the neck, bilaterally, and voltage sensing electrodes 216 at the lower thorax at the xiphoid level, bilaterally (see FIG. 7). In contrast, the transbrachial technique uses a segment of brachium between voltage sensing electrodes 112/114 positioned proximate the axilla (arm pit) and junction of the upper and lower arm at the level of the olecranon process of the elbow. (see FIG. 1).
  • FIG. 1 schematically shows one apparatus embodiment according to the present invention, and its electrical interface with a subject 100. Signal acquisition from the upper arm 102 (brachium) requires application of a constant magnitude alternating current (A.C.) 104 of high frequency and small amplitude to electrodes 106, 108 that are spaced apart, with one or more electrodes affixed to the skin of the axilla, as well as one or more electrodes placed medially at the level of the antecubital fossa creating a current field. In this embodiment, the electrodes are applied to the subject's left arm. In other embodiments, the electrodes may be positioned on the right arm.
  • With the current field thus generated, the potential difference between the current injecting electrodes or alternating current flow electrodes 106, 108 is measured by a voltmeter 110 connected to the voltage sensing electrodes 112, 114 placed within the current field (see FIG. 1). A baseline impedance between the voltage sensing electrodes 112, 114, as well as a change in impedance, ΔZ(t) can be measured transbrachially. When the ΔZ(t) signal is electronically differentiated, dZ/dt(brachium) results, its peak systolic magnitude being dZ/dtmax(brachium). For the purposes of the invention disclosed herein, dZ/dtmax (Ω/s2) is equivalent to the nadir, or peak negative value of the rate of change of the impedance variation, dZ/dtmin. It is also understood and stipulated that −dZ/dtmax=+dZ/dtmax=dZ/dtmin. When dZ/dt undergoes electronic differentiation, d2Z/dt2(brachium) (Ω/s3) results, its peak magnitude being d2Z/dt2 max(brachium). For the purposes of the invention disclosed herein, d2Z/dt2 and d2Z/dt2 max are equivalent to the rate of change and peak negative rate of change of dZ/dt and dZ/dtmin, respectively. In the context of dZ/dtmax being equivalent to dZ/dtmin, the interpretation of the peak magnitudes of their respective derivatives should be understood and furthermore stipulated as equivalent (d2Z/dt2 max=d2Z/dt2 min). Many different methods of applying the electrodes or electrode arrays to the arm are envisioned, such as spot electrodes, arm band(s) both circumferential and non-circumferential, adhesive strips or other attachment means known in the art. In one embodiment, an 8 spot-electrode array can be implemented. Alternatively, in another embodiment, a 4 spot-electrode array, placed on the inner, or medial aspect of the upper arm, proximate the brachial artery, can be implemented. Alternatively, 4 non-circumferential band (strip) electrodes, embedded in an adhesive carrier, may be affixed to the brachium, medially, and used in lieu of spot electrodes.
  • The voltages measured by the Voltmeter 110 not only contains a signal caused by the AC applied, but may also contain a signal component from which an electrocardiogram (ECG) can be derived. The application of filters separates the AC related and ECG related signal components. In another embodiment, EKG 116 may also be measured by placing EKG electrodes 118 on the patient 100. In the figure, a 3-lead EKG is shown and EKG is measured by known means. The magnitude of the alternating current (A.C.) 104 and voltmeter 110 may be components of an apparatus 120. The apparatus 120 may also include an input device and a processor. The input device may be any suitable device that provides information to the apparatus, such as a keyboard. The input device may also receive information from external sources, such as the EKG 116. The processor is in communication with the data input device, the alternating current source 104 and electrodes 106, 108, and the voltmeter 110 and electrodes 112, 114. The processor is capable of receiving the information and calculating the stroke volume and cardiac output of the patient 100. The stroke volume and cardiac output of the patient may be displayed on a screen or be sent to other devices via a data output device of the apparatus.
  • FIG. 2 shows an example of ECG, ΔZ(t) and dZ/dt waveforms from a human subject 100, where TRR=the R-R interval, or the time for one cardiac cycle; Q=onset of ventricular depolarization; ---------=maximum systolic upslope extrapolation of ΔZ(t); B=aortic valve opening; C=peak rate of change of the thoracic cardiogenic impedance variation, dZ/dtmax; X=aortic valve closure; Y=pulmonic valve closure; O=rapid ventricular filling wave; Q-B interval=pre-ejection period, TPE; B-C interval=time-to-peak dZ/dt, TTP; B-X interval=left ventricular ejection period, TLVE. dZ/dt waveform to the right shows dZ/dtmax remaining constant throughout the ejection interval, TLVE.
  • Rationale for use of the brachium as an appropriate anatomic site for SV measurement by the bioimpedance technique is as follows. When A.C. (I) is injected through a segment of upper arm, otherwise known as the brachium, the boundaries of which are the deltoid muscles of the shoulder and axilla, proximally, and the elbow and antecubital fossa, distally, a quasi-static voltage, U0, and voltage change, ΔU(t), can be measured between the current injecting electrodes. Embedded within the brachial musculature and connective tissue, and anatomically situated medial to the brachial bone, otherwise known as the humerus, is the brachial artery. The brachial artery is a large artery, continuous with both the subclavian and axillary arteries, and, whereas the left subclavian artery is a major branch of the arch of the thoracic aorta, the right subclavian artery is a branch of the brachiocephalic artery. The contents of the upper arm, including connective tissue, bone, nervous tissue, veins, and the brachial artery, comprise an impedance (Z) to current flow. The passage of A.C. across the brachium generates a quasi-static voltage, U0, and, concordant with every pressure pulse of the brachial artery, a time-dependent drop in measured transbrachial voltage, ΔU(t), this following shortly after the onset of left ventricular ejection. The magnitude of the time delay (Δt, ms) between the brachial artery pressure pulse and the onset of left ventricular ejection is a function of pulse wave velocity. Transthoracically, the peak rate of change of impedance, dZ/dtmax, resulting from electronic differentiation of ΔZ(t), corresponds in time with peak aortic blood acceleration, dv/dtmax (cm/s2). Thus, in the preferred embodiment of the invention, dZ/dtmax(brachial) represents the ohmic analog of peak blood acceleration in the brachial artery. Chemla et al. (Fundam Clin Pharmacol 1996; 10:393-399) showed that the measured acceleration of blood in the brachial artery is highly correlated (r=0.79) and linearly proportional with blood acceleration in the ascending aorta. Moreover, whereas the magnitude of brachial artery blood velocity is affected by downstream peripheral vasoactivity (vasodilation, vasoconstriction), the magnitude of brachial artery blood acceleration is modulated only by beta (β) adrenergic stimulation or depression of the cardiac adrenoceptors (Chemla D, et al. Am J Cardiol 1990; 65:494-500). As extrapolated from Visser (Ann Biomed Eng 1989; 17:463-463), when flowing blood is interrogated by a field of alternating current (A.C.), the acceleration of blood in the aorta is measured as the aortic reduced average blood acceleration which is the mean aortic acceleration divided by the vessel radius: (dv/dt(mean)/R). When [(dv/dt(mean)/R]max, or peak aortic reduced average blood acceleration (1/s2), undergoes square root transformation, peak aortic reduced average blood velocity (1/s) results. Likewise, when dZ/dtmax/Z0 (1/s2) undergoes square root transformation, ohmic mean velocity results, ΔZv(t)max/Z0 (1/s). In the context of the present invention, where dZ/dtmax represents ohmic mean acceleration, dZ/dtmax/Z0 is herein referred to as the acceleration index (ACI). Because of the high correlation of dv/dtmax measured in the aorta, with that of the brachial artery, it is claimed by the present invention that SV can also be obtained from the brachial artery. It has been discovered that processing the mean value of the second time-derivative of the impedance change provides results equivalent to processing the peak value of the first time-derivative, but with better accuracy. Thus, in the preferred embodiment, SV TB = V C · [ ( 2 Z / t mean 2 Z 0 ) · 10 - 2 ] 3 · T LVE equation 1
    where SVTB=transbrachial SV (mL); Vc(brachium) equals the volume conductor (mL); Z0 equals the quasi-static transbrachial base impedance (Ohm, Ω); Tlve equals left ventricular ejection time (s), and d2Z/dt2 mean equals the mean value of the second time-derivative of the cardiogenically induced transbrachial impedance variation (Ω/s3). In a first embodiment:
    d 2 Z/dt 2 mean=(dZ/dt max)/TTP m,  equation 1a
    where TTPm=measured time to peak dZ/dt, defined as the temporal interval from point B to point C (seconds, s) in FIG. 2. In a second embodiment:
    d 2 Z/dt 2 mean=(dZ/dt max)/TTP m,  equation 1b
    where 0.01 s≦TTP≦0.1 s, such as TTPb=0.06 s. In a third embodiment:
    d 2 Z/dt 2 mean=(dZ/dt max)/TTP c,  equation 1c
    where TTPc=corrected rise time, or time to peak dZ/dt, which is the corrected temporal interval from point B to point C (s). Here, TTPc=(10−a)×[1/(ACI)b=s; where [1/(ACI)b]=ACIb; “a” is a negative exponent such as 0>a≧5 or even a=−2; ACI=acceleration index=dZ/dtmax/Z0=1/s2; “b” is an exponent such as 0.1≦−b≦1.0 or even b=0.5; where −b (minus b) is a negative exponent such as 0>−b≧10 or even −b (minus b)=−0.5; and 0.01 s≦TTPc≦0.1 s. For example, in a preferred embodiment, TTPc=(10−2)×[1/(ACI)0.5]=s. In a fourth embodiment:
    d 2 Z/dt 2 mean =dZ/dt max /TTP d,  equation 1d
    where TTPd(s)=(TTPc+TTPm)/A; and 0<A≦5 such as A=2; and 0.01 s≦TTPd≦0.1 s. In a fifth embodiment:
    d 2 Z/dt 2 mean =dZ/dt max /TTP e,  equation 1e
    where TTPc(s)=(TTPc+TTPb)/A; and 0<A≦5 such as A=2; and 0.01 s≦TTPc≦0.1 s. In a sixth embodiment:
    d 2 Z/dt 2 mean =dZ/dt max /TTP f,  equation 1f
    where TTPf=(TTPm+TTPb)/A; and 0<A≦5 such as A=2; and 0.01 s≦TTPf≦0.1 s. In a seventh embodiment:
    d 2 Z/dt 2 mean =d 2 Z/dt 2 max /B;  equation 1g
    where d2Z/dt2 max (Ω/s3)=the maximal systolic upslope extrapolation of dZ/dt, peak first time-derivative of dZ/dt, or systolic peak rate of change of dZ/dt; and 0<B≦10 such as B=2.
  • In one exemplary embodiment:
      • Vc(brachium)=C1·[W·C2];
      • C1=0<C1≦50,000,
      • W=weight in kilograms (kg)
      • C2=C3/(BMIn)y
      • BMIn=BMIp/C4
        • a. 35≦C3≦100 (mL/kg), wherein C3 in the preferred embodiment=70 mL/kg;
        • b. BMIn=normalized body mass index (dimensionless), wherein 0.5≦BMIn≦1.0, wherein the preferred embodiment, BMIn=1
        • c. BMIp=a person's body mass index=Weight (kg)/Height (meters)2 (kg/m2) where W=a person's weight (kg), and H=a person's height (m).
        • d. 15≦C4≦40 kg/m2, wherein C4=ideal body mass index=24 kg/m2 in the preferred embodiment.
        • e. 0.25≦y≦1.0, wherein y=0.5 in the preferred embodiment.
        • f. 35≦C2≦100, wherein the preferred embodiment, C2=70 mL/kg
  • Unlike previously described bioimpedance techniques, which broadly assume a plethysmographic or volumetric origin for the cardiogenic impedance change, ΔZ(t), and its peak first time derivative, dZ(t)/dtmax, the present technique assumes dZ(t)/dtmax to represent the ohmic equivalent of the peak acceleration of red blood cells. Thus, when the first time-derivative of ΔZ(t) is taken, dZ/dt, its peak magnitude, dZ/dtmax, can be shown to coincide in time with the peak red blood cell acceleration, dv/dtmax (cm/s2), and not with the peak rate of change of volume, dV/dtmax (mL/s). Consequently, to obtain ohmic mean velocity, dZ/dtmax/Z0 (1/s2) must undergo square root transformation. This transformation is to be known as square root Acceleration Step-down Transformation: √{square root over ([)}(dZ/dtmax)/Z0] (1/s). It can also be shown that the cube root of the normalized mean value of the second time-derivative of the impedance change times 10−2, that is, [(d2Z/dt2 mean/Z0)×10−2]0.333, yields ohmic mean velocity equivalent to the taking the square root of the peak value of the first time-derivative. Signal processing the mean value of the second time-derivative of ΔZ(t) has certain advantages over processing the first time-derivative (e.g. time to peak values calculated more accurately because of using cube root calculations over square root calculations), which are provided in another preferred embodiment using the transthoracic approach. To obtain brachial artery ohmic mean velocity, the cube root transformation is thus implemented in the preferred embodiment. Thus, the signal processing technique, comprising part of the invention, implies that the proper designation for the transbrachial approach is correctly stated as Transbrachial Electrical Bioimpedance Cardiovelocimetry or, simply, Transbrachial bioimpedance velocimetry.
  • Impedance Measurement Techniques in the Preferred Embodiment of the Invention
  • FIGS. 3 a and 3 b show the relationship between the dZ/dt curve and either the dP/dt or d(SpO2)/dt curves. FIG. 3 a further shows an example where points B and X are apparent on the dZ/dt curve and FIG. 3 b shows an example where point B is not detectable, but point X is detectable on the dZ/dt curve. Thus, determination of left ventricular ejection time, (Tlve), onset of flow (point B), and the ohmic equivalent of peak brachial artery reduced average blood acceleration, (dZ/dtmax/Z0 (brachium)), while ideally measured directly from the dZ/dt curve, are supplemented obligatorily by alternative means. The said alternative/obligatory means for determining Tlve are those obtained from means such as from the waveform corresponding to the photoplethysmographic pulse oximetry waveform, ΔSpO2(t), or its first time-derivative, d(SpO2)/dt, and/or by the waveform obtained from a non-invasive applanated radial arterial pressure pulse waveform, ΔP(t)(radial), or its first time-derivative, dP/dt(radial). (see FIG. 4)
  • The said means for determining point B on the transbrachial dZ/dt curve are those methods used for determining Tlve when point X on the transbrachial dZ/dt curve, or its first time derivative (d2Z/dt2), are identifiable by those skilled in the art of bioimpedance curve analysis. When point X is not identifiable on the transbrachial dZ/dt curve, or its first time-derivative, d2Z/dt2, then alternative means for point B detection are necessary. In the absence of an identifiable point X by those skilled in the art of dZ/dt curve analysis, said means for point B detection include use of the first time-derivative of the applanated radial pressure waveform tracing, dP/dt(radial). It will be clear to those skilled in the art of curve analysis, why the aforementioned said means are superior to those disclosed by others, and most recently by Baura et al. (U.S. Pat. No. 6,561,986 B2).
  • Methods for Determination of Left Ventricular Ejection Time (Tlve, LVET):
      • 1. dZ/dt waveform analysis: Tlve (LVET) measured across the brachium is defined as the temporal interval from point B, which corresponds to aortic valve opening (AVO), albeit with a time delay, to point X, which coincides in time, albeit with a short time delay, to aortic valve closure (AVC), these time delays dictated by pulse wave velocity.
      • 2. d2Z/dt2 waveform analysis: TLVE (LVET) measured from the peak second time-derivative of the transbrachial or transthoracic impedance pulse variation is defined as the temporal interval proximate the peak of the earliest positive systolic deflection of d2Z/dt2 (i.e. d2Z/dt2 max), corresponding proximate in time with point B on the dZ/dt curve and the nadir of the rising foot of ΔZ(t), to, in the usual case, the second zero baseline impedance crossing of d2Z/dt2. To the ejection interval just described, 20 ms (i.e. 0.02 seconds) should be added. The second Zero crossing, following the first zero crossing corresponding to point C (dZ/dtmax), corresponds temporally with point X on the dZ/dt curve. These relationships are better appreciated by inspection of FIG. 8.
      • 3. Pulse Oximetry waveform (ΔSpO2(t)): LVET is defined as the temporal interval (seconds) from the onset of the oximetric pulse at zero baseline, signifying the onset of ejection, albeit with a time delay, to the oximetry wave equivalent of the dicrotic notch, which signifies aortic valve closure, albeit with a time delay, and the end of ejection. The oximetry waveform can be obtained from any appropriate site on, or within the human body, but, in the preferred embodiment, the distal digit of the human finger is deemed most appropriate.
      • 4. Applanation Tonometry Pressure Pulse waveform (ΔP(t)(radial)): LVET is defined as the temporal interval (seconds) from the onset of the pressure pulse at zero baseline, signifying the onset of ejection, albeit with a time delay, to the dicrotic notch equivalent, which signifies aortic valve closure, albeit with a time delay, and the end of ejection. In the preferred embodiment, the pressure pulse waveform is obtained from the radial artery at the wrist, but may be obtained from any site on the arm, specifically from either brachial artery.
      • 5. Regression Equations for Tlve versus Heart Rate (HR): LVET is determined by Weissler's regression equations: Male: Tlve=−0.0017·HR+0.413; Female: Tlve=−0.0016·HR+0.418.
        It should be noted that any of the above methods of LVET determinations can be used solely, or in combination with each other (e.g. average multiple LVET determinations, use one or more determinations that appear to produce better results, etc.).
  • With exemplary dZ/dt waveforms, such as those shown in FIG. 3 a, point B and point X are readily distinguishable by one skilled in the art of curve analysis. However, these fiducial landmarks are frequently distorted by motion and ventilation artifacts (especially using the transthoracic approach), as well as by certain disease processes. LVET may be more accurately measured by curve analysis of the pulse oximetry and applanation tonometry waveforms (FIG. 3 b), or their first time-derivatives. In one embodiment of the invention, either or both methods may be implemented. Of these techniques, applanation tonometry is most likely to demonstrate a dicrotic notch, and, therefore, is considered the preferred technique. Furthermore, for those skilled in the art of computer waveform analysis, the points coinciding with the beginning and end of ejection can be readily identified from the first time-derivative curves of both the oximetry and applanation tonometry waveforms; namely, d(SpO2)/dt and dP/dt(radial). In the preferred embodiment, the best method constitutes computer analysis of the first time-derivatives. In another embodiment, regression equations for heart rate versus LVET may be implemented.
  • Methods for Point B Detection on the dZ/dt Waveform:
      • 1. Methods for determining point B when point X is readily identifiable by one skilled in the art of curve analysis (see FIG. 3 a).
  • Point B on the transbrachial dZ/dt waveform is known to coincide with aortic valve opening, albeit with a time delay. Exemplary dZ/dt waveforms demonstrate a distinct change in slope at, or not uncommonly above the zero baseline impedance, followed by a steep, positive linear segment ending at point C, or dZ/dtmax. When a distinct change in slope leading to point C is detected at or above the baseline, one skilled in the art of curve analysis can readily identify point B. However, as demonstrated by Debski T T et al. (Biol Psychol 1993; 36:63-74) using the transthoracic method, despite using fiducial landmarks on the time-derivatives of dZ/dt (i.e. d2Z/dt2 and d3Z/dt3) to identify this change in slope, detection of point B can be problematic. This inability to correctly identify point B is obvious to those skilled in the art of curve analysis, and especially curve analysis of dZ/dt, by inspection of FIG. 3 b. The method disclosed herein provides a new and innovative solution for point B detection. The new method employs one, or a combination of methods disclosed under determination of LVET; namely, ΔSpO2(t) and/or ΔP(t)(radial) (as shown in FIG. 4), or, respectively, their time derivatives, d(SpO2)/dt and/or dP/dt(radial) (as shown in FIG. 5). The technique of point B detection, as disclosed herein as a preferred embodiment, involves computerized curve fitting and alignment in time of temporal landmark X on the transbrachial dZ/dt curve with the dicrotic notch equivalent of one or both of the measured aforementioned oximetry and pressure curves, and/or preferably with one or both of their first time-derivatives. Independently, or in concert, one or both of the first time-derivative curves can be aligned in time with the transbrachial dZ/dt curve, such that the temporal point of the termination of flow, or aortic valve closure (AVC) equivalent on the first derivative oximetry or pressure curves, can be aligned in time with point X of the transbrachial dZ/dt curve. Point B, coinciding with aortic valve opening (AVO), and the beginning of flow, albeit with a time delay, is identified by determining the temporal point on the transbrachial dZ/dt curve, intersecting, and coinciding in time with the point of onset of flow/pressure on the ΔSpO2(t)/AP(t) curves, and/or on their first time-derivatives. This temporal point is identified as a discreet point at the baseline occurring before the first positive maximum upslope measured from foot of the respective baselines of the ΔSpO2(t) and/or ΔP(t)(radial) curves, and/or from their first time-derivatives (FIGS. 3 a, 3 b, 4, 5, 6).
      • 2. Method for determining point B when point X is not readily identifiable by one skilled in the art of curve analysis (see FIG. 6):
  • When point X is not readily identifiable by one skilled in the art of curve analysis, then alternative means must be applied. Requiring alternative means, for example, would be the inability to identify the first zero crossing at baseline impedance after the zero crossing of point C (dZ/dtmax) on the second time-derivative curve of ΔZ(t)(i.e., d2Z/dt2), where said zero crossing corresponds in time to point X and AVC. Said alternative means requires application of the first time-derivative of the applanation tonometry curve, dP/dt(radial) (FIG. 6). For one skilled in the art of curve analysis, said means requires alignment in time of the earliest maximum positive peak of dP/dt (dP/dtmax) with point C of the transbrachial dZ/dt curve. With point dP/dtmax and point C aligned in time, point B can be identified by applying a perpendicular through, and coinciding in time with the onset at baseline of the first positive deflection of dP/dt, where said perpendicular line must intersect the dZ/dt curve at or above baseline impedance. The point of intersection of the perpendicular with the transbrachial dZ/dt curve is designated point B. When the above methods are unavailable, or fail to supply waveforms with fiducial landmarks necessary for point B detection, as assessed by pre-determined criteria, then, as default methods, a point on the transbrachial dZ/dt curve occurring 55 ms prior to point C, but obligatorily at or above baseline impedance, or alternatively, a point 15% above baseline impedance on the dZ/dt curve, is taken as point B.
  • Method for determining the maximum systolic upslope of transbrachial dZ/dt, otherwise known as transbrachial d2Z/dt2 max: Employing one or a combination of the techniques described herein for point B detection, transbrachial d2Z/dt2 max is the measured peak positive deflection of the d2Z/dt2 (d2Z/dt2 max) curve usually occurring temporally proximate point B on the transbrachial dZ/dt curve.
  • In one embodiment, external calibration of the SV/CO by means of the transbrachial approach:
      • 1. External calibration of Vc(brachium) by means of the transthoracic method: Determination of Vc(cal).
  • Because of the high correlation of dv/dtmax measured in the aorta with that of the brachial artery, it is claimed that: V C ( thorax ) = [ ( 2 Z / t mean 2 Z 0 ) thorax · 10 - 2 ] 3 · T LVE = V C ( brachium ) · [ ( 2 Z / t mean 2 Z 0 ) transbrachial · 10 - 2 ] 3 · T LVE ( equation 2 )
  • Since Tlve is equivalent for both sides of equation 2, Vc(brachium) can be found thusly, V C ( cal ) = V C ( brachium ) = V C ( thorax ) [ ( 2 Z / t mean 2 Z 0 ) thorax · 10 - 2 ] 3 [ ( 2 Z / t mean 2 Z 0 ) transbrachial · 10 - 2 ] 3 equation 3 Where V c ( thorax ) = C 1 ( thorax ) · [ W ( kg ) · C 2 ] equation 4
    Where 0.10≦C1(thorax)≦0.75, wherein C1=0.25 in the preferred embodiment,
  • Thus, SV by the transbrachial method, externally calibrated from the transthoracic approach is given as, SV transbrachial = V C ( cal ) · [ ( 2 Z / t mean 2 Z 0 ) transbrachial · 10 - 2 ] 3 · T LVE equation 5
  • 2. Determination of SV from the transbrachial approach by means of auto-calibration: A priori determination of C1(brachium) as a mean value for a population, n.
  • In order to satisfy the requirements of equation 2, Vc(brachium) is found by determining Vc(cal) from the solution of equation 3. This operation requires insertion of Vc(thorax) as determined from equation 4. Therefore, Vc(cal) in equation 3 can be determined thusly;
    V c(cal) =C 1(brachium) ·[W(kg)·C 2]  equation 6
    Where, C1 is thus,
    C 1(brachium) =V c(cal) /[W(kg)·C 2]  equation 7
    where, 0≦C1≦50,000, wherein the preferred embodiment C1 is proprietary.
  • By solving equation 7 for a population, n, determining Vc(cal) from equation 3, the mean value of the constant, C1(brachium), can be found for the general population as follows;
    C 1(brachium)(mean)=[(C 1-1 +C 1-2+C1-4+ . . . C1-n)/n]  equation 8
    Where C1(mean) ideally=C1-1 through C1-n. Thus, SV determination by the transbrachial approach by auto-calibration is given as, SV transbrachial = [ C 1 ( mean , brachial ) · W · C 2 ] · [ ( 2 Z / t mean 2 Z 0 ) transbrachial · 10 - 2 ] 3 · T LVE equation 9
  • In another embodiment of the present invention, the SV of the left ventricle can be determined using the transthoracic approach and the mean value of the second time-derivative of the cardiogenically induced transthoracic impedance variation (Ω/s3). As implemented by inspection of FIG. 7 and related description, a tetrapolar spot electrode array can be applied to a person's body. The description of signal acquisition and processing is precisely that described above with respect to FIG. 2. Stroke volume determination by means of the transthoracic (TT) application is implemented by means of the following equation (where, in the general embodiment, the stroke volume (SV) equation is given as): SV TT = V C ( TT ) · [ ( 2 Z / t mean 2 Z 0 ) TT · 10 - 2 ] 3 · T LVE equation 10
    Where:
      • SV(TT)=stroke volume by the transthoracic approach (mL).
      • VC(TT)=Volume conductor (mL), otherwise known as the volume of electrically participating thoracic tissue (VEPT). Ohmic mean velocity ( 1 / sec ) = [ ( 2 Z / t mean 2 Z 0 ) TT · 10 - 2 ] 3 equation 11
  • d2Z/dt2 mean is determined identically as described above with regard to equations 1 and 1a-1g.
  • TLVE=left ventricular ejection time (seconds, sec, s).
  • Evaluation, ranges and preferred embodiments of the input variables of equation 10 are as follows: V C ( TT ) = ζ · ( C 2 · W · C 3 ) ζ ( zeta ) = ( C 1 - C 1 Z 0 Z c + Z 0 2 Z c 2 ) ( dimensionless ) equation 12
      • 1.0≦C1≦8.0, wherein the preferred embodiment C1=4.
      • Z0 (Ohm, Ω)=the transthoracic base impedance, or static D.C. component of the total transthoracic impedance, Z(t).
      • Zc (Ω)=the critical level of base impedance; wherein, 15 Ω≦Zc≦25 Ω, and wherein the preferred embodiment, Zc=20 Ω. For all values of Z0<20 Ω, ζ>1.0, and for all values of Z0<20 Ω, ζ=1.0.
      • 0.10≦C2≦0.75, wherein the preferred embodiment, C2=0.25
      • W=weight in kilograms (kg) C 3 = C 4 ( BMI N ) y = mL / kg equation 13
      • 35≦C4≦100=mL/kg, wherein the preferred embodiment, C4=70 mL/kg.
      • BMIN=normalized body mass index (dimensionless), where 0.5≦BMIN≦5.0, wherein the preferred embodiment BMIN=1.0.
      • BMIN=BMIp/C5
  • BMIP=a person's body mass index (W(kg)/H(m2), i.e. kg/m2), where W=a person's weight in kilograms (kg), and H=a person's height in meters (m).
      • C5=a person's ideal body mass index (kg/m2), where 10≦C5≦100, wherein the preferred embodiment, C5=24 kg/m2.
      • y is an exponent, 0.25≦y≦1.0, wherein the preferred embodiment, y=0.5, where 1/y=equivalent root function, where, 1/y√{square root over (BMIN)}=(BMIN)y.
      • 35≦C3≦100, wherein the preferred embodiment, C3=70 mL/kg. V C ( TT ) ( mL ) = [ ( C 1 - C 1 Z 0 Z c + Z 0 2 Z c 2 ) · ( C 2 · W · C 3 ) ] equation 14
  • The preferred embodiment the product of C2 and C3 is a constant K, (i.e. C2·C3=K), where 10≦C2·C3≦35, wherein the preferred embodiment, C2·C3=17.5 mL/kg
  • Where equation 11 in the broadest definition is given as, Ohmic mean velocity ( 1 / sec ) = [ ( n 2 Z / t mean n 2 Z 0 ) TT · 10 - n 3 ] n 1 = [ ( n 2 Z / t mean n 2 Z 0 ) TT · 10 - n 3 ] n 4 equation 15
    Where:
      • 1.0≦n1≦10, wherein the preferred embodiment, n1=3.
      • 1.0≦n2≦10, wherein the preferred embodiment, n2=2.
      • 0<n3≦5, wherein the preferred embodiment, n3=2.
      • 0.1≦n4≦1.0, wherein the preferred embodiment, n4=0.333. n 4 = 1 n 1
      • TLVE=left ventricular ejection time (seconds, sec).
  • Where the stroke volume equation for the transthoracic application in its broadest definition is given as, SV TT = [ ( C 1 - C 1 Z 0 Z c + Z 0 2 Z c 2 ) · ( C 2 · W · C 3 ) ] · [ ( n 2 Z / t mean n 2 Z 0 ) · 10 - n 3 ] n 1 · T LVE equation 16
    Wherein as operationally implemented in the preferred embodiment, the stroke volume equation for the transthoracic application is given as, SV TT = [ ( C 1 - C 1 Z 0 Z c + Z 0 2 Z c 2 ) · ( C 2 · W · C 3 ) ] · [ ( 2 Z / t mean 2 Z 0 ) TT · 10 - 2 ] 3 · T LVE equation 17
  • It is to be understood that the present invention is not limited to the embodiment(s) described above and illustrated herein, but encompasses any and all variations falling within the scope of the appended claims.

Claims (34)

1. An apparatus for determining stroke volume by bioimpedance from a patient, comprising:
two or more spaced apart alternating current flow electrodes positionable on a patient;
two or more spaced apart voltage sensing electrodes positionable on the patient and between the alternating current flow electrodes;
an alternating current source electrically connected to the alternating current flow electrodes;
a voltmeter electrically connected to the voltage sensing electrodes; and
a processing unit in communication with the voltage sensing electrodes, wherein the processing unit is capable of using a voltage sensed by the voltage sensing electrodes to calculate a mean value of a second time-derivative of a cardiogenically induced impedance variation of the patient and to determine therefrom a stroke volume of the patient.
2. The apparatus of claim 1, wherein the alternating current flow electrodes and the voltage sensing electrodes are positionable on an arm of the patient proximate a brachial artery.
3. The apparatus of claim 1, wherein the alternating current flow electrodes and the voltage sensing electrodes are positionable on a base of the patient's neck and on a lower thorax of the patient.
4. The apparatus of claim 1, wherein the processing unit calculates the mean value of the second time-derivative of the cardiogenically induced impedance variation as a peak rate of change of the cardiogenically induced impedance variation divided by a time to peak value of the cardiogenically induced impedance.
5. The apparatus of claim 4, wherein the time to peak value is a measured time period between an aortic valve of the patient opening and a peak rate of change of the cardiogenically induced impedance variation.
6. The apparatus of claim 4, wherein the time to peak value is a selected time value at least as great as 0.01 seconds and no greater than 0.1 seconds.
7. The apparatus of claim 4, wherein the time to peak value is a corrected rise time between an aortic valve of the patient opening and a peak rate of change of the cardiogenically induced impedance variation.
8. The apparatus of claim 4, wherein the time to peak value is calculated as (TTPc+TTPm)/A, where TTPm is a measured time period between an aortic valve of the patient opening and a peak rate of change of the cardiogenically induced impedance variation, TTPc is a corrected rise time between an aortic valve of the patient opening and a peak rate of change of the cardiogenically induced impedance variation, and A is a value greater than zero but not greater than 5.
9. The apparatus of claim 4, wherein the time to peak value is calculated as (TTPc+TTPb)/A, where TTPb is a selected time value at least as great as 0.01 seconds and no greater than 0.1 seconds, and TTPc is a corrected rise time between an aortic valve of the patient opening and a peak rate of change of the cardiogenically induced impedance variation, and A is a value greater than zero but not greater than 5.
10. The apparatus of claim 4, wherein the time to peak value is calculated as (TTPm+TTPb)/A, where TTPb is a selected time value at least as great as 0.01 seconds and no greater than 0.1 seconds, and TTPm is a measured time period between an aortic valve of the patient opening and a peak rate of change of the cardiogenically induced impedance variation, and A is a value greater than zero but not greater than 5.
11. The apparatus of claim 1, wherein the processing unit calculates the mean value of the second time-derivative of the cardiogenically induced impedance variation as a maximum value of the second time-derivative of the cardiogenically induced impedance variation divided by value that is greater than zero but no greater than 10.
12. The apparatus of claim 1, wherein the processing unit uses the equation:
SV = V C · [ ( 2 Z / t mean 2 Z 0 ) · 10 - 2 ] 3 · T LVE
to determine the stroke volume of the patient, wherein Vc is a volume of electrically participating tissue of the patient, d2Z/dt2 mean is a mean value of a second time-derivative of a cardiogenically induced impedance variation of the patient, Z0 is a quasi-static base impedance, and TLVE is a left ventricular ejection time of the patient.
13. The apparatus of claim 12, wherein TLVE is obtained from a dZ/dt waveform.
14. The apparatus of claim 13, wherein a trigger for initiating processing of the dZ/dt waveform is obtained from an R wave of an antecedent ECG waveform or a C wave of an antecedent dZ/dt waveform.
15. The apparatus of claim 12, wherein TLVE is obtained from a pulse oximetry waveform (ΔSpO2(t)), or its first time-derivative, dSpO2(t)/dt.
16. The apparatus of claim 12, wherein TLVE is obtained from an applanation tonometry (pressure) waveform (ΔP(t)), or its first time-derivative dP(t)/dt.
17. The apparatus of claim 12, wherein TLVE is obtained from regression equations.
18. A method of determining stroke volume by bioimpedance from a patient, comprising:
positioning two or more spaced apart alternating current flow electrodes on a patient;
positioning two or more spaced apart voltage sensing electrodes on the patient and between the alternating current flow electrodes;
providing an alternating current flow (I(t)) through the electrically conductive electrodes creating a current field;
measuring a voltage (U(t)) between the voltage sensing electrodes within the current field; and
calculating a mean value of a second time-derivative of a cardiogenically induced impedance variation of the patient using the measured voltage (U(t)); and
calculating the stroke volume (SV) of the patient using the calculated second time-derivative mean value.
19. The method of claim 18, wherein the positioning of the alternating current flow electrodes includes positioning the alternating current flow electrodes on an arm of the patient proximate a brachial artery, and wherein the positioning of the voltage sensing electrodes including positioning the voltage sensing electrodes on the arm of the patient proximate the brachial artery.
20. The method of claim 18, wherein the positioning of the alternating current flow electrodes includes positioning the alternating current flow electrodes on a base of the patient's neck and on a lower thorax of the patient, and wherein the positioning of the voltage sensing electrodes including positioning the voltage sensing electrodes on the base of the patient's neck and on the lower thorax of the patient
21. The method of claim 18, wherein the calculating of the second time-derivative mean value includes dividing a peak rate of change of the cardiogenically induced impedance variation by a time to peak value of the cardiogenically induced impedance.
22. The method of claim 21, wherein the time to peak value is determined by measuring a time period between an aortic valve of the patient opening and a peak rate of change of the cardiogenically induced impedance variation.
23. The method of claim 21, wherein the time to peak value is determined by selecting a time value at least as great as 0.01 seconds and no greater than 0.1 seconds.
24. The method of claim 21, wherein the time to peak value is determined by correcting a rise time between an aortic valve of the patient opening and a peak rate of change of the cardiogenically induced impedance variation.
25. The method of claim 21, wherein the time to peak value is determined by calculating (TTPc+TTPm)/A, where TTPm is a measured time period between an aortic valve of the patient opening and a peak rate of change of the cardiogenically induced impedance variation, TTPc is a corrected rise time between an aortic valve of the patient opening and a peak rate of change of the cardiogenically induced impedance variation, and A is a value greater than zero but not greater than 5.
26. The method of claim 21, wherein the time to peak value is determined by calculating (TTPc+TTPb)/A, where TTPb is a selected time value at least as great as 0.01 seconds and no greater than 0.1 seconds, and TTPc is a corrected rise time between an aortic valve of the patient opening and a peak rate of change of the cardiogenically induced impedance variation, and A is a value greater than zero but not greater than 5.
27. The method of claim 21, wherein the time to peak value is determined by calculating (TTPm+TTPb)/A, where TTPb is a selected time value at least as great as 0.01 seconds and no greater than 0.1 seconds, and TTPm is a measured time period between an aortic valve of the patient opening and a peak rate of change of the cardiogenically induced impedance variation, and A is a value greater than zero but not greater than 5.
28. The method of claim 18, wherein mean value of the second time-derivative of the cardiogenically induced impedance variation is calculated as a maximum value of the second time-derivative of the cardiogenically induced impedance variation divided by value that is greater than zero but no greater than 10.
29. The method of claim 18, wherein the calculating the stroke volume (SV) includes:
determining a volume of electrically participating tissue Vc of the patient;
determining a quasi-static base impedance Z0;
determining a left ventricular ejection time TLVE of the patient; and
calculating the stroke volume (SV) of the patient using the equation:
SV = V C · [ ( 2 Z / t mean 2 Z 0 ) · 10 - 2 ] 3 · T LVE
wherein d2Z/dt2 mean is the second time-derivative mean value.
30. The method of claim 29, wherein TLVE is determined from a dZ/dt waveform.
31. The method of claim 30, wherein a trigger for initiating a processing of the dZ/dt waveform is obtained from an R wave or a C wave of an antecedent ECG waveform.
32. The method of claim 29, wherein TLVE is determined from a pulse oximetry waveform (ΔSpO2(t)), or its first time-derivative, dSpO2(t)/dt.
33. The method of claim 29, wherein TLVE is determined from an applanation tonometry (pressure) waveform (ΔP(t)), or its first time-derivative dP(t)/dt.
34. The method of claim 29, wherein TLVE is determined from regression equations.
US11/158,521 2004-06-16 2005-06-21 Apparatus and method for determination of stroke volume using the brachial artery Active - Reinstated 2028-06-19 US7806830B2 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US11/158,521 US7806830B2 (en) 2004-06-16 2005-06-21 Apparatus and method for determination of stroke volume using the brachial artery
PCT/US2005/044671 WO2006063255A2 (en) 2004-12-08 2005-12-08 Determining of stroke volume using the brachial artery
US12/489,335 US7740590B2 (en) 2004-06-16 2009-06-22 Apparatus and method for determination of stroke volume using the brachial artery

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US10/870,281 US7261697B2 (en) 2004-06-16 2004-06-16 Apparatus for determination of stroke volume using the brachial artery
US63461604P 2004-12-08 2004-12-08
US11/158,521 US7806830B2 (en) 2004-06-16 2005-06-21 Apparatus and method for determination of stroke volume using the brachial artery

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US10/870,281 Continuation-In-Part US7261697B2 (en) 2004-06-16 2004-06-16 Apparatus for determination of stroke volume using the brachial artery

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US12/489,335 Division US7740590B2 (en) 2004-06-16 2009-06-22 Apparatus and method for determination of stroke volume using the brachial artery

Publications (2)

Publication Number Publication Date
US20060009710A1 true US20060009710A1 (en) 2006-01-12
US7806830B2 US7806830B2 (en) 2010-10-05

Family

ID=36578625

Family Applications (2)

Application Number Title Priority Date Filing Date
US11/158,521 Active - Reinstated 2028-06-19 US7806830B2 (en) 2004-06-16 2005-06-21 Apparatus and method for determination of stroke volume using the brachial artery
US12/489,335 Expired - Lifetime US7740590B2 (en) 2004-06-16 2009-06-22 Apparatus and method for determination of stroke volume using the brachial artery

Family Applications After (1)

Application Number Title Priority Date Filing Date
US12/489,335 Expired - Lifetime US7740590B2 (en) 2004-06-16 2009-06-22 Apparatus and method for determination of stroke volume using the brachial artery

Country Status (2)

Country Link
US (2) US7806830B2 (en)
WO (1) WO2006063255A2 (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100081960A1 (en) * 2008-09-30 2010-04-01 Nellcor Puritan Bennett Llc Bioimpedance System and Sensor and Technique for Using the Same
US20140249440A1 (en) * 2010-12-28 2014-09-04 Matt Banet Body-worn system for continuous, noninvasive measurement of cardiac output, stroke volume, cardiac power, and blood pressure
US20150324977A1 (en) * 2014-05-09 2015-11-12 Siemens Medical Solutions Usa, Inc. Image-Based Waveform Parameter Estimation
US9451888B1 (en) * 2015-03-27 2016-09-27 Cordeus, Inc. Method and apparatus for determination of left ventricular stroke volume and cardiac output using the arteries of the forearm
US20190046069A1 (en) * 2015-07-10 2019-02-14 Bodyport Inc. Cardiovascular signal acquisition, fusion, and noise mitigation
EP2737847B1 (en) * 2011-07-29 2020-02-26 Universitat Politècnica De Catalunya Method and apparatus for obtaining cardiovascular information by measuring between two extremities
WO2022182975A1 (en) * 2021-02-26 2022-09-01 Massachusetts Institute Of Technology Methods to simulate metrics of vascular function from clinical data
US11696715B2 (en) 2015-07-10 2023-07-11 Bodyport Inc. Cardiovascular signal acquisition, fusion, and noise mitigation
CN117281494A (en) * 2023-11-27 2023-12-26 安徽通灵仿生科技有限公司 Method and device for identifying signal characteristic points of arterial blood pressure signals

Families Citing this family (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6850788B2 (en) 2002-03-25 2005-02-01 Masimo Corporation Physiological measurement communications adapter
US8840549B2 (en) 2006-09-22 2014-09-23 Masimo Corporation Modular patient monitor
US8241222B2 (en) 2008-07-31 2012-08-14 Medtronic, Inc. Monitoring hemodynamic status based on intracardiac or vascular impedance
US20110208015A1 (en) 2009-07-20 2011-08-25 Masimo Corporation Wireless patient monitoring system
US9153112B1 (en) 2009-12-21 2015-10-06 Masimo Corporation Modular patient monitor
DE102010016172A1 (en) 2010-03-26 2011-09-29 Medis Medizinische Messtechnik Gmbh Arrangement and method for non-invasive detection of haemodynamic parameters
US9084585B1 (en) * 2011-10-06 2015-07-21 Cerner Innovation, Inc. System and method to classify left ventricular and mitral valve function
WO2013056160A2 (en) 2011-10-13 2013-04-18 Masimo Corporation Medical monitoring hub
US9943269B2 (en) 2011-10-13 2018-04-17 Masimo Corporation System for displaying medical monitoring data
WO2013081586A1 (en) * 2011-11-29 2013-06-06 King Saud University Systems and methods to measure fluid in a body segment
US9895069B2 (en) 2011-11-29 2018-02-20 King Saud University Systems and methods to measure fluid in a body segment
US10149616B2 (en) 2012-02-09 2018-12-11 Masimo Corporation Wireless patient monitoring device
US9060745B2 (en) 2012-08-22 2015-06-23 Covidien Lp System and method for detecting fluid responsiveness of a patient
US8731649B2 (en) 2012-08-30 2014-05-20 Covidien Lp Systems and methods for analyzing changes in cardiac output
US9357937B2 (en) 2012-09-06 2016-06-07 Covidien Lp System and method for determining stroke volume of an individual
US9241646B2 (en) 2012-09-11 2016-01-26 Covidien Lp System and method for determining stroke volume of a patient
US20140081152A1 (en) 2012-09-14 2014-03-20 Nellcor Puritan Bennett Llc System and method for determining stability of cardiac output
US9749232B2 (en) 2012-09-20 2017-08-29 Masimo Corporation Intelligent medical network edge router
US8977348B2 (en) 2012-12-21 2015-03-10 Covidien Lp Systems and methods for determining cardiac output
AT514017B1 (en) 2013-02-22 2020-11-15 Dr Skrabal Falko Hemodynamic EKG
US10832818B2 (en) 2013-10-11 2020-11-10 Masimo Corporation Alarm notification system
AT516499B1 (en) 2015-04-22 2016-06-15 Skrabal Falko Dr Body impedance meter
JP6940483B2 (en) 2015-08-31 2021-09-29 マシモ・コーポレイション Wireless patient monitoring system and method
CN105207635A (en) * 2015-10-31 2015-12-30 惠州华阳通用电子有限公司 Method and device for automatic volume control
US10617302B2 (en) 2016-07-07 2020-04-14 Masimo Corporation Wearable pulse oximeter and respiration monitor
WO2018071715A1 (en) 2016-10-13 2018-04-19 Masimo Corporation Systems and methods for patient fall detection
US10524668B2 (en) 2018-02-05 2020-01-07 Aerobex, Inc. Method and apparatus for determination of left ventricular stroke volume and cardiac output using the arteries of the forearm by means of integration technique
WO2019204368A1 (en) 2018-04-19 2019-10-24 Masimo Corporation Mobile patient alarm display
JP2023518303A (en) 2020-03-20 2023-04-28 マシモ・コーポレイション Wearable device for non-invasive body temperature measurement
USD980091S1 (en) 2020-07-27 2023-03-07 Masimo Corporation Wearable temperature measurement device
USD974193S1 (en) 2020-07-27 2023-01-03 Masimo Corporation Wearable temperature measurement device
USD1000975S1 (en) 2021-09-22 2023-10-10 Masimo Corporation Wearable temperature measurement device
USD1048908S1 (en) 2022-10-04 2024-10-29 Masimo Corporation Wearable sensor

Citations (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US579349A (en) * 1897-03-23 Automatic movable dam ob sluiceway gate
US3340867A (en) * 1964-08-19 1967-09-12 Univ Minnesota Impedance plethysmograph
US4450527A (en) * 1982-06-29 1984-05-22 Bomed Medical Mfg. Ltd. Noninvasive continuous cardiac output monitor
US4548211A (en) * 1984-01-12 1985-10-22 Marks Lloyd A Computer assisted admittance plethysmograph
US4562843A (en) * 1980-09-29 1986-01-07 Ljubomir Djordjevich System for determining characteristics of blood flow
US4676253A (en) * 1985-07-18 1987-06-30 Doll Medical Research, Inc. Method and apparatus for noninvasive determination of cardiac output
US4807638A (en) * 1987-10-21 1989-02-28 Bomed Medical Manufacturing, Ltd. Noninvasive continuous mean arterial blood prssure monitor
US4836214A (en) * 1986-12-01 1989-06-06 Bomed Medical Manufacturing, Ltd. Esophageal electrode array for electrical bioimpedance measurement
US4953556A (en) * 1984-12-13 1990-09-04 Evans John M Method and apparatus for the measurement of thoracic field potentiometry
US5103828A (en) * 1988-07-14 1992-04-14 Bomed Medical Manufacturing, Ltd. System for therapeutic management of hemodynamic state of patient
US5178154A (en) * 1990-09-18 1993-01-12 Sorba Medical Systems, Inc. Impedance cardiograph and method of operation utilizing peak aligned ensemble averaging
US5309917A (en) * 1991-09-12 1994-05-10 Drexel University System and method of impedance cardiography and heartbeat determination
US5316004A (en) * 1990-12-28 1994-05-31 Regents Of The University Of Minnesota Method for vascular impedance measurement
US5423326A (en) * 1991-09-12 1995-06-13 Drexel University Apparatus and method for measuring cardiac output
US5469859A (en) * 1992-06-24 1995-11-28 N.I. Medical Ltd. Non-invasive method and device for collecting measurements representing body activity and determining cardiorespiratory parameters of the human body based upon the measurements collected
US5503157A (en) * 1995-03-17 1996-04-02 Sramek; Bohumir System for detection of electrical bioimpedance signals
US5505209A (en) * 1994-07-07 1996-04-09 Reining International, Ltd. Impedance cardiograph apparatus and method
US5685316A (en) * 1996-04-08 1997-11-11 Rheo-Graphic Pte Ltd. Non-invasive monitoring of hemodynamic parameters using impedance cardiography
US5782774A (en) * 1996-04-17 1998-07-21 Imagyn Medical Technologies California, Inc. Apparatus and method of bioelectrical impedance analysis of blood flow
US6016445A (en) * 1996-04-16 2000-01-18 Cardiotronics Method and apparatus for electrode and transthoracic impedance estimation
US6058325A (en) * 1996-04-16 2000-05-02 Cardiotronics Method and apparatus for high current electrode, transthoracic and transmyocardial impedance estimation
US6102869A (en) * 1995-09-12 2000-08-15 Heinemann & Gregori Gmbh Process and device for determining the cardiac output
US6186955B1 (en) * 1998-11-16 2001-02-13 Gail D. Baura Noninvasive continuous cardiac output monitor
US20020193689A1 (en) * 2001-04-03 2002-12-19 Osypka Medical Gmbh Apparatus and method for determining an approximation of the stroke volume and the cardiac output of the heart
US6561986B2 (en) * 2001-01-17 2003-05-13 Cardiodynamics International Corporation Method and apparatus for hemodynamic assessment including fiducial point detection
US6602201B1 (en) * 2000-07-10 2003-08-05 Cardiodynamics International Corporation Apparatus and method for determining cardiac output in a living subject
US20030153058A1 (en) * 1999-06-03 2003-08-14 Ajinomoto Co. Inc. Method for producing L-arginine
US6636754B1 (en) * 2000-07-10 2003-10-21 Cardiodynamics International Corporation Apparatus and method for determining cardiac output in a living subject
US6829501B2 (en) * 2001-12-20 2004-12-07 Ge Medical Systems Information Technologies, Inc. Patient monitor and method with non-invasive cardiac output monitoring

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6322518B1 (en) * 1993-12-06 2001-11-27 Heska Corporation Method and apparatus for measuring cardiac output
US5791349A (en) 1996-04-17 1998-08-11 Urohealth, Inc. Apparatus and method of bioelectrical impedance analysis of blood flow
US7822470B2 (en) 2001-10-11 2010-10-26 Osypka Medical Gmbh Method for determining the left-ventricular ejection time TLVE of a heart of a subject

Patent Citations (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US579349A (en) * 1897-03-23 Automatic movable dam ob sluiceway gate
US3340867A (en) * 1964-08-19 1967-09-12 Univ Minnesota Impedance plethysmograph
US4562843A (en) * 1980-09-29 1986-01-07 Ljubomir Djordjevich System for determining characteristics of blood flow
US4450527A (en) * 1982-06-29 1984-05-22 Bomed Medical Mfg. Ltd. Noninvasive continuous cardiac output monitor
US4548211A (en) * 1984-01-12 1985-10-22 Marks Lloyd A Computer assisted admittance plethysmograph
US4953556A (en) * 1984-12-13 1990-09-04 Evans John M Method and apparatus for the measurement of thoracic field potentiometry
US4676253A (en) * 1985-07-18 1987-06-30 Doll Medical Research, Inc. Method and apparatus for noninvasive determination of cardiac output
US4836214A (en) * 1986-12-01 1989-06-06 Bomed Medical Manufacturing, Ltd. Esophageal electrode array for electrical bioimpedance measurement
US4807638A (en) * 1987-10-21 1989-02-28 Bomed Medical Manufacturing, Ltd. Noninvasive continuous mean arterial blood prssure monitor
US5103828A (en) * 1988-07-14 1992-04-14 Bomed Medical Manufacturing, Ltd. System for therapeutic management of hemodynamic state of patient
US5178154A (en) * 1990-09-18 1993-01-12 Sorba Medical Systems, Inc. Impedance cardiograph and method of operation utilizing peak aligned ensemble averaging
US5316004A (en) * 1990-12-28 1994-05-31 Regents Of The University Of Minnesota Method for vascular impedance measurement
US5309917A (en) * 1991-09-12 1994-05-10 Drexel University System and method of impedance cardiography and heartbeat determination
US5423326A (en) * 1991-09-12 1995-06-13 Drexel University Apparatus and method for measuring cardiac output
US5469859A (en) * 1992-06-24 1995-11-28 N.I. Medical Ltd. Non-invasive method and device for collecting measurements representing body activity and determining cardiorespiratory parameters of the human body based upon the measurements collected
US5505209A (en) * 1994-07-07 1996-04-09 Reining International, Ltd. Impedance cardiograph apparatus and method
US5529072A (en) * 1995-03-17 1996-06-25 Sramek; Bohumir System for detection of electrical bioimpedance signals
US5503157A (en) * 1995-03-17 1996-04-02 Sramek; Bohumir System for detection of electrical bioimpedance signals
US6102869A (en) * 1995-09-12 2000-08-15 Heinemann & Gregori Gmbh Process and device for determining the cardiac output
US5685316A (en) * 1996-04-08 1997-11-11 Rheo-Graphic Pte Ltd. Non-invasive monitoring of hemodynamic parameters using impedance cardiography
US6161038A (en) * 1996-04-08 2000-12-12 Rheo-Graphic Pte Ltd. Non-invasive monitoring of hemodynamic parameters using impedance cardiography
US6016445A (en) * 1996-04-16 2000-01-18 Cardiotronics Method and apparatus for electrode and transthoracic impedance estimation
US6058325A (en) * 1996-04-16 2000-05-02 Cardiotronics Method and apparatus for high current electrode, transthoracic and transmyocardial impedance estimation
US6095987A (en) * 1996-04-17 2000-08-01 Imagyn Medical Techonologies California, Inc. Apparatus and methods of bioelectrical impedance analysis of blood flow
US5782774A (en) * 1996-04-17 1998-07-21 Imagyn Medical Technologies California, Inc. Apparatus and method of bioelectrical impedance analysis of blood flow
US6186955B1 (en) * 1998-11-16 2001-02-13 Gail D. Baura Noninvasive continuous cardiac output monitor
US20030153058A1 (en) * 1999-06-03 2003-08-14 Ajinomoto Co. Inc. Method for producing L-arginine
US6602201B1 (en) * 2000-07-10 2003-08-05 Cardiodynamics International Corporation Apparatus and method for determining cardiac output in a living subject
US6636754B1 (en) * 2000-07-10 2003-10-21 Cardiodynamics International Corporation Apparatus and method for determining cardiac output in a living subject
US6561986B2 (en) * 2001-01-17 2003-05-13 Cardiodynamics International Corporation Method and apparatus for hemodynamic assessment including fiducial point detection
US20020193689A1 (en) * 2001-04-03 2002-12-19 Osypka Medical Gmbh Apparatus and method for determining an approximation of the stroke volume and the cardiac output of the heart
US6511438B2 (en) * 2001-04-03 2003-01-28 Osypka Medical Gmbh Apparatus and method for determining an approximation of the stroke volume and the cardiac output of the heart
US6829501B2 (en) * 2001-12-20 2004-12-07 Ge Medical Systems Information Technologies, Inc. Patient monitor and method with non-invasive cardiac output monitoring

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100081960A1 (en) * 2008-09-30 2010-04-01 Nellcor Puritan Bennett Llc Bioimpedance System and Sensor and Technique for Using the Same
US8406865B2 (en) * 2008-09-30 2013-03-26 Covidien Lp Bioimpedance system and sensor and technique for using the same
US20140249440A1 (en) * 2010-12-28 2014-09-04 Matt Banet Body-worn system for continuous, noninvasive measurement of cardiac output, stroke volume, cardiac power, and blood pressure
US10856752B2 (en) * 2010-12-28 2020-12-08 Sotera Wireless, Inc. Body-worn system for continuous, noninvasive measurement of cardiac output, stroke volume, cardiac power, and blood pressure
EP2737847B1 (en) * 2011-07-29 2020-02-26 Universitat Politècnica De Catalunya Method and apparatus for obtaining cardiovascular information by measuring between two extremities
US20150324977A1 (en) * 2014-05-09 2015-11-12 Siemens Medical Solutions Usa, Inc. Image-Based Waveform Parameter Estimation
US9498140B2 (en) * 2014-05-09 2016-11-22 Siemens Medical Solutions Usa, Inc. Image-based waveform parameter estimation
US9451888B1 (en) * 2015-03-27 2016-09-27 Cordeus, Inc. Method and apparatus for determination of left ventricular stroke volume and cardiac output using the arteries of the forearm
US20190046069A1 (en) * 2015-07-10 2019-02-14 Bodyport Inc. Cardiovascular signal acquisition, fusion, and noise mitigation
US11696715B2 (en) 2015-07-10 2023-07-11 Bodyport Inc. Cardiovascular signal acquisition, fusion, and noise mitigation
WO2022182975A1 (en) * 2021-02-26 2022-09-01 Massachusetts Institute Of Technology Methods to simulate metrics of vascular function from clinical data
CN117281494A (en) * 2023-11-27 2023-12-26 安徽通灵仿生科技有限公司 Method and device for identifying signal characteristic points of arterial blood pressure signals

Also Published As

Publication number Publication date
WO2006063255A3 (en) 2007-04-05
US20090259132A1 (en) 2009-10-15
US7740590B2 (en) 2010-06-22
WO2006063255A2 (en) 2006-06-15
US7806830B2 (en) 2010-10-05

Similar Documents

Publication Publication Date Title
US7806830B2 (en) Apparatus and method for determination of stroke volume using the brachial artery
US7261697B2 (en) Apparatus for determination of stroke volume using the brachial artery
US6511438B2 (en) Apparatus and method for determining an approximation of the stroke volume and the cardiac output of the heart
Bernstein et al. Stroke volume equation for impedance cardiography
US9808168B2 (en) Method and system for non-invasive measurement of cardiac parameters
CN105832317B (en) System for measuring cardiac output, stroke volume, heart force and blood pressure
US6228033B1 (en) Apparatuses and methods for a noninvasive measurement of physiological parameters
US10617322B2 (en) System, method and apparatus for measuring blood flow and blood volume
US8414498B2 (en) System, method and apparatus for measuring blood flow and blood volume
Solà et al. Wearable PWV technologies to measure blood pressure: eliminating brachial cuffs
Shimazu et al. Noninvasive measurement of beat-to-beat vascular viscoelastic properties in human fingers and forearms
US9451888B1 (en) Method and apparatus for determination of left ventricular stroke volume and cardiac output using the arteries of the forearm
Nakonezny et al. New ambulatory impedance cardiograph validated against the Minnesota Impedance Cardiograph
US20170188963A1 (en) Physiological monitoring system featuring floormat and handheld sensor
Langer et al. Respiratory-induced hemodynamic changes measured by whole-body multichannel impedance plethysmography.
US20170188890A1 (en) Physiological monitoring system featuring floormat and handheld sensor
US20170188964A1 (en) Physiological monitoring system featuring floormat and handheld sensor
US10524668B2 (en) Method and apparatus for determination of left ventricular stroke volume and cardiac output using the arteries of the forearm by means of integration technique
US20170188944A1 (en) Physiological monitoring system featuring floormat and handheld sensor
Steiner Hidden information in three-axial ECG data of normal subjects: Fractal dimensions of corresponding points from successive QRS loops as a potential sport and age dependent marker
Solà et al. Wearable PWV technologies to measure Blood Pressure: getting rid of brachial cuffs
Karpov et al. Preeclampsia Amplitude-time Characteristics for ΔZ/Δt Curve
Buell Impedance cardiography and plethysmography
SENGTHIPPHANY HEART RATE VARIABILITY ANALYSIS FROM PHOTOPLETHYSMOGRAPHY DURING RESTING AND EXERCISE

Legal Events

Date Code Title Description
AS Assignment

Owner name: CORDEUS, INC., CALIFORNIA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:BERNSTEIN, DONALD P.;REEL/FRAME:016648/0162

Effective date: 20050904

FEPP Fee payment procedure

Free format text: PETITION RELATED TO MAINTENANCE FEES GRANTED (ORIGINAL EVENT CODE: PMFG); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY

Free format text: PETITION RELATED TO MAINTENANCE FEES FILED (ORIGINAL EVENT CODE: PMFP); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY

REMI Maintenance fee reminder mailed
LAPS Lapse for failure to pay maintenance fees
REIN Reinstatement after maintenance fee payment confirmed
PRDP Patent reinstated due to the acceptance of a late maintenance fee

Effective date: 20141125

FP Lapsed due to failure to pay maintenance fee

Effective date: 20141005

FPAY Fee payment

Year of fee payment: 4

STCF Information on status: patent grant

Free format text: PATENTED CASE

FEPP Fee payment procedure

Free format text: 7.5 YR SURCHARGE - LATE PMT W/IN 6 MO, SMALL ENTITY (ORIGINAL EVENT CODE: M2555)

MAFP Maintenance fee payment

Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2552)

Year of fee payment: 8

MAFP Maintenance fee payment

Free format text: PAYMENT OF MAINTENANCE FEE, 12TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2553); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY

Year of fee payment: 12