Identification of controlled-complexity thermal therapy models derived from magnetic resonance thermometry images

Medical imaging provides information valuable in diagnosis, planning, and control of therapies. In this paper, we develop a method that uses a specific type of imaging--the magnetic resonance thermometry--to identify accurate and computationally efficient site and patient-specific computer models for thermal therapies, such as focused ultrasound surgery, hyperthermia, and thermally triggered targeted drug delivery. The developed method uses a sequence of acquired MR thermometry images to identify a treatment model describing the deposition and dissipation of thermal energy in tissues. The proper orthogonal decomposition of thermal images is first used to identify a set of empirical eigenfunctions, which captures spatial correlations in the thermal response of tissues. Using the reduced subset of eigenfunction as a functional basis, low-dimensional thermal response and the ultrasound specific absorption rate models are then identified. Once identified, the treatment models can be used to plan, optimize, and control the treatment. The developed approach is validated experimentally using the results of MR thermal imaging of a tissue phantom during focused ultrasound sonication. The validation demonstrates that our approach produces accurate low-dimensional treatment models and provides a convenient tool for balancing the accuracy of model predictions and the computational complexity of the treatment models.