A unifying probabilistic Bayesian approach to derive electron density from MRI for radiation therapy treatment planning

MRI significantly improves the accuracy and reliability of target delineation in radiation therapy for certain tumors due to its superior soft tissue contrast compared to CT. A treatment planning process with MRI as the sole imaging modality will eliminate systematic CT/MRI co-registration errors, reduce cost and radiation exposure, and simplify clinical workflow. However, MRI lacks the key electron density information necessary for accurate dose calculation and generating reference images for patient setup. The purpose of this work is to develop a unifying method to derive electron density from standard T1-weighted MRI. We propose to combine both intensity and geometry information into a unifying probabilistic Bayesian framework for electron density mapping. For each voxel, we compute two conditional probability density functions (PDFs) of electron density given its: (1) T1-weighted MRI intensity, and (2) geometry in a reference anatomy, obtained by deformable image registration between the MRI of the atlas and test patient. The two conditional PDFs containing intensity and geometry information are combined into a unifying posterior PDF, whose mean value corresponds to the optimal electron density value under the mean-square error criterion. We evaluated the algorithm's accuracy of electron density mapping and its ability to detect bone in the head for eight patients, using an additional patient as the atlas or template. Mean absolute HU error between the estimated and true CT, as well as receiver operating characteristics for bone detection (HU > 200) were calculated. The performance was compared with a global intensity approach based on T1 and no density correction (set whole head to water). The proposed technique significantly reduced the errors in electron density estimation, with a mean absolute HU error of 126, compared with 139 for deformable registration (p = 2 x 10(-4)), 283 for the intensity approach (p = 2 x 10(-6)) and 282 without density correction (p = 5 x 10(-6)). For 90% sensitivity in bone detection, the proposed method achieved a specificity of 86%, compared with 80, 11 and 10% using deformable registration, intensity and without density correction, respectively. Notably, the Bayesian approach was more robust against anatomical differences between patients, with a specificity of 62% in the worst case (patient), compared to 30% specificity in registration-based approach. In conclusion, the proposed unifying Bayesian method provides accurate electron density estimation and bone detection from MRI of the head with highly heterogeneous anatomy.