1
|
Amini S, Ahmadikia H. Treatment Modelling of a 3D Tumour in Brain by Laser-Induced Interstitial Thermotherapy. J Lasers Med Sci 2023; 14:e26. [PMID: 37744014 PMCID: PMC10517584 DOI: 10.34172/jlms.2023.26] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Accepted: 03/13/2023] [Indexed: 09/26/2023]
Abstract
Introduction: There are some ways to examine heat transfer in tumor tissue, which is an important issue in bioengineering. One of these ways uses the bioheat equation, proposed by Pennes, in a continuous medium. Another one uses a porous medium to model heat transfer in living tissues. The objective of this paper was to study an approach to modelling the temperature distribution and tumour ablation in brain tissue and compare results to Pennes' approach. Methods: This approach presents and uses a porous medium as the tissue instead of a continuous medium. In addition, the two approaches (simulation in continuous and porous medium) are compared in terms of temperature simulation and amount of cell ablation. The density, heat conduction factor, and blood perfusion rate are considered functions of temperature. Results: In these approaches, after an 85-second treatment, the temperature increases to about 90°C. The temperature increase of the porous medium is relatively the same as that of the continuous medium and for this reason, the amount of cancerous cells that are ablated in a porous medium is approximately the same as that in a continuous medium. The volume of cell ablation is about 6500 mm3 for two ideas. In addition, the degree of damage, computed from the Arrhenius integral method, and the ablated volume of the tumour endorse equality at the end of treatment. According to the results, similar to the continuous approach, the porous approach predicts the temperature and amount of volume of damaged cells. Conclusion: Therefore, it is possible to use the porous approach instead of the Pennes approach for tumour treatment.
Collapse
Affiliation(s)
- Sadegh Amini
- Mechanical Engineering Department, University of Isfahan, Isfahan, Iran
| | - Hossein Ahmadikia
- Mechanical Engineering Department, University of Isfahan, Isfahan, Iran
| |
Collapse
|
2
|
Faridi P, Keselman P, Fallahi H, Prakash P. Experimental assessment of microwave ablation computational modeling with MR thermometry. Med Phys 2020; 47:3777-3788. [PMID: 32506550 DOI: 10.1002/mp.14318] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Revised: 05/22/2020] [Accepted: 05/24/2020] [Indexed: 12/23/2022] Open
Abstract
PURPOSE Computational models are widely used during the design and characterization of microwave ablation (MWA) devices, and have been proposed for pretreatment planning. Our objective was to assess three-dimensional (3D) transient temperature and ablation profiles predicted by MWA computational models with temperature profiles measured experimentally using magnetic resonance (MR) thermometry in ex vivo bovine liver. MATERIALS AND METHODS We performed MWA in ex vivo tissue under MR guidance using a custom, 2.45 GHz water-cooled applicator. MR thermometry data were acquired for 2 min prior to heating, during 5-10 min microwave exposures, and for 3 min following heating. Fiber-optic temperature sensors were used to validate the accuracy of MR temperature measurements. A total of 13 ablation experiments were conducted using 30-50 W applied power at the applicator input. MWA computational models were implemented using the finite element method, and incorporated temperature-dependent changes in tissue physical properties. Model-predicted ablation zone extents were compared against MRI-derived Arrhenius thermal damage maps using the Dice similarity coefficient (DSC). RESULTS Prior to heating, the observed standard deviation of MR temperature data was in the range of 0.3-0.7°C. Mean absolute error between MR temperature measurements and fiber-optic temperature probes during heating was in the range of 0.5-2.8°C. The mean DSC between model-predicted ablation zones and MRI-derived Arrhenius thermal damage maps for 13 experimental set-ups was 0.95. When comparing simulated and experimentally (i.e. using MRI) measured temperatures, the mean absolute error (MAE %) relative to maximum temperature change was in the range 5%-8.5%. CONCLUSION We developed a system for characterizing 3D transient temperature and ablation profiles with MR thermometry during MWA in ex vivo liver tissue, and applied the system for experimental validation of MWA computational models.
Collapse
Affiliation(s)
- Pegah Faridi
- Mike Wiegers Department of Electrical and Computer Engineering, Kansas State University, Manhattan, KS, 66506, USA
| | - Paul Keselman
- Hoglund Brain Imaging Center, University of Kansas Medical Center, Kansas City, KS, 66160, USA
| | - Hojjatollah Fallahi
- Mike Wiegers Department of Electrical and Computer Engineering, Kansas State University, Manhattan, KS, 66506, USA
| | - Punit Prakash
- Mike Wiegers Department of Electrical and Computer Engineering, Kansas State University, Manhattan, KS, 66506, USA
| |
Collapse
|
3
|
Fuentes D, Fahrenholtz SJ, Guo C, MacLellan CJ, Layman RR, Rivière B, Stafford RJ, Cressman E. Mathematical modeling of mass and energy transport for thermoembolization. Int J Hyperthermia 2020; 37:356-365. [PMID: 32308071 PMCID: PMC10558277 DOI: 10.1080/02656736.2020.1749317] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2019] [Revised: 03/16/2020] [Accepted: 03/23/2020] [Indexed: 12/27/2022] Open
Abstract
Background: Thermoembolization presents a unique treatment alternative for patients diagnosed with hepatocellular carcinoma. The approach delivers a reagent that undergoes an exothermic chemical reaction and combines the benefits of embolic as well as thermal- and chemical-ablative therapy modalities. The target tissue and vascular bed are subjected to simultaneous hyperthermia, ischemia, and chemical denaturation in a single procedure. To guide optimal delivery, we developed a mathematical model for understanding the competing diffusive and convective effects observed in thermoembolization delivery protocols.Methods: A mixture theory formulation was used to mathematically model thermoembolization as chemically reacting transport of an electrophile, dichloroacetyl chloride (DCACl), within porous living tissue. Mass and energy transport of each relevant constituent are considered. Specifically, DCACl is injected into the vessels and exothermically reacts with water in the blood or tissue to form dichloroacetic acid and hydrochloric acid. Neutralization reactions are assumed instantaneous in this approach. We validated the mathematical model predictions of temperature using MR thermometry of the thermoembolization procedure performed in ex vivo kidney.Results: Mathematical modeling predictions of tissue death were highly dependent on the vascular geometry, injection pressure, and intrinsic amount of exothermic energy released from the chemical species, and were able to recapitulate the temperature distributions observed in MR thermometry.Conclusion: These efforts present a first step toward formalizing a mathematical model for thermoembolization and are promising for providing insight for delivery protocol optimization. While our approach captured the observed experimental temperature measurements, larger-scale experimental validation is needed to prioritize additional model complexity and fidelity.
Collapse
Affiliation(s)
- David Fuentes
- Department of Imaging Physics, Anderson Cancer Center, The University of Texas MD, Houston, TX, USA
- Department of Computational and Applied Mathematics, Rice University, Houston, TX, USA
| | - Samuel J. Fahrenholtz
- Department of Imaging Physics, Anderson Cancer Center, The University of Texas MD, Houston, TX, USA
| | - Chunxiao Guo
- Department of Interventional Radiology, Anderson Cancer Center, The University of Texas MD, Houston, TX, USA
| | - Christopher J. MacLellan
- Department of Imaging Physics, Anderson Cancer Center, The University of Texas MD, Houston, TX, USA
| | - Rick R. Layman
- Department of Imaging Physics, Anderson Cancer Center, The University of Texas MD, Houston, TX, USA
| | - Beatrice Rivière
- Department of Computational and Applied Mathematics, Rice University, Houston, TX, USA
| | - R. Jason Stafford
- Department of Imaging Physics, Anderson Cancer Center, The University of Texas MD, Houston, TX, USA
| | - Erik Cressman
- Department of Interventional Radiology, Anderson Cancer Center, The University of Texas MD, Houston, TX, USA
| |
Collapse
|
4
|
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
|
5
|
Peng Y, Zou C, Qiao Y, Tie C, Wan Q, Jiang R, Cheng C, Liang D, Zheng H, Li F, Liu X. Fast MR thermometry using an echo-shifted sequence with simultaneous multi-slice imaging. MAGNETIC RESONANCE MATERIALS IN PHYSICS BIOLOGY AND MEDICINE 2018; 31:771-779. [PMID: 29948236 DOI: 10.1007/s10334-018-0692-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/03/2018] [Revised: 05/15/2018] [Accepted: 06/04/2018] [Indexed: 12/20/2022]
Abstract
PURPOSE Real-time monitoring is important for the safety and effectiveness of high-intensity focused ultrasound (HIFU) therapy. Magnetic resonance imaging is the preferred imaging modality for HIFU monitoring, with its unique capability of temperature imaging. For real-time temperature imaging, higher temporal resolution and larger spatial coverage are needed. In this study, a sequence based on the echo-shifted RF-spoiled gradient echo (GRE) with simultaneous multi-slice (SMS) imaging was designed for fast temperature imaging. METHODS A phantom experiment was conducted to evaluate the accuracy of the echo-shifted sequence using a fluorescent fiber thermometer as reference. The temperature uncertainty of the echo-shifted sequence was compared with the traditional GRE sequence at room temperature through the ex vivo porcine muscle. Finally, the ex vivo porcine liver tissue experiment using HIFU heating was performed to demonstrate that the spatial coverage was increased without decreasing temporal resolution. RESULTS The echo-shifted sequence had a better temperature uncertainty performance compared with the traditional GRE sequence with the same temporal resolution. The ex vivo heating experiment confirmed that by combining the SMS technique and echo-shifted sequence, the spatial coverage was increased without decreasing the temporal resolution while maintaining high temperature measurement precision. CONCLUSION The proposed technique was validated as an effective real-time method for monitoring HIFU therapy.
Collapse
Affiliation(s)
- Yuhong Peng
- State Key Laboratory of Ultrasound Engineering in Medicine, Chongqing Key Laboratory of Biomedical Engineering, Co-Founded by Chongqing and the Ministry of Science and Technology, College of Biomedical Engineering, Chongqing Medical University, 153 Box, 1 Yixueyuan Road, Yuzhong District, Chongqing, 400016, China.,Paul C Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, 1068 Xueyuan Blvd, Nanshan District, Shenzhen, 518055, Guangdong, China.,Chongqing Collaborative Innovation Center for Minimally-invasive and Noninvasive Medicine, Chongqing, 400016, China
| | - Chao Zou
- Paul C Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, 1068 Xueyuan Blvd, Nanshan District, Shenzhen, 518055, Guangdong, China
| | - Yangzi Qiao
- Paul C Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, 1068 Xueyuan Blvd, Nanshan District, Shenzhen, 518055, Guangdong, China
| | - Changjun Tie
- Paul C Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, 1068 Xueyuan Blvd, Nanshan District, Shenzhen, 518055, Guangdong, China
| | - Qian Wan
- Paul C Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, 1068 Xueyuan Blvd, Nanshan District, Shenzhen, 518055, Guangdong, China
| | - Rui Jiang
- Paul C Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, 1068 Xueyuan Blvd, Nanshan District, Shenzhen, 518055, Guangdong, China.,Shenzhen College of Advanced Technology, University of Chinese Academy of Sciences, Shenzhen, Guangdong, 518055, China
| | - Chuanli Cheng
- Paul C Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, 1068 Xueyuan Blvd, Nanshan District, Shenzhen, 518055, Guangdong, China.,Shenzhen College of Advanced Technology, University of Chinese Academy of Sciences, Shenzhen, Guangdong, 518055, China
| | - Dong Liang
- Paul C Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, 1068 Xueyuan Blvd, Nanshan District, Shenzhen, 518055, Guangdong, China
| | - Hairong Zheng
- Paul C Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, 1068 Xueyuan Blvd, Nanshan District, Shenzhen, 518055, Guangdong, China
| | - Faqi Li
- State Key Laboratory of Ultrasound Engineering in Medicine, Chongqing Key Laboratory of Biomedical Engineering, Co-Founded by Chongqing and the Ministry of Science and Technology, College of Biomedical Engineering, Chongqing Medical University, 153 Box, 1 Yixueyuan Road, Yuzhong District, Chongqing, 400016, China. .,Chongqing Collaborative Innovation Center for Minimally-invasive and Noninvasive Medicine, Chongqing, 400016, China.
| | - Xin Liu
- Paul C Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, 1068 Xueyuan Blvd, Nanshan District, Shenzhen, 518055, Guangdong, China. .,Chongqing Collaborative Innovation Center for Minimally-invasive and Noninvasive Medicine, Chongqing, 400016, China.
| |
Collapse
|
6
|
Ebbini ES, Simon C, Liu D. Real-time Ultrasound Thermography and Thermometry. IEEE SIGNAL PROCESSING MAGAZINE 2018; 35:166-174. [PMID: 30283214 PMCID: PMC6167021 DOI: 10.1109/msp.2017.2773338] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
|
7
|
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
|
8
|
Zhou M, Melancon M, Stafford RJ, Li J, Nick AM, Tian M, Sood AK, Li C. Precision Nanomedicine Using Dual PET and MR Temperature Imaging-Guided Photothermal Therapy. J Nucl Med 2016; 57:1778-1783. [PMID: 27283932 DOI: 10.2967/jnumed.116.172775] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2016] [Accepted: 04/26/2016] [Indexed: 12/28/2022] Open
Abstract
Imaging-based techniques have enabled the direct integration of noninvasive imaging with minimally invasive interventions such as photothermal therapy (PTT) to improve the precision of treatment. METHODS We investigated the feasibility of PTT for ovarian cancer under the guidance of PET and MR temperature imaging using copper sulfide nanoparticles (CuS NPs). The tumor distribution of the CuS NPs after systemic administration was assessed using highly sensitive, quantifiable PET imaging. Two wavelengths of near-infrared (NIR) lasers-808 and 980 nm-were tested for PTT using noninvasive MR temperature imaging real-time monitoring. RESULTS The in vivo studies revealed that the 980-nm NIR laser had better photothermal effects than the 808-nm NIR laser. These results were in accord with the histologic findings. In vivo PTT using CuS NPs combined with 980-nm laser irradiation achieved significant tumor ablation compared with no treatment control in both subcutaneous (P = 0.007) and orthotopic (P < 0.001) models of ovarian cancer with regard to the percentage of necrotic damage. CONCLUSION Our results indicate that real-time monitoring of the accuracy of PTT is a promising approach for future clinical translation of this emerging thermal ablation technique.
Collapse
Affiliation(s)
- Min Zhou
- Department of Cancer Systems Imaging, The University of Texas MD Anderson Cancer Center, Houston, Texas.,Department of Nuclear Medicine and PET Center, The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China.,Institute of Translational Medicine, School of Medicine, Zhejiang University, Hangzhou, China
| | - Marites Melancon
- Department of Interventional Radiology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - R Jason Stafford
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Junjie Li
- Department of Cancer Systems Imaging, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Alpa M Nick
- Department of Gynecologic Oncology and Reproductive Medicine, The University of Texas MD Anderson Cancer Center, Houston, Texas; and
| | - Mei Tian
- Department of Nuclear Medicine and PET Center, The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
| | - Anil K Sood
- Department of Gynecologic Oncology and Reproductive Medicine, The University of Texas MD Anderson Cancer Center, Houston, Texas; and.,Center for RNAi and Non-Coding RNA, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Chun Li
- Department of Cancer Systems Imaging, The University of Texas MD Anderson Cancer Center, Houston, Texas
| |
Collapse
|
9
|
Closed-Bore Interventional MRI: Percutaneous Biopsies and Ablations. AJR Am J Roentgenol 2015; 205:W400-10. [DOI: 10.2214/ajr.15.14732] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
10
|
Fahrenholtz SJ, Moon TY, Franco M, Medina D, Danish S, Gowda A, Shetty A, Maier F, Hazle JD, Stafford RJ, Warburton T, Fuentes D. A model evaluation study for treatment planning of laser-induced thermal therapy. Int J Hyperthermia 2015; 31:705-14. [PMID: 26368014 DOI: 10.3109/02656736.2015.1055831] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
A cross-validation analysis evaluating computer model prediction accuracy for a priori planning magnetic resonance-guided laser-induced thermal therapy (MRgLITT) procedures in treating focal diseased brain tissue is presented. Two mathematical models are considered. (1) A spectral element discretisation of the transient Pennes bioheat transfer equation is implemented to predict the laser-induced heating in perfused tissue. (2) A closed-form algorithm for predicting the steady-state heat transfer from a linear superposition of analytic point source heating functions is also considered. Prediction accuracy is retrospectively evaluated via leave-one-out cross-validation (LOOCV). Modelling predictions are quantitatively evaluated in terms of a Dice similarity coefficient (DSC) between the simulated thermal dose and thermal dose information contained within N = 22 MR thermometry datasets. During LOOCV analysis, the transient model's DSC mean and median are 0.7323 and 0.8001 respectively, with 15 of 22 DSC values exceeding the success criterion of DSC ≥ 0.7. The steady-state model's DSC mean and median are 0.6431 and 0.6770 respectively, with 10 of 22 passing. A one-sample, one-sided Wilcoxon signed-rank test indicates that the transient finite element method model achieves the prediction success criteria, DSC ≥ 0.7, at a statistically significant level.
Collapse
Affiliation(s)
- Samuel J Fahrenholtz
- a Department of Imaging Physics , M.D. Anderson Cancer Center, University of Texas , Houston , Texas , USA .,b Graduate School of Biomedical Sciences, University of Texas , Houston , Texas , USA
| | - Tim Y Moon
- c Department of Computational and Applied Mathematics , Rice University , Houston , Texas , USA
| | - Michael Franco
- c Department of Computational and Applied Mathematics , Rice University , Houston , Texas , USA
| | - David Medina
- c Department of Computational and Applied Mathematics , Rice University , Houston , Texas , USA
| | - Shabbar Danish
- d Department of Neurosurgery , Robert Wood Johnson Hospital , New Brunswick, New Jersey , USA , and
| | | | | | - Florian Maier
- a Department of Imaging Physics , M.D. Anderson Cancer Center, University of Texas , Houston , Texas , USA
| | - John D Hazle
- a Department of Imaging Physics , M.D. Anderson Cancer Center, University of Texas , Houston , Texas , USA .,b Graduate School of Biomedical Sciences, University of Texas , Houston , Texas , USA
| | - Roger J Stafford
- a Department of Imaging Physics , M.D. Anderson Cancer Center, University of Texas , Houston , Texas , USA .,b Graduate School of Biomedical Sciences, University of Texas , Houston , Texas , USA
| | - Tim Warburton
- c Department of Computational and Applied Mathematics , Rice University , Houston , Texas , USA
| | - David Fuentes
- a Department of Imaging Physics , M.D. Anderson Cancer Center, University of Texas , Houston , Texas , USA .,b Graduate School of Biomedical Sciences, University of Texas , Houston , Texas , USA
| |
Collapse
|
11
|
Meshless method with operator splitting technique for transient nonlinear bioheat transfer in two-dimensional skin tissues. Int J Mol Sci 2015; 16:2001-19. [PMID: 25603180 PMCID: PMC4307345 DOI: 10.3390/ijms16012001] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2014] [Accepted: 01/07/2015] [Indexed: 11/21/2022] Open
Abstract
A meshless numerical scheme combining the operator splitting method (OSM), the radial basis function (RBF) interpolation, and the method of fundamental solutions (MFS) is developed for solving transient nonlinear bioheat problems in two-dimensional (2D) skin tissues. In the numerical scheme, the nonlinearity caused by linear and exponential relationships of temperature-dependent blood perfusion rate (TDBPR) is taken into consideration. In the analysis, the OSM is used first to separate the Laplacian operator and the nonlinear source term, and then the second-order time-stepping schemes are employed for approximating two splitting operators to convert the original governing equation into a linear nonhomogeneous Helmholtz-type governing equation (NHGE) at each time step. Subsequently, the RBF interpolation and the MFS involving the fundamental solution of the Laplace equation are respectively employed to obtain approximated particular and homogeneous solutions of the nonhomogeneous Helmholtz-type governing equation. Finally, the full fields consisting of the particular and homogeneous solutions are enforced to fit the NHGE at interpolation points and the boundary conditions at boundary collocations for determining unknowns at each time step. The proposed method is verified by comparison of other methods. Furthermore, the sensitivity of the coefficients in the cases of a linear and an exponential relationship of TDBPR is investigated to reveal their bioheat effect on the skin tissue.
Collapse
|
12
|
Wichmann JL, Beeres M, Borchard BM, Naguib NNN, Bodelle B, Lee C, Zangos S, Vogl TJ, Mack MG, Eichler K. Evaluation of MRI T1-based treatment monitoring during laser-induced thermotherapy of liver metastases for necrotic size prediction. Int J Hyperthermia 2013; 30:19-26. [DOI: 10.3109/02656736.2013.854931] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
|
13
|
Fernando R, Downs J, Maples D, Ranjan A. MRI-guided monitoring of thermal dose and targeted drug delivery for cancer therapy. Pharm Res 2013; 30:2709-17. [PMID: 23780716 DOI: 10.1007/s11095-013-1110-8] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2013] [Accepted: 06/04/2013] [Indexed: 12/14/2022]
Abstract
Application of localized hyperthermia treatment for solid tumor therapy is under active clinical investigation. The success of this treatment methodology, whether for tumor ablation or drug delivery, requires accurate target localization and real-time temperature mapping of the targeted region. Magnetic Resonance Imaging (MRI) can monitor temperature elevations in tissues in real-time during tumor therapy. MRI can also be applied in concert with methods such as High Intensity Focused Ultrasound (HIFU) to enable image-guided drug delivery (IGDD) from temperature sensitive nanocarriers, by exploiting not only its anatomic resolution, but its ability to detect and measure drug release using markers co-loaded with drugs within the nanocarriers. We review this rapidly emerging technology, providing an overview of MRI-guided tissue thermal dose monitoring for HIFU and Laser therapy, its role in targeted drug delivery and its future potential for clinical translation.
Collapse
Affiliation(s)
- Ruchika Fernando
- Laboratory of Nanomedicine & Targeted Therapy Department of Physiological Sciences, Oklahoma State University, Stillwater, Oklahoma, USA
| | | | | | | |
Collapse
|
14
|
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
|
15
|
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
|
16
|
Focal laser ablation of prostate cancer: definition, needs, and future. Adv Urol 2012; 2012:589160. [PMID: 22666240 PMCID: PMC3362007 DOI: 10.1155/2012/589160] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2012] [Revised: 03/13/2012] [Accepted: 03/14/2012] [Indexed: 01/13/2023] Open
Abstract
Current challenges and innovations in prostate cancer management concern the development of focal therapies that allow the treatment of only the cancer areas sparing the rest of the gland to minimize the potential morbidity. Among these techniques, focal laser ablation (FLA) appears as a potential candidate to reach the goal of focusing energy delivery on the identified targets. The aim of this study is to perform an up-to-date review of this new therapeutic modality. Relevant literature was identified using MEDLINE database with no language restrictions (entries: focal therapy, laser interstitial thermotherapy, prostate cancer, FLA) and by cross-referencing from previously identified studies. Precision, real-time monitoring, MRI compatibility, and low cost of integrated system are principal advantages of FLA. Feasibility and safety of this technique have been reported in phase I assays. FLA might eventually prove to be a middle ground between active surveillance and radical treatment. In conclusion, FLA may have found a role in the management of prostate cancer. However, further trials are required to demonstrate the oncologic effectiveness in the long term.
Collapse
|
17
|
Fuentes D, Yung J, Hazle JD, Weinberg JS, Stafford RJ. Kalman filtered MR temperature imaging for laser induced thermal therapies. IEEE TRANSACTIONS ON MEDICAL IMAGING 2012; 31:984-94. [PMID: 22203706 PMCID: PMC3873725 DOI: 10.1109/tmi.2011.2181185] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
The feasibility of using a stochastic form of Pennes bioheat model within a 3-D finite element based Kalman filter (KF) algorithm is critically evaluated for the ability to provide temperature field estimates in the event of magnetic resonance temperature imaging (MRTI) data loss during laser induced thermal therapy (LITT). The ability to recover missing MRTI data was analyzed by systematically removing spatiotemporal information from a clinical MR-guided LITT procedure in human brain and comparing predictions in these regions to the original measurements. Performance was quantitatively evaluated in terms of a dimensionless L(2) (RMS) norm of the temperature error weighted by acquisition uncertainty. During periods of no data corruption, observed error histories demonstrate that the Kalman algorithm does not alter the high quality temperature measurement provided by MR thermal imaging. The KF-MRTI implementation considered is seen to predict the bioheat transfer with RMS error < 4 for a short period of time, ∆t < 10 s, until the data corruption subsides. In its present form, the KF-MRTI method currently fails to compensate for consecutive for consecutive time periods of data loss ∆t > 10 sec.
Collapse
Affiliation(s)
- D. Fuentes
- The University of Texas M.D Anderson Cancer Center, Department of Imaging Physics, Houston TX 77030, USA
| | - J. Yung
- The University of Texas M.D Anderson Cancer Center, Department of Imaging Physics, Houston TX 77030, USA
| | - J. D. Hazle
- The University of Texas M.D Anderson Cancer Center, Department of Imaging Physics, Houston TX 77030, USA
| | - J. S. Weinberg
- The University of Texas M.D Anderson Cancer Center, Department of Neurosurgery, Houston TX 77030, USA
| | - R. J. Stafford
- The University of Texas M.D Anderson Cancer Center, Department of Imaging Physics, Houston TX 77030, USA
| |
Collapse
|
18
|
Hirata A, Masuda H, Kanai Y, Asai R, Fujiwara O, Arima T, Kawai H, Watanabe S, Lagroye I, Veyret B. Computational modeling of temperature elevation and thermoregulatory response in the brains of anesthetized rats locally exposed at 1.5 GHz. Phys Med Biol 2012; 56:7639-57. [PMID: 22086327 DOI: 10.1088/0031-9155/56/23/019] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
The dominant effect of human exposures to microwaves is caused by temperature elevation ('thermal effect'). In the safety guidelines/standards, the specific absorption rate averaged over a specific volume is used as a metric for human protection from localized exposure. Further investigation on the use of this metric is required, especially in terms of thermophysiology. The World Health Organization (2006 RF research agenda) has given high priority to research into the extent and consequences of microwave-induced temperature elevation in children. In this study, an electromagnetic-thermal computational code was developed to model electromagnetic power absorption and resulting temperature elevation leading to changes in active blood flow in response to localized 1.457 GHz exposure in rat heads. Both juvenile (4 week old) and young adult (8 week old) rats were considered. The computational code was validated against measurements for 4 and 8 week old rats. Our computational results suggest that the blood flow rate depends on both brain and core temperature elevations. No significant difference was observed between thermophysiological responses in 4 and 8 week old rats under these exposure conditions. The computational model developed herein is thus applicable to set exposure conditions for rats in laboratory investigations, as well as in planning treatment protocols in the thermal therapy.
Collapse
Affiliation(s)
- Akimasa Hirata
- Department of Computer Science and Engineering, Nagoya Institute of Technology, Nagoya 466-8555, Japan.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
19
|
Feng Y, Fuentes D. Model-based planning and real-time predictive control for laser-induced thermal therapy. Int J Hyperthermia 2011; 27:751-61. [PMID: 22098360 PMCID: PMC3930104 DOI: 10.3109/02656736.2011.611962] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
In this article, the major idea and mathematical aspects of model-based planning and real-time predictive control for laser-induced thermal therapy (LITT) are presented. In particular, a computational framework and its major components developed by authors in recent years are reviewed. The framework provides the backbone for not only treatment planning but also real-time surgical monitoring and control with a focus on MR thermometry enabled predictive control and applications to image-guided LITT, or MRgLITT. Although this computational framework is designed for LITT in treating prostate cancer, it is further applicable to other thermal therapies in focal lesions induced by radio-frequency (RF), microwave and high-intensity-focused ultrasound (HIFU). Moreover, the model-based dynamic closed-loop predictive control algorithms in the framework, facilitated by the coupling of mathematical modelling and computer simulation with real-time imaging feedback, has great potential to enable a novel methodology in thermal medicine. Such technology could dramatically increase treatment efficacy and reduce morbidity.
Collapse
Affiliation(s)
- Yusheng Feng
- Computational Bioengineering and Control Lab, The University of Texas at San Antonio, USA.
| | | |
Collapse
|