Comparison of data-driven and general temporal constraints on compressed sensing for breast DCE MRI.

PURPOSE: Current breast DCE-MRI strategies provide high sensitivity for cancer detection but are known to be insufficient in fully capturing rapidly changing contrast kinetics at high spatial resolution across both breasts. Advanced acquisition and reconstruction strategies aim to improve spatial and temporal resolution and increase specificity for disease characterization. In this work, we evaluate the spatial and temporal fidelity of a modified data-driven low-rank-based model (known as MOCCO, model consistency condition) compressed-sensing (CS) reconstruction compared to CS with temporal total variation with radial acquisition for high spatial-temporal breast DCE MRI. METHODS: Reconstruction performance was characterized using numerical simulations of a golden-angle stack-of-stars breast DCE-MRI acquisition at 5-second temporal resolution. Specifically, MOCCO was compared with CS total variation and conventional SENSE reconstructions. The temporal model for MOCCO was prelearned over the source data, whereas CS total variation was performed using a first-order temporal gradient sparsity transform. RESULTS: The MOCCO reconstruction was able to capture rapid lesion kinetics while providing high image quality across a range of optimal regularization values. It also recovered kinetics in small lesions (1.5 mm) in line-profile analysis and error images, whereas g-factor maps showed relatively low and constant values with no significant artifacts. The CS-TV method demonstrated either recovery of high spatial resolution with reduced temporal accuracy using large regularization values, or recovery of rapid lesion kinetics with reduced image quality using low regularization values. CONCLUSION: Simulations demonstrated that MOCCO with radial acquisition provides a robust imaging technique for improving temporal fidelity, while maintaining high spatial resolution and image quality in the setting of bilateral breast DCE MRI.

Utility of Stack-of-stars Acquisition for Hepatobiliary Phase Imaging without Breath-holding.

PURPOSE: Post-contrast liver magnetic resonance imaging is typically performed with breath-hold 3D gradient echo sequences. However, breath-holding for >10 s is difficult for some patients. In this study, we compared the quality of hepatobiliary phase (HBP) imaging without breath-holding using the prototype pulse sequences stack-of-stars liver acquisition with volume acceleration (LAVA) (LAVA Star) with or without navigator echoes (LAVA Starnavi+ and LAVA Starnavi-) and Cartesian LAVA with navigator echoes (Cartesian LAVAnavi+). METHODS: Seventy-two patients were included in this single-center, retrospective, cross-sectional study. HBP imaging using the three LAVA sequences (Cartesian LAVAnavi+, LAVA Starnavi-, and LAVA Starnavi+) without breath-holding was performed for all patients using a 3T magnetic resonance system. Two independent radiologists qualitatively analyzed (overall image quality, liver edge sharpness, hepatic vein clarity, streak artifacts, and respiratory motion/pulsation artifacts) HBP images taken by the three sequences using a five-point scale. Quantitative evaluations were also performed by calculating the liver-to-spleen, -lesion, and -portal vein (PV) signal intensity ratios. The results were compared between the three sequences using the Friedman test. RESULTS: LAVA Starnavi+ showed the best image quality and hepatic vein clarity (P < 0.0001). LAVA Starnavi- showed the lowest image quality (P < 0.0001-0.0106). LAVA Starnavi+ images showed fewer streak artifacts than LAVA Starnavi- images (P < 0.0001), while Cartesian LAVAnavi+ images showed no streak artifacts. Cartesian LAVAnavi+ images showed stronger respiratory motion/pulsation artifacts than the others (P < 0.0001). LAVA Starnavi- images showed the highest liver-to-spleen ratios (P < 0.0001-0.0005). Cartesian LAVAnavi+ images showed the lowest liver-to-lesion and -PV ratios (P < 0.0001-0.0108). CONCLUSION: In terms of image quality, the combination of stack-of-stars acquisition and navigator echoes is the best for HBP imaging without breath-holding.

Undersampled linogram trajectory for fast imaging (ULTI): experiments at 3 T and 7 T.

In this study, the performance of linogram acquisition was investigated for the reconstruction of images from undersampled data using parallel imaging methods. The point spread function (PSF) of linogram sampling was analyzed for image sharpness and artifacts. Generalized auto-calibrating partially parallel acquisition was implemented for this new sampling scheme, and images were reconstructed with high acceleration rates. The results were compared with conventional radial sampling methods using simulations and phantom experiments at 3 T. Additionally, a human volunteer was scanned at 7 T. The results demonstrated that the PSF was sharper and the mean artifact power was lower for linogram sampling compared with radial sampling. Results of simulations and phantom experiments were in accord with the findings of the PSF analysis. In simulations, errors in the reconstructed images were lower for linogram sampling. In phantom experiments, fine details and sharp edges were preserved for linogram sampling, while details were blurred for radial sampling. The in vivo human study demonstrated that linogram sampling could provide high quality images of anatomy, even at high acceleration rates. Linogram sampling not only possesses the advantages of radial sampling, such as reduced sensitivity to motion and higher acceleration rates, but it also provides sharper images with fewer artifacts. Moreover, it is less prone to off-resonance artifacts compared with radial sampling. Copyright © 2016 John Wiley & Sons, Ltd.

Wavelet Domain Radiofrequency Pulse Design Applied to Magnetic Resonance Imaging.

A new method for designing radiofrequency (RF) pulses with numerical optimization in the wavelet domain is presented. Numerical optimization may yield solutions that might otherwise have not been discovered with analytic techniques alone. Further, processing in the wavelet domain reduces the number of unknowns through compression properties inherent in wavelet transforms, providing a more tractable optimization problem. This algorithm is demonstrated with simultaneous multi-slice (SMS) spin echo refocusing pulses because reduced peak RF power is necessary for SMS diffusion imaging with high acceleration factors. An iterative, nonlinear, constrained numerical minimization algorithm was developed to generate an optimized RF pulse waveform. Wavelet domain coefficients were modulated while iteratively running a Bloch equation simulator to generate the intermediate slice profile of the net magnetization. The algorithm minimizes the L2-norm of the slice profile with additional terms to penalize rejection band ripple and maximize the net transverse magnetization across each slice. Simulations and human brain imaging were used to demonstrate a new RF pulse design that yields an optimized slice profile and reduced peak energy deposition when applied to a multiband single-shot echo planar diffusion acquisition. This method may be used to optimize factors such as magnitude and phase spectral profiles and peak RF pulse power for multiband simultaneous multi-slice (SMS) acquisitions. Wavelet-based RF pulse optimization provides a useful design method to achieve a pulse waveform with beneficial amplitude reduction while preserving appropriate magnetization response for magnetic resonance imaging.

Quantitative analysis of the efficacy of gradient table correction on improving the accuracy of fiber tractography.

PURPOSE: Several groups suggested that the gradient table of a DTI data set should be reoriented to compensate for head motion. Although the effects of this correction were demonstrated qualitatively, its efficacy was not demonstrated quantitatively to date. The main goal of this study was to investigate the efficacy of gradient table correction on improving the accuracy of fiber tractography. METHODS: First, the effects of gradient table correction on the estimation of fractional anisotropy (FA) maps and the primary diffusion direction were quantified and compared with the inherent uncertainty in the estimation process. Then, the effects of gradient table correction on tractography were quantified. RESULTS: The corrections in FA values were only a fraction of the typical values seen in major fasciculi and inter-subject variance. The corrections to the primary diffusion direction were also much smaller than the uncertainty inherent in the estimation of its direction. However, the directional estimates were biased due to head motion and deviated fiber tracking. CONCLUSIONS: Corrections to FA values were negligible and are not expected to affect group comparisons. However, a small but consistent bias was introduced to the estimates of primary diffusion direction, which might affect brain connectivity analyses based on fiber tracking.

Highly accelerated projection imaging with coil sensitivity encoding for rapid MRI.

PURPOSE: Rapid magnetic resonance imaging (MRI) acquisition is typically achieved by acquiring all or most lines of k-space after one radio frequency (RF) excitation. Parallel imaging techniques can further accelerate data acquisition by acquiring fewer phase-encoded lines and utilizing the spatial sensitivity information of the RF coil arrays. The goal of this study was to develop a new MRI data acquisition and reconstruction technique that is capable of reconstructing a two-dimensional (2D) image using highly undersampled k-space data without any special hardware. Such a technique would be very efficient, as it would significantly reduce the time wasted during multiple RF excitations or phase encoding and gradient switching periods. METHODS: The essence of this new technique is to densely sample a small number of projections, which should be acquired at an angle other than 0° or multiples of 45°. This results in multiple rays passing through a voxel and provides new and independent measurements for each voxel. Then the images are reconstructed using the unique information coming from these projections combined with RF coil sensitivity profiles. The feasibility of this new technique was investigated with realistic simulations and experimental studies using a phantom and compared with conventional nonuniform fast Fourier transform technique. Eigenvalue analysis and error calculations were conducted to find optimal projection angles and minimum requirements for dense sampling. RESULTS: Reconstruction of 64 × 64 images was done using a single projection from simulated data under different noise levels. Simulated reconstruction was also tested with two projections to assess the improvement. Experimental phantom images were reconstructed at higher resolution using 4, 8, and 16 projections. Cross-sectional profiles illustrate that the new technique resolved compartment boundaries clearly. CONCLUSIONS: Simulations demonstrated that only a single k-space line might be sufficient to reconstruct a 2D image using this new technique. Experimental results showed that this is a promising new technique for fast imaging. Using the information from the simulations and fast imaging parameters published in the literature, it could be predicted that a two-dimensional image could be acquired in about 10 ms. One of the major advantages of this new technique is that it does not require any additional hardware and can be implemented on a conventional scanner with an eight-channel coil.