Golden-angle radial sparse parallel MRI: combination of compressed sensing, parallel imaging, and golden-angle radial sampling for fast and flexible dynamic volumetric MRI

3D printing from MRI Data: Harnessing strengths and minimizing weaknesses

3D printing facilitates the creation of accurate physical models of patient-specific anatomy from medical imaging datasets. While the majority of models to date are created from computed tomography (CT) data, there is increasing interest in creating models from other datasets, such as ultrasound and magnetic resonance imaging (MRI). MRI, in particular, holds great potential for 3D printing, given its excellent tissue characterization and lack of ionizing radiation. There are, however, challenges to 3D printing from MRI data as well. Here we review the basics of 3D printing, explore the current strengths and weaknesses of printing from MRI data as they pertain to model accuracy, and discuss considerations in the design of MRI sequences for 3D printing. Finally, we explore the future of 3D printing and MRI, including creative applications and new materials.

LEVEL OF EVIDENCE

5 J. Magn. Reson. Imaging 2017;45:635-645.

Direct estimation of tracer-kinetic parameter maps from highly undersampled brain dynamic contrast enhanced MRI

PURPOSE

The purpose of this work was to develop and evaluate a T₁ -weighted dynamic contrast enhanced (DCE) MRI methodology where tracer-kinetic (TK) parameter maps are directly estimated from undersampled (k,t)-space data.

THEORY AND METHODS

The proposed reconstruction involves solving a nonlinear least squares optimization problem that includes explicit use of a full forward model to convert parameter maps to (k,t)-space, utilizing the Patlak TK model. The proposed scheme is compared against an indirect method that creates intermediate images by parallel imaging and compressed sensing before to TK modeling. Thirteen fully sampled brain tumor DCE-MRI scans with 5-second temporal resolution are retrospectively undersampled at rates R = 20, 40, 60, 80, and 100 for each dynamic frame. TK maps are quantitatively compared based on root mean-squared-error (rMSE) and Bland-Altman analysis. The approach is also applied to four prospectively R = 30 undersampled whole-brain DCE-MRI data sets.

RESULTS

In the retrospective study, the proposed method performed statistically better than indirect method at R ≥ 80 for all 13 cases. This approach provided restoration of TK parameter values with less errors in tumor regions of interest, an improvement compared to a state-of-the-art indirect method. Applied prospectively, the proposed method provided whole-brain, high-resolution TK maps with good image quality.

CONCLUSION

Model-based direct estimation of TK maps from k,t-space DCE-MRI data is feasible and is compatible up to 100-fold undersampling. Magn Reson Med 78:1566-1578, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Comparison of conventional DCE-MRI and a novel golden-angle radial multicoil compressed sensing method for the evaluation of breast lesion conspicuity

PURPOSE

To compare a novel multicoil compressed sensing technique with flexible temporal resolution, golden-angle radial sparse parallel (GRASP), to conventional fat-suppressed spoiled three-dimensional (3D) gradient-echo (volumetric interpolated breath-hold examination, VIBE) MRI in evaluating the conspicuity of benign and malignant breast lesions.

MATERIALS AND METHODS

Between March and August 2015, 121 women (24-84 years; mean, 49.7 years) with 180 biopsy-proven benign and malignant lesions were imaged consecutively at 3.0 Tesla in a dynamic contrast-enhanced (DCE) MRI exam using sagittal T1-weighted fat-suppressed 3D VIBE in this Health Insurance Portability and Accountability Act-compliant, retrospective study. Subjects underwent MRI-guided breast biopsy (mean, 13 days [1-95 days]) using GRASP DCE-MRI, a fat-suppressed radial "stack-of-stars" 3D FLASH sequence with golden-angle ordering. Three readers independently evaluated breast lesions on both sequences. Statistical analysis included mixed models with generalized estimating equations, kappa-weighted coefficients and Fisher's exact test.

RESULTS

All lesions demonstrated good conspicuity on VIBE and GRASP sequences (4.28 ± 0.81 versus 3.65 ± 1.22), with no significant difference in lesion detection (P = 0.248). VIBE had slightly higher lesion conspicuity than GRASP for all lesions, with VIBE 12.6% (0.63/5.0) more conspicuous (P < 0.001). Masses and nonmass enhancement (NME) were more conspicuous on VIBE (P < 0.001), with a larger difference for NME (14.2% versus 9.4% more conspicuous). Malignant lesions were more conspicuous than benign lesions (P < 0.001) on both sequences.

CONCLUSION

GRASP DCE-MRI, a multicoil compressed sensing technique with high spatial resolution and flexible temporal resolution, has near-comparable performance to conventional VIBE imaging for breast lesion evaluation.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 3 J. MAGN. RESON. IMAGING 2017;45:1746-1752.

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

PURPOSE

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

THEORY AND METHODS

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

RESULTS

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

CONCLUSION

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

Hepatic function imaging using dynamic Gd-EOB-DTPA enhanced MRI and pharmacokinetic modeling

PURPOSE

To determine whether pharmacokinetic modeling parameters with different output assumptions of dynamic contrast-enhanced MRI (DCE-MRI) using Gd-EOB-DTPA correlate with serum-based liver function tests, and compare the goodness of fit of the different output assumptions.

METHODS

A 6-min DCE-MRI protocol was performed in 38 patients. Four dual-input two-compartment models with different output assumptions and a published one-compartment model were used to calculate hepatic function parameters. The Akaike information criterion fitting error was used to evaluate the goodness of fit. Imaging-based hepatic function parameters were compared with blood chemistry using correlation with multiple comparison correction.

RESULTS

The dual-input two-compartment model assuming venous flow equals arterial flow plus portal venous flow and no bile duct output better described the liver tissue enhancement with low fitting error and high correlation with blood chemistry. The relative uptake rate Kir derived from this model was found to be significantly correlated with direct bilirubin (r = -0.52, P = 0.015), prealbumin concentration (r = 0.58, P = 0.015), and prothrombin time (r = -0.51, P = 0.026).

CONCLUSION

It is feasible to evaluate hepatic function by proper output assumptions. The relative uptake rate has the potential to serve as a biomarker of function. Magn Reson Med 78:1488-1495, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Voxelwise analysis of simultaneously acquired and spatially correlated 18 F-fluorodeoxyglucose (FDG)-PET and intravoxel incoherent motion metrics in breast cancer

PURPOSE

Diffusion-weighted imaging (DWI) and ¹⁸ F-fluorodeoxyglucose-positron emission tomography (¹⁸ F-FDG-PET) independently correlate with malignancy in breast cancer, but the relationship between their structural and metabolic metrics is not completely understood. This study spatially correlates diffusion, perfusion, and glucose avidity in breast cancer with simultaneous PET/MR imaging and compares correlations with clinical prognostics.

METHODS

In this Health Insurance Portability and Accountability Act-compliant prospective study, with written informed consent and approval of the institutional review board and using simultaneously acquired FDG-PET and DWI, tissue diffusion (D_t ), and perfusion fraction (f_p ) from intravoxel incoherent motion (IVIM) analysis were registered to FDG-PET within 14 locally advanced breast cancers. Lesions were analyzed using 2D histograms and correlation coefficients between D_t , f_p , and standardized uptake value (SUV). Correlations were compared with prognostics from biopsy, metastatic burden from whole-body PET, and treatment history.

RESULTS

SUV||D_t correlation coefficient significantly distinguished treated (0.11 ± 0.24) from nontreated (-0.33 ± 0.26) patients (P = 0.005). SUV||f_p correlations were on average negative for the whole cohort (-0.17 ± 0.13).

CONCLUSION

Simultaneously acquired and registered FDG-PET/DWI allowed quantifiable descriptions of breast cancer microenvironments that may provide a framework for monitoring and predicting response to treatment. Magn Reson Med 78:1147-1156, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Separation of benign and malignant breast lesions using dynamic contrast enhanced MRI in a biopsy cohort

PURPOSE

To assess the diagnostic utility of contrast kinetic analysis for breast lesions and background parenchyma of women undergoing MRI-guided biopsies, for whom standard clinical analysis had failed to separate benign and malignant lesions.

MATERIALS AND METHODS

This study included 115 women who had indeterminate lesions based on routine diagnostic breast MRI exams and underwent an MRI (3 Tesla) -guided biopsy of one or more lesions suspicious for breast cancer. Breast dynamic contrast-enhanced (DCE)-MRI was performed using a radial stack-of-stars three-dimensional spoiled gradient echo pulse sequence and modified k-space weighted image contrast image reconstruction. Contrast kinetic model analysis was conducted to characterize the contrast enhancement patterns measured in lesions and background parenchyma (BP). The transfer rate (K^trans ), interstitial volume fraction (v_e ), and vascular volume fraction (v_p ) estimated from the lesion and BP were used to separate malignant from benign lesions.

RESULTS

The patients with malignant lesions had significantly (P < 0.05) higher median lesion-K^trans (0.081 min^-1 ), higher median BP-K^trans (0.032 min^-1 ), and BP-v_p (0.020) than those without malignant lesions (0.056 min^-1 , 0.017 min^-1 and 0.012, respectively). The area under the receiver operating characteristic curve (AUC) of the BP-K^trans (0.687) was highest among the single parameters and higher than that of the lesion-K^trans (0.664). The combined logistic regression model of lesion-K^trans , lesion-v_e , BP-K^trans , BP-v_e , and BP-v_p had the highest AUC of 0.730.

CONCLUSION

Our results suggest that the contrast kinetic analysis of DCE-MRI data can be used to differentiate the malignant lesions from the benign and high-risk lesions among the indeterminate breast lesions recommended for MRI-guided biopsy exams.

LEVEL OF EVIDENCE

3 J. MAGN. RESON. IMAGING 2017;45:1385-1393.

Motion correction based reconstruction method for compressively sampled cardiac MR imaging

Respiratory motion during Magnetic Resonance (MR) acquisition causes strong blurring artifacts in the reconstructed images. These artifacts become more pronounced when used with the fast imaging reconstruction techniques like compressed sensing (CS). Recently, an MR reconstruction technique has been done with the help of compressed sensing (CS), to provide good quality sparse images from the highly under-sampled k-space data. In order to maximize the benefits of CS, it is obvious to use CS with the motion corrected samples. In this paper, we propose a novel CS based motion corrected image reconstruction technique. First, k-space data have been assigned to different respiratory state with the help of frequency domain phase correlation method. Then, multiple sparsity constraints has been used to provide good quality reconstructed cardiac cine images with the highly under-sampled k-space data. The proposed method exploits the multiple sparsity constraints, in combination with demon based registration technique and a novel reconstruction technique to provide the final motion free images. The proposed method is very simple to implement in clinical settings as compared to existing motion corrected methods. The performance of the proposed method is examined using simulated data and clinical data. Results show that this method performs better than the reconstruction of CS based method of cardiac cine images. Different acceleration rates have been used to show the performance of the proposed method.

Dynamic contrast-enhanced MRI of the prostate: An intraindividual assessment of the effect of temporal resolution on qualitative detection and quantitative analysis of histopathologically proven prostate cancer

PURPOSE

To assess the effects of temporal resolution (R_T ) in dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) on qualitative tumor detection and quantitative pharmacokinetic parameters in prostate cancer.

MATERIALS AND METHODS

This retrospective Institutional Review Board (IRB)-approved study included 58 men (64 ± 7 years). They underwent 3T prostate MRI showing dominant peripheral zone (PZ) tumors (24 with Gleason ≥ 4 + 3), prior to prostatectomy. Continuously acquired DCE utilizing GRASP (Golden-angle RAdial Sparse Parallel) was retrospectively reconstructed at R_T of 1.4 sec, 3.7 sec, 6.0 sec, 9.7 sec, and 14.9 sec. A reader placed volumes-of-interest on dominant tumors and benign PZ, generating quantitative pharmacokinetic parameters (k^trans , v_e ) at each R_T . Two blinded readers assessed each R_T for lesion presence, location, conspicuity, and reader confidence on a 5-point scale. Data were assessed by mixed-model analysis of variance (ANOVA), generalized estimating equation (GEE), and receiver operating characteristic (ROC) analysis.

RESULTS

R_T did not affect sensitivity (R1_all : 69.0%-72.4%, all P_adj  = 1.000; R1_{GS≥4 + 3} : 83.3-91.7%, all P_adj  = 1.000; R2_all : 60.3-69.0%, all P_adj  = 1.000; R2_{GS≥4 + 3} : 58.3%-79.2%, all P_adj  = 1.000). R1 reported greater conspicuity of GS ≥ 4 + 3 tumors at R_T of 1.4 sec vs. 14.9 sec (4.29 ± 1.23 vs. 3.46 ± 1.44; P_adj  = 0.029). No other tumor conspicuity pairwise comparison reached significance (R1_all : 2.98-3.43, all P_adj ≥ 0.205; R2_all : 2.57-3.19, all P_adj ≥ 0.059; R1_{GS≥4 + 3} : 3.46-4.29, all other P_adj ≥ 0.156; R2_{GS≥4 + 3} : 2.92-3.71, all P_adj ≥ 0.439). There was no effect of R_T on reader confidence (R1_all : 3.17-3.34, all P_adj  = 1.000; R2_all : 2.83-3.19, all P_adj ≥ 0.801; R1_{GS≥4 + 3} : 3.79-4.21, all P_adj  = 1.000; R2_{GS≥4 + 3} : 3.13-3.79, all P_adj  = 1.000). k^trans and v_e of tumor and benign tissue did not differ across R_T (all adjusted P values [P_adj ] = 1.000). R_T did not significantly affect area under the curve (AUC) of K^trans or v_e for differentiating tumor from benign (all P_adj  = 1.000).

CONCLUSION

Current PI-RADS recommendations for R_T of 10 seconds may be sufficient, with further reduction to the stated PI-RADS preference of R_T ≤ 7 seconds offering no benefit in tumor detection or quantitative analysis.

LEVEL OF EVIDENCE

3 J. MAGN. RESON. IMAGING 2017;45:1464-1475.

Free-breathing volumetric fat/water separation by combining radial sampling, compressed sensing, and parallel imaging

PURPOSE

Conventional fat/water separation techniques require that patients hold breath during abdominal acquisitions, which often fails and limits the achievable spatial resolution and anatomic coverage. This work presents a novel approach for free-breathing volumetric fat/water separation.

METHODS

Multiecho data are acquired using a motion-robust radial stack-of-stars three-dimensional GRE sequence with bipolar readout. To obtain fat/water maps, a model-based reconstruction is used that accounts for the off-resonant blurring of fat and integrates both compressed sensing and parallel imaging. The approach additionally enables generation of respiration-resolved fat/water maps by detecting motion from k-space data and reconstructing different respiration states. Furthermore, an extension is described for dynamic contrast-enhanced fat-water-separated measurements.

RESULTS

Uniform and robust fat/water separation is demonstrated in several clinical applications, including free-breathing noncontrast abdominal examination of adults and a pediatric subject with both motion-averaged and motion-resolved reconstructions, as well as in a noncontrast breast exam. Furthermore, dynamic contrast-enhanced fat/water imaging with high temporal resolution is demonstrated in the abdomen and breast.

CONCLUSION

The described framework provides a viable approach for motion-robust fat/water separation and promises particular value for clinical applications that are currently limited by the breath-holding capacity or cooperation of patients. Magn Reson Med 78:565-576, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

A Projection Algorithm for Gradient Waveforms Design in Magnetic Resonance Imaging

Collecting the maximal amount of information in a given scanning time is a major concern in magnetic resonance imaging (MRI) to speed up image acquisition. The hardware constraints (gradient magnitude, slew rate, etc.), physical distortions (e.g., off-resonance effects) and sampling theorems (Shannon, compressed sensing) must be taken into account simultaneously, which makes this problem extremely challenging. To date, the main approach to design gradient waveform has consisted of selecting an initial shape (e.g., spiral, radial lines, etc.) and then traversing it as fast as possible using optimal control. In this paper, we propose an alternative solution which first consists of defining a desired parameterization of the trajectory and then of optimizing for minimal deviation of the sampling points within gradient constraints. This method has various advantages. First, it better preserves the density of the input curve which is critical in sampling theory. Second, it allows to smooth high curvature areas making the acquisition time shorter in some cases. Third, it can be used both in the Shannon and CS sampling theories. Last, the optimized trajectory is computed as the solution of an efficient iterative algorithm based on convex programming. For piecewise linear trajectories, as compared to optimal control reparameterization, our approach generates a gain in scanning time of 10% in echo planar imaging while improving image quality in terms of signal-to-noise ratio (SNR) by more than 6 dB. We also investigate original trajectories relying on traveling salesman problem solutions. In this context, the sampling patterns obtained using the proposed projection algorithm are shown to provide significantly better reconstructions (more than 6 dB) while lasting the same scanning time.

Self-gated MRI of multiple beat morphologies in the presence of arrhythmias

PURPOSE

Develop self-gated MRI for distinct heartbeat morphologies in subjects with arrhythmias.

METHODS

Golden angle radial data was obtained in seven sinus and eight arrhythmias subjects. An image-based cardiac navigator was derived from single-shot images, distinct beat types were identified, and images were reconstructed for repeated morphologies. Image sharpness, contrast, and volume variation were quantified and compared with self-gated MRI. Images were scored for image quality and artifacts. Hemodynamic parameters were computed for each distinct beat morphology in bigeminy and trigeminy subjects and for sinus beats in patients with infrequent premature ventricular contractions.

RESULTS

Images of distinct beat types were reconstructed except for two patients with infrequent premature ventricular contractions. Image contrast and sharpness were similar to sinus self-gated images (contrast = 0.45 ± 0.13 and 0.43 ± 0.15; sharpness = 0.21 ± 0.11 and 0.20 ± 0.05). Visual scoring was highest in self-gated images (4.1 ± 0.3) compared with real-time (3.9 ± 0.4) and ECG-gated cine (3.4 ± 1.5). ECG-gated cine had less artifacts than self-gating (2.3 ± 0.7 and 2.1 ± 0.2), but was affected by misgating in two subjects. Among arrhythmia subjects, post-extrasystole/sinus (58.1 ± 8.6 mL) and interrupted sinus (61.4 ± 5.9 mL) stroke volume was higher than extrasystole (32.0 ± 16.5 mL; P < 0.02).

CONCLUSION

Self-gated imaging can reconstruct images during ectopy and allowed for quantification of hemodynamic function of different beat morphologies. Magn Reson Med 78:678-688, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Self-gated 4D multiphase, steady-state imaging with contrast enhancement (MUSIC) using rotating cartesian K-space (ROCK): Validation in children with congenital heart disease

PURPOSE

To develop and validate a cardiac-respiratory self-gating strategy for the recently proposed multiphase steady-state imaging with contrast enhancement (MUSIC) technique.

METHODS

The proposed SG strategy uses the ROtating Cartesian K-space (ROCK) sampling, which allows for retrospective k-space binning based on motion surrogates derived from k-space center line. The k-space bins are reconstructed using a compressed sensing algorithm. Ten pediatric patients underwent cardiac MRI for clinical reasons. The original MUSIC and 2D-CINE images were acquired as a part of the clinical protocol, followed by the ROCK-MUSIC acquisition, all under steady-state intravascular distribution of ferumoxytol. Subjective scores and image sharpness were used to compare the images of ROCK-MUSIC and original MUSIC.

RESULTS

All scans were completed successfully without complications. The ROCK-MUSIC acquisition took 5 ± 1 min, compared to 8 ± 2 min for the original MUSIC. Image scores of ROCK-MUSIC were significantly better than original MUSIC at the ventricular outflow tracts (3.9 ± 0.3 vs. 3.3 ± 0.6, P < 0.05). There was a strong trend toward superior image scores for ROCK-MUSIC in the other anatomic locations.

CONCLUSION

ROCK-MUSIC provided images of equal or superior image quality compared to original MUSIC, and this was achievable with 40% savings in scan time and without the need for physiologic signal. Magn Reson Med 78:472-483, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Radial simultaneous multi-slice CAIPI for ungated myocardial perfusion

OBJECTIVE

Simultaneous multi-slice (SMS) imaging is a slice acceleration technique that acquires multiple slices in the same time as a single slice. Radial controlled aliasing in parallel imaging results in higher acceleration (radial CAIPIRINHA or CAIPI) is a promising SMS method with less severe slice aliasing artifacts as compared to its Cartesian counterpart. Here we use radial CAIPI with data undersampling and constrained reconstruction to improve the utility of ungated cardiac perfusion acquisitions. We test the proposed framework with a traditional saturation recovery fast low-angle shot (turboFLASH) sequence and also without saturation recovery as a steady-state spoiled gradient echo (SPGR) sequence on animal and human studies.

METHODS

Simulations and phantom studies were performed for both the turboFLASH and the SPGR radial CAIPI methods. Ungated undersampled golden ratio radial CAIPI data with saturation recovery were acquired in 8 dogs and 2 human subjects. The CAIPI data without saturation pulses were acquired in 4 human subjects. For both methods, slice acceleration factors of two and three were used. A new spatio-temporal reconstruction using total variation and patch-based low rank constraints was used to jointly reconstruct the multi-slice multi-coil images.

RESULTS

Phantom scans and computer simulations showed that ungated SPGR generally provides better contrast to noise ratio (CNR) than the saturation recovery sequence if the saturation recovery time is less than 100ms. Both of the ungated radial CAIPI methods demonstrated promising image quality in terms of preserving dynamics of the contrast agent and maintaining anatomical structures, even with three slices acquired simultaneously.

CONCLUSION

Ungated simultaneous multi-slice acquisitions with either a saturation recovery turboFLASH sequence or a steady-state gradient echo SPGR sequence are feasible and provide increased slice coverage without loss of temporal resolution. Compared with a sensitivity encoding (SENSE) SMS reconstruction, the constrained reconstruction method provides better image quality for undersampled radial CAIPI data.

Surrogate arterial phase imaging using a long duration (≈1.5 min) radial acquisition T1-weighted sequence: an alternative in patients unable to breath-hold

BACKGROUND

Pediatric and adult patients unable to suspend respiration generally undergo magnetic resonance (MR) examinations that lack arterial phase imaging, which is a phase that provides substantial information on disease processes. An MR strategy that provides this type of information may be of considerable value.

PURPOSE

To describe and assess the feasibility and enhancement quality of early-phase imaging utilizing long-duration radial 3D-GRE imaging by initiating the sequence prior to starting contrast injection.

MATERIAL AND METHODS

Thirty-three consecutive patients (10 men, 23 women; 50.7 ± 25.5 years) underwent free-breathing gadolinium-enhanced radial 3D-GRE, with sequence initiation 30 s prior to contrast injection. Late hepatic arterial (LHA) phase was chosen for comparison. Images were evaluated for enhancement and overall image quality. Organ enhancement was calculated. Sub-group analysis was performed.

RESULTS

Twenty-two examinations of radial 3D-GRE sequences were acquired during the LHA phase. Organ enhancement scores were of satisfactory to good quality (range, 3.32-3.82). There was a significant trend of superior overall enhancement quality scores in pediatrics and examinations performed at 3 T (P = 0.0225 and 0.0001, respectively).

CONCLUSION

Arterial phase abdominal MR imaging is feasible using conventional radial 3D-GRE by adopting this simplistic proposed approach, which may allow arterial-phase imaging in patients unable to breath-hold.

Respiratory Motion-Resolved Compressed Sensing Reconstruction of Free-Breathing Radial Acquisition for Dynamic Liver Magnetic Resonance Imaging

OBJECTIVE

This study aimed to demonstrate feasibility of free-breathing radial acquisition with respiratory motion-resolved compressed sensing reconstruction [extra-dimensional golden-angle radial sparse parallel imaging (XD-GRASP)] for multiphase dynamic gadolinium ethoxybenzyl diethylenetriamine pentaacetic acid (Gd-EOB-DTPA)-enhanced liver imaging, and to compare image quality to compressed sensing reconstruction with respiratory motion-averaging (GRASP) and prior conventional breath-held Cartesian-sampled data sets [BH volume interpolated breath-hold examination (VIBE)] in same patients.

SUBJECTS AND METHODS

In this Health Insurance Portability and Accountability Act-compliant prospective study, 16 subjects underwent free-breathing continuous radial acquisition during Gd-EOB-DTPA injection and had prior BH-VIBE available. Acquired data were reconstructed using motion-averaging GRASP approach in which consecutive 84 spokes were grouped in each contrast-enhanced phase for a temporal resolution of approximately 14 seconds. Additionally, respiratory motion-resolved reconstruction was performed from the same k-space data by sorting each contrast-enhanced phase into multiple respiratory motion states using compressed sensing algorithm named XD-GRASP, which exploits sparsity along both the contrast-enhancement and respiratory-state dimensions.Contrast-enhanced dynamic multiphase XD-GRASP, GRASP, and BH-VIBE images were anonymized, pooled together in a random order, and presented to 2 board-certified radiologists for independent evaluation of image quality, with higher score indicating more optimal examination.

RESULTS

The XD-GRASP reconstructions had significantly (all P < 0.05) higher overall image quality scores compared to GRASP for early arterial (reader 1: 4.3 ± 0.6 vs 3.31 ± 0.6; reader 2: 3.81 ± 0.8 vs 3.38 ± 0.9) and late arterial (reader 1: 4.5 ± 0.6 vs 3.63 ± 0.6; reader 2: 3.56 ± 0.5 vs 2.88 ± 0.7) phases of enhancement for both readers. The XD-GRASP also had higher overall image quality score in portal venous phase, which was significant for reader 1 (4.44 ± 0.5 vs 3.75 ± 0.8; P = 0.002). In addition, the XD-GRASP had higher overall image quality score compared to BH-VIBE for early (reader 1: 4.3 ± 0.6 vs 3.88 ± 0.6; reader 2: 3.81 ± 0.8 vs 3.50 ± 1.0) and late (reader 1: 4.5 ± 0.6 vs 3.44 ± 0.6; reader 2: 3.56 ± 0.5 vs 2.94 ± 0.9) arterial phases.

CONCLUSION

Free-breathing motion-resolved XD-GRASP reconstructions provide diagnostic high-quality multiphase images in patients undergoing Gd-EOB-DTPA-enhanced liver examination.

Evaluation of upper airway collapsibility using real-time MRI

PURPOSE

To develop and demonstrate a real-time MRI method for assessing upper airway collapsibility in sleep apnea.

MATERIALS AND METHODS

Data were acquired on a clinical 3 Tesla scanner using a radial CAIPIRIHNA sequence with modified golden angle view ordering and reconstructed using parallel imaging and compressed sensing with temporal finite difference sparsity constraint. Segmented airway areas together with synchronized facemask pressure were used to calculate airway compliance and projected closing pressure, Pclose , at four axial locations along the upper airway. This technique was demonstrated in five adolescent obstructive sleep apnea (OSA) patients, three adult OSA patients and four healthy volunteers. Heart rate, oxygen saturation, facemask pressure, and abdominal/chest movements were monitored in real-time during the experiments to determine sleep/wakefulness.

RESULTS

Student's t-tests showed that both compliance and Pclose were significantly different between healthy controls and OSA patients (P < 0.001). The results also suggested that a narrower airway site does not always correspond to higher collapsibility.

CONCLUSION

With the proposed methods, both compliance and Pclose can be calculated and used to quantify airway collapsibility in OSA with an awake scan of 30 min total scan room time. J. Magn. Reson. Imaging 2016;44:158-167.

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.

Compressed sensing for longitudinal MRI: An adaptive-weighted approach

PURPOSE

Repeated brain MRI scans are performed in many clinical scenarios, such as follow up of patients with tumors and therapy response assessment. In this paper, the authors show an approach to utilize former scans of the patient for the acceleration of repeated MRI scans.

METHODS

The proposed approach utilizes the possible similarity of the repeated scans in longitudinal MRI studies. Since similarity is not guaranteed, sampling and reconstruction are adjusted during acquisition to match the actual similarity between the scans. The baseline MR scan is utilized both in the sampling stage, via adaptive sampling, and in the reconstruction stage, with weighted reconstruction. In adaptive sampling, k-space sampling locations are optimized during acquisition. Weighted reconstruction uses the locations of the nonzero coefficients in the sparse domains as a prior in the recovery process. The approach was tested on 2D and 3D MRI scans of patients with brain tumors.

RESULTS

The longitudinal adaptive compressed sensing MRI (LACS-MRI) scheme provides reconstruction quality which outperforms other CS-based approaches for rapid MRI. Examples are shown on patients with brain tumors and demonstrate improved spatial resolution. Compared with data sampled at the Nyquist rate, LACS-MRI exhibits signal-to-error ratio (SER) of 24.8 dB with undersampling factor of 16.6 in 3D MRI.

CONCLUSIONS

The authors presented an adaptive method for image reconstruction utilizing similarity of scans in longitudinal MRI studies, where possible. The proposed approach can significantly reduce scanning time in many applications that consist of disease follow-up and monitoring of longitudinal changes in brain MRI.

Fast 3D ultrashort echo-time spiral projection imaging using golden-angle: A flexible protocol for in vivo mouse imaging at high magnetic field

PURPOSE

To develop a fast three-dimensional (3D) k-space encoding method based on spiral projection imaging (SPI) with an interleaved golden-angle approach and to validate this novel sequence on small animal models.

METHODS

A disk-like trajectory, in which each disk contained spirals, was developed. The 3D encoding was performed by tilting the disks with a golden angle. The sharpness was first calculated at different T2* values. Then, the sharpness was measured on phantom using variable undersampling ratios. Finally, the sampling method was validated by whole brain time-of-flight angiography and ultrasmall superparamagnetic iron oxide (USPIO) enhanced free-breathing liver angiography on mouse.

RESULTS

The in vitro results demonstrated the robustness of the method for short T2* and high undersampling ratios. In vivo experiments showed the ability to properly detect small vessels in the brain with an acquisition time shorter than 1 min. Free-breathing mice liver angiography showed the insensitivity of this protocol toward motions and flow artifacts, and enabled the visualization of liver motion during breathing.

CONCLUSIONS

The method implemented here allowed fast 3D k-space sampling with a high undersampling ratio. Combining the advantages of center-out spirals with the flexibility of the golden angle approach could have major implications for real-time imaging. Magn Reson Med 77:1831-1840, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Minimum Field Strength Simulator for Proton Density Weighted MRI

OBJECTIVE

To develop and evaluate a framework for simulating low-field proton-density weighted MRI acquisitions based on high-field acquisitions, which could be used to predict the minimum B0 field strength requirements for MRI techniques. This framework would be particularly useful in the evaluation of de-noising and constrained reconstruction techniques.

MATERIALS AND METHODS

Given MRI raw data, lower field MRI acquisitions can be simulated based on the signal and noise scaling with field strength. Certain assumptions are imposed for the simulation and their validity is discussed. A validation experiment was performed using a standard resolution phantom imaged at 0.35 T, 1.5 T, 3 T, and 7 T. This framework was then applied to two sample proton-density weighted MRI applications that demonstrated estimation of minimum field strength requirements: real-time upper airway imaging and liver proton-density fat fraction measurement.

RESULTS

The phantom experiment showed good agreement between simulated and measured images. The SNR difference between simulated and measured was ≤ 8% for the 1.5T, 3T, and 7T cases which utilized scanners with the same geometry and from the same vendor. The measured SNR at 0.35T was 1.8- to 2.5-fold less than predicted likely due to unaccounted differences in the RF receive chain. The predicted minimum field strength requirements for the two sample applications were 0.2 T and 0.3 T, respectively.

CONCLUSIONS

Under certain assumptions, low-field MRI acquisitions can be simulated from high-field MRI data. This enables prediction of the minimum field strength requirements for a broad range of MRI techniques.

Accelerated MRI for the assessment of cardiac function

Heart disease is a worldwide public health problem; assessment of cardiac function is an important part of the diagnosis and management of heart disease. MRI of the heart can provide clinically useful information on cardiac function, although it is still not routinely used in clinical practice, in part because of limited imaging speed. New accelerated methods for performing cardiovascular MRI (CMR) have the potential to provide both increased imaging speed and robustness to CMR, as well as access to increased functional information. In this review, we will briefly discuss the main methods currently employed to accelerate CMR methods, such as parallel imaging, k-t undersampling and compressed sensing, as well as new approaches that extend the idea of compressed sensing and exploit sparsity to provide richer information of potential use in clinical practice.

On the influence of respiratory motion in radial tissue phase mapping cardiac MRI

PURPOSE

To investigate the impact of respiratory motion on radial tissue phase mapping (TPM) measurements, and to improve image quality and scan efficiency without compromising velocity fidelity by increasing the respiratory acceptance window with and without motion correction.

MATERIALS AND METHODS

A radial golden angle TPM sequence was measured in 10 healthy volunteers in three short axis slices at 3T. Ungated ( CFREE), self-gated with a single acceptance window ( CREF), motion-corrected averaging using all ( CMCall), or selected ( CMC) data reconstructions were compared by means of various image quality measures and resulting velocities.

RESULTS

Using all data ( CFREE) resulted in significantly higher perceived signal-to-noise ratio (SNR) (P < 0.001), but significantly reduced sharpness (P < 0.001) and contrast (P = 0.02), when compared to CREF. Coefficient of variation (CV) and perceived sharpness were not significantly different (P > 0.05). With motion-correction, perceived sharpness could be significantly improved ( CMC: P = 0.002; CMCall: P = 0.002) in comparison to CFREE. Velocity peaks of CFREE were significantly reduced compared to CREF (all peaks: P < 0.001; except the longitudinal "E" peak: P = 0.03). The peak velocities in CMC and CMCall were not significantly different from CREF (all peaks: P > 0.08; except longitudinal "E"/"A" peaks: P > 0.01).

CONCLUSION

Free-breathing reconstruction results in good perceived image sharpness and velocity information with slightly, but significantly, reduced peak velocities. For achieving velocities and image quality comparable to data from a single acceptance window, but higher gating efficiency, selected motion-corrected TPM (CMC) can be applied. J. Magn. Reson. Imaging 2016;44:1218-1228.

Quantification of Left Ventricular Function With Premature Ventricular Complexes Reveals Variable Hemodynamics

BACKGROUND

Premature ventricular complexes (PVCs) are prevalent in the general population and are sometimes associated with reduced ventricular function. Current echocardiographic and cardiovascular magnetic resonance imaging techniques do not adequately address the effect of PVCs on left ventricular function.

METHODS AND RESULTS

Fifteen subjects with a history of frequent PVCs undergoing cardiovascular magnetic resonance imaging had real-time slice volume quantification performed using a 2-dimensional (2D) real-time cardiovascular magnetic resonance imaging technique. Synchronization of 2D real-time imaging with patient ECG allowed for different beats to be categorized by the loading beat RR duration and beat RR duration. For each beat type, global volumes were quantified via summation over all slices covering the entire ventricle. Different patterns of ectopy, including isolated PVCs, bigeminy, trigeminy, and interpolated PVCs, were observed. Global functional measurement of the different beat types based on timing demonstrated differences in preload, stroke volume, and ejection fraction. An average of hemodynamic function was quantified for each subject depending on the frequency of each observed beat type.

CONCLUSIONS

Application of real-time cardiovascular magnetic resonance imaging in patients with PVCs revealed differential contribution of PVCs to hemodynamics.

Pulmonary imaging using respiratory motion compensated simultaneous PET/MR

PURPOSE

Pulmonary positron emission tomography (PET) imaging is confounded by blurring artifacts caused by respiratory motion. These artifacts degrade both image quality and quantitative accuracy. In this paper, the authors present a complete data acquisition and processing framework for respiratory motion compensated image reconstruction (MCIR) using simultaneous whole body PET/magnetic resonance (MR) and validate it through simulation and clinical patient studies.

METHODS

The authors have developed an MCIR framework based on maximum a posteriori or MAP estimation. For fast acquisition of high quality 4D MR images, the authors developed a novel Golden-angle RAdial Navigated Gradient Echo (GRANGE) pulse sequence and used it in conjunction with sparsity-enforcing k-t FOCUSS reconstruction. The authors use a 1D slice-projection navigator signal encapsulated within this pulse sequence along with a histogram-based gate assignment technique to retrospectively sort the MR and PET data into individual gates. The authors compute deformation fields for each gate via nonrigid registration. The deformation fields are incorporated into the PET data model as well as utilized for generating dynamic attenuation maps. The framework was validated using simulation studies on the 4D XCAT phantom and three clinical patient studies that were performed on the Biograph mMR, a simultaneous whole body PET/MR scanner.

RESULTS

The authors compared MCIR (MC) results with ungated (UG) and one-gate (OG) reconstruction results. The XCAT study revealed contrast-to-noise ratio (CNR) improvements for MC relative to UG in the range of 21%-107% for 14 mm diameter lung lesions and 39%-120% for 10 mm diameter lung lesions. A strategy for regularization parameter selection was proposed, validated using XCAT simulations, and applied to the clinical studies. The authors' results show that the MC image yields 19%-190% increase in the CNR of high-intensity features of interest affected by respiratory motion relative to UG and a 6%-51% increase relative to OG.

CONCLUSIONS

Standalone MR is not the traditional choice for lung scans due to the low proton density, high magnetic susceptibility, and low T2 (∗) relaxation time in the lungs. By developing and validating this PET/MR pulmonary imaging framework, the authors show that simultaneous PET/MR, unique in its capability of combining structural information from MR with functional information from PET, shows promise in pulmonary imaging.

4D respiratory motion-compensated image reconstruction of free-breathing radial MR data with very high undersampling

PURPOSE

To develop four-dimensional (4D) respiratory time-resolved MRI based on free-breathing acquisition of radial MR data with very high undersampling.

METHODS

We propose the 4D joint motion-compensated high-dimensional total variation (4D joint MoCo-HDTV) algorithm, which alternates between motion-compensated image reconstruction and artifact-robust motion estimation at multiple resolution levels. The algorithm is applied to radial MR data of the thorax and upper abdomen of 12 free-breathing subjects with acquisition times between 37 and 41 s and undersampling factors of 16.8. Resulting images are compared with compressed sensing-based 4D motion-adaptive spatio-temporal regularization (MASTeR) and 4D high-dimensional total variation (HDTV) reconstructions.

RESULTS

For all subjects, 4D joint MoCo-HDTV achieves higher similarity in terms of normalized mutual information and cross-correlation than 4D MASTeR and 4D HDTV when compared with reference 4D gated gridding reconstructions with 8.4 ± 1.1 times longer acquisition times. In a qualitative assessment of artifact level and image sharpness by two radiologists, 4D joint MoCo-HDTV reveals higher scores (P < 0.05) than 4D HDTV and 4D MASTeR at the same undersampling factor and the reference 4D gated gridding reconstructions, respectively.

CONCLUSIONS

4D joint MoCo-HDTV enables time-resolved image reconstruction of free-breathing radial MR data with undersampling factors of 16.8 while achieving low-streak artifact levels and high image sharpness. Magn Reson Med 77:1170-1183, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Three-dimensional MR Cholangiopancreatography in a Breath Hold with Sparsity-based Reconstruction of Highly Undersampled Data

Purpose To develop a three-dimensional breath-hold (BH) magnetic resonance (MR) cholangiopancreatographic protocol with sampling perfection with application-optimized contrast using different flip-angle evolutions (SPACE) acquisition and sparsity-based iterative reconstruction (SPARSE) of prospectively sampled 5% k-space data and to compare the results with conventional respiratory-triggered (RT) acquisition. Materials and Methods This HIPAA-compliant prospective study was institutional review board approved. Twenty-nine patients underwent conventional RT SPACE and BH-accelerated SPACE acquisition with 5% k-space sampling at 3 T. Spatial resolution and other parameters were matched when possible. BH SPACE images were reconstructed by enforcing joint multicoil sparsity in the wavelet domain (SPARSE-SPACE). Two board-certified radiologists independently evaluated BH SPARSE-SPACE and RT SPACE images for image quality parameters in the pancreatic duct and common bile duct by using a five-point scale. The Wilcoxon signed-rank test was used to compare BH SPARSE-SPACE and RT SPACE images. Results Acquisition time for BH SPARSE-SPACE was 20 seconds, which was significantly (P < .001) shorter than that for RT SPACE (mean ± standard deviation, 338.8 sec ± 69.1). Overall image quality scores were higher for BH SPARSE-SPACE than for RT SPACE images for both readers for the proximal, middle, and distal pancreatic duct, but the difference was not statistically significant (P > .05). For reader 1, distal common bile duct scores were significantly higher with BH SPARSE-SPACE acquisition (P = .036). More patients had acceptable or better overall image quality (scores ≥ 3) with BH SPARSE-SPACE than with RT SPACE acquisition, respectively, for the proximal (23 of 29 [79%] vs 22 of 29 [76%]), middle (22 of 29 [76%] vs 18 of 29 [62%]), and distal (20 of 29 [69%] vs 13 of 29 [45%]) pancreatic duct and the proximal (25 of 28 [89%] vs 22 of 28 [79%]) and distal (25 of 28 [89%] vs 24 of 28 [86%]) common bile duct. Conclusion BH SPARSE-SPACE showed similar or superior image quality for the pancreatic and common duct compared with that of RT SPACE despite 17-fold shorter acquisition time. (©) RSNA, 2016.

Development of Real-Time Magnetic Resonance Imaging of Mouse Hearts at 9.4 Tesla--Simulations and First Application

A novel method for real-time magnetic resonance imaging for the assessment of cardiac function in mice at 9.4 T is proposed. The technique combines a highly undersampled radial gradient echo acquisition with an image reconstruction utilizing both parallel imaging and compressed sensing. Simulations on an in silico phantom were performed to determine the achievable acceleration factor and to optimize regularization parameters. Several parameters characterizing the quality of the reconstructed images (such as spatial and temporal image sharpness or compartment areas) were calculated for this purpose. Subsequently, double-gated segmented cine data as well as non-gated undersampled real-time data using only six projections per timeframe (temporal resolution  ∼ 10 ms) were acquired in a mid-ventricular slice of four normal mouse hearts in vivo. The highly accelerated data sets were then subjected to the introduced reconstruction technique and results were validated against the fully sampled references. Functional parameters obtained from real-time and fully sampled data agreed well with a comparable accuracy for left-ventricular volumes and a slightly larger scatter for mass. This study introduces and validates a real-time cine-MRI technique, which significantly reduces scan time in preclinical cardiac functional imaging and has the potential to investigate mouse models with abnormal heart rhythm.

Accelerating chemical exchange saturation transfer (CEST) MRI by combining compressed sensing and sensitivity encoding techniques

PURPOSE

To evaluate the feasibility of accelerated chemical-exchange-saturation-transfer (CEST) imaging using a combination of compressed sensing (CS) and sensitivity encoding (SENSE) at 3 Tesla.

THEORY AND METHODS

Two healthy volunteers and six high-grade glioma patients were recruited. Raw CEST image k-space data were acquired (with varied radiofrequency saturation power levels for the healthy volunteer study), and a sequential CS and SENSE reconstruction (CS-SENSE) was assessed. The MTR_asym (3.5 ppm) signals were compared with varied CS-SENSE acceleration factors.

RESULTS

In the healthy volunteer study, a CS-SENSE acceleration factor of R = 2 × 2 (CS × SENSE) was achieved without compromising the reconstructed MTR_asym (3.5 ppm) image quality. The MTR_asym (3.5 ppm) signals obtained from the CS-SENSE reconstruction with R = 2 × 2 were well preserved compared with the reference image (R = 2 for only SENSE). In the glioma patient study, the MTR_asym (3.5 ppm) signals were significantly higher in the tumor region (Gd-enhancing tumor core) than in the normal-appearing white matter (P < 0.001). There was no significant MTR_asym (3.5 ppm) difference between the reference image and CS-SENSE-reconstructed image in the acceleration factor of R = 2 × 2.

CONCLUSION

Combining the SENSE technique with CS (R = 2 × 2) enables considerable acceleration of CEST image acquisition and potentially has a wide range of clinical applications. Magn Reson Med 77:779-786, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Free-breathing liver perfusion imaging using 3-dimensional through-time spiral generalized autocalibrating partially parallel acquisition acceleration

OBJECTIVES

The goal of this study was to develop free-breathing high-spatiotemporal resolution dynamic contrast-enhanced liver magnetic resonance imaging using non-Cartesian parallel imaging acceleration, and quantitative liver perfusion mapping.

MATERIALS AND METHODS

This study was approved by the local institutional review board and written informed consent was obtained from all participants. Ten healthy subjects and 5 patients were scanned on a Siemens 3-T Skyra scanner. A stack-of-spirals trajectory was undersampled in-plane with a reduction factor of 6 and reconstructed using 3-dimensional (3D) through-time non-Cartesian generalized autocalibrating partially parallel acquisition. High-resolution 3D images were acquired with a true temporal resolution of 1.6 to 1.9 seconds while the subjects were breathing freely. A dual-input single-compartment model was used to retrieve liver perfusion parameters from dynamic contrast-enhanced magnetic resonance imaging data, which were coregistered using an algorithm designed to reduce the effects of dynamic contrast changes on registration. Image quality evaluation was performed on spiral images and conventional images from 5 healthy subjects.

RESULTS

Images with a spatial resolution of 1.9 × 1.9 × 3 mm3 were obtained with whole-liver coverage. With an imaging speed of better than 2 s/vol, free-breathing scans were achieved and dynamic changes in enhancement were captured. The overall image quality of free-breathing spiral images was slightly lower than that of conventional long breath-hold Cartesian images, but it provided clinically acceptable or better image quality. The free-breathing 3D images were registered with almost no residual motion in liver tissue. After the registration, quantitative whole-liver 3D perfusion maps were obtained and the perfusion parameters are all in good agreement with the literature.

CONCLUSIONS

This high-spatiotemporal resolution free-breathing 3D liver imaging technique allows voxelwise quantification of liver perfusion.

XD-GRASP: Golden-angle radial MRI with reconstruction of extra motion-state dimensions using compressed sensing

PURPOSE

To develop a novel framework for free-breathing MRI called XD-GRASP, which sorts dynamic data into extra motion-state dimensions using the self-navigation properties of radial imaging and reconstructs the multidimensional dataset using compressed sensing.

METHODS

Radial k-space data are continuously acquired using the golden-angle sampling scheme and sorted into multiple motion-states based on respiratory and/or cardiac motion signals derived directly from the data. The resulting undersampled multidimensional dataset is reconstructed using a compressed sensing approach that exploits sparsity along the new dynamic dimensions. The performance of XD-GRASP is demonstrated for free-breathing three-dimensional (3D) abdominal imaging, two-dimensional (2D) cardiac cine imaging and 3D dynamic contrast-enhanced (DCE) MRI of the liver, comparing against reconstructions without motion sorting in both healthy volunteers and patients.

RESULTS

XD-GRASP separates respiratory motion from cardiac motion in cardiac imaging, and respiratory motion from contrast enhancement in liver DCE-MRI, which improves image quality and reduces motion-blurring artifacts.

CONCLUSION

XD-GRASP represents a new use of sparsity for motion compensation and a novel way to handle motions in the context of a continuous acquisition paradigm. Instead of removing or correcting motion, extra motion-state dimensions are reconstructed, which improves image quality and also offers new physiological information of potential clinical value.

Multiparametric MR Imaging in Abdominal Malignancies

Modern MR imaging protocols can yield both anatomic and functional information for the assessment of hepatobiliary and pancreatic malignancies. Diffusion-weighted imaging is fully integrated into state-of-the-art protocols for tumor detection, characterization, and therapy monitoring. Hepatobiliary contrast agents have gained ground in the evaluation of focal liver lesions during the last years. Perfusion MR imaging is expected to have a central role for monitoring therapy in body tumors treated with antivascular drugs. Approaches such as Magnetic resonance (MR) elastography and (1)H-MR spectroscopy are still confined to research centers, but with the potential to grow in a short time frame.

A fast and flexible MRI system for the study of dynamic vocal tract shaping

PURPOSE

The aim of this work was to develop and evaluate an MRI-based system for study of dynamic vocal tract shaping during speech production, which provides high spatial and temporal resolution.

METHODS

The proposed system utilizes (a) custom eight-channel upper airway coils that have high sensitivity to upper airway regions of interest, (b) two-dimensional golden angle spiral gradient echo acquisition, (c) on-the-fly view-sharing reconstruction, and (d) off-line temporal finite difference constrained reconstruction. The system also provides simultaneous noise-cancelled and temporally aligned audio. The system is evaluated in 3 healthy volunteers, and 1 tongue cancer patient, with a broad range of speech tasks.

RESULTS

We report spatiotemporal resolutions of 2.4 × 2.4 mm² every 12 ms for single-slice imaging, and 2.4 × 2.4 mm² every 36 ms for three-slice imaging, which reflects roughly 7-fold acceleration over Nyquist sampling. This system demonstrates improved temporal fidelity in capturing rapid vocal tract shaping for tasks, such as producing consonant clusters in speech, and beat-boxing sounds. Novel acoustic-articulatory analysis was also demonstrated.

CONCLUSION

A synergistic combination of custom coils, spiral acquisitions, and constrained reconstruction enables visualization of rapid speech with high spatiotemporal resolution in multiple planes. Magn Reson Med 77:112-125, 2017. © 2016 Wiley Periodicals, Inc.

Recommendations for real-time speech MRI

Real-time magnetic resonance imaging (RT-MRI) is being increasingly used for speech and vocal production research studies. Several imaging protocols have emerged based on advances in RT-MRI acquisition, reconstruction, and audio-processing methods. This review summarizes the state-of-the-art, discusses technical considerations, and provides specific guidance for new groups entering this field. We provide recommendations for performing RT-MRI of the upper airway. This is a consensus statement stemming from the ISMRM-endorsed Speech MRI summit held in Los Angeles, February 2014. A major unmet need identified at the summit was the need for consensus on protocols that can be easily adapted by researchers equipped with conventional MRI systems. To this end, we provide a discussion of tradeoffs in RT-MRI in terms of acquisition requirements, a priori assumptions, artifacts, computational load, and performance for different speech tasks. We provide four recommended protocols and identify appropriate acquisition and reconstruction tools. We list pointers to open-source software that facilitate implementation. We conclude by discussing current open challenges in the methodological aspects of RT-MRI of speech.

Performance of simultaneous high temporal resolution quantitative perfusion imaging of bladder tumors and conventional multi-phase urography using a novel free-breathing continuously acquired radial compressed-sensing MRI sequence

PURPOSE

To investigate the feasibility of high temporal resolution quantitative perfusion imaging of bladder tumors performed simultaneously with conventional multi-phase MR urography (MRU) using a novel free-breathing continuously acquired radial MRI sequence with compressed-sensing reconstruction.

METHODS

22 patients with bladder lesions underwent MRU using GRASP (Golden-angle RAdial Sparse Parallel) acquisition. Multi-phase contrast-enhanced abdominopelvic GRASP was performed during free-breathing (1.4×1.4×3.0mm(3) voxel size; 3:44min acquisition). Two dynamic datasets were retrospectively reconstructed by combining different numbers of sequentially acquired spokes into each dynamic frame: 110 spokes per frame for 25-s temporal resolution (serving as conventional MRU for clinical interpretation) and 8 spokes per frame for 1.7-s resolution. Using 1.7-s resolution images, ROIs were placed within bladder lesions and normal bladder wall, a femoral artery arterial input function was generated, and the Generalized Kinetic Model was applied.

RESULTS

Biopsy/cystectomy demonstrated 16 bladder tumors (13 stage≥T2, 3 stage≤T1) and 6 benign lesions. All lesions were well visualized using 25-s clinical multi-phase images. Using 1.7-s resolution images, K(trans) was significantly higher in tumors (0.38±0.24) than normal bladder (0.12±0.02=8, p<0.001) or benign lesions (0.15±0.04, p=0.033). Ratio between K(trans) of lesions and normal bladder was nearly double for tumors than benign lesions (4.3±3.4 vs. 2.2±1.6), and K(trans) was nearly double in stage≥T2 than stage≤T1 tumors (0.44±0.24 vs. 0.24±0.24), although these did not approach significance (p=0.180-0.209), possibly related to small sample size.

CONCLUSION

GRASP allows simultaneous quantitative high temporal resolution perfusion of bladder lesions during clinical MRU examinations using only one contrast injection and without additional scan time.

GOCART: GOlden-angle CArtesian randomized time-resolved 3D MRI

PURPOSE

To develop and evaluate a novel 3D Cartesian sampling scheme which is well suited for time-resolved 3D MRI using parallel imaging and compressed sensing.

METHODS

The proposed sampling scheme, termed GOlden-angle CArtesian Randomized Time-resolved (GOCART) 3D MRI, is based on golden angle (GA) Cartesian sampling, with random sampling of the ky-kz phase encode locations along each Cartesian radial spoke. This method was evaluated in conjunction with constrained reconstruction of retrospectively and prospectively undersampled in-vivo dynamic contrast enhanced (DCE) MRI data and simulated phantom data.

RESULTS

In in-vivo retrospective studies and phantom simulations, images reconstructed from phase encodes defined by GOCART were equal to or superior to those with Poisson disc or GA sampling schemes. Typical GOCART sampling tables were generated in <100ms. GOCART has also been successfully utilized prospectively to produce clinically valuable whole-brain DCE-MRI images.

CONCLUSION

GOCART is a practical and efficient sampling scheme for time-resolved 3D MRI. It shows great potential for highly accelerated DCE-MRI and is well suited to modern reconstruction methods such as parallel imaging and compressed sensing.

Gradient-Modulated PETRA MRI

Image blurring due to off-resonance and fast T 2* signal decay is a common issue in radial ultrashort echo time MRI sequences. One solution is to use a higher readout bandwidth, but this may be impractical for some techniques like pointwise encoding time reduction with radial acquisition (PETRA), which is a hybrid method of zero echo time and single point imaging techniques. Specifically, PETRA has severe specific absorption rate (SAR) and radiofrequency (RF) pulse peak power limitations when using higher bandwidths in human measurements. In this study, we introduce gradient modulation (GM) to PETRA to reduce image blurring artifacts while keeping SAR and RF peak power low. Tolerance of GM-PETRA to image blurring was evaluated in simulations and experiments by comparing with the conventional PETRA technique. We performed inner ear imaging of a healthy subject at 7T. GM-PETRA showed significantly less image blurring due to off-resonance and fast T2* signal decay compared to PETRA. In in vivo imaging, GM-PETRA nicely captured complex structures of the inner ear such as the cochlea and semicircular canals. Gradient modulation can improve the PETRA image quality and mitigate SAR and RF peak power limitations without special hardware modification in clinical scanners.

Compressed sensing reconstruction of cardiac cine MRI using golden angle spiral trajectories

In dynamic cardiac cine Magnetic Resonance Imaging (MRI), the spatiotemporal resolution is limited by the low imaging speed. Compressed sensing (CS) theory has been applied to improve the imaging speed and thus the spatiotemporal resolution. The purpose of this paper is to improve CS reconstruction of under sampled data by exploiting spatiotemporal sparsity and efficient spiral trajectories. We extend k-t sparse algorithm to spiral trajectories to achieve high spatio temporal resolutions in cardiac cine imaging. We have exploited spatiotemporal sparsity of cardiac cine MRI by applying a 2D+time wavelet-Fourier transform. For efficient coverage of k-space, we have used a modified version of multi shot (interleaved) spirals trajectories. In order to reduce incoherent aliasing artifact, we use different random undersampling pattern for each temporal frame. Finally, we have used nonuniform fast Fourier transform (NUFFT) algorithm to reconstruct the image from the non-uniformly acquired samples. The proposed approach was tested in simulated and cardiac cine MRI data. Results show that higher acceleration factors with improved image quality can be obtained with the proposed approach in comparison to the existing state-of-the-art method. The flexibility of the introduced method should allow it to be used not only for the challenging case of cardiac imaging, but also for other patient motion where the patient moves or breathes during acquisition.

Childhood extracranial neoplasms: the role of imaging in drug development and clinical trials

Cancer is the leading cause of death in children older than 1 year of age and new drugs are necessary to improve outcomes. Imaging is crucial to the drug development process and assessment of therapeutic response. In adults, tumours are often assessed with CT using size criteria. Unfortunately, techniques established in adults are not necessarily applicable in children due to differing pathophysiology, ability to cooperate and increased susceptibility to ionising radiation. MRI, in particular quantitative MRI, has to date not been fully utilised in children with extracranial neoplasms. The specific challenges of imaging in children, the potential for functional imaging techniques to inform upon and their inclusion in clinical trials are discussed.

Estimating liver perfusion from free-breathing continuously acquired dynamic gadolinium-ethoxybenzyl-diethylenetriamine pentaacetic acid-enhanced acquisition with compressed sensing reconstruction

OBJECTIVE

The purpose of this study was to estimate perfusion metrics in healthy and cirrhotic liver with pharmacokinetic modeling of high-temporal resolution reconstruction of continuously acquired free-breathing gadolinium-ethoxybenzyl-diethylenetriamine pentaacetic acid-enhanced acquisition in patients undergoing clinically indicated liver magnetic resonance imaging.

SUBJECTS AND METHODS

In this Health Insurance Portability and Accountability Act-compliant prospective study, 9 cirrhotic and 10 noncirrhotic patients underwent clinical magnetic resonance imaging, which included continuously acquired radial stack-of-stars 3-dimensional gradient recalled echo sequence with golden-angle ordering scheme in free breathing during contrast injection. A total of 1904 radial spokes were acquired continuously in 318 to 340 seconds. High-temporal resolution data sets were formed by grouping 13 spokes per frame for temporal resolution of 2.2 to 2.4 seconds, which were reconstructed using the golden-angle radial sparse parallel technique that combines compressed sensing and parallel imaging. High-temporal resolution reconstructions were evaluated by a board-certified radiologist to generate gadolinium concentration-time curves in the aorta (arterial input function), portal vein (venous input function), and liver, which were fitted to dual-input dual-compartment model to estimate liver perfusion metrics that were compared between cirrhotic and noncirrhotic livers.

RESULTS

The cirrhotic livers had significantly lower total plasma flow (70.1 ± 10.1 versus 103.1 ± 24.3 mL/min per 100 mL; P < 0.05), lower portal venous flow (33.4 ± 17.7 versus 89.9 ± 20.8 mL/min per 100 mL; P < 0.05), and higher arterial perfusion fraction (52.0% ± 23.4% versus 12.4% ± 7.1%; P < 0.05). The mean transit time was higher in the cirrhotic livers (24.4 ± 4.7 versus 15.7 ± 3.4 seconds; P < 0.05), and the hepatocellular uptake rate was lower (3.03 ± 2.1 versus 6.53 ± 2.4 100/min; P < 0.05).

CONCLUSIONS

Liver perfusion metrics can be estimated from free-breathing dynamic acquisition performed for every clinical examination without additional contrast injection or time. This is a novel paradigm for dynamic liver imaging.

Anisotropic field-of-view support for golden angle radial imaging

PURPOSE

To provide anisotropic field-of-view (FOV) support for golden angle radial imaging.

THEORY AND METHODS

In radial imaging, uniform spoke density leads to a circular FOV, which is excessive for objects with anisotropic dimensions. Larson et al previously showed that the angular k-space spoke density can be determined by the desired anisotropic FOV. We show that conventional golden angle sampling can be deployed in an angle-normalized space and transformed back to k-space such that the desired nonuniform spoke density is preserved for arbitrary temporal window length. Elliptical FOVs were used to illustrate this generalized mapping approach. Point-spread-function and spoke density analysis was performed. Phantom and in vivo cardiac images were acquired.

RESULTS

Simulations, phantom, and in vivo experiments confirmed that the proposed method is able to achieve anisotropic FOV while still maintaining the benefits of golden angle sampling. This approach requires 50% less spokes for elliptical FOV with major-to-minor-axis ratio of 1:0.3, when compared with isotropic FOV with the same undersampling factor.

CONCLUSION

We demonstrate a simple method for applying golden angle view ordering to anisotropic FOV radial imaging. This can reduce imaging time for objects with anisotropic dimensions while still allowing arbitrary temporal window selection. Magn Reson Med 76:229-236, 2016. © 2015 Wiley Periodicals, Inc.