A joint compressed-sensing and super-resolution approach for very high-resolution diffusion imaging

POSE: POSition Encoding for accelerated quantitative MRI

Quantitative MRI utilizes multiple acquisitions with varying sequence parameters to sufficiently characterize a biophysical model of interest, resulting in undesirable scan times. Here we propose, validate and demonstrate a new general strategy for accelerating MRI using subvoxel shifting as a source of encoding called POSition Encoding (POSE). The POSE framework applies unique subvoxel shifts along the acquisition parameter dimension, thereby creating an extra source of encoding. Combining with a biophysical signal model of interest, accelerated and enhanced resolution maps of biophysical parameters are obtained. This has been validated and demonstrated through numerical Bloch equation simulations, phantom experiments and in vivo experiments using the variable flip angle signal model in 3D acquisitions as an application example. Monte Carlo simulations were performed using in vivo data to investigate our method's noise performance. POSE quantification results from numerical Bloch equation simulations of both a numerical phantom and realistic digital brain phantom concur well with the reference method, validating our method both theoretically and for realistic situations. NIST phantom experiment results show excellent overall agreement with the reference method, confirming our method's applicability for a wide range of T1 values. In vivo results not only exhibit good agreement with the reference method, but also show g-factors that significantly outperforms conventional parallel imaging methods with identical acceleration. Furthermore, our results show that POSE can be combined with parallel imaging to further accelerate while maintaining superior noise performance over parallel imaging that uses lower acceleration factors.

MHAN: Multi-Stage Hybrid Attention Network for MRI reconstruction and super-resolution

High-quality magnetic resonance imaging (MRI) affords clear body tissue structure for reliable diagnosing. However, there is a principal problem of the trade-off between acquisition speed and image quality. Image reconstruction and super-resolution are crucial techniques to solve these problems. In the main field of MR image restoration, most researchers mainly focus on only one of these aspects, namely reconstruction or super-resolution. In this paper, we propose an efficient model called Multi-Stage Hybrid Attention Network (MHAN) that performs the multi-task of recovering high-resolution (HR) MR images from low-resolution (LR) under-sampled measurements. Our model is highlighted by three major modules: (i) an Amplified Spatial Attention Block (ASAB) capable of enhancing the differences in spatial information, (ii) a Self-Attention Block with a Data-Consistency Layer (DC-SAB), which can improve the accuracy of the extracted feature information, (iii) an Adaptive Local Residual Attention Block (ALRAB) that focuses on both spatial and channel information. MHAN employs an encoder-decoder architecture to deeply extract contextual information and a pipeline to provide spatial accuracy. Compared with the recent multi-task model T2Net, our MHAN improves by 2.759 dB in PSNR and 0.026 in SSIM with scaling factor ×2 and acceleration factor 4× on T2 modality.

Rapid Diffusion Magnetic Resonance Imaging Using Slice-Interleaved Encoding

In this paper, we present a robust reconstruction scheme for diffusion MRI (dMRI) data acquired using slice-interleaved diffusion encoding (SIDE). When combined with SIDE undersampling and simultaneous multi-slice (SMS) imaging, our reconstruction strategy is capable of significantly reducing the amount of data that needs to be acquired, enabling high-speed diffusion imaging for pediatric, elderly, and claustrophobic individuals. In contrast to the conventional approach of acquiring a full diffusion-weighted (DW) volume per diffusion wavevector, SIDE acquires in each repetition time (TR) a volume that consists of interleaved slice groups, each group corresponding to a different diffusion wavevector. This strategy allows SIDE to rapidly acquire data covering a large number of wavevectors within a short period of time. The proposed reconstruction method uses a diffusion spectrum model and multi-dimensional total variation to recover full DW images from DW volumes that are slice-undersampled due to unacquired SIDE volumes. We formulate an inverse problem that can be solved efficiently using the alternating direction method of multipliers (ADMM). Experiment results demonstrate that DW images can be reconstructed with high fidelity even when the acquisition is accelerated by 25 folds.

Through-Plane Super-Resolution With Autoencoders in Diffusion Magnetic Resonance Imaging of the Developing Human Brain

Fetal brain diffusion magnetic resonance images (MRI) are often acquired with a lower through-plane than in-plane resolution. This anisotropy is often overcome by classical upsampling methods such as linear or cubic interpolation. In this work, we employ an unsupervised learning algorithm using an autoencoder neural network for single-image through-plane super-resolution by leveraging a large amount of data. Our framework, which can also be used for slice outliers replacement, overperformed conventional interpolations quantitatively and qualitatively on pre-term newborns of the developing Human Connectome Project. The evaluation was performed on both the original diffusion-weighted signal and the estimated diffusion tensor maps. A byproduct of our autoencoder was its ability to act as a denoiser. The network was able to generalize fetal data with different levels of motions and we qualitatively showed its consistency, hence supporting the relevance of pre-term datasets to improve the processing of fetal brain images.

Automated cerebral cortex segmentation based solely on diffusion tensor imaging for investigating cortical anisotropy

To extract Diffusion Tensor Imaging (DTI) parameters from the human cortex, the inner and outer boundaries of the cortex are usually defined on 3D-T1-weighted images and then applied to the co-registered DTI. However, this analysis requires the acquisition of an additional high-resolution structural image that may not be practical in various imaging studies. Here an automatic cortical boundary segmentation method was developed to work directly only on the native DTI images by using fractional anisotropy (FA) maps and mean diffusion weighted images (DWI), the latter with acceptable gray-white matter image contrast. This new method was compared to the conventional cortical segmentations generated from high-resolution T1 structural images in 5 participants. In addition, the proposed method was applied to 15 healthy young adults (10 cross-sectional, 5 test-retest) to measure FA, MD, and radiality of the primary eigenvector across the cortex on whole-brain 1.5 mm isotropic images acquired in 3.5 min at 3T. The proposed method generated reasonable segmentations of the cortical boundaries for all individuals and large proportions of the proposed method segmentations (more than 85%) were within ±1 mm from those generated with the conventional approach on higher resolution T1 structural images. Both FA (~0.15) and MD (~0.77 × 10^-3 mm²/s) extracted halfway between the cortical boundaries were relatively stable across the cortex, although focal regions such as the posterior bank of the central sulcus, anterior insula, and medial temporal lobe showed higher FA. The primary eigenvectors were primarily oriented radially to the middle cortical surface, but there were tangential orientations in the sulcal fundi as well as in the posterior bank of the central sulcus. The proposed method demonstrates the feasibility and accuracy of cortical analysis in native DTI space while avoiding the acquisition of other imaging contrasts like 3D T1-weighted scans.

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.

High-fidelity, accelerated whole-brain submillimeter in vivo diffusion MRI using gSlider-spherical ridgelets (gSlider-SR)

PURPOSE

To develop an accelerated, robust, and accurate diffusion MRI acquisition and reconstruction technique for submillimeter whole human brain in vivo scan on a clinical scanner.

METHODS

We extend the ultra-high resolution diffusion MRI acquisition technique, gSlider, by allowing undersampling in q-space and radiofrequency (RF)-encoding space, thereby dramatically reducing the total acquisition time of conventional gSlider. The novel method, termed gSlider-SR, compensates for the lack of acquired information by exploiting redundancy in the dMRI data using a basis of spherical ridgelets (SR), while simultaneously enhancing the signal-to-noise ratio. Using Monte Carlo simulation with realistic noise levels and several acquisitions of in vivo human brain dMRI data (acquired on a Siemens Prisma 3T scanner), we demonstrate the efficacy of our method using several quantitative metrics.

RESULTS

For high-resolution dMRI data with realistic noise levels (synthetically added), we show that gSlider-SR can reconstruct high-quality dMRI data at different acceleration factors preserving both signal and angular information. With in vivo data, we demonstrate that gSlider-SR can accurately reconstruct 860 μm diffusion MRI data (64 diffusion directions at b = 2000 s / mm 2 ), at comparable quality as that obtained with conventional gSlider with four averages, thereby providing an eight-fold reduction in scan time (from 1 hour 20 to 10 minutes).

CONCLUSIONS

gSlider-SR enables whole-brain high angular resolution dMRI at a submillimeter spatial resolution with a dramatically reduced acquisition time, making it feasible to use the proposed scheme on existing clinical scanners.

Challenges for biophysical modeling of microstructure

The biophysical modeling efforts in diffusion MRI have grown considerably over the past 25 years. In this review, we dwell on the various challenges along the journey of bringing a biophysical model from initial design to clinical implementation, identifying both hurdles that have been already overcome and outstanding issues. First, we describe the critical initial task of selecting which features of tissue microstructure can be estimated using a model and which acquisition protocol needs to be implemented to make the estimation possible. The model performance should necessarily be tested in realistic numerical simulations and in experimental data - adapting the fitting strategy accordingly, and parameter estimates should be validated against complementary techniques, when/if available. Secondly, the model performance and validity should be explored in pathological conditions, and, if appropriate, dedicated models for pathology should be developed. We build on examples from tumors, ischemia and demyelinating diseases. We then discuss the challenges associated with clinical translation and added value. Finally, we single out four major unresolved challenges that are related to: the availability of a microstructural ground truth, the validation of model parameters which cannot be accessed with complementary techniques, the development of a generalized standard model for any brain region and pathology, and the seamless communication between different parties involved in the development and application of biophysical models of diffusion.

Fast submillimeter diffusion MRI using gSlider-SMS and SNR-enhancing joint reconstruction

PURPOSE

We evaluate a new approach for achieving diffusion MRI data with high spatial resolution, large volume coverage, and fast acquisition speed.

THEORY AND METHODS

A recent method called gSlider-SMS enables whole-brain submillimeter diffusion MRI with high signal-to-noise ratio (SNR) efficiency. However, despite the efficient acquisition, the resulting images can still suffer from low SNR due to the small size of the imaging voxels. This work proposes to mitigate the SNR problem by combining gSlider-SMS with a regularized SNR-enhancing reconstruction approach.

RESULTS

Illustrative results show that, from gSlider-SMS data acquired over a span of only 25 minutes on a 3T scanner, the proposed method is able to produce 71 MRI images (64 diffusion encoding orientations with b = 1500 s/ mm 2 , and 7 images without diffusion weighting) of the entire in vivo human brain with nominal 0.66 mm spatial resolution. Using data acquired from 75 minutes of acquisition as a gold standard reference, we demonstrate that the proposed SNR-enhancement procedure leads to substantial improvements in estimated diffusion parameters compared to conventional gSlider reconstruction. Results also demonstrate that the proposed method has advantages relative to denoising methods based on low-rank matrix modeling. A theoretical analysis of the trade-off between spatial resolution and SNR suggests that the proposed approach has high efficiency.

CONCLUSIONS

The combination of gSlider-SMS with advanced regularized reconstruction enables high-resolution quantitative diffusion MRI from a relatively fast acquisition.

Anatomical assessment of trigeminal nerve tractography using diffusion MRI: A comparison of acquisition b-values and single- and multi-fiber tracking strategies

BACKGROUND

The trigeminal nerve (TGN) is the largest cranial nerve and can be involved in multiple inflammatory, compressive, ischemic or other pathologies. Currently, imaging-based approaches to identify the TGN mostly rely on T2-weighted magnetic resonance imaging (MRI), which provides localization of the cisternal portion of the TGN where the contrast between nerve and cerebrospinal fluid (CSF) is high enough to allow differentiation. The course of the TGN within the brainstem as well as anterior to the cisternal portion, however, is more difficult to display on traditional imaging sequences. An advanced imaging technique, diffusion MRI (dMRI), enables tracking of the trajectory of TGN fibers and has the potential to visualize anatomical regions of the TGN not seen on T2-weighted imaging. This may allow a more comprehensive assessment of the nerve in the context of pathology. To date, most work in TGN tracking has used clinical dMRI acquisitions with a b-value of 1000 s/mm2 and conventional diffusion tensor MRI (DTI) tractography methods. Though higher b-value acquisitions and multi-tensor tractography methods are known to be beneficial for tracking brain white matter fiber tracts, there have been no studies conducted to evaluate the performance of these advanced approaches on nerve tracking of the TGN, in particular on tracking different anatomical regions of the TGN.

OBJECTIVE

We compare TGN tracking performance using dMRI data with different b-values, in combination with both single- and multi-tensor tractography methods. Our goal is to assess the advantages and limitations of these different strategies for identifying the anatomical regions of the TGN.

METHODS

We proposed seven anatomical rating criteria including true and false positive structures, and we performed an expert rating study of over 1000 TGN visualizations, as follows. We tracked the TGN using high-quality dMRI data from 100 healthy adult subjects from the Human Connectome Project (HCP). TGN tracking performance was compared across dMRI acquisitions with b = 1000 s/mm2, b = 2000 s/mm2 and b = 3000 s/mm2, using single-tensor (1T) and two-tensor (2T) unscented Kalman filter (UKF) tractography. This resulted in a total of six tracking strategies. The TGN was identified using an anatomical region-of-interest (ROI) selection approach. First, in a subset of the dataset we identified ROIs that provided good TGN tracking performance across all tracking strategies. Using these ROIs, the TGN was then tracked in all subjects using the six tracking strategies. An expert rater (GX) visually assessed and scored each TGN based on seven anatomical judgment criteria. These criteria included the presence of multiple expected anatomical segments of the TGN (true positive structures), specifically branch-like structures, cisternal portion, mesencephalic trigeminal tract, and spinal cord tract of the TGN. False positive criteria included the presence of any fibers entering the temporal lobe, the inferior cerebellar peduncle, or the middle cerebellar peduncle. Expert rating scores were analyzed to compare TGN tracking performance across the six tracking strategies. Intra- and inter-rater validation was performed to assess the reliability of the expert TGN rating result.

RESULTS

The TGN was selected using two anatomical ROIs (Meckel's Cave and cisternal portion of the TGN). The two-tensor tractography method had significantly better performance on identifying true positive structures, while generating more false positive streamlines in comparison to the single-tensor tractography method. TGN tracking performance was significantly different across the three b-values for almost all structures studied. Tracking performance was reported in terms of the percentage of subjects achieving each anatomical rating criterion. Tracking of the cisternal portion and branching structure of the TGN was generally successful, with the highest performance of over 98% using two-tensor tractography and b = 1000 or b = 2000. However, tracking the smaller mesencephalic and spinal cord tracts of the TGN was quite challenging (highest performance of 37.5% and 57.07%, using two-tensor tractography with b = 1000 and b = 2000, respectively). False positive connections to the temporal lobe (over 38% of subjects for all strategies) and cerebellar peduncles (100% of subjects for all strategies) were prevalent. High joint probability of agreement was obtained in the inter-rater (on average 83%) and intra-rater validation (on average 90%), showing a highly reliable expert rating result.

CONCLUSIONS

Overall, the results of the study suggest that researchers and clinicians may benefit from tailoring their acquisition and tracking methodology to the specific anatomical portion of the TGN that is of the greatest interest. For example, tracking of branching structures and TGN-T2 overlap can be best achieved with a two-tensor model and an acquisition using b = 1000 or b = 2000. In general, b = 1000 and b = 2000 acquisitions provided the best-rated tracking results. Further research is needed to improve both sensitivity and specificity of the depiction of the TGN anatomy using dMRI.

Globally Optimized Super-Resolution of Diffusion MRI Data via Fiber Continuity

In this paper, we introduce a technique for super-resolution reconstruction of diffusion MRI, harnessing fiber-continuity (FC) as a constraint in a global whole-brain optimization framework. FC is a biologically-motivated constraint that relates orientation information between neighboring voxels. We show that it can be used to effectively constrain the inverse problem of recovering high-resolution data from low-resolution data. Since voxels are inter-related by FC, we devise a global optimization framework that allows solutions pertaining to all voxels to be solved simultaneously. We demonstrate that the proposed super-resolution framework is effective for diffusion MRI data of a glioma patient, a healthy subject, and a macaque.

XQ-SR: Joint x-q space super-resolution with application to infant diffusion MRI

Diffusion MRI (DMRI) is a powerful tool for studying early brain development and disorders. However, the typically low spatio-angular resolution of DMRI diminishes structural details and limits quantitative analysis to simple diffusion models. This problem is aggravated for infant DMRI since (i) the infant brain is significantly smaller than that of an adult, demanding higher spatial resolution to capture subtle structures; and (ii) the typically limited scan time of unsedated infants poses significant challenges to DMRI acquisition with high spatio-angular resolution. Post-acquisition super-resolution (SR) is an important alternative for increasing the resolution of DMRI data without prolonging acquisition times. However, most existing methods focus on the SR of only either the spatial domain (x-space) or the diffusion wavevector domain (q-space). For more effective resolution enhancement, we propose a framework for joint SR in both spatial and wavevector domains. More specifically, we first establish the signal relationships in x-q space using a robust neighborhood matching technique. We then harness the signal relationships to regularize the ill-posed inverse problem associated with the recovery of high-resolution data from their low-resolution counterpart. Extensive experiments on synthetic, adult, and infant DMRI data demonstrate that our method is able to recover high-resolution DMRI data with remarkably improved quality.

Diffusion Acceleration with Gaussian process Estimated Reconstruction (DAGER)

PURPOSE

Image acceleration provides multiple benefits to diffusion MRI, with in-plane acceleration reducing distortion and slice-wise acceleration increasing the number of directions that can be acquired in a given scan time. However, as acceleration factors increase, the reconstruction problem becomes ill-conditioned, particularly when using both in-plane acceleration and simultaneous multislice imaging. In this work, we develop a novel reconstruction method for in vivo MRI acquisition that provides acceleration beyond what conventional techniques can achieve.

THEORY AND METHODS

We propose to constrain the reconstruction in the spatial (k) domain by incorporating information from the angular (q) domain. This approach exploits smoothness of the signal in q-space using Gaussian processes, as has previously been exploited in post-reconstruction analysis. We demonstrate in-plane undersampling exceeding the theoretical parallel imaging limits, and simultaneous multislice combined with in-plane undersampling at a total factor of 12. This reconstruction is cast within a Bayesian framework that incorporates estimation of smoothness hyper-parameters, with no need for manual tuning.

RESULTS

Simulations and in vivo results demonstrate superior performance of the proposed method compared with conventional parallel imaging methods. These improvements are achieved without loss of spatial or angular resolution and require only a minor modification to standard pulse sequences.

CONCLUSION

The proposed method provides improvements over existing methods for diffusion acceleration, particularly for high simultaneous multislice acceleration with in-plane undersampling.

Test-retest reproducibility of white matter parcellation using diffusion MRI tractography fiber clustering

There are two popular approaches for automated white matter parcellation using diffusion MRI tractography, including fiber clustering strategies that group white matter fibers according to their geometric trajectories and cortical-parcellation-based strategies that focus on the structural connectivity among different brain regions of interest. While multiple studies have assessed test-retest reproducibility of automated white matter parcellations using cortical-parcellation-based strategies, there are no existing studies of test-retest reproducibility of fiber clustering parcellation. In this work, we perform what we believe is the first study of fiber clustering white matter parcellation test-retest reproducibility. The assessment is performed on three test-retest diffusion MRI datasets including a total of 255 subjects across genders, a broad age range (5-82 years), health conditions (autism, Parkinson's disease and healthy subjects), and imaging acquisition protocols (three different sites). A comprehensive evaluation is conducted for a fiber clustering method that leverages an anatomically curated fiber clustering white matter atlas, with comparison to a popular cortical-parcellation-based method. The two methods are compared for the two main white matter parcellation applications of dividing the entire white matter into parcels (i.e., whole brain white matter parcellation) and identifying particular anatomical fiber tracts (i.e., anatomical fiber tract parcellation). Test-retest reproducibility is measured using both geometric and diffusion features, including volumetric overlap (wDice) and relative difference of fractional anisotropy. Our experimental results in general indicate that the fiber clustering method produced more reproducible white matter parcellations than the cortical-parcellation-based method.

Multifold Acceleration of Diffusion MRI via Deep Learning Reconstruction from Slice-Undersampled Data

Diffusion MRI (dMRI), while powerful for characterization of tissue microstructure, suffers from long acquisition time. In this paper, we present a method for effective diffusion MRI reconstruction from slice-undersampled data. Instead of full diffusion-weighted (DW) image volumes, only a subsample of equally-spaced slices need to be acquired. We show that complementary information from DW volumes corresponding to different diffusion wavevectors can be harnessed using graph convolutional neural networks for reconstruction of the full DW volumes. The experimental results indicate a high acceleration factor of up to 5 can be achieved with minimal information loss.

Super-resolution Diffusion Tensor Imaging for Delineating the Facial Nerve in Patients with Vestibular Schwannoma

Objectives  Predicting the course of cranial nerves (CNs) VII and VIII in the cerebellopontine angle on preoperative imaging for vestibular schwannoma (VS) may help guide surgical resection and reduce complications. Diffusion magnetic resonance imaging dMRI is commonly used for this purpose, but is limited by its resolution. We investigate the use of super-resolution reconstruction (SRR), where several different dMRIs are combined into one dataset. We hypothesize that SRR improves the visualization of the CN VII and VIII. Design  Retrospective case review. Setting  Tertiary referral center. SRR was performed on the basis of axial and parasagittal single-shot epiplanar diffusion tensor imaging on a 3.0-tesla MRI scanner. Participants  Seventeen adult patients with suspected neoplasms of the lateral skull base. Main Outcome Measures  We assessed separability of the two distinct nerves on fractional anisotropy (FA) maps, the tractography of the nerves through the cerebrospinal fluid (CSF), and FA in the CSF as a measure of noise. Results  SRR increases separability of the CN VII and VIII (16/17 vs. 0/17, p  = 0.008). Mean FA of CSF surrounding the nerves is significantly lower in SRRs (0.07 ± 0.02 vs. 0.13 ± 0.03 [axial images]/0.14 ± 0.05 [parasagittal images], p  = 0.00003/ p  = 0.00005). Combined scanning times (parasagittal and axial) used for SRR were shorter (8 minute 25 seconds) than a comparable high-resolution scan (15 minute 17 seconds). Conclusion  SRR improves the resolution of CN VII and VIII. The technique can be readily applied in the clinical setting, improving surgical counseling and planning in patients with VS.

Advances in microstructural diffusion neuroimaging for psychiatric disorders

Understanding the neuropathological underpinnings of mental disorders such as schizophrenia, major depression, and bipolar disorder is an essential step towards the development of targeted treatments. Diffusion MRI studies utilizing the diffusion tensor imaging (DTI) model have been extremely successful to date in identifying microstructural brain abnormalities in individuals suffering from mental illness, especially in regions of white matter, although identified abnormalities have been biologically non-specific. Building on DTI's success, in recent years more advanced diffusion MRI methods have been developed and applied to the study of psychiatric populations, with the aim of offering increased sensitivity to subtle neurological abnormalities, as well as improved specificity to candidate pathologies such as demyelination and neuroinflammation. These advanced methods, however, usually come at the cost of prolonged imaging sequences or reduced signal to noise, and they are more difficult to evaluate compared with the more simplified approach taken by the now common DTI model. To date, a limited number of advanced diffusion MRI methods have been employed to study schizophrenia, major depression and bipolar disorder populations. In this review we survey these studies, compare findings across diverse methods, discuss the main benefits and limitations of the different methods, and assess the extent to which the application of more advanced diffusion imaging approaches has led to novel and transformative information with regards to our ability to better understand the etiology and pathology of mental disorders.

Investigation into local white matter abnormality in emotional processing and sensorimotor areas using an automatically annotated fiber clustering in major depressive disorder

This work presents an automatically annotated fiber cluster (AAFC) method to enable identification of anatomically meaningful white matter structures from the whole brain tractography. The proposed method consists of 1) a study-specific whole brain white matter parcellation using a well-established data-driven groupwise fiber clustering pipeline to segment tractography into multiple fiber clusters, and 2) a novel cluster annotation method to automatically assign an anatomical tract annotation to each fiber cluster by employing cortical parcellation information across multiple subjects. The novelty of the AAFC method is that it leverages group-wise information about the fiber clusters, including their fiber geometry and cortical terminations, to compute a tract anatomical label for each cluster in an automated fashion. We demonstrate the proposed AAFC method in an application of investigating white matter abnormality in emotional processing and sensorimotor areas in major depressive disorder (MDD). Seven tracts of interest related to emotional processing and sensorimotor functions are automatically identified using the proposed AAFC method as well as a comparable method that uses a cortical parcellation alone. Experimental results indicate that our proposed method is more consistent in identifying the tracts across subjects and across hemispheres in terms of the number of fibers. In addition, we perform a between-group statistical analysis in 31 MDD patients and 62 healthy subjects on the identified tracts using our AAFC method. We find statistical differences in diffusion measures in local regions within a fiber tract (e.g. 4 fiber clusters within the identified left hemisphere cingulum bundle (consisting of 14 clusters) are significantly different between the two groups), suggesting the ability of our method in identifying potential abnormality specific to subdivisions of a white matter structure.

Joint spatial-angular sparse coding for dMRI with separable dictionaries

Diffusion MRI (dMRI) provides the ability to reconstruct neuronal fibers in the brain, in vivo, by measuring water diffusion along angular gradient directions in q-space. High angular resolution diffusion imaging (HARDI) can produce better estimates of fiber orientation than the popularly used diffusion tensor imaging, but the high number of samples needed to estimate diffusivity requires longer patient scan times. To accelerate dMRI, compressed sensing (CS) has been utilized by exploiting a sparse dictionary representation of the data, discovered through sparse coding. The sparser the representation, the fewer samples are needed to reconstruct a high resolution signal with limited information loss, and so an important area of research has focused on finding the sparsest possible representation of dMRI. Current reconstruction methods however, rely on an angular representation per voxel with added spatial regularization, and so, for non-zero signals, one is required to have at least one non-zero coefficient per voxel. This means that the global level of sparsity must be greater than the number of voxels. In contrast, we propose a joint spatial-angular representation of dMRI that will allow us to achieve levels of global sparsity that are below the number of voxels. A major challenge, however, is the computational complexity of solving a global sparse coding problem over large-scale dMRI. In this work, we present novel adaptations of popular sparse coding algorithms that become better suited for solving large-scale problems by exploiting spatial-angular separability. Our experiments show that our method achieves significantly sparser representations of HARDI than is possible by the state of the art.

A joint space-angle regularization approach for single 4D diffusion image super-resolution

PURPOSE

Low signal-to-noise-ratio and limited scan time of diffusion magnetic resonance imaging (dMRI) in current clinical settings impede obtaining images with high spatial and angular resolution (HSAR) for a reliable fiber reconstruction with fine anatomical details. To overcome this problem, we propose a joint space-angle regularization approach to reconstruct HSAR diffusion signals from a single 4D low resolution (LR) dMRI, which is down-sampled in both 3D-space and q-space.

METHODS

Different from the existing works which combine multiple 4D LR diffusion images acquired using specific acquisition protocols, the proposed method reconstructs HSAR dMRI from only a single 4D dMRI by exploring and integrating two key priors, that is, the nonlocal self-similarity in the spatial domain as a prior to increase spatial resolution and ridgelet approximations in the diffusion domain as another prior to increase the angular resolution of dMRI. To more effectively capture nonlocal self-similarity in the spatial domain, a novel 3D block-based nonlocal means filter is imposed as the 3D image space regularization term which is accurate in measuring the similarity and fast for 3D reconstruction. To reduce computational complexity, we use the L₂ -norm instead of sparsity constraint on the representation coefficients.

RESULTS

Experimental results demonstrate that the proposed method can obtain the HSAR dMRI efficiently with approximately 2% per-voxel root-mean-square error between the actual and reconstructed HSAR dMRI.

CONCLUSION

The proposed approach can effectively increase the spatial and angular resolution of the dMRI which is independent of the acquisition protocol, thus overcomes the inherent resolution limitation of imaging systems.

Motion-robust sub-millimeter isotropic diffusion imaging through motion corrected generalized slice dithered enhanced resolution (MC-gSlider) acquisition

PURPOSE

To develop an efficient MR technique for ultra-high resolution diffusion MRI (dMRI) in the presence of motion.

METHODS

gSlider is an SNR-efficient high-resolution dMRI acquisition technique. However, subject motion is inevitable during a prolonged scan for high spatial resolution, leading to potential image artifacts and blurring. In this study, an integrated technique termed Motion Corrected gSlider (MC-gSlider) is proposed to obtain high-quality, high-resolution dMRI in the presence of large in-plane and through-plane motion. A motion-aware reconstruction with spatially adaptive regularization is developed to optimize the conditioning of the image reconstruction under difficult through-plane motion cases. In addition, an approach for intra-volume motion estimation and correction is proposed to achieve motion correction at high temporal resolution.

RESULTS

Theoretical SNR and resolution analysis validated the efficiency of MC-gSlider with regularization, and aided in selection of reconstruction parameters. Simulations and in vivo experiments further demonstrated the ability of MC-gSlider to mitigate motion artifacts and recover detailed brain structures for dMRI at 860 μm isotropic resolution in the presence of motion with various ranges.

CONCLUSION

MC-gSlider provides motion-robust, high-resolution dMRI with a temporal motion correction sensitivity of 2 s, allowing for the recovery of fine detailed brain structures in the presence of large subject movements.

Incorporation of Prior Knowledge of Signal Behavior Into the Reconstruction to Accelerate the Acquisition of Diffusion MRI Data

Diffusion MRI data are generally acquired using hyperpolarized gases during patient breath-hold, which yields a compromise between achievable image resolution, lung coverage, and number of -values. In this paper, we propose a novel method that accelerates the acquisition of diffusion MRI data by undersampling in both the spatial and -value dimensions and incorporating knowledge about signal decay into the reconstruction (SIDER). SIDER is compared with total variation (TV) reconstruction by assessing its effect on both the recovery of ventilation images and the estimated mean alveolar dimensions (MADs). Both methods are assessed by retrospectively undersampling diffusion data sets ( =8) of healthy volunteers and patients with Chronic Obstructive Pulmonary Disease (COPD) for acceleration factors between x2 and x10. TV led to large errors and artifacts for acceleration factors equal to or larger than x5. SIDER improved TV, with a lower solution error and MAD histograms closer to those obtained from fully sampled data for acceleration factors up to x10. SIDER preserved image quality at all acceleration factors, although images were slightly smoothed and some details were lost at x10. In conclusion, we developed and validated a novel compressed sensing method for lung MRI imaging and achieved high acceleration factors, which can be used to increase the amount of data acquired during breath-hold. This methodology is expected to improve the accuracy of estimated lung microstructure dimensions and provide more options in the study of lung diseases with MRI.

Image formation in diffusion MRI: A review of recent technical developments

Diffusion magnetic resonance imaging (MRI) is a standard imaging tool in clinical neurology, and is becoming increasingly important for neuroscience studies due to its ability to depict complex neuroanatomy (eg, white matter connectivity). Single-shot echo-planar imaging is currently the predominant formation method for diffusion MRI, but suffers from blurring, distortion, and low spatial resolution. A number of methods have been proposed to address these limitations and improve diffusion MRI acquisition. Here, the recent technical developments for image formation in diffusion MRI are reviewed. We discuss three areas of advance in diffusion MRI: improving image fidelity, accelerating acquisition, and increasing the signal-to-noise ratio.

LEVEL OF EVIDENCE

5 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2017;46:646-662.

A New Sparse Representation Framework for Reconstruction of an Isotropic High Spatial Resolution MR Volume From Orthogonal Anisotropic Resolution Scans

In magnetic resonance (MR), hardware limitations, scan time constraints, and patient movement often result in the acquisition of anisotropic 3-D MR images with limited spatial resolution in the out-of-plane views. Our goal is to construct an isotropic high-resolution (HR) 3-D MR image through upsampling and fusion of orthogonal anisotropic input scans. We propose a multiframe super-resolution (SR) reconstruction technique based on sparse representation of MR images. Our proposed algorithm exploits the correspondence between the HR slices and the low-resolution (LR) sections of the orthogonal input scans as well as the self-similarity of each input scan to train pairs of overcomplete dictionaries that are used in a sparse-land local model to upsample the input scans. The upsampled images are then combined using wavelet fusion and error backprojection to reconstruct an image. Features are learned from the data and no extra training set is needed. Qualitative and quantitative analyses were conducted to evaluate the proposed algorithm using simulated and clinical MR scans. Experimental results show that the proposed algorithm achieves promising results in terms of peak signal-to-noise ratio, structural similarity image index, intensity profiles, and visualization of small structures obscured in the LR imaging process due to partial volume effects. Our novel SR algorithm outperforms the nonlocal means (NLM) method using self-similarity, NLM method using self-similarity and image prior, self-training dictionary learning-based SR method, averaging of upsampled scans, and the wavelet fusion method. Our SR algorithm can reduce through-plane partial volume artifact by combining multiple orthogonal MR scans, and thus can potentially improve medical image analysis, research, and clinical diagnosis.

Fast single image super-resolution using estimated low-frequency k-space data in MRI

PURPOSE

Single image super-resolution (SR) is highly desired in many fields but obtaining it is often technically limited in practice. The purpose of this study was to propose a simple, rapid and robust single image SR method in magnetic resonance (MR) imaging (MRI).

METHODS

The idea is based on the mathematical formulation of the intrinsic link in k-space between a given (modulus) low-resolution (LR) image and the desired SR image. The method consists of two steps: 1) estimating the low-frequency k-space data of the desired SR image from a single LR image; 2) reconstructing the SR image using the estimated low-frequency and zero-filled high-frequency k-space data. The method was evaluated on digital phantom images, physical phantom MR images and real brain MR images, and compared with existing SR methods.

RESULTS

The proposed SR method exhibited a good robustness by reaching a clearly higher PSNR (25.77dB) and SSIM (0.991) averaged over different noise levels in comparison with existing edge-guided nonlinear interpolation (EGNI) (PSNR=23.78dB, SSIM=0.983), zero-filling (ZF) (PSNR=24.09dB, SSIM=0.985) and total variation (TV) (PSNR=24.54dB, SSIM=0.987) methods while presenting the same order of computation time as the ZF method but being much faster than the EGNI or TV method. The average PSNR or SSIM over different slice images of the proposed method (PSNR=26.33 dB or SSIM=0.955) was also higher than the EGNI (PSNR=25.07dB or SSIM=0.952), ZF (PSNR=24.97dB or SSIM=0.950) and TV (PSNR=25.70dB or SSIM=0.953) methods, demonstrating its good robustness to variation in anatomical structure of the images. Meanwhile, the proposed method always produced less ringing artifacts than the ZF method, gave a clearer image than the EGNI method, and did not exhibit any blocking effect presented in the TV method. In addition, the proposed method yielded the highest spatial consistency in the inter-slice dimension among the four methods.

CONCLUSIONS

This study proposed a fast, robust and efficient single image SR method with high spatial consistency in the inter-slice dimension for clinical MR images by estimating the low-frequency k-space data of the desired SR image from a single spatial modulus LR image.

Accelerated radial diffusion spectrum imaging using a multi-echo stimulated echo diffusion sequence

PURPOSE

Diffusion spectrum imaging (DSI) provides us non-invasively and robustly with anatomical details of brain microstructure. To achieve sufficient angular resolution, DSI requires a large number of q-space samples, leading to long acquisition times. This need is mitigated here by combining the beneficial properties of Radial q-space sampling for DSI with a Multi-Echo Stimulated Echo Sequence (MESTIM).

METHODS

Full 2D k-spaces for each of several q-space samples, along the same radially outward line in q-space, are acquired in one readout train with one spin and three stimulated echoes. RF flip angles are carefully chosen to distribute spin magnetization over the echoes and the DSI reconstruction is adapted to account for differences in diffusion time among echoes.

RESULTS

Individual datasets and bootstrapped reproducibility analysis demonstrate image quality and SNR of the more-than-twofold-accelerated RDSI MESTIM sequence. Orientation distribution functions (ODF) and tractography results benefit from the longer diffusion times of the latter echoes in the echo train.

CONCLUSION

A MESTIM sequence can be used to shorten RDSI acquisition times significantly without loss of image or ODF quality. Further acceleration is possible by combination with simultaneous multi-slice techniques. Magn Reson Med 79:306-316, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Investigation into Deep Brain Stimulation Lead Designs: A Patient-Specific Simulation Study

New deep brain stimulation (DBS) electrode designs offer operation in voltage and current mode and capability to steer the electric field (EF). The aim of the study was to compare the EF distributions of four DBS leads at equivalent amplitudes (3 V and 3.4 mA). Finite element method (FEM) simulations (n = 38) around cylindrical contacts (leads 3389, 6148) or equivalent contact configurations (leads 6180, SureStim1) were performed using homogeneous and patient-specific (heterogeneous) brain tissue models. Steering effects of 6180 and SureStim1 were compared with symmetric stimulation fields. To make relative comparisons between simulations, an EF isolevel of 0.2 V/mm was chosen based on neuron model simulations (n = 832) applied before EF visualization and comparisons. The simulations show that the EF distribution is largely influenced by the heterogeneity of the tissue, and the operating mode. Equivalent contact configurations result in similar EF distributions. In steering configurations, larger EF volumes were achieved in current mode using equivalent amplitudes. The methodology was demonstrated in a patient-specific simulation around the zona incerta and a "virtual" ventral intermediate nucleus target. In conclusion, lead design differences are enhanced when using patient-specific tissue models and current stimulation mode.

Sparse Reconstruction Techniques in Magnetic Resonance Imaging: Methods, Applications, and Challenges to Clinical Adoption

The family of sparse reconstruction techniques, including the recently introduced compressed sensing framework, has been extensively explored to reduce scan times in magnetic resonance imaging (MRI). While there are many different methods that fall under the general umbrella of sparse reconstructions, they all rely on the idea that a priori information about the sparsity of MR images can be used to reconstruct full images from undersampled data. This review describes the basic ideas behind sparse reconstruction techniques, how they could be applied to improve MRI, and the open challenges to their general adoption in a clinical setting. The fundamental principles underlying different classes of sparse reconstructions techniques are examined, and the requirements that each make on the undersampled data outlined. Applications that could potentially benefit from the accelerations that sparse reconstructions could provide are described, and clinical studies using sparse reconstructions reviewed. Lastly, technical and clinical challenges to widespread implementation of sparse reconstruction techniques, including optimization, reconstruction times, artifact appearance, and comparison with current gold standards, are discussed.

Non Local Spatial and Angular Matching: Enabling higher spatial resolution diffusion MRI datasets through adaptive denoising

Diffusion magnetic resonance imaging (MRI) datasets suffer from low Signal-to-Noise Ratio (SNR), especially at high b-values. Acquiring data at high b-values contains relevant information and is now of great interest for microstructural and connectomics studies. High noise levels bias the measurements due to the non-Gaussian nature of the noise, which in turn can lead to a false and biased estimation of the diffusion parameters. Additionally, the usage of in-plane acceleration techniques during the acquisition leads to a spatially varying noise distribution, which depends on the parallel acceleration method implemented on the scanner. This paper proposes a novel diffusion MRI denoising technique that can be used on all existing data, without adding to the scanning time. We first apply a statistical framework to convert both stationary and non stationary Rician and non central Chi distributed noise to Gaussian distributed noise, effectively removing the bias. We then introduce a spatially and angular adaptive denoising technique, the Non Local Spatial and Angular Matching (NLSAM) algorithm. Each volume is first decomposed in small 4D overlapping patches, thus capturing the spatial and angular structure of the diffusion data, and a dictionary of atoms is learned on those patches. A local sparse decomposition is then found by bounding the reconstruction error with the local noise variance. We compare against three other state-of-the-art denoising methods and show quantitative local and connectivity results on a synthetic phantom and on an in-vivo high resolution dataset. Overall, our method restores perceptual information, removes the noise bias in common diffusion metrics, restores the extracted peaks coherence and improves reproducibility of tractography on the synthetic dataset. On the 1.2 mm high resolution in-vivo dataset, our denoising improves the visual quality of the data and reduces the number of spurious tracts when compared to the noisy acquisition. Our work paves the way for higher spatial resolution acquisition of diffusion MRI datasets, which could in turn reveal new anatomical details that are not discernible at the spatial resolution currently used by the diffusion MRI community.