Analytically exact correction scheme for signal extraction from noisy magnitude MR signals

Human Olivocochlear Effects: A Statistical Detection Approach Applied to the Cochlear Microphonic Evoked by Swept Tones

The human medial olivocochlear (MOC) reflex was assessed by observing the effects of contralateral acoustic stimulation (CAS) on the cochlear microphonic (CM) across a range of probe frequencies. A frequency-swept probe tone (125-4757 Hz, 90 dB SPL) was presented in two directions (up sweep and down sweep) to normal-hearing young adults. This study assessed MOC effects on the CM in individual participants using a statistical approach that calculated minimum detectable changes in magnitude and phase based on CM signal-to-noise ratio (SNR). Significant increases in CM magnitude, typically 1-2 dB in size, were observed for most participants from 354 to 1414 Hz, where the size and consistency of these effects depended on participant, probe frequency, sweep direction, and SNR. CAS-related phase lags were also observed, consistent with CM-based MOC studies in laboratory animals. Observed effects on CM magnitude and phase were in the opposite directions of reported effects on otoacoustic emissions (OAEs). OAEs are sensitive to changes in the motility of outer hair cells located near the peak region of the traveling wave, while the effects of CAS on the CM likely originate from MOC-related changes in the conductance of outer hair cells located in the basal tail of the traveling wave. Thus, MOC effects on the CM are complementary to those observed for OAEs.

New non-local mean methods for MRI denoising based on global self-similarity between values

Magnetic resonance imaging (MRI) is a non-invasive medical imaging technique that provides high-resolution 3D images and valuable insights into human tissue conditions. Even at present, the refinement of denoising methods for MRI remains a crucial concern for improving the quality of the images. This study aims to improve the prefiltered rotationally invariant non-local principal component analysis (PRI-NL-PCA) algorithm. We relaxed the original restrictions using particle swarm optimization to determine optimal parameters for the PCA part of the original algorithm. In addition, we adjusted the prefiltered rotationally invariant non-local mean (PRI-NLM) part by traversing the signal intensities of voxels instead of their spatial positions to reduce duplicate calculations and expand the search volume to the whole image when estimating voxels' signal intensities. The new method demonstrated superior denoising performance compared to the original approach. Moreover, in most cases, the new algorithm ran faster. Furthermore, our proposed method can also be applied to process Gaussian noise in natural images and has the potential to enhance other NLM-based denoising algorithms.

Assessment of Mechanically Induced Changes in Helical Fiber Microstructure Using Diffusion Tensor Imaging

Noninvasive methods to detect microstructural changes in collagen-based fibrous tissues are necessary to differentiate healthy from damaged tissues in vivo but are sparse. Diffusion Tensor Imaging (DTI) is a noninvasive imaging technique used to quantitatively infer tissue microstructure with previous work primarily focused in neuroimaging applications. Yet, it is still unclear how DTI metrics relate to fiber microstructure and function in musculoskeletal tissues such as ligament and tendon, in part because of the high heterogeneity inherent to such tissues. To address this limitation, we assessed the ability of DTI to detect microstructural changes caused by mechanical loading in tissue-mimicking helical fiber constructs of known structure. Using high-resolution optical and micro-computed tomography imaging, we found that static and fatigue loading resulted in decreased sample diameter and a re-alignment of the macro-scale fiber twist angle similar with the direction of loading. However, DTI and micro-computed tomography measurements suggest microstructural differences in the effect of static versus fatigue loading that were not apparent at the bulk level. Specifically, static load resulted in an increase in diffusion anisotropy and a decrease in radial diffusivity suggesting radially uniform fiber compaction. In contrast, fatigue loads resulted in increased diffusivity in all directions and a change in the alignment of the principal diffusion direction away from the constructs' main axis suggesting fiber compaction and microstructural disruptions in fiber architecture. These results provide quantitative evidence of the ability of DTI to detect mechanically induced changes in tissue microstructure that are not apparent at the bulk level, thus confirming its potential as a noninvasive measure of microstructure in helically architected collagen-based tissues, such as ligaments and tendons.

The effects of axonal beading and undulation on axonal diameter estimation from diffusion MRI: Insights from simulations in human axons segmented from three-dimensional electron microscopy

The increasing availability of high-performance gradient systems in human MRI scanners has generated great interest in diffusion microstructural imaging applications such as axonal diameter mapping. Practically, sensitivity to axon diameter in diffusion MRI is attained at strong diffusion weightings b , where the deviation from the expected 1 / b scaling in white matter yields a finite transverse diffusivity, which is then translated into an axon diameter estimate. While axons are usually modeled as perfectly straight, impermeable cylinders, local variations in diameter (caliber variation or beading) and direction (undulation) are known to influence axonal diameter estimates and have been observed in microscopy data of human axons. In this study, we performed Monte Carlo simulations of diffusion in axons reconstructed from three-dimensional electron microscopy of a human temporal lobe specimen using simulated sequence parameters matched to the maximal gradient strength of the next-generation Connectome 2.0 human MRI scanner ( ≲ 500 mT/m). We show that axon diameter estimation is accurate for nonbeaded, nonundulating fibers; however, in fibers with caliber variations and undulations, the axon diameter is heavily underestimated due to caliber variations, and this effect overshadows the known overestimation of the axon diameter due to undulations. This unexpected underestimation may originate from variations in the coarse-grained axial diffusivity due to caliber variations. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.

Impact of truncating diffusion MRI scans on diffusional kurtosis imaging

OBJECTIVE

Diffusional kurtosis imaging (DKI) extends diffusion tensor imaging (DTI), characterizing non-Gaussian diffusion effects but requires longer acquisition times. To ensure the robustness of DKI parameters, data acquisition ordering should be optimized allowing for scan interruptions or shortening. Three methodologies were used to examine how reduced diffusion MRI scans impact DKI histogram-metrics: 1) the electrostatic repulsion model (OptEEM); 2) spherical codes (OptSC); 3) random (RandomTRUNC).

MATERIALS AND METHODS

Pre-acquired diffusion multi-shell data from 14 female healthy volunteers (29±5 years) were used to generate reordered data. For each strategy, subsets containing different amounts of the full dataset were generated. The subsampling effects were assessed on histogram-based DKI metrics from tract-based spatial statistics (TBSS) skeletonized maps. To evaluate each subsampling method on simulated data at different SNRs and the influence of subsampling on in vivo data, we used a 3-way and 2-way repeated measures ANOVA, respectively.

RESULTS

Simulations showed that subsampling had different effects depending on DKI parameter, with fractional anisotropy the most stable (up to 5% error) and radial kurtosis the least stable (up to 26% error). RandomTRUNC performed the worst while the others showed comparable results. Furthermore, the impact of subsampling varied across distinct histogram characteristics, the peak value the least affected (OptEEM: up to 5% error; OptSC: up to 7% error) and peak height (OptEEM: up to 8% error; OptSC: up to 11% error) the most affected.

CONCLUSION

The impact of truncation depends on specific histogram-based DKI metrics. The use of a strategy for optimizing the acquisition order is advisable to improve DKI robustness to exam interruptions.

A translational MRI approach to validate acute axonal damage detection as an early event in multiple sclerosis

Axonal degeneration is a central pathological feature of multiple sclerosis and is closely associated with irreversible clinical disability. Current noninvasive methods to detect axonal damage in vivo are limited in their specificity and clinical applicability, and by the lack of proper validation. We aimed to validate an MRI framework based on multicompartment modeling of the diffusion signal (AxCaliber) in rats in the presence of axonal pathology, achieved through injection of a neurotoxin damaging the neuronal terminal of axons. We then applied the same MRI protocol to map axonal integrity in the brain of multiple sclerosis relapsing-remitting patients and age-matched healthy controls. AxCaliber is sensitive to acute axonal damage in rats, as demonstrated by a significant increase in the mean axonal caliber along the targeted tract, which correlated with neurofilament staining. Electron microscopy confirmed that increased mean axonal diameter is associated with acute axonal pathology. In humans with multiple sclerosis, we uncovered a diffuse increase in mean axonal caliber in most areas of the normal-appearing white matter, preferentially affecting patients with short disease duration. Our results demonstrate that MRI-based axonal diameter mapping is a sensitive and specific imaging biomarker that links noninvasive imaging contrasts with the underlying biological substrate, uncovering generalized axonal damage in multiple sclerosis as an early event.

DENOISING DIFFUSION MRI: CONSIDERATIONS AND IMPLICATIONS FOR ANALYSIS

Development of diffusion MRI (dMRI) denoising approaches has experienced considerable growth over the last years. As noise can inherently reduce accuracy and precision in measurements, its effects have been well characterised both in terms of uncertainty increase in dMRI-derived features and in terms of biases caused by the noise floor, the smallest measurable signal given the noise level. However, gaps in our knowledge still exist in objectively characterising dMRI denoising approaches in terms of both of these effects and assessing their efficacy. In this work, we reconsider what a denoising method should and should not do and we accordingly define criteria to characterise the performance. We propose a comprehensive set of evaluations, including i) benefits in improving signal quality and reducing noise variance, ii) gains in reducing biases and the noise floor and improving, iii) preservation of spatial resolution, iv) agreement of denoised data against a gold standard, v) gains in downstream parameter estimation (precision and accuracy), vi) efficacy in enabling noise-prone applications, such as ultra-high-resolution imaging. We further provide newly acquired complex datasets (magnitude and phase) with multiple repeats that sample different SNR regimes to highlight performance differences under different scenarios. Without loss of generality, we subsequently apply a number of exemplar patch-based denoising algorithms to these datasets, including Non-Local Means, Marchenko-Pastur PCA (MPPCA) in the magnitude and complex domain and NORDIC, and compare them with respect to the above criteria and against a gold standard complex average of multiple repeats. We demonstrate that all tested denoising approaches reduce noise-related variance, but not always biases from the elevated noise floor. They all induce a spatial resolution penalty, but its extent can vary depending on the method and the implementation. Some denoising approaches agree with the gold standard more than others and we demonstrate challenges in even defining such a standard. Overall, we show that dMRI denoising performed in the complex domain is advantageous to magnitude domain denoising with respect to all the above criteria.

Frequency-dependent diffusion kurtosis imaging in the human brain using an oscillating gradient spin echo sequence and a high-performance head-only gradient

Measuring the time/frequency dependence of diffusion MRI is a promising approach to distinguish between the effects of different tissue microenvironments, such as membrane restriction, tissue heterogeneity, and compartmental water exchange. In this study, we measure the frequency dependence of diffusivity (D) and kurtosis (K) with oscillating gradient diffusion encoding waveforms and a diffusion kurtosis imaging (DKI) model in human brains using a high-performance, head-only MAGNUS gradient system, with a combination of b-values, oscillating frequencies (f), and echo time that has not been achieved in human studies before. Frequency dependence of diffusivity and kurtosis are observed in both global and local white matter (WM) and gray matter (GM) regions and characterized with a power-law model ∼Λ*fθ. The frequency dependences of diffusivity and kurtosis (including changes between fmin and fmax, Λ, and θ) vary over different WM and GM regions, indicating potential microstructural differences between regions. A trend of decreasing kurtosis over frequency in the short-time limit is successfully captured for in vivo human brains. The effects of gradient nonlinearity (GNL) on frequency-dependent diffusivity and kurtosis measurements are investigated and corrected. Our results show that the GNL has prominent scaling effects on the measured diffusivity values (3.5∼5.5% difference in the global WM and 6∼8% difference in the global cortex) and subsequently affects the corresponding power-law parameters (Λ, θ) while having a marginal influence on the measured kurtosis values (<0.05% difference) and power-law parameters (Λ, θ). This study expands previous OGSE studies and further demonstrates the translatability of frequency-dependent diffusivity and kurtosis measurements to human brains, which may provide new opportunities to probe human brain microstructure in health and disease.

Mapping tissue microstructure across the human brain on a clinical scanner with soma and neurite density image metrics

Soma and neurite density image (SANDI) is an advanced diffusion magnetic resonance imaging biophysical signal model devised to probe in vivo microstructural information in the gray matter (GM). This model requires acquisitions that include b values that are at least six times higher than those used in clinical practice. Such high b values are required to disentangle the signal contribution of water diffusing in soma from that diffusing in neurites and extracellular space, while keeping the diffusion time as short as possible to minimize potential bias due to water exchange. These requirements have limited the use of SANDI only to preclinical or cutting-edge human scanners. Here, we investigate the potential impact of neglecting water exchange in the SANDI model and present a 10-min acquisition protocol that enables to characterize both GM and white matter (WM) on 3 T scanners. We implemented analytical simulations to (i) evaluate the stability of the fitting of SANDI parameters when diminishing the number of shells; (ii) estimate the bias due to potential exchange between neurites and extracellular space in such reduced acquisition scheme, comparing it with the bias due to experimental noise. Then, we demonstrated the feasibility and assessed the repeatability and reproducibility of our approach by computing microstructural metrics of SANDI with AMICO toolbox and other state-of-the-art models on five healthy subjects. Finally, we applied our protocol to five multiple sclerosis patients. Results suggest that SANDI is a practical method to characterize WM and GM tissues in vivo on performant clinical scanners.

The influence of axonal beading and undulation on axonal diameter mapping

We consider the effect of non-cylindrical axonal shape on axonal diameter mapping with diffusion MRI. Practical sensitivity to axon diameter is attained at strong diffusion weightings b , where the deviation from the 1 / b scaling yields the finite transverse diffusivity, which is then translated into axon diameter. While axons are usually modeled as perfectly straight, impermeable cylinders, the local variations in diameter (caliber variation or beading) and direction (undulation) have been observed in microscopy data of human axons. Here we quantify the influence of cellular-level features such as caliber variation and undulation on axon diameter estimation. For that, we simulate the diffusion MRI signal in realistic axons segmented from 3-dimensional electron microscopy of a human brain sample. We then create artificial fibers with the same features and tune the amplitude of their caliber variations and undulations. Numerical simulations of diffusion in fibers with such tunable features show that caliber variations and undulations result in under- and over-estimation of axon diameters, correspondingly; this bias can be as large as 100%. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.

Synthetic-to-real domain adaptation with deep learning for fitting the intravoxel incoherent motion model of diffusion-weighted imaging

BACKGROUND

Intravoxel incoherent motion (IVIM) is a type of diffusion-weighted imaging (DWI), and IVIM model parameters (water molecule diffusion rate D_t , pseudo-diffusion coefficient D_p , and tissue perfusion fraction F_p ) have been widely used in the diagnosis and characterization of malignant lesions.

PURPOSE

This study proposes a deep-learning model with synthetic-to-real domain adaptation to fit the IVIM model parameters of DWI.

METHODS

Ninety-eight consecutive patients diagnosed with hepatocellular carcinoma between January 2017 and September 2020 were included in the study, and routine IVIM-DWI serial examinations were performed using a 3.0 T magnetic resonance imaging system in preoperative MR imaging. The proposed method is mainly composed of two modules: a convolutional neural network-based IVIM model fitting network to map b-value images to the IVIM parameter maps and a domain discriminator to improve the accuracy of the IVIM parameter maps in the real data. The proposed method was compared with previously reported fitting methods, including the nonlinear least squares (NLSs), IVIM-NET_optim , and self-supervised U-network methods. The IVIM parameter-fitting performance was assessed by measuring the DWI reconstruction performance and testing the robustness of each method against noise using noise-corrupted data.

RESULTS

The DWI reconstruction performance demonstrates that the proposed method has better reconstruction accuracy for DWI with a low signal-to-noise ratio, which implies that the proposed method improves the fitting accuracy of the IVIM parameters. Noise-corrupt experiments show that the proposed method is more robust against noise-corrupted signals. With the proposed method, no outliers were found in D_t , and outliers were reduced for F_p in the abnormal regions (proposed method: 1.85%; NLS: 5.90%; IVIM-NET_optim : 6.61%; and self-U-net: 25.36%). Moreover, experiments show that the proposed method has a more stable parameter estimation performance than the existing methods in the absence of real data.

CONCLUSIONS

IVIM parameters can be estimated using a synthetic-to-real domain-adaptation framework with deep learning, and the proposed method outperforms previously reported methods.

Graph theoretical approach to brain remodeling in multiple sclerosis

Multiple sclerosis (MS) is a neuroinflammatory disorder damaging structural connectivity. Natural remodeling processes of the nervous system can, to some extent, restore the damage caused. However, there is a lack of biomarkers to evaluate remodeling in MS. Our objective is to evaluate graph theory metrics (especially modularity) as a biomarker of remodeling and cognition in MS. We recruited 60 relapsing-remitting MS and 26 healthy controls. Structural and diffusion MRI, plus cognitive and disability evaluations, were done. We calculated modularity and global efficiency from the tractography-derived connectivity matrices. Association of graph metrics with T2 lesion load, cognition, and disability was evaluated using general linear models adjusting for age, gender, and disease duration wherever applicable. We showed that MS subjects had higher modularity and lower global efficiency compared with controls. In the MS group, modularity was inversely associated with cognitive performance but positively associated with T2 lesion load. Our results indicate that modularity increase is due to the disruption of intermodular connections in MS because of the lesions, with no improvement or preserving of cognitive functions.

Estimating axial diffusivity in the NODDI model

To estimate microstructure-related parameters from diffusion MRI data, biophysical models make strong, simplifying assumptions about the underlying tissue. The extent to which many of these assumptions are valid remains an open research question. This study was inspired by the disparity between the estimated intra-axonal axial diffusivity from literature and that typically assumed by the Neurite Orientation Dispersion and Density Imaging (NODDI) model (d_∥=1.7μm²/ms). We first demonstrate how changing the assumed axial diffusivity results in considerably different NODDI parameter estimates. Second, we illustrate the ability to estimate axial diffusivity as a free parameter of the model using high b-value data and an adapted NODDI framework. Using both simulated and in vivo data we investigate the impact of fitting to either real-valued or magnitude data, with Gaussian and Rician noise characteristics respectively, and what happens if we get the noise assumptions wrong in this high b-value and thus low SNR regime. Our results from real-valued human data estimate intra-axonal axial diffusivities of ∼2-2.5μm²/ms, in line with current literature. Crucially, our results demonstrate the importance of accounting for both a rectified noise floor and/or a signal offset to avoid biased parameter estimates when dealing with low SNR data.

An unsupervised convolutional neural network method for estimation of intravoxel incoherent motion parameters

Objective. Intravoxel incoherent motion (IVIM) imaging obtained by fitting a biexponential model to multiple b-value diffusion-weighted magnetic resonance imaging (DW-MRI) has been shown to be a promising tool for different clinical applications. Recently, several deep neural network (DNN) methods were proposed to generate IVIM imaging. Approach. In this study, we proposed an unsupervised convolutional neural network (CNN) method for estimation of IVIM parameters. We used both simulated and real abdominal DW-MRI data to evaluate the performance of the proposed CNN-based method, and compared the results with those obtained from a non-linear least-squares fit (TRR, trust-region reflective algorithm) and a feed-forward backward-propagation DNN-based method. Main results. The simulation results showed that both the DNN- and CNN-based methods had lower coefficients of variation than the TRR method, but the CNN-based method provided more accurate parameter estimates. The results obtained from real DW-MRI data showed that the TRR method produced many biased IVIM parameter estimates that hit the upper and lower parameter bounds. In contrast, both the DNN- and CNN-based methods yielded less biased IVIM parameter estimates. Overall, the perfusion fraction and diffusion coefficient obtained from the DNN- and CNN-based methods were close to literature values. However, compared with the CNN-based method, both the TRR and DNN-based methods tended to yield increased pseudodiffusion coefficients (55%–180%). Significance. Our preliminary results suggest that it is feasible to estimate IVIM parameters using CNN.

SENSE EPI reconstruction with 2D phase error correction and channel-wise noise removal

PURPOSE

To develop a robust reconstruction pipeline for EPI data that enables 2D Nyquist phase error correction using sensitivity encoding without incurring major noise artifacts in low SNR data.

METHODS

SENSE with 2D phase error correction (PEC-SENSE) was combined with channel-wise noise removal using Marcenko-Pastur principal component analysis (MPPCA) to simultaneously eliminate Nyquist ghost artifacts in EPI data and mitigate the noise amplification associated with phase correction using parallel imaging. The proposed pipeline (coined SPECTRE) was validated in phantom DW-EPI data using the accuracy and precision of diffusion metrics; ground truth values were obtained from data acquired with a spin echo readout. Results from the SPECTRE pipeline were compared against PEC-SENSE reconstructions with three alternate denoising strategies: (i) no denoising; (ii) denoising of magnitude data after image formation; (iii) denoising of complex data after image formation. SPECTRE was then tested using highb $$ b $$ -value (i.e., low SNR) diffusion data (up tob = 3000 $$ b=3000 $$  s/mm2 $$ {}^2 $$ ) in four healthy subjects.

RESULTS

Noise amplification associated with phase error correction incurred a 23% bias in phantom mean diffusivity (MD) measurements. Phantom MD estimates using the SPECTRE pipeline were within 8% of the ground truth value. In healthy volunteers, the SPECTRE pipeline visibly corrected Nyquist ghost artifacts and reduced associated noise amplification in highb $$ b $$ -value data.

CONCLUSION

The proposed reconstruction pipeline is effective in correcting low SNR data, and improves the accuracy and precision of derived diffusion metrics.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

Mapping the Human Connectome using Diffusion MRI at 300 mT/m Gradient Strength: Methodological Advances and Scientific Impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in Continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength dMRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for dMRI and where the field is headed in the coming years.

Diffusion MRI signal cumulants and hepatocyte microstructure at fixed diffusion time: Insights from simulations, 9.4T imaging, and histology

Purpose

Relationships between diffusion‐weighted MRI signals and hepatocyte microstructure were investigated to inform liver diffusion MRI modeling, focusing on the following question: Can cell size and diffusivity be estimated at fixed diffusion time, realistic SNR, and negligible contribution from extracellular/extravascular water and exchange?

Methods

Monte Carlo simulations were performed within synthetic hepatocytes for varying cell size/diffusivity L/D0, and clinical protocols (single diffusion encoding; maximum b‐value: {1000, 1500, 2000} s/mm²; 5 unique gradient duration/separation pairs; SNR = {∞, 100, 80, 40, 20}), accounting for heterogeneity in (D0,L) and perfusion contamination. Diffusion (D) and kurtosis (K) coefficients were calculated, and relationships between (D0,L) and (D,K) were visualized. Functions mapping (D,K) to (D0,L) were computed to predict unseen (D0,L) values, tested for their ability to classify discrete cell‐size contrasts, and deployed on 9.4T ex vivo MRI‐histology data of fixed mouse livers

Results

Relationships between (D,K) and (D0,L) are complex and depend on the diffusion encoding. Functions mapping D,K to (D0,L) captures salient characteristics of D0(D,K) and L(D,K) dependencies. Mappings are not always accurate, but they enable just under 70% accuracy in a three‐class cell‐size classification task (for SNR = 20, bmax = 1500 s/mm², δ = 20 ms, and Δ = 75 ms). MRI detects cell‐size contrasts in the mouse livers that are confirmed by histology, but overestimates the largest cell sizes.

Conclusion

Salient information about liver cell size and diffusivity may be retrieved from minimal diffusion encodings at fixed diffusion time, in experimental conditions and pathological scenarios for which extracellular, extravascular water and exchange are negligible.

Diffusion time dependence, power-law scaling, and exchange in gray matter

Characterizing neural tissue microstructure is a critical goal for future neuroimaging. Diffusion MRI (dMRI) provides contrasts that reflect diffusing spins’ interactions with myriad microstructural features of biological systems. However, the specificity of dMRI remains limited due to the ambiguity of its signals vis-à-vis the underlying microstructure. To improve specificity, biophysical models of white matter (WM) typically express dMRI signals according to the Standard Model (SM) and have more recently in gray matter (GM) taken spherical compartments into account (the SANDI model) in attempts to represent cell soma. The validity of the assumptions underlying these models, however, remains largely undetermined, especially in GM. To validate these assumptions experimentally, observing their unique, functional properties, such as the b^−1/2 power-law associated with one-dimensional diffusion, has emerged as a fruitful strategy. The absence of this signature in GM, in turn, has been explained by neurite water exchange, non-linear morphology, and/or by obscuring soma signal contributions. Here, we present diffusion simulations in realistic neurons demonstrating that curvature and branching does not destroy the stick power-law behavior in impermeable neurites, but also that their signal is drowned by the soma signal under typical experimental conditions. Nevertheless, by studying the GM dMRI signal's behavior as a function of diffusion weighting as well as time, we identify an attainable experimental regime in which the neurite signal dominates. Furthermore, we find that exchange-driven time dependence produces a signal behavior opposite to that which would be expected from restricted diffusion, thereby providing a functional signature that disambiguates the two effects. We present data from dMRI experiments in ex vivo rat brain at ultrahigh field of 16.4T and observe a time dependence that is consistent with substantial exchange but also with a GM stick power-law. The first finding suggests significant water exchange between neurites and the extracellular space while the second suggests a small sub-population of impermeable neurites. To quantify these observations, we harness the Kärger exchange model and incorporate the corresponding signal time dependence in the SM and SANDI models.

Nonparametric D-R1-R2 distribution MRI of the living human brain

Diffusion-relaxation correlation NMR can simultaneously characterize both the microstructure and the local chemical composition of complex samples that contain multiple populations of water. Recent developments on tensor-valued diffusion encoding and Monte Carlo inversion algorithms have made it possible to transfer diffusion-relaxation correlation NMR from small-bore scanners to clinical MRI systems. Initial studies on clinical MRI systems employed 5D D-R₁ and D-R₂ correlation to characterize healthy brain in vivo. However, these methods are subject to an inherent bias that originates from not including R₂ or R₁ in the analysis, respectively. This drawback can be remedied by extending the concept to 6D D-R₁-R₂ correlation. In this work, we present a sparse acquisition protocol that records all data necessary for in vivo 6D D-R₁-R₂ correlation MRI across 633 individual measurements within 25 min-a time frame comparable to previous lower-dimensional acquisition protocols. The data were processed with a Monte Carlo inversion algorithm to obtain nonparametric 6D D-R₁-R₂ distributions. We validated the reproducibility of the method in repeated measurements of healthy volunteers. For a post-therapy glioblastoma case featuring cysts, edema, and partially necrotic remains of tumor, we present representative single-voxel 6D distributions, parameter maps, and artificial contrasts over a wide range of diffusion-, R₁-, and R₂-weightings based on the rich information contained in the D-R₁-R₂ distributions.

Cross-scanner reproducibility and harmonization of a diffusion MRI structural brain network: A traveling subject study of multi-b acquisition

Characterization of brain networks by diffusion MRI (dMRI) has rapidly evolved, and there are ongoing movements toward data sharing and multi-center studies. To extract meaningful information from multi-center data, methods to correct for the bias caused by scanner differences, that is, harmonization, are urgently needed. In this work, we report the cross-scanner differences in structural network analyses using data from nine traveling subjects (four males and five females, 21-49 years-old) who underwent scanning using four 3T scanners (public database available from the Brain/MINDS Beyond Human Brain MRI project (http://mriportal.umin.jp/)). The reliability and reproducibility were compared to those of data from another set of four subjects (all males, 29-42 years-old) who underwent scan-rescan (interval, 105-147 days) with the same scanner as well as scan-rescan data from the Human Connectome Project database. The results demonstrated that the reliability of the edge weights and graph theory metrics was lower for data including different scanners, compared to the scan-rescan with the same scanner. Besides, systematic differences between scanners were observed, indicating the risk of bias in comparing networks obtained from different scanners directly. We further demonstrate that it is feasible to reduce inter-scanner variabilities while preserving the inter-subject differences among healthy individuals by modeling the scanner effects at the level of network matrices, when traveling-subject data are available for calibration between scanners. The present data and results are expected to serve as a basis for developing and evaluating novel harmonization methods.

Improved unsupervised physics-informed deep learning for intravoxel incoherent motion modeling and evaluation in pancreatic cancer patients

PURPOSE

Earlier work showed that IVIM-NETorig , an unsupervised physics-informed deep neural network, was faster and more accurate than other state-of-the-art intravoxel-incoherent motion (IVIM) fitting approaches to diffusion-weighted imaging (DWI). This study presents a substantially improved version, IVIM-NEToptim , and characterizes its superior performance in pancreatic cancer patients.

METHOD

In simulations (signal-to-noise ratio [SNR] = 20), the accuracy, independence, and consistency of IVIM-NET were evaluated for combinations of hyperparameters (fit S0, constraints, network architecture, number of hidden layers, dropout, batch normalization, learning rate), by calculating the normalized root-mean-square error (NRMSE), Spearman's ρ, and the coefficient of variation (CVNET ), respectively. The best performing network, IVIM-NEToptim was compared to least squares (LS) and a Bayesian approach at different SNRs. IVIM-NEToptim 's performance was evaluated in an independent dataset of 23 patients with pancreatic ductal adenocarcinoma. Fourteen of the patients received no treatment between two repeated scan sessions and nine received chemoradiotherapy between the repeated sessions. Intersession within-subject standard deviations (wSD) and treatment-induced changes were assessed.

RESULTS

In simulations (SNR = 20), IVIM-NEToptim outperformed IVIM-NETorig in accuracy (NRMSE(D) = 0.177 vs 0.196; NMRSE(f) = 0.220 vs 0.267; NMRSE(D*) = 0.386 vs 0.393), independence (ρ(D*, f) = 0.22 vs 0.74), and consistency (CVNET (D) = 0.013 vs 0.104; CVNET (f) = 0.020 vs 0.054; CVNET (D*) = 0.036 vs 0.110). IVIM-NEToptim showed superior performance to the LS and Bayesian approaches at SNRs < 50. In vivo, IVIM-NEToptim showed significantly less noisy parameter maps with lower wSD for D and f than the alternatives. In the treated cohort, IVIM-NEToptim detected the most individual patients with significant parameter changes compared to day-to-day variations.

CONCLUSION

IVIM-NEToptim is recommended for accurate, informative, and consistent IVIM fitting to DWI data.

Training data distribution significantly impacts the estimation of tissue microstructure with machine learning

PURPOSE

Supervised machine learning (ML) provides a compelling alternative to traditional model fitting for parameter mapping in quantitative MRI. The aim of this work is to demonstrate and quantify the effect of different training data distributions on the accuracy and precision of parameter estimates when supervised ML is used for fitting.

METHODS

We fit a two- and three-compartment biophysical model to diffusion measurements from in-vivo human brain, as well as simulated diffusion data, using both traditional model fitting and supervised ML. For supervised ML, we train several artificial neural networks, as well as random forest regressors, on different distributions of ground truth parameters. We compare the accuracy and precision of parameter estimates obtained from the different estimation approaches using synthetic test data.

RESULTS

When the distribution of parameter combinations in the training set matches those observed in healthy human data sets, we observe high precision, but inaccurate estimates for atypical parameter combinations. In contrast, when training data is sampled uniformly from the entire plausible parameter space, estimates tend to be more accurate for atypical parameter combinations but may have lower precision for typical parameter combinations.

CONCLUSION

This work highlights that estimation of model parameters using supervised ML depends strongly on the training-set distribution. We show that high precision obtained using ML may mask strong bias, and visual assessment of the parameter maps is not sufficient for evaluating the quality of the estimates.

Optimized bias and signal inference in diffusion-weighted image analysis (OBSIDIAN)

PURPOSE

Correction of Rician signal bias in magnitude MR images.

METHODS

A model-based, iterative fitting procedure is used to simultaneously estimate true signal and underlying Gaussian noise with standard deviation σ g on a pixel-by-pixel basis in magnitude MR images. A precomputed function that relates absolute residuals between measured signals and model fit to σ g is used to iteratively estimate σ g . The feasibility of the method is evaluated and compared to maximum likelihood estimation (MLE) for diffusion signal decay simulations and diffusion-weighted images of the prostate considering 21 linearly spaced b-values from 0 to 3000 s/mm² . A multidirectional analysis was performed with publically available brain data.

RESULTS

Model simulations show that the Rician bias correction algorithm is fast, with an accuracy and precision that is on par to model-based MLE and direct fitting in the case of pure Gaussian noise. Increased accuracy in parameter prediction in a low signal-to-noise ratio (SNR) scenario is ideally achieved by using a composite of multiple signal decays from neighboring voxels as input for the algorithm. For patient data, good agreement with high SNR reference data of diffusion in prostate is achieved.

CONCLUSIONS

OBSIDIAN is a novel, alternative, simple to implement approach for rapid Rician bias correction applicable in any case where differences between true signal decay and underlying model function can be considered negligible in comparison to noise. The proposed composite fitting approach permits accurate parameter estimation even in typical clinical scenarios with low SNR, which significantly simplifies comparison of complex diffusion parameters among studies.

Removal of partial Fourier-induced Gibbs (RPG) ringing artifacts in MRI

PURPOSE

To investigate and remove Gibbs-ringing artifacts caused by partial Fourier (PF) acquisition and zero filling interpolation in MRI data.

THEORY AND METHODS

Gibbs ringing of fully sampled data, leading to oscillations around tissue boundaries, is caused by the symmetric truncation of k-space. Such ringing can be removed by conventional methods, with the local subvoxel shifts method being the state-of-the-art. However, the asymmetric truncation of k-space in routinely used PF acquisitions leads to additional ringings of wider intervals in the PF sampling dimension that cannot be corrected solely based on magnitude images reconstructed via zero filling. Here, we develop a pipeline for the Removal of PF-induced Gibbs ringing (RPG) to remove ringing patterns of different periods by applying the conventional method twice. The proposed pipeline is validated on numerical phantoms, demonstrated on in vivo diffusion MRI measurements, and compared with the conventional method and neural network-based approach.

RESULTS

For PF = 7/8 and 6/8, Gibbs-ringings and subsequent bias in diffusion metrics induced by PF acquisition and zero filling are robustly removed by using the proposed RPG pipeline. For PF = 5/8, however, ringing removal via RPG leads to excessive image blurring due to the interplay of image phase and convolution kernel.

CONCLUSIONS

RPG corrects Gibbs-ringing artifacts in magnitude images of PF acquired data and reduces the bias in quantitative MR metrics. Considering the benefit of PF acquisition and the feasibility of ringing removal, we suggest applying PF = 6/8 when PF acquisition is necessary.

PreQual: An automated pipeline for integrated preprocessing and quality assurance of diffusion weighted MRI images

PURPOSE

Diffusion weighted MRI imaging (DWI) is often subject to low signal-to-noise ratios (SNRs) and artifacts. Recent work has produced software tools that can correct individual problems, but these tools have not been combined with each other and with quality assurance (QA). A single integrated pipeline is proposed to perform DWI preprocessing with a spectrum of tools and produce an intuitive QA document.

METHODS

The proposed pipeline, built around the FSL, MRTrix3, and ANTs software packages, performs DWI denoising; inter-scan intensity normalization; susceptibility-, eddy current-, and motion-induced artifact correction; and slice-wise signal drop-out imputation. To perform QA on the raw and preprocessed data and each preprocessing operation, the pipeline documents qualitative visualizations, quantitative plots, gradient verifications, and tensor goodness-of-fit and fractional anisotropy analyses.

RESULTS

Raw DWI data were preprocessed and quality checked with the proposed pipeline and demonstrated improved SNRs; physiologic intensity ratios; corrected susceptibility-, eddy current-, and motion-induced artifacts; imputed signal-lost slices; and improved tensor fits. The pipeline identified incorrect gradient configurations and file-type conversion errors and was shown to be effective on externally available datasets.

CONCLUSIONS

The proposed pipeline is a single integrated pipeline that combines established diffusion preprocessing tools from major MRI-focused software packages with intuitive QA.

Rapid microstructural plasticity in the cortical semantic network following a short language learning session

Despite the clear importance of language in our life, our vital ability to quickly and effectively learn new words and meanings is neurobiologically poorly understood. Conventional knowledge maintains that language learning—especially in adulthood—is slow and laborious. Furthermore, its structural basis remains unclear. Even though behavioural manifestations of learning are evident near instantly, previous neuroimaging work across a range of semantic categories has largely studied neural changes associated with months or years of practice. Here, we address rapid neuroanatomical plasticity accompanying new lexicon acquisition, specifically focussing on the learning of action-related language, which has been linked to the brain’s motor systems. Our results show that it is possible to measure and to externally modulate (using transcranial magnetic stimulation (TMS) of motor cortex) cortical microanatomic reorganisation after mere minutes of new word learning. Learning-induced microstructural changes, as measured by diffusion kurtosis imaging (DKI) and machine learning-based analysis, were evident in prefrontal, temporal, and parietal neocortical sites, likely reflecting integrative lexico-semantic processing and formation of new memory circuits immediately during the learning tasks. These results suggest a structural basis for the rapid neocortical word encoding mechanism and reveal the causally interactive relationship of modal and associative brain regions in supporting learning and word acquisition.

This combined neuroimaging and brain stimulation study reveals rapid and distributed microstructural plasticity after a single immersive language learning session, demonstrating the causal relevance of the motor cortex in encoding the meaning of novel action words.

Beyond the diffusion standard model in fixed rat spinal cord with combined linear and planar encoding

Information about tissue on the microscopic and mesoscopic scales can be accessed by modelling diffusion MRI signals, with the aim of extracting microstructure-specific biomarkers. The standard model (SM) of diffusion, currently the most broadly adopted microstructural model, describes diffusion in white matter (WM) tissues by two Gaussian components, one of which has zero radial diffusivity, to represent diffusion in intra- and extra-axonal water, respectively. Here, we reappraise these SM assumptions by collecting comprehensive double diffusion encoded (DDE) MRI data with both linear and planar encodings, which was recently shown to substantially enhance the ability to estimate SM parameters. We find however, that the SM is unable to account for data recorded in fixed rat spinal cord at an ultrahigh field of 16.4 T, suggesting that its underlying assumptions are violated in our experimental data. We offer three model extensions to mitigate this problem: first, we generalize the SM to accommodate finite radii (axons) by releasing the constraint of zero radial diffusivity in the intra-axonal compartment. Second, we include intracompartmental kurtosis to account for non-Gaussian behaviour. Third, we introduce an additional (third) compartment. The ability of these models to account for our experimental data are compared based on parameter feasibility and Bayesian information criterion. Our analysis identifies the three-compartment description as the optimal model. The third compartment exhibits slow diffusion with a minor but non-negligible signal fraction (∼12%). We demonstrate how failure to take the presence of such a compartment into account severely misguides inferences about WM microstructure. Our findings bear significance for microstructural modelling at large and can impact the interpretation of biomarkers extracted from the standard model of diffusion.

Heterogeneity of multiple sclerosis lesions in fast diffusional kurtosis imaging

Background

Mean kurtosis (MK), one of the parameters derived from diffusion kurtosis imaging (DKI), has shown increased sensitivity to tissue microstructure damage in several neurological disorders.

Methods

Thirty-seven patients with relapsing-remitting MS and eleven healthy controls (HC) received brain imaging on a 3T MR scanner, including a fast DKI sequence. MK and mean diffusivity (MD) were measured in the white matter of HC, normal-appearing white matter (NAWM) of MS patients, contrast-enhancing lesions (CE-L), FLAIR lesions (FLAIR-L) and black holes (BH).

Results

Overall 1529 lesions were analyzed, including 30 CE-L, 832 FLAIR-L and 667 BH. Highest MK values were obtained in the white matter of HC (0.814 ± 0.129), followed by NAWM (0.724 ± 0.137), CE-L (0.619 ± 0.096), FLAIR-L (0.565 ± 0.123) and BH (0.549 ± 0.12). Lowest MD values were obtained in the white matter of HC (0.747 ± 0.068 10⁻³mm²/sec), followed by NAWM (0.808 ± 0.163 10⁻³mm²/sec), CE-L (0.853 ± 0.211 10⁻³mm²/sec), BH (0.957 ± 0.304 10⁻³mm²/sec) and FLAIR-L (0.976 ± 0.35 10⁻³mm²/sec). While MK differed significantly between CE-L and non-enhancing lesions, MD did not.

Conclusion

MK adds predictive value to differentiate between MS lesions and might provide further information about diffuse white matter injury and lesion microstructure.

Data Preparation Protocol for Low Signal-to-Noise Ratio Fluorine-19 MRI

Fluorine-19 MRI shows great promise for a wide range of applications including renal imaging, yet the typically low signal-to-noise ratios and sparse signal distribution necessitate a thorough data preparation.This chapter describes a general data preparation workflow for fluorine MRI experiments. The main processing steps are: (1) estimation of noise level, (2) correction of noise-induced bias and (3) background subtraction. The protocol is supplemented by an example script and toolbox available online.This chapter is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers. This analysis protocol chapter is complemented by two separate chapters describing the basic concept and experimental procedure.

The sensitivity of diffusion MRI to microstructural properties and experimental factors

Diffusion MRI is a non-invasive technique to study brain microstructure. Differences in the microstructural properties of tissue, including size and anisotropy, can be represented in the signal if the appropriate method of acquisition is used. However, to depict the underlying properties, special care must be taken when designing the acquisition protocol as any changes in the procedure might impact on quantitative measurements. This work reviews state-of-the-art methods for studying brain microstructure using diffusion MRI and their sensitivity to microstructural differences and various experimental factors. Microstructural properties of the tissue at a micrometer scale can be linked to the diffusion signal at a millimeter-scale using modeling. In this paper, we first give an introduction to diffusion MRI and different encoding schemes. Then, signal representation-based methods and multi-compartment models are explained briefly. The sensitivity of the diffusion MRI signal to the microstructural components and the effects of curvedness of axonal trajectories on the diffusion signal are reviewed. Factors that impact on the quality (accuracy and precision) of derived metrics are then reviewed, including the impact of random noise, and variations in the acquisition parameters (i.e., number of sampled signals, b-value and number of acquisition shells). Finally, yet importantly, typical approaches to deal with experimental factors are depicted, including unbiased measures and harmonization. We conclude the review with some future directions and recommendations on this topic.

In vivo observation and biophysical interpretation of time-dependent diffusion in human cortical gray matter

The dependence of the diffusion MRI signal on the diffusion time t is a hallmark of tissue microstructure at the scale of the diffusion length. Here we measure the time-dependence of the mean diffusivity D(t) and mean kurtosis K(t) in cortical gray matter and in 25 ​gray matter sub-regions, in 10 healthy subjects. Significant diffusivity and kurtosis time-dependence is observed for t=21.2-100 ​ms, and is characterized by a power-law tail ∼t-ϑ with dynamical exponent ϑ. To interpret our measurements, we systematize the relevant scenarios and mechanisms for diffusion time-dependence in the brain. Using the effective medium theory formalism, we derive an exact relation between the power-law tails in D(t) and K(t). The estimated dynamical exponent ϑ≃1/2 in both D(t) and K(t) is consistent with one-dimensional diffusion in the presence of randomly positioned restrictions along neurites. We analyze the short-range disordered statistics of synapses on axon collaterals in the cortex, and perform one-dimensional Monte Carlo simulations of diffusion restricted by permeable barriers with a similar randomness in their placement, to confirm the ϑ=1/2 exponent. In contrast, the Kärger model of exchange is less consistent with the data since it does not capture the diffusivity time-dependence, and the estimated exchange time from K(t) falls below our measured t-range. Although we cannot exclude exchange as a contributing factor, we argue that structural disorder along neurites is mainly responsible for the observed time-dependence of diffusivity and kurtosis. Our observation and theoretical interpretation of the t-1/2 tail in D(t) and K(t) altogether establish the sensitivity of a macroscopic MRI signal to micrometer-scale structural heterogeneities along neurites in human gray matter in vivo.

Multi-parametric quantitative in vivo spinal cord MRI with unified signal readout and image denoising

Multi-parametric quantitative MRI (qMRI) of the spinal cord is a promising non-invasive tool to probe early microstructural damage in neurological disorders. It is usually performed in vivo by combining acquisitions with multiple signal readouts, which exhibit different thermal noise levels, geometrical distortions and susceptibility to physiological noise. This ultimately hinders joint multi-contrast modelling and makes the geometric correspondence of parametric maps challenging. We propose an approach to overcome these limitations, by implementing state-of-the-art microstructural MRI of the spinal cord with a unified signal readout in vivo (i.e. with matched spatial encoding parameters across a range of imaging contrasts). We base our acquisition on single-shot echo planar imaging with reduced field-of-view, and obtain data from two different vendors (vendor 1: Philips Achieva; vendor 2: Siemens Prisma). Importantly, the unified acquisition allows us to compare signal and noise across contrasts, thus enabling overall quality enhancement via multi-contrast image denoising methods. As a proof-of-concept, here we provide a demonstration with one such method, known as Marchenko-Pastur (MP) Principal Component Analysis (PCA) denoising. MP-PCA is a singular value (SV) decomposition truncation approach that relies on redundant acquisitions, i.e. such that the number of measurements is large compared to the number of components that are maintained in the truncated SV decomposition. Here we used in vivo and synthetic data to test whether a unified readout enables more efficient MP-PCA denoising of less redundant acquisitions, since these can be denoised jointly with more redundant ones. We demonstrate that a unified readout provides robust multi-parametric maps, including diffusion and kurtosis tensors from diffusion MRI, myelin metrics from two-pool magnetisation transfer, and T1 and T2 from relaxometry. Moreover, we show that MP-PCA improves the quality of our multi-contrast acquisitions, since it reduces the coefficient of variation (i.e. variability) by up to 17% for mean kurtosis, 8% for bound pool fraction (myelin-sensitive), and 13% for T1, while enabling more efficient denoising of modalities limited in redundancy (e.g. relaxometry). In conclusion, multi-parametric spinal cord qMRI with unified readout is feasible and provides robust microstructural metrics with matched resolution and distortions, whose quality benefits from multi-contrast denoising methods such as MP-PCA.

Automated characterization of noise distributions in diffusion MRI data

Knowledge of the noise distribution in magnitude diffusion MRI images is the centerpiece to quantify uncertainties arising from the acquisition process. The use of parallel imaging methods, the number of receiver coils and imaging filters applied by the scanner, amongst other factors, dictate the resulting signal distribution. Accurate estimation beyond textbook Rician or noncentral chi distributions often requires information about the acquisition process (e.g., coils sensitivity maps or reconstruction coefficients), which is usually not available. We introduce two new automated methods using the moments and maximum likelihood equations of the Gamma distribution to estimate noise distributions as they explicitly depend on the number of coils, making it possible to estimate all unknown parameters using only the magnitude data. A rejection step is used to make the framework automatic and robust to artifacts. Simulations using stationary and spatially varying noncentral chi noise distributions were created for two diffusion weightings with SENSE or GRAPPA reconstruction and 8, 12 or 32 receiver coils. Furthermore, MRI data of a water phantom with different combinations of parallel imaging were acquired on a 3T Philips scanner along with noise-only measurements. Finally, experiments on freely available datasets from a single subject acquired on a 3T GE scanner are used to assess reproducibility when limited information about the acquisition protocol is available. Additionally, we demonstrated the applicability of the proposed methods for a bias correction and denoising task on an in vivo dataset acquired on a 3T Siemens scanner. A generalized version of the bias correction framework for non integer degrees of freedom is also introduced. The proposed framework is compared with three other algorithms with datasets from three vendors, employing different reconstruction methods. Simulations showed that assuming a Rician distribution can lead to misestimation of the noise distribution in parallel imaging. Results on the acquired datasets showed that signal leakage in multiband can also lead to a misestimation of the noise distribution. Repeated acquisitions of in vivo datasets show that the estimated parameters are stable and have lower variability than compared methods. Results for the bias correction and denoising task show that the proposed methods reduce the appearance of noise at high b-value. The proposed algorithms herein can estimate both parameters of the noise distribution automatically, are robust to signal leakage artifacts and perform best when used on acquired noise maps.

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.

Denoise magnitude diffusion magnetic resonance images via variance-stabilizing transformation and optimal singular-value manipulation

Although shown to have a great utility for a wide range of neuroscientific and clinical applications, diffusion-weighted magnetic resonance imaging (dMRI) faces a major challenge of low signal-to-noise ratio (SNR), especially when pushing the spatial resolution for improved delineation of brain's fine structure or increasing the diffusion weighting for increased angular contrast or both. Here, we introduce a comprehensive denoising framework for denoising magnitude dMRI. The framework synergistically combines the variance stabilizing transform (VST) with optimal singular value manipulation. The purpose of VST is to transform the Rician data to Gaussian-like data so that an asymptotically optimal singular value manipulation strategy tailored for Gaussian data can be used. The output of the framework is the estimated underlying diffusion signal for each voxel in the image domain. The usefulness of the proposed framework for denoising magnitude dMRI is demonstrated using both simulation and real-data experiments. Our results show that the proposed denoising framework can significantly improve SNR across the entire brain, leading to substantially enhanced performances for estimating diffusion tensor related indices and for resolving crossing fibers when compared to another competing method. More encouragingly, the proposed method when used to denoise a single average of 7 ​Tesla Human Connectome Project-style diffusion acquisition provided comparable performances relative to those achievable with ten averages for resolving multiple fiber populations across the brain. As such, the proposed denoising method is expected to have a great utility for high-quality, high-resolution whole-brain dMRI, desirable for many neuroscientific and clinical applications.

Nonivasive quantification of axon radii using diffusion MRI

Axon caliber plays a crucial role in determining conduction velocity and, consequently, in the timing and synchronization of neural activation. Noninvasive measurement of axon radii could have significant impact on the understanding of healthy and diseased neural processes. Until now, accurate axon radius mapping has eluded in vivo neuroimaging, mainly due to a lack of sensitivity of the MRI signal to micron-sized axons. Here, we show how - when confounding factors such as extra-axonal water and axonal orientation dispersion are eliminated - heavily diffusion-weighted MRI signals become sensitive to axon radii. However, diffusion MRI is only capable of estimating a single metric, the effective radius, representing the entire axon radius distribution within a voxel that emphasizes the larger axons. Our findings, both in rodents and humans, enable noninvasive mapping of critical information on axon radii, as well as resolve the long-standing debate on whether axon radii can be quantified.

An enhanced adaptive non-local means algorithm for Rician noise reduction in magnetic resonance brain images

Background

The Rician noise formed in magnetic resonance (MR) imaging greatly reduced the accuracy and reliability of subsequent analysis, and most of the existing denoising methods are suitable for Gaussian noise rather than Rician noise. Aiming to solve this problem, we proposed fuzzy c-means and adaptive non-local means (FANLM), which combined the adaptive non-local means (NLM) with fuzzy c-means (FCM), as a novel method to reduce noise in the study.

Method

The algorithm chose the optimal size of search window automatically based on the noise variance which was estimated by the improved estimator of the median absolute deviation (MAD) for Rician noise. Meanwhile, it solved the problem that the traditional NLM algorithm had to use a fixed size of search window. Considering the distribution characteristics for each pixel, we designed three types of search window sizes as large, medium and small instead of using a fixed size. In addition, the combination with the FCM algorithm helped to achieve better denoising effect since the improved the FCM algorithm divided the membership degrees of images and introduced the morphological reconstruction to preserve the image details.

Results

The experimental results showed that the proposed algorithm (FANLM) can effectively remove the noise. Moreover, it had the highest peak signal-noise ratio (PSNR) and structural similarity (SSIM), compared with other three methods: non-local means (NLM), linear minimum mean square error (LMMSE) and undecimated wavelet transform (UWT). Using the FANLM method, the image details can be well preserved with the noise being mostly removed.

Conclusion

Compared with the traditional denoising methods, the experimental results showed that the proposed approach effectively suppressed the noise and the edge details were well retained. However, the FANLM method took an average of 13 s throughout the experiment, and its computational cost was not the shortest. Addressing these can be part of our future research.

Denoising of Diffusion MRI Data via Graph Framelet Matching in x-q Space

Diffusion magnetic resonance imaging (DMRI) suffers from lower signal-to-noise-ratio (SNR) due to MR signal attenuation associated with the motion of water molecules. To improve SNR, the non-local means (NLM) algorithm has demonstrated state-of-the-art performance in noise reduction. However, existing NLM algorithms do not take into account explicitly the fact that DMRI signal can vary significantly with local fiber orientations. Applying NLM naïvely can hence blur subtle structures and aggravate partial volume effects. To overcome this limitation, we improve NLM by performing neighborhood matching in non-flat domains and removing noise with information from both x -space (spatial domain) and q -space (wavevector domain). Specifically, we first encode the q -space sampling domain using a graph. We then perform graph framelet transforms to extract robust rotation-invariant features for each sampling point in x-q space. The resulting features are employed for robust neighborhood matching to locate recurrent information. Finally, we remove noise via an NLM framework. To adapt to the various types of noise in multi-coil MR imaging, we transform the signal before denoising so that it is Gaussian-distributed, allowing noise removal to be carried out in an unbiased manner. Our method is able to more effectively locate recurrent information in white matter structures with different orientations, avoiding the blurring effects caused by naïvely applying NLM. Experiments on synthetic, repetitively-acquired, and infant DMRI data demonstrate that our method is able to preserve subtle structures while effectively removing noise.

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.

Towards an optimised processing pipeline for diffusion magnetic resonance imaging data: Effects of artefact corrections on diffusion metrics and their age associations in UK Biobank

Increasing interest in the structural and functional organisation of the human brain encourages the acquisition of big data sets comprising multiple neuroimaging modalities, often accompanied by additional information obtained from health records, cognitive tests, biomarkers and genotypes. Diffusion weighted magnetic resonance imaging data enables a range of promising imaging phenotypes probing structural connections as well as macroanatomical and microstructural properties of the brain. The reliability and biological sensitivity and specificity of diffusion data depend on processing pipeline. A state-of-the-art framework for data processing facilitates cross-study harmonisation and reduces pipeline-related variability. Using diffusion magnetic resonance imaging (MRI) data from 218 individuals in the UK Biobank, we evaluate the effects of different processing steps that have been suggested to reduce imaging artefacts and improve reliability of diffusion metrics. In lack of a ground truth, we compared diffusion metric sensitivity to age between pipelines. By comparing distributions and age sensitivity of the resulting diffusion metrics based on different approaches (diffusion tensor imaging, diffusion kurtosis imaging and white matter tract integrity), we evaluate a general pipeline comprising seven postprocessing blocks: noise correction; Gibbs ringing correction; evaluation of field distortions; susceptibility, eddy-current and motion-induced distortion corrections; bias field correction; spatial smoothing and final diffusion metric estimations. Based on this evaluation, we suggest an optimised processing pipeline for diffusion weighted MRI data.

Improved accuracy of apparent diffusion coefficient quantification using a fully automatic noise bias compensation method: Preliminary evaluation in prostate diffusion weighted imaging

Noise in diffusion magnetic resonance imaging can introduce bias in apparent diffusion coefficient (ADC) quantification. Previous studies proposed methods that are site-specific techniques as research tools with limited availability and typically require manual intervention, not completely ready to use in the clinical environment. The purpose of this study was to develop a fully automatic computational method to correct noise bias in ADC quantification and perform a preliminary evaluation in the clinical prostate diffusion weighted imaging (DWI). Using a pseudo replica approach for the noise map calculation as well as a direct mapping and a stepwise Chebychev polynomial modelling approach for the ADC fitting, a fully automatic noise-bias-compensated ADC calculation method was proposed and implemented both on the scanner and offline. The proposed method was validated in a computer simulation and a standardized diffusion phantom with ground-truth values. Two in vivo studies were performed to evaluate the proposed method in the clinical environment. The first in vivo study performed acquisitions using a clinically routine prostate DWI protocol on 29 subjects to evaluate the consistency between simulated and empirical results. In the second in vivo study, prostate ADC values of 14 subjects were compared between data acquired with external coils only and reconstructed with the proposed method vs. acquired with external combined with endorectal coils and reconstructed with the conventional method. In statistical analyses, p < 0.05 was regarded as significantly different. In the computer simulation, the proposed method showed smaller error percentage than the other methods and was significantly different (p < 2.2 × 10-16). With low signal-to-noise ratio (SNR), the conventional method underestimated ADC values compared to the ground truth values of the diffusion phantom, while the results of the proposed method were more consistent with the ground truth values. Statistical analyses showed no significant differences between measured and simulated results in the first in vivo study (p = 0.5618). Data from the second in vivo study showed that agreement between ADC measured with external coils only and combined coils was improved for the proposed method (mean bias: 0.04 × 10-3 mm2/s, 95% confidence interval (CI) = [-0.01, 0.09] × 10-3 mm2/s, p = 0.187), compared to the conventional method (mean bias: -0.12 × 10-3 mm2/s, 95% CI = [-0.17, -0.06] × 10-3 mm2/s, p < 0.0001). The proposed method compensates noise bias in low-SNR diffusion-weighted acquisitions and results show improved ADC quantification accuracy in the prostate. This method may be suitable for both clinical imaging and research utilizing ADC quantification.