Dynamic magnetic resonance imaging method based on golden-ratio cartesian sampling and compressed sensing.
Ontology highlight
ABSTRACT:
SUBMITTER: Li S
PROVIDER: S-EPMC5790254 | biostudies-literature | 2018
REPOSITORIES: biostudies-literature
Publications
Dynamic magnetic resonance imaging method based on golden-ratio cartesian sampling and compressed sensing.
Li Shuo S Zhu Yanchun Y Xie Yaoqin Y Gao Song S
PloS one 20180130 1
Dynamic magnetic resonance imaging (DMRI) is used to noninvasively trace the movements of organs and the process of drug delivery. The results can provide quantitative or semiquantitative pathology-related parameters, thus giving DMRI great potential for clinical applications. However, conventional DMRI techniques suffer from low temporal resolution and long scan time owing to the limitations of the k-space sampling scheme and image reconstruction algorithm. In this paper, we propose a novel DMR ...[more]
Similar Datasets
Project description:To develop a fast and flexible free-breathing dynamic volumetric MRI technique, iterative Golden-angle RAdial Sparse Parallel MRI (iGRASP), that combines compressed sensing, parallel imaging, and golden-angle radial sampling.Radial k-space data are acquired continuously using the golden-angle scheme and sorted into time series by grouping an arbitrary number of consecutive spokes into temporal frames. An iterative reconstruction procedure is then performed on the undersampled time series where joint multicoil sparsity is enforced by applying a total-variation constraint along the temporal dimension. Required coil-sensitivity profiles are obtained from the time-averaged data.iGRASP achieved higher acceleration capability than either parallel imaging or coil-by-coil compressed sensing alone. It enabled dynamic volumetric imaging with high spatial and temporal resolution for various clinical applications, including free-breathing dynamic contrast-enhanced imaging in the abdomen of both adult and pediatric patients, and in the breast and neck of adult patients.The high performance and flexibility provided by iGRASP can improve clinical studies that require robustness to motion and simultaneous high spatial and temporal resolution. Magn Reson Med 72:707-717, 2014. © 2013 Wiley Periodicals, Inc.
| S-EPMC3991777 | biostudies-literature
Project description:Diffusion Magnetic Resonance Imaging (dMRI) has shown great potential in probing tissue microstructure and structural connectivity in the brain but is often limited by the lengthy scan time needed to sample the diffusion profile by acquiring multiple diffusion weighted images (DWIs). Although parallel imaging technique has improved the speed of dMRI acquisition, attaining high resolution three dimensional (3D) dMRI on preclinical MRI systems remained still time consuming. In this paper, kernel principal component analysis, a machine learning approach, was employed to estimate the correlation among DWIs. We demonstrated the feasibility of such correlation estimation from low-resolution training DWIs and used the correlation as a constraint to reconstruct high-resolution DWIs from highly under-sampled k-space data, which significantly reduced the scan time. Using full k-space 3D dMRI data of post-mortem mouse brains, we retrospectively compared the performance of the so-called kernel low rank (KLR) method with a conventional compressed sensing (CS) method in terms of image quality and ability to resolve complex fiber orientations and connectivity. The results demonstrated that the KLR-CS method outperformed the conventional CS method for acceleration factors up to 8 and was likely to enhance our ability to investigate brain microstructure and connectivity using high-resolution 3D dMRI.
| S-EPMC7359413 | biostudies-literature
Project description:PURPOSE:Conventional non-Cartesian compressed sensing requires multiple nonuniform Fourier transforms every iteration, which is computationally expensive. Accordingly, time-consuming reconstructions have slowed the adoption of undersampled 3D non-Cartesian acquisitions into clinical protocols. In this work we investigate several approaches to minimize reconstruction times without sacrificing accuracy. METHODS:The reconstruction problem can be reformatted to exploit the Toeplitz structure of matrices that are evaluated every iteration, but it requires larger oversampling than what is strictly required by nonuniform Fourier transforms. Accordingly, we investigate relative speeds of the two approaches for various nonuniform Fourier transform kernel sizes and oversampling for both GPU and CPU implementations. Second, we introduce a method to minimize matrix sizes by estimating the image support. Finally, density compensation weights have been used as a preconditioning matrix to improve convergence, but this increases noise. We propose a more general approach to preconditioning that allows a trade-off between accuracy and convergence speed. RESULTS:When using a GPU, the Toeplitz approach was faster for all practical parameters. Second, it was found that properly accounting for image support can prevent aliasing errors with minimal impact on reconstruction time. Third, the proposed preconditioning scheme improved convergence rates by an order of magnitude with negligible impact on noise. CONCLUSION:With the proposed methods, 3D non-Cartesian compressed sensing with clinically relevant reconstruction times (<2 min) is feasible using practical computer resources. Magn Reson Med 79:2685-2692, 2018. © 2017 International Society for Magnetic Resonance in Medicine.
| S-EPMC6755916 | biostudies-literature
Project description:BACKGROUND:4D flow cardiovascular magnetic resonance (CMR) enables visualization of complex blood flow and quantification of biomarkers for vessel wall disease, such as wall shear stress (WSS). Because of the inherently long acquisition times, many efforts have been made to accelerate 4D flow acquisitions, however, no detailed analysis has been made on the effect of Cartesian compressed sensing accelerated 4D flow CMR at different undersampling rates on quantitative flow parameters and WSS. METHODS:We implemented a retrospectively triggered 4D flow CMR acquisition with pseudo-spiral Cartesian k-space filling, which results in incoherent undersampling of k-t space. Additionally, this strategy leads to small jumps in k-space thereby minimizing eddy current related artifacts. The pseudo-spirals were rotated in a tiny golden-angle fashion, which provides optimal incoherence and a variable density sampling pattern with a fully sampled center. We evaluated this 4D flow protocol in a carotid flow phantom with accelerations of R = 2-20, as well as in carotids of 7 healthy subjects (27 ± 2 years, 4 male) for R = 10-30. Fully sampled 2D flow CMR served as a flow reference. Arteries were manually segmented and registered to enable voxel-wise comparisons of both velocity and WSS using a Bland-Altman analysis. RESULTS:Magnitude images, velocity images, and pathline reconstructions from phantom and in vivo scans were similar for all accelerations. For the phantom data, mean differences at peak systole for the entire vessel volume in comparison to R = 2 ranged from - 2.3 to - 5.3% (WSS) and - 2.4 to - 2.2% (velocity) for acceleration factors R = 4-20. For the in vivo data, mean differences for the entire vessel volume at peak systole in comparison to R = 10 were - 9.9, - 13.4, and - 16.9% (WSS) and - 8.4, - 10.8, and - 14.0% (velocity), for R = 20, 25, and 30, respectively. Compared to single slice 2D flow CMR acquisitions, peak systolic flow rates of the phantom showed no differences, whereas peak systolic flow rates in the carotid artery in vivo became increasingly underestimated with increasing acceleration. CONCLUSION:Acquisition of 4D flow CMR of the carotid arteries can be highly accelerated by pseudo-spiral k-space sampling and compressed sensing reconstruction, with consistent data quality facilitating velocity pathline reconstructions, as well as quantitative flow rate and WSS estimations. At an acceleration factor of R = 20 the underestimation of peak velocity and peak WSS was acceptable (< 10%) in comparison to an R = 10 accelerated 4D flow CMR reference scan. Peak flow rates were underestimated in comparison with 2D flow CMR and decreased systematically with higher acceleration factors.
| S-EPMC6971939 | biostudies-literature
Project description:BACKGROUND:Coronary magnetic resonance angiography (CMRA) is a promising technique for assessing the coronary arteries. However, a disadvantage of CMRA is the comparatively long acquisition time. Compressed sensing (CS) can considerably reduce the scan time. The aim of this study was to verify the feasibility of CS CMRA scanning during the waiting time between contrast injection and late gadolinium enhancement (LGE) scan in a clinical protocol. METHODS:Fifty clinical patients underwent contrast-enhanced CS CMRA and conventional CMRA on a 3?T CMR scanner. After contrast injection, CS CMRA was scanned during the waiting time for LGE CMR. A conventional CMRA scan was performed after LGE CMR. We assessed acquisition times and coronary artery image quality for each segment on a 4-point scale. Visible vessel length, sharpness and diameter of right (RCA), left anterior descending (LAD), and left circumflex (LCX) coronary arteries were also quantitatively compared among the scans. RESULTS:All CS CMRA scans were successfully performed within the LGE waiting time. The median total scan time was 207?s (163, 259?s) for CS and 785?s (698, 975?s) for conventional CMRA (p?<?0.001). No significant differences were observed in image quality scores, vessel length measurements, sharpness, and diameter between CS and conventional CMRA. CONCLUSIONS:We could achieve all CS CMRA scans within the LGE waiting time. Contrast-enhanced CS CMRA could considerably shorten the scan time while maintaining image quality compared with conventional CMRA.
| S-EPMC7017458 | biostudies-literature
Project description:BackgroundEndogenous contrast T1? cardiovascular magnetic resonance (CMR) can detect scar or infiltrative fibrosis in patients with ischemic or non-ischemic cardiomyopathy. Existing 2D T1? techniques have limited spatial coverage or require multiple breath-holds. The purpose of this project was to develop an accelerated, free-breathing 3D T1? mapping sequence with whole left ventricle coverage using a multicoil, compressed sensing (CS) reconstruction technique for rapid reconstruction of undersampled k-space data.MethodsWe developed a cardiac- and respiratory-gated, free-breathing 3D T1? sequence and acquired data using a variable-density k-space sampling pattern (A?=?3). The effect of the transient magnetization trajectory, incomplete recovery of magnetization between T1?-preparations (heart rate dependence), and k-space sampling pattern on T1? relaxation time error and edge blurring was analyzed using Bloch simulations for normal and chronically infarcted myocardium. Sequence accuracy and repeatability was evaluated using MnCl2 phantoms with different T1? relaxation times and compared to 2D measurements. We further assessed accuracy and repeatability in healthy subjects and compared these results to 2D breath-held measurements.ResultsThe error in T1? due to incomplete recovery of magnetization between T1?-preparations was T1?healthy =?6.1% and T1?infarct =?10.8% at 60?bpm and T1?healthy =?13.2% and T1?infarct =?19.6% at 90?bpm. At a heart rate of 60?bpm, error from the combined effects of readout-dependent magnetization transients, k-space undersampling and reordering was T1?healthy =?12.6% and T1?infarct =?5.8%. CS reconstructions had improved edge sharpness (blur metric?=?0.15) compared to inverse Fourier transform reconstructions (blur metric?=?0.48). There was strong agreement between the mean T1? estimated from the 2D and accelerated 3D data (R2 =?0.99; P <?0.05) acquired on the MnCl2 phantoms. The mean R1? estimated from the accelerated 3D sequence was highly correlated with MnCl2 concentration (R2 =?0.99; P?<?0.05). 3D T1? acquisitions were successful in all human subjects. There was no significant bias between undersampled 3D T1? and breath-held 2D T1? (mean bias?=?0.87) and the measurements had good repeatability (COV2D =?6.4% and COV3D =?7.1%).ConclusionsThis is the first report of an accelerated, free-breathing 3D T1? mapping of the left ventricle. This technique may improve non-contrast myocardial tissue characterization in patients with heart disease in a scan time appropriate for patients.
| S-EPMC6327532 | biostudies-literature
Project description:Compressed sensing has shown to be promising to accelerate magnetic resonance imaging. In this new technology, magnetic resonance images are usually reconstructed by enforcing its sparsity in sparse image reconstruction models, including both synthesis and analysis models. The synthesis model assumes that an image is a sparse combination of atom signals while the analysis model assumes that an image is sparse after the application of an analysis operator. Balanced model is a new sparse model that bridges analysis and synthesis models by introducing a penalty term on the distance of frame coefficients to the range of the analysis operator. In this paper, we study the performance of the balanced model in tight frame based compressed sensing magnetic resonance imaging and propose a new efficient numerical algorithm to solve the optimization problem. By tuning the balancing parameter, the new model achieves solutions of three models. It is found that the balanced model has a comparable performance with the analysis model. Besides, both of them achieve better results than the synthesis model no matter what value the balancing parameter is. Experiment shows that our proposed numerical algorithm constrained split augmented Lagrangian shrinkage algorithm for balanced model (C-SALSA-B) converges faster than previously proposed algorithms accelerated proximal algorithm (APG) and alternating directional method of multipliers for balanced model (ADMM-B).
| S-EPMC4388626 | biostudies-literature
Project description:Purpose To assess the effect of a combination of compressed sensing and SENSitivity Encoding (SENSE) acceleration techniques on radiation therapy magnetic resonance imaging (MRI) simulation workflows. Methods and Materials Thirty-seven acquisitions were performed with both SENSE-only (SENSE) and combined compressed sensing and SENSE (CS) techniques in 24 patients receiving radiation therapy MRI simulation for a wide range of disease sites. The anatomic field of view prescription and image resolution were identical for both SENSE and CS acquisitions to ensure fair comparison. The acquisition time of all images was recorded to assess time savings. For each image pair, image quality, and ability to contour were assessed by 2 radiation oncologists. Aside from direct image pair comparisons, the feasibility of using CS to improve MRI simulation protocols by increasing image resolution, field of view, and reducing motion artifacts was also evaluated. Results CS resulted in an average reduction of 27% in scan time with negligible changes in image quality and the ability to contour structures for RT treatment planning compared with SENSE. Physician scoring of image quality and ability to contour shows that while SENSE still has slightly better image quality compared with CS, this observed difference in image quality did not affect the ability to contour. In addition, the higher acceleration capability of CS enabled use of superior-inferior direction phase encoding in a sagittal 3-dimensional T2-weighted scan for substantially improved visibility of the prostatic urethra, which eliminated the need for a Foley catheter in most patients. Conclusions The combination of compressed sensing and parallel imaging resulted in marked improvements in the MRI Simulation workflow. The scan time was reduced without significantly affecting image quality in the context of ability to contour. The acceleration capabilities allowed for increased image resolution under similar scanning times as well as significantly improved urethra visualization in prostate simulations.
| S-EPMC8569477 | biostudies-literature
Project description:BACKGROUND:Three-dimensional time-resolved phase-contrast cardiovascular magnetic resonance (4D flow CMR) enables the quantification and visualisation of blood flow, but its clinical applicability remains hampered by its long scan time. The aim of this study was to evaluate the use of compressed sensing (CS) with on-line reconstruction to accelerate the acquisition and reconstruction of 4D flow CMR of the thoracic aorta. METHODS:4D flow CMR of the thoracic aorta was acquired in 20 healthy subjects using CS with acceleration factors ranging from 4 to 10. As a reference, conventional parallel imaging (SENSE) with acceleration factor 2 was used. Flow curves, net flows, peak flows and peak velocities were extracted from six contours along the aorta. To measure internal data consistency, a quantitative particle trace analysis was performed. Additionally, scan-rescan, inter- and intraobserver reproducibility were assessed. Subsequently, 4D flow CMR with CS factor 6 was acquired in 3 patients with differing aortopathies. The flow patterns resulting from particle trace visualisation were qualitatively analysed. RESULTS:All collected data were successfully acquired and reconstructed on-line. The average acquisition time including respiratory navigator efficiency with CS factor 6 was 5:02?±?2:23?min while reconstruction took approximately 9?min. For CS factors of 8 or less, mean differences in net flow, peak flow and peak velocity as compared to SENSE were below 2.2?±?7.8?ml/cycle, 4.6?±?25.2?ml/s and?-?7.9?±?13.0?cm/s, respectively. For a CS factor of 10 differences reached 5.4?±?8.0?ml/cycle, 14.4?±?28.3?ml/s and?-?4.0?±?12.2?cm/s. Scan-rescan analysis yielded mean differences in net flow of -?0.7?±?4.9?ml/cycle for SENSE and?-?0.2?±?8.5?ml/cycle for CS factor of 6. CONCLUSIONS:A six- to eightfold acceleration of 4D flow CMR using CS is feasible. Up to a CS acceleration rate of 6, no statistically significant differences in measured flow parameters could be observed with respect to the reference technique. Acquisitions in patients with aortopathies confirm the potential to integrate the proposed method in a clinical routine setting, whereby its main benefits are scan-time savings and direct on-line reconstruction.
| S-EPMC6802342 | biostudies-literature
Project description:PurposeDesign of a preconditioner for fast and efficient parallel imaging (PI) and compressed sensing (CS) reconstructions for Cartesian trajectories.TheoryPI and CS reconstructions become time consuming when the problem size or the number of coils is large, due to the large linear system of equations that has to be solved in ? 1 and ? 2 -norm based reconstruction algorithms. Such linear systems can be solved efficiently using effective preconditioning techniques.MethodsIn this article we construct such a preconditioner by approximating the system matrix of the linear system, which comprises the data fidelity and includes total variation and wavelet regularization, by a matrix that is block circulant with circulant blocks. Due to this structure, the preconditioner can be constructed quickly and its inverse can be evaluated fast using only two fast Fourier transformations. We test the performance of the preconditioner for the conjugate gradient method as the linear solver, integrated into the well-established Split Bregman algorithm.ResultsThe designed circulant preconditioner reduces the number of iterations required in the conjugate gradient method by almost a factor of 5. The speed up results in a total acceleration factor of approximately 2.5 for the entire reconstruction algorithm when implemented in MATLAB, while the initialization time of the preconditioner is negligible.ConclusionThe proposed preconditioner reduces the reconstruction time for PI and CS in a Split Bregman implementation without compromising reconstruction stability and can easily handle large systems since it is Fourier-based, allowing for efficient computations.
| S-EPMC6283050 | biostudies-literature