Robust multishot diffusion-weighted imaging of the abdomen with region-based shot rejection.

PURPOSE: Diffusion-weighted (DW) imaging provides a useful clinical contrast, but is susceptible to motion-induced dephasing caused by the application of strong diffusion gradients. Phase navigators are commonly used to resolve shot-to-shot motion-induced phase in multishot reconstructions, but poor phase estimates result in signal dropout and Apparent Diffusion Coefficient (ADC) overestimation. These artifacts are prominent in the abdomen, a region prone to involuntary cardiac and respiratory motion. To improve the robustness of DW imaging in the abdomen, region-based shot rejection schemes that selectively weight regions where the shot-to-shot phase is poorly estimated were evaluated. METHODS: Spatially varying weights for each shot, reflecting both the accuracy of the estimated phase and the degree of subvoxel dephasing, were estimated from the phase navigator magnitude images. The weighting was integrated into a multishot reconstruction using different formulations and phase navigator resolutions and tested with different phase navigator resolutions in multishot DW-echo Planar Imaging acquisitions of the liver and pancreas, using conventional monopolar and velocity-compensated diffusion encoding. Reconstructed images and ADC estimates were compared qualitatively. RESULTS: The proposed region-based shot rejection reduces banding and signal dropout artifacts caused by physiological motion in the liver and pancreas. Shot rejection allows conventional monopolar diffusion encoding to achieve median ADCs in the pancreas comparable to motion-compensated diffusion encoding, albeit with a greater spread of ADCs. CONCLUSION: Region-based shot rejection is a linear reconstruction that improves the motion robustness of multi-shot DWI and requires no sequence modifications.

Comprehensive assessment of nonuniform image quality: Application to imaging near metal.

PURPOSE: Comprehensive assessment of image quality requires accounting for spatial variations in (i) intensity artifact, (ii) geometric distortion, (iii) signal-to-noise ratio (SNR), and (iv) spatial resolution, among other factors. This work presents an ensemble of methods to meet this need, from phantom design to image analysis, and applies it to the scenario of imaging near metal. METHODS: A modular phantom design employing a gyroid lattice is developed to enable the co-registered volumetric quantitation of image quality near a metallic hip implant. A method for measuring spatial resolution by means of local point spread function (PSF) estimation is presented and the relative fitness of gyroid and cubic lattices is examined. Intensity artifact, geometric distortion, and SNR maps are also computed. Results are demonstrated with 2D-FSE and MAVRIC-SL scan protocols on a 3T MRI scanner. RESULTS: The spatial resolution method demonstrates a worst-case error of 0.17 pixels for measuring in-plane blurring up to 3 pixels (full width at half maximum). The gyroid outperforms a cubic lattice design for the local PSF estimation task. The phantom supports four configurations toggling the presence/absence of both metal and structure with good spatial correspondence for co-registered analysis of the four quality factors. The marginal scan time to evaluate one scan protocol amounts to five repetitions. The phantom design can be fabricated in 2 days at negligible material cost. CONCLUSION: The phantom and associated analysis methods can elucidate complex image quality trade-offs involving intensity artifact, geometric distortion, SNR, and spatial resolution. The ensemble of methods is suitable for benchmarking imaging performance near metal.

Distortion-free water-fat separated diffusion-weighted imaging using spatiotemporal joint reconstruction.

PURPOSE: Diffusion-weighted imaging (DWI) suffers from geometric distortion and chemical shift artifacts due to the commonly used Echo Planar Imaging (EPI) trajectory. Even with fat suppression in DWI, severe B0 and B1 variations can result in residual fat, which becomes both a source of image artifacts and a confounding factor in diffusion-weighted contrast in distinguishing benign and malignant tissues. This work presents a method for acquiring distortion-free diffusion-weighted images using spatiotemporal acquisition and joint reconstruction. Water-fat separation is performed by chemical-shift encoding. METHODS: Spatiotemporal acquisition is employed to obtain distortion-free images at a series of echo times. Chemical-shift encoding is used for water-fat separation. Reconstruction and separation are performed jointly in the spat-spectral domain. To address the shot-to-shot motion-induced phase in DWI, an Fast Spin Echo (FSE)-based phase navigator is incorporated into the sequence to obtain distortion-free phase information. The proposed method was validated in phantoms and in vivo for the brain, head and neck, and breast. RESULTS: The proposed method enables the acquisition of distortion-free diffusion-weighted images in the presence of B0 field inhomogenieties commonly observed in the body. Water and fat components are separated with no obvious spectral leakage artifacts. The estimated Apparent Diffusion Coefficient (ADC) is comparable to that of multishot DW-EPI. CONCLUSION: Distortion-free, water-fat separated diffusion-weighted images in body can be obtained through the utilization of spatiotemporal acquisition and joint reconstruction methods.

A diffusion-prepared reduced FOV sequence for prostate MRI near metallic implants.

PURPOSE: To enable diffusion weighted imaging in prostate patients with metallic total hip replacements in clinically feasible scan times for prostate cancer screening, and avoid distortion and dropout artifacts present in the conventionally used Echo Planar Imaging (EPI). METHODS: A reduced field of view (FOV) diffusion-prepared sequence that is robust to the B 0 $$ {\kern0em }_0 $$ inhomogeneities produced by total hip replacements was achieved using high radiofrequency (RF) bandwidth pulses and manipulation for stimulated echo pathways. The reduced FOV along the A/P direction was obtained using slice-select gradient reversal, and the prepared magnetization was imaged with a three-dimensional RF-spoiled gradient echo readout. The sequence was validated in phantom experiments, in vivo in healthy volunteers with and without total hip replacements, and in vivo in patients undergoing a standard MRI prostate exam. RESULTS: The proposed sequence is robust to shading and distortion artifacts that are encountered by standard diffusion-weighted EPI in the presence of moderate off-resonance. Apparent diffusion coefficient estimates obtained by the proposed sequence were comparable to those obtained with diffusion-weighted EPI. CONCLUSION: Acquisition of distortionless diffusion weighted images of the prostate is feasible in patients with total hip replacements on conventional, whole-body 3T MRI, using a b-value of 800 s / mm 2 $$ \mathrm{s}/{\mathrm{mm}}^2 $$ and nominal resolution of 1.7 × $$ \times $$ 1.7 × $$ \times $$ 4 mm3 in scan times of 6 min.

Distortionless, free-breathing, and respiratory resolved 3D diffusion weighted imaging of the abdomen.

PURPOSE: Abdominal imaging is frequently performed with breath holds or respiratory triggering to reduce the effects of respiratory motion. Diffusion weighted sequences provide a useful clinical contrast but have prolonged scan times due to low signal-to-noise ratio (SNR), and cannot be completed in a single breath hold. Echo-planar imaging (EPI) is the most commonly used trajectory for diffusion weighted imaging but it is susceptible to off-resonance artifacts. A respiratory resolved, three-dimensional (3D) diffusion prepared sequence that obtains distortionless diffusion weighted images during free-breathing is presented. Techniques to address the myriad of challenges including: 3D shot-to-shot phase correction, respiratory binning, diffusion encoding during free-breathing, and robustness to off-resonance are described. METHODS: A twice-refocused, M1-nulled diffusion preparation was combined with an RF-spoiled gradient echo readout and respiratory resolved reconstruction to obtain free-breathing diffusion weighted images in the abdomen. Cartesian sampling permits a sampling density that enables 3D shot-to-shot phase navigation and reduction of transient fat artifacts. Theoretical properties of a region-based shot rejection are described. The region-based shot rejection method was evaluated with free-breathing (normal and exaggerated breathing), and respiratory triggering. The proposed sequence was compared in vivo with multishot DW-EPI. RESULTS: The proposed sequence exhibits no evident distortion in vivo when compared to multishot DW-EPI, robustness to B0 and B1 field inhomogeneities, and robustness to motion from different respiratory patterns. CONCLUSION: Acquisition of distortionless, diffusion weighted images is feasible during free-breathing with a b-value of 500 s/mm2, scan time of 6 min, and a clinically viable reconstruction time.

Reproducibility of Quantitative Double-Echo Steady-State T2 Mapping of Knee Cartilage.

BACKGROUND: Cartilage T2 can detect joints at risk of developing osteoarthritis. The quantitative double-echo steady state (qDESS) sequence is attractive for knee cartilage T2 mapping because of its acquisition time of under 5 minutes. Understanding the reproducibility errors associated with qDESS T2 is essential to profiling the technical performance of this biomarker. PURPOSE: To examine the combined acquisition and segmentation reproducibility of knee cartilage qDESS T2 using two different regional analysis schemes: 1) manual segmentation of subregions loaded during common activities and 2) automatic subregional segmentation. STUDY TYPE: Prospective. SUBJECTS: 11 uninjured participants (age: 28 ± 3 years; 8 (73%) female). FIELD STRENGTH/SEQUENCE: 3-T, qDESS. ASSESSMENT: Test-retest T2 maps were acquired twice on the same day and with a 1-week interval between scans. For each acquisition, average cartilage T2 was calculated in four manually segmented regions encompassing tibiofemoral contact areas during common activities and 12 automatically segmented regions from the deep-learning open-source framework for musculoskeletal MRI analysis (DOSMA) encompassing medial and lateral anterior, central, and posterior tibiofemoral regions. Test-retest T2 values from matching regions were used to evaluate reproducibility. STATISTICAL TESTS: Coefficients of variation (%CV), root-mean-square-average-CV (%RMSA-CV), and intraclass correlation coefficients (ICCs) assessed test-retest T2 reproducibility. The median of test-retest standard deviations was used for T2 precision. Bland-Altman (BA) analyses examined test-retest biases. The smallest detectable difference (SDD) was defined as the BA limit of agreement of largest magnitude. Significance was accepted for P < 0.05. RESULTS: All cartilage regions across both segmentation schemes demonstrated intraday and interday qDESS T2 CVs and RMSA-CVs of ≤5%. T2 ICC values >0.75 were observed in the majority of regions but were more variable in interday tibial comparisons. Test-retest T2 precision was <1.3 msec. The T2 SDD was 3.8 msec. DATA CONCLUSION: Excellent CV and RMSA-CV reproducibility may suggest that qDESS T2 increases or decreases >5% (3.8 msec) could represent changes to cartilage composition. TECHNICAL EFFICACY: Stage 2.

A method for measuring B0 field inhomogeneity using quantitative double-echo in steady-state.

PURPOSE: To develop and validate a method for B 0 $$ {B}_0 $$ mapping for knee imaging using the quantitative Double-Echo in Steady-State (qDESS) exploiting the phase difference ( Δ Î¸ $$ \Delta \theta $$ ) between the two echoes acquired. Contrary to a two-gradient-echo (2-GRE) method, Δ Î¸ $$ \Delta \theta $$ depends only on the first echo time. METHODS: Bloch simulations were applied to investigate robustness to noise of the proposed methodology and all imaging studies were validated with phantoms and in vivo simultaneous bilateral knee acquisitions. Two phantoms and five healthy subjects were scanned using qDESS, water saturation shift referencing (WASSR), and multi-GRE sequences. Δ B 0 $$ \Delta {B}_0 $$ maps were calculated with the qDESS and the 2-GRE methods and compared against those obtained with WASSR. The comparison was quantitatively assessed exploiting pixel-wise difference maps, Bland-Altman (BA) analysis, and Lin's concordance coefficient ( ρ c $$ {\rho}_c $$ ). For in vivo subjects, the comparison was assessed in cartilage using average values in six subregions. RESULTS: The proposed method for measuring Δ B 0 $$ \Delta {B}_0 $$ inhomogeneities from a qDESS acquisition provided Δ B 0 $$ \Delta {B}_0 $$ maps that were in good agreement with those obtained using WASSR. Δ B 0 $$ \Delta {B}_0 $$ ρ c $$ {\rho}_c $$ values were ≥ $$ \ge $$ 0.98 and 0.90 in phantoms and in vivo, respectively. The agreement between qDESS and WASSR was comparable to that of a 2-GRE method. CONCLUSION: The proposed method may allow B0 correction for qDESS T 2 $$ {T}_2 $$ mapping using an inherently co-registered Δ B 0 $$ \Delta {B}_0 $$ map without requiring an additional B0 measurement sequence. More generally, the method may help shorten knee imaging protocols that require an auxiliary Δ B 0 $$ \Delta {B}_0 $$ map by exploiting a qDESS acquisition that also provides T 2 $$ {T}_2 $$ measurements and high-quality morphological imaging.

Noise2Recon: Enabling SNR-robust MRI reconstruction with semi-supervised and self-supervised learning.

PURPOSE: To develop a method for building MRI reconstruction neural networks robust to changes in signal-to-noise ratio (SNR) and trainable with a limited number of fully sampled scans. METHODS: We propose Noise2Recon, a consistency training method for SNR-robust accelerated MRI reconstruction that can use both fully sampled (labeled) and undersampled (unlabeled) scans. Noise2Recon uses unlabeled data by enforcing consistency between model reconstructions of undersampled scans and their noise-augmented counterparts. Noise2Recon was compared to compressed sensing and both supervised and self-supervised deep learning baselines. Experiments were conducted using retrospectively accelerated data from the mridata three-dimensional fast-spin-echo knee and two-dimensional fastMRI brain datasets. All methods were evaluated in label-limited settings and among out-of-distribution (OOD) shifts, including changes in SNR, acceleration factors, and datasets. An extensive ablation study was conducted to characterize the sensitivity of Noise2Recon to hyperparameter choices. RESULTS: In label-limited settings, Noise2Recon achieved better structural similarity, peak signal-to-noise ratio, and normalized-RMS error than all baselines and matched performance of supervised models, which were trained with 14 × $$ 14\times $$ more fully sampled scans. Noise2Recon outperformed all baselines, including state-of-the-art fine-tuning and augmentation techniques, among low-SNR scans and when generalizing to OOD acceleration factors. Augmentation extent and loss weighting hyperparameters had negligible impact on Noise2Recon compared to supervised methods, which may indicate increased training stability. CONCLUSION: Noise2Recon is a label-efficient reconstruction method that is robust to distribution shifts, such as changes in SNR, acceleration factors, and others, with limited or no fully sampled training data.

Multishot Diffusion-Weighted MRI of the Breasts in the Supine vs. Prone Position.

BACKGROUND: Diffusion-weighted imaging (DWI) may allow for breast cancer screening MRI without a contrast injection. Multishot methods improve prone DWI of the breasts but face different challenges in the supine position. PURPOSE: To establish a multishot DWI (msDWI) protocol for supine breast MRI and to evaluate the performance of supine vs. prone msDWI. STUDY TYPE: Prospective. POPULATION: Protocol optimization: 10 healthy women (ages 22-56), supine vs. prone: 24 healthy women (ages 22-62) and five women (ages 29-61) with breast tumors. FIELD STRENGTH/SEQUENCE: 3-T, protocol optimization msDWI: free-breathing (FB) 2-shots, FB 4-shots, respiratory-triggered (RT) 2-shots, RT 4-shots, supine vs. prone: RT 4-shot msDWI, T2-weighted fast-spin echo. ASSESSMENT: Protocol optimization and supine vs. prone: three observers performed an image quality assessment of sharpness, aliasing, distortion (vs. T2), perceived SNR, and overall image quality (scale of 1-5). Apparent diffusion coefficients (ADCs) in fibroglandular tissue (FGT) and breast tumors were measured. STATISTICAL TESTS: Effect of study variables on dichotomized ratings (4/5 vs. 1/2/3) and FGT ADCs were assessed with mixed-effects logistic regression. Interobserver agreement utilized Gwet's agreement coefficient (AC). Lesion ADCs were assessed by Bland-Altman analysis and concordance correlation (ρc ). P value <0.05 was considered statistically significant. RESULTS: Protocol optimization: 4-shots significantly improved sharpness and distortion; RT significantly improved sharpness, aliasing, perceived SNR, and overall image quality. FGT ADCs were not significantly different between shots (P = 0.812), FB vs. RT (P = 0.591), or side (P = 0.574). Supine vs. prone: supine images were rated significantly higher for sharpness, aliasing, and overall image quality. FGT ADCs were significantly higher supine; lesion ADCs were highly correlated (ρc  = 0.92). DATA CONCLUSION: Based on image quality, supine msDWI outperformed prone msDWI. Lesion ADCs were highly correlated between the two positions, while FGT ADCs were higher in the supine position. EVIDENCE LEVEL: 2. TECHNICAL EFFICACY: Stage 1.

Generalizability of Deep Learning Segmentation Algorithms for Automated Assessment of Cartilage Morphology and MRI Relaxometry.

BACKGROUND: Deep learning (DL)-based automatic segmentation models can expedite manual segmentation yet require resource-intensive fine-tuning before deployment on new datasets. The generalizability of DL methods to new datasets without fine-tuning is not well characterized. PURPOSE: Evaluate the generalizability of DL-based models by deploying pretrained models on independent datasets varying by MR scanner, acquisition parameters, and subject population. STUDY TYPE: Retrospective based on prospectively acquired data. POPULATION: Overall test dataset: 59 subjects (26 females); Study 1: 5 healthy subjects (zero females), Study 2: 8 healthy subjects (eight females), Study 3: 10 subjects with osteoarthritis (eight females), Study 4: 36 subjects with various knee pathology (10 females). FIELD STRENGTH/SEQUENCE: A 3-T, quantitative double-echo steady state (qDESS). ASSESSMENT: Four annotators manually segmented knee cartilage. Each reader segmented one of four qDESS datasets in the test dataset. Two DL models, one trained on qDESS data and another on Osteoarthritis Initiative (OAI)-DESS data, were assessed. Manual and automatic segmentations were compared by quantifying variations in segmentation accuracy, volume, and T2 relaxation times for superficial and deep cartilage. STATISTICAL TESTS: Dice similarity coefficient (DSC) for segmentation accuracy. Lin's concordance correlation coefficient (CCC), Wilcoxon rank-sum tests, root-mean-squared error-coefficient-of-variation to quantify manual vs. automatic T2 and volume variations. Bland-Altman plots for manual vs. automatic T2 agreement. A P value < 0.05 was considered statistically significant. RESULTS: DSCs for the qDESS-trained model, 0.79-0.93, were higher than those for the OAI-DESS-trained model, 0.59-0.79. T2 and volume CCCs for the qDESS-trained model, 0.75-0.98 and 0.47-0.95, were higher than respective CCCs for the OAI-DESS-trained model, 0.35-0.90 and 0.13-0.84. Bland-Altman 95% limits of agreement for superficial and deep cartilage T2 were lower for the qDESS-trained model, ±2.4 msec and ±4.0 msec, than the OAI-DESS-trained model, ±4.4 msec and ±5.2 msec. DATA CONCLUSION: The qDESS-trained model may generalize well to independent qDESS datasets regardless of MR scanner, acquisition parameters, and subject population. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 1.