A simple method for rectified noise floor suppression: Phase-corrected real data reconstruction with application to diffusion-weighted imaging

Improved quantitative parameter estimation for prostate T2 relaxometry using convolutional neural networks

OBJECTIVE

Quantitative parameter mapping conventionally relies on curve fitting techniques to estimate parameters from magnetic resonance image series. This study compares conventional curve fitting techniques to methods using neural networks (NN) for measuring T2 in the prostate.

MATERIALS AND METHODS

Large physics-based synthetic datasets simulating T2 mapping acquisitions were generated for training NNs and for quantitative performance comparisons. Four combinations of different NN architectures and training corpora were implemented and compared with four different curve fitting strategies. All methods were compared quantitatively using synthetic data with known ground truth, and further compared on in vivo test data, with and without noise augmentation, to evaluate feasibility and noise robustness.

RESULTS

In the evaluation on synthetic data, a convolutional neural network (CNN), trained in a supervised fashion using synthetic data generated from naturalistic images, showed the highest overall accuracy and precision amongst the methods. On in vivo data, this best performing method produced low-noise T2 maps and showed the least deterioration with increasing input noise levels.

DISCUSSION

This study showed that a CNN, trained with synthetic data in a supervised manner, may provide superior T2 estimation performance compared to conventional curve fitting, especially in low signal-to-noise regions.

Reduced cerebrospinal fluid motion in patients with Parkinson's disease revealed by magnetic resonance imaging with low b-value diffusion weighted imaging

BACKGROUND

Parkinson's disease is characterized by dopamine-responsive symptoms as well as aggregation of α-synuclein protofibrils. New diagnostic methods assess α-synuclein aggregation characteristics from cerebrospinal fluid (CSF) and recent pathophysiologic mechanisms suggest that CSF circulation disruptions may precipitate α-synuclein retention. Here, diffusion-weighted MRI with low-to-intermediate diffusion-weightings was applied to test the hypothesis that CSF motion is reduced in Parkinson's disease relative to healthy participants.

METHODS

Multi-shell diffusion weighted MRI (spatial resolution = 1.8 × 1.8 × 4.0 mm) with low-to-intermediate diffusion weightings (b-values = 0, 50, 100, 200, 300, 700, and 1000 s/mm2) was applied over the approximate kinetic range of suprasellar cistern fluid motion at 3 Tesla in Parkinson's disease (n = 27; age = 66 ± 6.7 years) and non-Parkinson's control (n = 32; age = 68 ± 8.9 years) participants. Wilcoxon rank-sum tests were applied to test the primary hypothesis that the noise floor-corrected decay rate of CSF signal as a function of b-value, which reflects increasing fluid motion, is reduced within the suprasellar cistern of persons with versus without Parkinson's disease and inversely relates to choroid plexus activity assessed from perfusion-weighted MRI (significance-criteria: p < 0.05).

RESULTS

Consistent with the primary hypothesis, CSF decay rates were higher in healthy (D = 0.00673 ± 0.00213 mm2/s) relative to Parkinson's disease (D = 0.00517 ± 0.00110 mm2/s) participants. This finding was preserved after controlling for age and sex and was observed in the posterior region of the suprasellar cistern (p < 0.001). An inverse correlation between choroid plexus perfusion and decay rate in the voxels within the suprasellar cistern (Spearman's-r=-0.312; p = 0.019) was observed.

CONCLUSIONS

Multi-shell diffusion MRI was applied to identify reduced CSF motion at the level of the suprasellar cistern in adults with versus without Parkinson's disease; the strengths and limitations of this methodology are discussed in the context of the growing literature on CSF flow.

Evaluation of deep learning reconstruction on diffusion-weighted imaging quality and apparent diffusion coefficient using an ice-water phantom

This study assessed the influence of deep learning reconstruction (DLR) on the quality of diffusion-weighted images (DWI) and apparent diffusion coefficient (ADC) using an ice-water phantom. An ice-water phantom with known diffusion properties (true ADC = 1.1 × 10-3 mm2/s at 0 °C) was imaged at various b-values (0, 1000, 2000, and 4000 s/mm2) using a 3 T magnetic resonance imaging scanner with slice thicknesses of 1.5 and 3.0 mm. All DWIs were reconstructed with or without DLR. ADC maps were generated using combinations of b-values 0 and 1000, 0 and 2000, and 0 and 4000 s/mm2. Based on the quantitative imaging biomarker alliance profile, the signal-to-noise ratio (SNRs) in DWIs was calculated, and the accuracy, precision, and within-subject parameter variance (wCV) of the ADCs were evaluated. DLR improved the SNR in DWIs with b-values ranging from 0 to 2000s/mm2; however, its effectiveness was diminished at 4000 s/mm2. There was no noticeable difference in the ADCs of images generated with or without implementing DLR. For a slice thickness of 1.5 mm and combined b-values of 0 and 4000 s/mm2, the ADC values were 0.97 × 10-3and 0.98 × 10-3mm2/s with and without DLR, respectively, both being lower than the true ADC value. Furthermore, DLR enhanced the precision and wCV of the ADC measurements. DLR can enhance the SNR, repeatability, and precision of ADC measurements; however, it does not improve their accuracies.

Recommendations for improved reproducibility of ADC derivation on behalf of the Elekta MRI-linac consortium image analysis working group

BACKGROUND AND PURPOSE

The apparent diffusion coefficient (ADC), a potential imaging biomarker for radiotherapy response, needs to be reproducible before translation into clinical use. The aim of this study was to evaluate the multi-centre delineation- and calculation-related ADC variation and give recommendations to minimize it.

MATERIALS AND METHODS

Nine centres received identical diffusion-weighted and anatomical magnetic resonance images of different cancerous tumours (adrenal gland, pelvic oligo metastasis, pancreas, and prostate). All centres delineated the gross tumour volume (GTV), clinical target volume (CTV), and viable tumour volume (VTV), and calculated ADCs using both their local calculation methods and each of the following calculation conditions: b-values 0-500 vs. 150-500 s/mm2, region-of-interest (ROI)-based vs. voxel-based calculation, and mean vs. median. ADC variation was assessed using the mean coefficient of variation across delineations (CVD) and calculation methods (CVC). Absolute ADC differences between calculation conditions were evaluated using Friedman's test. Recommendations for ADC calculation were formulated based on observations and discussions within the Elekta MRI-linac consortium image analysis working group.

RESULTS

The median (range) CVD and CVC were 0.06 (0.02-0.32) and 0.17 (0.08-0.26), respectively. The ADC estimates differed 18% between b-value sets and 4% between ROI/voxel-based calculation (p-values < 0.01). No significant difference was observed between mean and median (p = 0.64). Aligning calculation conditions between centres reduced CVC to 0.04 (0.01-0.16). CVD was comparable between ROI types.

CONCLUSION

Overall, calculation methods had a larger impact on ADC reproducibility compared to delineation. Based on the results, significant sources of variation were identified, which should be considered when initiating new studies, in particular multi-centre investigations.

Axon fiber orientation as the source of T1 relaxation anisotropy in white matter: A study on corpus callosum in vivo and ex vivo

PURPOSE

Recent studies indicate that T₁ in white matter (WM) is influenced by fiber orientation in B₀ . The purpose of the study was to investigate the interrelationships between axon fiber orientation in corpus callosum (CC) and T₁ relaxation time in humans in vivo as well as in rat brain ex vivo.

METHODS

Volunteers were scanned for relaxometric and diffusion MRI at 3 T and 7 T. Angular T₁ plots from WM were computed using fractional anisotropy and fiber-to-field-angle maps. T₁ and fiber-to-field angle were measured in five sections of CC to estimate the effects of inherently varying fiber orientations on T₁ within the same tracts in vivo. Ex vivo rat-brain preparation encompassing posterior CC was rotated in B₀ and T₁ , and diffusion MRI images acquired at 9.4 T. T₁ angular plots were determined at several rotation angles in B₀ .

RESULTS

Angular T₁ plots from global WM provided reference for estimated fiber orientation-linked T₁ changes within CC. In anterior midbody of CC in vivo, where small axons are dominantly present, a shift in axon orientation is accompanied by a change in T₁ , matching that estimated from WM T₁ data. In CC, where large and giant axons are numerous, the measured T₁ change is about 2-fold greater than the estimated one. Ex vivo rotation of the same midsagittal CC region of interest produced angular T₁ plots at 9.4 T, matching those observed at 7 T in vivo.

CONCLUSION

These data causally link axon fiber orientation in B₀ to the T₁ relaxation anisotropy in WM.

Improved Quantitative Parameter Estimation for Prostate T2 Relaxometry using Convolutional Neural Networks

This work seeks to evaluate multiple methods for quantitative parameter estimation from standard T2 mapping acquisitions in the prostate. The T2 estimation performance of methods based on neural networks (NN) was quantitatively compared to that of conventional curve fitting techniques. Large physics-based synthetic datasets simulating T2 mapping acquisitions were generated for training NNs and for quantitative performance comparisons. Ten combinations of different NN architectures, training strategies, and training corpora were implemented and compared with four different curve fitting strategies. All methods were compared quantitatively using synthetic data with known ground truth, and further compared on in vivo test data, with and without noise augmentation, to evaluate feasibility and noise robustness. In the evaluation on synthetic data, a convolutional neural network (CNN), trained in a supervised fashion using synthetic data generated from naturalistic images, showed the highest overall accuracy and precision amongst all the methods. On in vivo data, this best-performing method produced low-noise T2 maps and showed the least deterioration with increasing input noise levels. This study showed that a CNN, trained with synthetic data in a supervised manner, may provide superior T2 estimation performance compared to conventional curve fitting, especially in low signal-to-noise regions.

White matter microstructure and longitudinal relaxation time anisotropy in human brain at 3 and 7 T

A high degree of structural order by white matter (WM) fibre tracts creates a physicochemical environment where water relaxations are rendered anisotropic. Recently, angularly dependent longitudinal relaxation has been reported in human WM. We have characterised interrelationships between T1 relaxation and diffusion MRI microstructural indices at 3 and 7 T. Eleven volunteers consented to participate in the study. Multishell diffusion MR images were acquired with b-values of 0/1500/3000 and 0/1000/2000 s/mm² at 1.5 and 1.05 mm³ isotropic resolutions at 3 and 7 T, respectively. DTIFIT was used to compute DTI indices; the fibre-to-field angle (θ_FB ) maps were obtained using the principal eigenvector images. The orientations and volume fractions of multiple fibre populations were estimated using BedpostX in FSL, and the orientation dispersion index (ODI) was estimated using the NODDI protocol. MP2RAGE was used to acquire images for T1 maps at 1.0 and 0.9 mm³ isotropic resolutions at 3 and 7 T, respectively. At 3 T, T1 as a function of θ_FB in WM with high fractional anisotropy and one-fibre orientation volume fraction or low ODI shows a broad peak centred at 50^o , but a flat baseline at 0^o and 90^o . The broad peak amounted up to 7% of the mean T1. At 7 T, the broad peak appeared at 40^o and T1 in fibres running parallel to B0 was longer by up to 75 ms (8.3% of the mean T1) than in those perpendicular to the field. The peak at 40^o was approximately 5% of mean T1 (i.e., proportionally smaller than that at 54^o at 3 T). The data demonstrate T1 anisotropy in WM with high microstructural order at both fields. The angular patterns are indicative of the B0-dependency of T1 anisotropy. Thus myelinated WM fibres influence T1 contrast both by acting as a T1 contrast agent and rendering T1 dependent on fibre orientation with B0.

Does perfect filtering really guarantee perfect phase correction for diffusion MRI data?

Owing to its merit of avoiding noise-floor, phase correction is recently used to reconstruct real-valued diffusion MRI data by employing an image filter to estimate the noise-free background phase. However, several studies report an unexpected signal-loss issue for their reconstruction results, with its causing reason still remaining unclear. Although phase correction has achieved promising results in mitigating the signal-loss issue via improving the employed image filter, we have observed counterintuitive results that an advanced filter generates severe artifacts in our previous work. Considering the potential issues with phase correction procedures, in this paper, we argue that even a perfect image filter is insufficient to produce perfect phase correction. To point out the reason why phase correction introduces signal-loss and address this issue, we first propose a complex polar coordinate system (CPCS) to analyze its procedures in detail; second, based on CPCS, we find that phase correction has not sufficiently utilized the background phase, and thus propose a quantitative criterion to fully exploit the background phase; eventually, we propose a phase calibration procedure to remedy current phase correction. Extensive experimental results, including those on synthetic and real diffusion MRI data, demonstrate that our proposed method significantly reduces signal-loss and also eliminates artifacts in FA maps, particularly with improved accuracy on FA.

Diffusion MRI with pulsed and free gradient waveforms: Effects of restricted diffusion and exchange

Monitoring time dependence with diffusion MRI yields observables sensitive to compartment sizes (restricted diffusion) and membrane permeability (water exchange). However, restricted diffusion and exchange have opposite effects on the diffusion-weighted signal, which can lead to errors in parameter estimates. In this work, we propose a signal representation that incorporates the effects of both restricted diffusion and exchange up to second order in b-value and is compatible with gradient waveforms of arbitrary shape. The representation features mappings from a gradient waveform to two scalars that separately control the sensitivity to restriction and exchange. We demonstrate that these scalars span a two-dimensional space that can be used to choose waveforms that selectively probe restricted diffusion or exchange, eliminating the correlation between the two phenomena. We found that waveforms with specific but unconventional shapes provide an advantage over conventional pulsed and oscillating gradient acquisitions. We also show that parametrization of waveforms into a two-dimensional space can be used to understand protocols from other approaches that probe restricted diffusion and exchange. For example, we found that the variation of mixing time in filter-exchange imaging corresponds to variation of our exchange-weighting scalar at a fixed value of the restriction-weighting scalar. The proposed signal representation was evaluated using Monte Carlo simulations in identical parallel cylinders with hexagonal and random packing as well as parallel cylinders with gamma-distributed radii. Results showed that the approach is sensitive to sizes in the interval 4-12 μm and exchange rates in the simulated range of 0 to 20 s - 1 , but also that there is a sensitivity to the extracellular geometry. The presented theory constitutes a simple and intuitive description of how restricted diffusion and exchange influence the signal as well as a guide to protocol design capable of separating the two effects.

Accuracy and precision in super-resolution MRI: Enabling spherical tensor diffusion encoding at ultra-high b-values and high resolution

Diffusion MRI (dMRI) can probe the tissue microstructure but suffers from low signal-to-noise ratio (SNR) whenever high resolution is combined with high diffusion encoding strengths. Low SNR leads to poor precision as well as poor accuracy of the diffusion-weighted signal; the latter is caused by the rectified noise floor and can be observed as a positive bias in magnitude signal. Super-resolution techniques may facilitate a beneficial tradeoff between bias and resolution by allowing acquisition at low spatial resolution and high SNR, whereafter high spatial resolution is recovered by image reconstruction. In this work, we describe a super-resolution reconstruction framework for dMRI and investigate its performance with respect to signal accuracy and precision. Using phantom experiments and numerical simulations, we show that the super-resolution approach improves accuracy by facilitating a more beneficial trade-off between spatial resolution and diffusion encoding strength before the noise floor affects the signal. By contrast, precision is shown to have a less straightforward dependency on acquisition, reconstruction, and intrinsic tissue parameters. Indeed, we find a gain in precision from super-resolution reconstruction is substantial only when some spatial resolution is sacrificed. Finally, we deployed super-resolution reconstruction in a healthy brain for the challenging combination of spherical b-tensor encoding at ultra-high b-values and high spatial resolution—a configuration that produces a unique contrast that emphasizes tissue in which diffusion is restricted in all directions. This demonstration showcased that super-resolution reconstruction enables a vastly superior image contrast compared to conventional imaging, facilitating investigations that would otherwise have prohibitively low SNR, resolution or require non-conventional MRI hardware.

Quantitative diffusion MRI of the abdomen and pelvis

Diffusion MRI has enormous potential and utility in the evaluation of various abdominal and pelvic disease processes including cancer and noncancer imaging of the liver, prostate, and other organs. Quantitative diffusion MRI is based on acquisitions with multiple diffusion encodings followed by quantitative mapping of diffusion parameters that are sensitive to tissue microstructure. Compared to qualitative diffusion-weighted MRI, quantitative diffusion MRI can improve standardization of tissue characterization as needed for disease detection, staging, and treatment monitoring. However, similar to many other quantitative MRI methods, diffusion MRI faces multiple challenges including acquisition artifacts, signal modeling limitations, and biological variability. In abdominal and pelvic diffusion MRI, technical acquisition challenges include physiologic motion (respiratory, peristaltic, and pulsatile), image distortions, and low signal-to-noise ratio. If unaddressed, these challenges lead to poor technical performance (bias and precision) and clinical outcomes of quantitative diffusion MRI. Emerging and novel technical developments seek to address these challenges and may enable reliable quantitative diffusion MRI of the abdomen and pelvis. Through systematic validation in phantoms, volunteers, and patients, including multicenter studies to assess reproducibility, these emerging techniques may finally demonstrate the potential of quantitative diffusion MRI for abdominal and pelvic imaging applications.

Gaussianization of Diffusion MRI Data Using Spatially Adaptive Filtering

Diffusion MRI magnitude data, typically Rician or noncentral χ distributed, is affected by the noise floor, which falsely elevates signal, reduces image contrast, and biases estimation of diffusion parameters. Noise floor can be avoided by extracting real-valued Gaussian-distributed data from complex diffusion-weighted images via phase correction, which is performed by rotating each complex diffusion-weighted image based on its phase so that the actual image content resides in the real part. The imaginary part can then be discarded, leaving only the real part to form a Gaussian-noise image that is not confounded by the noise floor. The effectiveness of phase correction depends on the estimation of the background phase associated with factors such as brain motion, cardiac pulsation, perfusion, and respiration. Most existing smoothing techniques, applied to the real and imaginary images for phase estimation, assume spatially-stationary noise. This assumption does not necessarily hold in real data. In this paper, we introduce an adaptive filtering approach, called multi-kernel filter (MKF), for image smoothing catering to spatially-varying noise. Inspired by the mechanisms of human vision, MKF employs a bilateral filter with spatially-varying kernels. Extensive experiments demonstrate that MKF significantly improves spatial adaptivity and outperforms various state-of-the-art filters in signal Gaussianization.

Probing in vivo cortical myeloarchitecture in humans via line-scan diffusion acquisitions at 7 T with 250-500 micron radial resolution

PURPOSE

The goal of this study was to measure diffusion signals within the cerebral cortex using the line-scan technique to achieve extremely high resolution in the radial direction (ie, perpendicular to the cortical surface) and to demonstrate the utility of these measurements for investigating laminar architecture in the living human brain.

METHODS

Line-scan diffusion data with 250-500 micron radial resolution were acquired at 7 T on 8 healthy volunteers, with each line prescribed perpendicularly to primary somatosensory cortex (S1) and primary motor cortex (M1). Apparent diffusion coefficients, fractional anisotropy values, and radiality indices were measured as a function of cortical depth.

RESULTS

In the deep layers of S1, we found evidence for high anisotropy and predominantly tangential diffusion, with low anisotropy observed in superficial S1. In M1, moderate anisotropy and predominantly radial diffusion was seen at almost all cortical depths. These patterns were consistent across subjects and were conspicuous without averaging data across different locations on the cortical sheet.

CONCLUSION

Our results are in accord with the myeloarchitecture of S1 and M1, known from prior histology studies: in S1, dense bands of tangential myelinated fibers run through the deep layers but not the superficial ones, and in M1, radial myelinated fibers are prominent at most cortical depths. This work therefore provides support for the idea that high-resolution diffusion signals, measured with the line-scan technique and receiving a boost in SNR at 7 T, may serve as a sensitive probe of in vivo laminar architecture.

Removing rician bias in diffusional kurtosis of the prostate using real-data reconstruction

PURPOSE

To compare prostate diffusional kurtosis imaging (DKI) metrics generated using phase-corrected real data with those generated using magnitude data with and without noise compensation (NC).

METHODS

Diffusion-weighted images were acquired at 3T in 16 prostate cancer patients, measuring 6 b-values (0-1500 s/mm2 ), each acquired with 6 signal averages along 3 diffusion directions, with noise-only images acquired to allow NC. In addition to conventional magnitude averaging, phase-corrected real data were averaged in an attempt to reduce rician noise-bias, with a range of phase-correction low-pass filter (LPF) sizes (8-128 pixels) tested. Each method was also tested using simulations. Pixelwise maps of apparent diffusion (D) and apparent kurtosis (K) were calculated for magnitude data with and without NC and phase-corrected real data. Average values were compared in tumor, normal transition zone (NTZ), and normal peripheral zone (NPZ).

RESULTS

Simulations indicated LPF size can strongly affect K metrics, where 64-pixel LPFs produced accurate metrics. Relative to metrics estimated from magnitude data without NC, median NC K were lower (P < 0.0001) by 6/11/8% in tumor/NPZ/NTZ, 64-LPF real-data K were lower (P < 0.0001) by 4/10/7%, respectively.

CONCLUSION

Compared with magnitude data with NC, phase-corrected real data can produce similar K, although the choice of phase-correction LPF should be chosen carefully.

Adaptive phase correction of diffusion-weighted images

Phase correction (PC) is a preprocessing technique that exploits the phase of images acquired in Magnetic Resonance Imaging (MRI) to obtain real-valued images containing tissue contrast with additive Gaussian noise, as opposed to magnitude images which follow a non-Gaussian distribution, e.g. Rician. PC finds its natural application to diffusion-weighted images (DWIs) due to their inherent low signal-to-noise ratio and consequent non-Gaussianity that induces a signal overestimation bias that propagates to the calculated diffusion indices. PC effectiveness depends upon the quality of the phase estimation, which is often performed via a regularization procedure. We show that a suboptimal regularization can produce alterations of the true image contrast in the real-valued phase-corrected images. We propose adaptive phase correction (APC), a method where the phase is estimated by using MRI noise information to perform a complex-valued image regularization that accounts for the local variance of the noise. We show, on synthetic and acquired data, that APC leads to phase-corrected real-valued DWIs that present a reduced number of alterations and a reduced bias. The substantial absence of parameters for which human input is required favors a straightforward integration of APC in MRI processing pipelines.

Accelerated Segmented Diffusion-Weighted Prostate Imaging for Higher Resolution, Higher Geometric Fidelity, and Multi-b Perfusion Estimation

PURPOSE

The aim of this study was to improve the geometric fidelity and spatial resolution of multi-b diffusion-weighted magnetic resonance imaging of the prostate.

MATERIALS AND METHODS

An accelerated segmented diffusion imaging sequence was developed and evaluated in 25 patients undergoing multiparametric magnetic resonance imaging examinations of the prostate. A reduced field of view was acquired using an endorectal coil. The number of sampled diffusion weightings, or b-factors, was increased to allow estimation of tissue perfusion based on the intravoxel incoherent motion (IVIM) model. Apparent diffusion coefficients measured with the proposed segmented method were compared with those obtained with conventional single-shot echo-planar imaging (EPI).

RESULTS

Compared with single-shot EPI, the segmented method resulted in faster acquisition with 2-fold improvement in spatial resolution and a greater than 3-fold improvement in geometric fidelity. Apparent diffusion coefficient values measured with the novel sequence demonstrated excellent agreement with those obtained from the conventional scan (R = 0.91 for bmax = 500 s/mm and R = 0.89 for bmax = 1400 s/mm). The IVIM perfusion fraction was 4.0% ± 2.7% for normal peripheral zone, 6.6% ± 3.6% for normal transition zone, and 4.4% ± 2.9% for suspected tumor lesions.

CONCLUSIONS

The proposed accelerated segmented prostate diffusion imaging sequence achieved improvements in both spatial resolution and geometric fidelity, along with concurrent quantification of IVIM perfusion.

Complex diffusion-weighted image estimation via matrix recovery under general noise models

We propose a patch-based singular value shrinkage method for diffusion magnetic resonance image estimation targeted at low signal to noise ratio and accelerated acquisitions. It operates on the complex data resulting from a sensitivity encoding reconstruction, where asymptotically optimal signal recovery guarantees can be attained by modeling the noise propagation in the reconstruction and subsequently simulating or calculating the limit singular value spectrum. Simple strategies are presented to deal with phase inconsistencies and optimize patch construction. The pertinence of our contributions is quantitatively validated on synthetic data, an in vivo adult example, and challenging neonatal and fetal cohorts. Our methodology is compared with related approaches, which generally operate on magnitude-only data and use data-based noise level estimation and singular value truncation. Visual examples are provided to illustrate effectiveness in generating denoised and debiased diffusion estimates with well preserved spatial and diffusion detail.

Some new angles on the magic angle: what MSK radiologists know and don't know about this phenomenon

PURPOSE

Magic angle effects (MAE) are well-recognized in musculoskeletal (MSK) MRI. With short TE acquisitions, the signal intensity of tendons, ligaments, and menisci depend on their orientation relative to the main magnetic field (B₀). An interactive resident physics teaching module simulating MR imaging of a tendon forced us to identify and correct several misconceptions we had about MAE. We suspected these misconceptions were shared by other MSK radiologists.

MATERIALS AND METHODS

We surveyed members of the Society of Academic Bone Radiologists (SABR) regarding which pulse sequences, acquisition parameters, tissues and angles relative to B₀ were most likely to produce MAE.

RESULTS

Survey respondents knew that MAE strongly depend on TE and commonly appear on T1W, FSE and PD sequences, but were less aware that MAE may also appear on T2W, STIR and DWI sequences. They knew of MAE effects in tendons, ligaments and cartilage, but were less aware of those in entheses, peripheral nerves and intervertebral discs. Respondents underestimated the wide angular range (full-width at half-maximum ≈ 40^∘) over which significant MAE can be seen with short TE.

CONCLUSIONS

Collagen-containing tissues with parallel molecular alignment exhibit increased signal intensity when oriented at 55^∘ relative to B₀. Experienced MSK radiologists were found to underestimate the combinations of image parameters, pulse sequences, tissues and collagen orientations in which significant MAE may be seen. Our survey results highlight the need for ongoing MR physics education for practicing radiologists.

Computation of exact g-factor maps in 3D GRAPPA reconstructions

PURPOSE

To characterize the noise distributions in 3D-MRI accelerated acquisitions reconstructed with GRAPPA using an exact noise propagation analysis that operates directly in k-space.

THEORY AND METHODS

We exploit the extensive symmetries and separability in the reconstruction steps to account for the correlation between all the acquired k-space samples. Monte Carlo simulations and multi-repetition phantom experiments were conducted to test both the accuracy and feasibility of the proposed method; a high-resolution in-vivo experiment was performed to assess the applicability of our method to clinical scenarios.

RESULTS

Our theoretical derivation shows that the direct k-space analysis renders an exact noise characterization under the assumptions of stationarity and uncorrelation in the original k-space. Simulations and phantom experiments provide empirical support to the theoretical proof. Finally, the high-resolution in-vivo experiment demonstrates the ability of the proposed method to assess the impact of the sub-sampling pattern on the overall noise behavior.

CONCLUSIONS

By operating directly in the k-space, the proposed method is able to provide an exact characterization of noise for any Cartesian pattern sub-sampled along the two phase-encoding directions. Exploitation of the symmetries and separability into independent blocks through the image reconstruction procedure allows us to overcome the computational challenges related to the very large size of the covariance matrices involved.

Evaluation of Head and Neck Tumors with Functional MR Imaging

Head and neck cancer is one of the most common cancers worldwide. MR imaging-based diffusion and perfusion techniques enable the noninvasive assessment of tumor biology and physiology, which supplement information obtained from standard structural scans. Diffusion and perfusion MR imaging techniques provide novel biomarkers that can aid monitoring in pretreatment, during treatment, and posttreatment stages to improve patient selection for therapeutic strategies; provide evidence for change of therapy regime; and evaluate treatment response. This review discusses pertinent aspects of the role of diffusion and perfusion MR imaging and computational analysis methods in studying head and neck cancer.

Multi-institutional validation of a novel textural analysis tool for preoperative stratification of suspected thyroid tumors on diffusion-weighted MRI

PURPOSE

Ultrasound-guided fine needle aspirate cytology fails to diagnose many malignant thyroid nodules; consequently, patients may undergo diagnostic lobectomy. This study assessed whether textural analysis (TA) could noninvasively stratify thyroid nodules accurately using diffusion-weighted MRI (DW-MRI).

METHODS

This multi-institutional study examined 3T DW-MRI images obtained with spin echo echo planar imaging sequences. The training data set included 26 patients from Cambridge, United Kingdom, and the test data set included 18 thyroid cancer patients from Memorial Sloan Kettering Cancer Center (New York, New York, USA). Apparent diffusion coefficients (ADCs) were compared over regions of interest (ROIs) defined on thyroid nodules. TA, linear discriminant analysis (LDA), and feature reduction were performed using the 21 MaZda-generated texture parameters that best distinguished benign and malignant ROIs.

RESULTS

Training data set mean ADC values were significantly different for benign and malignant nodules (P = 0.02) with a sensitivity and specificity of 70% and 63%, respectively, and a receiver operator characteristic (ROC) area under the curve (AUC) of 0.73. The LDA model of the top 21 textural features correctly classified 89/94 DW-MRI ROIs with 92% sensitivity, 96% specificity, and an AUC of 0.97. This algorithm correctly classified 16/18 (89%) patients in the independently obtained test set of thyroid DW-MRI scans.

CONCLUSION

TA classifies thyroid nodules with high sensitivity and specificity on multi-institutional DW-MRI data sets. This method requires further validation in a larger prospective study. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance.

Using diffusion-weighted MRI to predict aggressive histological features in papillary thyroid carcinoma: a novel tool for pre-operative risk stratification in thyroid cancer

BACKGROUND

Initial management recommendations of papillary thyroid carcinoma (PTC) are very dependent on preoperative studies designed to evaluate the presence of PTC with aggressive features. The purpose of this study was to evaluate whether diffusion-weighted magnetic resonance imaging (DW-MRI) before surgery can be used as a tool to stratify tumor aggressiveness in patients with PTC.

METHODS

In this prospective study, 28 patients with PTC underwent DW-MRI studies on a three Tesla MR scanner prior to thyroidectomy. Due to image quality, 21 patients were finally suitable for further analysis. Apparent diffusion coefficients (ADCs) of normal thyroid tissues and PTCs for 21 patients were calculated. Tumor aggressiveness was defined by surgical histopathology. The Mann-Whitney U test was used to compare the difference in ADCs among groups of normal thyroid tissues and PTCs with and without features of tumor aggressiveness. Receiver operating characteristic (ROC) analysis was performed to assess the discriminative specificity, sensitivity, and accuracy of and determine the cutoff value for the ADC in stratifying PTCs with tumor aggressiveness.

RESULTS

There was no significant difference in ADC values between normal thyroid tissues and PTCs. However, ADC values of PTCs with extrathyroidal extension (ETE; 1.53±0.25×10(-3) mm2/s) were significantly lower than corresponding values from PTCs without ETE (2.37±0.67×10(-3) mm2/s; p<0.005). ADC values identified 3 papillary carcinoma patients with extrathyroidal extension that would have otherwise been candidates for observation based on ultrasound evaluations. The cutoff value of ADC to discriminate PTCs with and without ETE was determined at 1.85×10(-3) mm2/s with a sensitivity of 85%, specificity of 85%, and ROC curve area of 0.85.

CONCLUSION

ADC value derived from DW-MRI before surgery has the potential to stratify ETE in patients with PTCs.

Influence of noise correction on intra- and inter-subject variability of quantitative metrics in diffusion kurtosis imaging

Diffusion kurtosis imaging (DKI) is a promising extension of diffusion tensor imaging, giving new insights into the white matter microstructure and providing new biomarkers. Given the rapidly increasing number of studies, DKI has a potential to establish itself as a valuable tool in brain diagnostics. However, to become a routine procedure, DKI still needs to be improved in terms of robustness, reliability, and reproducibility. As it requires acquisitions at higher diffusion weightings, results are more affected by noise than in diffusion tensor imaging. The lack of standard procedures for post-processing, especially for noise correction, might become a significant obstacle for the use of DKI in clinical routine limiting its application. We considered two noise correction schemes accounting for the noise properties of multichannel phased-array coils, in order to improve the data quality at signal-to-noise ratio (SNR) typical for DKI. The SNR dependence of estimated DKI metrics such as mean kurtosis (MK), mean diffusivity (MD) and fractional anisotropy (FA) is investigated for these noise correction approaches in Monte Carlo simulations and in in vivo human studies. The intra-subject reproducibility is investigated in a single subject study by varying the SNR level and SNR spatial distribution. Then the impact of the noise correction on inter-subject variability is evaluated in a homogeneous sample of 25 healthy volunteers. Results show a strong impact of noise correction on the MK estimate, while the estimation of FA and MD was affected to a lesser extent. Both intra- and inter-subject SNR-related variability of the MK estimate is considerably reduced after correction for the noise bias, providing more accurate and reproducible measures. In this work, we have proposed a straightforward method that improves accuracy of DKI metrics. This should contribute to standardization of DKI applications in clinical studies making valuable inferences in group analysis and longitudinal studies.

Spatially variable Rician noise in magnetic resonance imaging

Magnetic resonance images tend to be influenced by various random factors usually referred to as "noise". The principal sources of noise and related artefacts can be divided into two types: arising from hardware (acquisition coil arrays, gradient coils, field inhomogeneity); and arising from the subject (physiological noise including body motion, cardiac pulsation or respiratory motion). These factors negatively affect the resolution and reproducibility of the images. Therefore, a proper noise treatment is important for improving the performance of clinical and research investigations. Noise reduction becomes especially critical for the images with a low signal-to-noise ratio, such as those typically acquired in diffusion tensor imaging at high diffusion weightings. The standard methods of signal correction usually assume a uniform distribution of the standard deviation of the noise across the image and evaluate a single correction parameter for the whole image. We pursue a more advanced approach based on the assumption of an inhomogeneous distribution of noise in space and evaluate correction factors for each voxel individually. The Rician nature of the underlying noise is considered for low and high signal-to-noise ratios. The approach developed here has been examined using numerical simulations and in vivo brain diffusion tensor imaging experiments. The efficacy and usefulness of this approach is demonstrated here and the resultant effective tool is described.