Accelerated MR Thermometry in the Presence of Uncertainties
Reza Madankan1, Wolfgang Stefan1, Christopher MacLellan1, Samuel Fahrenholtz1, Drew Mitchell1, R.J. Stafford1, John Hazle1, and David Fuentes1

1Imaging Physics, MD Anderson Cancer Center, Houston, TX, United States

Synopsis

Compressive sensing and sparse image reconstruction has received significant attention and has demonstrated potential in reduction of acquisition times. However, in many methods, under-sampling strategies are heuristically chosen and empirically validated. This often leads to a relatively larger number of k-space samples than needed for a particular application. The presented work develops a mathematically rigorous and quantitative methodology for k-space under-sampling with respect to model-based reconstruction of MR thermometry. The key idea of the proposed approach is to detect the useful samples of k-space in order to refine the model, and then the refined mathematical model is utilized to reconstruct the image.

Target Audience

Researchers interested in identifying and quantifying the optimal k-space locations for model-based reconstruction of MR thermometry data.

Purpose

Compressive sensing and sparse image reconstruction has received significant attention and has demonstrated potential in reduction of acquisition times. However, in many methods, under-sampling strategies are heuristically chosen and empirically validated. This often leads to a relatively larger number of k-space samples than needed for a particular application. The presented work develops a mathematically rigorous and quantitative methodology for k-space under-sampling with respect to model-based reconstruction of MR thermometry. Hence, the key question here is that whether there exist optimal locations in k-space such that they provide maximum information regarding the image? And if they exist, how one can develop a general framework to find these locations on k-space?

Methods

An accelerated model-based information theoretic approach is developed to perform the task of image reconstruction from observed k-space samples. The key idea of the proposed approach is to optimally detect the useful samples of k-space in order to refine the model, and then the refined mathematical model is utilized to reconstruct the image1. Uncertainty quantification techniques2 along with information theoretic concepts, i.e, entropy, are used to locate k-space regions that demonstrate the highest sensitivity with respect to a composite heat transfer and MR physics. The key contribution of this method is to optimally select the samples on k-space such that they result in the best estimate of the corresponding temperature map. Optimally selected observations on k-space are utilized to refine the mathematical model and the refined model is then used to generate the temperature field corresponding to the observations. Fig.1 illustrates schematic view for implementation of the proposed algorithm for brain tumor laser ablation. Note that only small fraction of lines in k-space are used for data acquisition (bottom left). Acquired measurement data is then merged with model prediction data (top left) in a minimum variance framework to estimate the temperature filed over the tissue (top right). Corresponding error between the estimated and the true temperature field is also shown in Fig.1 (bottom right).

Results

An in-silico simulation is compared to MR temperature imaging from a human brain ablation. Fig.2 shows the position of the tumor in brain. Performance of the proposed technique is retrospectively validated in fully sampled planar temperature imaging acquired during the thermal ablation process. Fully sampled k-space data acquired in-vivo3 are used as a control for the computed sub-sampling strategies. MR temperature imaging was performed on a 1.5 T MRI (ExciteHD®, General Electric, Milwaukee, WI) GE scanner at every Δt =15 seconds (TR/FA = 38 ms/30º, frequency × phase = 256 × 128, FOV = 26 cm × 26 cm, BW = 100 Hz/pixel, slice thickness 5 mm). Fig.3 represents acquired k-space data points based on the proposed technique and other conventional sampling methods. The error in temperature estimate, while using different k-space data acquisition techniques, is also shown in Fig. 4. As expected, proposed approach results in minimum error in comparison with other subsampling schemes.

Conclusion

The novelty of the proposed approach demonstrates that locations of high-information content with respect to a model based reconstruction of MR thermometry may be quantitatively identified. The information locations identified are a consequence of the physics based modeling of the laser induced heating and an intelligent characterization of the measurement uncertainties of the acquisition system. In particular, efficient k-space points are selected by variance maximization of model-based prediction of the MR signal. Our mathematical model of the MR thermometry acquisition is refined using this set of judiciously selected data observations to reconstruct the thermal image. Intuitively, the precision and confidence of the thermal image reconstruction depends on multiple factors like accuracy, information content, and the number of data observations. Optimally selected k-space samples allow reconstruction of high quality MR thermometry images with significantly less data compared to the fully sampled control data.

Acknowledgements

No acknowledgement found.

References

1. R. Madankan, W. Stefan, S. Fahrenholtz, C. MacLellan, J. Hazle, J. Stafford, J. S. Weinberg, G. Rao, and D. Fuentes, Accelerated Magnetic Resonance Thermometry in Presence of Uncertainties, Arxiv Preprint, http://arxiv.org/abs/1510.08875, 2015.

2. Madankan, R., et al. Computation of probabilistic hazard maps and source parameter estimation for volcanic ash transport and dispersion. J. of Computational Physics 271 (2014): 39-59.

3. Fahrenholtz, Samuel J., et al. "A model evaluation study for treatment planning of laser-induced thermal therapy." International Journal of Hyperthermia 31.7 (2015): 705-714.

Figures

Fig.1: Schematic view of information theory based data assimilation. Model predictions of temperature based on prior data (top left) and judiciously selected MRI data (bottom left) are merged in a minimum variance framework (center) to provide posterior estimate of the temperature (top right). The error between temperature estimate and its actual value is shown in bottom right.

Fig.2: Position of the tumor and the temperature field resulting from MRgLITT (FOV = 26 cm2). Black box and green line inside represent the region of interest and laser fiber, respectively.

Fig. 3: Planar Image Reconstruction: a) The 20 readout lines used for data acquisition, obtained from proposed approach. b) Uniformly distributed readout lines based on rectilinear undersampling. c) Readout lines based on variable-density Poisson disk undersampling.

Fig. 4: The error between the posterior estimate and actual measurement of temperature field, obtained by using a) the proposed technique, b) rectilinear undersampling, and c) variable-density Poisson disk undersampling. Note the significant amount of error in center of ROI for cases a and b. High value of the error in borders of the images is due to presence of artifacts (ROI = 55 mm 60 mm).



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
2107