Validation of deep learning techniques for quality augmentation in diffusion MRI for clinical studies.

The objective of this study is to evaluate the efficacy of deep learning (DL) techniques in improving the quality of diffusion MRI (dMRI) data in clinical applications. The study aims to determine whether the use of artificial intelligence (AI) methods in medical images may result in the loss of critical clinical information and/or the appearance of false information. To assess this, the focus was on the angular resolution of dMRI and a clinical trial was conducted on migraine, specifically between episodic and chronic migraine patients. The number of gradient directions had an impact on white matter analysis results, with statistically significant differences between groups being drastically reduced when using 21 gradient directions instead of the original 61. Fourteen teams from different institutions were tasked to use DL to enhance three diffusion metrics (FA, AD and MD) calculated from data acquired with 21 gradient directions and a b-value of 1000 s/mm2. The goal was to produce results that were comparable to those calculated from 61 gradient directions. The results were evaluated using both standard image quality metrics and Tract-Based Spatial Statistics (TBSS) to compare episodic and chronic migraine patients. The study results suggest that while most DL techniques improved the ability to detect statistical differences between groups, they also led to an increase in false positive. The results showed that there was a constant growth rate of false positives linearly proportional to the new true positives, which highlights the risk of generalization of AI-based tasks when assessing diverse clinical cohorts and training using data from a single group. The methods also showed divergent performance when replicating the original distribution of the data and some exhibited significant bias. In conclusion, extreme caution should be exercised when using AI methods for harmonization or synthesis in clinical studies when processing heterogeneous data in clinical studies, as important information may be altered, even when global metrics such as structural similarity or peak signal-to-noise ratio appear to suggest otherwise.

Spherical means-based free-water volume fraction from diffusion MRI increases non-linearly with age in the white matter of the healthy human brain.

The term free-water volume fraction (FWVF) refers to the signal fraction that could be found as the cerebrospinal fluid of the brain, which has been demonstrated as a sensitive measure that correlates with cognitive performance and various neuropathological processes. It can be quantified by properly fitting the isotropic component of the magnetic resonance (MR) signal in diffusion-sensitized sequences. Using N=287 healthy subjects (178F/109M) aged 25-94, this study examines in detail the evolution of the FWVF obtained with the spherical means technique from multi-shell acquisitions in the human brain white matter across the adult lifespan, which has been previously reported to exhibit a positive trend when estimated from single-shell data using the bi-tensor signal representation. We found evidence of a noticeably non-linear gain after the sixth decade of life, with a region-specific variate and varying change rate of the spherical means-based multi-shell FWVF parameter with age, at the same time, a heteroskedastic pattern across the adult lifespan is suggested. On the other hand, the FW corrected diffusion tensor imaging (DTI) leads to a region-dependent flattened age-related evolution of the mean diffusivity (MD) and fractional anisotropy (FA), along with a considerable reduction in their variability, as compared to the studies conducted over the standard (single-component) DTI. This way, our study provides a new perspective on the trajectory-based assessment of the brain and explains the conceivable reason for the variations observed in FA and MD parameters across the lifespan with previous studies under the standard diffusion tensor imaging.

Viability of AMURA biomarkers from single-shell diffusion MRI in clinical studies.

Diffusion Tensor Imaging (DTI) is the most employed method to assess white matter properties using quantitative parameters derived from diffusion MRI, but it presents known limitations that restrict the evaluation of complex structures. The objective of this study was to validate the reliability and robustness of complementary diffusion measures extracted with a novel approach, Apparent Measures Using Reduced Acquisitions (AMURA), with a typical diffusion MRI acquisition from a clinical context in comparison with DTI with application to clinical studies. Fifty healthy controls, 51 episodic migraine and 56 chronic migraine patients underwent single-shell diffusion MRI. Four DTI-based and eight AMURA-based parameters were compared between groups with tract-based spatial statistics to establish reference results. On the other hand, following a region-based analysis, the measures were assessed for multiple subsamples with diverse reduced sample sizes and their stability was evaluated with the coefficient of quartile variation. To assess the discrimination power of the diffusion measures, we repeated the statistical comparisons with a region-based analysis employing reduced sample sizes with diverse subsets, decreasing 10 subjects per group for consecutive reductions, and using 5,001 different random subsamples. For each sample size, the stability of the diffusion descriptors was evaluated with the coefficient of quartile variation. AMURA measures showed a greater number of statistically significant differences in the reference comparisons between episodic migraine patients and controls compared to DTI. In contrast, a higher number of differences was found with DTI parameters compared to AMURA in the comparisons between both migraine groups. Regarding the assessments reducing the sample size, the AMURA parameters showed a more stable behavior than DTI, showing a lower decrease for each reduced sample size or a higher number of regions with significant differences. However, most AMURA parameters showed lower stability in relation to higher coefficient of quartile variation values than the DTI descriptors, although two AMURA measures showed similar values to DTI. For the synthetic signals, there were AMURA measures with similar quantification to DTI, while other showed similar behavior. These findings suggest that AMURA presents favorable characteristics to identify differences of specific microstructural properties between clinical groups in regions with complex fiber architecture and lower dependency on the sample size or assessing technique than DTI.

Diffusion sampling schemes: A generalized methodology with nongeometric criteria.

PURPOSE: The aim of this paper is to show that geometrical criteria for designing multishell q $$ q $$ -space sampling procedures do not necessarily translate into reconstruction matrices with high figures of merit commonly used in the compressed sensing theory. In addition, we show that a well-known method for visiting k-space in radial three-dimensional acquisitions, namely, the Spiral Phyllotaxis, is a competitive initialization for the optimization of our nonconvex objective function. THEORY AND METHODS: We propose the gradient design method WISH (WeIghting SHells) which uses an objective function that accounts for weighted distances between gradients within M-tuples of consecutive shells, with M $$ M $$ ranging between 1 and the maximum number of shells S $$ S $$ . All the M $$ M $$ -tuples share the same weight ω M $$ {\omega}_M $$ . The objective function is optimized for a sample of these weights, using Spiral Phyllotaxis as initialization. State-of-the-art General Electrostatic Energy Minimization (GEEM) and Spherical Codes (SC) were used for comparison. For the three methods, reconstruction matrices of the attenuation signal using MAP-MRI were tested using figures of merit borrowed from the Compressed Sensing theory (namely, Restricted Isometry Property -RIP- and Coherence); we also tested the gradient design using a geometric criterion based on Voronoi cells. RESULTS: For RIP and Coherence, WISH got better results in at least one combination of weights, whilst the criterion based on Voronoi cells showed an unrelated pattern. CONCLUSION: The versatility provided by WISH is supported by better results. Optimization in the weight parameter space is likely to provide additional improvements. For a practical design with an intermediate number of gradients, our results recommend to carry out the methodology here used to determine the appropriate gradient table.

HYDI-DSI revisited: Constrained non-parametric EAP imaging without q-space re-gridding.

Hybrid Diffusion Imaging (HYDI) was one of the first attempts to use multi-shell samplings of the q-space to infer diffusion properties beyond Diffusion Tensor Imaging (DTI) or High Angular Resolution Diffusion Imaging (HARDI). HYDI was intended as a flexible protocol embedding both DTI (for lower b-values) and HARDI (for higher b-values) processing, as well as Diffusion Spectrum Imaging (DSI) when the entire data set was exploited. In the latter case, the spherical sampling of the q-space is re-gridded by interpolation to a Cartesian lattice whose extent covers the range of acquired b-values, hence being acquisition-dependent. The Discrete Fourier Transform (DFT) is afterwards used to compute the corresponding Cartesian sampling of the Ensemble Average Propagator (EAP) in an entirely non-parametric way. From this lattice, diffusion markers such as the Return To Origin Probability (RTOP) or the Mean Squared Displacement (MSD) can be numerically estimated. We aim at re-formulating this scheme by means of a Fourier Transform encoding matrix that eliminates the need for q-space re-gridding at the same time it preserves the non-parametric nature of HYDI-DSI. The encoding matrix is adaptively designed at each voxel according to the underlying DTI approximation, so that an optimal sampling of the EAP can be pursued without being conditioned by the particular acquisition protocol. The estimation of the EAP is afterwards carried out as a regularized Quadratic Programming (QP) problem, which allows to impose positivity constraints that cannot be trivially embedded within the conventional HYDI-DSI. We demonstrate that the definition of the encoding matrix in the adaptive space allows to analytically (as opposed to numerically) compute several popular descriptors of diffusion with the unique source of error being the cropping of high frequency harmonics in the Fourier analysis of the attenuation signal. They include not only RTOP and MSD, but also Return to Axis/Plane Probabilities (RTAP/RTPP), which are defined in terms of specific spatial directions and are not available with the former HYDI-DSI. We report extensive experiments that suggest the benefits of our proposal in terms of accuracy, robustness and computational efficiency, especially when only standard, non-dedicated q-space samplings are available.

Efficient estimation of propagator anisotropy and non-Gaussianity in multishell diffusion MRI with micro-structure adaptive convolution kernels and dual Fourier integral transforms.

PURPOSE: We seek to reformulate the so-called Propagator Anisotropy (PA) and Non-Gaussianity (NG), originally conceived for the Mean Apparent Propagator diffusion MRI (MAP-MRI), to the Micro-Structure adaptive convolution kernels and dual Fourier Integral Transforms (MiSFIT). These measures describe relevant normalized features of the Ensemble Average Propagator (EAP). THEORY AND METHODS: First, the indices, which are defined as the EAP's dissimilarity from an isotropic (PA) or a Gaussian (NG) one, are analytically reformulated within the MiSFIT framework. Then a comparison between the resulting maps is drawn by means of a visual analysis, a quantitative assessment via numerical simulations, a test-retest study across the MICRA dataset (6 subjects scanned five times) and, finally, a computational time evaluation. RESULTS: Findings illustrate the visual similarity between the indices computed with either technique. Evaluation against synthetic ground truth data, however, demonstrates MiSFIT's improved accuracy. In addition, the test-retest study reveals MiSFIT's higher degree of reliability in most of white matter regions. Finally, the computational time evaluation shows MiSFIT's time reduction up to two orders of magnitude. CONCLUSIONS: Despite being a direct development on the MAP-MRI representation, the PA and the NG can be reliably and efficiently computed within MiSFIT's framework. This, together with the previous findings in the original MiSFIT's article, could mean the difference that definitely qualifies diffusion MRI to be incorporated into regular clinical settings.

Influencia de la densidad de energía de ondas de choque focalizadas en el tratamiento de la fascitis plantar / Influence of energy density on the effectiveness of the treatment of plantar fasciitis with focused extracorporeal shock waves

Resumen La fascitis plantar (FP) es una patología frecuente e invalidante que puede tratarse con ondas de choque focalizadas. El objetivo principal del estudio fue valorar la eficacia del tratamiento con ondas de choque focalizadas en la FP según la densidad de energía utilizada. Se incluyeron 82 pacientes con diagnóstico clínico de FP que fueron asignados mediante muestreo aleatorio simple a dos grupos de tratamiento: densidad de energía media- alta (0,59mJ/mm2) y densidad de energía media-baja (0,27mJ/mm2). Se evaluaron el dolor y la funcionalidad, mediante las escalas EVA (Escala Visual Analógica) y AOFAS (American Orthopedic Foot and Ankle Society Ankle-Hindfoot Scale) respectivamente, al inicio del estudio (consulta base), y al primer, tercer y sexto mes tras el tratamiento. Por último, se evaluó el grado de satisfacción de los pacientes mediante la escala de Roles y Maudsley. Se compararon los resultados de las escalas en las revisiones posteriores al tratamiento, obteniéndose significación estadística para las variables principales del estudio (dolor y funcionalidad) en cada grupo de intervención. Aunque los niveles de dolor y la funcionalidad mejoraron en ambos grupos de estudio, se obtuvo una respuesta analgésica y funcional mayor y más precoz en el grupo tratado con densidad de energía media-alta.

Abstract Plantar fasciitis (FP) is a frequent and disabling condition that can be treated with focused extracorporeal shock waves. The main objective of this study was to assess the effectiveness of focused extracorporeal shockwave treatment in FP according to the energy density used. Eighty-two patients with a clinical diagnosis of FP were included and assigned, by simple random sampling, to two treatment groups: medium-high energy density (0.59mJ/mm2) and low-medium energy density (0.27mJ/mm2). Pain and functionality were assessed using the VAS (Visual Analogical Scale) and AOFAS (American Orthopedic Foot and Ankle Society Ankle-Hindfoot Scale) scales, respectively, at the start of the study (baseline consultation), and at the first, third and sixth month post-treatment. Finally, the degree of patient satisfaction was evaluated using the Roles and Maudsley score. The results of the scales in the post-treatment reviews were compared, and statistical significance was obtained for the main study variables (pain and functionality) in each intervention group. Although pain levels and functionality improved in both study groups after treatment, a greater and earlier analgesic and functional response was obtained for the medium-high energy density group.

Anisotropy measure from three diffusion-encoding gradient directions.

We propose a method that can provide information about the anisotropy and orientation of diffusion in the brain from only 3 orthogonal gradient directions without imposing additional assumptions. The method is based on the Diffusion Anisotropy (DiA) that measures the distance from a diffusion signal to its isotropic equivalent. The original formulation based on a Spherical Harmonics basis allows to go down to only 3 orthogonal directions in order to estimate the measure. In addition, an alternative simplification and a color-coding representation are also proposed. Acquisitions from a publicly available database are used to test the viability of the proposal. The DiA succeeded in providing anisotropy information from the white matter using only 3 diffusion-encoding directions. The price to pay for such reduced acquisition is an increment in the variability of the data and a subestimation of the metric on those tracts not aligned with the acquired directions. Nevertheless, the calculation of anisotropy information from DMRI is feasible using fewer than 6 gradient directions by using DiA. The method is totally compatible with existing acquisition protocols, and it may provide complementary information about orientation in fast diffusion acquisitions.

Moment-based representation of the diffusion inside the brain from reduced DMRI acquisitions: Generalized AMURA.

AMURA (Apparent Measures Using Reduced Acquisitions) was originally proposed as a method to infer micro-structural information from single-shell acquisitions in diffusion MRI. It reduces the number of samples needed and the computational complexity of the estimation of diffusion properties of tissues by assuming the diffusion anisotropy is roughly independent on the b-value. This simplification allows the computation of simplified expressions and makes it compatible with standard acquisition protocols commonly used even in clinical practice. The present work proposes an extension of AMURA that allows the calculation of general moments of the diffusion signals that can be applied to describe the diffusion process with higher accuracy. We provide simplified expressions to analytically compute a set of scalar indices as moments of arbitrary orders over either the whole 3-D space, particular directions, or particular planes. The existing metrics previously proposed for AMURA (RTOP, RTPP and RTAP) are now special cases of this generalization. An extensive set of experiments is performed on public data and a clinical clase acquired with a standard type acquisition. The new metrics provide additional information about the diffusion processes inside the brain.

Accurate free-water estimation in white matter from fast diffusion MRI acquisitions using the spherical means technique.

PURPOSE: To accurately estimate the partial volume fraction of free water in the white matter from diffusion MRI acquisitions not demanding strong sensitizing gradients and/or large collections of different b-values. Data sets considered comprise ∼ 32-64 gradients near b=1000s/mm2 plus ∼ 6 gradients near b=500s/mm2 . THEORY AND METHODS: The spherical means of each diffusion MRI set with the same b-value are computed. These means are related to the inherent diffusion parameters within the voxel (free- and cellular-water fractions; cellular-water diffusivity), which are solved by constrained nonlinear least squares regression. RESULTS: The proposed method outperforms those based on mixtures of two Gaussians for the kind of data sets considered. W.r.t. the accuracy, the former does not introduce significant biases in the scenarios of interest, while the latter can reach a bias of 5%-7% if fiber crossings are present. W.r.t. the precision, a variance near 10% , compared to 15%, can be attained for usual configurations. CONCLUSION: It is possible to compute reliable estimates of the free-water fraction inside the white matter by complementing typical DTI acquisitions with few gradients at a lowb-value. It can be done voxel-by-voxel, without imposing spatial regularity constraints.

Apparent propagator anisotropy from single-shell diffusion MRI acquisitions.

PURPOSE: The apparent propagator anisotropy (APA) is a new diffusion MRI metric that, while drawing on the benefits of the ensemble averaged propagator anisotropy (PA) compared to the fractional anisotropy (FA), can be estimated from single-shell data. THEORY AND METHODS: Computation of the full PA requires acquisition of large datasets with many diffusion directions and different b-values, and results in extremely long processing times. This has hindered adoption of the PA by the community, despite evidence that it provides meaningful information beyond the FA. Calculation of the complete propagator can be avoided under the hypothesis that a similar sensitivity/specificity may be achieved from apparent measurements at a given shell. Assuming that diffusion anisotropy (DiA) is nondependent on the b-value, a closed-form expression using information from one single shell (ie, b-value) is reported. RESULTS: Publicly available databases with healthy and diseased subjects are used to compare the APA against other anisotropy measures. The structural information provided by the APA correlates with that provided by the PA for healthy subjects, while it also reveals statistically relevant differences in white matter regions for two pathologies, with a higher reliability than the FA. Additionally, APA has a computational complexity similar to the FA, with processing-times several orders of magnitude below the PA. CONCLUSIONS: The APA can extract more relevant white matter information than the FA, without any additional demands on data acquisition. This makes APA an attractive option for adoption into existing diffusion MRI analysis pipelines.

Efficient and accurate EAP imaging from multi-shell dMRI with micro-structure adaptive convolution kernels and dual Fourier Integral Transforms (MiSFIT).

A number of computational techniques have been lately devised to image the Ensemble Average Propagator (EAP) within the white matter of the brain, propelled by the deployment of multi-shell acquisition protocols and databases: approaches like Mean Apparent Propagator Imaging (MAP-MRI) and its Laplacian-regularized version (MAPL) aim at describing the low frequency spectrum of the EAP (limited by the maximum b-value acquired) and afterwards computing scalar indices that embed useful descriptions of the white matter, e. g. the Return-to-Origin, Plane, or Axis Probabilities (RTOP, RTPP, RTAP). These methods resort to a non-parametric, bandwidth limited representation of the EAP that implies fitting a set of 3-D basis functions in a large-scale optimization problem. We propose a semi-parametric approach inspired by signal theory: the EAP is approximated as the spherical convolution of a Micro-Structure adaptive Gaussian kernel with a non-parametric orientation histogram, which aims at representing the low-frequency response of an ensemble of coherent sets of fiber bundles at the white matter. This way, the optimization involves just the 2 to 3 parameters that describe the kernel, making our approach far more efficient than the related state of the art. We devise dual Fourier domains Integral Transforms to analytically compute RTxP-like scalar indices as moments of arbitrary orders over either the whole 3-D space, particular directions, or particular planes. The so-called MiSFIT is both time efficient (a typical multi-shell data set can be processed in roughly one minute) and accurate: it provides estimates of widely validated indices like RTOP, RTPP, and RTAP comparable to MAPL for a wide variety of white matter configurations.

Micro-structure diffusion scalar measures from reduced MRI acquisitions.

In diffusion MRI, the Ensemble Average diffusion Propagator (EAP) provides relevant micro-structural information and meaningful descriptive maps of the white matter previously obscured by traditional techniques like Diffusion Tensor Imaging (DTI). The direct estimation of the EAP, however, requires a dense sampling of the Cartesian q-space involving a huge amount of samples (diffusion gradients) for proper reconstruction. A collection of more efficient techniques have been proposed in the last decade based on parametric representations of the EAP, but they still imply acquiring a large number of diffusion gradients with different b-values (shells). Paradoxically, this has come together with an effort to find scalar measures gathering all the q-space micro-structural information probed in one single index or set of indices. Among them, the return-to-origin (RTOP), return-to-plane (RTPP), and return-to-axis (RTAP) probabilities have rapidly gained popularity. In this work, we propose the so-called "Apparent Measures Using Reduced Acquisitions" (AMURA) aimed at computing scalar indices that can mimic the sensitivity of state of the art EAP-based measures to micro-structural changes. AMURA drastically reduces both the number of samples needed and the computational complexity of the estimation of diffusion properties by assuming the diffusion anisotropy is roughly independent from the radial direction. This simplification allows us to compute closed-form expressions from single-shell information, so that AMURA remains compatible with standard acquisition protocols commonly used even in clinical practice. Additionally, the analytical form of AMURA-based measures, as opposed to the iterative, non-linear reconstruction ubiquitous to full EAP techniques, turns the newly introduced apparent RTOP, RTPP, and RTAP both robust and efficient to compute.

Scalar diffusion-MRI measures invariant to acquisition parameters: A first step towards imaging biomarkers.

An imaging biomarker is a biologic feature in an image that is relevant to a patient's diagnosis or prognosis. In order to qualify as a biomarker, a measure must be robust and reproducible. However, the usual scalar measures derived from diffusion tensor imaging are known to be highly dependent on the variation of the acquisition parameters, which prevents their possible use as biomarkers. In this work, we propose a new set of quantitative measures based on diffusion magnetic resonance imaging from single-shell acquisitions that are designed to be robust to the variations of several acquisition parameters (number of gradient directions, b-value and SNR) while keeping a high discrimination power on differences in the diffusion characteristics of the tissue. These new scalar measures are analytically obtained from a generic diffusion function that does not require the calculation of a diffusion tensor. This way, on one hand, we avoid the use of a specific diffusion model and, on the other hand, we make easier the statistical characterization of the measures. Accordingly, the analysis of the measures bias is carried out and it is used to minimize their dependency with respect to the acquisition noise for different SNRs. The robustness and discrimination power of the measures are tested for different number of gradients, b-values and SNRs using a realistic phantom and three real datasets: (1) 13 control subjects and different acquisition parameters; (2) a public data set from a single subject acquired using multiple shells and (3) 32 schizophrenia patients and 32 age and sex-matched healthy controls with a varying number of gradient directions. The proposed quantitative measures exhibit low variability to the changes of the acquisition parameters, while at the same time they preserve a discrimination power that is able to detect significant changes in the anisotropy of the diffusion.

Impact of MR Acquisition Parameters on DTI Scalar Indexes: A Tractography Based Approach.

Acquisition parameters play a crucial role in Diffusion Tensor Imaging (DTI), as they have a major impact on the values of scalar measures such as Fractional Anisotropy (FA) or Mean Diffusivity (MD) that are usually the focus of clinical studies based on white matter analysis. This paper presents an analysis on the impact of the variation of several acquisition parameters on these scalar measures with a novel double focus. First, a tractography-based approach is employed, motivated by the significant number of clinical studies that are carried out using this technique. Second, the consequences of simultaneous changes in multiple parameters are analyzed: number of gradient directions, b-value and voxel resolution. Results indicate that the FA is most affected by changes in the number of gradients and voxel resolution, while MD is specially influenced by variations in the b-value. Even if the choice of a tractography algorithm has an effect on the numerical values of the final scalar measures, the evolution of these measures when acquisition parameters are modified is parallel.

Efficient and Robust Image Restoration Using Multiple-Feature L2-Relaxed Sparse Analysis Priors.

We propose a novel formulation for relaxed analysis-based sparsity in multiple dictionaries as a general type of prior for images, and apply it for Bayesian estimation in image restoration problems. Our formulation of a ℓ2-relaxed ℓ0 pseudo-norm prior allows for an especially simple maximum a posteriori estimation iterative marginal optimization algorithm, whose convergence we prove. We achieve a significant speedup over the direct (static) solution by using dynamically evolving parameters through the estimation loop. As an added heuristic twist, we fix in advance the number of iterations, and then empirically optimize the involved parameters according to two performance benchmarks. The resulting constrained dynamic method is not just fast and effective, it is also highly robust and flexible. First, it is able to provide an outstanding tradeoff between computational load and performance, in visual and objective, mean square error and structural similarity terms, for a large variety of degradation tests, using the same set of parameter values for all tests. Second, the performance benchmark can be easily adapted to specific types of degradation, image classes, and even performance criteria. Third, it allows for using simultaneously several dictionaries with complementary features. This unique combination makes ours a highly practical deconvolution method.

The DTI Challenge: Toward Standardized Evaluation of Diffusion Tensor Imaging Tractography for Neurosurgery.

BACKGROUND AND PURPOSE: Diffusion tensor imaging (DTI) tractography reconstruction of white matter pathways can help guide brain tumor resection. However, DTI tracts are complex mathematical objects and the validity of tractography-derived information in clinical settings has yet to be fully established. To address this issue, we initiated the DTI Challenge, an international working group of clinicians and scientists whose goal was to provide standardized evaluation of tractography methods for neurosurgery. The purpose of this empirical study was to evaluate different tractography techniques in the first DTI Challenge workshop. METHODS: Eight international teams from leading institutions reconstructed the pyramidal tract in four neurosurgical cases presenting with a glioma near the motor cortex. Tractography methods included deterministic, probabilistic, filtered, and global approaches. Standardized evaluation of the tracts consisted in the qualitative review of the pyramidal pathways by a panel of neurosurgeons and DTI experts and the quantitative evaluation of the degree of agreement among methods. RESULTS: The evaluation of tractography reconstructions showed a great interalgorithm variability. Although most methods found projections of the pyramidal tract from the medial portion of the motor strip, only a few algorithms could trace the lateral projections from the hand, face, and tongue area. In addition, the structure of disagreement among methods was similar across hemispheres despite the anatomical distortions caused by pathological tissues. CONCLUSIONS: The DTI Challenge provides a benchmark for the standardized evaluation of tractography methods on neurosurgical data. This study suggests that there are still limitations to the clinical use of tractography for neurosurgical decision making.

Improving GRAPPA reconstruction by frequency discrimination in the ACS lines.

PURPOSE: GRAPPA is a well-known parallel imaging method that recovers the MR magnitude image from aliasing by using a weighted interpolation of the data in k-space. To estimate the optimal reconstruction weights, GRAPPA uses a band along the center of the k-space where the signal is sampled at the Nyquist rate, the so-called autocalibrated (ACS) lines. However, while the subsampled lines usually belong to the medium- to high-frequency areas of the spectrum, the ACS lines include the low-frequency areas around the DC component. The use for estimation and reconstruction of areas of the k-space with very different features may negatively affect the final reconstruction quality. We propose a simple, yet powerful method to eliminate reconstruction artifacts, based on the discrimination of the low-frequency spectrum. METHODS: The proposal to improve the estimation of the weights lays on a proper selection of the coefficients within the ACS lines, which advises discarding those points around the DC component. A simple approach is the elimination of a square window in the center of the k-space, although more developed approaches can be used. RESULTS: The method is tested using real multiple-coil MRI acquisitions. We empirically show this approach achieves great enhancement rates, while keeping the same complexity of the original GRAPPA and reducing the g-factor. The reconstruction is even more accurate when combined with other reconstruction methods. Improvement rates of 35% are achieved for 32 ACS and acceleration rate of 3. CONCLUSIONS: The method proposed highly improves the accuracy of the GRAPPA coefficients and therefore the final image reconstruction. The method is fully compatible with the original GRAPPA formulation and with other optimization methods proposed in literature, and it can be easily implemented into the commercial scanning software.

Noise estimation in parallel MRI: GRAPPA and SENSE.

Parallel imaging methods allow to increase the acquisition rate via subsampled acquisitions of the k-space. SENSE and GRAPPA are the most popular reconstruction methods proposed in order to suppress the artifacts created by this subsampling. The reconstruction process carried out by both methods yields to a variance of noise value which is dependent on the position within the final image. Hence, the traditional noise estimation methods - based on a single noise level for the whole image - fail. In this paper we propose a novel methodology to estimate the spatial dependent pattern of the variance of noise in SENSE and GRAPPA reconstructed images. In both cases, some additional information must be known beforehand: the sensitivity maps of each receiver coil in the SENSE case and the reconstruction coefficients for GRAPPA.

Merging squared-magnitude approaches to DWI denoising: An adaptive Wiener filter tuned to the anatomical contents of the image.

We present a new method for denoising of Diffusion Weighted Images (DWI) that shares several desirable features of state-of-the-art proposals: 1) it works with the squared-magnitude signal, allowing for a closed-form formulation as a Linear Minimum Mean Squared Error (LMMSE) estimator, a.k.a. Wiener filter; 2) it jointly accounts for the DWI channels altogether, being able to unveil anatomical structures that remain hidden in each separated channel; 3) it uses a Non-Local Means (NLM)-like scheme to discriminate voxels corresponding to different fiber bundles, being able to enhance the anatomical structures of the DWI. We report extensive experiments evidencing the new approach outperforms several related methods for all the range of input signal-to-noise ratios (SNR). An open-source C++ implementation of the algorithm is also provided for the sake of reproducibility.