Efficient parallel reconstruction for high resolution multishot spiral diffusion data with low rank constraint

High-fidelity mesoscale in-vivo diffusion MRI through gSlider-BUDA and circular EPI with S-LORAKS reconstruction

PURPOSE

To develop a high-fidelity diffusion MRI acquisition and reconstruction framework with reduced echo-train-length for less T2* image blurring compared to typical highly accelerated echo-planar imaging (EPI) acquisitions at sub-millimeter isotropic resolution.

METHODS

We first proposed a circular-EPI trajectory with partial Fourier sampling on both the readout and phase-encoding directions to minimize the echo-train-length and echo time. We then utilized this trajectory in an interleaved two-shot EPI acquisition with reversed phase-encoding polarity, to aid in the correction of off-resonance-induced image distortions and provide complementary k-space coverage in the missing partial Fourier regions. Using model-based reconstruction with structured low-rank constraint and smooth phase prior, we corrected the shot-to-shot phase variations across the two shots and recover the missing k-space data. Finally, we combined the proposed acquisition/reconstruction framework with an SNR-efficient RF-encoded simultaneous multi-slab technique, termed gSlider, to achieve high-fidelity 720 µm and 500 µm isotropic resolution in-vivo diffusion MRI.

RESULTS

Both simulation and in-vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide distortion-corrected diffusion imaging at the mesoscale with markedly reduced T2*-blurring. The in-vivo results of 720 µm and 500 µm datasets show high-fidelity diffusion images with reduced image blurring and echo time using the proposed approaches.

CONCLUSIONS

The proposed method provides high-quality distortion-corrected diffusion-weighted images with ∼40% reduction in the echo-train-length and T2* blurring at 500µm-isotropic-resolution compared to standard multi-shot EPI.

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.

Anomaly detection of complex magnetic measurements using structured Hankel low-rank modeling and singular value decomposition

The magnetic anomalies generated by the ferromagnetic targets are usually buried within uncontrollable interference sources, such as the power frequency and random noises. In particular, the variability of the geomagnetic field and the low signal-to-noise ratio (SNR) of the magnetic anomalies cannot be avoided. In this paper, to improve the performance of magnetic anomaly detection (MAD) with a low SNR, we propose a novel structured low-rank (SLR) decomposition-based MAD method. In addition, a new framework based on the SLR and singular value decomposition (SVD) is constructed, dubbed SLR-SVD, and the corresponding working principle and implemented strategy are elaborated. Through comparing the SLR-SVD with two state-of-the-art methods, including principal component analysis and SVD, the results demonstrate that the proposed SLR-SVD can not only suppress the noise sufficiently, i.e., improving 55.26% approximately of the SNR, but also retain more boundary information of magnetic anomalies, i.e., decreasing approximately 68.05% of the mean squared error and improving approximately 28.47% of the structural similarity index.

Distortion-Free Diffusion Imaging Using Self-Navigated Cartesian Echo-Planar Time Resolved Acquisition and Joint Magnitude and Phase Constrained Reconstruction

Echo-planar time resolved imaging (EPTI) is an effective approach for acquiring high-quality distortion-free images with a multi-shot EPI (ms-EPI) readout. As with traditional ms-EPI acquisitions, inter-shot phase variations present a main challenge when incorporating EPTI into a diffusion-prepared pulse sequence. The aim of this study is to develop a self-navigated Cartesian EPTI-based (scEPTI) acquisition together with a magnitude and phase constrained reconstruction for distortion-free diffusion imaging. A self-navigated Cartesian EPTI-based diffusion-prepared pulse sequence is designed. The different phase components in EPTI diffusion signal are analyzed and an approach to synthesize a fully phase-matched navigator for the inter-shot phase correction is demonstrated. Lastly, EPTI contains richer magnitude and phase information than conventional ms-EPI, such as the magnitude and phase correlations along the temporal dimension. The potential of these magnitude and phase correlations to enhance the reconstruction is explored. The reconstruction results with and without phase matching and with and without phase or magnitude constraints are compared. Compared with reconstruction without phase matching, the proposed phase matching method can improve the accuracy of inter-shot phase correction and reduce signal corruption in the final diffusion images. Magnitude constraints further improve image quality by suppressing the background noise and thereby increasing SNR, while phase constraints can mitigate possible image blurring from adding magnitude constraints. The high-quality distortion-free diffusion images and simultaneous diffusion-relaxometry imaging capacity provided by the proposed EPTI design represent a highly valuable tool for both clinical and neuroscientific assessments of tissue microstructure.

A modular magneto-inductive sensor for low vector magnetic field measurements

The low magnetic field measurement has been utilized since ancient times in order to find economic resources, to detect magnetic anomalies, etc. In this case, the vector magnetic survey can simultaneously obtain the modulus and direction information of the magnetic field, which can contribute to obtaining more precise information and characteristics of magnetic field resources. This paper is concerned with the potential to exploit the signals of vector magnetic field measurement with a magneto-inductive (MI) sensor. To evaluate the capability of the MI sensor, a test platform is set up and its performance, including the noise floor, the resolution, and the sensitivity, is comprehensively characterized. Furthermore, a comparative geomagnetic observation and magnetic anomaly detection among the proposed MI sensor, a high-precision Overhauser sensor, and a commonly used and accepted commercial MI sensor are conducted. The experimental results identify the capability of the proposed MI sensor in weak magnetic detection.

Efficient T2 mapping with blip-up/down EPI and gSlider-SMS (T2 -BUDA-gSlider)

PURPOSE

To rapidly obtain high isotropic-resolution T₂ maps with whole-brain coverage and high geometric fidelity.

METHODS

A T₂ blip-up/down EPI acquisition with generalized slice-dithered enhanced resolution (T₂ -BUDA-gSlider) is proposed. A RF-encoded multi-slab spin-echo (SE) EPI acquisition with multiple TEs was developed to obtain high SNR efficiency with reduced TR. This was combined with an interleaved 2-shot EPI acquisition using blip-up/down phase encoding. An estimated field map was incorporated into the joint multi-shot EPI reconstruction with a structured low rank constraint to achieve distortion-free and robust reconstruction for each slab without navigation. A Bloch simulated subspace model was integrated into gSlider reconstruction and used for T₂ quantification.

RESULTS

In vivo results demonstrated that the T₂ values estimated by the proposed method were consistent with gold standard spin-echo acquisition. Compared to the reference 3D fast spin echo (FSE) images, distortion caused by off-resonance and eddy current effects were effectively mitigated.

CONCLUSION

BUDA-gSlider SE-EPI acquisition and gSlider-subspace joint reconstruction enabled distortion-free whole-brain T₂ mapping in 2 min at ~1 mm³ isotropic resolution, which could bring significant benefits to related clinical and neuroscience applications.

Distortion-free, high-isotropic-resolution diffusion MRI with gSlider BUDA-EPI and multicoil dynamic B0 shimming

PURPOSE

We combine SNR-efficient acquisition and model-based reconstruction strategies with newly available hardware instrumentation to achieve distortion-free in vivo diffusion MRI of the brain at submillimeter-isotropic resolution with high fidelity and sensitivity on a clinical 3T scanner.

METHODS

We propose blip-up/down acquisition (BUDA) for multishot EPI using interleaved blip-up/blip-down phase encoding and incorporate B₀ forward-modeling into structured low-rank reconstruction to enable distortion-free and navigator-free diffusion MRI. We further combine BUDA-EPI with an SNR-efficient simultaneous multislab acquisition (generalized slice-dithered enhanced resolution ["gSlider"]), to achieve high-isotropic-resolution diffusion MRI. To validate gSlider BUDA-EPI, whole-brain diffusion data at 860-μm and 780-μm data sets were acquired. Finally, to improve the conditioning and minimize noise penalty in BUDA reconstruction at very high resolutions where B₀ inhomogeneity can have a detrimental effect, the level of B₀ inhomogeneity was reduced by incorporating slab-by-slab dynamic shimming with a 32-channel AC/DC coil into the acquisition. Whole-brain 600-μm diffusion data were then acquired with this combined approach of gSlider BUDA-EPI with dynamic shimming.

RESULTS

The results of 860-μm and 780-μm datasets show high geometry fidelity with gSlider BUDA-EPI. With dynamic shimming, the BUDA reconstruction's noise penalty was further alleviated. This enables whole-brain 600-μm isotropic resolution diffusion imaging with high image quality.

CONCLUSIONS

The gSlider BUDA-EPI method enables high-quality, distortion-free diffusion imaging across the whole brain at submillimeter resolution, where the use of multicoil dynamic B₀ shimming further improves reconstruction performance, which can be particularly useful at very high resolutions.

Modeling and reduction of the gradient effect of a magnetic sensor for an aeromagnetic compensator

The traditional aeromagnetic compensation model does not consider the negative effects of the gradient of the geomagnetic background field. Especially in emergency aeromagnetic surveys, this gradient effect is obvious and can reduce the performance of magnetic compensation. To solve this problem, we propose a method of reducing the gradient effect of the magnetic sensor for aeromagnetic compensation. A magnetic gradient model of the Earth's background field is established for the aeromagnetic compensation, and the gradient effect is then suppressed using altitude data. An experimental platform was built, and flight tests were conducted. Comparison results show that the adoption of the proposed method increased the improvement ratio (IR) of our compensator on two flights from 20.4 and 21.3 to 25.9 and 26.0, with the growth rates of the IR, respectively, reaching 27% and 22%.

Model-Based Deep Learning for Reconstruction of Joint k-q Under-sampled High Resolution Diffusion MRI

We propose a model-based deep learning architecture for the reconstruction of highly accelerated diffusion magnetic resonance imaging (MRI) that enables high resolution imaging. The proposed reconstruction jointly recovers all the diffusion weighted images in a single step from a joint k-q under-sampled acquisition in a parallel MRI setting. We propose the novel use of a pre-trained denoiser as a regularizer in a model-based reconstruction for the recovery of highly under-sampled data. Specifically, we designed the denoiser based on a general diffusion MRI tissue microstructure model for multi-compartmental modeling. By using a wide range of biologically plausible parameter values for the multi-compartmental microstructure model, we simulated diffusion signal that spans the entire microstructure parameter space. A neural network was trained in an unsupervised manner using an autoencoder to learn the diffusion MRI signal subspace. We employed the autoencoder in a model-based reconstruction and show that the autoencoder provides a strong denoising prior to recover the q-space signal. We show reconstruction results on a simulated brain dataset that shows high acceleration capabilities of the proposed method.

Efficient noise reduction for the free induction decay signal from a proton precession magnetometer with time-frequency peak filtering

The proton precession magnetometer (PPM) is a commonly used device to measure the varying magnetic field. Since the frequency of the PPM sensing free induction decay (FID) signal is proportional to the magnetic field, the signal-to-noise ratio (SNR) is always a critical issue that influences the measurement accuracy severely due to the external interferences such as harmonic noise and random noise. In this study, to boost the SNR of the FID signal, an effective filtering algorithm based on time-frequency peak filtering (TFPF) analyzed with pseudo-Wigner-Ville distribution (PWVD) is proposed. Through pre-treating the collected noisy FID signal with frequency modulation and instantaneous frequency estimation using the peak value of the time-frequency characterization, the embedded noise can be decorrelated and the relative pure FID signal can be detected regardless of the impact of varying noise levels. The superiority of the proposed synaptic noise reduction framework, namely, TFPF-PWVD, was found by comparing it with state-of-the-art approaches under the same conditions. The results illustrated that even though in a strong-noisy scenario, the proposed TFPF-PWVD based approach still achieved the best SNR for the yielded sensing FID and the minimum standard deviation for the observed magnetic field data, which can enhance the geomagnetic measuring performance of a PPM.

Towed Overhauser marine magnetometer for weak magnetic anomaly detection in severe ocean conditions

A towed Overhauser marine geomagnetic magnetometer used for weak magnetic anomaly detection in severe ocean conditions is studied to investigate means to reduce the negative effect of dynamic behavior and magnetic noise associated with ocean waves. For the dynamic effect, a continuous polarization workflow is proposed to enhance the free-induction-decay signal, and then, a multi-angle pickup coil and a self-tracking programmable amplifier are used to further reduce the adverse effect caused by uncontrollable changes in the towfish attitude on the signal quality. Furthermore, to achieve adaptive suppression of magnetic noise in different ocean conditions and areas, a modified adaptive Kalman algorithm is assessed. In addition, an optimized Overhauser sensor and a towfish were developed. Overall, the experimental results show that the sensor can effectively suppress the dynamic effect and magnetic noise. Regarding the magnetic sensitivity, uncertainty and range are 12 pT/Hz^1/2@1Hz and 0.21 nT and 20 000 nT-100 000 nT, respectively. Moreover, underwater testing was performed to verify the function and the detection of the magnetic anomaly.

Suppression of the negative effect of abnormal data based on the Hough transform and application to the magnetic compensation of airborne optically-pumped magnetometer data

A magnetic sensor of the optically pumped magnetometer may enter a dead zone during an aeromagnetic survey, resulting in unavoidable abnormal data and seriously affecting the effect of aeromagnetic compensation. We propose a fast discrimination and culling method based on the Hough transform to prevent abnormal data from participating in the coefficient estimation. In the proposed method, the parameter space is partitioned into small buckets and the most frequently passed region of sinusoidal curves is detected to eliminate abnormal data. Although the conventional method performs similarly to the proposed method when there is only 1% abnormal data, it is theoretically shown that the proposed method has a better goodness of fit of 0.9518, compared with a value of 0.1956 for the conventional method, in the presence of 45% abnormal data. Furthermore, we construct an experimental platform and conduct a flight test in which the proposed method has an improvement ratio of 4.11 compared with a value of 0.34 for the conventional method.

High-fidelity, high-isotropic-resolution diffusion imaging through gSlider acquisition with B 1 + and T1 corrections and integrated ΔB0 /Rx shim array

PURPOSE

B 1 + and T₁ corrections and dynamic multicoil shimming approaches were proposed to improve the fidelity of high-isotropic-resolution generalized slice-dithered enhanced resolution (gSlider) diffusion imaging.

METHODS

An extended reconstruction incorporating B 1 + inhomogeneity and T₁ recovery information was developed to mitigate slab-boundary artifacts in short-repetition time (TR) gSlider acquisitions. Slab-by-slab dynamic B₀ shimming using a multicoil integrated ΔB₀ /Rx shim array and high in-plane acceleration (R_inplane = 4) achieved with virtual-coil GRAPPA were also incorporated into a 1-mm isotropic resolution gSlider acquisition/reconstruction framework to achieve a significant reduction in geometric distortion compared to single-shot echo planar imaging (EPI).

RESULTS

The slab-boundary artifacts were alleviated by the proposed B 1 + and T₁ corrections compared to the standard gSlider reconstruction pipeline for short-TR acquisitions. Dynamic shimming provided >50% reduction in geometric distortion compared to conventional global second-order shimming. One-millimeter isotropic resolution diffusion data show that the typically problematic temporal and frontal lobes of the brain can be imaged with high geometric fidelity using dynamic shimming.

CONCLUSIONS

The proposed B 1 + and T₁ corrections and local-field control substantially improved the fidelity of high-isotropic-resolution diffusion imaging, with reduced slab-boundary artifacts and geometric distortion compared to conventional gSlider acquisition and reconstruction. This enabled high-fidelity whole-brain 1-mm isotropic diffusion imaging with 64 diffusion directions in 20 min using a 3T clinical scanner.

A fusion of principal component analysis and singular value decomposition based multivariate denoising algorithm for free induction decay transversal data

The free induction decay (FID) transversal data determines the measurement accuracy of time-dependent geomagnetic fields, whereas the conservation of clean components and removal of noise cannot be easily achieved for this kind of data. Even though numerous techniques have been proven to be effective in improving the signal-to-noise ratio by filtering out frequency bands, how to efficiently reduce noise is still a crucial issue due to several restrictions, e.g., prior information requirement, stationary data assumption. To end this, a new multivariate algorithm based on the fusion of principal component analysis (PCA) and singular value decomposition (SVD), namely, principal component analysis and decomposition (PCAD), was presented. This novel algorithm aims to reduce noise as well as cancel the interference of FID transversal data. Specifically, the PCAD algorithm is able to obtain the dominant principal components of the FID and that of the noise floor by PCA, in which an optimal number of subspaces could be retained via a cumulative percent of variance criterion. Furthermore, the PCA was combined with an SVD filter whose singular values corresponding to the interferences were identified, and then the noise was suppressed by nulling the corresponding singular values, which was able to achieve an optimum trade-off between the preservation of pure FID data and the denoising efficiency. Our proposed PCAD algorithm was compared with the widely used filter methods via extensive experiments on synthetic and real FID transversal data under different noise levels. The results demonstrated that this method can preserve the FID transversal data better and shows a significant improvement in noise suppression.

Motion-robust reconstruction of multishot diffusion-weighted images without phase estimation through locally low-rank regularization

PURPOSE

The goal of this work is to propose a motion robust reconstruction method for diffusion-weighted MRI that resolves shot-to-shot phase mismatches without using phase estimation.

METHODS

Assuming that shot-to-shot phase variations are slowly varying, spatial-shot matrices can be formed using a local group of pixels to form columns, in which each column is from a different shot (excitation). A convex model with a locally low-rank constraint on the spatial-shot matrices is proposed. In vivo brain and breast experiments were performed to evaluate the performance of the proposed method.

RESULTS

The proposed method shows significant benefits when the motion is severe, such as for breast imaging. Furthermore, the resulting images can be used for reliable phase estimation in the context of phase-estimation-based methods to achieve even higher image quality.

CONCLUSION

We introduced the shot-locally low-rank method, a reconstruction technique for multishot diffusion-weighted MRI without explicit phase estimation. In addition, its motion robustness can be beneficial to neuroimaging and body imaging.

Combining SENSE and reduced field-of-view for high-resolution diffusion weighted magnetic resonance imaging

BACKGROUND

In diffusion-weighted magnetic resonance imaging (DWI) using single-shot echo planar imaging (ss-EPI), both reduced field-of-view (FOV) excitation and sensitivity encoding (SENSE) alone can increase in-plane resolution to some degree. However, when the two techniques are combined to further increase resolution without pronounced geometric distortion, the resulted images are often corrupted by high level of noise and artifact due to the numerical restriction in SENSE. Hence, this study is aimed to provide a reconstruction method to deal with this problem.

METHODS

The proposed reconstruction method was developed and implemented to deal with the high level of noise and artifact in the combination of reduced FOV imaging and traditional SENSE, in which all the imaging data were considered jointly by incorporating the motion induced phase variations among excitations. The in vivo human spine diffusion images from ten subjects were acquired at 1.5 T and reconstructed using the proposed method, and compared with SENSE magnitude average results for a range of reduction factors in reduced FOV. These images were evaluated by two radiologists using visual scores (considering distortion, noise and artifact levels) from 1 to 10.

RESULTS

The proposed method was able to reconstruct images with greatly reduced noise and artifact compared to SENSE magnitude average. The mean g-factors were maintained close to 1 along with enhanced signal-to-noise ratio efficiency. The image quality scores of the proposed method were significantly higher (P < 0.01) than SENSE magnitude average for all the evaluated reduction factors.

CONCLUSION

The proposed method can improve the combination of SENSE and reduced FOV for high-resolution ss-EPI DWI with reduced noise and artifact.

A diffusion-matched principal component analysis (DM-PCA) based two-channel denoising procedure for high-resolution diffusion-weighted MRI

Over the past several years, significant efforts have been made to improve the spatial resolution of diffusion-weighted imaging (DWI), aiming at better detecting subtle lesions and more reliably resolving white-matter fiber tracts. A major concern with high-resolution DWI is the limited signal-to-noise ratio (SNR), which may significantly offset the advantages of high spatial resolution. Although the SNR of DWI data can be improved by denoising in post-processing, existing denoising procedures may potentially reduce the anatomic resolvability of high-resolution imaging data. Additionally, non-Gaussian noise induced signal bias in low-SNR DWI data may not always be corrected with existing denoising approaches. Here we report an improved denoising procedure, termed diffusion-matched principal component analysis (DM-PCA), which comprises 1) identifying a group of (not necessarily neighboring) voxels that demonstrate very similar magnitude signal variation patterns along the diffusion dimension, 2) correcting low-frequency phase variations in complex-valued DWI data, 3) performing PCA along the diffusion dimension for real- and imaginary-components (in two separate channels) of phase-corrected DWI voxels with matched diffusion properties, 4) suppressing the noisy PCA components in real- and imaginary-components, separately, of phase-corrected DWI data, and 5) combining real- and imaginary-components of denoised DWI data. Our data show that the new two-channel (i.e., for real- and imaginary-components) DM-PCA denoising procedure performs reliably without noticeably compromising anatomic resolvability. Non-Gaussian noise induced signal bias could also be reduced with the new denoising method. The DM-PCA based denoising procedure should prove highly valuable for high-resolution DWI studies in research and clinical uses.

Robust sliding-window reconstruction for Accelerating the acquisition of MR fingerprinting

PURPOSE

To develop a method for accelerated and robust MR fingerprinting (MRF) with improved image reconstruction and parameter matching processes.

THEORY AND METHODS

A sliding-window (SW) strategy was applied to MRF, in which signal and dictionary matching was conducted between fingerprints consisting of mixed-contrast image series reconstructed from consecutive data frames segmented by a sliding window, and a precalculated mixed-contrast dictionary. The effectiveness and performance of this new method, dubbed SW-MRF, was evaluated in both phantom and in vivo. Error quantifications were conducted on results obtained with various settings of SW reconstruction parameters.

RESULTS

Compared with the original MRF strategy, the results of both phantom and in vivo experiments demonstrate that the proposed SW-MRF strategy either provided similar accuracy with reduced acquisition time, or improved accuracy with equal acquisition time. Parametric maps of T₁ , T₂ , and proton density of comparable quality could be achieved with a two-fold or more reduction in acquisition time. The effect of sliding-window width on dictionary sensitivity was also estimated.

CONCLUSION

The novel SW-MRF recovers high quality image frames from highly undersampled MRF data, which enables more robust dictionary matching with reduced numbers of data frames. This time efficiency may facilitate MRF applications in time-critical clinical settings. Magn Reson Med 78:1579-1588, 2017. © 2016 International Society for Magnetic Resonance in Medicine.