Bayesian estimation of relaxation times T1 in MR images of irradiated Fricke-agarose gels

The authors present a novel method for processing T1-weighted images acquired with Inversion-Recovery (IR) sequence. The method, developed within the Bayesian framework, takes into account a priori knowledge about the spatial regularity of the parameters to be estimated. Inference is drawn by means of Markov Chains Monte Carlo algorithms. The method has been applied to the processing of IR images from irradiated Fricke-agarose gels, proposed in the past as relative dosimeter to verify radiotherapeutic treatment planning systems. Comparison with results obtained from a standard approach shows that signal-to noise ratio (SNR) is strongly enhanced when the estimation of the longitudinal relaxation rate (R1) is performed with the newly proposed statistical approach. Furthermore, the method allows the use of more complex models of the signal. Finally, an appreciable reduction of total acquisition time can be obtained due to the possibility of using a reduced number of images. The method can also be applied to T1 mapping of other systems.