A multi-scale residual network for accelerated radial MR parameter mapping

The Role of Dynamic Contrast Enhanced Magnetic Resonance Imaging in Evaluating Prostate Adenocarcinoma: A Partially-Blinded Retrospective Study of a Prostatectomy Patient Cohort With Whole Gland Histopathology Correlation and Application of PI-RADS or TNM Staging

BACKGROUND

Dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI) in the current Prostate Imaging-Reporting and Data System version 2.1 (PI-RADS v2.1) is considered optional, with primary scoring based on T2-weighted imaging (T2WI) and diffusion weighted imaging (DWI). Our study is designed to assess the relative contribution of DCE MRI in a patient-cohort with whole mount prostate histopathology and spatially-mapped prostate adenocarcinoma (PCa) for reference.

METHODS

We performed a partially-blinded retrospective review of 47 prostatectomy patients with recent multi-parametric MRI (mpMRI). Scans included T2WI, DWI with apparent diffusion coefficient (ADC) mapping, and DCE imaging. Lesion conspicuity was scored on a 10-point scale with ≥ 6 considered "positive," and image quality was assessed on a 4-point scale for each sequence. The diagnostic contribution of DCE images was evaluated on a 4-point scale. The mpMRI studies were assigned PI-RADS scores and tumor, node, metastasis (TNM) T-stage with blinded comparison to spatially-mapped whole-mount pathology. Results were compared to the prospective clinical reports, which used standardized PI-RADS templates that emphasize T2WI, DWI and ADC.

RESULTS

Per lesion sensitivity for PCa was 93.5%, 82.6%, 63.0%, and 58.7% on T2WI, DCE, ADC and DWI, respectively. Mean lesion conspicuity was 8.5, 7.9, 6.2, and 6.1, on T2W, DCE, ADC and DWI, respectively. The higher values on T2WI and DCE imaging were not significantly different from each other but were both significantly different from DWI and ADC (p < 0.001). DCE scans were determined to have a marked diagnostic contribution in 83% of patients, with the most common diagnostic yield being detection of contralateral peripheral zone tumor or delineating presence/absence of extra-prostatic extension (EPE), contributing to more accurate PCa staging by PI-RADS or TNM, as compared to histopathology.

CONCLUSION

We demonstrate that DCE may contribute to lesion detection and local staging as compared to T2WI plus DWI-ADC alone and that lesion conspicuity using DCE is markedly improved as compared to DWI-ADC. These findings support modification of PI-RADS v2.1 to include use of DCE acquisitions and that a TNM staging is feasible on mpMRI as compared to surgical pathology.

HDNLS: Hybrid Deep-Learning and Non-Linear Least Squares-Based Method for Fast Multi-Component T1ρ Mapping in the Knee Joint

Non-linear least squares (NLS) methods are commonly used for quantitative magnetic resonance imaging (MRI), especially for multi-exponential T1ρ mapping, which provides precise parameter estimation for different relaxation models in tissues, such as mono-exponential (ME), bi-exponential (BE), and stretched-exponential (SE) models. However, NLS may suffer from problems like sensitivity to initial guesses, slow convergence speed, and high computational cost. While deep learning (DL)-based T1ρ fitting methods offer faster alternatives, they often face challenges such as noise sensitivity and reliance on NLS-generated reference data for training. To address these limitations of both approaches, we propose the HDNLS, a hybrid model for fast multi-component parameter mapping, particularly targeted for T1ρ mapping in the knee joint. HDNLS combines voxel-wise DL, trained with synthetic data, with a few iterations of NLS to accelerate the fitting process, thus eliminating the need for reference MRI data for training. Due to the inverse-problem nature of the parameter mapping, certain parameters in a specific model may be more sensitive to noise, such as the short component in the BE model. To address this, the number of NLS iterations in HDNLS can act as a regularization, stabilizing the estimation to obtain meaningful solutions. Thus, in this work, we conducted a comprehensive analysis of the impact of NLS iterations on HDNLS performance and proposed four variants that balance estimation accuracy and computational speed. These variants are Ultrafast-NLS, Superfast-HDNLS, HDNLS, and Relaxed-HDNLS. These methods allow users to select a suitable configuration based on their specific speed and performance requirements. Among these, HDNLS emerges as the optimal trade-off between performance and fitting time. Extensive experiments on synthetic data demonstrate that HDNLS achieves comparable performance to NLS and regularized-NLS (RNLS) with a minimum of a 13-fold improvement in speed. HDNLS is just a little slower than DL-based methods; however, it significantly improves estimation quality, offering a solution for T1ρ fitting that is fast and reliable.

Denoising and uncertainty estimation in parameter mapping with approximate Bayesian deep image priors

PURPOSE

To mitigate the problem of noisy parameter maps with high uncertainties by casting parameter mapping as a denoising task based on Deep Image Priors.

METHODS

We extend the concept of denoising with Deep Image Prior (DIP) into parameter mapping by treating the output of an image-generating network as a parametrization of tissue parameter maps. The method implicitly denoises the parameter mapping process by filtering low-level image features with an untrained convolutional neural network (CNN). Our implementation includes uncertainty estimation from Bernoulli approximate variational inference, implemented with MC dropout, which provides model uncertainty in each voxel of the denoised parameter maps. The method is modular, so the specifics of different applications (e.g., T1 mapping) separate into application-specific signal equation blocks. We evaluate the method on variable flip angle T1 mapping, multi-echo T2 mapping, and apparent diffusion coefficient mapping.

RESULTS

We found that deep image prior adapts successfully to several applications in parameter mapping. In all evaluations, the method produces noise-reduced parameter maps with decreased uncertainty compared to conventional methods. The downsides of the proposed method are the long computational time and the introduction of some bias from the denoising prior.

CONCLUSION

DIP successfully denoise the parameter mapping process and applies to several applications with limited hyperparameter tuning. Further, it is easy to implement since DIP methods do not use network training data. Although time-consuming, uncertainty information from MC dropout makes the method more robust and provides useful information when properly calibrated.

Uncertainty-aware physics-driven deep learning network for free-breathing liver fat and R2 * quantification using self-gated stack-of-radial MRI

PURPOSE

To develop a deep learning-based method for rapid liver proton-density fat fraction (PDFF) and R2 * quantification with built-in uncertainty estimation using self-gated free-breathing stack-of-radial MRI.

METHODS

This work developed an uncertainty-aware physics-driven deep learning network (UP-Net) to (1) suppress radial streaking artifacts because of undersampling after self-gating, (2) calculate accurate quantitative maps, and (3) provide pixel-wise uncertainty maps. UP-Net incorporated a phase augmentation strategy, generative adversarial network architecture, and an MRI physics loss term based on a fat-water and R2 * signal model. UP-Net was trained and tested using free-breathing multi-echo stack-of-radial MRI data from 105 subjects. UP-Net uncertainty scores were calibrated in a validation dataset and used to predict quantification errors for liver PDFF and R2 * in a testing dataset.

RESULTS

Compared with images reconstructed using compressed sensing (CS), UP-Net achieved structural similarity index >0.87 and normalized root mean squared error <0.18. Compared with reference quantitative maps generated using CS and graph-cut (GC) algorithms, UP-Net achieved low mean differences (MD) for liver PDFF (-0.36%) and R2 * (-0.37 s-1 ). Compared with breath-holding Cartesian MRI results, UP-Net achieved low MD for liver PDFF (0.53%) and R2 * (6.75 s-1 ). UP-Net uncertainty scores predicted absolute liver PDFF and R2 * errors with low MD of 0.27% and 0.12 s-1 compared to CS + GC results. The computational time for UP-Net was 79 ms/slice, whereas CS + GC required 3.2 min/slice.

CONCLUSION

UP-Net rapidly calculates accurate liver PDFF and R2 * maps from self-gated free-breathing stack-of-radial MRI. The pixel-wise uncertainty maps from UP-Net predict quantification errors in the liver.

SuperMAP: Deep ultrafast MR relaxometry with joint spatiotemporal undersampling

PURPOSE

To develop an ultrafast and robust MR parameter mapping network using deep learning.

THEORY AND METHODS

We design a deep learning framework called SuperMAP that directly converts a series of undersampled (both in k-space and parameter-space) parameter-weighted images into several quantitative maps, bypassing the conventional exponential fitting procedure. We also present a novel technique to simultaneously reconstruct T1rho and T2 relaxation maps within a single scan. Full data were acquired and retrospectively undersampled for training and testing using traditional and state-of-the-art techniques for comparison. Prospective data were also collected to evaluate the trained network. The performance of all methods is evaluated using the parameter qualification errors and other metrics in the segmented regions of interest.

RESULTS

SuperMAP achieved accurate T1rho and T2 mapping with high acceleration factors (R = 24 and R = 32). It exploited both spatial and temporal information and yielded low error (normalized mean square error of 2.7% at R = 24 and 2.8% at R = 32) and high resemblance (structural similarity of 97% at R = 24 and 96% at R = 32) to the gold standard. The network trained with retrospectively undersampled data also works well for the prospective data (with a slightly lower acceleration factor). SuperMAP is also superior to conventional methods.

CONCLUSION

Our results demonstrate the feasibility of generating superfast MR parameter maps through very few undersampled parameter-weighted images. SuperMAP can simultaneously generate T1rho and T2 relaxation maps in a short scan time.

Highly accelerated MR parametric mapping by undersampling the k-space and reducing the contrast number simultaneously with deep learning

Introduction. To propose a novel deep learning-based method called RG-Net (reconstruction and generation network) for highly accelerated MR parametric mapping by undersampling k-space and reducing the acquired contrast number simultaneously. Methods. The proposed framework consists of a reconstruction module and a generative module. The reconstruction module reconstructs MR images from the acquired few undersampled k-space data with the help of a data prior. The generative module then synthesizes the remaining multi-contrast images from the reconstructed images, where the exponential model is implicitly incorporated into the image generation through the supervision of fully sampled labels. The RG-Net was trained and tested on the T1ρ mapping data from 8 volunteers at net acceleration rates of 17, respectively. Regional T1ρ analysis for cartilage and the brain was performed to assess the performance of RG-Net. Results. RG-Net yields a high-quality T1ρ map at a high acceleration rate of 17. Compared with the competing methods that only undersample k-space, our framework achieves better performance in T1ρ value analysis. Conclusion. The proposed RG-Net can achieve a high acceleration rate while maintaining good reconstruction quality by undersampling k-space and reducing the contrast number simultaneously for fast MR parametric mapping. The generative module of our framework can also be used as an insertable module in other fast MR parametric mapping methods.

Convolutional Neural Network-Based Computer-Assisted Diagnosis of Hashimoto's Thyroiditis on Ultrasound

PURPOSE

This study investigates the efficiency of deep learning models in the automated diagnosis of Hashimoto's thyroiditis (HT) using real-world ultrasound data from ultrasound examinations by computer-assisted diagnosis (CAD) with artificial intelligence.

METHODS

We retrospectively collected ultrasound images from patients with and without HT from 2 hospitals in China between September 2008 and February 2018. Images were divided into a training set (80%) and a validation set (20%). We ensembled 9 convolutional neural networks (CNNs) as the final model (CAD-HT) for HT classification. The model's diagnostic performance was validated and compared to 2 hospital validation sets. We also compared the accuracy of CAD-HT against seniors/junior radiologists. Subgroup analysis of CAD-HT performance for different thyroid hormone levels (hyperthyroidism, hypothyroidism, and euthyroidism) was also evaluated.

RESULTS

39 280 ultrasound images from 21 118 patients were included in this study. The accuracy, sensitivity, and specificity of the HT-CAD model were 0.892, 0.890, and 0.895, respectively. HT-CAD performance between 2 hospitals was not significantly different. The HT-CAD model achieved a higher performance (P < 0.001) when compared to senior radiologists, with a nearly 9% accuracy improvement. HT-CAD had almost similar accuracy (range 0.87-0.894) for the 3 subgroups based on thyroid hormone level.

CONCLUSION

The HT-CAD strategy based on CNN significantly improved the radiologists' diagnostic accuracy of HT. Our model demonstrates good performance and robustness in different hospitals and for different thyroid hormone levels.

Abdominal T2-Weighted Imaging and T2 Mapping Using a Variable Flip Angle Radial Turbo Spin-Echo Technique

BACKGROUND

T2 mapping is of great interest in abdominal imaging but current methods are limited by low resolution, slice coverage, motion sensitivity, or lengthy acquisitions.

PURPOSE

Develop a radial turbo spin-echo technique with refocusing variable flip angles (RADTSE-VFA) for high spatiotemporal T2 mapping and efficient slice coverage within a breath-hold and compare to the constant flip angle counterpart (RADTSE-CFA).

STUDY TYPE

Prospective technical efficacy.

SUBJECTS

Testing performed on agarose phantoms and 12 patients. Focal liver lesion classification tested on malignant (N = 24) and benign (N = 11) lesions.

FIELD STRENGTH/SEQUENCE

1.5 T/RADTSE-VFA, RADTSE-CFA.

ASSESSMENT

A constrained objective function was used to optimize the refocusing flip angles. Phantom and/or in vivo data were used to assess relative contrast, T2 estimation, specific absorption rate (SAR), and focal liver lesion classification. STATISTICAL TESTS: t-Tests or Mann-Whitney Rank Sum tests were used.

RESULTS

Phantom data did not show significant differences in mean relative contrast (P = 0.10) and T2 accuracy (P = 0.99) between RADTSE-VFA and RADTSE-CFA. Adding noise caused T2 overestimation predominantly for RADTSE-CFA and low T2 values. In vivo results did not show significant differences in mean spleen-to-liver (P = 0.62) and kidney-to-liver (P = 0.49) relative contrast between RADTSE-VFA and RADTSE-CFA. Mean T2 values were not significantly different between the two techniques for spleen (T2VFA  = 109.2 ± 12.3 msec; T2CFA  = 110.7 ± 11.1 msec; P = 0.78) and kidney-medulla (T2VFA  = 113.0 ± 8.7 msec; T2CFA  = 114.0 ± 8.6 msec; P = 0.79). Liver T2 was significantly higher for RADTSE-CFA (T2VFA  = 52.6 ± 6.6 msec; T2CFA  = 60.4 ± 8.0 msec) consistent with T2 overestimation in the phantom study. Focal liver lesion classification had comparable T2 distributions for RADTSE-VFA and RADTSE-CFA for malignancies (P = 1.0) and benign lesions (P = 0.39). RADTSE-VFA had significantly lower SAR than RADTSE-CFA increasing slice coverage by 1.5.

DATA CONCLUSION

RADTSE-VFA provided noise-robust T2 estimation compared to the constant flip angle counterpart while generating T2-weighted images with comparable contrast. The VFA scheme minimized SAR improving slice efficiency for breath-hold imaging.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 1.

Synergistic multi-contrast cardiac magnetic resonance image reconstruction

Cardiac magnetic resonance imaging (CMR) is an important tool for the non-invasive diagnosis of a variety of cardiovascular diseases. Parametric mapping with multi-contrast CMR is able to quantify tissue alterations in myocardial disease and promises to improve patient care. However, magnetic resonance imaging is an inherently slow imaging modality, resulting in long acquisition times for parametric mapping which acquires a series of cardiac images with different contrasts for signal fitting or dictionary matching. Furthermore, extra efforts to deal with respiratory and cardiac motion by triggering and gating further increase the scan time. Several techniques have been developed to speed up CMR acquisitions, which usually acquire less data than that required by the Nyquist-Shannon sampling theorem, followed by regularized reconstruction to mitigate undersampling artefacts. Recent advances in CMR parametric mapping speed up CMR by synergistically exploiting spatial-temporal and contrast redundancies. In this article, we will review the recent developments in multi-contrast CMR image reconstruction for parametric mapping with special focus on low-rank and model-based reconstructions. Deep learning-based multi-contrast reconstruction has recently been proposed in other magnetic resonance applications. These developments will be covered to introduce the general methodology. Current technical limitations and potential future directions are discussed. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 1'.