Accelerated MRI with CIRcular Cartesian UnderSampling (CIRCUS): a variable density Cartesian sampling strategy for compressed sensing and parallel imaging

A novel automatic segmentation method directly based on magnetic resonance imaging K-space data for auxiliary diagnosis of glioma

Background

The use of segmentation architectures in medical imaging, particularly for glioma diagnosis, marks a significant advancement in the field. Traditional methods often rely on post-processed images; however, key details can be lost during the fast Fourier transformation (FFT) process. Given the limitations of these techniques, there is a growing interest in exploring more direct approaches. The adaption of segmentation architectures originally designed for road extraction for medical imaging represents an innovative step in this direction. By employing K-space data as the modal input, this method completely eliminates the information loss inherent in FFT, thereby potentially enhancing the precision and effectiveness of glioma diagnosis.

Methods

In the study, a novel architecture based on a deep-residual U-net was developed to accomplish the challenging task of automatically segmenting brain tumors from K-space data. Brain tumors from K-space data with different under-sampling rates were also segmented to verify the clinical application of our method.

Results

Compared to the benchmarks set in the 2018 Brain Tumor Segmentation (BraTS) Challenge, our proposed architecture had superior performance, achieving Dice scores of 0.8573, 0.8789, and 0.7765 for the whole tumor (WT), tumor core (TC), and enhanced tumor (ET) regions, respectively. The corresponding Hausdorff distances were 2.5649, 1.6146, and 2.7187 for the WT, TC, and ET regions, respectively. Notably, compared to traditional image-based approaches, the architecture also exhibited an improvement of approximately 10% in segmentation accuracy on the K-space data at different under-sampling rates.

Conclusions

These results show the superiority of our method compared to previous methods. The direct performance of lesion segmentation based on K-space data eliminates the time-consuming and tedious image reconstruction process, thus enabling the segmentation task to be accomplished more efficiently.

Cardiac MR: From Theory to Practice

Cardiovascular disease (CVD) is the leading single cause of morbidity and mortality, causing over 17. 9 million deaths worldwide per year with associated costs of over $800 billion. Improving prevention, diagnosis, and treatment of CVD is therefore a global priority. Cardiovascular magnetic resonance (CMR) has emerged as a clinically important technique for the assessment of cardiovascular anatomy, function, perfusion, and viability. However, diversity and complexity of imaging, reconstruction and analysis methods pose some limitations to the widespread use of CMR. Especially in view of recent developments in the field of machine learning that provide novel solutions to address existing problems, it is necessary to bridge the gap between the clinical and scientific communities. This review covers five essential aspects of CMR to provide a comprehensive overview ranging from CVDs to CMR pulse sequence design, acquisition protocols, motion handling, image reconstruction and quantitative analysis of the obtained data. (1) The basic MR physics of CMR is introduced. Basic pulse sequence building blocks that are commonly used in CMR imaging are presented. Sequences containing these building blocks are formed for parametric mapping and functional imaging techniques. Commonly perceived artifacts and potential countermeasures are discussed for these methods. (2) CMR methods for identifying CVDs are illustrated. Basic anatomy and functional processes are described to understand the cardiac pathologies and how they can be captured by CMR imaging. (3) The planning and conduct of a complete CMR exam which is targeted for the respective pathology is shown. Building blocks are illustrated to create an efficient and patient-centered workflow. Further strategies to cope with challenging patients are discussed. (4) Imaging acceleration and reconstruction techniques are presented that enable acquisition of spatial, temporal, and parametric dynamics of the cardiac cycle. The handling of respiratory and cardiac motion strategies as well as their integration into the reconstruction processes is showcased. (5) Recent advances on deep learning-based reconstructions for this purpose are summarized. Furthermore, an overview of novel deep learning image segmentation and analysis methods is provided with a focus on automatic, fast and reliable extraction of biomarkers and parameters of clinical relevance.

Mapping of fatty acid composition with free-breathing MR spectroscopic imaging and compressed sensing

Nonalcoholic fatty liver disease (NAFLD) is a growing health problem, and a major challenge in NAFLD management is identifying which patients are at risk of progression to more serious disease. Simple measurements of liver fat content are not strong predictors of clinical outcome, but biomarkers related to fatty acid composition (ie, saturated vs. unsaturated fat) may be more effective. MR spectroscopic imaging (MRSI) methods allow spatially resolved, whole-liver measurements of chemical composition but are traditionally limited by slow acquisition times. In this work we present an accelerated MRSI acquisition based on spin echo single point imaging (SE-SPI), which, using appropriate sampling and compressed sensing reconstruction, allows free-breathing acquisition in a mouse model of fatty liver disease. After validating the technique's performance in oil/water phantoms, we imaged mice that had received a normal diet or a methionine and choline deficient (MCD) diet, some of which also received supplemental injections of iron to mimic hepatic iron overload. SE-SPI was more resistant to the line-broadening effects of iron than single-voxel spectroscopy measurements, and was consistently able to measure the amplitudes of low-intensity spectral peaks that are important to characterizing fatty acid composition. In particular, in the mice receiving the MCD diet, SE-SPI showed a significant decrease in a metric associated with unsaturated fat, which is consistent with the literature. This or other related metrics may therefore offer more a specific biomarker of liver health than fat content alone. This preclinical study is an important precursor to clinical testing of the proposed method. MR-based quantification of fatty acid composition may allow for improved characterization of non-alcoholic fatty liver disease. A spectroscopic imaging method with appropriate sampling strategy allows whole-liver mapping of fat composition metrics in a free-breathing mouse model. Changes in metrics like the surrogate unsaturation index (UIs) are visible in mice receiving a diet which induces fat accumulation in the liver, as compared to a normal diet; such metrics may prove useful in future clinical studies of liver disease.

Fetal Cerebral Oxygenation Is Impaired in Congenital Heart Disease and Shows Variable Response to Maternal Hyperoxia

Background

Impairments in fetal oxygen delivery have been implicated in brain dysmaturation seen in congenital heart disease (CHD), suggesting a role for in utero transplacental oxygen therapy. We applied a novel imaging tool to quantify fetal cerebral oxygenation by measuring T2* decay. We compared T2* in fetuses with CHD with controls with a focus on cardiovascular physiologies (transposition or left‐sided obstruction) and described the effect of brief administration of maternal hyperoxia on T2* decay.

Methods and Results

This is a prospective study performed on pregnant mothers with a prenatal diagnosis of CHD compared with controls in the third trimester. Participants underwent a fetal brain magnetic resonance imaging scan including a T2* sequence before and after maternal hyperoxia. Comparisons were made between control and CHD fetuses including subgroup analyses by cardiac physiology. Forty‐four mothers (CHD=24, control=20) participated. Fetuses with CHD had lower total brain volume (238.2 mm³, 95% CI, 224.6–251.9) compared with controls (262.4 mm³, 95% CI, 245.0–279.8, P=0.04). T2* decay time was faster in CHD compared with controls (beta=−14.4, 95% CI, −23.3 to −5.6, P=0.002). The magnitude of change in T2* with maternal hyperoxia was higher in fetuses with transposition compared with controls (increase of 8.4 ms, 95% CI, 0.5–14.3, P=0.01), though between‐subject variability was noted.

Conclusions

Cerebral tissue oxygenation is lower in fetuses with complex CHD. There was variability in the response to maternal hyperoxia by CHD subgroup that can be tested in future larger studies. Cardiovascular physiology is critical when designing neuroprotective clinical trials in the fetus with CHD.

Evaluation of Golden-Angle-Sampled Dynamic Contrast-Enhanced MRI Reconstruction Using Objective Image Quality Measures: A Simulated Phantom Study

We aim to extend the use of image quality metrics (IQMs) from static magnetic resonance imaging (MRI) applications to dynamic MRI studies. We assessed the use of 2 IQMs, the root mean square error and structural similarity index, in evaluating the reconstruction of quantitative dynamic contrast-enhanced (DCE) MRI data acquired using golden-angle sampling and compressed sensing (CS). To address the difficulty of obtaining ground-truth knowledge of parameters describing dynamics in real patient data, we developed a Matlab simulation framework to assess quantitative CS-DCE-MRI. We began by validating the response of each IQM to the CS-MRI reconstruction process using static data and the performance of our simulation framework with simple dynamic data. We then extended the simulations to the more realistic extended Tofts model. When assessing the Tofts model, we tested 4 different methods of selecting a reference image for the IQMs. Results from the retrospective static CS-MRI reconstructions showed that each IQM is responsive to the CS-MRI reconstruction process. Simulations of a simple contrast evolution model validated the performance of our framework. Despite the complexity of the Tofts model, both IQM scores correlated well with the recovery accuracy of a central model parameter for all reference cases studied. This finding may form the basis of algorithms for automated selection of image reconstruction aspects, such as temporal resolution, in golden-angle-sampled CS-DCE-MRI. These further suggest that objective measures of image quality may find use in general dynamic MRI applications.

Non-contrast coronary magnetic resonance angiography: current frontiers and future horizons

Coronary magnetic resonance angiography (coronary MRA) is advantageous in its ability to assess coronary artery morphology and function without ionizing radiation or contrast media. However, technical limitations including reduced spatial resolution, long acquisition times, and low signal-to-noise ratios prevent it from clinical routine utilization. Nonetheless, each of these limitations can be specifically addressed by a combination of novel technologies including super-resolution imaging, compressed sensing, and deep-learning reconstruction. In this paper, we first review the current clinical use and motivations for non-contrast coronary MRA, discuss currently available coronary MRA techniques, and highlight current technical developments that hold unique potential to optimize coronary MRA image acquisition and post-processing. In the final section, we examine the various research-based coronary MRA methods and metrics that can be leveraged to assess coronary stenosis severity, physiological function, and atherosclerotic plaque characterization. We specifically discuss how such technologies may contribute to the clinical translation of coronary MRA into a robust modality for routine clinical use.

Reconstruction of undersampled 3D non-Cartesian image-based navigators for coronary MRA using an unrolled deep learning model

PURPOSE

To rapidly reconstruct undersampled 3D non-Cartesian image-based navigators (iNAVs) using an unrolled deep learning (DL) model, enabling nonrigid motion correction in coronary magnetic resonance angiography (CMRA).

METHODS

An end-to-end unrolled network is trained to reconstruct beat-to-beat 3D iNAVs acquired during a CMRA sequence. The unrolled model incorporates a nonuniform FFT operator in TensorFlow to perform the data-consistency operation, and the regularization term is learned by a convolutional neural network (CNN) based on the proximal gradient descent algorithm. The training set includes 6,000 3D iNAVs acquired from 7 different subjects and 11 scans using a variable-density (VD) cones trajectory. For testing, 3D iNAVs from 4 additional subjects are reconstructed using the unrolled model. To validate reconstruction accuracy, global and localized motion estimates from DL model-based 3D iNAVs are compared with those extracted from 3D iNAVs reconstructed with l 1 -ESPIRiT. Then, the high-resolution coronary MRA images motion corrected with autofocusing using the l 1 -ESPIRiT and DL model-based 3D iNAVs are assessed for differences.

RESULTS

3D iNAVs reconstructed using the DL model-based approach and conventional l 1 -ESPIRiT generate similar global and localized motion estimates and provide equivalent coronary image quality. Reconstruction with the unrolled network completes in a fraction of the time compared to CPU and GPU implementations of l 1 -ESPIRiT (20× and 3× speed increases, respectively).

CONCLUSIONS

We have developed a deep neural network architecture to reconstruct undersampled 3D non-Cartesian VD cones iNAVs. Our approach decreases reconstruction time for 3D iNAVs, while preserving the accuracy of nonrigid motion information offered by them for correction.

Incorporating prior knowledge via volumetric deep residual network to optimize the reconstruction of sparsely sampled MRI

For sparse sampling that accelerates magnetic resonance (MR) image acquisition, non-linear reconstruction algorithms have been developed, which incorporated patient specific a prior information. More generic a prior information could be acquired via deep learning and utilized for image reconstruction. In this study, we developed a volumetric hierarchical deep residual convolutional neural network, referred to as T-Net, to provide a data-driven end-to-end mapping from sparsely sampled MR images to fully sampled MR images, where cartilage MR images were acquired using an Ultra-short TE sequence and retrospectively undersampled using pseudo-random Cartesian and radial acquisition schemes. The network had a hierarchical architecture that promoted the sparsity of feature maps and increased the receptive field, which were valuable for signal synthesis and artifact suppression. Relatively dense local connections and global shortcuts were established to facilitate residual learning and compensate for details lost in hierarchical processing. Additionally, volumetric processing was adopted to fully exploit spatial continuity in three-dimensional space. Data consistency was further enforced. The network was trained with 336 three-dimensional images (each consisting of 32 slices) and tested by 24 images. The incorporation of a priori information acquired via deep learning facilitated high acceleration factors (as high as 8) while maintaining high image fidelity (quantitatively evaluated using the structural similarity index measurement). The proposed T-Net had an improved performance as compared to several state-of-the-art networks.

Highly accelerated 4D flow cardiovascular magnetic resonance using a pseudo-spiral Cartesian acquisition and compressed sensing reconstruction for carotid flow and wall shear stress

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.

Technical Note: Simultaneous segmentation and relaxometry for MRI through multitask learning

PURPOSE

This study demonstrated a magnetic resonance (MR) signal multitask learning method for three-dimensional (3D) simultaneous segmentation and relaxometry of human brain tissues.

MATERIALS AND METHODS

A 3D inversion-prepared balanced steady-state free precession sequence was used for acquiring in vivo multicontrast brain images. The deep neural network contained three residual blocks, and each block had 8 fully connected layers with sigmoid activation, layer norm, and 256 neurons in each layer. Online-synthesized MR signal evolutions and labels were used to train the neural network batch-by-batch. Empirically defined ranges of T1 and T2 values for the normal gray matter, white matter, and cerebrospinal fluid (CSF) were used as the prior knowledge. MRI brain experiments were performed on three healthy volunteers. The mean and standard deviation for the T1 and T2 values in vivo were reported and compared to literature values. Additional animal (N = 6) and prostate patient (N = 1) experiments were performed to compare the estimated T1 and T2 values with those from gold standard methods and to demonstrate clinical applications of the proposed method.

RESULTS

In animal validation experiment, the differences/errors (mean difference ± standard deviation of difference) between the T1 and T2 values estimated from the proposed method and the ground truth were 113 ± 486 and 154 ± 512 ms for T1, and 5 ± 33 and 7 ± 41 ms for T2, respectively. In healthy volunteer experiments (N = 3), whole brain segmentation and relaxometry were finished within ~ 5 s. The estimated apparent T1 and T2 maps were in accordance with known brain anatomy, and not affected by coil sensitivity variation. Gray matter, white matter, and CSF were successfully segmented. The deep neural network can also generate synthetic T1- and T2-weighted images.

CONCLUSION

The proposed multitask learning method can directly generate brain apparent T1 and T2 maps, as well as synthetic T1- and T2-weighted images, in conjunction with segmentation of gray matter, white matter, and CSF.

Extending Cardiac Functional Assessment with Respiratory-Resolved 3D Cine MRI

This study aimed to develop a cardiorespiratory-resolved 3D magnetic resonance imaging (5D MRI: x-y-z-cardiac-respiratory) approach based on 3D motion tracking for investigating the influence of respiration on cardiac ventricular function. A highly-accelerated 2.5-minute sparse MR protocol was developed for a continuous acquisition of cardiac images through multiple cardiac and respiratory cycles. The heart displacement along respiration was extracted using a 3D image deformation algorithm, and this information was used to cluster the acquired data into multiple respiratory phases. The proposed approach was tested in 15 healthy volunteers (7 females). Cardiac function parameters, including the end-systolic volume (ESV), end-diastolic volume (EDV), stroke volume (SV), and ejection fraction (EF), were measured for the left and right ventricle in both end-expiration and end-inspiration. Although with the proposed 5D cardiac MRI, there were no significant differences (p > 0.05, t-test) between end-expiration and end-inspiration measurements of the cardiac function in volunteers, incremental respiratory motion parameters that were derived from 3D motion tracking, such as the depth, expiration and inspiration distribution, correlated (p < 0.05, correlation coefficient, Mann-Whitney) with those volume-based parameters of cardiac function and varied between genders. The obtained initial results suggested that this new approach allows evaluation of cardiac function during specific respiratory phases. Thus, it can enable investigation of effects related to respiratory variability and better assessment of cardiac function for studying respiratory and/or cardiac dysfunction.

Cortical bone vessel identification and quantification on contrast-enhanced MR images

Background

Cortical bone porosity is a major determinant of bone strength. Despite the biomechanical importance of cortical bone porosity, the biological drivers of cortical porosity are unknown. The content of cortical pore space can indicate pore expansion mechanisms; both of the primary components of pore space, vessels and adipocytes, have been implicated in pore expansion. Dynamic contrast-enhanced MRI (DCE-MRI) is widely used in vessel detection in cardiovascular studies, but has not been applied to visualize vessels within cortical bone. In this study, we have developed a multimodal DCE-MRI and high resolution peripheral QCT (HR-pQCT) acquisition and image processing pipeline to detect vessel-filled cortical bone pores.

Methods

For this in vivo human study, 19 volunteers (10 males and 9 females; mean age =63±5) were recruited. Both distal and ultra-distal regions of the non-dominant tibia were imaged by HR-pQCT (82 µm nominal resolution) for bone structure segmentation and by 3T DCE-MRI (Gadavist; 9 min scan time; temporal resolution =30 sec; voxel size 230×230×500 µm³) for vessel visualization. The DCE-MRI was registered to the HR-pQCT volume and the voxels within the MRI cortical bone region were extracted. Features of the DCE data were calculated and voxels were categorized by a 2-stage hierarchical kmeans clustering algorithm to determine which voxels represent vessels. Vessel volume fraction (volume ratio of vessels to cortical bone), vessel density (average vessel count per cortical bone volume), and average vessel volume (mean volume of vessels) were calculated to quantify the status of vessel-filled pores in cortical bone. To examine spatial resolution and perform validation, a virtual phantom with 5 channel sizes and an applied pseudo enhancement curve was processed through the proposed image processing pipeline. Overlap volume ratio and Dice coefficient was calculated to measure the similarity between the detected vessel map and ground truth.

Results

In the human study, mean vessel volume fraction was 2.2%±1.0%, mean vessel density was 0.68±0.27 vessel/mm³, and mean average vessel volume was 0.032±0.012 mm³/vessel. Signal intensity for detected vessel voxels increased during the scan, while signal for non-vessel voxels within pores did not enhance. In the validation phantom, channels with diameter 250 µm or greater were detected successfully, with volume ratio equal to 1 and Dice coefficient above 0.6. Both statistics decreased dramatically for channel sizes less than 250 µm.

Conclusions

We have a developed a multi-modal image acquisition and processing pipeline that successfully detects vessels within cortical bone pores. The performance of this technique degrades for vessel diameters below the in-plane spatial resolution of the DCE-MRI acquisition. This approach can be applied to investigate the biological systems associated with cortical pore expansion.

Self-Attention Convolutional Neural Network for Improved MR Image Reconstruction

MRI is an advanced imaging modality with the unfortunate disadvantage of long data acquisition time. To accelerate MR image acquisition while maintaining high image quality, extensive investigations have been conducted on image reconstruction of sparsely sampled MRI. Recently, deep convolutional neural networks have achieved promising results, yet the local receptive field in convolution neural network raises concerns regarding signal synthesis and artifact compensation. In this study, we proposed a deep learning-based reconstruction framework to provide improved image fidelity for accelerated MRI. We integrated the self-attention mechanism, which captured long-range dependencies across image regions, into a volumetric hierarchical deep residual convolutional neural network. Basically, a self-attention module was integrated to every convolutional layer, where signal at a position was calculated as a weighted sum of the features at all positions. Furthermore, relatively dense shortcut connections were employed, and data consistency was enforced. The proposed network, referred to as SAT-Net, was applied on cartilage MRI acquired using an ultrashort TE sequence and retrospectively undersampled in a pseudo-random Cartesian pattern. The network was trained using 336 three dimensional images (each containing 32 slices) and tested with 24 images that yielded improved outcome. The framework is generic and can be extended to various applications.

Multiresolution spatiotemporal mechanical model of the heart as a prior to constrain the solution for 4D models of the heart

In several nuclear cardiac imaging applications (SPECT and PET), images are formed by reconstructing tomographic data using an iterative reconstruction algorithm with corrections for physical factors involved in the imaging detection process and with corrections for cardiac and respiratory motion. The physical factors are modeled as coefficients in the matrix of a system of linear equations and include attenuation, scatter, and spatially varying geometric response. The solution to the tomographic problem involves solving the inverse of this system matrix. This requires the design of an iterative reconstruction algorithm with a statistical model that best fits the data acquisition. The most appropriate model is based on a Poisson distribution. Using Bayes Theorem, an iterative reconstruction algorithm is designed to determine the maximum a posteriori estimate of the reconstructed image with constraints that maximizes the Bayesian likelihood function for the Poisson statistical model. The a priori distribution is formulated as the joint entropy (JE) to measure the similarity between the gated cardiac PET image and the cardiac MRI cine image modeled as a FE mechanical model. The developed algorithm shows the potential of using a FE mechanical model of the heart derived from a cardiac MRI cine scan to constrain solutions of gated cardiac PET images.

Developing an efficient phase-matched attenuation correction method for quiescent period PET in abdominal PET/MRI

Respiratory motion causes misalignments between positron emission tomography (PET) and magnetic resonance (MR)-derived attenuation maps (µ-maps) in addition to artifacts on both PET and MR images in simultaneous PET/MRI for organs such as liver that can experience motion of several centimeters. To address this problem, we developed an efficient MR-based attenuation correction (MRAC) method to generate phase-matched µ-maps for quiescent period PET (PET_Q) in abdominal PET/MRI. MRAC data was acquired with CIRcular Cartesian UnderSampling (CIRCUS) sampling during 100 s in free-breathing as an accelerated data acquisition strategy for phase-matched MRAC (MRAC_PM-CIRCUS). For comparison, MRAC data with raster (Default) k-space sampling was also acquired during 100 s in free-breathing (MRAC_PM-Default), and used to evaluate MRAC_PM-CIRCUS as well as un-matched MRAC (MRAC_UM) that was un-gated. We purposefully oversampled the MRAC_PM data to ensure we had enough information to capture all respiratory phases to make this comparison as robust as possible. The proposed MRAC_PM-CIRCUS was evaluated in 17 patients with ⁶⁸Ga-DOTA-TOC PET/MRI exams, suspected of having neuroendocrine tumors or liver metastases. Effects of CIRCUS sampling for accelerating a data acquisition were evaluated by simulating the data acquisition time retrospectively in increments of 5 s. Effects of MRAC_PM-CIRCUS on PET_Q were evaluated using uptake differences in the liver lesions (n  =  35), compared to PET_Q with MRAC_PM-Default and MRAC_UM. A Wilcoxon signed-rank test was performed to compare lesion uptakes between the MRAC methods. MRAC_PM-CIRCUS showed higher image quality compared to MRAC_PM-Default for the same acquisition times, demonstrating that a data acquisition time of 30 s was reasonable to achieve phase-matched µ-maps. Lesion update differences between MRAC_PM-CIRCUS (30 s) versus MRAC_PM-Default (reference, 100 s) were 0.1%  ±  1.4% (range of  -2.7% to 3.2%) and not significant (P  >  .05); while, the differences between MRAC_UM versus MRAC_PM-Default were 0.6%  ±  11.4% with a large variation (range of  -37% to 20%) and significant (P  <  .05). In conclusion, we demonstrated that a data acquisition of 30 s achieved phase-matched µ-maps when using specialized CIRCUS data sampling and phase-matched µ-maps improved PET_Q quantification significantly.

Highly accelerated intracranial 4D flow MRI: evaluation of healthy volunteers and patients with intracranial aneurysms

OBJECTIVES

To evaluate an accelerated 4D flow MRI method that provides high temporal resolution in a clinically feasible acquisition time for intracranial velocity imaging.

MATERIALS AND METHODS

Accelerated 4D flow MRI was developed by using a pseudo-random variable-density Cartesian undersampling strategy (CIRCUS) with the combination of k-t, parallel imaging and compressed sensing image reconstruction techniques (k-t SPARSE-SENSE). Four-dimensional flow data were acquired on five healthy volunteers and eight patients with intracranial aneurysms using CIRCUS (acceleration factor of R = 4, termed CIRCUS4) and GRAPPA (R = 2, termed GRAPPA2) as the reference method. Images with three times higher temporal resolution (R = 12, CIRCUS12) were also reconstructed from the same acquisition as CIRCUS4. Qualitative and quantitative image assessment was performed on the images acquired with different methods, and complex flow patterns in the aneurysms were identified and compared.

RESULTS

Four-dimensional flow MRI with CIRCUS was achieved in 5 min and allowed further improved temporal resolution of <30 ms. Volunteer studies showed similar qualitative and quantitative evaluation obtained with the proposed approach compared to the reference (overall image scores: GRAPPA2 3.2 ± 0.6; CIRCUS4 3.1 ± 0.7; CIRCUS12 3.3 ± 0.4; difference of the peak velocities: -3.83 ± 7.72 cm/s between CIRCUS4 and GRAPPA2, -1.72 ± 8.41 cm/s between CIRCUS12 and GRAPPA2). In patients with intracranial aneurysms, the higher temporal resolution improved capturing of the flow features in intracranial aneurysms (pathline visualization scores: GRAPPA2 2.2 ± 0.2; CIRCUS4 2.5 ± 0.5; CIRCUS12 2.7 ± 0.6).

CONCLUSION

The proposed rapid 4D flow MRI with a high temporal resolution is a promising tool for evaluating intracranial aneurysms in a clinically feasible acquisition time.

5D whole-heart sparse MRI

PURPOSE

A 5D whole-heart sparse imaging framework is proposed for simultaneous assessment of myocardial function and high-resolution cardiac and respiratory motion-resolved whole-heart anatomy in a single continuous noncontrast MR scan.

METHODS

A non-electrocardiograph (ECG)-triggered 3D golden-angle radial balanced steady-state free precession sequence was used for data acquisition. The acquired 3D k-space data were sorted into a 5D dataset containing separated cardiac and respiratory dimensions using a self-extracted respiratory motion signal and a recorded ECG signal. Images were then reconstructed using XD-GRASP, a multidimensional compressed sensing technique exploiting correlations/sparsity along cardiac and respiratory dimensions. 5D whole-heart imaging was compared with respiratory motion-corrected 3D and 4D whole-heart imaging in nine volunteers for evaluation of the myocardium, great vessels, and coronary arteries. It was also compared with breath-held, ECG-gated 2D cardiac cine imaging for validation of cardiac function quantification.

RESULTS

5D whole-heart images received systematic higher quality scores in the myocardium, great vessels and coronary arteries. Quantitative coronary sharpness and length were always better for the 5D images. Good agreement was obtained for quantification of cardiac function compared with 2D cine imaging.

CONCLUSION

5D whole-heart sparse imaging represents a robust and promising framework for simplified comprehensive cardiac MRI without the need for breath-hold and motion correction. Magn Reson Med 79:826-838, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Highly-accelerated self-gated free-breathing 3D cardiac cine MRI: validation in assessment of left ventricular function

OBJECTIVE

This work presents a highly-accelerated, self-gated, free-breathing 3D cardiac cine MRI method for cardiac function assessment.

MATERIALS AND METHODS

A golden-ratio profile based variable-density, pseudo-random, Cartesian undersampling scheme was implemented for continuous 3D data acquisition. Respiratory self-gating was achieved by deriving motion signal from the acquired MRI data. A multi-coil compressed sensing technique was employed to reconstruct 4D images (3D+time). 3D cardiac cine imaging with self-gating was compared to bellows gating and the clinical standard breath-held 2D cine imaging for evaluation of self-gating accuracy, image quality, and cardiac function in eight volunteers. Reproducibility of 3D imaging was assessed.

RESULTS

Self-gated 3D imaging provided an image quality score of 3.4 ± 0.7 vs 4.0 ± 0 with the 2D method (p = 0.06). It determined left ventricular end-systolic volume as 42.4 ± 11.5 mL, end-diastolic volume as 111.1 ± 24.7 mL, and ejection fraction as 62.0 ± 3.1%, which were comparable to the 2D method, with bias ± 1.96 × SD of -0.8 ± 7.5 mL (p = 0.90), 2.6 ± 3.3 mL (p = 0.84) and 1.4 ± 6.4% (p = 0.45), respectively.

CONCLUSION

The proposed 3D cardiac cine imaging method enables reliable respiratory self-gating performance with good reproducibility, and provides comparable image quality and functional measurements to 2D imaging, suggesting that self-gated, free-breathing 3D cardiac cine MRI framework is promising for improved patient comfort and cardiac MRI scan efficiency.

Free breathing whole-heart 3D CINE MRI with self-gated Cartesian trajectory

Purpose

To present a method that uses a novel free-running self-gated acquisition to achieve isotropic resolution in whole heart 3D Cartesian cardiac CINE MRI.

Material and methods

3D cardiac CINE MRI using navigator gating results in long acquisition times. Recently, several frameworks based on self-gated non-Cartesian trajectories have been proposed to accelerate this acquisition. However, non-Cartesian reconstructions are computationally expensive due to gridding, particularly in 3D. In this work, we propose a novel highly efficient self-gated Cartesian approach for 3D cardiac CINE MRI. Acquisition is performed using CArtesian trajectory with Spiral PRofile ordering and Tiny golden angle step for eddy current reduction (so called here CASPR-Tiger). Data is acquired continuously under free breathing (retrospective ECG gating, no preparation pulses interruption) for 4–5 min and 4D whole-heart volumes (3D + cardiac phases) with isotropic spatial resolution are reconstructed from all available data using a soft gating technique combined with temporal total variation (TV) constrained iterative SENSE reconstruction.

Results

For data acquired on eight healthy subjects and three patients, the reconstructed images using the proposed method had good contrast and spatio-temporal variations, correctly recovering diastolic and systolic cardiac phases. Non-significant differences (P > 0.05) were observed in cardiac functional measurements obtained with proposed 3D approach and gold standard 2D multi-slice breath-hold acquisition.

Conclusion

The proposed approach enables isotropic 3D whole heart Cartesian cardiac CINE MRI in 4 to 5 min free breathing acquisition.

•

A novel self-gated 3D Cartesian acquisition is proposed for free breathing whole-heart cardiac MRI

•

The proposed framework has efficient k-space sampling, better eddy current performance and high computational efficiency

•

The Proposed method is able to achieve high spatio-temporal resolution 3D cardiac CINE

•

The proposed method only requires four to five minute free breathing scan

Compressed sensing for body MRI

The introduction of compressed sensing for increasing imaging speed in magnetic resonance imaging (MRI) has raised significant interest among researchers and clinicians, and has initiated a large body of research across multiple clinical applications over the last decade. Compressed sensing aims to reconstruct unaliased images from fewer measurements than are traditionally required in MRI by exploiting image compressibility or sparsity. Moreover, appropriate combinations of compressed sensing with previously introduced fast imaging approaches, such as parallel imaging, have demonstrated further improved performance. The advent of compressed sensing marks the prelude to a new era of rapid MRI, where the focus of data acquisition has changed from sampling based on the nominal number of voxels and/or frames to sampling based on the desired information content. This article presents a brief overview of the application of compressed sensing techniques in body MRI, where imaging speed is crucial due to the presence of respiratory motion along with stringent constraints on spatial and temporal resolution. The first section provides an overview of the basic compressed sensing methodology, including the notion of sparsity, incoherence, and nonlinear reconstruction. The second section reviews state-of-the-art compressed sensing techniques that have been demonstrated for various clinical body MRI applications. In the final section, the article discusses current challenges and future opportunities.

LEVEL OF EVIDENCE

5 J. Magn. Reson. Imaging 2017;45:966-987.

Sparse Reconstruction Techniques in Magnetic Resonance Imaging: Methods, Applications, and Challenges to Clinical Adoption

The family of sparse reconstruction techniques, including the recently introduced compressed sensing framework, has been extensively explored to reduce scan times in magnetic resonance imaging (MRI). While there are many different methods that fall under the general umbrella of sparse reconstructions, they all rely on the idea that a priori information about the sparsity of MR images can be used to reconstruct full images from undersampled data. This review describes the basic ideas behind sparse reconstruction techniques, how they could be applied to improve MRI, and the open challenges to their general adoption in a clinical setting. The fundamental principles underlying different classes of sparse reconstructions techniques are examined, and the requirements that each make on the undersampled data outlined. Applications that could potentially benefit from the accelerations that sparse reconstructions could provide are described, and clinical studies using sparse reconstructions reviewed. Lastly, technical and clinical challenges to widespread implementation of sparse reconstruction techniques, including optimization, reconstruction times, artifact appearance, and comparison with current gold standards, are discussed.

High-temporospatial-resolution dynamic contrast-enhanced (DCE) wrist MRI with variable-density pseudo-random circular Cartesian undersampling (CIRCUS) acquisition: evaluation of perfusion in rheumatoid arthritis patients

This study is to evaluate highly accelerated three-dimensional (3D) dynamic contrast-enhanced (DCE) wrist MRI for assessment of perfusion in rheumatoid arthritis (RA) patients. A pseudo-random variable-density undersampling strategy, circular Cartesian undersampling (CIRCUS), was combined with k-t SPARSE-SENSE reconstruction to achieve a highly accelerated 3D DCE wrist MRI. Two healthy volunteers and 10 RA patients were studied. Two patients were on methotrexate (MTX) only (Group I) and the other eight were treated with a combination therapy of MTX and anti-tumor necrosis factor (TNF) therapy (Group II). Patients were scanned at baseline and 3 month follow-up. DCE MR images were used to evaluate perfusion in synovitis and bone marrow edema pattern in the RA wrist joints. A series of perfusion parameters was derived and compared with clinical disease activity scores of 28 joints (DAS28). 3D DCE wrist MR images were obtained with a spatial resolution of 0.3 × 0.3 × 1.5 mm(3) and temporal resolution of 5 s (with an acceleration factor of 20). The derived perfusion parameters, most notably transition time (dT) of synovitis, showed significant negative correlations with DAS28-ESR (r = -0.80, p < 0.05) and DAS28-CRP (r = -0.87, p < 0.05) at baseline and also correlated significantly with treatment responses evaluated by clinical score changes between baseline and 3 month follow-up (with DAS28-ESR r = -0.79, p < 0.05, and DAS28-CRP r = -0.82, p < 0.05). Highly accelerated 3D DCE wrist MRI with improved temporospatial resolution has been achieved in RA patients and provides accurate assessment of neovascularization and perfusion in RA joints, showing promise as a potential tool for evaluating treatment responses.