Archana Vadiraj Malagi1, Arjunlokesh Netaji2, Virendra Kumar3, Esha Baidya Kayal1, Kedar Khare4, Chandan Jyoti Das2, Fernando Calamante5, Amit Mehndiratta6,7. 1. Centre for Biomedical Engineering, Indian Institute of Technology Delhi, Hauz Khas, New Delhi, 110016, India. 2. Department of Radio-Diagnosis, All India Institute of Medical Sciences, New Delhi, India. 3. Department of Nuclear Magnetic Resonance, All India Institute of Medical Sciences, New Delhi, India. 4. Department of Physics, Indian Institute of Technology Delhi, New Delhi, India. 5. Sydney Imaging and School of Biomedical Engineering, University of Sydney, Sydney, Australia. 6. Centre for Biomedical Engineering, Indian Institute of Technology Delhi, Hauz Khas, New Delhi, 110016, India. amit.mehndiratta@keble.oxon.org. 7. Department of Biomedical Engineering, All India Institute of Medical Sciences, New Delhi, India. amit.mehndiratta@keble.oxon.org.
Abstract
OBJECTIVE: To implement an advanced spatial penalty-based reconstruction to constrain the intravoxel incoherent motion (IVIM)-diffusion kurtosis imaging (DKI) model and investigate whether it provides a suitable alternative at 1.5 T to the traditional IVIM-DKI model at 3 T for clinical characterization of prostate cancer (PCa) and benign prostatic hyperplasia (BPH). MATERIALS AND METHODS: Thirty-two patients with biopsy-proven PCa were recruited for MRI examination (n = 16 scanned at 1.5 T, n = 16 scanned at 3 T). Diffusion-weighted imaging (DWI) with 13 b values (b = 0 to 2000 s/mm2 up to 3 averages, 1.5 T: TR = 5.774 s, TE = 81 ms and 3 T: TR = 4.899 s, TE = 100 ms), T2-weighted, and T1-weighted imaging were used on the 1.5 T and 3 T MRI scanner, respectively. The IVIM-DKI signal was modeled using the traditional IVIM-DKI model and a novel model in which the total variation (TV) penalty function was combined with the traditional model to optimize non-physiological variations. Paired and unpaired t-tests were used to compare intra-scanner and scanner group differences in IVIM-DKI parameters obtained using the novel and the traditional models. Analysis of variance with post hoc test and receiver operating characteristic (ROC) curve analysis were used to assess the ability of parameters obtained using the novel model (at 1.5 T) and the traditional model (at 3 T) to characterize prostate lesions. RESULTS: IVIM-DKI modeled using novel model with TV spatial penalty function at 1.5 T, produced parameter maps with 50-78% lower coefficient of variation (CV) than traditional model at 3 T. Novel model estimated higher D with lower D*, f and k values at both field strengths compared to traditional model. For scanner differences, the novel model at 1.5 T estimated lower D* and f values as compared to traditional model at 3 T. At 1.5 T, D and f values were significantly lower with k values significantly higher in tumor than BPH and healthy tissue. D (AUC: 0.98), f (AUC: 0.82), and k (AUC: 0.91) parameters estimated using novel model showed high diagnostic performance in cancer lesion detection at 1.5 T. DISCUSSION: In comparison with the IVIM-DKI model at 3 T, IVIM-DKI signal modeled with the TV penalty function at 1.5 T showed lower estimation errors. The proposed novel model can be utilized for improved detection of prostate lesions.
OBJECTIVE: To implement an advanced spatial penalty-based reconstruction to constrain the intravoxel incoherent motion (IVIM)-diffusion kurtosis imaging (DKI) model and investigate whether it provides a suitable alternative at 1.5 T to the traditional IVIM-DKI model at 3 T for clinical characterization of prostate cancer (PCa) and benign prostatic hyperplasia (BPH). MATERIALS AND METHODS: Thirty-two patients with biopsy-proven PCa were recruited for MRI examination (n = 16 scanned at 1.5 T, n = 16 scanned at 3 T). Diffusion-weighted imaging (DWI) with 13 b values (b = 0 to 2000 s/mm2 up to 3 averages, 1.5 T: TR = 5.774 s, TE = 81 ms and 3 T: TR = 4.899 s, TE = 100 ms), T2-weighted, and T1-weighted imaging were used on the 1.5 T and 3 T MRI scanner, respectively. The IVIM-DKI signal was modeled using the traditional IVIM-DKI model and a novel model in which the total variation (TV) penalty function was combined with the traditional model to optimize non-physiological variations. Paired and unpaired t-tests were used to compare intra-scanner and scanner group differences in IVIM-DKI parameters obtained using the novel and the traditional models. Analysis of variance with post hoc test and receiver operating characteristic (ROC) curve analysis were used to assess the ability of parameters obtained using the novel model (at 1.5 T) and the traditional model (at 3 T) to characterize prostate lesions. RESULTS: IVIM-DKI modeled using novel model with TV spatial penalty function at 1.5 T, produced parameter maps with 50-78% lower coefficient of variation (CV) than traditional model at 3 T. Novel model estimated higher D with lower D*, f and k values at both field strengths compared to traditional model. For scanner differences, the novel model at 1.5 T estimated lower D* and f values as compared to traditional model at 3 T. At 1.5 T, D and f values were significantly lower with k values significantly higher in tumor than BPH and healthy tissue. D (AUC: 0.98), f (AUC: 0.82), and k (AUC: 0.91) parameters estimated using novel model showed high diagnostic performance in cancer lesion detection at 1.5 T. DISCUSSION: In comparison with the IVIM-DKI model at 3 T, IVIM-DKI signal modeled with the TV penalty function at 1.5 T showed lower estimation errors. The proposed novel model can be utilized for improved detection of prostate lesions.