Evaluation of a numerical simulation for cryoablation - comparison with bench data, clinical kidney and lung cases.

PURPOSE: The accuracy of a numerical simulation of cryoablation ice balls was evaluated in gel phantom data as well as clinical kidney and lung cases. MATERIALS AND METHODS: To evaluate the accuracy, 64 experimental single-needle cryoablations and 12 multi-needle cryoablations in gel phantoms were re-simulated with the corresponding freeze-thaw-freeze cycles. The simulated temperatures were compared over time with the measurements of thermocouples. For single needles, temperature values were compared at each thermocouple location. For multiple needles, Euclidean distances between simulated and measured isotherms (10 °C, 0 °C, -20 °C, -40 °C) were computed. Furthermore, surface and volume of simulated 0 °C isotherms were compared to cryoablation-induced ice balls in 14 kidney and 13 lung patients. For this purpose, needle positions and relevant anatomical structures defining material parameters (kidney/lung, tumor) were reconstructed from pre-ablation CT images and fused with postablation CT images (from which ice balls were extracted by manual delineation). RESULTS: The single-needle gel phantom cases showed less than 5 °C prediction error on average. Over all multiple needle experiments in gel, the mean and maximum isotherm distance were less than 2.3 mm and 4.1 mm, respectively. Average Dice coefficients of 0.82/0.63 (kidney/lung) and mean surface distances of 2.59/3.12 mm quantify the prediction performance of the numerical simulation. However, maximum surface distances of 10.57/10.8 mm indicate that locally larger errors have to be expected. CONCLUSION: A very good agreement of the numerical simulations for gel experiments was measured and a satisfactory agreement of the numerical simulations with measured ice balls in patient data was shown.

MRI-based computational hemodynamics in patients with aortic coarctation using the lattice Boltzmann methods: Clinical validation study.

PURPOSE: To introduce a scheme based on a recent technique in computational hemodynamics, known as the lattice Boltzmann methods (LBM), to noninvasively measure pressure gradients in patients with a coarctation of the aorta (CoA). To provide evidence on the accuracy of the proposed scheme, the computed pressure drop values are compared against those obtained using the reference standard method of catheterization. MATERIALS AND METHODS: Pre- and posttreatment LBM-based pressure gradients for 12 patients with CoA were simulated for the time point of peak systole using the open source library OpenLB. Four-dimensional (4D) flow-sensitive phase-contrast MRI at 1.5 Tesla was used to acquire flow and to setup the simulation. The vascular geometry was reconstructed using 3D whole-heart MRI. Patients underwent pre- and postinterventional pressure catheterization as a reference standard. RESULTS: There is a significant linear correlation between the pretreatment catheter pressure drops and those computed based on the LBM simulation, r=.85, P<.001. The bias was -0.58 ± 4.1 mmHg and was not significant ( P=0.64) with a 95% confidence interval (CI) of -3.22 to 2.06. For the posttreatment results, the bias was larger and at -2.54 ± 3.53 mmHg with a 95% CI of -0.17 to -4.91 mmHg. CONCLUSION: The results indicate a reasonable agreement between the simulation results and the catheter measurements. LBM-based computational hemodynamics can be considered as an alternative to more traditional computational fluid dynamics schemes for noninvasive pressure calculations and can assist in diagnosis and therapy planning. LEVEL OF EVIDENCE: 3 J. Magn. Reson. Imaging 2017;45:139-146.

Algorithmically generated rodent hepatic vascular trees in arbitrary detail.

Physiologically realistic geometric models of the vasculature in the liver are indispensable for modelling hepatic blood flow, the main connection between the liver and the organism. Current in vivo imaging techniques do not provide sufficiently detailed vascular trees for many simulation applications, so it is necessary to use algorithmic refinement methods. The method of Constrained Constructive Optimization (CCO) (Schreiner et al., 2006) is well suited for this purpose. Its results after calibration have been previously compared to experimentally acquired human vascular trees (Schwen and Preusser, 2012). The goal of this paper is to extend this calibration to the case of rodents (mice and rats), the most commonly used animal models in liver research. Based on in vivo and ex vivo micro-CT scans of rodent livers and their vasculature, we performed an analysis of various geometric features of the vascular trees. Starting from pruned versions of the original vascular trees, we applied the CCO procedure and compared these algorithmic results to the original vascular trees using a suitable similarity measure. The calibration of the postprocessing improved the algorithmic results compared to those obtained using standard CCO. In terms of angular features, the average similarity increased from 0.27 to 0.61, improving the total similarity from 0.28 to 0.40. Finally, we applied the calibrated algorithm to refine measured vascular trees to the (higher) level of detail desired for specific applications. Having successfully adapted the CCO algorithm to the rodent model organism, the resulting individual-specific refined hepatic vascular trees can now be used for advanced modeling involving, e.g., detailed blood flow simulations.

Spatio-temporal simulation of first pass drug perfusion in the liver.

The liver is the central organ for detoxification of xenobiotics in the body. In pharmacokinetic modeling, hepatic metabolization capacity is typically quantified as hepatic clearance computed as degradation in well-stirred compartments. This is an accurate mechanistic description once a quasi-equilibrium between blood and surrounding tissue is established. However, this model structure cannot be used to simulate spatio-temporal distribution during the first instants after drug injection. In this paper, we introduce a new spatially resolved model to simulate first pass perfusion of compounds within the naive liver. The model is based on vascular structures obtained from computed tomography as well as physiologically based mass transfer descriptions obtained from pharmacokinetic modeling. The physiological architecture of hepatic tissue in our model is governed by both vascular geometry and the composition of the connecting hepatic tissue. In particular, we here consider locally distributed mass flow in liver tissue instead of considering well-stirred compartments. Experimentally, the model structure corresponds to an isolated perfused liver and provides an ideal platform to address first pass effects and questions of hepatic heterogeneity. The model was evaluated for three exemplary compounds covering key aspects of perfusion, distribution and metabolization within the liver. As pathophysiological states we considered the influence of steatosis and carbon tetrachloride-induced liver necrosis on total hepatic distribution and metabolic capacity. Notably, we found that our computational predictions are in qualitative agreement with previously published experimental data. The simulation results provide an unprecedented level of detail in compound concentration profiles during first pass perfusion, both spatio-temporally in liver tissue itself and temporally in the outflowing blood. We expect our model to be the foundation of further spatially resolved models of the liver in the future.

An integrated model-based software for FUS in moving abdominal organs.

Focused ultrasound surgery (FUS) is a non-invasive method for tissue ablation that has the potential for complete and controlled local tumour destruction with minimal side effects. The treatment of abdominal organs such as the liver, however, requires particular technological support in order to enable a safe, efficient and effective treatment. As FUS is applied from outside the patient's body, suitable imaging methods, such as magnetic resonance imaging or diagnostic ultrasound, are needed to guide and track the procedure. To facilitate an efficient FUS procedure in the liver, the organ motion during breathing and the partial occlusion by the rib cage need to be taken into account in real time, demanding a continuous patient-specific adaptation of the treatment configuration. Modelling the patient's respiratory motion and combining this with tracking data improves the accuracy of motion predictions. Modelling and simulation of the FUS effects within the body allows the use of treatment planning and has the potential to be used within therapy to increase knowledge about the patient status. This article describes integrated model-based software for patient-specific modelling and prediction for FUS treatments of moving abdominal organs.

Toward focused ultrasound liver surgery under free breathing.

Focused ultrasound surgery is an outstanding novel technique for cancer treatment because it is completely noninvasive with the potential for complete and controlled local tumor destruction. Because focused ultrasound surgery is applied from outside of the patient's body, imaging such as ultrasound or magnetic resonance imaging is required to plan and monitor the intervention. For the treatment of liver tumors, several complexities have to be taken into account, including accessibility of the target and protection of structures at risk. To allow for safe and efficient treatment under free respiration, in which the liver moves significantly, both planning and execution have to be performed specifically according to the patient's individual breathing. This article reviews the state of the art of liver applications, the tremendous challenges of this field, and approaches to overcome these challenges. This includes modeling of the patient-individual breathing cycle, detection of and adaptation to the actual breathing, and simulation and monitoring of the therapy.

Instrument localisation for endovascular aneurysm repair: Comparison of two methods based on tracking systems or using imaging.

BACKGROUND: In endovascular aneuysm repair (EVAR) procedures, medical instruments are currently navigated with a two-dimensional imaging based guidance requiring X-rays and contrast agent. METHODS: Novel approaches for obtaining the three-dimensional instrument positions are introduced. Firstly, a method based on fibre optical shape sensing, one electromagnetic sensor and a preoperative computed tomography (CT) scan is described. Secondly, an approach based on image processing using one 2D fluoroscopic image and a preoperative CT scan is introduced. RESULTS: For the tracking based method, average errors from 1.81 to 3.13 mm and maximum errors from 3.21 to 5.46 mm were measured. For the image-based approach, average errors from 3.07 to 6.02 mm and maximum errors from 8.05 to 15.75 mm were measured. CONCLUSION: The tracking based method is promising for usage in EVAR procedures. For the image-based approach are applications in smaller vessels more suitable, since its errors increase with the vessel diameter.

State of the art in computer-assisted planning, intervention, and assessment of liver-tumor ablation.

Percutaneous, image-guided thermal tumor ablation procedures are used increasingly for minimally invasive, local treatment of tumors in the liver. The planning of these procedures; the support of targeting, monitoring, and controlling during the intervention itself; and the assessment of the treatment response can all benefit significantly from computer assistance. The outcome can be optimized by supporting the physician in the process of determining an intervention strategy that enables complete destruction of the targeted tumor while reducing the danger of complications. During the intervention, computer-assisted methods can be used to guide the physician in the implementation of the intended strategy by providing planning information. Assessment of the intervention result is carried out by comparison of the achieved coagulation with the target tumor volume. Supporting this comparison facilitates the early detection of potential recurrences. This report provides an overview of state-of-the-art computer-assisted methods for the support of thermal tumor ablations in the liver. Proper approaches for image segmentation, access-path determination, simulation, visualization, interventional guidance, and post-interventional assessment, as well as integrated work flow-oriented solutions, are reviewed with respect to technical aspects and applicability in the clinical setting.

Landmark-based evaluation of a deformable motion correction for DCE-MRI of the liver.

PURPOSE: Annotation of meaningful landmark ground truth on DCE-MRI is difficult and laborious. Motion correction methods applied to DCE-MRI of the liver are thus mostly evaluated using qualitative or indirect measures. We propose a novel landmark annotation scheme that facilitates the generation of landmark ground truth on larger clinical datasets. METHODS: In our annotation scheme, landmarks are equally distributed over all time points of all available dataset cases and annotated by multiple observers on a per-pair basis. The scheme is used to annotate 26 DCE-MRI of the liver. A subset of the ground truth is used to optimize parameters of a deformable motion correction. Several variants of the motion correction are evaluated on the remaining cases with respect to distances of corresponding landmarks after registration, deformation field properties, and qualitative measures. RESULTS: A landmark ground truth on 26 cases could be generated in under 12 h per observer with a mean inter-observer distance below the mean voxel diagonal. Furthermore, the landmarks are spatially well distributed within the liver. Parameter optimization significantly improves the performance of the motion correction, and landmark distance after registration is 2 mm. Qualitative evaluation of the motion correction reflects the quantitative results. CONCLUSIONS: The annotation scheme makes a landmark-based evaluation of motion corrections for hepatic DCE-MRI practically feasible for larger clinical datasets. The comparably large number of cases enables both optimization and evaluation of motion correction methods.

Multiscale optimization of the probe placement for radiofrequency ablation.

RATIONALE AND OBJECTIVES: We present a model for the optimal placement of mono- and bipolar probes in radiofrequency (RF) ablation. The model is based on a system of partial differential equations that describe the electric potential of the probe and the steady state of the induced heat distribution. MATERIALS AND METHODS: To optimize the probe placement we minimize a temperature-based objective function under the constraining system of partial differential equations. Further, the extension of the resulting optimality system for the use of multiple coupled RF probes is discussed. We choose a multiscale gradient descent approach to solve the optimality system. RESULTS: This article describes the discretization and implementation of the approach with finite elements on three-dimensional hexahedral grids. CONCLUSION: Applications of the optimization to artificial test scenarios as well as a comparison to a real RF ablation show the usefulness of the approach.

Fast Numerical Simulation of Focused Ultrasound Treatments During Respiratory Motion With Discontinuous Motion Boundaries.

OBJECTIVE: Focused ultrasound (FUS) is rapidly gaining clinical acceptance for several target tissues in the human body. Yet, treating liver targets is not clinically applied due to a high complexity of the procedure (noninvasiveness, target motion, complex anatomy, blood cooling effects, shielding by ribs, and limited image-based monitoring). To reduce the complexity, numerical FUS simulations can be utilized for both treatment planning and execution. These use-cases demand highly accurate and computationally efficient simulations. METHODS: We propose a numerical method for the simulation of abdominal FUS treatments during respiratory motion of the organs and target. Especially, a novel approach is proposed to simulate the heating during motion by solving Pennes' bioheat equation in a computational reference space, i.e., the equation is mathematically transformed to the reference. The approach allows for motion discontinuities, e.g., the sliding of the liver along the abdominal wall. RESULTS: Implementing the solver completely on the graphics processing unit and combining it with an atlas-based ultrasound simulation approach yields a simulation performance faster than real time (less than 50-s computing time for 100 s of treatment time) on a modern off-the-shelf laptop. The simulation method is incorporated into a treatment planning demonstration application that allows to simulate real patient cases including respiratory motion. CONCLUSION: The high performance of the presented simulation method opens the door to clinical applications. SIGNIFICANCE: The methods bear the potential to enable the application of FUS for moving organs.

Methodology on quantification of sonication duration for safe application of MR guided focused ultrasound for liver tumour ablation.

BACKGROUND AND OBJECTIVE: Magnetic Resonance Guided Focused Ultrasound (MRgFUS) for liver tumour ablation is a challenging task due to motion caused by breathing and occlusion due the ribcage between the transducer and the tumour. To overcome these challenges, a novel system for liver tumour ablation during free breathing has been designed. METHODS: The novel TRANS-FUSIMO Treatment System (TTS, EUFP7) interacts with a Magnetic Resonance (MR) scanner and a focused ultrasound transducer to sonicate to a moving target in liver. To meet the requirements of ISO 13485; a quality management system for medical device design, the system needs to be tested for certain process parameters. The duration of sonication and, the delay after the sonication button is activated, are among the parameters that need to be quantified for efficient and safe ablation of tumour tissue. A novel methodology is developed to quantify these process parameters. A computerised scope is programmed in LabVIEW to collect data via hydrophone; where the coordinates of fiber-optic sensor assembly was fed into the TRANS-FUSIMO treatment software via Magnetic Resonance Imaging (MRI) to sonicate to the tip of the sensor, which is synchronised with the clock of the scope, embedded in a degassed water tank via sensor assembly holder. The sonications were executed for 50 W, 100 W, 150 W for 10 s to quantify the actual sonication duration and the delay after the emergency stop by two independent operators for thirty times. The deviation of the system from the predefined specs was calculated. Student's-T test was used to investigate the user dependency. RESULTS: The duration of sonication and the delay after the sonication were quantified successfully with the developed method. TTS can sonicate with a maximum deviation of 0.16 s (Std 0.32) from the planned duration and with a delay of 14 ms (Std 0.14) for the emergency stop. Student's T tests indicate that the results do not depend on operators (p > .05). CONCLUSION: The evidence obtained via this protocol is crucial for translation- of-research into the clinics for safe application of MRgFUS. The developed protocol could be used for system maintenance in compliance with quality systems in clinics for daily quality assurance routines.

A focused ultrasound treatment system for moving targets (part I): generic system design and in-silico first-stage evaluation.

BACKGROUND: Focused ultrasound (FUS) is entering clinical routine as a treatment option. Currently, no clinically available FUS treatment system features automated respiratory motion compensation. The required quality standards make developing such a system challenging. METHODS: A novel FUS treatment system with motion compensation is described, developed with the goal of clinical use. The system comprises a clinically available MR device and FUS transducer system. The controller is very generic and could use any suitable MR or FUS device. MR image sequences (echo planar imaging) are acquired for both motion observation and thermometry. Based on anatomical feature tracking, motion predictions are estimated to compensate for processing delays. FUS control parameters are computed repeatedly and sent to the hardware to steer the focus to the (estimated) target position. All involved calculations produce individually known errors, yet their impact on therapy outcome is unclear. This is solved by defining an intuitive quality measure that compares the achieved temperature to the static scenario, resulting in an overall efficiency with respect to temperature rise. To allow for extensive testing of the system over wide ranges of parameters and algorithmic choices, we replace the actual MR and FUS devices by a virtual system. It emulates the hardware and, using numerical simulations of FUS during motion, predicts the local temperature rise in the tissue resulting from the controls it receives. RESULTS: With a clinically available monitoring image rate of 6.67 Hz and 20 FUS control updates per second, normal respiratory motion is estimated to be compensable with an estimated efficiency of 80%. This reduces to about 70% for motion scaled by 1.5. Extensive testing (6347 simulated sonications) over wide ranges of parameters shows that the main source of error is the temporal motion prediction. A history-based motion prediction method performs better than a simple linear extrapolator. CONCLUSIONS: The estimated efficiency of the new treatment system is already suited for clinical applications. The simulation-based in-silico testing as a first-stage validation reduces the efforts of real-world testing. Due to the extensible modular design, the described approach might lead to faster translations from research to clinical practice.

Zonated quantification of steatosis in an entire mouse liver.

Many physiological processes and pathological conditions in livers are spatially heterogeneous, forming patterns at the lobular length scale or varying across the organ. Steatosis, a common liver disease characterized by lipids accumulating in hepatocytes, exhibits heterogeneity at both these spatial scales. The main goal of the present study was to provide a method for zonated quantification of the steatosis patterns found in an entire mouse liver. As an example application, the results were employed in a pharmacokinetics simulation. For the analysis, an automatic detection of the lipid vacuoles was used in multiple slides of histological serial sections covering an entire mouse liver. Lobuli were determined semi-automatically and zones were defined within the lobuli. Subsequently, the lipid content of each zone was computed. The steatosis patterns were found to be predominantly periportal, with a notable organ-scale heterogeneity. The analysis provides a quantitative description of the extent of steatosis in unprecedented detail. The resulting steatosis patterns were successfully used as a perturbation to the liver as part of an exemplary whole-body pharmacokinetics simulation for the antitussive drug dextromethorphan. The zonated quantification is also applicable to other pathological conditions that can be detected in histological images. Besides being a descriptive research tool, this quantification could perspectively complement diagnosis based on visual assessment of histological images.

Interactive multi-criteria planning for radiofrequency ablation.

PURPOSE: Image-guided radiofrequency ablation (RFA) is a broadly used minimally invasive method for the thermal destruction of focal liver malignancies using needle-shaped instruments. The established planning workflow is based on examination of 2D slices and manual definition of the access path. During that process, multiple criteria for all possible trajectories have to be taken into account. Hence, it demands considerable experience and constitutes a significant mental task. METHODS: An access path determination method based on image processing and numerical optimization is proposed. Fast GPU-based simulation approximation is utilized to incorporate the heat distribution including realistic cooling effects from nearby blood vessels. A user interface for intuitive exploration of the optimization results is introduced. RESULTS: The proposed methods are integrated into a clinical software assistant. To evaluate the suitability of the interactive optimization approach for the identification of meaningful therapy strategies, a retrospective study has been carried out. The system is able to propose clinically relevant trajectories to the target by incorporating multiple criteria. CONCLUSIONS: A novel method for planning of image-guided radiofrequency ablation by means of interactive access path determination based on optimization is presented. A first retrospective study indicates that the method is suited to improve the classical planning of RFA.

Representative Sinusoids for Hepatic Four-Scale Pharmacokinetics Simulations.

The mammalian liver plays a key role for metabolism and detoxification of xenobiotics in the body. The corresponding biochemical processes are typically subject to spatial variations at different length scales. Zonal enzyme expression along sinusoids leads to zonated metabolization already in the healthy state. Pathological states of the liver may involve liver cells affected in a zonated manner or heterogeneously across the whole organ. This spatial heterogeneity, however, cannot be described by most computational models which usually consider the liver as a homogeneous, well-stirred organ. The goal of this article is to present a methodology to extend whole-body pharmacokinetics models by a detailed liver model, combining different modeling approaches from the literature. This approach results in an integrated four-scale model, from single cells via sinusoids and the organ to the whole organism, capable of mechanistically representing metabolization inhomogeneity in livers at different spatial scales. Moreover, the model shows circulatory mixing effects due to a delayed recirculation through the surrounding organism. To show that this approach is generally applicable for different physiological processes, we show three applications as proofs of concept, covering a range of species, compounds, and diseased states: clearance of midazolam in steatotic human livers, clearance of caffeine in mouse livers regenerating from necrosis, and a parameter study on the impact of different cell entities on insulin uptake in mouse livers. The examples illustrate how variations only discernible at the local scale influence substance distribution in the plasma at the whole-body level. In particular, our results show that simultaneously considering variations at all relevant spatial scales may be necessary to understand their impact on observations at the organism scale.

A detailed physiologically based model to simulate the pharmacokinetics and hormonal pharmacodynamics of enalapril on the circulating endocrine Renin-Angiotensin-aldosterone system.

The renin-angiotensin-aldosterone system (RAAS) plays a key role in the pathogenesis of cardiovascular disorders including hypertension and is one of the most important targets for drugs. A whole body physiologically based pharmacokinetic (wb PBPK) model integrating this hormone circulation system and its inhibition can be used to explore the influence of drugs that interfere with this system, and thus to improve the understanding of interactions between drugs and the target system. In this study, we describe the development of a mechanistic RAAS model and exemplify drug action by a simulation of enalapril administration. Enalapril and its metabolite enalaprilat are potent inhibitors of the angiotensin-converting-enzyme (ACE). To this end, a coupled dynamic parent-metabolite PBPK model was developed and linked with the RAAS model that consists of seven coupled PBPK models for aldosterone, ACE, angiotensin 1, angiotensin 2, angiotensin 2 receptor type 1, renin, and prorenin. The results indicate that the model represents the interactions in the RAAS in response to the pharmacokinetics (PK) and pharmacodynamics (PD) of enalapril and enalaprilat in an accurate manner. The full set of RAAS-hormone profiles and interactions are consistently described at pre- and post-administration steady state as well as during their dynamic transition and show a good agreement with literature data. The model allows a simultaneous representation of the parent-metabolite conversion to the active form as well as the effect of the drug on the hormone levels, offering a detailed mechanistic insight into the hormone cascade and its inhibition. This model constitutes a first major step to establish a PBPK-PD-model including the PK and the mode of action (MoA) of a drug acting on a dynamic RAAS that can be further used to link to clinical endpoints such as blood pressure.

Segmentation of stochastic images with a stochastic random walker method.

We present an extension of the random walker segmentation to images with uncertain gray values. Such gray-value uncertainty may result from noise or other imaging artifacts or more general from measurement errors in the image acquisition process. The purpose is to quantify the influence of the gray-value uncertainty onto the result when using random walker segmentation. In random walker segmentation, a weighted graph is built from the image, where the edge weights depend on the image gradient between the pixels. For given seed regions, the probability is evaluated for a random walk on this graph starting at a pixel to end in one of the seed regions. Here, we extend this method to images with uncertain gray values. To this end, we consider the pixel values to be random variables (RVs), thus introducing the notion of stochastic images. We end up with stochastic weights for the graph in random walker segmentation and a stochastic partial differential equation (PDE) that has to be solved. We discretize the RVs and the stochastic PDE by the method of generalized polynomial chaos, combining the recent developments in numerical methods for the discretization of stochastic PDEs and an interactive segmentation algorithm. The resulting algorithm allows for the detection of regions where the segmentation result is highly influenced by the uncertain pixel values. Thus, it gives a reliability estimate for the resulting segmentation, and it furthermore allows determining the probability density function of the segmented object volume.