Development and benchmarking diffusion magnetic resonance imaging analysis for integration into radiation treatment planning

Purpose: The rising promise in the utility of advanced multi-parametric magnetic resonance (MR) imaging in radiotherapy (RT) treatment planning creates a necessity for testing and enhancing the accuracy of quantitative imaging analysis. Standardizing the analysis of diffusion weighted imaging (DWI) and diffusion tensor imaging (DTI) to generate meaningful and reproducible apparent diffusion coefficient (ADC) and fractional anisotropy (FA) lays the requisite needed for clinical integration. The aim of the demonstrated work is to benchmark the generation of the ADC and FA parametric map analyses using integrated tools in a commercial treatment planning system against the currently used ones. Methods: Three software packages were used for generating ADC and FA maps in this study; one tool was developed within a commercial treatment planning system, another by the Functional Magnetic Resonance Imaging of the Brain (FMRIB) Software Library FSL (Analysis Group, FMRIB, Oxford, United Kingdom), and an in-house tool developed at the M.D. Anderson Cancer Center. The ADC and FA maps generated by all three packages for 35 subjects were subtracted from one another, and the standard deviation of the images’ differences was used to compare the reproducibility. The reproducibility of the ADC maps was compared with the Quantitative Imaging Biomarkers Alliance (QIBA) protocol, while that of the FA maps was compared to data in published literature. Results: Results show that the discrepancies between the ADC maps calculated for each patient using the three different software algorithms are less than 2% which meets the 3.6% recommended QIBA requirement. Except for a small number of isolated points, the majority of differences in FA maps for each patient produced by the three methods did not exceed 0.02 which is 10 times lower than the differences seen in healthy gray and white matter. The results were also compared to the maps generated by existing MR Imaging consoles and showed that the robustness of console generated ADC and FA maps is largely dependent on the correct application of scaling factors, that only if correctly placed; the differences between the three tested methods and the console generated values were within the recommended QIBA guidelines. Conclusions: Cross-comparison difference maps demonstrated that quantitative reproducibility of ADC and FA metrics generated using our tested commercial treatment planning system were comparable to in-house and established tools as benchmarks. This integrated approach facilitates the clinical utility of diffusion imaging in radiation treatment planning workflow.