Search

Purpose: Parallel imaging and compressed sensing reconstructions of large MRI datasets often have a prohibitive computational cost that bottlenecks clinical deployment, especially for 3D non-Cartesian acquisitions. One common approach is to reduce the number of coil channels actively used during reconstruction as in coil compression. While effective for Cartesian imaging, coil compression inherently loses signal energy, producing shading artifacts that compromise image quality for 3D non-Cartesian imaging. We propose coil sketching, a general and versatile method for computationally-efficient iterative MR image reconstruction. Theory and Methods: We based our method on randomized sketching algorithms, a type of large-scale optimization algorithms well established in the fields of machine learning and big data analysis. We adapt the sketching theory to the MRI reconstruction problem via a structured sketching matrix that, similar to coil compression, considers high-energy virtual coils obtained from principal component analysis. But, unlike coil compression, it also considers random linear combinations of the remaining low-energy coils, effectively leveraging information from all coils. Results: First, we performed ablation experiments to validate the sketching matrix design on both Cartesian and non-Cartesian datasets. The resulting design yielded both improved computational efficiency and preserved signal-to-noise ratio (SNR) as measured by the inverse g-factor. Then, we verified the efficacy of our approach on high-dimensional non-Cartesian 3D cones datasets, where coil sketching yielded up to three-fold faster reconstructions with equivalent image quality. Conclusion: Coil sketching is a general and versatile reconstruction framework for computationally fast and memory-efficient reconstruction., Comment: 19 pages, 7 figures, 3 tables

Accelerated MRI protocols routinely involve a predefined sampling pattern that undersamples the k-space. Finding an optimal pattern can enhance the reconstruction quality, however this optimization is a challenging task. To address this challenge, we introduce a novel deep learning framework, AutoSamp, based on variational information maximization that enables joint optimization of sampling pattern and reconstruction of MRI scans. We represent the encoder as a non-uniform Fast Fourier Transform that allows continuous optimization of k-space sample locations on a non-Cartesian plane, and the decoder as a deep reconstruction network. Experiments on public MRI datasets show improved reconstruction quality of the proposed AutoSamp method over the prevailing variable density and variable density Poisson disc sampling. We demonstrate that our data-driven sampling optimization method achieves 4.4dB, 2.0dB, 0.75dB, 0.7dB PSNR improvements over reconstruction with Poisson Disc masks for acceleration factors of R = 5, 10, 15, 25, respectively. Furthermore, we analyze the characteristics of the learned sampling patterns with respect to changes in acceleration factor, measurement noise, underlying anatomy, and coil sensitivities. We show that all these factors contribute to the optimization result by affecting the sampling density, k-space coverage and point spread functions of the learned sampling patterns.

BackgroundNon-invasive assessment of the hemodynamic changes of cirrhosis might help guide management of patients with liver disease but are currently limited.PurposeTo determine whether free-breathing 4D flow MRI can be used to quantify the hemodynamic effects of cirrhosis and introduce hydraulic circuit indexes of severity.Study typeRetrospective.PopulationForty-seven patients including 26 with cirrhosis.Field strength/sequence3 T/free-breathing 4D flow MRI with soft gating and golden-angle view ordering.AssessmentMeasurements of the supra-celiac abdominal aorta, supra-renal abdominal aorta (SRA), celiac trunk (CeT), superior mesenteric artery (SMA), splenic artery (SpA), common hepatic artery (CHA), portal vein (PV), and supra-renal inferior vena cava (IVC) were made by two radiologists. Measures of hepatic vascular resistance (hepatic arterial relative resistance [HARR]; portal resistive index [PRI]) were proposed and calculated.Statistical analysisBland-Altman, Pearson's correlation, Tukey's multiple comparison, and Cohen's kappa. P

Many real-world signal sources are complex-valued, having real and imaginary components. However, the vast majority of existing deep learning platforms and network architectures do not support the use of complex-valued data. MRI data is inherently complex-valued, so existing approaches discard the richer algebraic structure of the complex data. In this work, we investigate end-to-end complex-valued convolutional neural networks - specifically, for image reconstruction in lieu of two-channel real-valued networks. We apply this to magnetic resonance imaging reconstruction for the purpose of accelerating scan times and determine the performance of various promising complex-valued activation functions. We find that complex-valued CNNs with complex-valued convolutions provide superior reconstructions compared to real-valued convolutions with the same number of trainable parameters, over a variety of network architectures and datasets.

Purpose: To develop a framework to reconstruct large-scale volumetric dynamic MRI from rapid continuous and non-gated acquisitions, with applications to pulmonary and dynamic contrast enhanced (DCE) imaging. Theory and Methods: The problem considered here requires recovering hundred-gigabytes of dynamic volumetric image data from a few gigabytes of k-space data, acquired continuously over several minutes. This reconstruction is vastly under-determined, heavily stressing computing resources as well as memory management and storage. To overcome these challenges, we leverage intrinsic three dimensional (3D) trajectories, such as 3D radial and 3D cones, with ordering that incoherently cover time and k-space over the entire acquisition. We then propose two innovations: (1) A compressed representation using multi-scale low rank matrix factorization that constrains the reconstruction problem, and reduces its memory footprint. (2) Stochastic optimization to reduce computation, improve memory locality, and minimize communications between threads and processors. We demonstrate the feasibility of the proposed method on DCE imaging acquired with a golden-angle ordered 3D cones trajectory and pulmonary imaging acquired with a bit-reversed ordered 3D radial trajectory. We compare it with "soft-gated" dynamic reconstruction for DCE and respiratory resolved reconstruction for pulmonary imaging. Results: The proposed technique shows transient dynamics that are not seen in gating based methods. When applied to datasets with irregular, or non-repetitive motions, the proposed method displays sharper image features. Conclusion: We demonstrated a method that can reconstruct massive 3D dynamic image series in the extreme undersampling and extreme computation setting., Comment: Accepted to Magnetic Resonance in Medicine

Compressed sensing in MRI enables high subsampling factors while maintaining diagnostic image quality. This technique enables shortened scan durations and/or improved image resolution. Further, compressed sensing can increase the diagnostic information and value from each scan performed. Overall, compressed sensing has significant clinical impact in improving the diagnostic quality and patient experience for imaging exams. However, a number of challenges exist when moving compressed sensing from research to the clinic. These challenges include hand-crafted image priors, sensitive tuning parameters, and long reconstruction times. Data-driven learning provides a solution to address these challenges. As a result, compressed sensing can have greater clinical impact. In this tutorial, we will review the compressed sensing formulation and outline steps needed to transform this formulation to a deep learning framework. Supplementary open source code in python will be used to demonstrate this approach with open databases. Further, we will discuss considerations in applying data-driven compressed sensing in the clinical setting., Comment: Submitted to the Special Issue on Computational MRI: Compressed Sensing and Beyond in the IEEE Signal Processing Magazine

BACKGROUND:In patients with mitral or tricuspid valve regurgitation, evaluation of regurgitant severity is essential for determining the need for surgery. While transthoracic echocardiography is widely accessible, it has limited reproducibility for grading inlet valve regurgitation. Multiplanar cardiac MRI is the quantitative standard but requires specialized local expertise, and is thus not widely available. Volumetric 4D flow MRI has potential for quantitatively grading the severity of inlet valve regurgitation in adult patients. PURPOSE:To evaluate the accuracy and reproducibility of volumetric 4D flow MRI for quantification of inlet valvular regurgitation compared to conventional multiplanar MRI, which may simplify and improve accessibility of cardiac MRI. STUDY TYPE:This retrospective, HIPAA-compliant imaging-based comparison study was conducted at a single institution. SUBJECTS:Twenty-one patients who underwent concurrent multiplanar and 4D flow cardiac MRI between April 2015 and January 2017. FIELD STRENGTH/SEQUENCES:3T; steady-state free-precession (SSFP), 2D phase contrast (2D-PC), and postcontrast 4D flow. ASSESSMENT:We evaluated the intertechnique (4D flow vs. 2D-PC), intermethod (direct vs. indirect measurement), interobserver and intraobserver reproducibility of measurements of regurgitant flow volume (RFV), fraction (RF), and volume (RVol). STATISTICAL TESTS:Statistical analysis included Pearson correlation, Bland-Altman statistics, and intraclass correlation coefficients. RESULTS:There was high concordance between 4D flow and multiplanar MRI, whether using direct or indirect methods of quantifying regurgitation (r = 0.813-0.985). Direct interrogation of the regurgitant jet with 4D flow showed high intraobserver consistency (r = 0.976-0.999) and interobserver consistency (r = 0.861-0.992), and correlated well with traditional indirect measurements obtained as the difference between stroke volume and forward outlet valve flow. DATA CONCLUSION:4D flow MRI provides highly reproducible measurements of mitral and tricuspid regurgitant volume, and may be used in place of conventional multiplanar MRI. LEVEL OF EVIDENCE:4 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018;48:1147-1158.

Purpose: Off-resonance artifact correction by deep-learning, to facilitate rapid pediatric body imaging with a scan time efficient 3D cones trajectory. Methods: A residual convolutional neural network to correct off-resonance artifacts (Off-ResNet) was trained with a prospective study of 30 pediatric magnetic resonance angiography exams. Each exam acquired a short-readout scan (1.18 ms +- 0.38) and a long-readout scan (3.35 ms +- 0.74) at 3T. Short-readout scans, with longer scan times but negligible off-resonance blurring, were used as reference images and augmented with additional off-resonance for supervised training examples. Long-readout scans, with greater off-resonance artifacts but shorter scan time, were corrected by autofocus and Off-ResNet and compared to short-readout scans by normalized root-mean-square error (NRMSE), structural similarity index (SSIM), and peak signal-to-noise ratio (PSNR). Scans were also compared by scoring on eight anatomical features by two radiologists, using analysis of variance with post-hoc Tukey's test. Reader agreement was determined with intraclass correlation. Results: Long-readout scans were on average 59.3% shorter than short-readout scans. Images from Off-ResNet had superior NRMSE, SSIM, and PSNR compared to uncorrected images across +-1kHz off-resonance (P

PurposeA new acquisition and reconstruction method called T2 Shuffling is presented for volumetric fast spin-echo (three-dimensional [3D] FSE) imaging. T2 Shuffling reduces blurring and recovers many images at multiple T2 contrasts from a single acquisition at clinically feasible scan times (6-7 min).Theory and methodsThe parallel imaging forward model is modified to account for temporal signal relaxation during the echo train. Scan efficiency is improved by acquiring data during the transient signal decay and by increasing echo train lengths without loss in signal-to-noise ratio (SNR). By (1) randomly shuffling the phase encode view ordering, (2) constraining the temporal signal evolution to a low-dimensional subspace, and (3) promoting spatio-temporal correlations through locally low rank regularization, a time series of virtual echo time images is recovered from a single scan. A convex formulation is presented that is robust to partial voluming and radiofrequency field inhomogeneity.ResultsRetrospective undersampling and in vivo scans confirm the increase in sharpness afforded by T2 Shuffling. Multiple image contrasts are recovered and used to highlight pathology in pediatric patients. A proof-of-principle method is integrated into a clinical musculoskeletal imaging workflow.ConclusionThe proposed T2 Shuffling method improves the diagnostic utility of 3D FSE by reducing blurring and producing multiple image contrasts from a single scan. Magn Reson Med 77:180-195, 2017. © 2016 Wiley Periodicals, Inc.

We develop a novel next generation light-weight highly flexible pediatric coil array, combine it with a high-density pediatric posterior array or conventional posterior phased array, and determine feasibility of pediatric clinical use. A highly flexible 16 element MRI receiver coil was constructed with low-profile noise controlling preamplifiers that minimized reactive and resistive coupling. Element decoupling was assessed in flat and highly flexed states. With IRB approval and informed consent and assent, 24 consecutive subjects undergoing torso or extremity MRI were prospectively recruited. Care team members were surveyed on preference for the coil versus conventional coils and diagnostic acceptability of the images was recorded. Confidence interval of proportion of diagnostic exams was calculated. The array without cable weighed 480 grams, demonstrated good flexibility while maintaining element decoupling. The coil was preferred by all nurses and anesthesiologists involved in the care of the patients. Technologists preferred the coil in 96 percent of cases, and 23 of 24 exams were diagnostically adequate with 85 percent confidence interval of 90-100 percent. Light-weight highly flexible coil arrays can be constructed that maintain element decoupling. Pediatric clinical image quality is likely to be diagnostic, with acceptance by members of the care team., Comment: 22 pages, 7 figures, 2 tables

PurposeTo investigate four-dimensional flow denoising using the divergence-free wavelet (DFW) transform and compare its performance with existing techniques.Theory and methodsDFW is a vector-wavelet that provides a sparse representation of flow in a generally divergence-free field and can be used to enforce "soft" divergence-free conditions when discretization and partial voluming result in numerical nondivergence-free components. Efficient denoising is achieved by appropriate shrinkage of divergence-free wavelet and nondivergence-free coefficients. SureShrink and cycle spinning are investigated to further improve denoising performance.ResultsDFW denoising was compared with existing methods on simulated and phantom data and was shown to yield better noise reduction overall while being robust to segmentation errors. The processing was applied to in vivo data and was demonstrated to improve visualization while preserving quantifications of flow data.ConclusionDFW denoising of four-dimensional flow data was shown to reduce noise levels in flow data both quantitatively and visually.

PurposeTo determine whether it is feasible to use solely an accelerated 4D phase-contrast magnetic resonance imaging (4D-PC MRI) acquisition to quantify net and regurgitant flow volume through each of the cardiac valves.Materials and methodsAccelerated, 4D-PC MRI examinations performed between March 2010 through June 2011 as part of routine MRI examinations for congenital, structural heart disease were retrospectively reviewed and analyzed using valve-tracking visualization and quantification algorithms developed in Java and OpenGL. Excluding patients with transposition or single ventricle physiology, a total of 34 consecutive pediatric patients (19 male, 15 female; mean age 6.9 years; age range 10 months to 15 years) were identified. 4D-PC flow measurements were compared at each valve and against routine measurements from conventional cardiac MRI using Bland-Altman and Pearson correlation analysis.ResultsInlet and outlet valve net flow were highly correlated between all valves (P = 0.940-0.985). The sum of forward flow at the outlet valve and regurgitant flow at the inlet valve were consistent with volumetric displacements in each ventricle (P = 0.939-0.948). These were also highly consistent with conventional planar MRI measurements with net flow (P = 0.923-0.935) and regurgitant fractions (P = 0.917-0.972) at the outlet valve and ventricular volumes (P = 0.925-0.965).ConclusionIt is possible to obtain consistent measurements of net and regurgitant blood flow across the inlet and outlet valves relying solely on accelerated 4D-PC. This may facilitate more efficient clinical quantification of valvular regurgitation.