Convolutional network denoising for acceleration of multi-shot diffusion MRI.

Multi-shot echo planar imaging is a promising technique to reduce geometric distortions and increase spatial resolution in diffusion-weighted MRI (DWI), at the expense of increased scan time. Moreover, performing DWI in the body requires multiple repetitions to obtain sufficient signal-to-noise ratio, which further increases the scan time. This work proposes to reduce the number of repetitions and perform denoising of high b-value images using a convolutional network denoising trained on single-shot DWI to accelerate the acquisition of multi-shot DWI. Convolutional network denoising is demonstrated to accelerate the acquisition of 2-shot DWI by a factor of 4 compared to the clinical standard on patients with rectal cancer. Image quality was evaluated using qualitative scores from expert body radiologists between accelerated and non-accelerated acquisition. Additionally, the effect of convolutional network denoising on each image quality score was analyzed using a Wilcoxon signed-rank test. Convolutional network denoising would enable to increase the number of shots without increasing scan time for significant geometric artifact reduction and spatial resolution increase.

Free-breathing high isotropic resolution quantitative susceptibility mapping (QSM) of liver using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction.

PURPOSE: To enable free-breathing and high isotropic resolution liver quantitative susceptibility mapping (QSM) using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction. METHODS: Using 3D multi-echo UTE cones MRI, a respiratory motion was estimated from the k-space center of the imaging data. After sorting the k-space data with estimated motion, respiratory motion state-resolved reconstruction was performed for multi-echo data followed by nonlinear least-squares fitting for proton density fat fraction (PDFF), R 2 * $$ {\mathrm{R}}_2^{\ast } $$ , and fat-corrected B0 field maps. PDFF and B0 field maps were subsequently used for QSM reconstruction. The proposed method was compared with motion-averaged (gridding) reconstruction and conventional 3D multi-echo Cartesian MRI in moving gadolinium phantom and in vivo studies. Region of interest (ROI)-based linear regression analysis was performed on these methods to investigate correlations between gadolinium concentration and QSM in the phantom study and between R 2 * $$ {\mathrm{R}}_2^{\ast } $$ and QSM in in vivo study. RESULTS: Cones with motion-resolved reconstruction showed sharper image quality compared to motion-averaged reconstruction with a substantial reduction of motion artifacts in both moving phantom and in vivo studies. For ROI-based linear regression analysis of the phantom study, susceptibility values from cones with motion-resolved reconstruction ( QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.31 × gadolinium mM + $$ \times {\mathrm{gadolinium}}_{\mathrm{mM}}+ $$ 0.05, R 2 $$ {R}^2 $$ = 0.999) and Cartesian without motion ( QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.32 × gadolinium mM + $$ \times {\mathrm{gadolinium}}_{\mathrm{mM}}+ $$ 0.04, R 2 $$ {R}^2 $$ = 1.000) showed linear relationships with gadolinium concentrations and showed good agreement with each other. For in vivo, motion-resolved reconstruction showed higher goodness of fit ( QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.00261 × R 2 s - 1 * - $$ \times {\mathrm{R}}_{2_{{\mathrm{s}}^{-1}}}^{\ast }- $$ 0.524, R 2 $$ {R}^2 $$ = 0.977) compared to motion-averaged reconstruction ( QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.0021 × R 2 s - 1 * - $$ \times {\mathrm{R}}_{2_{{\mathrm{s}}^{-1}}}^{\ast }- $$ 0.572, R 2 $$ {R}^2 $$ = 0.723) in ROI-based linear regression analysis between R 2 * $$ {\mathrm{R}}_2^{\ast } $$ and QSM. CONCLUSION: Feasibility of free-breathing liver QSM was demonstrated with motion-resolved 3D multi-echo UTE cones MRI, achieving high isotropic resolution currently unachievable in conventional Cartesian MRI.

Accelerated Diffusion-Weighted MRI of Rectal Cancer Using a Residual Convolutional Network.

This work presents a deep-learning-based denoising technique to accelerate the acquisition of high b-value diffusion-weighted MRI for rectal cancer. A denoising convolutional neural network (DCNN) with a combined L1-L2 loss function was developed to denoise high b-value diffusion-weighted MRI data acquired with fewer repetitions (NEX: number of excitations) using the low b-value image as an anatomical guide. DCNN was trained using 85 datasets acquired on patients with rectal cancer and tested on 20 different datasets with NEX = 1, 2, and 4, corresponding to acceleration factors of 16, 8, and 4, respectively. Image quality was assessed qualitatively by expert body radiologists. Reader 1 scored similar overall image quality between denoised images with NEX = 1 and NEX = 2, which were slightly lower than the reference. Reader 2 scored similar quality between NEX = 1 and the reference, while better quality for NEX = 2. Denoised images with fourfold acceleration (NEX = 4) received even higher scores than the reference, which is due in part to the effect of gas-related motion in the rectum, which affects longer acquisitions. The proposed deep learning denoising technique can enable eightfold acceleration with similar image quality (average image quality = 2.8 ± 0.5) and fourfold acceleration with higher image quality (3.0 ± 0.6) than the clinical standard (2.5 ± 0.8) for improved diagnosis of rectal cancer.

Free-breathing R2∗ mapping of hepatic iron overload in children using 3D multi-echo UTE cones MRI.

PURPOSE: To enable motion-robust, ungated, free-breathing R2∗ mapping of hepatic iron overload in children with 3D multi-echo UTE cones MRI. METHODS: A golden-ratio re-ordered 3D multi-echo UTE cones acquisition was developed with chemical-shift encoding (CSE). Multi-echo complex-valued source images were reconstructed via gridding and coil combination, followed by confounder-corrected R2∗ (=1/ T2∗ ) mapping. A phantom containing 15 different concentrations of gadolinium solution (0-300 mM) was imaged at 3T. 3D multi-echo UTE cones with an initial TE of 0.036 ms and Cartesian CSE-MRI (IDEAL-IQ) sequences were performed. With institutional review board approval, 85 subjects (81 pediatric patients with iron overload + 4 healthy volunteers) were imaged at 3T using 3D multi-echo UTE cones with free breathing (FB cones), IDEAL-IQ with breath holding (BH Cartesian), and free breathing (FB Cartesian). Overall image quality of R2∗ maps was scored by 2 blinded experts and compared by a Wilcoxon rank-sum test. For each pediatric subject, the paired R2∗ maps were assessed to determine if a corresponding artifact-free 15 mm region-of-interest (ROI) could be identified at a mid-liver level on both images. Agreement between resulting R2∗ quantification from FB cones and BH/FB Cartesian was assessed with Bland-Altman and linear correlation analyses. RESULTS: ROI-based regression analysis showed a linear relationship between gadolinium concentration and R2∗ in IDEAL-IQ (y = 8.83x - 52.10, R2 = 0.995) as well as in cones (y = 9.19x - 64.16, R2 = 0.992). ROI-based Bland-Altman analysis showed that the mean difference (MD) was 0.15% and the SD was 5.78%. However, IDEAL-IQ R2∗ measurements beyond 200 mM substantially deviated from a linear relationship for IDEAL-IQ (y = 5.85x + 127.61, R2 = 0.827), as opposed to cones (y = 10.87x - 166.96, R2 = 0.984). In vivo, FB cones R2∗ had similar image quality with BH and FB Cartesian in 15 and 42 cases, respectively. FB cones R2∗ had better image quality scores than BH and FB Cartesian in 3 and 21 cases, respectively, where BH/FB Cartesian exhibited severe ghosting artifacts. ROI-based Bland-Altman analyses were 2.23% (MD) and 6.59% (SD) between FB cones and BH Cartesian and were 0.21% (MD) and 7.02% (SD) between FB cones and FB Cartesian, suggesting a good agreement between FB cones and BH (FB) Cartesian R2∗ . Strong linear relationships were observed between BH Cartesian and FB cones (y = 1.00x + 1.07, R2 = 0.996) and FB Cartesian and FB cones (y = 0.98x + 1.68, R2 = 0.999). CONCLUSION: Golden-ratio re-ordered 3D multi-echo UTE Cones MRI enabled motion-robust, ungated, and free-breathing R2∗ mapping of hepatic iron overload, with comparable R2∗ measurements and image quality to BH Cartesian, and better image quality than FB Cartesian.

Dipole modeling of multispectral signal for detecting metallic biopsy markers during MRI-guided breast biopsy: a pilot study.

PURPOSE: During MRI-guided breast biopsy, a metallic biopsy marker is deployed at the biopsy site to guide future interventions. Conventional MRI during biopsy cannot distinguish such markers from biopsy site air, and a post-biopsy mammogram is therefore performed to localize marker placement. The purpose of this pilot study is to develop dipole modeling of multispectral signal (DIMMS) as an MRI alternative to eliminate the cost, inefficiency, inconvenience, and ionizing radiation of a mammogram for biopsy marker localization. METHODS: DIMMS detects and localizes the biopsy marker by fitting the measured multispectral imaging (MSI) signal to the MRI signal model and marker properties. MSI was performed on phantoms containing titanium biopsy markers and air to illustrate the clinical challenge that DIMMS addresses and on 20 patients undergoing MRI-guided breast biopsy to assess DIMMS feasibility for marker detection. DIMMS was compared to conventional MSI field map thresholding, using the post-procedure mammogram as the reference standard. RESULTS: Biopsy markers were detected and localized in 20 of 20 cases using MSI with automated DIMMS post-processing (using a threshold of 0.7) and in 18 of 20 cases using MSI field mapping (using a threshold of 0.65 kHz). CONCLUSION: MSI with DIMMS post-processing is a feasible technique for biopsy marker detection and localization during MRI-guided breast biopsy. With a 2-min MSI scan, DIMMS is a promising MRI alternative to the standard-of-care post-biopsy mammogram.

Cluster analysis of time evolution (CAT) for quantitative susceptibility mapping (QSM) and quantitative blood oxygen level-dependent magnitude (qBOLD)-based oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO2 ) mapping.

PURPOSE: To improve the accuracy of QSM plus quantitative blood oxygen level-dependent magnitude (QSM + qBOLD or QQ)-based mapping of the oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO2 ) using cluster analysis of time evolution (CAT). METHODS: 3D multi-echo gradient echo and arterial spin labeling images were acquired in 11 healthy subjects and 5 ischemic stroke patients. DWI was also carried out on patients. CAT was developed for analyzing signal evolution over TE. QQ-based OEF and CMRO2 were reconstructed with and without CAT, and results were compared using region of interest analysis and a paired t-test. RESULTS: Simulations demonstrated that CAT substantially reduced noise error in QQ-based OEF. In healthy subjects, QQ-based OEF appeared less noisy and more uniform with CAT than without CAT; average OEF with and without CAT in cortical gray matter was 32.7 ± 4.0% and 37.9 ± 4.5%, with corresponding CMRO2 of 148.4 ± 23.8 and 171.4 ± 22.4 µmol/100 g/min, respectively. In patients, regions of low OEF were confined within the ischemic lesions defined on DWI when using CAT, which was not observed without CAT. CONCLUSION: The cluster analysis of time evolution (CAT) significantly improves the robustness of QQ-based OEF against noise.

Cerebral metabolic rate of oxygen (CMRO2 ) mapping by combining quantitative susceptibility mapping (QSM) and quantitative blood oxygenation level-dependent imaging (qBOLD).

PURPOSE: To map the cerebral metabolic rate of oxygen (CMRO2 ) by estimating the oxygen extraction fraction (OEF) from gradient echo imaging (GRE) using phase and magnitude of the GRE data. THEORY AND METHODS: 3D multi-echo gradient echo imaging and perfusion imaging with arterial spin labeling were performed in 11 healthy subjects. CMRO2 and OEF maps were reconstructed by joint quantitative susceptibility mapping (QSM) to process GRE phases and quantitative blood oxygen level-dependent (qBOLD) modeling to process GRE magnitudes. Comparisons with QSM and qBOLD alone were performed using ROI analysis, paired t-tests, and Bland-Altman plot. RESULTS: The average CMRO2 value in cortical gray matter across subjects were 140.4 ± 14.9, 134.1 ± 12.5, and 184.6 ± 17.9 µmol/100 g/min, with corresponding OEFs of 30.9 ± 3.4%, 30.0 ± 1.8%, and 40.9 ± 2.4% for methods based on QSM, qBOLD, and QSM+qBOLD, respectively. QSM+qBOLD provided the highest CMRO2 contrast between gray and white matter, more uniform OEF than QSM, and less noisy OEF than qBOLD. CONCLUSION: Quantitative CMRO2 mapping that fits the entire complex GRE data is feasible by combining QSM analysis of phase and qBOLD analysis of magnitude.

Cardiac quantitative susceptibility mapping (QSM) for heart chamber oxygenation.

PURPOSE: To demonstrate the feasibility of in vivo quantitative susceptibility mapping (QSM) in cardiac MRI and to show that mixed-venous oxygen saturation (SvO2 ) can be measured non-invasively using QSM. METHODS: Electrocardiographic-gated multi-echo 2D gradient echo data were collected at 1.5 T from 14 healthy volunteers during successive breath-holds. Phase wraps and fat chemical shift were removed using a graph-cut-based phase analysis and IDEAL in an iterative approach. The large susceptibility range from air in the lungs to blood in the heart was addressed by using the preconditioning approach in the dipole field inversion. SvO2 was calculated based on the difference in blood susceptibility between the right ventricle (RV) and left ventricle (LV). Cardiac QSM quality was assessed by two independent readers. RESULTS: Nine out of fourteen volunteers (64%) yielded interpretable cardiac QSM. QSM maps showed strong differential contrast between RV and LV blood with RV blood having higher susceptibility values (291.5 ± 32.4 ppb), which correspond to 78.3 ± 2.3% SvO2 . CONCLUSION: In vivo cardiac QSM is feasible and can be used to measure SvO2 , but improvements in data acquisition are needed. Magn Reson Med 79:1545-1552, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Coherence enhancement in quantitative susceptibility mapping by means of anisotropic weighting in morphology enabled dipole inversion.

PURPOSE: To investigate an anisotropic structural prior in morphology enabled dipole inversion (MEDI) for improving accuracy in quantitative susceptibility mapping (QSM). THEORY AND METHODS: Anisotropic weighting (AW) was devised and implemented to incorporate orientation information into the edge agreement in the MEDI method. AW performance was compared with isotropic weighting by testing and validating on in vivo brain multiple orientation MRI data using COSMOS and the (33) component of the susceptibility tensor as reference. RESULTS: Suppressing streaking artifacts, AW improved not only QSM image quality but also accuracy in terms of RMSE (root mean square error), HFEN (high frequency error norm), SSIM (structural similarity index), and GDA (gradient direction agreement). In addition, it outperformed isotropic weighting in region of interest-based analysis. From a computational perspective, AW was as fast as isotropic weighting, taking approximately the same central processing unit times. CONCLUSION: Using AW in MEDI improves QSM accuracy compared with isotropic weighting. Magn Reson Med 79:1172-1180, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Quantitative Susceptibility Mapping (QSM) Algorithms: Mathematical Rationale and Computational Implementations.

Quantitative susceptibility mapping (QSM) solves the magnetic field-to-magnetization (tissue susceptibility) inverse problem under conditions of noisy and incomplete field data acquired using magnetic resonance imaging. Therefore, sophisticated algorithms are necessary to treat the ill-posed nature of the problem and are reviewed here. The forward problem is typically presented as an integral form, where the field is the convolution of the dipole kernel and tissue susceptibility distribution. This integral form can be equivalently written as a partial differential equation (PDE). Algorithmic challenges are to reduce streaking and shadow artifacts characterized by the fundamental solution of the PDE. Bayesian maximum a posteriori estimation can be employed to solve the inverse problem, where morphological and relevant biomedical knowledge (specific to the imaging situation) are used as priors. As the cost functions in Bayesian QSM framework are typically convex, solutions can be robustly computed using a gradient-based optimization algorithm. Moreover, one can not only accelerate Bayesian QSM, but also increase its effectiveness at reducing shadows using prior knowledge based preconditioners. Improving the efficiency of QSM is under active development, and a rigorous analysis of preconditioning needs to be carried out for further investigation.Quantitative susceptibility mapping (QSM) solves the magnetic field-to-magnetization (tissue susceptibility) inverse problem under conditions of noisy and incomplete field data acquired using magnetic resonance imaging. Therefore, sophisticated algorithms are necessary to treat the ill-posed nature of the problem and are reviewed here. The forward problem is typically presented as an integral form, where the field is the convolution of the dipole kernel and tissue susceptibility distribution. This integral form can be equivalently written as a partial differential equation (PDE). Algorithmic challenges are to reduce streaking and shadow artifacts characterized by the fundamental solution of the PDE. Bayesian maximum a posteriori estimation can be employed to solve the inverse problem, where morphological and relevant biomedical knowledge (specific to the imaging situation) are used as priors. As the cost functions in Bayesian QSM framework are typically convex, solutions can be robustly computed using a gradient-based optimization algorithm. Moreover, one can not only accelerate Bayesian QSM, but also increase its effectiveness at reducing shadows using prior knowledge based preconditioners. Improving the efficiency of QSM is under active development, and a rigorous analysis of preconditioning needs to be carried out for further investigation.

Primal-dual and forward gradient implementation for quantitative susceptibility mapping.

PURPOSE: To investigate the computational aspects of the prior term in quantitative susceptibility mapping (QSM) by (i) comparing the Gauss-Newton conjugate gradient (GNCG) algorithm that uses numerical conditioning (ie, modifies the prior term) with a primal-dual (PD) formulation that avoids this, and (ii) carrying out a comparison between a central and forward difference scheme for the discretization of the prior term. THEORY AND METHODS: A spatially continuous formulation of the regularized QSM inversion problem and its PD formulation were derived. The Chambolle-Pock algorithm for PD was implemented and its convergence behavior was compared with that of GNCG for the original QSM. Forward and central difference schemes were compared in terms of the presence of checkerboard artifacts. All methods were tested and validated on a gadolinium phantom, ex vivo brain blocks, and in vivo brain MRI data with respect to COSMOS. RESULTS: The PD approach provided a faster convergence rate than GNCG. The GNCG convergence rate slowed considerably with smaller (more accurate) values of the conditioning parameter. Using a forward difference suppressed the checkerboard artifacts in QSM, as compared with the central difference. The accuracy of PD and GNCG were validated based on excellent correlation with COSMOS. CONCLUSIONS: The PD approach with forward difference for the gradient showed improved convergence and accuracy over the GNCG method using central difference. Magn Reson Med 78:2416-2427, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Preconditioned total field inversion (TFI) method for quantitative susceptibility mapping.

PURPOSE: To investigate systematic errors in traditional quantitative susceptibility mapping (QSM) where background field removal and local field inversion (LFI) are performed sequentially, to develop a total field inversion (TFI) QSM method to reduce these errors, and to improve QSM quality in the presence of large susceptibility differences. THEORY AND METHODS: The proposed TFI is a single optimization problem which simultaneously estimates the background and local fields, preventing error propagation from background field removal to QSM. To increase the computational speed, a new preconditioner is introduced and analyzed. TFI is compared with the traditional combination of background field removal and LFI in a numerical simulation and in phantom, 5 healthy subjects, and 18 patients with intracerebral hemorrhage. RESULTS: Compared with the traditional method projection onto dipole fields+LFI, preconditioned TFI substantially reduced error in QSM along the air-tissue boundaries in simulation, generated high-quality in vivo QSM within similar processing time, and suppressed streaking artifacts in intracerebral hemorrhage QSM. Moreover, preconditioned TFI was capable of generating QSM for the entire head including the brain, air-filled sinus, skull, and fat. CONCLUSION: Preconditioned total field inversion improves the accuracy of QSM over the traditional method where background and local fields are separately estimated. Magn Reson Med 78:303-315, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Adaptive comb filtering for motion artifact reduction from PPG with a structure of adaptive lattice IIR notch filter.

Photoplethysmographic (PPG) signal can provide important information about cardiovascular and respiratory conditions of individuals in a hospital or daily life. However, PPG can be distorted by motion artifacts significantly. Therefore, the reduction of the effects of motion artifacts is very important procedure for monitoring cardio-respiratory system by PPG. There have been many adaptive techniques to reduce motion artifacts from PPG signal including normalized least mean squares (NLMS) method, recursive least squares (RLS) filter, and Kalman filter. In the present study, we propose the adaptive comb filter (ACF) for reducing the effects of motion artifacts from PPG signal. ACF with adaptive lattice infinite impulse response (IIR) notch filter (ALNF) successfully reduced the motion artifacts from the quasi-periodic PPG signal.