Identification of ultrasound contrast agent dilution systems for ejection fraction measurements

Contrast-Enhanced Ultrasound Quantification: From Kinetic Modeling to Machine Learning

Ultrasound contrast agents (UCAs) have opened up immense diagnostic possibilities by combined use of indicator dilution principles and dynamic contrast-enhanced ultrasound (DCE-US) imaging. UCAs are microbubbles encapsulated in a biocompatible shell. With a rheology comparable to that of red blood cells, UCAs provide an intravascular indicator for functional imaging of the (micro)vasculature by quantitative DCE-US. Several models of the UCA intravascular kinetics have been proposed to provide functional quantitative maps, aiding diagnosis of different pathological conditions. This article is a comprehensive review of the available methods for quantitative DCE-US imaging based on temporal, spatial and spatiotemporal analysis of the UCA kinetics. The recent introduction of novel UCAs that are targeted to specific vascular receptors has advanced DCE-US to a molecular imaging modality. In parallel, new kinetic models of increased complexity have been developed. The extraction of multiple quantitative maps, reflecting complementary variables of the underlying physiological processes, requires an integrative approach to their interpretation. A probabilistic framework based on emerging machine-learning methods represents nowadays the ultimate approach, improving the diagnostic accuracy of DCE-US imaging by optimal combination of the extracted complementary information. The current value and future perspective of all these advances are critically discussed.

Absolute ultrasound perfusion parameter quantification of a tissue-mimicking phantom using bolus tracking [Correspondence]

This study presents three methods for absolute quantification in ultrasound perfusion analysis based on bolus tracking. The first two methods deconvolve the perfusion time sequence with a measured AIF, using a nonparametric or a parametric model of the tissue residue function, respectively. The third method is a simplified approach avoiding deconvolution by assuming a narrow AIF. A phantom with a dialyzer filter as a tissue-mimicking model was used for evaluation. Estimated mean transit times and blood volumes were compared with the theoretical values. A match with a maximum error of 12% was achieved.

Ultrasound perfusion analysis combining bolus-tracking and burst-replenishment

A new signal model and processing method for quantitative ultrasound perfusion analysis is presented, called bolus-and-burst. The method has the potential to provide absolute values of blood flow, blood volume, and mean transit time. Furthermore, it provides an estimate of the local arterial input function which characterizes the arterial tree, allowing accurate estimation of the bolus arrival time. The method combines two approaches to ultrasound perfusion analysis: bolus-tracking and burst-replenishment. A pharmacokinetic model based on the concept of arterial input functions and tissue residue functions is used to model both the bolus and replenishment parts of the recording. The pharmacokinetic model is fitted to the data using blind deconvolution. A preliminary assessment of the new perfusion-analysis method is presented on clinical recordings.

Comparison and evaluation of indicator dilution models for bolus of ultrasound contrast agents

Dynamic contrast-enhanced ultrasound (DCE-US) imaging is a promising diagnostic method, which enables the evaluation of tissue perfusion via different parameters. The mean transit time and time-to-peak parameters are the main time parameters and their values depend on the model used for the approximation of the noisy perfusion curves. In this paper, we described a new comparison of different perfusion models using a tissue mimicking phantom. The following models were compared: log-normal, lagged, Erlang, Gamma and the local density random walk model. We discovered that the mean-square error is not the best criterion for model evaluation. More important is the comparison between the estimated time perfusion parameters and the physical parameters of the developed tissue mimicking phantom. Based on the statistical analysis, we can suggest that for the DCE-US perfusion analysis more models should be used, excluding the log-normal model, which gives the highest error of mean transit time value.

A method to expedite data acquisition for multiple spatial-temporal analyses of tissue perfusion by contrast-enhanced ultrasound

For semiquantitative analyses of tissue perfusion using contrast-enhanced ultrasound the acquisition and processing of time intensity curves (TIC) is required. These TICs can be computed for each pixel of an image plane, yielding parametric images of classification numbers like "blood volume" and "flow rate." The expenditure of time for data acquisition and analysis typically limits semiquantitative perfusion imaging to a single image plane in 2-D. 3-D techniques, however, provide a higher diagnostic value since more information (e.g., of an entire lesion) is obtained. Moreover, spatial compounding, being a 2-D-technique where an object is imaged from different viewing angles, is known to improve image quality by reducing artifacts and speckle noise. Both techniques, 3-D and compounding, call for optimized acquisition and processing of TICs in several image planes (3-D) or in several (overlapping) sections of the same image plane (compounding) to decrease the time needed for data acquisition. Here, an approach of interleaved imaging is presented which is applicable, among others, to contrast perfusion imaging using the replenishment method. The total acquisition time is decreased by sequentially scanning image planes twice for short time spans - first, immediately after microbubble destruction to record the initial rise of the TICs, and second, a sufficient time thereafter to assess final values of the TIC. Data from both periods are combined to fit a model function from which parameters are extracted such as perfusion rate and blood volume. This approach was evaluated by in vitro measurements on a perfusion-mimicking phantom for both, individual images such as would be used for volume reconstruction in 3-D and compound images obtained from full angle spatial compounding (FASC, 360 degrees ). An error analysis is conducted to derive the deviation of the extracted parameters of the proposed method compared with the conventional one. These deviations are entailed by a reduction in acquisition time of the proposed method, which can be adjusted by several parameters, depending on the prevailing flow. Optimization strategies are proposed to find optimal values for those settings.

Transpulmonary dilution system identification for pulmonary blood volume measurements by contrast echocardiography

Pulmonary blood volume (PBV) is an important parameter to assess the condition of the lungs and the transpulmonary circulation. Unfortunately, its measurement is limited by the need for central catheterization. Contrast ultrasonography allows overcoming this problem. A small bolus of ultrasound contrast agent is injected in a peripheral vein and its passage through the right ventricle and left atrium is detected non-invasively by an ultrasound transducer. The PBV is then given by the product of the mean transit time (MTT) of the contrast between the two measurement sites times the cardiac output. The MTT is estimated by specific model interpolation of the measured dilution curves. In this paper we present a new method for PBV measurements based on a system identification approach. This method identifies the parameters of the model that represents the dilution system impulse response. No subsequent model interpolation is needed. The local density random walk model is adopted to represent the transpulmonary dilution system. Volume measurements show accurate in-vitro results with a correlation coefficient higher than 0.99. The clinical feasibility is confirmed by 70 measurements in patients. Beyond an accurate quantification of pulmonary blood volume, the proposed method also permits the characterization of the transpulmonary hemodynamics, possibly adding novel diagnostic value to the measurement.

Identification of cardiovascular dilution systems by contrast ultrasound

Indicator dilution techniques permit accurate measurements of important cardiovascular parameters, such as pulmonary blood volume (PBV) and ejection fraction (EF). However, their use is limited by the need for central catheterization. Contrast ultrasonography allows overcoming this problem. PBV and EF can be measured by a dilution system identification algorithm after detection of multiple dilution curves by an ultrasound scanner. In this paper, we present a system identification method that exploits the a priori knowledge on the dilution system and finds the optimum parameters for the parametric model representing the dilution system impulse response. No subsequent model interpolation is needed. Volume measurements show accurate in-vitro results and clinical feasibility, while 50 EF measurements in patients show a 0.88 correlation coefficient with echocardiographic biplane estimates. In conclusion, adding a priori knowledge to the system identification algorithm leads to increased accuracy and robustness of the method for PBV and EF measurements.

Contrast ultrasound methods for left-ventricle ejection fraction measurements

The left-ventricle ejection fraction is an important cardiac-efficiency measure that is regularly used in cardiology. Standard estimations are based on time-consuming geometrical analysis and modelling, which requires experienced cardiologists. Alternative methods are very invasive due to the need for cardiac catheterization. In this paper we present and study a minimally-invasive indicator dilution technique for ejection fraction quantification that has recently been developed. It is based on a peripheral injection of an ultrasound contrast agent bolus. Left-atrium and left-ventricle acoustic intensities are recorded versus time by transthoracic echocardiography during contrast bolus passage. The measured curves are corrected for attenuation distortion, filtered to suppress the measurement noise, and processed by an adaptive Wiener deconvolution algorithm for the estimation of the left-ventricle impulse response. The estimated impulse response is interpolated by a mono-compartment exponential model for the ejection fraction assessment. An adaptive search of the interval for the model fitting is also included. The feasibility of the method is tested on 52 measurements in patients with left-ventricle ejection fractions between 10% and 80%. The results are promising and show a 0.83 correlation coefficient with echographic biplane ejection fraction measurements.

Bubble dynamics involved in ultrasonic imaging

In clinical ultrasound, blood cells cannot be differentiated from surrounding tissue, due to the low acoustic impedance difference between blood cells and their surroundings. Resonant gas bubbles introduced in the bloodstream are ideal markers, if rapid dissolution can be prevented. Ultrasound contrast agents consist of microscopically small bubbles encapsulated by an elastic shell. These microbubbles oscillate upon ultrasound insonification. Microbubbles with thin lipid shells have demonstrated highly nonlinear behavior. To enhance diagnostic ultrasound imaging techniques and to explore therapeutic applications, these medical microbubbles have been modeled. Several detection techniques have been proposed to improve the detectability of the microbubbles. A new generation of contrast agents, with special targeting ligands attached to the shells, may assist the imaging of nonphysical properties of target tissue. Owing to microbubble-based contrast agents, ultrasound is becoming an even more important technique in clinical diagnostics.

Quantification in echocardiography

Until recently, more than 2200 Swan Ganz catheters were used annually in the operating rooms (OR) and intensive care unit (ICU) of the Catharina Hospital in Eindhoven, The Netherlands. After cardiologists who were specialists in echocardiography (ECHO) trained anesthesiologists in ECHO, the need for these catheters in cardiac and noncardiac surgery was reduced. Initially intended as a local teaching project, an ECHO teaching compact disk (CD) was produced during the training and distributed later worldwide, thanks to a positive review in a major anesthesiology publication. By reducing the number of Swan Ganz catheters, the hospital could finance and acquire two echocardiography machines for the OR and ICU. The availability of these machines resulted in a further reduction of the number of Swan Ganz catheters. However, the need for quantification (eg, measurements of cardiac output) remained. During the creation of the ECHO teaching CD, the idea was born to apply indicator-dilution principles on injected echo contrast. This study was performed in cooperation with the Signal Processing Department of the Eindhoven University of Technology. Advanced signal processing and modelling were used to develop algorithms to enable quantification of intrapulmonary blood volume, ejection-fraction, and flow from the transesophageal echocardiography approach. These quantitative measurements, which can be performed on an outpatient basis, may become a real asset in cardiology, anesthesiology, and intensive care.