Maximum spherical mean value filtering for whole-brain QSM.

PURPOSE: To develop a tissue field-filtering algorithm, called maximum spherical mean value (mSMV), for reducing shadow artifacts in QSM of the brain without requiring brain-tissue erosion. THEORY AND METHODS: Residual background field is a major source of shadow artifacts in QSM. The mSMV algorithm filters large field-magnitude values near the border, where the maximum value of the harmonic background field is located. The effectiveness of mSMV for artifact removal was evaluated by comparing existing QSM algorithms in numerical brain simulation as well as using in vivo human data acquired from 11 healthy volunteers and 93 patients. RESULTS: Numerical simulation showed that mSMV reduces shadow artifacts and improves QSM accuracy. Better shadow reduction, as demonstrated by lower QSM variation in the gray matter and higher QSM image quality score, was also observed in healthy subjects and in patients with hemorrhages, stroke, and multiple sclerosis. CONCLUSION: The mSMV algorithm allows QSM maps that are substantially equivalent to those obtained using SMV-filtered dipole inversion without eroding the volume of interest.

Microstructure-Informed Myelin Mapping (MIMM) from Gradient Echo MRI using Stochastic Matching Pursuit.

Quantification of the myelin content of the white matter is important for studying demyelination in neurodegenerative diseases such as Multiple Sclerosis (MS), particularly for longitudinal monitoring. A novel noninvasive MRI method, called Microstructure-Informed Myelin Mapping (MIMM), is developed to quantify the myelin volume fraction (MVF) by utilizing a multi gradient echo sequence (mGRE) and a detailed biophysical model of tissue microstructure. Myelin is modeled as anisotropic negative susceptibility source based on the Hollow Cylindrical Fiber Model (HCFM), and iron as isotropic positive susceptibility source in the extracellular region. Voxels with a range of biophysical parameters are simulated to create a dictionary of MR echo time magnitude signals and total susceptibility values. MRI signals measured using a mGRE sequence are then matched voxel-by-voxel to the created dictionary to obtain the spatial distributions of myelin and iron. Three different MIMM versions are presented to deal with the fiber orientation dependent susceptibility effects of the myelin sheaths: a basic variation, which assumes fiber orientation is an unknown to fit, two orientation informed variations, which assume the fiber orientation distribution is available either from a separate diffusion tensor imaging (DTI) acquisition or from a DTI atlas based fiber orientation map. While all showed a significant linear correlation with the reference method based on T2-relaxometry (p < 0.0001), DTI orientation informed and atlas orientation informed variations reduced overestimation at white matter tracts compared to the basic variation. Finally, the implications and usefulness of attaining an additional iron susceptibility distribution map are discussed. Highlights: novel stochastic matching pursuit algorithm called microstructure-informed myelin mapping (MIMM) is developed to quantify Myelin Volume Fraction (MVF) using Magnetic Resonance Imaging (MRI) and microstructural modeling.utilizes a detailed biophysical model to capture the susceptibility effects on both magnitude and phase to quantify myelin and iron.matter fiber orientation effects are considered for the improved MVF quantification in the major fiber tracts.acquired myelin and iron maps may be utilized to monitor longitudinal disease progress.

SNR and total acquisition time analysis of multi-echo FLASH pulse sequence for current density imaging.

Magnetic Resonance Current Density Imaging (MRCDI) is an imaging modality providing cross-sectional current density (J¯) information inside the body. The clinical applicability of MRCDI is highly dependent on the sensitivity of the acquired noisy current-induced magnetic flux density (B∼z) distributions. Here, a novel analysis is developed to investigate the combined effect of relevant parameters of the RF spoiled gradient echo (FLASH) pulse sequence on the SNR level and the total acquisition time (TAT) of the acquired B∼z images. The proposed analysis then is expanded for a multi-echo FLASH (ME-FLASH) pulse sequence to take advantage of combining the multiple echoes to achieve B∼zcomb distribution with a higher SNR than the one achievable with a single echo acquisition. The optimized sequence parameters to acquire a B∼z distribution with the highest possible SNR for a given acquisition time or the desired SNR in the shortest scan time are estimated using the proposed analysis. The analysis also provides different sets of sequence parameters to acquire B∼z distributions with the same SNR at almost the same TAT. Furthermore, the effects of intensive utilization of the gradients and the magnetohydrodynamic (MHD) flow velocity on the acquired B∼z distribution in MRCDI experiments is investigated. The analytical results of the proposed analysis are validated experimentally using an imaging phantom having the conductivity and the relaxation parameters of the brain white matter tissue.

Low-frequency conductivity tensor imaging with a single current injection using DT-MREIT.

Diffusion tensor-magnetic resonance electrical impedance tomography (DT-MREIT) is an imaging modality to obtain low-frequency anisotropic conductivity distribution employing diffusion tensor imaging and MREIT techniques. DT-MREIT is based on the linear relationship between the conductivity and water self-diffusion tensors in a porous medium, like the brain white matter. Several DT-MREIT studies in the literature provide cross-sectional anisotropic conductivity images of tissue phantoms, canine brain, and the human brain. In these studies, the conductivity tensor images are reconstructed using the diffusion tensor and current density data acquired by injecting two linearly independent current patterns. In this study, a novel reconstruction algorithm is devised for DT-MREIT to reconstruct the conductivity tensor images using a single current injection. Therefore, the clinical applicability of DT-MREIT can be improved by reducing the total acquisition time, the number of current injection cables, and contact electrodes to half by decreasing the number of current injection patterns to one. The proposed method is evaluated utilizing simulated measurements and physical experiments. The results obtained show the successful reconstruction of the anisotropic conductivity distribution using the proposed single current DT-MREIT.