ShapeMed-Knee: A Dataset and Neural Shape Model Benchmark for Modeling 3D Femurs.

Analyzing anatomic shapes of tissues and organs is pivotal for accurate disease diagnostics and clinical decision-making. One prominent disease that depends on anatomic shape analysis is osteoarthritis, which affects 30 million Americans. To advance osteoarthritis diagnostics and prognostics, we introduce ShapeMed-Knee, a 3D shape dataset with 9,376 high-resolution, medical-imaging-based 3D shapes of both femur bone and cartilage. Besides data, ShapeMed-Knee includes two benchmarks for assessing reconstruction accuracy and five clinical prediction tasks that assess the utility of learned shape representations. Leveraging ShapeMed-Knee, we develop and evaluate a novel hybrid explicit-implicit neural shape model which achieves up to 40% better reconstruction accuracy than a statistical shape model and implicit neural shape model. Our hybrid models achieve state-of-the-art performance for preserving cartilage biomarkers; they're also the first models to successfully predict localized structural features of osteoarthritis, outperforming shape models and convolutional neural networks applied to raw magnetic resonance images and segmentations. The ShapeMed-Knee dataset provides medical evaluations to reconstruct multiple anatomic surfaces and embed meaningful disease-specific information. ShapeMed-Knee reduces barriers to applying 3D modeling in medicine, and our benchmarks highlight that advancements in 3D modeling can enhance the diagnosis and risk stratification for complex diseases. The dataset, code, and benchmarks will be made freely accessible.

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.

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.

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.

Clinical evaluation of isotropic MAVRIC-SL for symptomatic hip arthroplasties at 3 T MRI.

BACKGROUND: 3D multi-spectral imaging (MSI) of metal implants necessitates relatively long scan times. OBJECTIVE: We implemented a fast isotropic 3D MSI technique at 3 T and compared its image quality and clinical utility to non-isotropic MSI in the evaluation of hip implants. METHODS: Two musculoskeletal radiologists scored images from coronal proton density-weighted conventional MAVRIC-SL and an isotropic MAVRIC-SL sequence accelerated with robust-component-analysis on a 3-point scale (3: diagnostic, 2: moderately diagnostic, 1: non-diagnostic) for overall image quality, metal artifact, and visualization around femoral and acetabular components. Grades were compared using a signed Wilcoxon test. Images were evaluated for effusion, synovitis, osteolysis, loosening, pseudotumor, fracture, and gluteal tendon abnormalities. Reformatted axial and sagittal images for both sequences were subsequently generated and compared for image quality with the Wilcoxon test. Whether these reformats increased diagnostic confidence or revealed additional pathology, including findings unrelated to arthroplasty that may contribute to hip pain, was also compared using the McNemar test. Inter-rater agreement was measured by Cohen's kappa. RESULTS: 39 symptomatic patients with a total of 59 hip prostheses were imaged (mean age, 70 years ±9, 14 males, 25 females). Comparison scores between coronal images showed no significant difference in image quality, metal artifact, or visualization of the femur and acetabulum. Except for loosening, reviewers identified more positive cases of pathology on the original coronally-acquired isotropic sequence. In comparison of reformatted axial and sagittal images, the isotropic sequence scored significantly (p < 0.01) higher for overall image quality (3.0 vs 2.0) and produced significantly (p < 0.01) more cases of increased diagnostic confidence (42.4% vs 7.6%) or additional diagnoses (50.8% vs 22.9%). Inter-rater agreement was substantial (k = 0.798) for image quality. Mean scan times were 4.2 mins (isotropic) and 7.1 mins (non-isotropic). CONCLUSION: Compared to the non-isotropic sequence, isotropic 3D MSI was acquired in less time while maintaining diagnostically acceptable image quality. It identified more pathology, including postoperative complications and potential pain-generating pathology unrelated to arthroplasty. This fast isotropic 3D MSI sequence demonstrates promise for improving diagnostic evaluation of symptomatic hip prostheses at 3 T while simultaneously reducing scan time.

The Impact of Occlusion on Depth Perception at Arm's Length.

This paper investigates the accuracy of Augmented Reality (AR) technologies, particularly commercially available optical see-through displays, in depicting virtual content inside the human body for surgical planning. Their inherent limitations result in inaccuracies in perceived object positioning. We examine how occlusion, specifically with opaque surfaces, affects perceived depth of virtual objects at arm's length working distances. A custom apparatus with a half-silvered mirror was developed, providing accurate depth cues excluding occlusion, differing from commercial displays. We carried out a study, contrasting our apparatus with a HoloLens 2, involving a depth estimation task under varied surface complexities and illuminations. In addition, we explored the effects of creating a virtual "hole" in the surface. Subjects' depth estimation accuracy and confidence were a ssessed. Results showed more depth estimation variation with HoloLens and significant depth error beneath complex occluding surfaces. However, creating a virtual hole significantly reduced depth errors and increased subjects' confidence, irrespective of accuracy enhancement. These findings have important implications for the design and use of mixed-reality technologies in surgical applications, and industrial applications such as using virtual content to guide maintenance or repair of components hidden beneath the opaque outer surface of equipment. A free copy of this paper and all supplemental materials are available at https://bit.ly/3YbkwjU.

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.

Stereoscopic calibration for augmented reality visualization in microscopic surgery.

PURPOSE: Middle and inner ear procedures target hearing loss, infections, and tumors of the temporal bone and lateral skull base. Despite the advances in surgical techniques, these procedures remain challenging due to limited haptic and visual feedback. Augmented reality (AR) may improve operative safety by allowing the 3D visualization of anatomical structures from preoperative computed tomography (CT) scans on real intraoperative microscope video feed. The purpose of this work was to develop a real-time CT-augmented stereo microscope system using camera calibration and electromagnetic (EM) tracking. METHODS: A 3D printed and electromagnetically tracked calibration board was used to compute the intrinsic and extrinsic parameters of the surgical stereo microscope. These parameters were used to establish a transformation between the EM tracker coordinate system and the stereo microscope image space such that any tracked 3D point can be projected onto the left and right images of the microscope video stream. This allowed the augmentation of the microscope feed of a 3D printed temporal bone with its corresponding CT-derived virtual model. Finally, the calibration board was also used for evaluating the accuracy of the calibration. RESULTS: We evaluated the accuracy of the system by calculating the registration error (RE) in 2D and 3D in a microsurgical laboratory setting. Our calibration workflow achieved a RE of 0.11 ± 0.06 mm in 2D and 0.98 ± 0.13 mm in 3D. In addition, we overlaid a 3D CT model on the microscope feed of a 3D resin printed model of a segmented temporal bone. The system exhibited small latency and good registration accuracy. CONCLUSION: We present the calibration of an electromagnetically tracked surgical stereo microscope for augmented reality visualization. The calibration method achieved accuracy within a range suitable for otologic procedures. The AR process introduces enhanced visualization of the surgical field while allowing depth perception.

[Formula: see text] Field inhomogeneity correction for qDESS [Formula: see text] mapping: application to rapid bilateral knee imaging.

PURPOSE: [Formula: see text] mapping is a powerful tool for studying osteoarthritis (OA) changes and bilateral imaging may be useful in investigating the role of between-knee asymmetry in OA onset and progression. The quantitative double-echo in steady-state (qDESS) can provide fast simultaneous bilateral knee [Formula: see text] and high-resolution morphometry for cartilage and meniscus. The qDESS uses an analytical signal model to compute [Formula: see text] relaxometry maps, which require knowledge of the flip angle (FA). In the presence of [Formula: see text] inhomogeneities, inconsistencies between the nominal and actual FA can affect the accuracy of [Formula: see text] measurements. We propose a pixel-wise [Formula: see text] correction method for qDESS [Formula: see text] mapping exploiting an auxiliary [Formula: see text] map to compute the actual FA used in the model. METHODS: The technique was validated in a phantom and in vivo with simultaneous bilateral knee imaging. [Formula: see text] measurements of femoral cartilage (FC) of both knees of six healthy participants were repeated longitudinally to investigate the association between [Formula: see text] variation and [Formula: see text]. RESULTS: The results showed that applying the [Formula: see text] correction mitigated [Formula: see text] variations that were driven by [Formula: see text] inhomogeneities. Specifically, [Formula: see text] left-right symmetry increased following the [Formula: see text] correction ([Formula: see text] = 0.74 > [Formula: see text] = 0.69). Without the [Formula: see text] correction, [Formula: see text] values showed a linear dependence with [Formula: see text]. The linear coefficient decreased using the [Formula: see text] correction (from 24.3 ± 1.6 ms to 4.1 ± 1.8) and the correlation was not statistically significant after the application of the Bonferroni correction (p value > 0.01). CONCLUSION: The study showed that [Formula: see text] correction could mitigate variations driven by the sensitivity of the qDESS [Formula: see text] mapping method to [Formula: see text], therefore, increasing the sensitivity to detect real biological changes. The proposed method may improve the robustness of bilateral qDESS [Formula: see text] mapping, allowing for an accurate and more efficient evaluation of OA pathways and pathophysiology through longitudinal and cross-sectional studies.

Improving Data-Efficiency and Robustness of Medical Imaging Segmentation Using Inpainting-Based Self-Supervised Learning.

We systematically evaluate the training methodology and efficacy of two inpainting-based pretext tasks of context prediction and context restoration for medical image segmentation using self-supervised learning (SSL). Multiple versions of self-supervised U-Net models were trained to segment MRI and CT datasets, each using a different combination of design choices and pretext tasks to determine the effect of these design choices on segmentation performance. The optimal design choices were used to train SSL models that were then compared with baseline supervised models for computing clinically-relevant metrics in label-limited scenarios. We observed that SSL pretraining with context restoration using 32 × 32 patches and Poission-disc sampling, transferring only the pretrained encoder weights, and fine-tuning immediately with an initial learning rate of 1 × 10-3 provided the most benefit over supervised learning for MRI and CT tissue segmentation accuracy (p < 0.001). For both datasets and most label-limited scenarios, scaling the size of unlabeled pretraining data resulted in improved segmentation performance. SSL models pretrained with this amount of data outperformed baseline supervised models in the computation of clinically-relevant metrics, especially when the performance of supervised learning was low. Our results demonstrate that SSL pretraining using inpainting-based pretext tasks can help increase the robustness of models in label-limited scenarios and reduce worst-case errors that occur with supervised learning.

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.

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.

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.

A joint linear reconstruction for multishot diffusion weighted non-Carr-Purcell-Meiboom-Gill fast spin echo with full signal.

PURPOSE: Diffusion weighted Fast Spin Echo (DW-FSE) is a promising approach for distortionless DW imaging that is robust to system imperfections such as eddy currents and off-resonance. Due to non-Carr-Purcell-Meiboom-Gill (CPMG) magnetization, most DW-FSE sequences discard a large fraction of the signal ( 2 - 2 × $$ \sqrt{2}-2\times $$ ), reducing signal-to-noise ratio (SNR) efficiency compared to DW-EPI. The full FSE signal can be preserved by quadratically incrementing the transmit phase of the refocusing pulses, but this method of resolving non-CPMG magnetization has only been applied to single-shot DW-FSE due to challenges associated with image reconstruction. We present a joint linear reconstruction for multishot quadratic phase increment data that addresses these challenges and corrects ghosting from both shot-to-shot phase and intrashot signal oscillations. Multishot imaging reduces T2 blur and joint reconstruction of shots improves g-factor performance. A thorough analysis on the condition number of the proposed linear system is described. METHODS: A joint multishot reconstruction is derived from the non-CPMG signal model. Multishot quadratic phase increment DW-FSE was tested in a standardized diffusion phantom and compared to single-shot DW-FSE and DW-EPI in vivo in the brain, cervical spine, and prostate. The pseudo multiple replica technique was applied to generate g-factor and SNR maps. RESULTS: The proposed joint shot reconstruction eliminates ghosting from shot-to-shot phase and intrashot oscillations. g-factor performance is improved compared to previously proposed reconstructions, permitting efficient multishot imaging. apparent diffusion coefficient estimates in phantom experiments and in vivo are comparable to those obtained with conventional methods. CONCLUSION: Multi-shot non-CPMG DW-FSE data with full signal can be jointly reconstructed using a linear model.

Volumetric and multispectral DWI near metallic implants using a non-linear phase Carr-Purcell-Meiboom-Gill diffusion preparation.

PURPOSE: DWI near metal implants has not been widely explored due to substantial challenges associated with through-slice and in-plane distortions, the increased encoding requirement of different spectral bins, and limited SNR. There is no widely adopted clinical protocol for DWI near metal since the commonly used EPI trajectory fails completely due to distortion from extreme off-resonance ranging from 2 to 20 kHz. We present a sequence that achieves DWI near metal with moderate b-values (400-500 s/mm2 ) and volumetric coverage in clinically feasible scan times. THEORY AND METHODS: Multispectral excitation with Cartesian sampling, view angle tilting, and kz phase encoding reduce in-plane and through-plane off-resonance artifacts, and Carr-Purcell-Meiboom-Gill (CPMG) spin-echo refocusing trains counteract T2* effects. The effect of random phase on the refocusing train is eliminated using a stimulated echo diffusion preparation. Root-flipped Shinnar-Le Roux refocusing pulses permits preparation of a high spectral bandwidth, which improves imaging times by reducing the number of excitations required to cover the desired spectral range. B1 sensitivity is reduced by using an excitation that satisfies the CPMG condition in the preparation. A method for ADC quantification insensitive to background gradients is presented. RESULTS: Non-linear phase refocusing pulses reduces the peak B1 by 46% which allows RF bandwidth to be doubled. Simulations and phantom experiments show that a non-linear phase CPMG pulse pair reduces B1 sensitivity. Application in vivo demonstrates complementary contrast to conventional multispectral acquisitions and improved visualization compared to DW-EPI. CONCLUSION: Volumetric and multispectral DW imaging near metal can be achieved with a 3D encoded sequence.

Circumventing the curse of dimensionality in magnetic resonance fingerprinting through a deep learning approach.

Magnetic resonance fingerprinting (MRF) is a rapidly developing approach for fast quantitative MRI. A typical drawback of dictionary-based MRF is an explosion of the dictionary size as a function of the number of reconstructed parameters, according to the "curse of dimensionality", which determines an explosion of resource requirements. Neural networks (NNs) have been proposed as a feasible alternative, but this approach is still in its infancy. In this work, we design a deep learning approach to MRF using a fully connected network (FCN). In the first part we investigate, by means of simulations, how the NN performance scales with the number of parameters to be retrieved in comparison with the standard dictionary approach. Four MRF sequences were considered: IR-FISP, bSSFP, IR-FISP-B1 , and IR-bSSFP-B1 , the latter two designed to be more specific for B1+ parameter encoding. Estimation accuracy, memory usage, and computational time required to perform the estimation task were considered to compare the scalability capabilities of the dictionary-based and the NN approaches. In the second part we study optimal training procedures by including different data augmentation and preprocessing strategies during training to achieve better accuracy and robustness to noise and undersampling artifacts. The study is conducted using the IR-FISP MRF sequence exploiting both simulations and in vivo acquisitions. Results demonstrate that the NN approach outperforms the dictionary-based approach in terms of scalability capabilities. Results also allow us to heuristically determine the optimal training strategy to make an FCN able to predict T1 , T2 , and M0 maps that are in good agreement with those obtained with the original dictionary approach. k-SVD denoising is proposed and found to be critical as a preprocessing step to handle undersampled data.

Clinical utility of accelerated MAVRIC-SL with robust-PCA compared to conventional MAVRIC-SL in evaluation of total hip arthroplasties.

OBJECTIVE: To compare the diagnostic performance of a conventional metal artifact suppression sequence MAVRIC-SL (multi-acquisition variable-resonance image combination selective) and a novel 2.6-fold faster sequence employing robust principal component analysis (RPCA), in the MR evaluation of hip implants at 3 T. MATERIALS AND METHODS: Thirty-six total hip implants in 25 patients were scanned at 3 T using a conventional MAVRIC-SL proton density-weighted sequence and an RPCA MAVRIC-SL proton density-weighted sequence. Comparison was made of image quality, geometric distortion, visualization around acetabular and femoral components, and conspicuity of abnormal imaging findings using the Wilcoxon signed-rank test and a non-inferiority test. Abnormal findings were correlated with subsequent clinical management and intraoperative findings if the patient underwent subsequent surgery. RESULTS: Mean scores for conventional MAVRIC-SL were better than RPCA MAVRIC-SL for all qualitative parameters (p < 0.05), although the probability of RPCA MAVRIC-SL being clinically useful was non-inferior to conventional MAVRIC-SL (within our accepted 10% difference, p < 0.05), except for visualization around the acetabular component. Abnormal imaging findings were seen in 25 hips, and either equally visible or visible but less conspicuous on RPCA MAVRIC-SL in 21 out of 25 cases. In 4 cases, a small joint effusion was queried on MAVRIC-SL but not RPCA MAVRIC-SL, but the presence or absence of a small effusion did not affect subsequent clinical management and patient outcome. CONCLUSION: While the overall image quality is reduced, RPCA MAVRIC-SL allows for significantly reduced scan time and maintains almost equal diagnostic performance.

Biopsy marker localization with thermo-acoustic ultrasound for lumpectomy guidance.

PURPOSE: Almost one in four lumpectomies fails to fully remove cancerous tissue from the breast, requiring reoperation. This high failure rate suggests that existing lumpectomy guidance methods are inadequate for allowing surgeons to consistently identify the proper volume of tissue for excision. Current guidance techniques either provide little information about the tumor position or require surgeons to frequently switch between making incisions and manually probing for a marker placed at the lesion site. This article explores the feasibility of thermo-acoustic ultrasound (TAUS) to enable hands-free localization of metallic biopsy markers throughout surgery, which would allow for continuous visualization of the lesion site in the breast without the interruption of surgery. In a TAUS-based localization system, microwave excitations would be transmitted into the breast, and the amplification in microwave absorption around the metallic markers would generate acoustic signals from the marker sites through the thermo-acoustic effect. Detection and ranging of these signals by multiple acoustic receivers on the breast could then enable marker localization through acoustic multilateration. METHODS: Physics simulations were used to characterize the TAUS signals generated from different markers by microwave excitations. First, electromagnetic simulations determined the spatial pattern of the amplification in microwave absorption around the markers. Then, acoustic simulations characterized the acoustic fields generated from these markers at various acoustic frequencies. TAUS-based one-dimensional (1D) ranging of two metallic markers-including a biopsy marker that is FDA-approved for clinical use-immersed in saline was also performed using a bench-top setup. To perform TAUS acquisitions, a microwave applicator was driven by 2.66 GHz microwave signals that were amplitude-modulated by chirps at the desired acoustic excitation frequencies, and the resulting TAUS signal from the markers was detected by an ultrasonic transducer. RESULTS: The simulation results show that the geometry of the marker strongly impacts the quantity and spatial pattern of both the microwave absorption around the marker and the resulting TAUS signal generated from the marker. The simulated TAUS signal maps and acoustic frequency responses also make clear that the marker geometry plays an important role in determining the overall system response. Using the bench-top setup, TAUS detection and 1D localization of the markers were successfully demonstrated for multiple different combinations of microwave applicator and metallic marker. These initial results indicate that TAUS-based localization of biopsy markers is feasible. CONCLUSIONS: Through microwave excitations and acoustic detection, TAUS can be used to localize metallic biopsy markers. With further development, TAUS opens new avenues to enable a more intuitive lumpectomy guidance system that could help to achieve better lumpectomy outcomes.

The International Workshop on Osteoarthritis Imaging Knee MRI Segmentation Challenge: A Multi-Institute Evaluation and Analysis Framework on a Standardized Dataset.

PURPOSE: To organize a multi-institute knee MRI segmentation challenge for characterizing the semantic and clinical efficacy of automatic segmentation methods relevant for monitoring osteoarthritis progression. MATERIALS AND METHODS: A dataset partition consisting of three-dimensional knee MRI from 88 retrospective patients at two time points (baseline and 1-year follow-up) with ground truth articular (femoral, tibial, and patellar) cartilage and meniscus segmentations was standardized. Challenge submissions and a majority-vote ensemble were evaluated against ground truth segmentations using Dice score, average symmetric surface distance, volumetric overlap error, and coefficient of variation on a holdout test set. Similarities in automated segmentations were measured using pairwise Dice coefficient correlations. Articular cartilage thickness was computed longitudinally and with scans. Correlation between thickness error and segmentation metrics was measured using the Pearson correlation coefficient. Two empirical upper bounds for ensemble performance were computed using combinations of model outputs that consolidated true positives and true negatives. RESULTS: Six teams (T 1-T 6) submitted entries for the challenge. No differences were observed across any segmentation metrics for any tissues (P = .99) among the four top-performing networks (T 2, T 3, T 4, T 6). Dice coefficient correlations between network pairs were high (> 0.85). Per-scan thickness errors were negligible among networks T 1-T 4 (P = .99), and longitudinal changes showed minimal bias (< 0.03 mm). Low correlations (ρ < 0.41) were observed between segmentation metrics and thickness error. The majority-vote ensemble was comparable to top-performing networks (P = .99). Empirical upper-bound performances were similar for both combinations (P = .99). CONCLUSION: Diverse networks learned to segment the knee similarly, where high segmentation accuracy did not correlate with cartilage thickness accuracy and voting ensembles did not exceed individual network performance.See also the commentary by Elhalawani and Mak in this issue.Keywords: Cartilage, Knee, MR-Imaging, Segmentation © RSNA, 2020Supplemental material is available for this article.

Zero echo time pediatric musculoskeletal magnetic resonance imaging: initial experience.

BACKGROUND: Projection radiography (XR) is often supplemented by both CT (to evaluate osseous structures with ionizing radiation) and MRI (for marrow and soft-tissue assessment). Zero echo time (ZTE) MR imaging produces a "CT-like" osseous contrast that might obviate CT. OBJECTIVE: This study investigated our institution's initial experience in implementing an isotropic ZTE MR imaging sequence for pediatric musculoskeletal examinations. MATERIALS AND METHODS: Pediatric patients referred for extremity MRI at 3 tesla (T) underwent ZTE MR imaging to yield images with contrast similar to that of CT. A radiograph-like image was also created with ray-sum image processing. We assessed ZTE-CT/XR anatomical image quality (Sanat) from 0 (nondiagnostic) to 5 (outstanding). Further, we made image comparisons on a 5-point scale (Scomp) (range of -2 = conventional CT/XR greater anatomical delineation to +2 = ZTE-CT/XR greater anatomical delineation; 0=same) for three cohorts: (1) ZTE-XR to conventional radiography, (2) ZTE-CT to conventional CT and (3) pathological lesion assessment on ZTE-XR to conventional radiography. We measured cortical thickness of ZTE-XR and ZTE-CT and compared these with conventional imaging. We calculated confidence interval of proportions, Wilcoxon rank sum test and intraclass correlation coefficients for inter-reader agreement. RESULTS: Cohorts 1, 2 and 3 consisted of 40, 20 and 35 cases, respectively (age range 0.6-23.0 years). ZTE-CT versus CT and ZTE-XR versus radiography of cortical thicknesses were not significantly different (P=0.55 and P=0.31, respectively). Cortical delineation was rated diagnostic or better (score of 3, 4 or 5) in all cases (confidence interval of proportions = 100%) for ZTE-CT/XR. Similarly, intramedullary cavity delineation was rated diagnostic or better in all cases for ZTE-CT, and ZTE-XR was at least diagnostic in 58-63% of cases. For cohort 2, cortex and intramedullary cavity Scomp for ZTE-CT was comparable to those of conventional CT, with confidence interval of proportion (sum of score of -1 to +2) of 93-100% and 95%, respectively. Pathology visualized on ZTE-CT/XR was comparable; Scomp confidence interval of proportions was 95%/97-100%, with improved delineation of non-displaced fractures on ZTE-XR. Readers had moderate to near-perfect intraclass correlation coefficient (range=0.60-0.93). CONCLUSION: Implementation of a diagnostic-quality ZTE MRI sequence in the pediatric population is feasible and can be performed as a complementary pulse sequence to enhance musculoskeletal MRI studies. Compared to conventional CT, ZTE has comparable cortical delineation, intramedullary cavity and pathology visualization. While not intended as a replacement for conventional radiography, ZTE-XR provides similar visualization of pathology.