The present invention relates to a dynamic anthropomorphic cardiovascular phantom, model or system and, particularly, a dynamic anthropomorphic cardiovascular phantom for simulating injectable fluid propagation.
The following information is provided to assist the reader to understand the invention disclosed below and the environment in which it will typically be used. The terms used herein are not intended to be limited to any particular narrow interpretation unless clearly stated otherwise in this document. References set forth herein may facilitate understanding of the present invention or the background of the present invention. The disclosure of all references cited herein is incorporated by reference.
Computed tomography (CT) is a medical imaging technique in which digital geometry processing is used to generate a three-dimensional imaging of a region of interest (ROI) within a body from a series of two-dimensional X-ray images taken around a single axis of rotation. A contrast medium injected into the patients venous system is typically used to enhance the contrast of various structures and/or fluids within the region of interest. The development of subsecond rotation and multi-slice CT of MSCT (for example, up to 64-slice) has provided both high resolution and high speed. Moreover, images with even higher temporal resolution can be generated using electrocardiogram (ECG) gating to correlate CT data with corresponding phases of cardiac contraction.
As MSCT enables the acquisition of sub-millimeter slices and near isotropic voxels, anatomical territories previously unacquirable are now almost routinely imaged. In particular, dynamic structures such as the coronary vasculature, are able to be imaged as a result of the MSCT scanner's ability to acquire a volumetric data set in 5-20 seconds, well within the breath-hold interval of most patients. Furthermore, the latest generation of MSCT also enables cardiac CT even in the presence of high heart-rates and arrhythmia. Cardiac CT Angiography (CCTA) is therefore a demanding imaging regime in which the clinician must perfect her technique to maximize image quality.
Intrinsic to the quality of a CCTA exam is the proper dosing, delivery, and timing of the iodinated contrast bolus required for image contrast enhancement. Because CCTA is only concerned with arterial imagery during the first-pass of the contrast bolus, the timing of the scanner acquisition relative to the peak contrast enhancement in the cardiac anatomy is important. There are many published works demonstrating the benefits of personalizing and optimizing the delivery of contrast material to the individual and the image procedure. In the case of CCTA with MSCT scanners, the optimization and timing of the contrast bolus is crucial.
The widespread clinical adoption of multi-slice CT has lead to challenges in adapting imaging and contrast delivery techniques which were developed in connection with single-slice, helical CT scanners. Development of optimal contrast injection techniques to generate ideal contrast enhancement has been the subject of numerous CT studies. The outcomes of such studies, however, are often not comparable and interpretation controversy arises because of variations in injection technique, parameters, and contrast medium properties. A consensus for optimal injection parameters and choice of contrast media for intravenous contrast enhancement in CT-scanning has yet to be reached and further investigation is ongoing and mandatory as scanning technology advances. To best compare different injection protocols, hemodynamic values like blood pressure, blood volume and cardiac output should ideally be held constant, as these factors significantly influence contrast enhancement. With non-uniform distribution of these factors, comparison of the contrast application is difficult to assess. Furthermore, measurements of time-enhancement curves at defined anatomic sites are required to ensure an exact analysis of contrast enhancement. To obtain a time-enhancement curve, serial CT scans at one anatomic level are necessary. However, such serial scans are not feasible in a patient study as a result of the increased radiation burden of the long acquisition time.
In an attempt to avoid the increased radiation burden and expense associated with CT patient studies, use of a flow phantom to simulate or emulate convective transport properties of the cardiovascular system has been studied. Repeated injections of contrast demonstrated some utility of the model in replicating the contrast enhancement pattern of the abdominal aorta in MSCT. Although utility was demonstrated for such a flow phantom, a number of significant limitations hinder the use of the flow phantom in developing improved CT injection protocols for contrast media.
It is, therefore, desirable to develop improved cardiovascular flow phantoms, models or systems for use in studies of the propagation of injectable fluids, including, for example, contrast media and/or other drugs.
In one aspect, a cardiovascular flow system includes: a cardiovascular model system, a pump system in fluid connection with the cardiovascular model system, and an ECG simulator in communicative connection with the pump system. The ECG simulator system is adapted to create and transmit a simulated ECG signal. The ECG simulator system uses a signal received from the pump system to adjust the simulated ECG signal transmitted from the ECG simulator system. The cardiovascular flow system further includes an injection port adapted to be placed in fluid connection with an injector to inject at least one fluid into the system.
The cardiovascular model system of the cardiovascular flow system can simulate mammalian physiologic hemodynamic parameters. For example, the hemodynamic parameters can include, but are not limited to, systolic pressure, diastolic pressure and convective transport of an injected fluid from the injection port to at least one region of interest within the cardiovascular model system.
The cardiovascular model system can additionally or alternatively simulate mammalian anatomy for at least one region of interest. The at least one region of interest can, for example, include the coronary artery tree and the thoracic aorta. The three-dimensional conformation of the region of interest can be simulated. Furthermore, movement within the at least one region of interest (for example, resulting from the heart's pulsatile pumping action) can be simulated. In that regard, motion within the region of interest can occur as a result of pulsatile flow of fluid therethrough occurs.
The ECG system can, for example, include a processing system to adjust the ECG signal to synchronize to the ECG signal to the signal received from the pump system. The ECG system can include a controller adapter to adjust timing between the ECG signal and signal received from the pump system. In several embodiments, the pump system includes a sensor (for example, a microswitch) in operative connection therewith that is adapted to provide the signal from the pump system. In a number of embodiments, the ECG system interfaces with an imaging system or scanner to transmit the simulated ECG signal to the scanner.
The cardiovascular model system can further include a lung sink (simulating flow through the lungs) and a body sink (simulating flow through several components of the body/system circulation) in fluid connection within the cardiovascular model system.
One or more interconnections within the cardiovascular model system can be formed to simulate at least one of mammalian physiology and anatomy. In several embodiments, the interconnections are selected to simulate physiologic pressure and/or convective transport.
The cardiovascular flow system can further include a kidney simulation system in fluid connection with the cardiovascular model system. The kidney system can, for example, include a first pump system to pump fresh fluid into the cardiovascular model system and a second pump system to pump fluid out of the cardiovascular system. The kidney simulation system can, for example, be controlled to control at least one of fluid volume and pressure within the cardiovascular model system.
The cardiovascular flow system can further include a pressure relief system to prevent pressure within the system from exceeding a predetermined pressure. The pressure relief system can, for example, be adapted to maintain pressure within the cardiovascular model system within limits of physiologic pressure.
In several embodiments, the cardiovascular model system includes a cardiac vasculature model including a pulmonary artery crossing the cardiac anatomy.
The cardiovascular flow system can further include the injector.
In another aspect, a method of simulating flow of an injection fluid through a mammalian cardiovascular system includes: providing a cardiovascular model, the cardiovascular model comprising an injection port adapted to be placed in fluid connection with an injector to inject the injection fluid into the cardiovascular model; pumping fluid through the cardiovascular model; configuring cardiac output, heart rate, stroke volume in the cardiovascular model; providing a simulated ECG signal, wherein the ECG signal is synchronized with pumping pulse or pumping pressure; injecting the injection fluid into the cardiovascular model via the injection port; and detecting the injection fluid in at least one region of interest within the cardiovascular model.
The synchronization can, for example, be based on an aortic pressure.
The method can further include scanning the at least one region of interest of the cardiovascular model with an imaging system to produce an image. The injection fluid can, for example, include a contrast medium (or contrast enhancement medium).
The cardiovascular model can include a coronary vasculature model, a body sink and a lung sink in fluid connection. The coronary vasculature model, the body sink and the lung sink can, for example, be in fluid connection via tubing adapted to simulate (or to enable simulation of) physiological pressures and/or convective transport.
In another aspect, a cardiovascular flow system includes: a cardiovascular model system, a pump system in fluid connection with the cardiovascular model system, an injection port adapted to be placed in fluid connection with an injector to inject at least one fluid into the system; and a kidney system or kidney simulation system in fluid connection with the cardiovascular model system. The kidney system includes a first pump system to pump fresh fluid into the cardiovascular model system and a second pump system to pump fluid out of the cardiovascular system. The first pump system and the second pump system can, for example, be controlled to prevent fluid volume increase within the cardiovascular model system upon injection of fluid into the cardiovascular model system via the injection port. The kidney system can be controlled to control pressure within the cardiovascular model system.
In a further aspect, a cardiovascular flow system includes: a cardiovascular model system, a pump system in fluid connection with the cardiovascular model system, an injection port adapted to be placed in fluid connection with an injector to inject at least one fluid into the system; and a pressure relief system to prevent pressure within the cardiovascular model system from exceeding a predetermined pressure. The pressure relief system can, for example, be adapted to maintain pressure within the cardiovascular model system within limits of physiologic pressure.
In still a further aspect, a cardiovascular flow system includes: a cardiovascular model system for simulating anatomy and physiology, a pump system in fluid connection with the cardiovascular model system, and an ECG simulator in communicative connection with the pump system. The ECG simulator system is adapted to create and transmit a simulated ECG signal. The ECG simulator system uses a signal received from the pump system to adjust the simulated ECG signal transmitted from the ECG simulator system. The cardiovascular flow system further includes an injection port adapted to be placed in fluid connection with an injector to inject at least one fluid into the system.
In several embodiments, the flow models, phantoms or systems described herein simulate or mimic the delicate anatomy of the coronary vasculature (including, for example, the coronary artery tree and the thoracic aorta) as well as the pulsatile nature of the heart and the motion resulting therefrom. The crossing of the coronary arteries by the pulmonary artery is also simulated or mimicked.
In several embodiments, an ECG simulation system interfaces to the ECG triggering or gating system of an imaging system or scanner. The flow models, phantoms or systems of the present invention also simulate or mimic physiologic hemodynamic parameters including, but not limited to, blood pressure (diastolic blood pressure, systolic blood pressure) and convective transport of an injection fluid to one or more regions of interest within the flow models, phantoms or systems. Users can also configure the system (that is, adjust one or more system variables) to achieve desired physiologic conditions.
The systems of the present invention improve characterization of the dynamics related to early-phase contrast enhancement in the cardiac and/or other anatomy to assist, for example, in the development and validation of novel, optimal contrast delivery control approaches. The systems of the present invention provide, among other things, configurable cardiac output, heart rate, stroke volume, and blood volume, accurate transport delays of an injected bolus, accurate dilution of the contrast media through a cardiopulmonary circuit, the ability to interface an ECG signal to the scanner, accurate motion of a region of interest such as the coronary tree, pressures, flow rates and other conditions within the physiologic realm, which are responsive to changes in the driving function, and pulsatile blood flow and heart rates.
The present invention, along with the attributes and attendant advantages thereof, will best be appreciated and understood in view of the following detailed description taken in conjunction with the accompanying drawings.
As used herein and in the appended claims, the singular forms “a,” “an”, and “the” include plural references unless the content clearly dictates otherwise. Thus, for example, reference to “a pressure transducer” includes a plurality of such pressure transducers and equivalents thereof known to those skilled in the art, and so forth, and reference to “the pressure transducer” is a reference to one or more such pressure transducers and equivalents thereof known to those skilled in the art, and so forth.
Although the anthropomorphic cardiovascular flow systems, cardiovascular phantom/model systems or flow phantoms described herein are discussed in connection with the use of contrast media in CT studies, and particularly cardiac CT angiography studies, one skilled in the art appreciates that the cardiovascular flow systems can be used to study the propagation of virtually any type of fluid that can be detected using one or more detection systems to virtually any region of interest within the cardiovascular flow system.
Coronary tree model 100 includes a modification that allows a model of the pulmonary artery 160 to cross the cardiac anatomy. The pulmonary artery cross enables acquisition of single-level, axial scans mimicking the contrast enhancement pattern at the pulmonary trunk in the human vascular system. System 10 incorporates the novel aspects of the pulmonary artery and aortas to generate time enhancement curves of the contrast propagation through the cardiopulmonary circuit similar to that in humans.
A ventricular pump system 200 is placed in fluid communication with coronary tree model 100 and is used too simulate the heart in system 10. A representative example of a ventricular pump suitable for use in system 10 (and used in the studies set forth herein) is the model BS4 553305 pulsatile blood pump available from Harvard Apparatus Inc. of Holliston Mass., USA (sometimes referred to as the Harvard pump). Pump system 200 can, for example, be configured to deliver stroke volumes of 15 to 100 ml, heart rates from 10 to 100 beats per minute (BPM), and systolic/diastolic ratios from 35% to 50% full cycle. Pump system 100 is designed to safely administer blood or a blood emulant, however distilled water may also be used in place of the blood emulant (that is, as a first or carrier fluid simulating blood and different from a second or injection fluid to be propagated through system 10).
In several embodiments, a trigger system is connected to pump system 200 to generated a trigger signal at the end of full displacement. The trigger system includes a sensor such as a microswitch 300 (see
Coronary tree model 100, body sink 600, and a lung sink 620 are connected via tubing, which is sized (for example, in diameter and length) to simulate corresponding physiological characteristics. In several studied embodiments, connecting tubing (fabricated from TYGON® plastic tubing) with approximately 1.0 cm to 2.5 cm (O.D.) was used. The lengths of a number of the connecting tubing segments are based on the constraint of the CT scanner gantry. As, for example, illustrated in
The tubing dimensions also provide physiologic pressures in system 10. For example, Tables 1 and 2 set forth the physical dimensions of the various tubing segments associated with system as illustrated in
System 10 is calibrated so that the central fluid volume provides a simulation of the mammalian (for example, human) cardiovascular circuit. The fluid volume may, for example, be set between 3 and 8 L. A number of components of system 10 are described below in connection with filling and draining procedures for system 10.
During filling of system 10, a reservoir tank 630 is filled with 3 liters (L) of water (more than 3 L may not allow sufficient pressure to be generated in system 10). Reservoir tank 620 is connected by a tube to a fill port 640. Fill port 640 and a vent 604 on body sink 600 are opened. Tubing connected to body vent 604 is placed into an empty container. A vent 624 on lung sink 620 is closed. A drain port 650 is set to a normal setting. At this point, reservoir water from reservoir tank 630 is pumped into system 10 until empty (using, for example, a pump handle 632 of reservoir tank 630). Pressure within reservoir tank 630 is then released (for example, by unscrewing a center piston in connection with handle 632 as known in the art for pressure tanks). Fill port 640 is then closed.
Reservoir tank 630 is then filled with 3 liters of water. Reservoir tank 630 is reconnected to fill port 640, which is then opened. Reservoir tank 630 is then pumped until body sink 600 is full. Fill port 640 and body vent port 604 are then closed. Tubing from lung vent port 624 is then placed into an empty container and lung vent port 624 is opened. Fill port 640 is opened and water is pumped into system 10 until lung sink 620 is full. At this point, fill port 640 is closed and lung vent port 624 is closed.
Pump system 200 is then, for example, run at 30 beats per minute (BPM), with a stroke volume of 90 ml, and a Systolic/Diastolic % of 50/50. Fill port 640 is opened. From this point forward, pressure should be monitored via one or more pressure transducer out connectors (see, for example,
At this point, water can be added to reservoir tank 630 if, for example, an excess container has more than 400 ml of water therein. The additional water (ml) to be added can be set to equal the volume of water in the excess container (ml) minus 400 ml. Fill port 640 is opened, and lung vent 624 is opened to evacuate air. Reservoir tank 630 is pumped until system 10 has been evacuated of air. After air evacuation, fill port 640 is closed. Lung vent 624 is also closed. The total water volume in system 10 is now 5.5 liters. This volume can be adjusted as desired.
Tubing from drain port 650 is then placed into the excess container (not shown). Reservoir tank 630 is emptied. Empty reservoir tank 630 is connected to body sink vent 604. Drain port 640 is now set to “drain”. Body sink vent 620 is opened and reservoir tank 630 is pumped to effect a desired system volume and to allow for additional injected volume (via an injection port 642). A graduated scale can, for example, be provided on a side of body sink 600. Drain port 640 is then once again set to normal and body sink vent 604 is closed.
To begin a drain procedure for system 10, drain tubing from drain/fill port 640 is placed into the empty excess container, which is positioned below the level of system 10. Empty reservoir tank 630 is connected to body sink vent 604, which is then opened. A drain clamp (not shown) is then opened and drain/fill port 640 to set to drain. Reservoir tank 30 is pumped to evacuate water until body sink 600 is empty. Body sink vent 604 is then closed and empty reservoir tank 630 is connected to lung sink vent 624, which is then opened. Reservoir tank 630 is pumped to evacuate the remaining water. Subsequently, pump system 200 is disconnected via connection fittings as known in the art and any remaining water is removed.
As, for example, illustrated in
In several embodiments, a dilution, filter or kidney system 800 is used to simulate operation of the kidneys. For example, in-line pump system can be used to filter contrast or other fluid for which propagation is being studied from system 10 during and/or between experiments to simulate function of a kidney. Kidney system 800 performs the function of removing “waste products” from system 10 and regulating the fluid level therein. As described herein, a typical use of system 10 includes injecting a contrast dye into the system 10 via injection port 642. As the contrast dye is injected into system 10, the additional injected volume of the contrast fluid is added to the existing fluid in system 10. Kidney system 800 can, for example, be configured and/or controlled to remove a volume of fluid equal to the additional injected volume of fluid. In addition, as the contrast is injected into system 10, the contrast mixes with the original fluid in system 10. At some point, the mixture of the fluids reduces the usability of system 10. Kidney system 800 can be used to replenish system 10 with fresh fluid (for example, water) while also removing “waste” fluids from system 10 (thereby reducing the concentration of injected fluid within system 10). A block diagram of kidney system 800 is set forth in
In the illustrated embodiment, kidney system 800 includes a first or fresh reservoir 810 containing fresh fluid (for example, a fluid having no concentration of the propagation fluid being studied—such as fresh water) to replenish system 10 and a second or waste reservoir 820 to collect waste fluid (that is, previously circulated fluid which can include some concentration of the propagation fluid) from system 10. Each of fresh reservoir 810 and waste reservoir 820 includes a fluid level sensor 812 and 822, respectively, to measure fluid levels therein. Fluid level sensor 812 is, for example, used to detect when fresh reservoir 910 is depleted and needs to be refilled. Fluid level sensor 822 is, for example, used to detect when waste reservoir 820 is full or nearly full and needs to be drained. A first fluid pumping system 830 is in fluid connection with fresh reservoir 810 to deliver fresh fluid into the system 10. A second fluid pumping system 840 is in fluid connection with waster reservoir 820 to remove excess fluid injected into system 10 and also to remove the waste fluid from the system 10. At least one motor drive 850 is in operative connection with pumping systems 830 and 840 to drive fluid pumping. Various connectors such as connecting luers and tubing as known in the art can be used to place the elements of kidney system 800 (and system 10 generally) in fluid connection to transport fluids therebetween. One or more check valves 860 can be used to ensure a desired direction of fluid flow. Kidney system 800 can be in communicative connection with a data acquisition/control system 900 of system 10, which can, for example, be used (with respect to kidney system 800) to read data from level sensors 812 and 822, to control pump motor speed and to effect on/off functions.
Pumping mechanisms suitable for use as pumping system 830 and pumping system 840 are described, for example, in U.S. Pat. Nos. 5,916,197 and 6,197,000, the disclosure of which are incorporated herein by reference. In one embodiment, pumping systems 830 and 840 were operatively connected to a common motor drive so that the amount of fluid pumped into system 10 by first pumping system 830 was always approximately equal or equal to the amount of fluid pumped out of system 10 by second pumping system 850. The flow rate of the pumping systems can be adjusted to control the filtering/dilution process. To be operative to adjust fluid volume within system 10 (for example, to remove a volume of fluid from system 10 equal to a volume injected via injection port 642 to maintain a generally constant volume of fluid in system 10) and/or to control pressure within system 10, first pumping system 830 and second pumping system 840 can be controlled separately (that is, to effect addition of or removal of a determined volume fluid from system 10).
In several embodiments, a manually adjustable automatic pressure relief system 680 was incorporated within system 10. In a number of studies, a commercially available model R-4104-10 pressure regulator available from Airtrol Components, Inc. of New Berlin, Wis. USA was used. See U.S. Pat. No. 4,315,520. Pressure relief and/or control can alternatively be accomplished by other means. For example, one or more electromechanically controlled valves can be used. Indeed, additional benefits can be realized with the use of an electromechanical valve. For example, a quicker response time to overpressure conditions can be achieved, as well as automatic control over the maximum pressure limit and pressure variability. Moreover, use of an electromechanical control valve in conjunction with data acquisition/control system 900 can provide the ability to control, monitor, store and analyze performance data.
In the embodiment illustrated in
Pressure relief system 680 alleviates several potential issues with the use of system 10. For example, as system 10 is being prepared for initial use (that is, during filling and priming) an overpressure condition can be encountered. Pressure relief system 680 assists in preventing overpressures that can potentially cause failures in fittings and/or damage to compliant coronary tree model 100.
Furthermore, without pressure relief system 680, contrast injected into system 10 during a power injection can increase the pressure to undesirable levels. Once again, fittings/connectors could fail and/or coronary tree model 100 could be damaged. Moreover, without the use of pressure relief system 680 during an injection, control of the system pressure to be within a physiological blood pressure range (for example, set to 120/80 mmHg) is difficult in that system pressure can increase significantly during and after the injection.
Several experiments were conducted to measure the effectiveness of pressure relief system 680.
Cardiac output can be adjusted from, for example, 2.5 L/min to 8 L/min based on the setting provided by (Harvard) pump system 200. Measurements of flow rate are made within system 10 by one or more flow rate sensors 690 such as in-line, ultrasonic transducers. Suitable flow rate sensors are, for example, available from Transonic Inc., if Ithaca, N.Y. USA. Once again, intravascular pressure is measured with pressure transducers such as pressure transducers 662 and 664. An example of a suitable pressure transducer is the model PX35D pressure sensor available from Omega Inc, of Stamford, Conn. USA.
Analog data from transducers or sensors, along with the ECG triggering pulse can, for example, be acquired with data acquisition (DAQ) hardware/software system which can from a part of data acquisition/control system 900. Suitable data acquisition systems are for example, available from National Instruments of Austin, Tex. USA (NIDAQ). Data acquisition/control system 900 can, for example, also include one or more processors (for example, one or more microprocessors) and associated memory, display etc. as known in the art. Data acquisition/control system 900 can, for example, include a desktop or laptop computer in communicative connection with data acquisition hardware.
The fluid level or volume, forced pressure and pump settings of system 10 are readily configurable to keep the fluid (blood emulant) pressure within typical physiologic limits, for example, 70 mmHg diastolic pressure to 120 mmHg systolic (see, for example,
The design of the phantom system was adjusted to mimic the convective transport delay of contrast material typical in humans for a particular region of interest. The dimensions and volumes of components of system 10 (for example, as set forth in Table 1) were determined by empirical study and mathematical modeling (using known engineering principles) to achieve desired transport delays of contrast in humans. It is typical, for example, for contrast to arrive in the pulmonary artery 8-13 seconds after the start of contrast injection. Further, contrast material typically arrives in the coronary arteries 15-22 seconds after the start of contrast injection. In several embodiments, the delay from the injection site to the simulated pulmonary artery in system 10 ranges from approximately 5-12 seconds, as is typical in humans. Furthermore, the delay of contrast from the simulated pulmonary artery to the simulated coronary arteries and thoracic artery ranges from approximately 6-18 seconds. It is important, when evaluating contrast injection protocols, for the arrival of contrast in such regions of interest to be representative of the delay encountered in the human cardio-pulmonary circuit. Scanning should be commensurate with the arrival of the contrast bolus in a desired region of interest. In light of the studies of system 10 (which has not been fully optimized), methods and protocols validated using system 10 are applicable for use on humans.
As described above, system 10 can be used to study the propagation of virtually any type of fluid that can be detected using one or more detection systems. In that regard, any fluid for which concentration (or a variable related to concentration) can be detected using a sensor can be studied in the present invention. The sensor can be embodied in a scanner system such as a CT scanner, an MRI scanner, an ultrasound scanner etc. as known in the medical imaging arts, which are suitable for measuring contrast enhancement (which is related to contrast concentration). Other types of sensors (for example, represented by contrast sensors and light imaging system 1100 in
ECG Enhancement:
As described above, system 10 includes ECG simulator system 100 to provide a synchronized ECG scanner-gating signal to scanner 500. Implementation of the ECG signal into system 10 was achieved in several embodiments using a signal modifier system in communication with simulator 400 (such as a BIO-TEK ECG Plus simulator) and trigger system 300 in operative connection with pump system 200. A typical ECG simulator transmits a non-synchronized, continuous ECG output waveform. In the case of system 10, electromechanical modifications provide an ECG gating implementation.
A first modification requires generating a synchronizing pulse signal from pump 200 and transmitting the signal to ECG simulator/system 400. Referring, for example, to
A second modification requires the addition of electronics to ECG simulator 400. Referring to
Referring to
A high level output signal of ECG simulator 400 can be obtained at + and − connectors 450 and 460, respectively, on the front panel of ECG simulator 400. The output signal is at a level sufficient to enable data acquisition. This signal can be used with a data acquisition system to monitor, store and analyze the ECG signal.
To demonstrate the ability of system 10 to reflect the effects on contrast bolus morphology and transport, several injection and imaging experiments were designed with a Siemens 64-Slice SOMATOM SENSATION 64 CT scanner. Bae et al demonstrated the linearity (additive) time invariant properties of contrast pharmacokinetics with respect to injection duration. See Bae, K. T., “Peak contrast enhancement in CT and MR angiography: when does it occur and why? Pharmacokinetic study in a porcine model,” Radiology, 2003. 227(3): p. 809-16. It is expected for a linear system, that the contrast enhancement pattern resulting from an injection of a long duration (measured at one-level in the vasculature) should be a linear superposition of injections with shorter durations.
The experiments of Bae et al were replicated by performing subsequent bolus injections into the phantom using a power injector (the STELLANT D injector available from MEDRAD, Inc. of Pittsburgh, Pa. USA) and measuring the resulting enhancement pattern at a fixed level in a region of interest of system/model 10a (that is, the pulmonary trunk). Bolus injections of iopromide 370 mgVml (ULTRAVIST 370 available from Bayer Schering Pharmaceutical of Berlin, Germany) lasting 5, 10, 15 and 20 seconds were repeated at injection flow rates of 4, 5 and 6 ml/s. The phantom/system parameters used in this experiment are set forth in Table 3. The scans were acquired every second from 2 seconds post injection to 40 seconds post injection.
Experimental Results
To test the linear predictive value of a test bolus, an imaging experiment was executed in which 20 ml of iopromide 370 mgI/ml was injected at 4 ml/s and followed by a saline flush of 50 ml. A time enhanced curve (TEC) was measured in the ascending aorta. Another injection of 80 ml at 4 ml/s was performed and the resulting TEC was measured.
Linearity Results
The results of the linearity injection for contrast at 6 ml/s (5, 10, 15, and 20 second) are illustrated in
Identification Experiments
The measurement from a 20 ml test bolus was used to construct a transfer function for the system. The resulting enhancement prediction (using 20 ml bolus as the input) is plotted against the actual signal in
Having an estimated transfer function, a predicted enhancement pattern was generated for an 80 ml injection. The resulting prediction is set forth in
Visual Example
An example of contrast enhancement visualization is presented in
The foregoing description and accompanying drawings set forth the preferred embodiments of the invention at the present time. Various modifications, additions and alternative designs will, of course, become apparent to those skilled in the art in light of the foregoing teachings without departing from the scope of the invention. The scope of the invention is indicated by the following claims rather than by the foregoing description. All changes and variations that fall within the meaning and range of equivalency of the claims are to be embraced within their scope.
This application claims benefit of U.S. Provisional Patent Application Ser. No. 61/033,761, filed Mar. 4, 2008, the disclosure of which is incorporated herein by reference.
Number | Name | Date | Kind |
---|---|---|---|
3701345 | Heilman et al. | Oct 1972 | A |
3812843 | Wootten et al. | May 1974 | A |
3888239 | Rubinstein | Jun 1975 | A |
3898983 | Elam | Aug 1975 | A |
3941126 | Dietrich et al. | Mar 1976 | A |
4090502 | Tajika | May 1978 | A |
4135247 | Gordon et al. | Jan 1979 | A |
4151845 | Clemens | May 1979 | A |
4187057 | Xanthopoulos | Feb 1980 | A |
4191183 | Mendelson | Mar 1980 | A |
4207871 | Jenkins | Jun 1980 | A |
4280494 | Cosgrove, Jr. et al. | Jul 1981 | A |
4315520 | Atkinson et al. | Feb 1982 | A |
4319568 | Tregoning | Mar 1982 | A |
4340153 | Spivey | Jul 1982 | A |
4392849 | Petre et al. | Jul 1983 | A |
4396385 | Kelly et al. | Aug 1983 | A |
4402310 | Kimura | Sep 1983 | A |
4409966 | Lambrecht et al. | Oct 1983 | A |
4444198 | Petre | Apr 1984 | A |
4447230 | Gula et al. | May 1984 | A |
4477923 | Baumann et al. | Oct 1984 | A |
4479761 | Bilstad et al. | Oct 1984 | A |
4512764 | Wunsch | Apr 1985 | A |
4544949 | Kurihara | Oct 1985 | A |
4551133 | Zegers de Beyl et al. | Nov 1985 | A |
4559036 | Wunsch | Dec 1985 | A |
4563175 | LaFond | Jan 1986 | A |
4585009 | Barker et al. | Apr 1986 | A |
4611340 | Okazaki | Sep 1986 | A |
4633307 | Honda | Dec 1986 | A |
4634426 | Kamen | Jan 1987 | A |
4636144 | Abe et al. | Jan 1987 | A |
4655197 | Atkinson | Apr 1987 | A |
4672651 | Horiba et al. | Jun 1987 | A |
4710166 | Thompson et al. | Dec 1987 | A |
4723261 | Janssen et al. | Feb 1988 | A |
4750643 | Wortrich | Jun 1988 | A |
4754786 | Roberts | Jul 1988 | A |
4783273 | Knutsson et al. | Nov 1988 | A |
4795429 | Feldstein | Jan 1989 | A |
4798590 | O'Leary et al. | Jan 1989 | A |
4838856 | Mulreany et al. | Jun 1989 | A |
4840620 | Kobayashi et al. | Jun 1989 | A |
4853521 | Claeys et al. | Aug 1989 | A |
4854324 | Hirschman et al. | Aug 1989 | A |
4857056 | Talonn | Aug 1989 | A |
4874359 | White et al. | Oct 1989 | A |
4880014 | Karowitz et al. | Nov 1989 | A |
4887554 | Whitford | Dec 1989 | A |
4903705 | Imamura et al. | Feb 1990 | A |
4913154 | Ermert et al. | Apr 1990 | A |
4925444 | Orkin et al. | May 1990 | A |
4943279 | Samiotes et al. | Jul 1990 | A |
4946439 | Eggers | Aug 1990 | A |
4966579 | Polaschegg | Oct 1990 | A |
5009654 | Minshall et al. | Apr 1991 | A |
5018173 | Komai et al. | May 1991 | A |
5032112 | Fairchild et al. | Jul 1991 | A |
5034987 | Fujimoto et al. | Jul 1991 | A |
5040537 | Katakura | Aug 1991 | A |
5078683 | Sancoff et al. | Jan 1992 | A |
5088981 | Howson et al. | Feb 1992 | A |
5104374 | Bishko et al. | Apr 1992 | A |
5111492 | Klausz | May 1992 | A |
5128121 | Berg et al. | Jul 1992 | A |
5135000 | Akselrod et al. | Aug 1992 | A |
5180896 | Gibby et al. | Jan 1993 | A |
5301672 | Kalender | Apr 1994 | A |
5310997 | Roach et al. | May 1994 | A |
5317506 | Coutre et al. | May 1994 | A |
5328463 | Barton et al. | Jul 1994 | A |
5339799 | Kami et al. | Aug 1994 | A |
5349625 | Born et al. | Sep 1994 | A |
5361761 | Van Lysel et al. | Nov 1994 | A |
5368562 | Blomqquist et al. | Nov 1994 | A |
5378231 | Johnson et al. | Jan 1995 | A |
5382232 | Hague et al. | Jan 1995 | A |
5383858 | Reilly et al. | Jan 1995 | A |
5385540 | Abbott et al. | Jan 1995 | A |
5400792 | Hoebel et al. | Mar 1995 | A |
5417213 | Prince | May 1995 | A |
5431627 | Pastrone et al. | Jul 1995 | A |
5445621 | Poli et al. | Aug 1995 | A |
5456255 | Abe et al. | Oct 1995 | A |
5458128 | Polanyi et al. | Oct 1995 | A |
5459769 | Brown | Oct 1995 | A |
5469769 | Sawada et al. | Nov 1995 | A |
5472403 | Cornacchia et al. | Dec 1995 | A |
5474683 | Bryant et al. | Dec 1995 | A |
5485831 | Holdsworth et al. | Jan 1996 | A |
5494036 | Uber, III et al. | Feb 1996 | A |
5494822 | Sadri | Feb 1996 | A |
5496273 | Pastrone et al. | Mar 1996 | A |
5507412 | Ebert et al. | Apr 1996 | A |
5515851 | Goldstein | May 1996 | A |
5522798 | Johnson et al. | Jun 1996 | A |
5531679 | Schulman et al. | Jul 1996 | A |
5569181 | Heilman et al. | Oct 1996 | A |
5573515 | Wilson et al. | Nov 1996 | A |
5583902 | Bae | Dec 1996 | A |
5632623 | Kolff et al. | May 1997 | A |
5687208 | Bae et al. | Nov 1997 | A |
5713358 | Mistretta et al. | Feb 1998 | A |
5724976 | Mine et al. | Mar 1998 | A |
5743266 | Levene et al. | Apr 1998 | A |
5768405 | Makram-Ebeid | Jun 1998 | A |
5840026 | Uber et al. | Nov 1998 | A |
5865744 | Lemelson | Feb 1999 | A |
5881124 | Giger et al. | Mar 1999 | A |
5916197 | Reilly et al. | Jun 1999 | A |
6009342 | Brasch et al. | Dec 1999 | A |
6055985 | Bae et al. | May 2000 | A |
6062866 | Prom | May 2000 | A |
6197000 | Reilly et al. | Mar 2001 | B1 |
6201889 | Vannah | Mar 2001 | B1 |
6205871 | Saloner et al. | Mar 2001 | B1 |
6236706 | Hsieh | May 2001 | B1 |
6317623 | Griffiths et al. | Nov 2001 | B1 |
6344030 | Duchon et al. | Feb 2002 | B1 |
6385483 | Uber, III et al. | May 2002 | B1 |
6397098 | Uber, III et al. | May 2002 | B1 |
6470889 | Bae et al. | Oct 2002 | B1 |
6471674 | Emig et al. | Oct 2002 | B1 |
6478735 | Pope et al. | Nov 2002 | B1 |
6522913 | Swanson et al. | Feb 2003 | B2 |
6556695 | Packer et al. | Apr 2003 | B1 |
6626862 | Duchon et al. | Sep 2003 | B1 |
6635030 | Bae et al. | Oct 2003 | B1 |
6643537 | Zatezalo et al. | Nov 2003 | B1 |
6656157 | Wilson et al. | Dec 2003 | B1 |
6673033 | Sciulli et al. | Jan 2004 | B1 |
6731971 | Evans, III et al. | May 2004 | B2 |
6776764 | Pinsky | Aug 2004 | B2 |
6901283 | Evans, III et al. | May 2005 | B2 |
6970735 | Uber, III et al. | Nov 2005 | B2 |
7094216 | Trombley, III et al. | Aug 2006 | B2 |
7326186 | Trombley, III et al. | Feb 2008 | B2 |
7925330 | Kalafut et al. | Apr 2011 | B2 |
7996381 | Uber, III et al. | Aug 2011 | B2 |
8016598 | Eggert et al. | Sep 2011 | B2 |
8055328 | Uber, III et al. | Nov 2011 | B2 |
8197437 | Kalafut et al. | Jun 2012 | B2 |
8295914 | Kalafut et al. | Oct 2012 | B2 |
20030135111 | Meaney et al. | Jul 2003 | A1 |
20030212364 | Mann et al. | Nov 2003 | A1 |
20030216683 | Shekalim | Nov 2003 | A1 |
20040010229 | Houde et al. | Jan 2004 | A1 |
20040025452 | McLean | Feb 2004 | A1 |
20040033477 | Ramphal et al. | Feb 2004 | A1 |
20040064040 | Masuda et al. | Apr 2004 | A1 |
20040097875 | Bae | May 2004 | A1 |
20040162484 | Nemoto | Aug 2004 | A1 |
20040215144 | Duchon et al. | Oct 2004 | A1 |
20040254525 | Uber et al. | Dec 2004 | A1 |
20070282263 | Kalafut et al. | Dec 2007 | A1 |
20080293025 | Zamierowsi et al. | Nov 2008 | A1 |
20080293029 | Wilkins et al. | Nov 2008 | A1 |
Number | Date | Country |
---|---|---|
2045070 | Feb 1992 | CA |
0619122 | Oct 1994 | EP |
60194934 | Oct 1985 | JP |
8500292 | Jan 1985 | WO |
8803815 | Jun 1988 | WO |
9415664 | Jul 1994 | WO |
9820919 | May 1998 | WO |
0061216 | Oct 2000 | WO |
Entry |
---|
Awai, K. et al, “Simulation of Aortic Peak Enhancement on MDCT using a Contrast Material Flow Phantom, feasibility study”, AJR Am J. Roentgenol, 2006. 186(2): pp. 379-385. |
Kalafut, J.F. et al, “An Anthropomorphic, Cardiovascular Flow Phantom for CT Contrast and Imaging Protocol Optimization”, ECR 2008, EPOS, Scientific Exhibit, pp. 13. |
Bae, K.T. et al, “Peak Contrast Enhancement in CT and MR Angiography: When Does it Occur and Why? Pharmacokinetic Study in a Porcine Model”, Radiology, Jun. 2003, pp. 809-816. |
Cademartini, F. and Luccichenti, G., et al (2004) “Sixteen-Row Multislice Computed Tomography: Basic Concepts, Protocols, and Enhanced Clinical Application,” Semir Ultrasound CT MR 25(1):2-16. |
K.T. Bae, J.P. Heiken and J.A. Brink, “Aortic and Hepatic Contrast Medium Enhancement at CT Part I. Prediction and a Computer Model.” Radiology, vol. 207, pp. 647-655, 1998. |
K.T. Bae, “Peak Contrast Enhancement in CT and MR Angiography: When Does It Occur and Why? Pharmacokinetic Study in a Porcine Model”, Radiology. vol. 227, pp. 809-816, 2003. |
K.T. Bae et al., “Multiphasic Injection Method for Uniform Prolonged Vascular Enhancement at CT Angiography: Pharmacokinetic Analysis and Experimental Porcine Method,” Radiology, vol. 216, pp. 872-880, 2000. |
D. Fleischmann and K. Hittmair, “Mathematical Analysis of Arterial Enhancement and Optimization of Bolus Geometry for CT Angiography Using the Discrete Fourier Transform”, J Comput Assist Tomogr, vol. 23, pp. 474-484, 1999. |
Fisher and Teo, “Optimal Insulin Infusion Resulting from a Mathematical Model of Blood Glucuose Dynamics,” IEEE Trans Biomed Eng. vol. 36(4). pp. 479-486, 1989. |
Jacobs, “Algorithm for Optimal Linear Model-Based Control with Application to Pharmacokinetic Model-Driven Drug Delivery”, IEEE Trans Biomed Eng. vol. 37(1), pp. 107-109, 1990. |
Wada and Ward, “The Hybrid Model: A New Pharmacokinetic Model for Computer-Controlled Infusion Pumps”, IEEE Trans. Biomed.Eng. vol. 41(2), pp. 134-142, 1994. |
Wada and Ward, “Open Loop Control of Multiple Drug Effects in Anesthesia”, IEEE Trans. Biomed Eng. vol. 42(7), pp. 666-677, 1995. |
Neatpisamvanit and Boston, “Estimation of Plasma Insulin from Plasma Glucose”, IEEE Trans Biomed Eng. vol. 49 (11), pp. 1253-1259, 2002. |
Gentilini et al., “A New Paradigm for the Closed-Loop Intraoperative Administration of Analgesics in Humans”, IEEE Tran Biomed Eng. vol. 49(4). pp. 289-2999, 2002. |
Garrett, J.S., Lanzer, P., et al., “Measurement of Cardiac Output by Cine Computed Tomography”, Am J Cardiol 56 (10): 657-61, 1985. |
Mahnken, A.H., Henzler, D, et al., “Determination of Cardiac Output with Multislice Spiral Computed Tomography: A Validation Study”, Invest Radiol 39(8): 451-4, 2004. |
Mahnken, A.H., Klotz, E, et al., “Measurement of Cardiac Output from a Test-Bolus Injection in Multislice Computed Tomography”, Eur Radiol 13(11): 2498-504, 2003. |
International Search Report and Written Opinion for counterpart PCT Application No. PCT/US2005/41913, 2004. |
Fleischmann, Dominik, “Present and Future Trends in Multiple Detector-Row CT Applications; CT Angiography”, Eur. Radiol. 12 (Suppl. 2) 2002. |
Korosec, Frank, “Basic Principles of Phase-contrast, Time-of-flight, and Contrast-enhanced MR Angiography”, 1999. |
Baker, Aaron; et al. “Fluid Mechanics Analysis of a Spring-Loaded Jet Injector.” IEEE Transactions on Biomedical Engineering, vol. 46, No. 2, Feb. 1899. |
Guytan, A.C., “Circulatory Physiology: Cardiac Output and Regulation”, Saunders, Philadelphia, p. 173, ISBN: 07216436004, Nov. 2004. |
Becker, C.R., et al., Optimal contrast application for cardiac 4-detector-row computed tomography. Invest Radiol, 2003. 38(11): p. 690-4. |
Schoepf, U.J., et al., Coronary CT angiography. Radiology, 2007. 244(1): p. 48-63. |
Johnson, T.R. et al., Optimization of contrast material administration for electrocardiogram-gated, 2007. |
Fleischmann, D. and K. Hittmair, Mathematical analysis of arterial enhancement and optimization of bolus geometry for CT angiography using the discrete fourier transform. J Com put Assist Tomogr, 1999. 23(3):p. 474-84. |
Awai, K., et al., Prediction of aortic peak enhancement in monophasic contrast injection protocols at multidetector CT; phantom and patient studies. Radiat Med, 2007.25(1 ):p. 14-21. |
Bae, K.T., J.P. Heiken, and J.A. Brink, Aortic and hepatic contrast medium enhancement at CT. Part II. Effect of reduced cardiac output in a porcine model. Radiology, 1998. 207(3):p. 657-62. |
Number | Date | Country | |
---|---|---|---|
20090226867 A1 | Sep 2009 | US |
Number | Date | Country | |
---|---|---|---|
61033761 | Mar 2008 | US |