1
|
Alpers J, Rötzer M, Gutberlet M, Wacker F, Hensen B, Hansen C. Adaptive simulation of 3D thermometry maps for interventional MR-guided tumor ablation using Pennes' bioheat equation and isotherms. Sci Rep 2022; 12:20356. [PMID: 36437405 PMCID: PMC9701800 DOI: 10.1038/s41598-022-24911-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2022] [Accepted: 11/22/2022] [Indexed: 11/28/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Julian Alpers
- grid.5807.a0000 0001 1018 4307Otto-von-Guericke University, Faculty of Computer Science, Magdeburg, 39106 Germany ,grid.5807.a0000 0001 1018 4307Otto-von-Guericke University, Research Campus STIMULATE, Magdeburg, 39106 Germany
| | - Maximilian Rötzer
- grid.5807.a0000 0001 1018 4307Otto-von-Guericke University, Faculty of Computer Science, Magdeburg, 39106 Germany ,grid.5807.a0000 0001 1018 4307Otto-von-Guericke University, Research Campus STIMULATE, Magdeburg, 39106 Germany
| | - Marcel Gutberlet
- grid.10423.340000 0000 9529 9877Hannover Medical School, Institute for Diagnostic and Interventional Radiology, Hannover, 30625 Germany ,grid.5807.a0000 0001 1018 4307Otto-von-Guericke University, Research Campus STIMULATE, Magdeburg, 39106 Germany
| | - Frank Wacker
- grid.10423.340000 0000 9529 9877Hannover Medical School, Institute for Diagnostic and Interventional Radiology, Hannover, 30625 Germany ,grid.5807.a0000 0001 1018 4307Otto-von-Guericke University, Research Campus STIMULATE, Magdeburg, 39106 Germany
| | - Bennet Hensen
- grid.10423.340000 0000 9529 9877Hannover Medical School, Institute for Diagnostic and Interventional Radiology, Hannover, 30625 Germany ,grid.5807.a0000 0001 1018 4307Otto-von-Guericke University, Research Campus STIMULATE, Magdeburg, 39106 Germany
| | - Christian Hansen
- grid.5807.a0000 0001 1018 4307Otto-von-Guericke University, Faculty of Computer Science, Magdeburg, 39106 Germany ,grid.5807.a0000 0001 1018 4307Otto-von-Guericke University, Research Campus STIMULATE, Magdeburg, 39106 Germany
| |
Collapse
|
2
|
Hübner F, Blauth S, Leithäuser C, Schreiner R, Siedow N, Vogl TJ. Validating a simulation model for laser-induced thermotherapy using MR thermometry. Int J Hyperthermia 2022; 39:1315-1326. [PMID: 36220179 DOI: 10.1080/02656736.2022.2129102] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022] Open
Abstract
OBJECTIVES We want to investigate whether temperature measurements obtained from MR thermometry are accurate and reliable enough to aid the development and validation of simulation models for Laser-induced interstitial thermotherapy (LITT). METHODS Laser-induced interstitial thermotherapy (LITT) is applied to ex-vivo porcine livers. An artificial blood vessel is used to study the cooling effect of large blood vessels in proximity to the ablation zone. The experimental setting is simulated using a model based on partial differential equations (PDEs) for temperature, radiation, and tissue damage. The simulated temperature distributions are compared to temperature data obtained from MR thermometry. RESULTS The overall agreement between measurement and simulation is good for two of our four test cases, while for the remaining cases drift problems with the thermometry data have been an issue. At higher temperatures local deviations between simulation and measurement occur in close proximity to the laser applicator and the vessel. This suggests that certain aspects of the model may need some refinement. CONCLUSION Thermometry data is well-suited for aiding the development of simulations models since it shows where refinements are necessary and enables the validation of such models.
Collapse
Affiliation(s)
- Frank Hübner
- Institute for Diagnostic and Interventional Radiology of the J.W. Goethe University Hospital, Frankfurt am Main, Germany
| | | | | | - Roland Schreiner
- Institute for Diagnostic and Interventional Radiology of the J.W. Goethe University Hospital, Frankfurt am Main, Germany
| | | | - Thomas J Vogl
- Institute for Diagnostic and Interventional Radiology of the J.W. Goethe University Hospital, Frankfurt am Main, Germany
| |
Collapse
|
3
|
Soltani-Sarvestani MA, Cotin S, Saccomandi P. Unscented Kalman Filtering for Real Time Thermometry During Laser Ablation Interventions. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2022; 2022:3485-3488. [PMID: 36085919 DOI: 10.1109/embc48229.2022.9871282] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
We present a data-assimilation Bayesian framework in the context of laser ablation for the treatment of cancer. For solving the nonlinear estimation of the tissue temperature evolving during the therapy, the Unscented Kalman Filter (UKF) predicts the next thermal status and controls the ablation process, based on sparse temperature information. The purpose of this paper is to study the outcome of the prediction model based on UKF and to assess the influence of different model settings on the framework performances. In particular, we analyze the effects of the time resolution of the filter and the number and the location of the observations. Clinical Relevance - The application of a data-assimilation approach based on limited temperature information allows to monitor and predict in real-time the thermal effects induced by thermal therapy for tumors.
Collapse
|
4
|
Schulmann N, Soltani-Sarvestani MA, De Landro M, Korganbayev S, Cotin S, Saccomandi P. Model-Based Thermometry for Laser Ablation Procedure Using Kalman Filters and Sparse Temperature Measurements. IEEE Trans Biomed Eng 2022; 69:2839-2849. [PMID: 35230944 DOI: 10.1109/tbme.2022.3155574] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
In this work, we implement a data-assimilation Bayesian framework for the reconstruction of the spatiotemporal profile of the tissue temperature during laser irradiation. The predictions of a physical model simulating the heat transfer in the tissue are associated with sparse temperature measurements, using an Unscented Kalman Filter. We compare a standard state-estimation filtering procedure with a joint-estimation (state and parameters) approach: whereas in the state-estimation only the temperature is evaluated, in the joint-estimation the filter corrects also uncertain model parameters (i.e., the medium thermal diffusivity, and laser beam properties). We have tested the method on synthetic temperature data, and on the temperature measured on agar-gel phantom and porcine liver with fiber optic sensors. The joint-estimation allows retrieving an accurate estimate of the temperature distribution with a maximal error < 1.5 C in both synthetic and liver 1D data, and < 2 C in phantom 2D data. Our approach allows also suggesting a strategy for optimizing the temperature estimation based on the positions of the sensors. Under the constraint of using only two sensors, optimal temperature estimations are obtained when one sensor is placed in proximity of the source, and the other one is in a non-symmetrical position.
Collapse
|
5
|
Fahrenholtz SJ, Guo C, MacLellan CJ, Yung JP, Hwang KP, Layman RR, Stafford RJ, Cressman E. Temperature mapping of exothermic in situ chemistry: imaging of thermoembolization via MR. Int J Hyperthermia 2020; 36:730-738. [PMID: 31362538 DOI: 10.1080/02656736.2019.1635274] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Purpose: MR temperature imaging (MRTI) was employed for visualizing the spatiotemporal evolution of the exotherm of thermoembolization, an investigative transarterial treatment for solid tumors. Materials and methods: Five explanted kidneys were injected with thermoembolic solutions, and monitored by MRTI. In three nonselective experiments, 5 ml of 4 mol/l dichloroacetyl chloride (DCA-Cl) solution in a hydrocarbon vehicle was injected via the main renal artery. For two of these three, MRTI temperature data were compared to fiber optic thermal probes. Another two kidneys received selective injections, treating only portions of the kidneys with 1 ml of 2 mol/l DCA-Cl. MRTI data were acquired and compared to changes in pre- and post-injection CT. Specimens were bisected and photographed for gross pathology 24 h post-procedure. Results: MRTI temperature estimates were within ±1 °C of the probes. In experiments without probes, MRTI measured increases of 30 °C. Some regions had not reached peak temperature by the end of the >18 min acquisition. MRTI indicated the initial heating occurred in the renal cortex, gradually spreading more proximally toward the main renal artery. Gross pathology showed the nonselective injection denatured the entire kidney whereas in the selective injections, only the treated territory was coagulated. Conclusion: The spatiotemporal evolution of thermoembolization was visualized for the first time using noninvasive MRTI, providing unique insight into the thermodynamics of thermoembolization. Précis Thermoembolization is being investigated as a novel transarterial treatment. In order to begin to characterize delivery of this novel treatment modality and aid translation from the laboratory to patients, we employ MR temperature imaging to visualize the spatiotemporal distribution of temperature from thermoembolization in ex vivo tissue.
Collapse
Affiliation(s)
- Samuel John Fahrenholtz
- a Department of Imaging Physics, The University of Texas MD Anderson Cancer Center , Houston , TX , USA
| | - Chunxiao Guo
- b Department of Interventional Radiology, The University of Texas MD Anderson Cancer Center , Houston , TX , USA
| | - Christopher J MacLellan
- a Department of Imaging Physics, The University of Texas MD Anderson Cancer Center , Houston , TX , USA
| | - Joshua P Yung
- a Department of Imaging Physics, The University of Texas MD Anderson Cancer Center , Houston , TX , USA
| | - Ken-Pin Hwang
- a Department of Imaging Physics, The University of Texas MD Anderson Cancer Center , Houston , TX , USA
| | - Rick R Layman
- a Department of Imaging Physics, The University of Texas MD Anderson Cancer Center , Houston , TX , USA
| | - R Jason Stafford
- a Department of Imaging Physics, The University of Texas MD Anderson Cancer Center , Houston , TX , USA
| | - Erik Cressman
- b Department of Interventional Radiology, The University of Texas MD Anderson Cancer Center , Houston , TX , USA
| |
Collapse
|
6
|
Odéen H, Parker DL. Magnetic resonance thermometry and its biological applications - Physical principles and practical considerations. PROGRESS IN NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY 2019; 110:34-61. [PMID: 30803693 PMCID: PMC6662927 DOI: 10.1016/j.pnmrs.2019.01.003] [Citation(s) in RCA: 70] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/11/2018] [Accepted: 01/23/2019] [Indexed: 05/25/2023]
Abstract
Most parameters that influence the magnetic resonance imaging (MRI) signal experience a temperature dependence. The fact that MRI can be used for non-invasive measurements of temperature and temperature change deep inside the human body has been known for over 30 years. Today, MR temperature imaging is widely used to monitor and evaluate thermal therapies such as radio frequency, microwave, laser, and focused ultrasound therapy. In this paper we cover the physical principles underlying the biological applications of MR temperature imaging and discuss practical considerations and remaining challenges. For biological tissue, the MR signal of interest comes mostly from hydrogen protons of water molecules but also from protons in, e.g., adipose tissue and various metabolites. Most of the discussed methods, such as those using the proton resonance frequency (PRF) shift, T1, T2, and diffusion only measure temperature change, but measurements of absolute temperatures are also possible using spectroscopic imaging methods (taking advantage of various metabolite signals as internal references) or various types of contrast agents. Currently, the PRF method is the most used clinically due to good sensitivity, excellent linearity with temperature, and because it is largely independent of tissue type. Because the PRF method does not work in adipose tissues, T1- and T2-based methods have recently gained interest for monitoring temperature change in areas with high fat content such as the breast and abdomen. Absolute temperature measurement methods using spectroscopic imaging and contrast agents often offer too low spatial and temporal resolution for accurate monitoring of ablative thermal procedures, but have shown great promise in monitoring the slower and usually less spatially localized temperature change observed during hyperthermia procedures. Much of the current research effort for ablative procedures is aimed at providing faster measurements, larger field-of-view coverage, simultaneous monitoring in aqueous and adipose tissues, and more motion-insensitive acquisitions for better precision measurements in organs such as the heart, liver, and kidneys. For hyperthermia applications, larger coverage, motion insensitivity, and simultaneous aqueous and adipose monitoring are also important, but great effort is also aimed at solving the problem of long-term field drift which gets interpreted as temperature change when using the PRF method.
Collapse
Affiliation(s)
- Henrik Odéen
- University of Utah, Utah Center for Advanced Imaging Research, Department of Radiology and Imaging Sciences, 729 Arapeen Drive, Salt Lake City, UT 84108-1217, USA.
| | - Dennis L Parker
- University of Utah, Utah Center for Advanced Imaging Research, Department of Radiology and Imaging Sciences, 729 Arapeen Drive, Salt Lake City, UT 84108-1217, USA.
| |
Collapse
|
7
|
Odéen H, Parker DL. Improved MR thermometry for laser interstitial thermotherapy. Lasers Surg Med 2019; 51:286-300. [PMID: 30645017 DOI: 10.1002/lsm.23049] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/28/2018] [Indexed: 12/24/2022]
Abstract
OBJECTIVES To develop, test and evaluate improved 2D and 3D protocols for proton resonance frequency shift magnetic resonance temperature imaging (MRTI) of laser interstitial thermal therapy (LITT). The objective was to develop improved MRTI protocols in terms of temperature measurement precision and volume coverage compared to the 2D MRTI protocol currently used with a commercially available LITT system. METHODS Four different 2D protocols and four different 3D protocols were investigated. The 2D protocols used multi-echo readouts to prolong the total MR sampling time and hence the MRTI precision, without prolonging the total acquisition time. The 3D protocols provided volumetric thermometry by acquiring a slab of 12 contiguous slices in the same acquisition time as the 2D protocols. The study only considered readily available pulse sequences (Cartesian 2D and 3D gradient recalled echo and echo planar imaging [EPI]) and methods (partial Fourier and parallel imaging) to ensure wide availability and rapid clinical implementation across vendors and field strengths. In vivo volunteer studies were performed to investigate and compare MRTI precision and image quality. Phantom experiments with LITT heating were performed to investigate and compare MRTI precision and accuracy. Different coil setups were used in the in vivo studies to assess precision differences between using local (such as flex and head coils) and non-local (i.e., body coil) receive coils. Studies were performed at both 1.5 T and 3 T. RESULTS The improved 2D protocols provide up to a factor of two improvement in the MRTI precision in the same acquisition time, compared to the currently used clinical protocol. The 3D echo planar imaging protocols provide comparable precision as the currently used 2D clinical protocol, but over a substantially larger field of view, without increasing the acquisition time. As expected, local receive coils perform substantially better than the body coil, and 3 T provides better MRTI accuracy and precision than 1.5 T. 3D data can be zero-filled interpolated in all three dimensions (as opposed to just two dimensions for 2D data), reducing partial volume effects and measuring higher maximum temperature rises. CONCLUSIONS With the presented protocols substantially improved MRTI precision (for 2D imaging) or greatly improved field of view coverage (for 3D imaging) can be achieved in the same acquisition time as the currently used protocol. Only widely available pulse sequences and acquisition methods were investigated, which should ensure quick translation to the clinic. Lasers Surg. Med. 51:286-300, 2019. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Henrik Odéen
- Utah Center for Advanced Imaging Research, Department of Radiology and Imaging Sciences, University of Utah, Salt Lake City, Utah
| | - Dennis L Parker
- Utah Center for Advanced Imaging Research, Department of Radiology and Imaging Sciences, University of Utah, Salt Lake City, Utah
| |
Collapse
|
8
|
Yung JP, Fuentes D, MacLellan CJ, Maier F, Liapis Y, Hazle JD, Stafford RJ. Referenceless magnetic resonance temperature imaging using Gaussian process modeling. Med Phys 2017; 44:3545-3555. [PMID: 28317125 DOI: 10.1002/mp.12231] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2016] [Revised: 12/15/2016] [Accepted: 01/09/2017] [Indexed: 11/12/2022] Open
Abstract
PURPOSE During magnetic resonance (MR)-guided thermal therapies, water proton resonance frequency shift (PRFS)-based MR temperature imaging can quantitatively monitor tissue temperature changes. It is widely known that the PRFS technique is easily perturbed by tissue motion, tissue susceptibility changes, magnetic field drift, and modality-dependent applicator-induced artifacts. Here, a referenceless Gaussian process modeling (GPM)-based estimation of the PRFS is investigated as a methodology to mitigate unwanted background field changes. The GPM offers a complementary trade-off between data fitting and smoothing and allows prior information to be used. The end result being the GPM provides a full probabilistic prediction and an estimate of the uncertainty. METHODS GPM was employed to estimate the covariance between the spatial position and MR phase measurements. The mean and variance provided by the statistical model extrapolated background phase values from nonheated neighboring voxels used to train the model. MR phase predictions in the heating ROI are computed using the spatial coordinates as the test input. The method is demonstrated in ex vivo rabbit liver tissue during focused ultrasound heating with manually introduced perturbations (n = 6) and in vivo during laser-induced interstitial thermal therapy to treat the human brain (n = 1) and liver (n = 1). RESULTS Temperature maps estimated using the GPM referenceless method demonstrated a RMS error of <0.8°C with artifact-induced reference-based MR thermometry during ex vivo heating using focused ultrasound. Nonheated surrounding areas were <0.5°C from the artifact-free MR measurements. The GPM referenceless MR temperature values and thermally damaged regions were within the 95% confidence interval during in vivo laser ablations. CONCLUSIONS A new approach to estimation for referenceless PRFS temperature imaging is introduced that allows for an accurate probabilistic extrapolation of the background phase. The technique demonstrated reliable temperature estimates in the presence of the background phase changes and was demonstrated useful in the in vivo brain and liver ablation scenarios presented.
Collapse
Affiliation(s)
- Joshua P Yung
- Unit 1902, Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Blvd., Houston, TX, 77030, USA.,The University of Texas Graduate School of Biomedical Sciences at Houston, 6767 Bertner Ave., Houston, TX, 77030, USA
| | - David Fuentes
- Unit 1902, Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Blvd., Houston, TX, 77030, USA.,The University of Texas Graduate School of Biomedical Sciences at Houston, 6767 Bertner Ave., Houston, TX, 77030, USA
| | - Christopher J MacLellan
- Unit 1902, Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Blvd., Houston, TX, 77030, USA.,The University of Texas Graduate School of Biomedical Sciences at Houston, 6767 Bertner Ave., Houston, TX, 77030, USA
| | - Florian Maier
- Unit 1902, Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Blvd., Houston, TX, 77030, USA
| | - Yannis Liapis
- Unit 1902, Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Blvd., Houston, TX, 77030, USA
| | - John D Hazle
- Unit 1902, Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Blvd., Houston, TX, 77030, USA.,The University of Texas Graduate School of Biomedical Sciences at Houston, 6767 Bertner Ave., Houston, TX, 77030, USA
| | - R Jason Stafford
- Unit 1902, Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Blvd., Houston, TX, 77030, USA.,The University of Texas Graduate School of Biomedical Sciences at Houston, 6767 Bertner Ave., Houston, TX, 77030, USA
| |
Collapse
|
9
|
Madankan R, Stefan W, Fahrenholtz SJ, MacLellan CJ, Hazle JD, Stafford RJ, Weinberg JS, Rao G, Fuentes D. Accelerated magnetic resonance thermometry in the presence of uncertainties. Phys Med Biol 2017; 62:214-245. [PMID: 27991449 DOI: 10.1088/1361-6560/62/1/214] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
A model-based information theoretic approach is presented to perform the task of magnetic resonance (MR) thermal image reconstruction from a limited number of observed samples on k-space. The key idea of the proposed approach is to optimally detect samples of k-space that are information-rich with respect to a model of the thermal data acquisition. These highly informative k-space samples can then be used to refine the mathematical model and efficiently reconstruct the image. The information theoretic reconstruction was demonstrated retrospectively in data acquired during MR-guided laser induced thermal therapy (MRgLITT) procedures. The approach demonstrates that locations with high-information content with respect to a model-based reconstruction of MR thermometry may be quantitatively identified. These information-rich k-space locations are demonstrated to be useful as a guide for k-space undersampling techniques. The effect of interactively increasing the predicted number of data points used in the subsampled model-based reconstruction was quantified using the L2-norm of the distance between the subsampled and fully sampled reconstruction. Performance of the proposed approach was also compared with uniform rectilinear subsampling and variable-density Poisson disk subsampling techniques. The proposed subsampling scheme resulted in accurate reconstructions using a small fraction of k-space points, suggesting that the reconstruction technique may be useful in improving the efficiency of thermometry data temporal resolution.
Collapse
Affiliation(s)
- R Madankan
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas, USA
| | | | | | | | | | | | | | | | | |
Collapse
|
10
|
Zhang Y, Chen S, Deng K, Chen B, Wei X, Yang J, Wang S, Ying K. Kalman Filtered Bio Heat Transfer Model Based Self-adaptive Hybrid Magnetic Resonance Thermometry. IEEE TRANSACTIONS ON MEDICAL IMAGING 2017; 36:194-202. [PMID: 27552745 DOI: 10.1109/tmi.2016.2601440] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
To develop a self-adaptive and fast thermometry method by combining the original hybrid magnetic resonance thermometry method and the bio heat transfer equation (BHTE) model. The proposed Kalman filtered Bio Heat Transfer Model Based Self-adaptive Hybrid Magnetic Resonance Thermometry, abbreviated as KalBHT hybrid method, introduced the BHTE model to synthesize a window on the regularization term of the hybrid algorithm, which leads to a self-adaptive regularization both spatially and temporally with change of temperature. Further, to decrease the sensitivity to accuracy of the BHTE model, Kalman filter is utilized to update the window at each iteration time. To investigate the effect of the proposed model, computer heating simulation, phantom microwave heating experiment and dynamic in-vivo model validation of liver and thoracic tumor were conducted in this study. The heating simulation indicates that the KalBHT hybrid algorithm achieves more accurate results without adjusting λ to a proper value in comparison to the hybrid algorithm. The results of the phantom heating experiment illustrate that the proposed model is able to follow temperature changes in the presence of motion and the temperature estimated also shows less noise in the background and surrounding the hot spot. The dynamic in-vivo model validation with heating simulation demonstrates that the proposed model has a higher convergence rate, more robustness to susceptibility problem surrounding the hot spot and more accuracy of temperature estimation. In the healthy liver experiment with heating simulation, the RMSE of the hot spot of the proposed model is reduced to about 50% compared to the RMSE of the original hybrid model and the convergence time becomes only about one fifth of the hybrid model. The proposed model is able to improve the accuracy of the original hybrid algorithm and accelerate the convergence rate of MR temperature estimation.
Collapse
|
11
|
Odéen H, Almquist S, de Bever J, Christensen DA, Parker DL. MR thermometry for focused ultrasound monitoring utilizing model predictive filtering and ultrasound beam modeling. J Ther Ultrasound 2016; 4:23. [PMID: 27688881 PMCID: PMC5032243 DOI: 10.1186/s40349-016-0067-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2016] [Accepted: 09/02/2016] [Indexed: 12/28/2022] Open
Abstract
Background A major challenge in using magnetic resonance temperature imaging (MRTI) to monitor focused ultrasound (FUS) applications is achieving high spatio-temporal resolution over a large field of view (FOV). This is important to accurately monitor all ultrasound (US) power depositions. Magnetic resonance (MR) subsampling in conjunction with thermal model-based reconstruction of the MRTI utilizing Pennes bioheat transfer equation (PBTE) is one promising approach. The thermal properties used in the thermal model are often estimated from a pre-treatment, low-power sonication. Methods In this proof-of-concept study we investigate the use of US simulations computed using the hybrid angular spectrum (HAS) method to estimate the US power deposition density Q, thereby avoiding the pre-treatment sonication and any potential tissue damage. MRTI reconstructions are performed using a thermal model-based reconstruction method called model predictive filtering (MPF). Experiments are performed in a homogeneous gelatin phantom and in a gelatin phantom with embedded plastic skull. MPF reconstructions are compared to separate sonications imaged with fully sampled data over a smaller FOV. Temperature root-mean-square errors (RMSE) and focal spot positions and shapes are evaluated. Results HAS simulations accurately predict the location of the focal spot (to within 1 mm) in both phantoms. Accurate temperature maps (RMSE below 1 °C), where the location of the focal spot agrees well with fully sampled “truth” (to within 1 mm), are also achieved in both phantoms. Conclusions HAS simulations can be used to accurately predict the focal spot location in homogeneous media and when focusing through an aberrating plastic skull. The HAS simulated power deposition (Q) patterns can be used in the MPF thermal model-based reconstruction to obtain accurate temperature maps with high spatio-temporal resolution over large FOVs.
Collapse
Affiliation(s)
- Henrik Odéen
- Utah Center for Advanced Imaging Research, Department of Radiology, University of Utah, Salt Lake City, UT USA
| | - Scott Almquist
- School of Computing, University of Utah, Salt Lake City, UT USA
| | - Joshua de Bever
- Utah Center for Advanced Imaging Research, Department of Radiology, University of Utah, Salt Lake City, UT USA
| | - Douglas A Christensen
- Department of Bioengineering, University of Utah, Salt Lake City, UT USA ; Department of Electrical and Computer Engineering, University of Utah, Salt Lake City, UT USA
| | - Dennis L Parker
- Utah Center for Advanced Imaging Research, Department of Radiology, University of Utah, Salt Lake City, UT USA
| |
Collapse
|
12
|
Odéen H, Todd N, Dillon C, Payne A, Parker DL. Model predictive filtering MR thermometry: Effects of model inaccuracies, k-space reduction factor, and temperature increase rate. Magn Reson Med 2015; 75:207-16. [PMID: 25726934 DOI: 10.1002/mrm.25622] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2014] [Revised: 12/24/2014] [Accepted: 12/29/2014] [Indexed: 01/20/2023]
Abstract
PURPOSE Evaluate effects of model parameter inaccuracies (thermal conductivity, k, and ultrasound power deposition density, Q), k-space reduction factor (R), and rate of temperature increase ( T˙) in a thermal model-based reconstruction for MR-thermometry during focused-ultrasound heating. METHODS Simulations and ex vivo experiments were performed to investigate the accuracy of the thermal model and the model predictive filtering (MPF) algorithm for varying R and T˙, and their sensitivity to errors in k and Q. Ex vivo data was acquired with a segmented EPI pulse sequence to achieve large field-of-view (192 × 162 × 96 mm) four-dimensional temperature maps with high spatiotemporal resolution (1.5 × 1.5 × 2.0 mm, 1.7 s). RESULTS In the simulations, 50% errors in k and Q resulted in maximum temperature root mean square errors (RMSE) of 6 °C for model only and 3 °C for MPF. Using recently developed methods, estimates of k and Q were accurate to within 3%. The RMSE between MPF and true temperature increased with R and T˙. In the ex vivo study the RMSE remained below 0.7 °C for R ranging from 4 to 12 and T˙ of 0.28-0.75 °C/s. CONCLUSION Errors in MPF temperatures occur due to errors in k and Q. These MPF temperature errors increase with increase in R and T˙, but are smaller than those obtained using the thermal model alone.
Collapse
Affiliation(s)
- Henrik Odéen
- Department of Radiology, University of Utah, Salt Lake City, Utah, USA.,Department of Physics and Astronomy, University of Utah, Salt Lake City, Utah, USA
| | - Nick Todd
- Department of Radiology, University of Utah, Salt Lake City, Utah, USA
| | - Christopher Dillon
- Department of Bioengineering, University of Utah, Salt Lake City, Utah, USA
| | - Allison Payne
- Department of Radiology, University of Utah, Salt Lake City, Utah, USA
| | - Dennis L Parker
- Department of Radiology, University of Utah, Salt Lake City, Utah, USA
| |
Collapse
|
13
|
Fahrenholtz SJ, Stafford RJ, Maier F, Hazle JD, Fuentes D. Generalised polynomial chaos-based uncertainty quantification for planning MRgLITT procedures. Int J Hyperthermia 2013; 29:324-35. [PMID: 23692295 DOI: 10.3109/02656736.2013.798036] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
PURPOSE A generalised polynomial chaos (gPC) method is used to incorporate constitutive parameter uncertainties within the Pennes representation of bioheat transfer phenomena. The stochastic temperature predictions of the mathematical model are critically evaluated against MR thermometry data for planning MR-guided laser-induced thermal therapies (MRgLITT). METHODS The Pennes bioheat transfer model coupled with a diffusion theory approximation of laser tissue interaction was implemented as the underlying deterministic kernel. A probabilistic sensitivity study was used to identify parameters that provide the most variance in temperature output. Confidence intervals of the temperature predictions are compared to MR temperature imaging (MRTI) obtained during phantom and in vivo canine (n = 4) MRgLITT experiments. The gPC predictions were quantitatively compared to MRTI data using probabilistic linear and temporal profiles as well as 2-D 60 °C isotherms. RESULTS Optical parameters provided the highest variance in the model output (peak standard deviation: anisotropy 3.51 °C, absorption 2.94 °C, scattering 1.84 °C, conductivity 1.43 °C, and perfusion 0.94 °C). Further, within the statistical sense considered, a non-linear model of the temperature and damage-dependent perfusion, absorption, and scattering is captured within the confidence intervals of the linear gPC method. Multivariate stochastic model predictions using parameters with the dominant sensitivities show good agreement with experimental MRTI data. CONCLUSIONS Given parameter uncertainties and mathematical modelling approximations of the Pennes bioheat model, the statistical framework demonstrates conservative estimates of the therapeutic heating and has potential for use as a computational prediction tool for thermal therapy planning.
Collapse
Affiliation(s)
- Samuel J Fahrenholtz
- Department of Imaging Physics, The University of Texas M.D. Anderson Cancer Center, Houston, TX 77054, USA
| | | | | | | | | |
Collapse
|
14
|
Fuentes D, Elliott A, Weinberg JS, Shetty A, Hazle JD, Stafford RJ. An inverse problem approach to recovery of in vivo nanoparticle concentrations from thermal image monitoring of MR-guided laser induced thermal therapy. Ann Biomed Eng 2012; 41:100-11. [PMID: 22918665 DOI: 10.1007/s10439-012-0638-9] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2012] [Accepted: 08/04/2012] [Indexed: 12/15/2022]
Abstract
Quantification of local variations in the optical properties of tumor tissue introduced by the presence of gold-silica nanoparticles (NP) presents significant opportunities in monitoring and control of NP-mediated laser induced thermal therapy (LITT) procedures. Finite element methods of inverse parameter recovery constrained by a Pennes bioheat transfer model were applied to estimate the optical parameters. Magnetic resonance temperature imaging (MRTI) acquired during a NP-mediated LITT of a canine transmissible venereal tumor in brain was used in the presented statistical inverse problem formulation. The maximum likelihood (ML) value of the optical parameters illustrated a marked change in the periphery of the tumor corresponding with the expected location of NP and area of selective heating observed on MRTI. Parameter recovery information became increasingly difficult to infer in distal regions of tissue where photon fluence had been significantly attenuated. Finite element temperature predictions using the ML parameter values obtained from the solution of the inverse problem are able to reproduce the NP selective heating within 5 °C of measured MRTI estimations along selected temperature profiles. Results indicate the ML solution found is able to sufficiently reproduce the selectivity of the NP mediated laser induced heating and therefore the ML solution is likely to return useful optical parameters within the region of significant laser fluence.
Collapse
Affiliation(s)
- D Fuentes
- Department of Imaging Physics, The University of Texas M.D. Anderson Cancer Center, Houston, TX 77030, USA.
| | | | | | | | | | | |
Collapse
|
15
|
Navab N, Taylor R, Yang GZ. Guest editorial: special issue on interventional imaging. IEEE TRANSACTIONS ON MEDICAL IMAGING 2012; 31:857-859. [PMID: 22582415 DOI: 10.1109/tmi.2012.2189153] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
|