Magnetic Resonance Temperature Imaging Research Articles

TU‐A‐9A‐04: Development of a Thermally Stable Phantom for Photoacoustic and Magnetic Resonance Temperature Imaging

Purpose:Photoacoustic‐ultrasonic (PAUS) imaging, which utilizes an ultrasound transducer to provide co‐registered photoacoustic and pulse‐echo ultrasound images, is capable of measuring temperature non‐invasively while simultaneously providing anatomical images. The sub‐millimeter resolution and centimeter‐order penetration depths achievable with PAUS imaging have the potential to deliver active monitoring of both a targeted tumor microenvironment and nearby healthy tissue during thermal ablation. These characteristics make PAUS imaging a promising new technique for guidance and monitoring during photothermal ablations of solid tumors. In order to assess the potential clinical role of PAUS imaging, the technique was validated against the clinically accepted magnetic resonance thermal imaging approach (MRTI).Methods:To facilitate co‐registration between modalities, the phantom had inclusions of gold nanoshells encapsulating super‐paramagnetic iron oxide (SPIO) particles, where gold enhances the PA signal and SPIOs provide negative contrast on MRI. Several phantom designs were assessed for resilience to heating. PA images were acquired on a Vevo LAZR (FUJIFILM VisualSonics Inc., Toronto, Ontario) PAultrasound small‐animal imaging system (21MHz) operating at 710nm. MRTI experiments were performed using a 6‐channel flex coil (GE Healthcare, Waukesha, WI) on a 3T MRI scanner (Discovery MR750, GE Healthcare, Waukesha, WI) using a fast multi gradient echo acquisition (16 echoes, 128×128 acquisition matrix, 25.6×25.6cm field of view, 3mm slice thickness, 60ms TR, 20° flip angle, 2.9ms minimum TE and 1.6ms echo spacing). The accuracy and spatio‐temporal resolution of PA thermography was cross‐validated with both MRTI and a fluoroptic temperature sensor (LumaSense Technologies, Santa Clara, CA) in the custom‐designed phantom.Results:A thermally stable, dual‐modality phantom was created for cross‐validation of photoacoustic thermography and MRTI. Axial and lateral resolutions of PA images were sub‐millimeter with a temporal resolution of 0.2s, which will accommodate precise real‐time guidance and monitoring.Conclusion:These results indicate that a PA thermography technique offers tremendous promise for real‐time thermal monitoring of ablative therapy.Funding support provided by Cancer Prevention Research Institute of Texas and Julia Jones Matthews family. No disclosures or conflicts of interest.

SU‐E‐J‐161: Inverse Problems for Optical Parameters in Laser Induced Thermal Therapy

Purpose:Magnetic resonance‐guided laser‐induced thermal therapy (MRgLITT) is investigated as a neurosurgical intervention for oncological applications throughout the body by active post market studies. Real‐time MR temperature imaging is used to monitor ablative thermal delivery in the clinic. Additionally, brain MRgLITT could improve through effective planning for laser fiber's placement. Mathematical bioheat models have been extensively investigated but require reliable patient specific physical parameter data, e.g. optical parameters. This abstract applies an inverse problem algorithm to characterize optical parameter data obtained from previous MRgLITT interventions.Methods:The implemented inverse problem has three primary components: a parameter‐space search algorithm, a physics model, and training data. First, the parameter‐space search algorithm uses a gradient‐based quasi‐Newton method to optimize the effective optical attenuation coefficient, μ_eff. A parameter reduction reduces the amount of optical parameter‐space the algorithm must search. Second, the physics model is a simplified bioheat model for homogeneous tissue where closed‐form Green's functions represent the exact solution. Third, the training data was temperature imaging data from 23 MRgLITT oncological brain ablations (980 nm wavelength) from seven different patients.Results:To three significant figures, the descriptive statistics for μ_eff were 1470 m−1 mean, 1360 m−1 median, 369 m−1 standard deviation, 933 m−1 minimum and 2260 m−1 maximum. The standard deviation normalized by the mean was 25.0%. The inverse problem took <30 minutes to optimize all 23 datasets.Conclusion:As expected, the inferred average is biased by underlying physics model. However, the standard deviation normalized by the mean is smaller than literature values and indicates an increased precision in the characterization of the optical parameters needed to plan MRgLITT procedures. This investigation demonstrates the potential for the optimization and validation of more sophisticated bioheat models that incorporate the uncertainty of the data into the predictions, e.g. stochastic finite element methods.

Robust phase unwrapping for MR temperature imaging using a magnitude-sorted list, multi-clustering algorithm.

Several methods in MRI use the phase information of the complex signal and require phase unwrapping (e.g., B0 field mapping, chemical shift imaging, and velocity measurements). In this work, an algorithm was developed focusing on the needs and requirements of MR temperature imaging applications. The proposed method performs fully automatic unwrapping using a list of all pixels sorted by magnitude in descending order and creates and merges clusters of unwrapped pixels until the entire image is unwrapped. The algorithm was evaluated using simulated phantom data and in vivo clinical temperature imaging data. The evaluation of the phantom data demonstrated no errors in regions with signal-to-noise ratios of at least 4.5. For the in vivo data, the algorithm did not fail at an average of more than one pixel for signal-to-noise ratios greater than 6.3. Processing times less than 30 ms per image were achieved by unwrapping pixels inside a region of interest (53 × 53 pixels) used for referenceless MR temperature imaging. The algorithm has been demonstrated to operate robustly with clinical in vivo data in this study. The processing time for common regions of interest in referenceless MR temperature imaging allows for online updates of temperature maps without noticeable delay.

Absolute MR thermometry using nanocarriers

Accurate time-resolved temperature mapping is crucial for the safe use of hyperthermia-mediated drug delivery. We here propose a magnetic resonance imaging temperature mapping method in which drug delivery systems serve not only to improve tumor targeting, but also as an accurate and absolute nano-thermometer. This method is based on the temperature-dependent chemical shift difference between water protons and the protons in different groups of drug delivery systems. We show that the chemical shift of the protons in the ethylene oxide group in polyethylene glycol (PEG) is temperature-independent, whereas the proton resonance of water decreases with increasing temperature. The frequency difference between both resonances is linear and does not depend on pH and physiological salt conditions. In addition, we show that the proton resonance of the methyl group in N-(2-hydroxypropyl)-methacrylamide (HPMA) is temperature-independent. Therefore, PEGylated liposomes, polymeric mPEG-b-pHPMAm-Lac2 micelles and HPMA copolymers can provide a temperature-independent reference frequency for absolute magnetic resonance (MR) thermometry. Subsequently, we show that multigradient echo MR imaging with PEGylated liposomes in situ allows accurate, time-resolved temperature mapping. In conclusion, nanocarrier materials may serve as highly versatile tools for tumor-targeted drug delivery, acting not only as hyperthermia-responsive drug delivery systems, but also as accurate and precise nano-thermometers.

Magnetic Resonance-Guided Laser Induced Thermal Therapy for Glioblastoma Multiforme: A Review

Magnetic resonance-guided laser induced thermotherapy (MRgLITT) has become an increasingly relevant therapy for tumor ablation due to its minimally invasive approach and broad applicability across many tissue types. The current state of the art applies laser irradiation via cooled optical fiber applicators in order to generate ablative heat and necrosis in tumor tissue. Magnetic resonance temperature imaging (MRTI) is used concurrently with this therapy to plan treatments and visualize tumor necrosis. Though application in neurosurgery remains in its infancy, MRgLITT has been found to be a promising therapy for many types of brain tumors. This review examines the current use of MRgLITT with regard to the special clinical challenge of glioblastoma multiforme and examines the potential applications of next-generation nanotherapy specific to the treatment of glioblastoma.

Estimating nanoparticle optical absorption with magnetic resonance temperature imaging and bioheat transfer simulation

Purpose: Optically activated nanoparticle-mediated heating for thermal therapy applications is an area of intense research. The ability to characterise the spatio-temporal heating potential of these particles for use in modelling under various exposure conditions can aid in the exploration of new approaches for therapy as well as more quantitative prospective approaches to treatment planning. The purpose of this research was to investigate an inverse solution to the heat equation using magnetic resonance temperature imaging (MRTI) feedback, for providing optical characterisation of two types of nanoparticles (gold–silica nanoshells and gold nanorods).Methods: The optical absorption of homogeneous nanoparticle–agar mixtures was measured during exposure to an 808 nm laser using real-time MRTI. A coupled finite element solution of heat transfer was registered with the data and used to solve the inverse problem. The L2 norm of the difference between the temperature increase in the model and MRTI was minimised using a pattern search algorithm by varying the absorption coefficient of the mixture.Results: Absorption fractions were within 10% of literature values for similar nanoparticles. Comparison of temporal and spatial profiles demonstrated good qualitative agreement between the model and the MRTI. The weighted root mean square error was <1.5 σMRTI and the average Dice similarity coefficient for ΔT = 5 °C isotherms was >0.9 over the measured time interval.Conclusion: This research demonstrates the feasibility of using an indirect method for making minimally invasive estimates of nanoparticle absorption that might be expanded to analyse a variety of geometries and particles of interest.

Catheter-based and endoluminal ultrasound applicators for magnetic resonance image-guided thermal therapy of pancreatic cancer: Preliminary investigations

Ultrasound devices are being investigated for endoluminal and intraductal access for targeted thermal ablation or hyperthermia of pancreas under MR guidance and temperature monitoring. Simulations using patient-specific 3D models were developed for applicator design and development of treatment delivery strategies. MR-compatible devices were constructed for endoluminal (3-5 MHz planar or lightly focused rectangular elements, 12-mm OD assembly, expandable balloon), transgastric interstitial and intraductal (6-8 MHz multi-sectored tubular elements, 2-mm catheter) deployment. Micro-coils were integrated for active MR tracking of position and alignment. The proof-of-concept devices were tested in phantoms, ex vivo tissues, cadaveric porcine models, and in vivo animal models under 3T MR temperature imaging (MRTI). Results indicate endoluminal devices could ablate 2-2.5 cm depth from gastric wall for tumors of the pancreatic head, and multi-sectored tubular intraductal and interstitial applicators could ablate 2.3-3.4 cm diameter targets with directional control. Intraductal applicators could produce effective hyperthermia (>40 C) extending 15 mm radial. Customized tracking sequences could be used to locate 3D position of the applicators. Endoluminal, interstitial, and intraductal ultrasound applicators show promise for ablation or hyperthermia of pancreatic tumors. MR guidance can be employed for positioning these devices with active tracking coils and real time temperature monitoring. (NIH-P01CA159992.)

Integrated ultrasound and magnetic resonance imaging for simultaneous temperature and cavitation monitoring during focused ultrasound therapies

Ultrasound can be used to noninvasively produce different bioeffects via viscous heating, acoustic cavitation, or their combination, and these effects can be exploited to develop a wide range of therapies for cancer and other disorders. In order to accurately localize and control these different effects, imaging methods are desired that can map both temperature changes and cavitation activity. To address these needs, the authors integrated an ultrasound imaging array into an MRI-guided focused ultrasound (MRgFUS) system to simultaneously visualize thermal and mechanical effects via passive acoustic mapping (PAM) and MR temperature imaging (MRTI), respectively. The system was tested with an MRgFUS system developed for transcranial sonication for brain tumor ablation in experiments with a tissue mimicking phantom and a phantom-filled ex vivo macaque skull. In experiments on cavitation-enhanced heating, 10 s continuous wave sonications were applied at increasing power levels (30-110 W) until broadband acoustic emissions (a signature for inertial cavitation) were evident. The presence or lack of signal in the PAM, as well as its magnitude and location, were compared to the focal heating in the MRTI. Additional experiments compared PAM with standard B-mode ultrasound imaging and tested the feasibility of the system to map cavitation activity produced during low-power (5 W) burst sonications in a channel filled with a microbubble ultrasound contrast agent. When inertial cavitation was evident, localized activity was present in PAM and a marked increase in heating was observed in MRTI. The location of the cavitation activity and heating agreed on average after registration of the two imaging modalities; the distance between the maximum cavitation activity and focal heating was -3.4 ± 2.1 mm and -0.1 ± 3.3 mm in the axial and transverse ultrasound array directions, respectively. Distortions and other MRI issues introduced small uncertainties in the PAM∕MRTI registration. Although there was substantial variation, a nonlinear relationship between the average intensity of the cavitation maps, which was relatively constant during sonication, and the peak temperature rise was evident. A fit to the data to an exponential had a correlation coefficient (R(2)) of 0.62. The system was also found to be capable of visualizing cavitation activity with B-mode imaging and of passively mapping cavitation activity transcranially during cavitation-enhanced heating and during low-power sonication with an ultrasound contrast agent. The authors have demonstrated the feasibility of integrating an ultrasound imaging array into an MRgFUS system to simultaneously map localized cavitation activity and temperature. The authors anticipate that this integrated approach can be utilized to develop controllers for cavitation-enhanced ablation and facilitate the optimization and development of this and other ultrasound therapies. The integrated system may also provide a useful tool to study the bioeffects of acoustic cavitation.

Approaches for modelling interstitial ultrasound ablation of tumours within or adjacent to bone: Theoretical and experimental evaluations

Purpose: The objectives of this study were to develop numerical models of interstitial ultrasound ablation of tumours within or adjacent to bone, to evaluate model performance through theoretical analysis, and to validate the models and approximations used through comparison to experiments. Methods: 3D transient biothermal and acoustic finite element models were developed, employing four approximations of 7-MHz ultrasound propagation at bone/soft tissue interfaces. The various approximations considered or excluded reflection, refraction, angle-dependence of transmission coefficients, shear mode conversion, and volumetric heat deposition. Simulations were performed for parametric and comparative studies. Experiments within ex vivo tissues and phantoms were performed to validate the models by comparison to simulations. Temperature measurements were conducted using needle thermocouples or magnetic resonance temperature imaging (MRTI). Finite element models representing heterogeneous tissue geometries were created based on segmented MR images. Results: High ultrasound absorption at bone/soft tissue interfaces increased the volumes of target tissue that could be ablated. Models using simplified approximations produced temperature profiles closely matching both more comprehensive models and experimental results, with good agreement between 3D calculations and MRTI. The correlation coefficients between simulated and measured temperature profiles in phantoms ranged from 0.852 to 0.967 (p-value < 0.01) for the four models. Conclusions: Models using approximations of interstitial ultrasound energy deposition around bone/soft tissue interfaces produced temperature distributions in close agreement with comprehensive simulations and experimental measurements. These models may be applied to accurately predict temperatures produced by interstitial ultrasound ablation of tumours near and within bone, with applications toward treatment planning.

Design and initial evaluation of a treatment planning software system for MRI-guided laser ablation in the brain

An open-source software system for planning magnetic resonance (MR)-guided laser-induced thermal therapy (MRgLITT) in brain is presented. The system was designed to provide a streamlined and operator-friendly graphical user interface (GUI) for simulating and visualizing potential outcomes of various treatment scenarios to aid in decisions on treatment approach or feasibility. A portable software module was developed on the 3D Slicer platform, an open-source medical imaging and visualization framework. The module introduces an interactive GUI for investigating different laser positions and power settings as well as the influence of patient-specific tissue properties for quickly creating and evaluating custom treatment options. It also provides a common treatment planning interface for use by both open-source and commercial finite element solvers. In this study, an open-source finite element solver for Pennes' bioheat equation is interfaced to the module to provide rapid 3D estimates of the steady-state temperature distribution and potential tissue damage in the presence of patient-specific tissue boundary conditions identified on segmented MR images. The total time to initialize and simulate an MRgLITT procedure using the GUI was [Formula: see text]5 min. Each independent simulation took [Formula: see text]30 s, including the time to visualize the results fused with the planning MRI. For demonstration purposes, a simulated steady-state isotherm contour [Formula: see text] was correlated with MR temperature imaging (N = 5). The mean Hausdorff distance between simulated and actual contours was 2.0 mm [Formula: see text], whereas the mean Dice similarity coefficient was 0.93 [Formula: see text]. We have designed, implemented, and conducted initial feasibility evaluations of a software tool for intuitive and rapid planning of MRgLITT in brain. The retrospective in vivo dataset presented herein illustrates the feasibility and potential of incorporating fast, image-based bioheat predictions into an interactive virtual planning environment for such procedures.

Effects of MRTI sampling characteristics on estimation of HIFU SAR and tissue thermal diffusivity

While the non-invasive and three-dimensional nature of magnetic-resonance temperature imaging (MRTI) makes it a valuable tool for high-intensity focused ultrasound (HIFU) treatments, random and systematic errors in MRTI measurements may propagate into temperature-based parameter estimates used for pretreatment planning. This study assesses the MRTI effects of zero-mean Gaussian noise (SD = 0.0–2.0 °C), temporal sampling (tacq = 1.0–8.0 s), and spatial averaging (Res = 0.5–2.0 mm isotropic) on HIFU temperature measurements and temperature-based estimates of the amplitude and full width half maximum (FWHM) of the HIFU specific absorption rate and of tissue thermal diffusivity. The ultrasound beam used in simulations and ex vivo pork loin experiments has lateral and axial FWHM dimensions of 1.4 mm and 7.9 mm respectively. For spatial averaging simulations, beams with lateral FWHM varying from 1.2–2.2 mm are also assessed. Under noisy conditions, parameter estimates are improved by fitting to data from larger voxel regions. Varying the temporal sampling results in minimal changes in measured temperatures (<2% change) and parameter estimates (<5% change). For the HIFU beams studied, a spatial resolution of 1 × 1 × 3 mm3 or smaller is required to keep errors in temperature and all estimated parameters less than 10%. By quantifying the errors associated with these sampling characteristics, this work provides researchers with appropriate MRTI conditions for obtaining estimates of parameters essential to pretreatment modeling of HIFU thermal therapies.

Multiparametric fat–water separation method for fast chemical‐shift imaging guidance of thermal therapies

A k-means-based classification algorithm is investigated to assess suitability for rapidly separating and classifying fat/water spectral peaks from a fast chemical shift imaging technique for magnetic resonance temperature imaging. Algorithm testing is performed in simulated mathematical phantoms and agar gel phantoms containing mixed fat/water regions. Proton resonance frequencies (PRFs), apparent spin-spin relaxation (T2*) times, and T1-weighted (T1-W) amplitude values were calculated for each voxel using a single-peak autoregressive moving average (ARMA) signal model. These parameters were then used as criteria for k-means sorting, with the results used to determine PRF ranges of each chemical species cluster for further classification. To detect the presence of secondary chemical species, spectral parameters were recalculated when needed using a two-peak ARMA signal model during the subsequent classification steps. Mathematical phantom simulations involved the modulation of signal-to-noise ratios (SNR), maximum PRF shift (MPS) values, analysis window sizes, and frequency expansion factor sizes in order to characterize the algorithm performance across a variety of conditions. In agar, images were collected on a 1.5T clinical MR scanner using acquisition parameters close to simulation, and algorithm performance was assessed by comparing classification results to manually segmented maps of the fat/water regions. Performance was characterized quantitatively using the Dice Similarity Coefficient (DSC), sensitivity, and specificity. The simulated mathematical phantom experiments demonstrated good fat/water separation depending on conditions, specifically high SNR, moderate MPS value, small analysis window size, and low but nonzero frequency expansion factor size. Physical phantom results demonstrated good identification for both water (0.997 ± 0.001, 0.999 ± 0.001, and 0.986 ± 0.001 for DSC, sensitivity, and specificity, respectively) and fat (0.763 ± 0.006, 0.980 ± 0.004, and 0.941 ± 0.002 for DSC, sensitivity, and specificity, respectively). Temperature uncertainties, based on PRF uncertainties from a 5 × 5-voxel ROI, were 0.342 and 0.351°C for pure and mixed fat/water regions, respectively. Algorithm speed was tested using 25 × 25-voxel and whole image ROIs containing both fat and water, resulting in average processing times per acquisition of 2.00 ± 0.07 s and 146 ± 1 s, respectively, using uncompiled MATLAB scripts running on a shared CPU server with eight Intel Xeon(TM) E5640 quad-core processors (2.66 GHz, 12 MB cache) and 12 GB RAM. Results from both the mathematical and physical phantom suggest the k-means-based classification algorithm could be useful for rapid, dynamic imaging in an ROI for thermal interventions. Successful separation of fat/water information would aid in reducing errors from the nontemperature sensitive fat PRF, as well as potentially facilitate using fat as an internal reference for PRF shift thermometry when appropriate. Additionally, the T1-W or R2* signals may be used for monitoring temperature in surrounding adipose tissue.

Magnetic Resonance Temperature Imaging of Laser-Induced Thermotherapy: Assessment of Fast Sequences in Ex Vivo Porcine Liver

To evaluate magnetic resonance sequences for T(1) and proton resonance frequency (PRF) thermometry during laser-induced thermotherapy (LITT) in liver tissue. During LITT (1064 nm; 30 W; 3-cm diffuser; 2-3 min) in ex vivo porcine liver, temperature was measured (25-70°C) utilizing a fiberoptic thermometer and MRI was performed with a 1.5-T scanner through the following sequences: segmented echo planar imaging (seg-EPI) for the PRF method; fast low-angle shot (FLASH), inversion-recovery turbo FLASH (IRTF), saturation-recovery turbo FLASH (SRTF) and true-fast imaging (TRUFI) for the T(1) method. Phase angle and signal amplitude (regarding PRF/T(1)) was recorded in regions of interest, on images under fiberoptic probe tips. Sequences' thermal coefficients were determined by calibrating phase angle and signal amplitude against temperature and subsequently validated. Coefficients of -0.0089 ± 0.0003 ppm °C(-1) (seg-EPI) and -0.917 ± 0.046, -1.166 ± 0.058, -1.038 ± 0.054 and -1.443 ± 0.118°C(-1) (FLASH, IRTF, SRTF and TRUFI, respectively) were obtained. Precisions of 0.71, 1.34, 2.07, 2.44 and 3.21°C and, through Bland-Altman analysis, accuracies of -0.67, 0.79, 1.65, 1.57 and 2.13°C (seg-EPI, FLASH, IRTF, SRTF and TRUFI, respectively) were determined. The PRF method with seg-EPI sequence is preferred for thermometry during LITT owing to higher precision and accuracy. Among T(1)-method sequences, FLASH showed higher accuracy and robustness.

Generalised polynomial chaos-based uncertainty quantification for planning MRgLITT procedures

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.

Temperature imaging of laser-induced thermotherapy (LITT) by MRI: evaluation of different sequences in phantom

The purpose of this study was to evaluate magnetic resonance (MR) temperature imaging of the laser-induced thermotherapy (LITT) comparing the proton resonance frequency (PRF) and T 1 thermometry methods. LITT was applied to a liver-mimicking acrylamide gel phantom. Temperature rise up to 70 °C was measured using a MR-compatible fiber-optic thermometer. MR imaging was performed by a 1.5-T scanner utilizing fast gradient echo sequences including a segmented echo planar imaging (seg-EPI) sequence for PRF and the following sequences for T 1 method: fast low-angle shot (FLASH), inversion recovery turbo flash (IRTF), saturation recovery turbo flash (SRTF), and true fast imaging (TRUFI). Temperature-induced change of the pixel values in circular regions of interest, selected on images under the temperature probe tip, was recorded. For each sequence, a calibration constant could be determined to be -0.0088 ± 0.0002 ppm °C(-1) (EPI), -1.15 ± 0.03 °C(-1) (FLASH), -1.49 ± 0.03 °C(-1) (IRTF), -1.21 ± 0.03 °C(-1) (SRTF), and -2.52 ± 0.12 °C(-1) (TRUFI). These constants were evaluated in further LITT experiments in phantom comparing the calculated temperatures with the fiber optic-measured ones; temperature precisions of 0.60 °C (EPI), 0.81 °C (FLASH), 1.85 °C (IRTF), 1.95 °C (SRTF), and 3.36 °C (TRUFI) were obtained. Furthermore, performing the Bland-Altman analysis, temperature accuracy was determined to be 0.23 °C (EPI), 0.31 °C (FLASH), 1.66 °C (IRTF), 1.19 °C (SRTF), and 3.20 °C (TRUFI). In conclusion, the seg-EPI sequence was found to be more convenient for MR temperature imaging of LITT due to its relatively high precision and accuracy. Among the T 1 method sequences, FLASH showed the highest accuracy and robustness.

Extended Kalman Filtering for Continuous Volumetric MR-Temperature Imaging

Real time magnetic resonance (MR) thermometry has evolved into the method of choice for the guidance of high-intensity focused ultrasound (HIFU) interventions. For this role, MR-thermometry should preferably have a high temporal and spatial resolution and allow observing the temperature over the entire targeted area and its vicinity with a high accuracy. In addition, the precision of real time MR-thermometry for therapy guidance is generally limited by the available signal-to-noise ratio (SNR) and the influence of physiological noise. MR-guided HIFU would benefit of the large coverage volumetric temperature maps, including characterization of volumetric heating trajectories as well as near- and far-field heating. In this paper, continuous volumetric MR-temperature monitoring was obtained as follows. The targeted area was continuously scanned during the heating process by a multi-slice sequence. Measured data and a priori knowledge of 3-D data derived from a forecast based on a physical model were combined using an extended Kalman filter (EKF). The proposed reconstruction improved the temperature measurement resolution and precision while maintaining guaranteed output accuracy. The method was evaluated experimentally ex vivo on a phantom, and in vivo on a porcine kidney, using HIFU heating. On the in vivo experiment, it allowed the reconstruction from a spatio-temporally under-sampled data set (with an update rate for each voxel of 1.143 s) to a 3-D dataset covering a field of view of 142.5×285×54 mm(3) with a voxel size of 3×3×6 mm(3) and a temporal resolution of 0.127 s. The method also provided noise reduction, while having a minimal impact on accuracy and latency.

An Inverse Problem Approach to Recovery of In Vivo Nanoparticle Concentrations from Thermal Image Monitoring of MR-Guided Laser Induced Thermal Therapy

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.

鎖状メチレン基プロトンのスピン格子緩和時間と水プロトンの共鳴周波数を用いた脂肪・水混在組織の非侵襲磁気共鳴温度分布画像化の可能性

脂肪・高含水組織混在部位のための磁気共鳴を用いた非侵襲温度分布可視化技術が提案された. この方法は鎖状メチレン基 (CH2) プロトンのスピン格子緩和時間 (T1) と水 (H2O) プロトンの磁気共鳴周波数の温度依存性を利用する. 多点ディクソン法と多フリップ角法を組み合わせて鎖状メチレン基のT1 と水の共鳴周波数を同時に測定する . 2つの容器に入れたマヨネーズの一方をマイクロ波加温した後, 室温に保ったものと共にガドリニウム水溶液に浸し, 冷却過程において3テスラ磁気共鳴画像化装置にてSPGR (SPoiled Gradient Recalled acquisition in steady state) で撮像した. FOV, 32×32 cm2 ; マトリクス, 64×64 ; 並列画像化因子, 2 ; TR, 36 ms ; TE, 1.15 ms ; フリップ角, 20, 50, 70度とした. 各フリップ角におけるデータをIDEAL (Iterative Decomposition of water and fat with Echo Asymmetry and Least squares estimation) 法のアルゴリズムで処理し, 水, 鎖状メチレン基, 終端メチル基の各成分を分離した. 他の含有量の少ない脂肪プロトン成分は簡単のため無視した. 脂肪の温度分布は3つの異なるフリップ角で得られたメチレン基信号のT1 から求めた. 水の温度分布は水とメチレン基の間の位相差の変化, ないしは水位相の相対変化から得た. 得られたこれらの温度画像をボクセル毎の脂肪・水含有率に応じて重み付け平均し, 最終的な温度画像を得た. この結果, プローブを用いて測定した温度上昇とよく相関する脂肪・水混在領域における温度画像が得られ, 本法の有用性が示された.

Kalman Filtered MR Temperature Imaging for Laser Induced Thermal Therapies

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.

MR temperature imaging of nanoshell mediated laser ablation

Minimally invasive thermal therapy using high-power diode lasers is an active area of clinical research. Gold nanoshells (AuNS) can be tuned to absorb light in the range used for laser ablation and may facilitate more conformal tumor heating and sparing of normal tissue via enhanced tumor specific heating. This concept was investigated in a xenograft model of prostate cancer (PC-3) using MR temperature imaging (MRTI) in a 1.5T scanner to characterize the spatiotemporal temperature distribution resulting from nanoparticle mediated heating. Tumors with and without intravenously injected AuNS were exposed to an external laser tuned to 808 nm for 180 sec at 4 W/cm2 under real-time monitoring with proton resonance frequency shift based MRTI. Microscopy indicated that these nanoparticles (140–150 nm) accumulated passively in the tumor and remained close to the tumor microvasculature. MRTI measured a statistically significant (p < 0.001) increase in maximum temperature in the tumor cortex (mean = 21 ± 7°C) in +AuNS tumors versus control tumors. Analysis of the temperature maps helped demonstrate that the overall distribution of temperature within +AuNS tumors was demonstrably higher versus control, and resulted in damage visible on histopathology. This research demonstrates that passive uptake of intravenously injected AuNS in PC-3 xenografts converts the tumor vasculature into a potent heating source for nanoparticle mediated ablation at power levels which do not generate significant damage in normal tissue. When used in conjunction with MRTI, this has implications for development and validation of more conformal delivery of therapy for interstitial laser ablations.