MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation

Denoising of diffusion MRI using random matrix theory

Tractography based on non-invasive diffusion imaging is central to the study of human brain connectivity. To date, the approach has not been systematically validated in ground truth studies. Based on a simulated human brain data set with ground truth tracts, we organized an open international tractography challenge, which resulted in 96 distinct submissions from 20 research groups. Here, we report the encouraging finding that most state-of-the-art algorithms produce tractograms containing 90% of the ground truth bundles (to at least some extent). However, the same tractograms contain many more invalid than valid bundles, and half of these invalid bundles occur systematically across research groups. Taken together, our results demonstrate and confirm fundamental ambiguities inherent in tract reconstruction based on orientation information alone, which need to be considered when interpreting tractography and connectivity results. Our approach provides a novel framework for estimating reliability of tractography and encourages innovation to address its current limitations.

/pdf/the-challenge-of-mapping-the-human-connectome-based-on-xi6eaf30xo.pdf

The challenge of mapping the human connectome based on diffusion tractography

Purpose To develop a fast and stable method for correcting the gibbs-ringing artifact. Methods Gibbs-ringing is a well-known artifact which manifests itself as spurious oscillations in the vicinity of sharp image gradients at tissue boundaries. The origin can be seen in the truncation of k-space during MRI data-acquisition. Correction techniques like Gegenbauer reconstruction or extrapolation methods aim at recovering these missing data. Here, we present a simple and robust method which exploits a different view on the Gibbs-phenomenon: The truncation in k-space can be interpreted as a convolution of the underlying image with a sinc-function. As the image is reconstructed on a discretized grid, the severity of the ringing artifacts depends on how this grid is located with respect to the edge and the oscillation pattern of the function. We propose to reinterpolate the image based on local, subvoxel-shifts to sample the ringing pattern at the zero-crossings of the oscillating sinc-function. Results With the proposed method, the artifact can simply, effectively, and robustly be removed with a minimal amount of image smoothing. Conclusions The robustness of the method suggests it as a suitable candidate for an implementation in the standard image processing pipeline in clinical routine. Magn Reson Med 76:1574-1581, 2016. © 2015 International Society for Magnetic Resonance in Medicine.

/pdf/gibbs-ringing-artifact-removal-based-on-local-subvoxel-l9ymfwwgbm.pdf

Gibbs-ringing artifact removal based on local subvoxel-shifts.

Judgement, as one of the core tenets of medicine, relies upon the integration of multilayered data with nuanced decision making. Cancer offers a unique context for medical decisions given not only its variegated forms with evolution of disease but also the need to take into account the individual condition of patients, their ability to receive treatment, and their responses to treatment. Challenges remain in the accurate detection, characterization, and monitoring of cancers despite improved technologies. Radiographic assessment of disease most commonly relies upon visual evaluations, the interpretations of which may be augmented by advanced computational analyses. In particular, artificial intelligence (AI) promises to make great strides in the qualitative interpretation of cancer imaging by expert clinicians, including volumetric delineation of tumors over time, extrapolation of the tumor genotype and biological course from its radiographic phenotype, prediction of clinical outcome, and assessment of the impact of disease and treatment on adjacent organs. AI may automate processes in the initial interpretation of images and shift the clinical workflow of radiographic detection, management decisions on whether or not to administer an intervention, and subsequent observation to a yet to be envisioned paradigm. Here, the authors review the current state of AI as applied to medical imaging of cancer and describe advances in 4 tumor types (lung, brain, breast, and prostate) to illustrate how common clinical problems are being addressed. Although most studies evaluating AI applications in oncology to date have not been vigorously validated for reproducibility and generalizability, the results do highlight increasingly concerted efforts in pushing AI technology to clinical use and to impact future directions in cancer care.

https://acsjournals.onlinelibrary.wiley.com/doi/epdf/10.3322/caac.21552

Artificial intelligence in cancer imaging: Clinical challenges and applications.

PURPOSE To estimate the spatially varying noise map using a redundant series of magnitude MR images. METHODS We exploit redundancy in non-Gaussian distributed multidirectional diffusion MRI data by identifying its noise-only principal components, based on the theory of noisy covariance matrices. The bulk of principal component analysis eigenvalues, arising due to noise, is described by the universal Marchenko-Pastur distribution, parameterized by the noise level. This allows us to estimate noise level in a local neighborhood based on the singular value decomposition of a matrix combining neighborhood voxels and diffusion directions. RESULTS We present a model-independent local noise mapping method capable of estimating the noise level down to about 1% error. In contrast to current state-of-the-art techniques, the resultant noise maps do not show artifactual anatomical features that often reflect physiological noise, the presence of sharp edges, or a lack of adequate a priori knowledge of the expected form of MR signal. CONCLUSIONS Simulations and experiments show that typical diffusion MRI data exhibit sufficient redundancy that enables accurate, precise, and robust estimation of the local noise level by interpreting the principal component analysis eigenspectrum in terms of the Marchenko-Pastur distribution. Magn Reson Med 76:1582-1593, 2016. © 2015 International Society for Magnetic Resonance in Medicine.

/pdf/diffusion-mri-noise-mapping-using-random-matrix-theory-u2oncjex5k.pdf

Diffusion MRI noise mapping using random matrix theory.

Weighted linear least squares estimation of diffusion MRI parameters: Strengths, limitations, and pitfalls

Mean kurtosis and linearly related kurtosis tensor parameters are potential new noninvasive biomarkers in the grading of gliomas.

https://pubs.rsna.org/doi/pdf/10.1148/radiol.12110927

Gliomas: Diffusion Kurtosis MR Imaging in Grading

The ultimate promise of diffusion MRI (dMRI) models is specificity to neuronal microstructure, which may lead to distinct clinical biomarkers using noninvasive imaging. While multi-compartment models are a common approach to interpret water diffusion in the brain in vivo, the estimation of their parameters from the dMRI signal remains an unresolved problem. Practically, even when q space is highly oversampled, nonlinear fit outputs suffer from heavy bias and poor precision. So far, this has been alleviated by fixing some of the model parameters to a priori values, for improved precision at the expense of accuracy. Here we use a representative two-compartment model to show that fitting fails to determine the five model parameters from over 60 measurement points. For the first time, we identify the reasons for this poor performance. The first reason is the existence of two local minima in the parameter space for the objective function of the fitting procedure. These minima correspond to qualitatively different sets of parameters, yet they both lie within biophysically plausible ranges. We show that, at realistic signal-to-noise ratio values, choosing between the two minima based on the associated objective function values is essentially impossible. Second, there is an ensemble of very low objective function values around each of these minima in the form of a pipe. The existence of such a direction in parameter space, along which the objective function profile is very flat, explains the bias and large uncertainty in parameter estimation, and the spurious parameter correlations: in the presence of noise, the minimum can be randomly displaced by a very large amount along each pipe. Our results suggest that the biophysical interpretation of dMRI model parameters crucially depends on establishing which of the minima is closer to the biophysical reality and the size of the uncertainty associated with each parameter.

/pdf/degeneracy-in-model-parameter-estimation-for-multi-24p30l0fqx.pdf

Jelle Veraart

Papers

Denoising of diffusion MRI using random matrix theory

Diffusion MRI noise mapping using random matrix theory.

Weighted linear least squares estimation of diffusion MRI parameters: Strengths, limitations, and pitfalls

Gliomas: Diffusion Kurtosis MR Imaging in Grading

Degeneracy in model parameter estimation for multi‐compartmental diffusion in neuronal tissue