WO2009089183A1 - Non-invasive method and device for measuring cardiac output - Google Patents

Non-invasive method and device for measuring cardiac output Download PDF

Info

Publication number
WO2009089183A1
WO2009089183A1 PCT/US2009/030157 US2009030157W WO2009089183A1 WO 2009089183 A1 WO2009089183 A1 WO 2009089183A1 US 2009030157 W US2009030157 W US 2009030157W WO 2009089183 A1 WO2009089183 A1 WO 2009089183A1
Authority
WO
WIPO (PCT)
Prior art keywords
heart
signal
person
magnitude
frequency
Prior art date
Application number
PCT/US2009/030157
Other languages
French (fr)
Inventor
Stephen C. Jacobsen
Tomasz J. Petelenz
Stephen C. Peterson
Original Assignee
Raytheon Sarcos, Llc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Raytheon Sarcos, Llc filed Critical Raytheon Sarcos, Llc
Priority to EP09700309A priority Critical patent/EP2237719A1/en
Priority to CN200980107669.8A priority patent/CN101983030B/en
Priority to JP2010541585A priority patent/JP5406211B2/en
Priority to MX2010007414A priority patent/MX2010007414A/en
Publication of WO2009089183A1 publication Critical patent/WO2009089183A1/en

Links

Classifications

    • 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 
    • 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/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/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/0507Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  using microwaves or terahertz waves

Definitions

  • the present invention relates to medical devices, and more particularly to devices and methods of operation of said devices for non-invasively measuring the cardiac output of a patient.
  • Hemorrhage is a major cause of death of battlefield casualties who are not killed immediately and survive beyond the initial five minutes from injury. Half of battlefield deaths are the result of uncontrolled bleeding.
  • Successful treatment requires immediate medical attention that can be very limited in the battlefield until the casualty is evacuated to higher echelons of care.
  • Monitoring of the physiological status during transport and the end points of the resuscitative treatment is critical to both immediate survival as well as long-term treatment of a patient. It has been suggested that severely injured casualties should be evacuated in less than one hour to a front-line high echelon care unit and that during initial treatment and evacuation, physiological monitoring should be used to guide resuscitation efforts.
  • Cardiac output is an important indicator of hemodynamic status of the circulatory system and is important in diagnosis and treatment of heart disease and in guiding fluid resuscitation following major trauma with hemorrhage.
  • Cardiac output provides a direct measure of the heart's ability to pump blood and is determined as a product of stroke volume and heart rate.
  • As the hemorrhage progresses arterial blood pressure decreases and sympathetic compensatory mechanisms are activated that shunt the blood from the peripheral to central compartment in order to maintain blood pressure and thus oxygenation of vital organs. In healthy and physically trained individuals, these compensatory mechanisms maintain blood pressure disproportionately high within the body only to cause a rapid circulatory collapse as the hemorrhage causes massive loss of blood volume.
  • Cardiac output reacts faster to the progression of hemorrhagic shock and drops sooner than the arterial pressure in response to blood loss during hemorrhage.
  • Early intervention and resuscitation based on monitoring of cardiac output is the most complete way of capturing the physiological impact of the hemorrhage and resulting circulatory shock.
  • arterial blood pressure is not available under battlefield conditions and is commonly substituted by traditional, non-invasive (occlusive) blood pressure measurements in peripheral circulation.
  • the compensatory mechanisms that maintain blood supply to the vital organs do so by increasing peripheral resistance and diverting blood from the peripheral to central pool, causing the cessation of pressure pulse in peripheral circulation, and further difficulties in the determination of blood pressure.
  • the present invention seeks to overcome these by providing a method of measuring cardiac output of a person, the method comprising placing a device about the mid-sternal position on the person, said device comprising a RF signal generator coupled to an antenna, said antenna having a dielectric material disposed about the exterior of the antenna.
  • the method further comprises propagating a first signal having a predetermined frequency towards the heart of the person and receiving and measuring a portion of the first signal returned from the heart of the person with the device.
  • the method further comprises comparing the magnitude of the first signal propagated into the heart of the person to the magnitude of the portion of the first signal returned from the heart of the person and calculating a return loss of the signal.
  • the method comprises propagating an additional signal into a portion of the heart of the person using the device, wherein said additional signal has a frequency different from the first signal and receiving and measuring a portion of the additional signal returned from the heart of the person with the device.
  • the method further comprises comparing the magnitude of the additional signal propagated into the heart to the person to the magnitude of the additional signal returned from the heart of the person and calculating a return loss of the additional signal as well as comparing the return loss of the first signal to the return loss of the additional signal.
  • the method further comprises repeating certain steps referenced above while incrementally varying the frequency of the signal with each repetitive step and determining a maximum return loss value of the signals propagated into the heart.
  • the present invention improves upon deficiencies in the prior art by providing a system for measuring the cardiac output of the heart of a patient, comprising a hand-held device comprising a housing containing a signal generator coupled to an antenna and a dielectric material disposed about the antenna, wherein said hand-held device is adapted to generate a plurality of radio frequency signals and direct said signals towards the heart of the person and measure a magnitude of a signal returned from the heart of the person.
  • the hand-held device further comprises a processor contained within said housing, wherein said processor is adapted to compare differences between a magnitude of a signal propagated into the heart of the person and the magnitude of the signal returned off the heart.
  • the processor is further adapted to determine a signal frequency having a maximum return loss value based on the differences between the magnitude of the signal propagated into the heart of the person and the magnitude of the signal returned off of the heart of the person. Additionally, the processor is further adapted to estimate a change in the amplitude of motion of a portion of a wall of the heart based on the differences between the magnitude of the signal propagated by the device and the magnitude of the signal returned off of the portion of the heart.
  • FIG. 1 shows one embodiment of an RF device in accordance with one embodiment of the present invention
  • FIG. 2 shows a schematic block diagram of a RF device in accordance with one embodiment of the present invention
  • FIG. 5 shows a comparison of the sensor signal with and without motion artifact reduction in accordance with one embodiment of the present invention
  • FIG. 6 shows a record of pressure sensor stabilizer data in accordance with one embodiment of the present invention
  • FIG. 7 shows a comparison of an ECG trace synchronously recorded with signals using an RF system in accordance with one embodiment of the present invention
  • FIG. 8 illustrates an estimated depth of penetration of an EM wave in various tissues in accordance with one embodiment of the present invention
  • FIG. 9 is an illustration of a frequency tuning sweep and resulting work point in accordance with one embodiment of the present invention.
  • FIG. 10 shows a data plot of the sensor magnitude versus pulmonary artery stroke volume from experiment SW3 in accordance with one embodiment of the present invention
  • FIG. 1 1 shows the bleed summary of experiment SW3 in accordance with one embodiment of the present invention
  • FIG. 12 shows a comparison of stroke volume derived from PA flow test and sensor signal during full bleed and re-infusion during experiment SW3 in accordance with one embodiment of the present invention.
  • FIG. 13 shows calibration lines for PEEP experiments 4, 6, 13, and 14 in accordance with one embodiment of the present invention.
  • a method and device for measuring the cardiac output of a patient.
  • the method and device are based on the general operating principle that alternating current electrical signals used to produce and detect radio waves may be utilized to detect the status of biological systems within the body.
  • radio waves are propagated into the body. A certain portion of those waves are absorbed by the body. A certain portion of the waves not absorbed by the body are returned back near their point of origin. The difference between the magnitude of the waves sent into the body and the magnitude of the waves returned back near their point of origin (referred to herein as "return loss”) is utilized to estimate certain characteristics of biological tissue (e.g., cardiac tissues) within the body.
  • a device 10 for estimating changes in the cardiac output of a patient comprising a radio frequency (“RF") signal generator coupled to an antenna 15, said antenna 15 having a dielectric material 20 disposed about an exterior surface of the antenna 15.
  • the dielectric material 20 may be disposed about select portions of the exterior surface of the antenna 15 or substantially about the entire exterior surface of the antenna 15.
  • the RF signal has a near field component and a far field component.
  • the dielectric material 20 has a predetermined thickness and a predetermined dielectric constant such that a substantial portion of the near field component of the electronic signal returned from the patient is not received by the antenna 15.
  • the dielectric material 20 is designed such that substantially the entire near field component of the electronic signal returned from the patient is not received by the antenna 15.
  • the dielectric material 20 has a dielectric constant ranging from 5 to 25 and has a thickness ranging from 2 mm to 20 mm.
  • the design of the dielectric material 20 is dependant upon numerous characteristics which vary depending on the type of dielectric material used. That is, the thickness of the material is a function of its dielectric constant and thus it is specific for a given material. Accordingly, no limitation should be read into the exact specifications of the dielectric material used herein.
  • the RF signal is generated by a voltage controlled oscillator 25 ("VCO") contained within a device housing 11.
  • VCO voltage controlled oscillator
  • the VCO 25 generates the RF signal at differing frequencies by varying the tuning voltage.
  • VCO's 25 which may be used include the ROS-2150 made by Minicircuits, however any VCO 25 capable of generating an appropriate RF signal may be utilized.
  • a micro-controller 30 via a digital to analog (“D/A") converter 31 sweeps the frequency of the RF signal between approximately 1000 and 2000 MHz by causing the tuning voltage to vary between approximately zero and 20 volts.
  • the output of the VCO 25 is less than 4 mW. While any suitable operating frequency may be used, in one embodiment of the present invention, the VCO 25 operates at a frequency ranging from 0.5 GHz to 2.5 GHz.
  • the device 10 further comprises a network analyzer 12 such as that manufactured by Analog Devices (AD8302 IC).
  • the detector is capable of identifying the magnitude and phase of subject signals and has a 60 dB dynamic range and outputs a magnitude voltage of 30 mV per dB of return loss.
  • the phase output is 10 mV per degree over 180 degrees.
  • the device 10 further comprises microcontroller 30 (such as those manufactured by Cygnal) to read the network analyzer outputs, determine the minimum or maximum values of magnitude and phase during a frequency sweep and generate a digital sweep that enters the D/A converter 31 that then drives the VCO 25. This closed loop generates the VCO sweep, determines the frequency of maximum return loss (magnitude or phase), and then locks on to that frequency.
  • the microcontroller 30 tracks the magnitude and phase of the "organ's" signal.
  • a low pass filter 35 (“LPF") is used to prevent the
  • the transmitted RF energy goes to the antenna 15 where the coupled signals are routed to the network analyzer 12.
  • a voltage multiplier 32 and appropriately placed attenuators 33 may also be used in accordance with certain aspects of the invention.
  • appropriate programming control and display interfaces are incorporated on the face of the RF device.
  • the programming control and display interfaces are located on a remote programming device.
  • FIG. 3 a block diagram of the RF device with analog detection means in accordance with one embodiment of the present invention is illustrated. While similar to the digital detection system shown in FIG. 2, additional components are shown to achieve proper detection of an analog signal. Additional components include a sawtooth generator 36, a generator trigger signal 37 and an amplifier/signal conditioner/peak detector/analog digital converter device 41.
  • the RF components discussed herein are fitted to a 4 layer PCB with 50-ohm traces.
  • the RF board components are shielded and, in some embodiments, feed through capacitors are used to prevent RF signals from "leaking out” of the shield via the DC or low frequency signals.
  • the device is capable of measuring the signal transmitted from the VCO 25 to an antenna 15 where the signal has been modulated by the organ system of the patient.
  • the device 10 comprises an additional dipole antenna 15a, an additional LPF 35a, a single directional coupler and a second network analyzer 12a, one of whose inputs is from the VCO 25. Output from the network analyzer 12a is routed to the data collection system.
  • the device further comprises a switch sensor for activating and deactivating the device 10.
  • the switch sensor would activate the device 10 only when it is in full contact with the body of the patient and would deactivate the device 10 when not in full contact with the body.
  • errant RF radiation emitted into the surrounding environment would be minimized by using the switch sensor.
  • the switch sensor may comprise a touch-activated switch sensor or other desirable switch sensor as suits a particular application.
  • Coupling of the RF signal into the body is maximized when there is a good impedance match between the coupler and the body. It is believed that biological tissues present about a 50-ohm load for the RF signal. Appropriately placed dipoles and loop antennas provide as much as 50 dB return loss at frequencies when placed near the body, even through clothing. As the antenna 15 nears the body, the tissue loads the antenna 15 changing its resonance to below that of free-space resonance. In accordance with one aspect of the invention, the measurement of the return loss is performed at 1 millisecond intervals, resulting in I kilosample per second signal acquisition rate.
  • the core of the system is the network analyzer 12 which measures the ratio of the returned RF signal to the sampled outgoing signal generated by the VCO 25.
  • the optimal frequency of the VCO 25 is found by scanning the frequency range within approximately 1-2 GHz range and determining the frequency of the maximum return loss. The magnitude, phase, and frequency at the maximum return loss point is then digitized and displayed and/or stored in the data file. Additional functions of the device included acquisition of up to 4 analog channels (e.g., flow, impedance, ECG, other) of 0-5V range and synchronizing the analog data channels with the RF magnitude, phase and frequency.
  • a dipole antenna is utilized which comprises copper foil tape applied to a suitable backing substrate.
  • a 50-ohm coaxial cable connects the dipole antenna to a processor.
  • a circular loop antenna was employed for which the components of the electric and magnetic field as a function of distance from the antenna are given by numerous equations. These equations represent fields resulting from a sinusoidally-excited current loop with a magnetic moment "dm". Electric and magnetic components of the field are given by the relationships:
  • the results of tests of monitoring of the heart and breathing action in human volunteers in the presence and absence of the interfacial motion are shown in FIG. 5.
  • the output signal of the antenna 15 having the dielectric material 20 was attenuated as compared to the signal measured with the antenna 15 without the dielectric material 20.
  • the signal was much less sensitive to motion artifacts induced by moving the antenna 15 with respect to the surface of the patient.
  • the motion artifact-related signal originated in the near-field of the antenna 15 (i.e., at the interface between the antenna and the body).
  • the near-field of the antenna 15 i.e., at the interface between the antenna and the body.
  • complex electric field modulation is dominated by fast changing 1/r 2 function, whereas in the far field, these interactions result in a slower 1/r relationship.
  • Placing a dielectric material 20 on the surface of the antenna 15 causes shortening of the wavelength in the dielectric region by a factor of l/(sqrt( ⁇ r)).
  • a proper selection of the dielectric material 20 and its geometry (thickness) permits enclosing the near-field region entirely in the dielectric material, resulting in the significant reduction of the device sensitivity to interfacial motion as noted above.
  • any desirable dielectric material may be utilized as suits a particular application.
  • Example dielectric materials include ceramics, glass, plastics, and polymers. Fluid dielectric materials may also be used as desired including, but not limited to gels, liquids, polymeric fluids, and hydrocarbon-based fluids. Additionally, while specific reference has been made to a specific thickness, it is important to note that the thickness of the dielectric material 20 will vary according to the type of dielectric material used.
  • a pressure sensor is coupled to the device 10 to detect the application force and motion of the antenna 15 with reference to the body.
  • Sensor software can be created to permit valid data measurements only while the sensor is balanced and indicates "no-motion" conditions.
  • pressure sensors are disposed on all four corners of the antenna 15 and are connected to a quad bridge circuit.
  • three strain gages are built into the applicator shield and used in the same manner.
  • Software is programmed into the network analyzer 12 to visually show the strain on each gage and prevent measurement if the difference in strain exceeds an adjustable threshold.
  • the method comprises placing a device 10, such as that described in the preceding sections of this application, about the mid-sternal position of the person.
  • the device 10 comprises an RF signal generator 25 coupled to an antenna 15, said antenna 15 having a dielectric material 20 disposed about the exterior 16 of the antenna 15.
  • the method further comprises propagating a first signal having a predetermined frequency towards the heart of the person and receiving and measuring a portion of the first signal returned from the heart of the person with the device 10.
  • the method further comprises comparing the return loss of the first signal to the return loss of the additional signal and repeating certain of the above-referenced steps while incrementally varying the frequency of the signal with each repetitive step.
  • the method also comprises determining a maximum return loss value of the signals propagated into the heart, of the person over the utilized frequency range. As noted above, in one embodiment, the frequency of the signal is swept between 1000 and 2000 MHz in an effort to determined the maximum return loss for a particular patient at a particular frequency.
  • optimal results are obtained when the device 10 is positioned in the mid-sternal position. That is, the optimal results are obtained when the signals from the device are predominantly focused on the right ventricle of the patient's heart.
  • the cardiac cycle consists of two major periods, systole, during which the ventricles contract and, under high pressure, eject blood into the pulmonary artery and into the aorta, and diastole, during which the heart is in its low pressure state and the ventricles are filling with blood.
  • systole the ventricles contract and, under high pressure, eject blood into the pulmonary artery and into the aorta
  • diastole during which the heart is in its low pressure state and the ventricles are filling with blood.
  • the heart is maximally extended, while at the beginning of the diastole the extension of the heart muscle is minimal.
  • the systole starts with the R-wave, which is immediately followed by contraction of the left ventricle, and ends after the T-wave with the closing of the aortic valve.
  • tension developed in cardiac muscle fibers is proportional to their stretching and thinning, or to the volume of blood filling the heart, or End Diastolic Volume ("EDV").
  • EDV End Diastolic Volume
  • the EDV is one of the factors that determine the size of the subsequent stroke volume. If a signal is proportional to the volume of the right ventricle, the magnitude of the signal should vary between cardiac cycles with the venous return, or preload. Therefore, the device acts as a sensitive indicator of the physiological events that affect venous return.
  • both tricuspid and pulmonic valves are closed, and the initial ventricular contraction results in an isovolumetric intraventricular pressure build-up.
  • FIG. 7 based on the comparison of the left ventricular pressure wave and S-wave timing, components of a recorded RF signal can be identified that correspond to mechanical contraction of the heart and thus provide information about changes in the cardiac output.
  • an analysis begins with the R-wave, which initiates contraction of the ventricles. Upon contraction, with some delay related to spreading of the electrical wave, the ventricles contract isovolumetrically (both valves are still closed).
  • the RS trade rises initially slowly, then rapidly, indicating a thickness change of the ventricular wall as detected by the antenna, consistent with contraction of the left ventricle.
  • EDP End Diastolic Pressure
  • the aortic valve opens and blood is ejected from the left ventricle into the aorta, initially rapidly, and then slowly as the aortic pressure rises to the End Systolic Level, resulting in reduced low.
  • the left ventricle reaches its minimum volume (maximum point on the RS curve), and starts expanding.
  • the RS reaches its maximum. That is, the heart wall has reached its maximum thickness.
  • the aortic valve closes.
  • the left ventricle continues to relax, with both valves closed (isovolumic relaxation), resulting in reduced thickness of the ventricular wall, which appears on the RS as a section with significantly reduced slope.
  • the mitral valve opens, initiating the period of rapid filling.
  • the ventricular wall continues thinning as blood fills the ventricle.
  • Both mitral and aortic valves are closed, and the left ventricle continues its (isovolumic) relaxation.
  • the mitral valve opens, and the ventricle fills with blood, expanding its size and causing a thinning of the wall.
  • the minimum of the RS curve corresponds to a filled ventricle ready for ejection following atrial contraction (P-wave). Following the P-wave, the mitral valve closes, and the cycle begins again. Similar pressure - valve opening/closing relationships are also true for the right atrium, right ventricle closing and opening of the pulmonic valve during systole.
  • the device described in detail above utilizes the interaction of electromagnetic energy with biological tissues.
  • the RF device utilizes the interaction of electromagnetic energy with biological tissues.
  • the FCC has compiled a database of these properties at http://www.fcc.gov/fcc-bin/dielec.sh and this database was used extensively throughout the research and development of the RF device described herein.
  • EM electromagnetic
  • ⁇ + j ⁇ where gamma is the complex propagation constant, alpha is the attenuation constant and beta is the phase constant is helpful for understanding modeling of the device signals. It is important to note that these "constants" are constant for only one frequency and particular material. That is, as the signal frequency is varied and the subject tissue is changed, the constants are variable. Beta is equal to 2*pi/lambda, the wavelength in the material. Alpha has the units of Nepers/meter.
  • omega the radian frequency
  • the permeability of free space
  • epsilon the permittivity of the material
  • sigma is the material's conductivity (Siemens/meter).
  • the intrinsic impedance of a lossy dielectric is given by:
  • EM energy is attenuated by lossy dielectrics and returned at impedance boundaries.
  • a useful equation is the transmission line equation that describes the impedance seen by an EM wave as it traverses materials of differing impedance.
  • ⁇ in ⁇ (( ⁇ 2 +j ⁇ i*tanpi)/( ⁇ ,+j ⁇ 2*tan ⁇ l)) where 1 is the thickness of material 1.
  • the intrinsic impedances of muscle and blood are estimated to be: muscle (49.2+j 5.91) and blood (46.5+J7.69) (note that the real part of the impedances is very close to 50 ohms so the tissues present a good match to 50 ohm cable).
  • P (Zi-Zo)/ (Zj +Zo) where Zi and Zo are the complex impedances of the load (blood) and line (muscle).
  • p (- 0.0247, ⁇ .022) and about 2.5% of the energy is reflected at the muscle blood interface.
  • One objective of developing a simplified propagation model for calculating RF properties of the antenna as a function of cardiac output and presence/absence of hemo/pneumothorax was to determine whether, given the dielectric properties of the tissues at the frequency range between 1 and 2GHz, using a simplified 2-D slab configuration.
  • the model is based on the interaction of a planar electromagnetic wave with the body tissues modeled as slabs and has not been designed to provide predictive information, but to assist in understanding the behavior of the reflection coefficient as a function of the changes in thickness of muscle, blood, and air and thus help in selecting the frequency of the measurement and characteristic waveform of the recorded sensor signal as a function of the beating heart. It is believed that the model assists in assessing the mechanisms responsible for modulation of the return loss signal as a function of either ventricular volume and of the movement of the ventricular wall. Similarly, introduction of an air pocket, as in pneumothorax, would create an impedance change that may be predicted by the model and therefore measurable by a suitable RF device. Selection of the working frequency impacts the depth of penetration of the RF radiation into the tissues.
  • each layer corresponding to different tissue forms a load for the EM wave traveling through the immediately preceding layer of tissue, and calculations are completed iteratively for all the layers of tissues, over a range of frequencies and dimensions.
  • Zi is the complex impedance of the load
  • Zo is the complex input impedance "looking into” the line.
  • Zi is the characteristic impedance of the chest wall, which is assumed to be a 2cm thick layer with dielectric properties corresponding to those of the muscle, bone and fat in 1 : 1 : 1 proportion.
  • the input impedance at each interface can be computed using the formula:
  • is the phase constant parameter representing lossy (imaginary) part of the complex propagation constant
  • ⁇ + j ⁇
  • is the attenuation constant
  • Z 1 SQRT ( ⁇ / ⁇ r ⁇ o)*[l/SQRT(l-j ⁇ / ⁇ ⁇ r ⁇ 0 ))]
  • FIG. 9 is an illustration of a frequency tuning sweep and resulting work point in accordance with one embodiment of the present invention.
  • animal testing was conducted to establish metrological correlation between the measured sensor signal and the stroke volume/cardiac output tests that are impossible, or at best extremely dangerous to perform on human subjects.
  • the animal studies included (i) measurements of cardiac output in the pulmonary artery during hemorrhage, and (ii) measurements of cardiac output using the device described herein as a function of blood loss during controlled hemorrhage and in correlation with the Pulmonary Artery Flow ("PAF").
  • PAF Pulmonary Artery Flow
  • Sensor readings were recorded during: (i) a control period prior to hemorrhage, (ii) bleeding and re-infusion of blood at a rate of 100 ml/min, up to approximately thirty-three percent of the estimated blood volume of the animal, and (iii) terminal bleeding.
  • the following experimental variables were continuously recorded in a personal computer: (i) ECG using 3-lead standard ECG configuration, (H) breathing, (iii) sensor output, AC- coupled, (iv) sensor output, DC-coupled, (v) PA flow, instantaneous, (vi) PA flow, average, and (vii) frequency (i.e., the optimized maximum return loss frequency).
  • the sensor output signal contains two components: an additive, breathing and motion-modulated baseline at breathing and motion-artifact frequencies, and the stroke volume-related signal at heart beat frequency. The separation of these components is performed in the frequency domain, and is followed by subsequent interpretation of the filtered signal in time domain.
  • the digital filtering discussed herein includes the following steps: (i) calculating the FFT of the selected data segments, (ii) plotting the power spectrum, (iii) dividing the power spectrum of the signal into two separate spectra containing the breathing-related and (separately) the heart-beat-related frequency components of the original signal, (iv) re-creating the heart-related and breathing-related waveforms using inverse FFT function, and (v) plotting the peak amplitudes of the sensor signal versus the corresponding peak amplitudes of the PA flow or volume bled/re- infused.
  • the frequency bandwidth used for the analysis is
  • FIG. 10 shows one implementation of the above-described filter using MatLab data processing software for a selected data segment. Specifically, FIG. 10 shows a data plot of the sensor magnitude versus pulmonary artery stroke volume from experiment SW3 in accordance with one embodiment of the present invention. FIG. 11 provides a summary data from the bleed experiment in SW3 experiment.
  • FIG. 12 shows a bleed re-infusion cycle in the animal of SW3. Stroke volume derived by integrating PA instantaneous flow is plotted versus blood loss/gain. From the same data record, the bottom panel shows the sensor signal plotted similarly.
  • FIG. 13 provides correlation graph for pooled unsealed PEEP08 and PEEPl 4 data. The correlation coefficient, R — 0.91 16, indicated a strong correlation between stroke volume derived from COTB measurements, and that from the RF sensor signal.

Abstract

A system comprising a housing containing a signal generator coupled to an antenna and a dielectric material disposed about the antenna. The device is adapted to generate and direct a plurality of signals towards the heart of the person and measure a magnitude of a signal returned from the heart. The device further comprises a processor to compare differences between a magnitude of a signal propagated and the magnitude of the signal returned off the heart and determine a signal frequency having a maximum return loss value based on those differences. The processor also estimates a change in the amplitude of motion of a portion of a wall of the heart based on the differences between the magnitude of the signal propagated by the device and the magnitude of the signal returned off of the portion of the heart.

Description

NON-INVASIVE METHOD AND DEVICE FOR MEASURING CARDIAC
OUTPUT
PRIORITY CLAIM
This application claims the benefit of United States Provisional Application No. 61/010,035 filed on January 4, 2008 entitled "Non-Invasive Method and Device for Measuring Cardiac Output" the entirety of which is incorporated herein by reference.
GOVERNMENT
This invention was made with government support under Grant No. DAMD1700C0013 awarded by the Department of Defense. The government has certain rights in the invention.
FIELD OF THE INVENTION
The present invention relates to medical devices, and more particularly to devices and methods of operation of said devices for non-invasively measuring the cardiac output of a patient.
BACKGROUND
There are numerous parameters that are measured invasively and non-invasively to gain diagnostic information about the cardiovascular system. The major fiinction of the cardiovascular system is to supply sufficient amounts of oxygen and nutrients to meet the metabolic demands of the tissues. Cardiac output is a primary determinant of global oxygen transport from the heart to the body and cardiovascular insufficiencies are reflected in cardiac output. Measurement of cardiac output provides both diagnostic and prognostic information, and the means to monitor the adequacy of therapy, as well as to assess the condition of a patient following a major trauma with hemorrhage and during the ensuing fluid resuscitation.
Hemorrhage is a major cause of death of battlefield casualties who are not killed immediately and survive beyond the initial five minutes from injury. Half of battlefield deaths are the result of uncontrolled bleeding. Successful treatment requires immediate medical attention that can be very limited in the battlefield until the casualty is evacuated to higher echelons of care. Monitoring of the physiological status during transport and the end points of the resuscitative treatment is critical to both immediate survival as well as long-term treatment of a patient. It has been suggested that severely injured casualties should be evacuated in less than one hour to a front-line high echelon care unit and that during initial treatment and evacuation, physiological monitoring should be used to guide resuscitation efforts.
Cardiac output is an important indicator of hemodynamic status of the circulatory system and is important in diagnosis and treatment of heart disease and in guiding fluid resuscitation following major trauma with hemorrhage. Cardiac output provides a direct measure of the heart's ability to pump blood and is determined as a product of stroke volume and heart rate. As the hemorrhage progresses, arterial blood pressure decreases and sympathetic compensatory mechanisms are activated that shunt the blood from the peripheral to central compartment in order to maintain blood pressure and thus oxygenation of vital organs. In healthy and physically trained individuals, these compensatory mechanisms maintain blood pressure disproportionately high within the body only to cause a rapid circulatory collapse as the hemorrhage causes massive loss of blood volume. Cardiac output, in turn, reacts faster to the progression of hemorrhagic shock and drops sooner than the arterial pressure in response to blood loss during hemorrhage. Early intervention and resuscitation based on monitoring of cardiac output is the most complete way of capturing the physiological impact of the hemorrhage and resulting circulatory shock. In addition to the delay in blood pressure drop caused by the compensatory mechanisms, arterial blood pressure is not available under battlefield conditions and is commonly substituted by traditional, non-invasive (occlusive) blood pressure measurements in peripheral circulation. As mentioned above, the compensatory mechanisms that maintain blood supply to the vital organs do so by increasing peripheral resistance and diverting blood from the peripheral to central pool, causing the cessation of pressure pulse in peripheral circulation, and further difficulties in the determination of blood pressure. Although blood pressure is currently used to monitor the progression of hemorrhagic shock and the effectiveness of fluid resuscitation, direct measurement of cardiac output would provide a more effective, direct measure of both and, consequently, a more effective therapy. Clinicians have an increasing number of available sensors and instrumentation that support measurements of cardiac output. These include, without limitation, indicator dilution techniques with or without the use of a pulmonary artery catheter, arterial pulse contour techniques, aortic pulsed Doppler, both of the ascending and descending aortas, indirect measures using arterio- venous gas content differences and expired gas measures via the Fick's equation, and bio-impedance techniques. Due to the high costs and need for highly specialized medical personnel and well-equipped facilities, as well as potentially severe complications associated with invasive measurement of cardiac output, noninvasive techniques are highly desirable. For battlefield applications, only non- invasive methods that require minimal skills and that can be deployed under "far forward" and "casualty transport" conditions are practical. None of the above methods has as yet fulfilled these requirements.
SUMMARY OF THE INVENTION In light of the problems and deficiencies inherent in the prior art, the present invention seeks to overcome these by providing a method of measuring cardiac output of a person, the method comprising placing a device about the mid-sternal position on the person, said device comprising a RF signal generator coupled to an antenna, said antenna having a dielectric material disposed about the exterior of the antenna. The method further comprises propagating a first signal having a predetermined frequency towards the heart of the person and receiving and measuring a portion of the first signal returned from the heart of the person with the device. The method further comprises comparing the magnitude of the first signal propagated into the heart of the person to the magnitude of the portion of the first signal returned from the heart of the person and calculating a return loss of the signal. Additionally, the method comprises propagating an additional signal into a portion of the heart of the person using the device, wherein said additional signal has a frequency different from the first signal and receiving and measuring a portion of the additional signal returned from the heart of the person with the device. The method further comprises comparing the magnitude of the additional signal propagated into the heart to the person to the magnitude of the additional signal returned from the heart of the person and calculating a return loss of the additional signal as well as comparing the return loss of the first signal to the return loss of the additional signal. The method further comprises repeating certain steps referenced above while incrementally varying the frequency of the signal with each repetitive step and determining a maximum return loss value of the signals propagated into the heart.
In accordance with additional embodiments, the present invention improves upon deficiencies in the prior art by providing a system for measuring the cardiac output of the heart of a patient, comprising a hand-held device comprising a housing containing a signal generator coupled to an antenna and a dielectric material disposed about the antenna, wherein said hand-held device is adapted to generate a plurality of radio frequency signals and direct said signals towards the heart of the person and measure a magnitude of a signal returned from the heart of the person. The hand-held device further comprises a processor contained within said housing, wherein said processor is adapted to compare differences between a magnitude of a signal propagated into the heart of the person and the magnitude of the signal returned off the heart. The processor is further adapted to determine a signal frequency having a maximum return loss value based on the differences between the magnitude of the signal propagated into the heart of the person and the magnitude of the signal returned off of the heart of the person. Additionally, the processor is further adapted to estimate a change in the amplitude of motion of a portion of a wall of the heart based on the differences between the magnitude of the signal propagated by the device and the magnitude of the signal returned off of the portion of the heart.
BRIEF DESCRIPTION OF THE DRAWINGS
The present invention will become more fully apparent from the following description and appended claims, taken in conjunction with the accompanying drawings. Understanding that these drawings merely depict exemplary embodiments of the present invention they are, therefore, not to be considered limiting of its scope. It will be readily appreciated that the components of the present invention, as generally described and illustrated in the figures herein, could be arranged and designed in a wide variety of different configurations. Nonetheless, the invention will be described and explained with additional specificity and detail through the use of the accompanying drawings in which:
FIG. 1 shows one embodiment of an RF device in accordance with one embodiment of the present invention; FIG. 2 shows a schematic block diagram of a RF device in accordance with one embodiment of the present invention;
FIG. 3 shows a schematic block diagram of a RF system in accordance with one embodiment of the present invention; FIG. 4 illustrates calculated contributions to a complex electric field of certain terms in a near field region of a circular loop antenna in a material with a relative dielectric constant in accordance with one embodiment of the present invention;
FIG. 5 shows a comparison of the sensor signal with and without motion artifact reduction in accordance with one embodiment of the present invention; FIG. 6 shows a record of pressure sensor stabilizer data in accordance with one embodiment of the present invention;
FIG. 7 shows a comparison of an ECG trace synchronously recorded with signals using an RF system in accordance with one embodiment of the present invention;
FIG. 8 illustrates an estimated depth of penetration of an EM wave in various tissues in accordance with one embodiment of the present invention;
FIG. 9 is an illustration of a frequency tuning sweep and resulting work point in accordance with one embodiment of the present invention;
FIG. 10 shows a data plot of the sensor magnitude versus pulmonary artery stroke volume from experiment SW3 in accordance with one embodiment of the present invention;
FIG. 1 1 shows the bleed summary of experiment SW3 in accordance with one embodiment of the present invention;
FIG. 12 shows a comparison of stroke volume derived from PA flow test and sensor signal during full bleed and re-infusion during experiment SW3 in accordance with one embodiment of the present invention; and
FIG. 13 shows calibration lines for PEEP experiments 4, 6, 13, and 14 in accordance with one embodiment of the present invention.
DESCRIPTION OF EXAMPLE EMBODIMENT^ The following detailed description of exemplary embodiments of the invention makes reference to the accompanying drawings, which form a part hereof and in which are shown, by way of illustration, exemplary embodiments in which the invention may be practiced. While these exemplary embodiments are described in sufficient detail to enable those skilled in the art to practice the invention, it should be understood that other embodiments may be realized and that various changes to the invention may be made without departing from the spirit and scope of the present invention. Thus, the following more detailed description of the embodiments of the present invention is not intended to limit the scope of the invention, as claimed, but is presented for purposes of illustration only to describe the features and characteristics of the present invention, and to sufficiently enable one skilled in the art to practice the invention. Accordingly, the scope of the present invention is to be defined solely by the appended claims, The following detailed description and exemplary embodiments of the invention will be best understood by reference to the accompanying drawings, wherein the elements and features of the invention are designated by numerals throughout.
In accordance with one embodiment of the invention, a method and device are disclosed for measuring the cardiac output of a patient. The method and device are based on the general operating principle that alternating current electrical signals used to produce and detect radio waves may be utilized to detect the status of biological systems within the body. Generally speaking, radio waves are propagated into the body. A certain portion of those waves are absorbed by the body. A certain portion of the waves not absorbed by the body are returned back near their point of origin. The difference between the magnitude of the waves sent into the body and the magnitude of the waves returned back near their point of origin (referred to herein as "return loss") is utilized to estimate certain characteristics of biological tissue (e.g., cardiac tissues) within the body. Additional details regarding the signal propagating device and method of estimating cardiac output from data the device gathers are presented in turn below. The present invention provides several significant advantages over prior related medical device. Many advantages will be apparent in light of the detailed description set forth below, with reference to the accompanying drawings. These advantages are not meant to be limiting in any way. Indeed, one skilled in the art will appreciate that other advantages may be realized, other than those specifically recited herein, upon practicing the present invention. RF DEVICE
With reference now to FIGS. 1 through 3, and in accordance with one embodiment of the invention, a device 10 for estimating changes in the cardiac output of a patient is disclosed comprising a radio frequency ("RF") signal generator coupled to an antenna 15, said antenna 15 having a dielectric material 20 disposed about an exterior surface of the antenna 15. The dielectric material 20 may be disposed about select portions of the exterior surface of the antenna 15 or substantially about the entire exterior surface of the antenna 15.
The RF signal has a near field component and a far field component. In one aspect of the invention, the dielectric material 20 has a predetermined thickness and a predetermined dielectric constant such that a substantial portion of the near field component of the electronic signal returned from the patient is not received by the antenna 15. In one aspect of the invention, the dielectric material 20 is designed such that substantially the entire near field component of the electronic signal returned from the patient is not received by the antenna 15. In one embodiment, the dielectric material 20 has a dielectric constant ranging from 5 to 25 and has a thickness ranging from 2 mm to 20 mm. Importantly, the design of the dielectric material 20 is dependant upon numerous characteristics which vary depending on the type of dielectric material used. That is, the thickness of the material is a function of its dielectric constant and thus it is specific for a given material. Accordingly, no limitation should be read into the exact specifications of the dielectric material used herein.
With reference to FIG. 2, a block diagram of the RF device with digital detection means in accordance with one embodiment of the present invention is illustrated. In one aspect of the invention, the RF signal is generated by a voltage controlled oscillator 25 ("VCO") contained within a device housing 11. The VCO 25 generates the RF signal at differing frequencies by varying the tuning voltage. Examples of VCO's 25 which may be used include the ROS-2150 made by Minicircuits, however any VCO 25 capable of generating an appropriate RF signal may be utilized. As noted herein in more detail, initially, a micro-controller 30 via a digital to analog ("D/A") converter 31 sweeps the frequency of the RF signal between approximately 1000 and 2000 MHz by causing the tuning voltage to vary between approximately zero and 20 volts. In one aspect of the invention, the output of the VCO 25 is less than 4 mW. While any suitable operating frequency may be used, in one embodiment of the present invention, the VCO 25 operates at a frequency ranging from 0.5 GHz to 2.5 GHz.
In one aspect of the invention, the device 10 further comprises a network analyzer 12 such as that manufactured by Analog Devices (AD8302 IC). The detector is capable of identifying the magnitude and phase of subject signals and has a 60 dB dynamic range and outputs a magnitude voltage of 30 mV per dB of return loss. The phase output is 10 mV per degree over 180 degrees. The device 10 further comprises microcontroller 30 (such as those manufactured by Cygnal) to read the network analyzer outputs, determine the minimum or maximum values of magnitude and phase during a frequency sweep and generate a digital sweep that enters the D/A converter 31 that then drives the VCO 25. This closed loop generates the VCO sweep, determines the frequency of maximum return loss (magnitude or phase), and then locks on to that frequency. As the heart beats or lungs move, the microcontroller 30 tracks the magnitude and phase of the "organ's" signal. In one aspect of the invention, a low pass filter 35 ("LPF") is used to prevent the
VCO second harmonic from entering the network analyzer 12 (e.g., AD8302 magnitude and phase integrated circuit). In one aspect, the LPF 35 has a -3 dB "corner" frequency of 3700 MHz. In an additional aspect of the invention, the LPF output is routed to two directional couplers 38. The directional couplers 38, 39 sample the RF signal both in the forward (38) and reverse (39) directions. In one embodiment, the sampled signal is 20 dB down from the transmitted signal. In an additional aspect, directional couplers 38, 39 exhibit directivity. That is, they are relatively insensitive to signals coming in a reverse direction. The directional couplers 38, 39 noted herein have directivities of 20 to 34 dB depending on the frequency of the RF signal. The transmitted RF energy goes to the antenna 15 where the coupled signals are routed to the network analyzer 12. A voltage multiplier 32 and appropriately placed attenuators 33 may also be used in accordance with certain aspects of the invention. In one aspect, appropriate programming control and display interfaces are incorporated on the face of the RF device. In other aspects, the programming control and display interfaces are located on a remote programming device. With reference now to FIG. 3, a block diagram of the RF device with analog detection means in accordance with one embodiment of the present invention is illustrated. While similar to the digital detection system shown in FIG. 2, additional components are shown to achieve proper detection of an analog signal. Additional components include a sawtooth generator 36, a generator trigger signal 37 and an amplifier/signal conditioner/peak detector/analog digital converter device 41.
Referring generally again to FIGS. 1 through 3, the RF components discussed herein are fitted to a 4 layer PCB with 50-ohm traces. The RF board components are shielded and, in some embodiments, feed through capacitors are used to prevent RF signals from "leaking out" of the shield via the DC or low frequency signals. In accordance with one aspect of the invention, the device is capable of measuring the signal transmitted from the VCO 25 to an antenna 15 where the signal has been modulated by the organ system of the patient. In this aspect, the device 10 comprises an additional dipole antenna 15a, an additional LPF 35a, a single directional coupler and a second network analyzer 12a, one of whose inputs is from the VCO 25. Output from the network analyzer 12a is routed to the data collection system.
In one embodiment of the present invention, the device further comprises a switch sensor for activating and deactivating the device 10. The switch sensor would activate the device 10 only when it is in full contact with the body of the patient and would deactivate the device 10 when not in full contact with the body. Advantageously, errant RF radiation emitted into the surrounding environment would be minimized by using the switch sensor. The switch sensor may comprise a touch-activated switch sensor or other desirable switch sensor as suits a particular application.
Coupling of the RF signal into the body is maximized when there is a good impedance match between the coupler and the body. It is believed that biological tissues present about a 50-ohm load for the RF signal. Appropriately placed dipoles and loop antennas provide as much as 50 dB return loss at frequencies when placed near the body, even through clothing. As the antenna 15 nears the body, the tissue loads the antenna 15 changing its resonance to below that of free-space resonance. In accordance with one aspect of the invention, the measurement of the return loss is performed at 1 millisecond intervals, resulting in I kilosample per second signal acquisition rate. As noted above, the core of the system is the network analyzer 12 which measures the ratio of the returned RF signal to the sampled outgoing signal generated by the VCO 25. The optimal frequency of the VCO 25 is found by scanning the frequency range within approximately 1-2 GHz range and determining the frequency of the maximum return loss. The magnitude, phase, and frequency at the maximum return loss point is then digitized and displayed and/or stored in the data file. Additional functions of the device included acquisition of up to 4 analog channels (e.g., flow, impedance, ECG, other) of 0-5V range and synchronizing the analog data channels with the RF magnitude, phase and frequency. In accordance with one aspect of the invention, a dipole antenna is utilized which comprises copper foil tape applied to a suitable backing substrate. A 50-ohm coaxial cable connects the dipole antenna to a processor.
In accordance with an additional embodiment of the present invention, a full wave loop antenna placed in a cylindrical shield may be utilized. As is the case with many biomedical sensors, the output of the sensor (i.e., the antenna) is affected by multiple factors giving rise to common mode error, such as sensor-body interface, relative motion of the sensor and the tissues at the sensor placement site, RF absorption and reflection by internal body organs other than the heart, and by motion artifacts generated by both intrinsic (e.g., breathing) and external factors (e.g., sensor pressure, sensor motion, etc.). Interfacial motion generates signal components that are much larger than the signals related to the heart and therefore it is essential to provide signal conditioning and data processing that can separate the signal of interest. That is, signals that relate to the stroke volume of the heart need to be measured while also preserving the breathing and respiratory components of the measured signal, such as breathing rate and depth. The difficulty of achieving these tasks is compounded by small levels of the RF signal from the sensor and the resulting low Signal to Noise ratio (SNR), and by the high ratio of the DC to AC component of the return loss signal.
Generally speaking, interactions of electromagnetic fields with matter are a function of frequency, dielectric, conductive and magnetic properties of the material, geometry of the radiating element (antenna) and distance, and orientation of the antenna with respect to the observation point. Calculations of an antenna field have been extensively studied for different antenna configurations, sizes and materials, and fundamental equations approximating electric and magnetic fields for simple antenna geometries and homogeneous media, are available. However, in cases of either complex antenna shapes or non-homogeneous media, calculations are very difficult, if possible. Therefore, while the following analysis is an oversimplification of the actual setting, it provides significant insights into the origin of motion artifact and offers a practical way of reducing the influence of undesirable interference. In one aspect of the present invention, a circular loop antenna was employed for which the components of the electric and magnetic field as a function of distance from the antenna are given by numerous equations. These equations represent fields resulting from a sinusoidally-excited current loop with a magnetic moment "dm". Electric and magnetic components of the field are given by the relationships:
Eφ, =30β3 * dm * [1 / (β r) - j / (β r)2 ] sin(θ ) e( -j(β r)
The electric field contains components are proportional to the 1 / (β r) and 1 / (βr)2 and magnetic field components which are a function of 1 / (β r)3 where r is a distance from the radiating element, and (β is a wave number, (β = 2Pi / λ). Therefore, the behavior of the complex electric field as a function of distance is determined by the functions: 1 / (βr), and 1 / (βr)2, 1 / (β r)3 for magnetic components)
If the boundary between near and far electrical field is defined as a distance at which r / λ = I / 2Pi, then at that distance the contributions of each of the terms l/( βr)n becomes equal (i.e., for r = λ / 2Pi, 1 / (βr) ~ 1 / (βr)2). In the near field, the higher power term, 1 / (βr)2 dominates the complex field components, while beyond that limit, the contribution of this term becomes insignificant and the electric field is determined by a familiar relationship: E ~ 1 / r. Thus, it is believed that any interactions between the radiated EM wave and an object will be stronger in the near field region than in the far field. FIG. 4 illustrates the calculated contributions to the complex electric field, of the 1 / (βr) n terms (curves 1, 2) in the near field region of the circular loop antenna in the material with the relative dielectric constant εr = 15. In the material with such high dielectric constant, the wavelength, 1, at 1.5GHz is only about 5cm (λ=c/f * (l/sqrt(εr)), where c is the speed of light and f is the frequency. Curve (1+2) represents the sum contribution from the 1 / (βr)n components.
The above analysis suggests that at the frequency range of interest all the interactions between the antenna and the body surface of the patient occur within the near-field distance, and therefore tend to be dominated by higher-power terms of the complex field, while the impact of the heart (that is positioned under the tissues in the distance of more than 2 cm from the surface) is dominated by the 1 / (βr) term, which varies much slower than the quadratic term. Accordingly, it is believed that a circular loop antenna would be expected to exhibit much higher sensitivity to interfacial motion than to the motion/absorption at the tissues and organs at distances further away from the antenna surface, especially considering that the wavelength inside the body is even shorter (average dielectric constant of the tissues is about 50, therefore λ~ 2.4cm).
Accordingly, it is believed that if the near-field region is contained within the stationary (not subject to interfacial perturbations) antenna limit, and all the interactions between the antenna and the body could be limited to the far-field region where the complex field is dominated by the 1/r term only, the return loss sensitivity would be the same for all the interactions contributing to the return loss signal of the sensor. That is, the relative contribution from the surface motion artifact to the signal would be greatly reduced. As such, by disposing a dielectric material about the antenna, motion artifacts are greatly reduced.
In one embodiment of the present invention, ceramic disks of 8.5 mm and 5.5 mm thickness and εr =15 are attached to the front surface of the antenna 15. The results of tests of monitoring of the heart and breathing action in human volunteers in the presence and absence of the interfacial motion are shown in FIG. 5. Advantageously, the output signal of the antenna 15 having the dielectric material 20 was attenuated as compared to the signal measured with the antenna 15 without the dielectric material 20. However, the signal was much less sensitive to motion artifacts induced by moving the antenna 15 with respect to the surface of the patient. While the bandwidth of the signal detected by the device 10 after passing through the dielectric material 20 is reduced, thereby making the system less efficient, without the dielectric material 20, the motion artifact signal entirely obliterates the heart beat signal. For the frequencies used in the testing of the device 10 (e.g., 1000 MHz to 2000
MHz), the motion artifact-related signal originated in the near-field of the antenna 15 (i.e., at the interface between the antenna and the body). In this region, complex electric field modulation is dominated by fast changing 1/r2 function, whereas in the far field, these interactions result in a slower 1/r relationship. Placing a dielectric material 20 on the surface of the antenna 15 causes shortening of the wavelength in the dielectric region by a factor of l/(sqrt(εr)). Thus, a proper selection of the dielectric material 20 and its geometry (thickness), permits enclosing the near-field region entirely in the dielectric material, resulting in the significant reduction of the device sensitivity to interfacial motion as noted above.
While specific reference has been made herein to ceramic as a dielectric material 20, any desirable dielectric material may be utilized as suits a particular application. Example dielectric materials include ceramics, glass, plastics, and polymers. Fluid dielectric materials may also be used as desired including, but not limited to gels, liquids, polymeric fluids, and hydrocarbon-based fluids. Additionally, while specific reference has been made to a specific thickness, it is important to note that the thickness of the dielectric material 20 will vary according to the type of dielectric material used.
In an additional embodiment of the present invention, additional modifications to the antenna 15 were made to minimize motion artifact. As has been discussed herein, the antenna 15 is very sensitive to tangential (lateral), as well as normal motion of the antenna 15 with reference to the skin surface. The sensitivity includes changes in the force of the sensor on a surface of the body. These factors affected both the initial sensor coupling to the body (DC loading level, dB), and to the lesser extent, an AC component of the signal.
As noted above and in accordance with one embodiment of the present invention, a pressure sensor is coupled to the device 10 to detect the application force and motion of the antenna 15 with reference to the body. Sensor software can be created to permit valid data measurements only while the sensor is balanced and indicates "no-motion" conditions. In one aspect of the invention, pressure sensors are disposed on all four corners of the antenna 15 and are connected to a quad bridge circuit. In one embodiment of the antenna design, three strain gages are built into the applicator shield and used in the same manner. Software is programmed into the network analyzer 12 to visually show the strain on each gage and prevent measurement if the difference in strain exceeds an adjustable threshold. An indicator light is disposed on the device 10 so the user knows when the strain is equal on all three gages thus indicating the applicator is evenly positioned. An example force sensor-based stabilization signal is shown in FIG. 6. At the left of the graph, the sensor applicator is pressed against the chest and the variation in the sensor signal is illustrated. Towards the center of the graph, applicator motion becomes minimal and the sensor is stabilized. On the right hand side of the graph, sensor motion is again occurring and as such, data should not be recorded. Illustrated further in FIG. 5, are example sensor signals with and without motion artifact reduction. The signal at the left of the top panel graph is off-scale and cannot be interpreted while suppressing the motion-free signal segment. The bottom panel illustrates the result of sensor signal measurement with motion artifact reduction design.
SIGNAL PROCESSING AND NON-INVASIVE METHOD OF MEASURING CARDIAC OUTPUT In accordance with one embodiment of the present invention and as described in more detail through this application, a method for estimating certain parameters of a person's cardiac output is disclosed. Broadly speaking, the method comprises placing a device 10, such as that described in the preceding sections of this application, about the mid-sternal position of the person. The device 10 comprises an RF signal generator 25 coupled to an antenna 15, said antenna 15 having a dielectric material 20 disposed about the exterior 16 of the antenna 15. The method further comprises propagating a first signal having a predetermined frequency towards the heart of the person and receiving and measuring a portion of the first signal returned from the heart of the person with the device 10. The magnitude of the first signal propagated into the heart of the person is compared to the magnitude of the portion of the first signal returned from the heart of the person and a return loss of the signal is calculated. The method further comprises propagating an additional signal into a portion of the heart of the person using the device, wherein said additional signal has a frequency different from the first signal and receiving and measuring a portion of the additional signal returned from the heart of the person with the device 10. The magnitude of the additional signal propagated into the heart to the person is then compared to the magnitude of the additional signal returned from the heart of the person and a return loss of the additional signal is calculated. The method further comprises comparing the return loss of the first signal to the return loss of the additional signal and repeating certain of the above-referenced steps while incrementally varying the frequency of the signal with each repetitive step. The method also comprises determining a maximum return loss value of the signals propagated into the heart, of the person over the utilized frequency range. As noted above, in one embodiment, the frequency of the signal is swept between 1000 and 2000 MHz in an effort to determined the maximum return loss for a particular patient at a particular frequency.
It is believed that proper positioning of the device 10 on the patient is important for optimal operation of the device. In one aspect of the invention, optimal results are obtained when the device 10 is positioned in the mid-sternal position. That is, the optimal results are obtained when the signals from the device are predominantly focused on the right ventricle of the patient's heart.
In an effort to understand the physiological interpretation of the RF device signal, a time-domain analysis of the RF signal was recorded synchronously with a 3-lead ECG (RA-LA-LL). A comparison of that analysis is explained in more detail below and illustrated in the attached figures. Generally speaking, the cardiac cycle consists of two major periods, systole, during which the ventricles contract and, under high pressure, eject blood into the pulmonary artery and into the aorta, and diastole, during which the heart is in its low pressure state and the ventricles are filling with blood. At the beginning of the systole, the heart is maximally extended, while at the beginning of the diastole the extension of the heart muscle is minimal. The systole starts with the R-wave, which is immediately followed by contraction of the left ventricle, and ends after the T-wave with the closing of the aortic valve. According to Frank-Starling Law of the heart, tension developed in cardiac muscle fibers is proportional to their stretching and thinning, or to the volume of blood filling the heart, or End Diastolic Volume ("EDV"). The EDV is one of the factors that determine the size of the subsequent stroke volume. If a signal is proportional to the volume of the right ventricle, the magnitude of the signal should vary between cardiac cycles with the venous return, or preload. Therefore, the device acts as a sensitive indicator of the physiological events that affect venous return. Immediately following the R-wave, both tricuspid and pulmonic valves are closed, and the initial ventricular contraction results in an isovolumetric intraventricular pressure build-up. With general reference now to FIG. 7, based on the comparison of the left ventricular pressure wave and S-wave timing, components of a recorded RF signal can be identified that correspond to mechanical contraction of the heart and thus provide information about changes in the cardiac output. In one aspect of the invention, an analysis begins with the R-wave, which initiates contraction of the ventricles. Upon contraction, with some delay related to spreading of the electrical wave, the ventricles contract isovolumetrically (both valves are still closed). During isovolumetric contraction, the RS trade rises initially slowly, then rapidly, indicating a thickness change of the ventricular wall as detected by the antenna, consistent with contraction of the left ventricle. When the LVP exceeds the End Diastolic Pressure ("EDP") in the aorta, the aortic valve opens and blood is ejected from the left ventricle into the aorta, initially rapidly, and then slowly as the aortic pressure rises to the End Systolic Level, resulting in reduced low. During this period, the left ventricle reaches its minimum volume (maximum point on the RS curve), and starts expanding. Also during this period, the RS reaches its maximum. That is, the heart wall has reached its maximum thickness. As the LVP declines, the aortic valve closes. The left ventricle continues to relax, with both valves closed (isovolumic relaxation), resulting in reduced thickness of the ventricular wall, which appears on the RS as a section with significantly reduced slope. As the ventricular pressure drops below the level of atrial pressure, the mitral valve opens, initiating the period of rapid filling. The ventricular wall continues thinning as blood fills the ventricle. Both mitral and aortic valves are closed, and the left ventricle continues its (isovolumic) relaxation. When the LVP drops below the atrial pressure, the mitral valve opens, and the ventricle fills with blood, expanding its size and causing a thinning of the wall. The minimum of the RS curve corresponds to a filled ventricle ready for ejection following atrial contraction (P-wave). Following the P-wave, the mitral valve closes, and the cycle begins again. Similar pressure - valve opening/closing relationships are also true for the right atrium, right ventricle closing and opening of the pulmonic valve during systole.
Based on this relationship between data points shown in FIG. 7, it is believed that (i) the positive slope indicates increasing return loss (less return of the RF signal), (ii) the section curve corresponding to isovolumic relaxation is either more (1st trace) or less (last trace) pronounced, as the filling of the heart is influenced by the intrathoracic pressure during the breaching cycle, (iii) the negative slope of the return loss curve indicates decreasing return loss, (iv) the section of the Recorded Signal trace corresponding to the filling cycle is flat due to atrial "pull," and (v) the top section of the left intraventricular pressure curve is flat due to saturation of the LVP sensor. As a result, the correlation of return loss signal and movement of the ventricular wall enables estimation of the magnitude of movement of the heart wall and thereafter information regarding changes in stroke volume and cardiac output.
SIGNAL MODELING
According to one embodiment of the present invention, the device described in detail above utilizes the interaction of electromagnetic energy with biological tissues. The RF device utilizes the interaction of electromagnetic energy with biological tissues. Generally speaking, biological tissues are characterized by their conductivity and dielectric constant. It is believed that the magnetic permeability (μ) of tissues is similar to that of free space, and thus it is assumed that μ = 1. These properties are a function of frequency and tissue type. The FCC has compiled a database of these properties at http://www.fcc.gov/fcc-bin/dielec.sh and this database was used extensively throughout the research and development of the RF device described herein.
The interaction of electromagnetic ("EM") energy with materials can be described by several important equations. The equation λ = α + jβ where gamma is the complex propagation constant, alpha is the attenuation constant and beta is the phase constant is helpful for understanding modeling of the device signals. It is important to note that these "constants" are constant for only one frequency and particular material. That is, as the signal frequency is varied and the subject tissue is changed, the constants are variable. Beta is equal to 2*pi/lambda, the wavelength in the material. Alpha has the units of Nepers/meter. For lossy dielectrics the propagation constant can be written as λ =j oW(με)* V(l-j(σ/ωε)) where omega is the radian frequency, μ is the permeability of free space, epsilon is the permittivity of the material and sigma is the material's conductivity (Siemens/meter). The intrinsic impedance of a lossy dielectric is given by:
Figure imgf000019_0001
EM energy is attenuated by lossy dielectrics and returned at impedance boundaries. A useful equation is the transmission line equation that describes the impedance seen by an EM wave as it traverses materials of differing impedance. Thus the impedance looking into a series of slabs of material is: ηin = ηι((η2+jηi*tanpi)/(η,+jη2*tanβl)) where 1 is the thickness of material 1. To estimate the amount of the signal returned from the muscle/blood interface, the intrinsic impedances of muscle and blood are estimated to be: muscle (49.2+j 5.91) and blood (46.5+J7.69) (note that the real part of the impedances is very close to 50 ohms so the tissues present a good match to 50 ohm cable). To estimate the amount of energy that is returned from the blood/muscle interface it is assumed that P=(Zi-Zo)/ (Zj +Zo) where Zi and Zo are the complex impedances of the load (blood) and line (muscle). Thus p = (- 0.0247, γθ.022) and about 2.5% of the energy is reflected at the muscle blood interface. There is also a change in the phase of the returned wave referenced to the incoming wave of 0.022 radians. These background assumptions allow for the construction of more complex models where there are multiple slabs including air layers representing different tissue layers, such as heart wall, heart ventricle, etc.
Of the possible approaches to examine the coupling of RF EM energy into the body, two are discussed below. First, the body is a load presented to a transmission line via an "antenna coupling device" or second, an EM wave impinges on the body in the antenna's near field. Such an approach is an oversimplification but provides a qualitative assessment of the interaction of the applicator with the body.
One objective of developing a simplified propagation model for calculating RF properties of the antenna as a function of cardiac output and presence/absence of hemo/pneumothorax was to determine whether, given the dielectric properties of the tissues at the frequency range between 1 and 2GHz, using a simplified 2-D slab configuration.
The model is based on the interaction of a planar electromagnetic wave with the body tissues modeled as slabs and has not been designed to provide predictive information, but to assist in understanding the behavior of the reflection coefficient as a function of the changes in thickness of muscle, blood, and air and thus help in selecting the frequency of the measurement and characteristic waveform of the recorded sensor signal as a function of the beating heart. It is believed that the model assists in assessing the mechanisms responsible for modulation of the return loss signal as a function of either ventricular volume and of the movement of the ventricular wall. Similarly, introduction of an air pocket, as in pneumothorax, would create an impedance change that may be predicted by the model and therefore measurable by a suitable RF device. Selection of the working frequency impacts the depth of penetration of the RF radiation into the tissues. FIG. 8 is a schematic of the anticipated depth of penetration with the range of frequencies of interest to this project, based on the literature data. The model assumes that the body is represented by a sequence of layers, corresponding to chest wall, heart muscle, blood and bone. Based on a lumped parameter model of a transmission line, complex electrical impedance was calculated for each layer, and an impedance change (reflection) was determined at each inter layer interface. According to this model, each layer, corresponding to different tissue forms a load for the EM wave traveling through the immediately preceding layer of tissue, and calculations are completed iteratively for all the layers of tissues, over a range of frequencies and dimensions.
The reflection coefficient in a lossy terminated transmission line is given by:
Figure imgf000021_0001
where Zi is the complex impedance of the load, and Zo is the complex input impedance "looking into" the line. For example, for the first layer, Zi is the characteristic impedance of the chest wall, which is assumed to be a 2cm thick layer with dielectric properties corresponding to those of the muscle, bone and fat in 1 : 1 : 1 proportion. The input impedance at each interface can be computed using the formula:
Z1n = Zo*[(Z1+j*Zo*tanh(λl))/Zo+j*Z1*tanh(λl)] and tanh(λl) = (sinh(αl)*cosh(βl) + j*cosh(αl)*sinh(βl))/ (cosh(αl)*cosh(βx) +j*sinh(αl)*sinh(βl))
Where β = is the phase constant parameter representing lossy (imaginary) part of the complex propagation constant, λ = α + jβ and α is the attenuation constant. The intrinsic impedance is calculated from the formula:
Z1 = SQRT (μ/εrεo)*[l/SQRT(l-jσ/ω εrε0))]
FIG. 9 is an illustration of a frequency tuning sweep and resulting work point in accordance with one embodiment of the present invention. ANIMAL STUDIES
In the course of developing certain embodiments of the present invention, animal testing was conducted to establish metrological correlation between the measured sensor signal and the stroke volume/cardiac output tests that are impossible, or at best extremely dangerous to perform on human subjects. The animal studies included (i) measurements of cardiac output in the pulmonary artery during hemorrhage, and (ii) measurements of cardiac output using the device described herein as a function of blood loss during controlled hemorrhage and in correlation with the Pulmonary Artery Flow ("PAF"). Sensor readings were recorded during: (i) a control period prior to hemorrhage, (ii) bleeding and re-infusion of blood at a rate of 100 ml/min, up to approximately thirty-three percent of the estimated blood volume of the animal, and (iii) terminal bleeding. The following experimental variables were continuously recorded in a personal computer: (i) ECG using 3-lead standard ECG configuration, (H) breathing, (iii) sensor output, AC- coupled, (iv) sensor output, DC-coupled, (v) PA flow, instantaneous, (vi) PA flow, average, and (vii) frequency (i.e., the optimized maximum return loss frequency).
In order to separate the breathing-related baseline variations, selected segments of data are digitally filtered. In one aspect of the invention, the sensor output signal contains two components: an additive, breathing and motion-modulated baseline at breathing and motion-artifact frequencies, and the stroke volume-related signal at heart beat frequency. The separation of these components is performed in the frequency domain, and is followed by subsequent interpretation of the filtered signal in time domain.
Generally speaking, the digital filtering discussed herein includes the following steps: (i) calculating the FFT of the selected data segments, (ii) plotting the power spectrum, (iii) dividing the power spectrum of the signal into two separate spectra containing the breathing-related and (separately) the heart-beat-related frequency components of the original signal, (iv) re-creating the heart-related and breathing-related waveforms using inverse FFT function, and (v) plotting the peak amplitudes of the sensor signal versus the corresponding peak amplitudes of the PA flow or volume bled/re- infused. In one aspect of the invention, the frequency bandwidth used for the analysis is
3Hz (adjustable), and the breathing rate separation frequency is 0.75Hz, or 25% of frequency bandwidth. This routine resulted in zero-average signals that showed multiplicative (amplitude) modulation by breathing. This digital filter was implemented in MatLab and the calculations were performed off-line. However, this filter routine can be implemented for processing of the RF sensor data in real time on an embedded platform. This method of filtering the baseline drift and its additive contribution to the sensor signal is not adaptive and does not account for inter-subject variability. However, it provides a simple and effective method of separating the breathing and low frequency motion interferants. In dynamic situations, when static selection of filtration boundaries is not possible, the limits are adjusted dynamically as a function of heart and breathing frequencies.
FIG. 10 shows one implementation of the above-described filter using MatLab data processing software for a selected data segment. Specifically, FIG. 10 shows a data plot of the sensor magnitude versus pulmonary artery stroke volume from experiment SW3 in accordance with one embodiment of the present invention. FIG. 11 provides a summary data from the bleed experiment in SW3 experiment.
High correlation between the returned RF device signal with bled volume is practically the same as for the pulmonary artery stroke volume signal, indicating that the RF device can become a non-invasive, simple to use diagnostic tool for detection of hemorrhage. FIG. 12 (upper panel) shows a bleed re-infusion cycle in the animal of SW3. Stroke volume derived by integrating PA instantaneous flow is plotted versus blood loss/gain. From the same data record, the bottom panel shows the sensor signal plotted similarly. FIG. 13 provides correlation graph for pooled unsealed PEEP08 and PEEPl 4 data. The correlation coefficient, R — 0.91 16, indicated a strong correlation between stroke volume derived from COTB measurements, and that from the RF sensor signal. The foregoing detailed description describes the invention with reference to specific exemplary embodiments. However, it will be appreciated that various modifications and changes can be made without departing from the scope of the present invention as set forth in the appended claims. The detailed description and accompanying drawings are to be regarded as merely illustrative, rather than as restrictive, and all such modifications or changes, if any, are intended to fall within the scope of the present invention as described and set forth herein. More specifically, while illustrative exemplary embodiments of the invention have been described herein, the present invention is not limited to these embodiments, but includes any and all embodiments having modifications, omissions, combinations (e.g., of aspects across various embodiments), adaptations and/or alterations as would be appreciated by those in the art based on the foregoing detailed description. The limitations in the claims are to be interpreted broadly based on the language employed in the claims and not limited to examples described in the foregoing detailed description or during the prosecution of the application, which examples are to be construed as nonexclusive. For example, in the present disclosure, the term "preferably" is non-exclusive where it is intended to mean "preferably, but not limited to." Any steps recited in any method or process claims may be executed in any order and are not limited to the order presented in the claims. Means-plus-function or step-plus-function limitations will only be employed where for a specific claim limitation all of the following conditions are present in that limitation: a) "means for" or "step for" is expressly recited; and b) a corresponding function is expressly recited. The structure, material or acts that support the means-plus function are expressly recited in the description herein. Accordingly, the scope of the invention should be determined solely by the appended claims and their legal equivalents, rather than by the descriptions and examples given above. What is claimed and desired to be secured by Letters Patent is:

Claims

1. A method for measuring cardiac output of a person, comprising:
(a) placing a device about the center of the chest of the person at the mid-sternum position, said device comprising a RF signal generator coupled to an antenna, said antenna having a dielectric material disposed about the exterior of the antenna;
(b) propagating a first signal having a predetermined frequency towards the heart of the person;
(c) receiving and measuring a portion of the first signal returned from the heart of the person with the device; (d) comparing the magnitude of the first signal propagated into the heart of the person to the magnitude of the portion of the first signal returned from the heart of the person and calculating a return loss of the signal;
(e) propagating an additional signal into a portion of the heart of the person using the device, wherein said additional signal has a frequency different from the first signal; (f) receiving and measuring a portion of the signal returned from the heart of the person with the device;
(g) comparing the magnitude of the signal propagated into the heart to the person to the magnitude of the additional signal returned from the heart of the person and calculating a return loss of the additional signal; (h) comparing the return loss of the first signal to the return loss of the additional signal;
(i) repeating steps (e) through (h) while incrementally varying the frequency of the signal with each repetitive step; and
Q) determining a maximum return loss value of the signals propagated into the heart of the person over the utilized frequency range: and
(k) determining the frequency of the maximum return loss signal.
2. The method of claim 1, wherein the signal comprises a radio frequency signal.
3. The method of claim 1, wherein the dielectric material is selected from the group consisting of ceramics, glass, liquids, gels and other polymeric materials.
4. The method of claim 1 , further comprising estimating the amplitude of motion of the portion of a wall of the heart.
5. The method of claim 4, further comprising estimating a change in the amplitude of motion of the portion of the wall of the heart.
6. The method of claim 5, further comprising estimating a stroke volume of the heart based on the change in the amplitude of motion of the portion of the wall of the heart.
7. The method of claim 5, further comprising determining the timing of a QRS complex of the heart based on the change in the amplitude of motion of the wall of the heart.
8. The method of claim 5, further comprising determining the timing of a QRS complex of the heart based on the change in frequency of the maximum return loss.
9. The method of claim 1, wherein the portion of the heart comprises a right ventricle of the heart.
10. A system for measuring the cardiac output of the heart of a person, comprising: a device comprising a housing containing a signal generator coupled to an antenna and a dielectric material disposed about the antenna, wherein said hand-held device is adapted to generate a plurality of radio frequency signals and direct said signals towards the heart of the person and measure a magnitude of a signal returned off the heart of the person; a processor contained within said housing, said processor being adapted to compare differences between a magnitude of a signal propagated into the heart of the person and the magnitude of the signal reflected from the heart, said processor being further adapted to determine a signal frequency having a maximum return loss value based on the differences between the magnitude of the signal propagated into the heart of the person and the magnitude of the signal returned from the heart of the person; and said processor being further adapted to estimate a change in the amplitude of motion of a portion of a wall of the heart based on the differences between the magnitude of the signal propagated by the device and the magnitude of the signal returned off of the portion of the heart.
11. The system of claim 10, wherein said processor is further adapted to calculate the cardiac output of the person; and further comprising a display disposed about the handheld device adapted to display the measured signals and that can be calibrated to display the estimated stroke volume of the person.
12. The system of claim 11, wherein said processor is further adapted to calculate the cardiac output of the person based on an estimated stroke volume of the person.
13. The system of claim 10, wherein said device is adapted to generate a signal at a frequency and magnitude having a determined maximum return loss value into the portion of the heart and receive a portion of the signal propagated into the portion of the heart.
14. The system of claim 10, wherein the dielectric material is selected from the group consisting of ceramics, glass, liquids, gels and other polymeric materials..
15. The system of claim 12, wherein said processor is further adapted to estimate the stroke volume of the heart based on a change in of the amplitude of motion of the portion of the wall of the heart.
16. The system of claim 15, wherein said processor is adapted to further comprise the step of estimating a QRS complex of the heart based on the change of the amplitude of motion of the portion of the wall of the heart.
17. The system of claim 15, wherein said processor is further adapted to further comprise the step of estimating a QRS complex of the heart based on the change of the frequency of the maximum return loss of the signal propagated towards the portion of the wall of the heart.
18. The system of claim 10, wherein the dielectric material has a predetermined thickness and a predetermined dielectric constant such that a substantial portion of a near field component of the signal propagated towards the heart of the person is contained within said dielectric material.
19. The system of claim 10, wherein the dielectric material has a predetermined thickness and a predetermined dielectric constant such that substantially an entire near field component of the signal propagated towards the heart of the person is contained within said dielectric material.
PCT/US2009/030157 2008-01-04 2009-01-05 Non-invasive method and device for measuring cardiac output WO2009089183A1 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
EP09700309A EP2237719A1 (en) 2008-01-04 2009-01-05 Non-invasive method and device for measuring cardiac output
CN200980107669.8A CN101983030B (en) 2008-01-04 2009-01-05 Non-invasive method and device for measuring cardiac output
JP2010541585A JP5406211B2 (en) 2008-01-04 2009-01-05 Non-invasive method and non-invasive device for measuring cardiac output
MX2010007414A MX2010007414A (en) 2008-01-04 2009-01-05 Non-invasive method and device for measuring cardiac output.

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US1003508P 2008-01-04 2008-01-04
US61/010,035 2008-01-04

Publications (1)

Publication Number Publication Date
WO2009089183A1 true WO2009089183A1 (en) 2009-07-16

Family

ID=40473708

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2009/030157 WO2009089183A1 (en) 2008-01-04 2009-01-05 Non-invasive method and device for measuring cardiac output

Country Status (6)

Country Link
US (1) US8721559B2 (en)
EP (1) EP2237719A1 (en)
JP (1) JP5406211B2 (en)
CN (1) CN101983030B (en)
MX (1) MX2010007414A (en)
WO (1) WO2009089183A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013530741A (en) * 2010-05-28 2013-08-01 ザ ジェネラル ホスピタル コーポレーション Apparatus, system, method and computer-accessible medium for analyzing information on cardiovascular disease and cardiac function
US11457889B2 (en) 2021-01-25 2022-10-04 Lazaro Eduardo Hernandez System and method of non-invasive continuous echocardiographic monitoring

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8463361B2 (en) * 2007-05-24 2013-06-11 Lifewave, Inc. System and method for non-invasive instantaneous and continuous measurement of cardiac chamber volume
US8264342B2 (en) 2008-10-28 2012-09-11 RF Surgical Systems, Inc Method and apparatus to detect transponder tagged objects, for example during medical procedures
US8726911B2 (en) 2008-10-28 2014-05-20 Rf Surgical Systems, Inc. Wirelessly detectable objects for use in medical procedures and methods of making same
US9792408B2 (en) * 2009-07-02 2017-10-17 Covidien Lp Method and apparatus to detect transponder tagged objects and to communicate with medical telemetry devices, for example during medical procedures
US9226686B2 (en) 2009-11-23 2016-01-05 Rf Surgical Systems, Inc. Method and apparatus to account for transponder tagged objects used during medical procedures
US20120123231A1 (en) * 2010-11-11 2012-05-17 O'reilly Michael Monitoring cardiac output and vessel fluid volume
EP2851100A1 (en) * 2013-09-20 2015-03-25 Berlin Heart GmbH Blood pump control system and method for controlling a blood pump
CN104138255A (en) * 2013-12-24 2014-11-12 智群(香港)有限公司 Belt type heart electromagnetic wave collecting device and manufacturing method thereof
CN104382567A (en) * 2014-11-14 2015-03-04 电子科技大学 Near infrared spectrum human hemodynamics detection device and motion interference elimination method
US10687729B2 (en) * 2014-11-24 2020-06-23 Koninklijke Philips N.V. Apparatus and method for estimating a value of a physiological characteristic
US10856806B2 (en) * 2015-02-12 2020-12-08 University Of Hawaii Lung water content measurement system and calibration method
WO2017008020A1 (en) 2015-07-08 2017-01-12 The Johns Hopkins University Tissue ablation and assessment system and method of use thereof
US20210361180A1 (en) * 2015-11-16 2021-11-25 Jun-Chau Chien Gigahertz frequency fringing near-field biomedical sensor
US11083388B2 (en) * 2015-11-16 2021-08-10 Jun-Chau Chien Gigahertz frequency fringing near-field cardiovascular sensor
WO2018194093A1 (en) * 2017-04-19 2018-10-25 学校法人関西大学 Biological information estimation device
CN111031903B (en) 2017-06-16 2023-08-04 康奈尔大学 Method and system for electromagnetic near-field coherent sensing
US11392782B2 (en) 2018-05-14 2022-07-19 Cornell University Collaborative RFID reader using code divisional multiple access (CDMA) and methods for same
WO2020039533A1 (en) * 2018-08-23 2020-02-27 三菱電機株式会社 Biosensor device
GB201917582D0 (en) * 2019-12-02 2020-01-15 Cardiaccs As Estimating ventricular preload
US11925489B1 (en) 2020-11-20 2024-03-12 Stryker Corporation Manifold for filtering medical waste being drawn under vacuum into a medical waste collection system and related methods
US11786647B1 (en) 2022-01-31 2023-10-17 Stryker Corporation Medical waste collection systems, manifolds, and related methods

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5766208A (en) * 1994-08-09 1998-06-16 The Regents Of The University Of California Body monitoring and imaging apparatus and method
US20040249258A1 (en) * 2003-06-04 2004-12-09 Tupin Joe Paul System and method for extracting physiological data using ultra-wideband radar and improved signal processing techniques
US20050124901A1 (en) * 2003-12-05 2005-06-09 Misczynski Dale J. Method and apparatus for electrophysiological and hemodynamic real-time assessment of cardiovascular fitness of a user
WO2007033270A2 (en) * 2005-09-13 2007-03-22 Spectral Energetics, Inc. Device and method for a noninvasive cardiac monitor
WO2008036396A2 (en) * 2006-09-21 2008-03-27 Noninvasive Medical Technologies, Inc. Apparatus and method for non-invasive thoracic radio interrogation
WO2008148040A1 (en) * 2007-05-24 2008-12-04 Lifewave, Inc. System and method for non-invasive instantaneous and continuous measurement of cardiac chamber volume

Family Cites Families (77)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3572316A (en) * 1968-02-23 1971-03-23 Chromalloy American Corp Physiological signal monitoring system
US3646606A (en) * 1969-08-06 1972-02-29 Care Electronics Inc Physiological monitoring system
US3846704A (en) * 1972-12-13 1974-11-05 R Bessette Apparatus for evaluating athletic performance
US3972320A (en) * 1974-08-12 1976-08-03 Gabor Ujhelyi Kalman Patient monitoring system
US4129125A (en) * 1976-12-27 1978-12-12 Camin Research Corp. Patient monitoring system
US4204549A (en) * 1977-12-12 1980-05-27 Rca Corporation Coaxial applicator for microwave hyperthermia
US4270547A (en) * 1978-10-03 1981-06-02 University Patents, Inc. Vital signs monitoring system
US4240445A (en) * 1978-10-23 1980-12-23 University Of Utah Electromagnetic energy coupler/receiver apparatus and method
US4341227A (en) 1979-01-11 1982-07-27 Bsd Corporation System for irradiating living tissue or simulations thereof
US4312358A (en) * 1979-07-23 1982-01-26 Texas Instruments Incorporated Instrument for measuring and computing heart beat, body temperature and other physiological and exercise-related parameters
US4331154A (en) * 1979-10-15 1982-05-25 Tech Engineering & Design Blood pressure and heart rate measuring watch
US4494553A (en) * 1981-04-01 1985-01-22 F. William Carr Vital signs monitor
US4488559A (en) * 1981-06-30 1984-12-18 University Of Utah Apparatus and method for measuring lung water content
US4397313A (en) * 1981-08-03 1983-08-09 Clini-Therm Corporation Multiple microwave applicator system and method for microwave hyperthermia treatment
US4446874A (en) * 1981-12-30 1984-05-08 Clini-Therm Corporation Microwave applicator with discoupled input coupling and frequency tuning functions
WO1983003744A1 (en) * 1982-04-23 1983-11-10 Reinhold Herbert Edward Jr Ambulatory monitoring system with real time analysis and telephone transmission
US4608994A (en) * 1984-12-20 1986-09-02 Matsushita Electric Industrial Co., Ltd. Physiological monitoring system
US4841990A (en) * 1985-06-29 1989-06-27 Tokyo Keiki Co., Ltd. Applicator for use in hyperthermia
US4819860A (en) * 1986-01-09 1989-04-11 Lloyd D. Lillie Wrist-mounted vital functions monitor and emergency locator
US5027829A (en) * 1986-12-15 1991-07-02 Larsen Lawrence E Apparatus for diathermy treatment and control
US4926868A (en) * 1987-04-15 1990-05-22 Larsen Lawrence E Method and apparatus for cardiac hemodynamic monitor
IL86010A (en) * 1987-07-03 1993-02-21 Hi Bridge Kk Apparatus for measuring blood pressure
US4867442A (en) * 1987-10-09 1989-09-19 Matthews H Gerard Physical exercise aid
US4909260A (en) * 1987-12-03 1990-03-20 American Health Products, Inc. Portable belt monitor of physiological functions and sensors therefor
US4966154A (en) * 1988-02-04 1990-10-30 Jonni Cooper Multiple parameter monitoring system for hospital patients
JP2614887B2 (en) * 1988-02-18 1997-05-28 甚一 松田 Local heating device and three-dimensional resonator for local heating
US5178151A (en) * 1988-04-20 1993-01-12 Sackner Marvin A System for non-invasive detection of changes of cardiac volumes and aortic pulses
US4958638A (en) * 1988-06-30 1990-09-25 Georgia Tech Research Corporation Non-contact vital signs monitor
US5423322A (en) * 1988-12-29 1995-06-13 Medical Physics, Inc. Total compliance method and apparatus for noninvasive arterial blood pressure measurement
US5153584A (en) * 1989-03-17 1992-10-06 Cardiac Evaluation Center, Inc. Miniature multilead biotelemetry and patient location system
US5025795A (en) * 1989-06-28 1991-06-25 Kunig Horst E Non-invasive cardiac performance monitoring device and method
US5284144A (en) * 1989-11-22 1994-02-08 The United States Of America As Represented By The Secretary Of The Dept. Of Health & Human Services Apparatus for hyperthermia treatment of cancer
US5027824A (en) * 1989-12-01 1991-07-02 Edmond Dougherty Method and apparatus for detecting, analyzing and recording cardiac rhythm disturbances
US5022402A (en) * 1989-12-04 1991-06-11 Schieberl Daniel L Bladder device for monitoring pulse and respiration rate
US5101836A (en) * 1990-02-27 1992-04-07 The Board Of Trustees Of The Leland Stanford Junior University Flexible low profile microwave array for hyperthermia of superficially located tumors
US4991585A (en) * 1990-03-13 1991-02-12 Mmtc, Inc. Non-invasive respiration and/or heartbeat monitor or the like
US5062432A (en) * 1990-05-09 1991-11-05 Labelle Industries, Inc. Method and apparatus for monitoring personal core temperature
US5186181A (en) * 1990-07-27 1993-02-16 Cafiero Franconi Radio frequency thermotherapy
US5228449A (en) * 1991-01-22 1993-07-20 Athanasios G. Christ System and method for detecting out-of-hospital cardiac emergencies and summoning emergency assistance
US5148002A (en) * 1991-03-14 1992-09-15 Kuo David D Multi-functional garment system
US5335664A (en) * 1991-09-17 1994-08-09 Casio Computer Co., Ltd. Monitor system and biological signal transmitter therefor
FI92139C (en) * 1992-02-28 1994-10-10 Matti Myllymaeki Monitoring device for the health condition, which is attached to the wrist
US5263491A (en) * 1992-05-12 1993-11-23 William Thornton Ambulatory metabolic monitor
US5544651A (en) * 1992-09-08 1996-08-13 Wilk; Peter J. Medical system and associated method for automatic treatment
US5442369A (en) * 1992-12-15 1995-08-15 West Virginia University Toroidal antenna
US5654723A (en) * 1992-12-15 1997-08-05 West Virginia University Contrawound antenna
US5373301A (en) * 1993-01-04 1994-12-13 Checkpoint Systems, Inc. Transmit and receive antenna having angled crossover elements
US5404877A (en) * 1993-06-04 1995-04-11 Telectronics Pacing Systems, Inc. Leadless implantable sensor assembly and a cardiac emergency warning alarm
US5748002A (en) * 1996-01-26 1998-05-05 Phase Dynamics Inc. RF probe for montoring composition of substances
US5724025A (en) * 1993-10-21 1998-03-03 Tavori; Itzchak Portable vital signs monitor
US5416468A (en) * 1993-10-29 1995-05-16 Motorola, Inc. Two-tiered system and method for remote monitoring
HRP931478A2 (en) * 1993-12-06 1995-12-31 Ferek Petri Bo Idar An apparatus for cardial electrotherapy containing transmission line uring cardial contractions
US5829437A (en) * 1994-07-01 1998-11-03 Interstitial, Inc. Microwave method and system to detect and locate cancers in heterogenous tissues
US5704355A (en) * 1994-07-01 1998-01-06 Bridges; Jack E. Non-invasive system for breast cancer detection
US5778882A (en) * 1995-02-24 1998-07-14 Brigham And Women's Hospital Health monitoring system
US5683382A (en) * 1995-05-15 1997-11-04 Arrow International Investment Corp. Microwave antenna catheter
US5769879A (en) * 1995-06-07 1998-06-23 Medical Contouring Corporation Microwave applicator and method of operation
US5633875A (en) * 1995-06-07 1997-05-27 General Electric Company Protocol and mechanism for centralized asset tracking communications
US5734353A (en) * 1995-08-14 1998-03-31 Vortekx P.C. Contrawound toroidal helical antenna
US5841288A (en) * 1996-02-12 1998-11-24 Microwave Imaging System Technologies, Inc. Two-dimensional microwave imaging apparatus and methods
US5771001A (en) * 1996-11-18 1998-06-23 Cobb; Marlon J. Personal alarm system
US5861019A (en) * 1997-07-25 1999-01-19 Medtronic Inc. Implantable medical device microstrip telemetry antenna
US6064903A (en) * 1997-12-29 2000-05-16 Spectra Research, Inc. Electromagnetic detection of an embedded dielectric region within an ambient dielectric region
EP1067860A1 (en) * 1998-03-30 2001-01-17 Rush-Presbyterian-St.Luke's Medical Centre Method and apparatus for detection of air cavities in a body
US6368286B1 (en) * 1998-06-04 2002-04-09 Eric D. Whitman Pneumothorax detector
US6233479B1 (en) * 1998-09-15 2001-05-15 The Regents Of The University Of California Microwave hematoma detector
AUPR343401A0 (en) * 2001-02-28 2001-03-29 Nguyen, Hung Modelling and design for early warning systems using physiological responses
FI118403B (en) * 2001-06-01 2007-10-31 Pulse Finland Oy Dielectric antenna
US6522196B1 (en) * 2001-08-20 2003-02-18 Harris Corporation System and method for improving input return loss in RF amplifiers
US6829501B2 (en) * 2001-12-20 2004-12-07 Ge Medical Systems Information Technologies, Inc. Patient monitor and method with non-invasive cardiac output monitoring
US6648828B2 (en) * 2002-03-01 2003-11-18 Ge Medical Systems Information Technologies, Inc. Continuous, non-invasive technique for measuring blood pressure using impedance plethysmography
US7450746B2 (en) * 2002-06-07 2008-11-11 Verathon Inc. System and method for cardiac imaging
JP3880547B2 (en) * 2003-05-30 2007-02-14 独立行政法人科学技術振興機構 Electrode for measuring cardiac output and cardiac output measuring device
JP4742644B2 (en) * 2004-03-31 2011-08-10 日本光電工業株式会社 Blood volume measuring method, measuring apparatus and biological signal monitoring apparatus
US20090082679A1 (en) * 2004-06-21 2009-03-26 Impedance Cardiology Systems, Inc. Cardiac monitoring system
WO2008036404A2 (en) * 2006-09-21 2008-03-27 Noninvasive Medical Technologies, Inc. Antenna for thoracic radio interrogation
EP2224849B1 (en) * 2007-12-21 2014-07-23 St. Jude Medical AB Device for monitoring acute decompensated heart failure

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5766208A (en) * 1994-08-09 1998-06-16 The Regents Of The University Of California Body monitoring and imaging apparatus and method
US20040249258A1 (en) * 2003-06-04 2004-12-09 Tupin Joe Paul System and method for extracting physiological data using ultra-wideband radar and improved signal processing techniques
US20050124901A1 (en) * 2003-12-05 2005-06-09 Misczynski Dale J. Method and apparatus for electrophysiological and hemodynamic real-time assessment of cardiovascular fitness of a user
WO2007033270A2 (en) * 2005-09-13 2007-03-22 Spectral Energetics, Inc. Device and method for a noninvasive cardiac monitor
WO2008036396A2 (en) * 2006-09-21 2008-03-27 Noninvasive Medical Technologies, Inc. Apparatus and method for non-invasive thoracic radio interrogation
WO2008148040A1 (en) * 2007-05-24 2008-12-04 Lifewave, Inc. System and method for non-invasive instantaneous and continuous measurement of cardiac chamber volume

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ENRICO M STADERINI: "UWB Radars in Medicine", IEEE AEROSPACE AND ELECTRONIC SYSTEMS MAGAZINE, IEEE SERVICE CENTER, PISCATAWAY, NJ, US, vol. 17, no. 1, 1 January 2002 (2002-01-01), pages 13 - 18, XP011092427, ISSN: 0885-8985 *
NEWMAN DAVID G ET AL: "The non-invasive assessment of stroke volume and cardiac output by impedance cardiography: A review", AVIATION, SPACE AND ENVIRONMENTAL MEDICINE, AEROSPACE MEDICAL ASSOCIATION, ALEXANDRIA, VA, US, vol. 70, no. 8, 1 August 1999 (1999-08-01), pages 780 - 789, XP008093991, ISSN: 0095-6562 *
STADERINI E M ET AL: "Optimization criteria in the design of medical UWB radars in compliance with the regulatory masks", BIOMEDICAL CIRCUITS AND SYSTEMS CONFERENCE, 2007. BIOCAS 2007. IEEE, IEEE, PISCATAWAY, NJ, USA, 27 November 2007 (2007-11-27), pages 53 - 58, XP031237030, ISBN: 978-1-4244-1524-3 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013530741A (en) * 2010-05-28 2013-08-01 ザ ジェネラル ホスピタル コーポレーション Apparatus, system, method and computer-accessible medium for analyzing information on cardiovascular disease and cardiac function
US11457889B2 (en) 2021-01-25 2022-10-04 Lazaro Eduardo Hernandez System and method of non-invasive continuous echocardiographic monitoring

Also Published As

Publication number Publication date
CN101983030B (en) 2013-09-18
JP2011509123A (en) 2011-03-24
CN101983030A (en) 2011-03-02
US8721559B2 (en) 2014-05-13
MX2010007414A (en) 2011-02-23
US20100022900A1 (en) 2010-01-28
JP5406211B2 (en) 2014-02-05
EP2237719A1 (en) 2010-10-13

Similar Documents

Publication Publication Date Title
US8721559B2 (en) Non-invasive method and device for measuring cardiac output
US11806121B2 (en) Methods of identifying and locating tissue abnormalities in a biological tissue
US20230148891A1 (en) System and method for non-invasive instantaneous and continuous measurement of heart rate, stroke volume and ejection fraction
EP2506917B1 (en) Microwave monitoring of heart function
JP6081355B2 (en) Embedded radio frequency sensor
Iskander et al. Electromagnetic techniques for medical diagnosis: A review
JP5432254B2 (en) Automatic transmission output adjustment for Doppler radar
JP2008515548A (en) Continuous field tomography
JP2015511843A (en) Monitoring and diagnostic systems and methods
WO2008027509A2 (en) Remote-sensing method and device
CN110573069B (en) Biological information estimation device
Pedersen et al. An investigation of the use of microwave radiation for pulmonary diagnostics
Brovoll et al. Optimal frequency range for medical radar measurements of human heartbeats using body-contact radar
Libove et al. Wearable brain imager/BMI technology for structural, vascular and functional extraction
Lin Microwave noninvasive sensing of physiological signatures
JPH02109538A (en) Method and apparatus for detecting position of movable boundary between biological tissues

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200980107669.8

Country of ref document: CN

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 09700309

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2010541585

Country of ref document: JP

Ref document number: 2009700309

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: MX/A/2010/007414

Country of ref document: MX

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 1607/MUMNP/2010

Country of ref document: IN