Adaptive simulation of 3D thermometry maps for interventional MR-guided tumor ablation using Pennes' bioheat equation and isotherms.

Minimally-invasive thermal ablation procedures have become clinically accepted treatment options for tumors and metastases. Continuous and reliable monitoring of volumetric heat distribution promises to be an important condition for successful outcomes. In this work, an adaptive bioheat transfer simulation of 3D thermometry maps is presented. Pennes' equation model is updated according to temperature maps generated by uniformly distributed 2D MR phase images rotated around the main axis of the applicator. The volumetric heat diffusion and the resulting shape of the ablation zone can be modelled accurately without introducing a specific heat source term. Filtering the temperature maps by extracting isotherms reduces artefacts and noise, compresses information of the measured data and adds physical a priori knowledge. The inverse heat transfer for estimating values of the simulated tissue and heating parameters is done by reducing the sum squared error between these isotherms and the 3D simulation. The approach is evaluated on data sets consisting of 13 ex vivo bio protein phantoms, including six perfusion phantoms with simulated heat sink effects. Results show an overall average Dice score of 0.89 ± 0.04 (SEM < 0.01). The optimization of the parameters takes 1.05 ± 0.26 s for each acquired image. Future steps should consider the local optimization of the simulation parameters instead of a global one to better detect heat sinks without a priori knowledge. In addition, the use of a proper Kalman filter might increase robustness and accuracy if combined with our method.

Comparison study of reconstruction algorithms for volumetric necrosis maps from 2D multi-slice GRE thermometry images.

Cancer is a disease which requires a significant amount of careful medical attention. For minimally-invasive thermal ablation procedures, the monitoring of heat distribution is one of the biggest challenges. In this work, three approaches for volumetric heat map reconstruction (Delauney triangulation, minimum volume enclosing ellipsoids (MVEE) and splines) are presented based on uniformly distributed 2D MRI phase images rotated around the applicator's main axis. We compare them with our previous temperature interpolation method with respect to accuracy, robustness and adaptability. All approaches are evaluated during MWA treatment on the same data sets consisting of 13 ex vivo bio protein phantoms, including six phantoms with simulated heat sink effects. Regarding accuracy, the DSC similarity results show a strong trend towards the MVEE ([Formula: see text]) and the splines ([Formula: see text]) method compared to the Delauney triangulation ([Formula: see text]) or the temperature interpolation ([Formula: see text]). Robustness is increased for all three approaches and the adaptability shows a significant trend towards the initial interpolation method and the splines. To overcome local inhomogeneities in the acquired data, the use of adaptive simulations should be considered in the future. In addition, the transfer to in vivo animal experiments should be considered to test for clinical applicability.