An optimal radial profile order based on the Golden Ratio for time-resolved MRI

Oxygen enhanced lung MRI by simultaneous measurement of T1 and T2 * during free breathing using ultrashort TE

PURPOSE

To provide a robust method for the simultaneous quantification of T1 and T2 * in the human lung during free breathing. Breathing pure oxygen accelerates T1 and T2 * relaxation in the lung. While T1 shortening reflects an increased amount of dissolved molecular oxygen in lung tissue, T2 * shortening shows an increased concentration of oxygen in the alveolar gas. Therefore, both parameters reflect different aspects of the oxygen uptake and provide complementary lung functional information.

MATERIALS AND METHODS

A segmented inversion recovery Look-Locker multiecho sequence based on a multiecho 2D ultrashort TE (UTE) was employed for simultaneous T1 and T2 * quantification. The radial projections follow a modified golden angle ordering, allowing for respiratory self-gating and thus the reconstruction of a series of differently T1 and T2 *-weighted images in arbitrary breathing states. The method was evaluated in nine healthy volunteers while breathing room air and pure oxygen, with two volunteers examined at five oxygen concentrations.

RESULTS

Relative differences of ΔT1 between 7.9% and 12.7% and of ΔT2 * between 13.2% and 6.0% were found.

CONCLUSION

The proposed method provides inherently coregistered, quantitative T1 and T2 * maps in both expiration and inspiration from a single measurement acquired during free breathing and is thus well suited for clinical application.

Generating multiple contrasts using single-shot radial T1 sensitive and insensitive steady-state imaging

PURPOSE

Recently, the (Resolution Enhanced-) T1 insensitive steady-state imaging (TOSSI) approach has been proposed for the fast acquisition of T2 -weighted images. This has been achieved by balanced steady-state free precession (bSSFP) imaging between unequally spaced inversion pulses. The purpose of this work is to present an extension of this technique, considerably increasing both the efficiency and possibilities of TOSSI.

THEORY AND METHODS

A radial trajectory in combination with an appropriate view-sharing reconstruction is used. Because each projection traverses the contrast defining k-space center, several different contrasts can be extracted from a single-shot measurement. These contrasts include various T2 -weightings and T2 /T1 -weighting if an even number of inversion pulses is used, while an odd number allow the generation of several images with predefined tissue types cancelled.

RESULTS

The approach is validated for brain and abdominal imaging at 3.0 Tesla. Results are compared with RE-TOSSI, bSSFP, and turbo spin-echo images and are shown to provide similar contrasts in a fraction of scan time. Furthermore, the potential utility of the approach is illustrated by images obtained from a brain tumor patient.

CONCLUSION

Radial T1 sensitive and insensitive steady-state imaging is able to generate multiple contrasts out of one single-shot measurement in a short scan time.

Consideration of slice profiles in inversion recovery Look-Locker relaxation parameter mapping

PURPOSE

To include the flip angle distribution caused by the slice profile into the model used for describing the relaxation curves observed in inversion recovery Look-Locker FLASH T1 mapping for a more accurate determination of the relaxation parameters.

MATERIALS AND METHODS

For each inversion time, the flip angle dependent signal of the mono-exponential relaxation model is integrated across the slice profile. The resulting Consideration of Slice Profiles (CSP) relaxation curves are compared to the mono-exponential signal model in numerical simulations as well as in phantom and in-vivo experiments.

RESULTS

All measured relaxation curves showed systematic deviations from a mono-exponential curve increasing with flip angle and T1 but decreasing with repetition time. Additionally, the accuracy of T1 was found to be largely dependent on the temporal coverage of the relaxation curve. All these systematic errors were largely reduced by the CSP model.

CONCLUSION

The proposed CSP model represents a useful extension of the conventionally used mono-exponential relaxation model. Despite inherent model inaccuracies, the mono-exponential model was found to be sufficient for many T1 mapping situations. However, if only a poor temporal coverage of the relaxation process is achievable or a very precise modeling of the relaxation course is needed as in model-based techniques, the mono-exponential model leads to systematic errors and the CSP model should be used instead.

Ultrafast CEST imaging

We describe a new MR imaging method for the rapid characterization or screening of chemical exchange saturation transfer (CEST) contrast agents. It is based on encoding the chemical shift dimension with an additional gradient as proposed in previous ultrafast CEST spectroscopy approaches, but extends these with imaging capabilities. This allows us to investigate multiple compounds simultaneously with an arbitrary sample tube arrangement. The technique requires a fast multislice readout to ensure the saturation is not lost during data acquisition due to T1 relaxation. We therefore employ radial subsampling, acquiring only 10 projections per CEST image with a 128×128 matrix. To recover the images, we use a heuristic reconstruction algorithm that incorporates low rank and limited object support as prior knowledge. This way, we are able to acquire a spectral CEST data set consisting of 15 saturation offsets more than 16 times faster than compared with conventional CEST imaging.

Speech MRI: morphology and function

Magnetic Resonance Imaging (MRI) plays an increasing role in the study of speech. This article reviews the MRI literature of anatomical imaging, imaging for acoustic modelling and dynamic imaging. It describes existing imaging techniques attempting to meet the challenges of imaging the upper airway during speech and examines the remaining hurdles and future research directions.

Dynamic contrast-enhanced MRI of the prostate with high spatiotemporal resolution using compressed sensing, parallel imaging, and continuous golden-angle radial sampling: preliminary experience

PURPOSE

To demonstrate dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI) of the prostate with both high spatial and temporal resolution via a combination of golden-angle radial k-space sampling, compressed sensing, and parallel-imaging reconstruction (GRASP), and to compare image quality and lesion depiction between GRASP and conventional DCE in prostate cancer patients.

MATERIALS AND METHODS

Twenty prostate cancer patients underwent two 3T prostate MRI examinations on separate dates, one using standard DCE (spatial resolution 3.0 × 1.9 × 1.9 mm, temporal resolution 5.5 sec) and the other using GRASP (spatial resolution 3.0 × 1.1 × 1.1 mm, temporal resolution 2.3 sec). Two radiologists assessed measures of image quality and dominant lesion size. The experienced reader recorded differences in contrast arrival times between the dominant lesion and benign prostate.

RESULTS

Compared with standard DCE, GRASP demonstrated significantly better clarity of the capsule, peripheral/transition zone boundary, urethra, and periprostatic vessels; image sharpness; and lesion conspicuity for both readers (P < 0.001-0.020). GRASP showed improved interreader correlation for lesion size (GRASP: r = 0.691-0.824, standard: r = 0.495-0.542). In 8/20 cases, only GRASP showed earlier contrast arrival in tumor than benign; in no case did only standard DCE show earlier contrast arrival in tumor.

CONCLUSION

High spatiotemporal resolution prostate DCE is possible with GRASP, which has the potential to improve image quality and lesion depiction as compared with standard DCE.

Fast, free-breathing, in vivo fetal imaging using time-resolved 3D MRI technique: preliminary results

Fetal MR imaging is very challenging due to the movement of fetus and the breathing motion of the mother. Current clinical protocols involve quick 2D scouting scans to determine scan plane and often several attempts to reorient the scan plane when the fetus moves. This makes acquisition of fetal MR images clinically challenging and results in long scan times in order to obtain images that are of diagnostic quality. Compared to 2D imaging, 3D imaging of the fetus has many advantages such as higher SNR and ability to reformat images in multiple planes. However, it is more sensitive to motion and challenging for fetal imaging due to irregular fetal motion in addition to maternal breathing and cardiac motion. This aim of this study is to develop a fast 3D fetal imaging technique to resolve the challenge of imaging the moving fetus. This 3D imaging sequence has multi-echo radial sampling in-plane and conventional Cartesian encoding through plane, which provides motion robustness and high data acquisition efficiency. The utilization of a golden-ratio based projection profile allows flexible time-resolved image reconstruction with arbitrary temporal resolution at arbitrary time points as well as high signal-to-noise and contrast-to-noise ratio. The nice features of the developed image technique allow the 3D visualization of the movements occurring throughout the scan. In this study, we applied this technique to three human subjects for fetal MRI and achieved promising preliminary results of fetal brain, heart and lung imaging.

Multimodality multiparametric imaging of early tumor response to a novel antiangiogenic therapy based on anticalins

Anticalins are a novel class of targeted protein therapeutics. The PEGylated Anticalin Angiocal (PRS-050-PEG40) is directed against VEGF-A. The purpose of our study was to compare the performance of diffusion weighted imaging (DWI), dynamic contrast enhanced magnetic resonance imaging (DCE)-MRI and positron emission tomography with the tracer [18F]fluorodeoxyglucose (FDG-PET) for monitoring early response to antiangiogenic therapy with PRS-050-PEG40. 31 mice were implanted subcutaneously with A673 rhabdomyosarcoma xenografts and underwent DWI, DCE-MRI and FDG-PET before and 2 days after i.p. injection of PRS-050-PEG40 (n = 13), Avastin (n = 6) or PBS (n = 12). Tumor size was measured manually with a caliper. Imaging results were correlated with histopathology. In the results, the tumor size was not significantly different in the treatment groups when compared to the control group on day 2 after therapy onset (P = 0.09). In contrast the imaging modalities DWI, DCE-MRI and FDG-PET showed significant differences between the therapeutic compared to the control group as early as 2 days after therapy onset (P<0.001). There was a strong correlation of the early changes in DWI, DCE-MRI and FDG-PET at day 2 after therapy onset and the change in tumor size at the end of therapy (r = -0.58, 0.71 and 0.67 respectively). The imaging results were confirmed by histopathology, showing early necrosis and necroptosis in the tumors. Thus multimodality multiparametric imaging was able to predict therapeutic success of PRS-050-PEG40 and Avastin as early as 2 days after onset of therapy and thus promising for monitoring early response of antiangiogenic therapy.

Using self-consistency for an iterative trajectory adjustment (SCITA)

PURPOSE

To iteratively correct for deviations in radial trajectories with no need of additionally performed calibration scans.

THEORY AND METHODS

Radially acquired data sets-even when undersampled to a certain extend-inherently feature an oversampled area in the center of k-space. Thus, for a perfectly measured trajectory and neglecting noise, information is consistent between multiple measurements gridded to the same Cartesian position within this region. In the case of erroneous coordinates, this accordance-and therefore a correction of the trajectory-can be enforced by an algorithm iteratively shifting the projections with respect to each other by applying the GRAPPA operator. The method was validated in numerical simulations, as well as in radial acquisitions of a phantom and in vivo images at 3T. The results of the correction were compared to a previously proposed correction method.

RESULTS

The newly introduced technique allowed for a reliable trajectory correction in each of the presented examples. The method was able to remove artifacts as effectively as methods that are based on data from additional calibration scans.

CONCLUSION

The iterative technique introduced in this paper allows for a correction of trajectory errors in radial imaging with no need for additional calibration data.

Quantification of cerebral perfusion using dynamic quantitative susceptibility mapping

PURPOSE

The purpose of this study is to develop a dynamic quantitative susceptibility mapping (QSM) technique with sufficient temporal resolution to map contrast agent concentration in cerebral perfusion imaging.

METHODS

The dynamic QSM used a multiecho three-dimensional (3D) spoiled gradient echo golden angle interleaved spiral sequence during contrast bolus injection. Four-dimensional (4D) space-time resolved magnetic field reconstruction was performed using the temporal resolution acceleration with constrained evolution reconstruction method. Deconvolution of the gadolinium-induced field was performed at each time point with the morphology enabled dipole inversion method to generate a 4D gadolinium concentration map, from which three-dimensional spatial distributions of cerebral blood volume and cerebral blood flow were computed.

RESULTS

Initial in vivo brain imaging demonstrated the feasibility of using dynamic QSM for generating quantitative 4D contrast agent maps and imaging three-dimensional perfusion. The cerebral blood flow obtained with dynamic QSM agreed with that obtained using arterial spin labeling.

CONCLUSION

Dynamic QSM can be used to perform 4D mapping of contrast agent concentration in contrast-enhanced magnetic resonance imaging. The perfusion parameters derived from this 4D contrast agent concentration map were in good agreement with those obtained using arterial spin labeling.

Investigating the limitations of single breath-hold renal artery blood flow measurements using spiral phase contrast MR with R-R interval averaging

PURPOSE

1) To validate an R-R interval averaged golden angle spiral phase contrast magnetic resonance (RAGS PCMR) sequence against conventional cine PCMR for assessment of renal blood flow (RBF) in normal volunteers; and 2) To investigate the effects of motion and heart rate on the accuracy of flow measurements using an in silico simulation.

MATERIALS AND METHODS

In 20 healthy volunteers RAGS (∼6 sec breath-hold) and respiratory-navigated cine (∼5 min) PCMR were performed in both renal arteries to assess RBF. A simulation of RAGS PCMR was used to assess the effect of heart rate (30-105 bpm), vessel expandability (0-150%) and translational motion (x1.0-4.0) on the accuracy of RBF measurements.

RESULTS

There was good agreement between RAGS and cine PCMR in the volunteer study (bias: 0.01 L/min, limits of agreement: -0.04 to +0.06 L/min, P = 0.0001). The simulation demonstrated a positive linear relationship between heart rate and error (r = 0.9894, P < 0.0001), a negative linear relationship between vessel expansion and error (r = -0.9484, P < 0.0001), and a nonlinear, heart rate-dependent relationship between vessel translation and error.

CONCLUSION

We have demonstrated that RAGS PCMR accurately measures RBF in vivo. However, the simulation reveals limitations in this technique at extreme heart rates (<40 bpm, >100 bpm), or when there is significant motion (vessel expandability: >80%, vessel translation: >x2.2).

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

PURPOSE

This study proposes and evaluates a novel method for generating efficient undersampling patterns for 3D Cartesian acquisition with compressed sensing (CS) and parallel imaging (PI).

METHODS

Image quality achieved with schemes that accelerate data acquisition, including CS and PI, are sensitive to the design of the specific undersampling scheme used. Ideally random sampling is required to recover MR images from undersampled data with CS. In practice, pseudo-random sampling schemes are usually applied. Radial or spiral sampling either for Cartesian or non-Cartesian acquisitions has been using because of its favorable features such as interleaving flexibility. In this study, we propose to undersample data on the ky-kz plane of the 3D Cartesian acquisition by circularly selecting sampling points in a way that maintains the features of both random and radial or spiral sampling.

RESULTS

The proposed sampling scheme is shown to outperform conventional random and radial or spiral samplings for 3D Cartesian acquisition and is found to be comparable to advanced variable-density Poisson-Disk sampling (vPDS) while retaining interleaving flexibility for dynamic imaging, based on the results with retrospective undersampling. Our preliminary results with the prospective implementation of the proposed undersampling strategy demonstrated its favorable features.

CONCLUSIONS

The proposed undersampling patterns for 3D Cartesian acquisition possess the desirable properties of randomization and radial or spiral trajectories. It provides easy implementation, flexible sampling, and high accuracy of image reconstruction with CS and PI.

Real-time magnetic resonance imaging technique for determining left ventricle pressure-volume loops

BACKGROUND

Rapid determination of the left ventricular (LV) pressure-volume (PV) relationship as loading conditions are varied is the gold standard for assessment of LV function. Cine magnetic resonance imaging (MRI) does not have sufficient spatiotemporal resolution to assess beat-to-beat changes of the LV PV relationship required to measure the LV end-systolic elastance (EES) or preload-recruitable stroke work (PRSW). Our aim was to investigate real-time MRI and semiautomated LV measurement of LV volume to measure PV relations in large animals under normal and inotropically stressed physiologic conditions.

METHODS

We determined that PV relationships could be accurately measured using an image exposure time Tex less than 100 ms and frame rate Tfr less than 50 ms at elevated heart rates (∼140 beats per minute) using a golden angle radial MRI k-space trajectory and active contour segmentation.

RESULTS

With an optimized exposure time (Tex=95 ms and frame rate Tfr=2.8 ms), we found that there was no significant difference between cine and real-time MRI at rest in end-diastolic volume, end-systolic volume, ejection fraction, stroke volume, or cardiac output (n=5, p<0.05) at either normal or elevated heart rates. We found EES increased from 1.9±0.7 to 3.1±0.3 mm Hg/mL and PRSW increased from 6.2±1.2 to 9.1±0.9 mm Hg during continuous intravenous dobutamine infusion (n=5, p<0.05).

CONCLUSIONS

Real-time MRI can assess LV volumes, EES, and PRSW at baseline and elevated inotropic states.

Highly efficient respiratory motion compensated free-breathing coronary MRA using golden-step Cartesian acquisition

PURPOSE

To develop an efficient 3D affine respiratory motion compensation framework for Cartesian whole-heart coronary magnetic resonance angiography (MRA).

MATERIALS AND METHODS

The proposed method achieves 100% scan efficiency by estimating the affine respiratory motion from the data itself and correcting the acquired data in the reconstruction process. For this, a golden-step Cartesian sampling with spiral profile ordering was performed to enable reconstruction of respiratory resolved images at any breathing position and with different respiratory window size. Affine motion parameters were estimated from image-based registration of 3D undersampled respiratory resolved images reconstructed with iterative SENSE and motion correction was performed directly in the reconstruction using a multiple-coils generalized matrix formulation method. This approach was tested on healthy volunteers and compared against a conventional diaphragmatic navigator-gated acquisition using quantitative and qualitative image quality assessment.

RESULTS

The proposed approach achieved 47 ± 12% and 59 ± 6% vessel sharpness for the right (RCA) and left (LAD) coronary arteries, respectively. Also, good quality visual scores of 2.4 ± 0.74 and 2.44 ± 0.86 were observed for the RCA and LAD (scores from 0, no to 4, excellent coronary vessel delineation). A not statically significant difference (P = 0.05) was found between the proposed method and an 8-mm navigator-gated and tracked scan, although scan efficiency increased from 61 ± 10% to 100%.

CONCLUSION

We demonstrate the feasibility of a new 3D affine respiratory motion correction technique for Cartesian whole-heart CMRA that achieves 100% scan efficiency and therefore a predictable acquisition time. This approach yields image quality comparable to that of an 8-mm navigator-gated acquisition with lower scan efficiency. Further evaluation of this technique in patients is now warranted to determine its clinical use.

Self-gated radial MRI for respiratory motion compensation on hybrid PET/MR systems

Accurate localization and uptake quantification of lesions in the chest and abdomen using PET imaging is challenging due to the respiratory motion during the exam. The advent of hybrid PET/MR systems offers new ways to compensate for respiratory motion without exposing the patient to additional radiation. The use of self-gated reconstructions of a 3D radial stack-of-stars GRE acquisition is proposed to derive a high-resolution MRI motion model. The self-gating signal is used to perform respiratory binning of the simultaneously acquired PET raw data. Matching mu-maps are generated for every bin, and post-reconstruction registration is performed in order to obtain a motion-compensated PET volume from the individual gates. The proposed method is demonstrated in-vivo for three clinical patients. Motion-corrected reconstructions are compared against ungated and gated PET reconstructions. In all cases, motion-induced blurring of lesions in the liver and lung was substantially reduced, without compromising SNR as it is the case for gated reconstructions.

Augmented Lagrangian with variable splitting for faster non-Cartesian L1-SPIRiT MR image reconstruction

SPIRiT (iterative self-consistent parallel imaging reconstruction), and its sparsity-regularized variant L1-SPIRiT, are compatible with both Cartesian and non-Cartesian magnetic resonance imaging sampling trajectories. However, the non-Cartesian framework is more expensive computationally, involving a nonuniform Fourier transform with a nontrivial Gram matrix. We propose a novel implementation of the regularized reconstruction problem using variable splitting, alternating minimization of the augmented Lagrangian, and careful preconditioning. Our new method based on the alternating direction method of multipliers converges much faster than existing methods because of the preconditioners' heightened effectiveness. We demonstrate such rapid convergence substantially improves image quality for a fixed computation time. Our framework is a step forward towards rapid non-Cartesian L1-SPIRiT reconstructions.

Dynamically phase-cycled radial balanced SSFP imaging for efficient banding removal

PURPOSE

Balanced steady-state free precession (bSSFP) imaging suffers from banding artifacts due to its inherent sensitivity to inhomogeneities in the main magnetic field. These artifacts can be removed by the acquisition of multiple images at different frequency offsets. However, conventional phase-cycling is hindered by a long scan time. The purpose of this work is to present a novel approach for efficient banding removal in bSSFP imaging.

THEORY AND METHODS

To this end, the phase-cycle during a single-shot radial acquisition of an image was dynamically changed. Thus, each projection is acquired with a different frequency offset. Using conventional radial gridding, an artifact-free image can be reconstructed out of this dataset.

RESULTS

The approach is validated at clinical field strength [3.0 Tesla (T)] as well as at ultrahigh field (9.4T). Robust elimination of banding artifacts was obtained for different imaging regions, including brain imaging at ultrahigh field with an in-plane resolution of 0.25 × 0.25 mm(2). Besides banding artifact-free imaging, the applicability of the proposed technique for fat-water separation is demonstrated.

CONCLUSION

Dynamically phase-cycled radial bSSFP has the potential for banding-free bSSFP imaging in a short scan time, in the presence of severe field inhomogeneities and at high resolution.

Non-Cartesian parallel imaging reconstruction

Non-Cartesian parallel imaging has played an important role in reducing data acquisition time in MRI. The use of non-Cartesian trajectories can enable more efficient coverage of k-space, which can be leveraged to reduce scan times. These trajectories can be undersampled to achieve even faster scan times, but the resulting images may contain aliasing artifacts. Just as Cartesian parallel imaging can be used to reconstruct images from undersampled Cartesian data, non-Cartesian parallel imaging methods can mitigate aliasing artifacts by using additional spatial encoding information in the form of the nonhomogeneous sensitivities of multi-coil phased arrays. This review will begin with an overview of non-Cartesian k-space trajectories and their sampling properties, followed by an in-depth discussion of several selected non-Cartesian parallel imaging algorithms. Three representative non-Cartesian parallel imaging methods will be described, including Conjugate Gradient SENSE (CG SENSE), non-Cartesian generalized autocalibrating partially parallel acquisition (GRAPPA), and Iterative Self-Consistent Parallel Imaging Reconstruction (SPIRiT). After a discussion of these three techniques, several potential promising clinical applications of non-Cartesian parallel imaging will be covered.

Validation of a golden angle radial sequence (GOLD) for abdominal T1 mapping during free breathing: demonstrating clinical feasibility for quantifying gastric secretion and emptying

PURPOSE

To validate a magnetic resonance imaging sequence suitable for quantitative assessment of acid suppression by a proton pump inhibitor (PPI) on gastric secretion and emptying in clinical practice.

METHODS

A golden angle radial sequence (GOLD) was validated in a series of in vitro and in vivo experiments and clinical feasibility was shown in two studies. The impact of free breathing and image plane orientation on T1 values was evaluated in a controlled in vivo experiment. The free-breathing GOLD sequence was compared against a standard breath-hold gradient echo sequence for gastric half emptying time in 23 subjects during a gastric emptying study. Pilot data from five subjects assessed the sensitivity of the GOLD sequence to detect changes in acid secretion volume produced by PPI treatment.

RESULTS

The coronal free-breathing GOLD sequence and the axial breath-hold standard gradient echo sequence showed good agreement of the gastric half emptying time (6 ± 3 min, P = 0.053). The GOLD sequence demonstrated sensitivity to reduction of gastric secretion volumes induced by PPI treatment (55 ± 5 mL, P < 0.001).

CONCLUSION

The GOLD sequence allowed for free breathing, multislice, combined imaging and T1 mapping of the stomach content. GOLD presents a promising multipurpose, noninvasive imaging tool for monitoring gastric function in clinical studies.

Retrospective Rigid Motion Correction in k-Space for Segmented Radial MRI

Motion occurring during magnetic resonance imaging acquisition is a major factor of image quality degradation. Self-navigation can help reduce artefacts by estimating motion from the acquired data to enable motion correction. Popular self-navigation techniques rely on the availability of a fully-sampled motion-free reference to register the motion corrupted data with. In the proposed technique, rigid motion parameters are derived using the inherent correlation between radial segments in k-space. The registration is performed exclusively in k-space using the Phase Correlation Method, a popular registration technique in computer vision. Robust and accurate registration has been carried out from radial segments composed of as few as 32 profiles. Successful self-navigation has been performed on 2-D dynamic brain scans corrupted with continuous motion for six volunteers. Retrospective motion correction using the derived self-navigation parameters resulted in significant improvement of image quality compared to the conventional sliding window. This work also demonstrates the benefits of using a bit-reversed ordering scheme to limit undesirable effects specific to retrospective motion correction on radial trajectories. This method provides a fast and efficient mean of measuring rigid motion directly in k-space from dynamic radial data under continuous motion.

Noncontrast enhanced four-dimensional dynamic MRA with golden angle radial acquisition and K-space weighted image contrast (KWIC) reconstruction

PURPOSE

To explore the feasibility of 2D and 3D golden-angle radial acquisition strategies in conjunction with k-space weighted image contrast (KWIC) temporal filtering to achieve noncontrast enhanced dynamic MRA (dMRA) with high spatial resolution, low streaking artifacts and high temporal fidelity.

METHODS

Simulations and in vivo examinations in eight normal volunteers and an arteriovenous malformation patient were carried out. Both 2D and 3D golden angle radial sequences, preceded by spin tagging, were used for dMRA of the brain. The radial dMRA data were temporally filtered using the KWIC strategy and compared with matched standard Cartesian techniques.

RESULTS

The 2D and 3D dynamic MRA image series acquired with the proposed radial techniques demonstrated excellent image quality without discernible temporal blurring compared with standard Cartesian based approaches. The image quality of radial dMRA was equivalent to or higher than that of Cartesian dMRA by visual inspection. A reduction factor of up to 10 and 3 in scan time was achieved for 2D and 3D radial dMRA compared with the Cartesian-based counterparts.

CONCLUSION

The proposed 2D and 3D radial dMRA techniques demonstrated image quality comparable or even superior to those obtained with standard Cartesian methods, but within a fraction of the scan time.

Radial k-t SPIRiT: autocalibrated parallel imaging for generalized phase-contrast MRI

PURPOSE

To extend SPIRiT to additionally exploit temporal correlations for highly accelerated generalized phase-contrast MRI and to compare the performance of the proposed radial k-t SPIRiT method relative to frame-by-frame SPIRiT and radial k-t GRAPPA reconstruction for velocity and turbulence mapping in the aortic arch.

THEORY AND METHODS

Free-breathing navigator-gated two-dimensional radial cine imaging with three-directional multi-point velocity encoding was implemented and fully sampled data were obtained in the aortic arch of healthy volunteers. Velocities were encoded with three different first gradient moments per axis to permit quantification of mean velocity and turbulent kinetic energy. Velocity and turbulent kinetic energy maps from up to 14-fold undersampled data were compared for k-t SPIRiT, frame-by-frame SPIRiT, and k-t GRAPPA relative to the fully sampled reference.

RESULTS

Using k-t SPIRiT, improvements in magnitude and velocity reconstruction accuracy were found. Temporally resolved magnitude profiles revealed a reduction in spatial blurring with k-t SPIRiT compared with frame-by-frame SPIRiT and k-t GRAPPA for all velocity encodings, leading to improved estimates of turbulent kinetic energy.

CONCLUSION

k-t SPIRiT offers improved reconstruction accuracy at high radial undersampling factors and hence facilitates the use of generalized phase-contrast MRI for routine use.

Simultaneous static and cine nonenhanced MR angiography using radial sampling and highly constrained back projection reconstruction

PURPOSE

To describe a pulse sequence for simultaneous static and cine nonenhanced magnetic resonance angiography (NEMRA) of the peripheral arteries.

METHODS

The peripheral arteries of 10 volunteers and 6 patients with peripheral arterial disease (PAD) were imaged with the proposed cine NEMRA sequence on a 1.5 Tesla (T) system. The impact of multi-shot imaging and highly constrained back projection (HYPR) reconstruction was examined. The propagation rate of signal along the length of the arterial tree in the cine nonenhanced MR angiograms was quantified.

RESULTS

The cine NEMRA sequence simultaneously provided a static MR angiogram showing vascular anatomy as well as a cine display of arterial pulse wave propagation along the entire length of the peripheral arteries. Multi-shot cine NEMRA improved temporal resolution and reduced image artifacts. HYPR reconstruction improved image quality when temporal reconstruction footprints shorter than 100 ms were used (P < 0.001). Pulse wave propagation within the arterial tree as displayed by cine NEMRA was slower in patients with PAD than in volunteers.

CONCLUSION

Simultaneous static and cine NEMRA of the peripheral arteries is feasible. Multi-shot acquisition and HYPR reconstruction can be used to improve arterial conspicuity and temporal resolution.

Radial spectroscopic MRI of hyperpolarized [1-(13) C] pyruvate at 7 tesla

PURPOSE

The transient and nonrenewable signal from hyperpolarized metabolites necessitates extensive sequence optimization for encoding spatial, spectral, and dynamic information. In this work, we evaluate the utility of radial single-timepoint and cumulative spectroscopic MRI of hyperpolarized [1-(13) C] pyruvate and its metabolic products at 7 Tesla (T).

METHODS

Simulations of radial echo planar spectroscopic imaging (EPSI) and multiband frequency encoding (MBFE) acquisitions were performed to confirm feasibility and evaluate performance for HP (13) C imaging. Corresponding sequences were implemented on a 7T small-animal MRI system, tested in phantom, and demonstrated in a murine model of anaplastic thyroid cancer.

RESULTS

MBFE provides excellent spectral separation but is susceptible to blurring and T2 * signal loss inherent to using low readout gradients. The higher readout gradients and more flexible spectral encoding for EPSI result in good spatial resolution and spectral separation. Radial acquisition throughout HP signal evolution offers the flexibility for reconstructing spatial maps of mean metabolite distribution and global dynamic time courses of multiple metabolites.

CONCLUSION

Radial EPSI and MBFE acquisitions are well-suited for hyperpolarized (13) C MRI over short and long durations. Advantages to this approach include robustness to nonstationary magnetization, insensitivity to precise acquisition timing, and versatility for reconstructing dynamically acquired spectroscopic data.

Retrospective reconstruction of cardiac cine images from golden-ratio radial MRI using one-dimensional navigators

PURPOSE

To demonstrate radial golden-ratio-based cardiac cine imaging by using interspersed one-dimensional (1D) navigators.

MATERIALS AND METHODS

The 1D navigators were interspersed into the acquisition of radial spokes which were continuously rotated by an angle increment based on the golden-ratio. Performing correlation analysis between the 1D navigator projections, time points corresponding to the same cardiac motion phases were automatically identified and used to combine retrospectively golden-ratio rotated radial spokes from multiple data windows. Data windows were shifted consecutively for dynamic reconstruction of different cardiac motion frames. Experiments were performed during a single breathhold. By artificially reducing the amount of input data, signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) as well as artifact level was evaluated for different breathhold durations.

RESULTS

Analysis of the 1D navigator data provided a detailed correlation function revealing cardiac motion over time. Imaging results were comparable to images reconstructed based on a timely synchronized ECG. Cardiac cine images with a low artifact level and good image quality in terms of SNR and CNR were reconstructed from volunteer data achieving a CNR between the myocardium and the left ventricular cavity of 50 for the longest breathhold duration of 26 s. CNR maintained a value higher than 30 for acquisition times as low as 10 s.

CONCLUSION

Combining radial golden-ratio-based imaging with an intrinsic navigator is a promising and robust method for performing high quality cardiac cine imaging.

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

PURPOSE

To develop a fast and flexible free-breathing dynamic volumetric MRI technique, iterative Golden-angle RAdial Sparse Parallel MRI (iGRASP), that combines compressed sensing, parallel imaging, and golden-angle radial sampling.

METHODS

Radial k-space data are acquired continuously using the golden-angle scheme and sorted into time series by grouping an arbitrary number of consecutive spokes into temporal frames. An iterative reconstruction procedure is then performed on the undersampled time series where joint multicoil sparsity is enforced by applying a total-variation constraint along the temporal dimension. Required coil-sensitivity profiles are obtained from the time-averaged data.

RESULTS

iGRASP achieved higher acceleration capability than either parallel imaging or coil-by-coil compressed sensing alone. It enabled dynamic volumetric imaging with high spatial and temporal resolution for various clinical applications, including free-breathing dynamic contrast-enhanced imaging in the abdomen of both adult and pediatric patients, and in the breast and neck of adult patients.

CONCLUSION

The high performance and flexibility provided by iGRASP can improve clinical studies that require robustness to motion and simultaneous high spatial and temporal resolution. Magn Reson Med 72:707-717, 2014. © 2013 Wiley Periodicals, Inc.

Fast retrospectively triggered local pulse-wave velocity measurements in mice with CMR-microscopy using a radial trajectory

BACKGROUND

The aortic pulse-wave velocity (PWV) is an important indicator of cardiovascular risk. In recent studies MRI methods have been developed to measure this parameter noninvasively in mice. Present techniques require additional hardware for cardiac and respiratory gating. In this work a robust self-gated measurement of the local PWV in mice without the need of triggering probes is proposed.

METHODS

The local PWV of 6-months-old wild-type C57BL/6J mice (n=6) was measured in the abdominal aorta with a retrospectively triggered radial Phase Contrast (PC) MR sequence using the flow-area (QA) method. A navigator signal was extracted from the CMR data of highly asymmetric radial projections with short repetition time (TR=3 ms) and post-processed with high-pass and low-pass filters for retrospective cardiac and respiratory gating. The self-gating signal was used for a reconstruction of high-resolution Cine frames of the aortic motion. To assess the local PWV the volume flow Q and the cross-sectional area A of the aorta were determined. The results were compared with the values measured with a triggered Cartesian and an undersampled triggered radial PC-Cine sequence.

RESULTS

In all examined animals a self-gating signal could be extracted and used for retrospective breath-gating and PC-Cine reconstruction. With the non-triggered measurement PWV values of 2.3±0.2 m/s were determined. These values are in agreement with those measured with the triggered Cartesian (2.4±0.2 m/s) and the triggered radial (2.3±0.2 m/s) measurement. Due to the strong robustness of the radial trajectory against undersampling an acceleration of more than two relative to the prospectively triggered Cartesian sampling could be achieved with the retrospective method.

CONCLUSION

With the radial flow-encoding sequence the extraction of a self-gating signal is feasible. The retrospective method enables a robust and fast measurement of the local PWV without the need of additional trigger hardware.

Real-time cardiovascular magnetic resonance at 1.5 T using balanced SSFP and 40 ms resolution

BACKGROUND

While cardiovascular magnetic resonance (CMR) commonly employs ECG-synchronized cine acquisitions with balanced steady-state free precession (SSFP) contrast at 1.5 T, recent developments at 3 T demonstrate significant potential for T1-weighted real-time imaging at high spatiotemporal resolution using undersampled radial FLASH. The purpose of this work was to combine both ideas and to evaluate a corresponding real-time CMR method at 1.5 T with SSFP contrast.

METHODS

Radial gradient-echo sequences with fully balanced gradients and at least 15-fold undersampling were implemented on two CMR systems with different gradient performance. Image reconstruction by regularized nonlinear inversion (NLINV) was performed offline and resulted in real-time SSFP CMR images at a nominal resolution of 1.8 mm and with acquisition times of 40 ms.

RESULTS

Studies of healthy subjects demonstrated technical feasibility in terms of robustness and general image quality. Clinical applicability with access to quantitative evaluations (e.g., ejection fraction) was confirmed by preliminary applications to 27 patients with typical indications for CMR including arrhythmias and abnormal wall motion. Real-time image quality was slightly lower than for cine SSFP recordings, but considered diagnostic in all cases.

CONCLUSIONS

Extending conventional cine approaches, real-time radial SSFP CMR with NLINV reconstruction provides access to individual cardiac cycles and allows for studies of patients with irregular heartbeat.

Dynamic MRI of the TMJ under physical load

OBJECTIVES

The objective of this study was to examine the kinematics of structures of the temporomandibular joint (TMJ) under physiological load while masticating.

METHODS

Radial MRI was chosen as a fast imaging method to dynamically capture the motions of the joint's anatomy. The technique included a golden ratio-based increment angle and a sliding window reconstruction. The measurements were performed on 22 subjects with and without deformation/displacement of the intra-articular disc while they were biting on a cooled caramel toffee.

RESULTS

The reconstructed dynamic images provided sufficient information about the size and localization of the disc as well as the change of the intra-articular distance with and without loading.

CONCLUSIONS

The feasibility of the golden ratio-based radial MRI technique to dynamically capture the anatomy of the TMJ under physical load was demonstrated in this initial study.

Real-time 3D magnetic resonance imaging of the pharyngeal airway in sleep apnea

PURPOSE

To investigate the feasibility of real-time 3D magnetic resonance imaging (MRI) with simultaneous recording of physiological signals for identifying sites of airway obstruction during natural sleep in pediatric patients with sleep-disordered breathing.

METHODS

Experiments were performed using a three-dimensional Fourier transformation (3DFT) gradient echo sequence with prospective undersampling based on golden-angle radial spokes, and L1-norm regularized iterative self-consistent parallel imaging (L1-SPIRiT) reconstruction. This technique was demonstrated in three healthy adult volunteers and five pediatric patients with sleep-disordered breathing. External airway occlusion was used to induce partial collapse of the upper airway on inspiration and test the effectiveness of the proposed imaging method. Apneic events were identified using information available from synchronized recording of mask pressure and respiratory effort.

RESULTS

Acceptable image quality was obtained in seven of eight subjects. Temporary airway collapse induced via inspiratory loading was successfully imaged in all three volunteers, with average airway volume reductions of 63.3%, 52.5%, and 33.7%. Central apneic events and associated airway narrowing/closure were identified in two pediatric patients. During central apneic events, airway obstruction was observed in the retropalatal region in one pediatric patient.

CONCLUSION

Real-time 3D MRI of the pharyngeal airway with synchronized recording of physiological signals is feasible and may provide valuable information about the sites and nature of airway narrowing/collapse during natural sleep.

Free-breathing contrast-enhanced multiphase MRI of the liver using a combination of compressed sensing, parallel imaging, and golden-angle radial sampling

OBJECTIVE

The objectives of this study were to develop a new method for free-breathing contrast-enhanced multiphase liver magnetic resonance imaging (MRI) using a combination of compressed sensing, parallel imaging, and radial k-space sampling and to demonstrate the feasibility of this method by performing image quality comparison with breath-hold cartesian T1-weighted (conventional) postcontrast acquisitions in healthy participants.

MATERIALS AND METHODS

This Health Insurance Portability and Accountability Act-compliant prospective study received approval from the institutional review board. Eight participants underwent 3 separate contrast-enhanced fat-saturated T1-weighted gradient-echo MRI examinations with matching imaging parameters: conventional breath-hold examination with cartesian k-space sampling volumetric interpolate breath hold examination (BH-VIBE) and free-breathing acquisitions with interleaved angle-bisection and continuous golden-angle radial sampling schemes. Interleaved angle-bisection and golden-angle data from each 100 consecutive spokes were reconstructed using a combination of compressed sensing and parallel imaging (interleaved-angle radial sparse parallel [IARASP] and golden-angle radial sparse parallel [GRASP]) to generate multiple postcontrast phases.Arterial- and venous-phase BH-VIBE, IARASP, and GRASP reconstructions were evaluated by 2 radiologists in a blinded fashion. The readers independently assessed quality of enhancement (QE), overall image quality (IQ), and other parameters of image quality on a 5-point scale, with the highest score indicating the most desirable examination. Mixed model analysis of variance was used to compare each measure of image quality.

RESULTS

Images of BH-VIBE and GRASP had significantly higher QE and IQ values compared with IARASP for both phases (P < 0.05). The differences in QE between BH-VIBE and GRASP for the arterial and venous phases were not significant (P > 0.05). Although GRASP had lower IQ score compared with BH-VIBE for the arterial (3.9 vs 4.8; P < 0.0001) and venous (4.2 vs 4.8; P = 0.005) phases, GRASP received IQ scores of 3 or more in all participants, which was consistent with acceptable or better diagnostic image quality.

CONCLUSION

Contrast-enhanced multiphase liver MRI of diagnostic quality can be performed during free breathing using a combination of compressed sensing, parallel imaging, and golden-angle radial sampling.

4D flow MRI

Traditionally, magnetic resonance imaging (MRI) of flow using phase contrast (PC) methods is accomplished using methods that resolve single-directional flow in two spatial dimensions (2D) of an individual slice. More recently, three-dimensional (3D) spatial encoding combined with three-directional velocity-encoded phase contrast MRI (here termed 4D flow MRI) has drawn increased attention. 4D flow MRI offers the ability to measure and to visualize the temporal evolution of complex blood flow patterns within an acquired 3D volume. Various methodological improvements permit the acquisition of 4D flow MRI data encompassing individual vascular structures and entire vascular territories such as the heart, the adjacent aorta, the carotid arteries, abdominal, or peripheral vessels within reasonable scan times. To subsequently analyze the flow data by quantitative means and visualization of complex, three-directional blood flow patterns, various tools have been proposed. This review intends to introduce currently used 4D flow MRI methods, including Cartesian and radial data acquisition, approaches for accelerated data acquisition, cardiac gating, and respiration control. Based on these developments, an overview is provided over the potential this new imaging technique has in different parts of the body from the head to the peripheral arteries.

Self-navigated tissue phase mapping using a golden-angle spiral acquisition-proof of concept in patients with pulmonary hypertension

PURPOSE

To create a high temporal- and spatial-resolution retrospectively cardiac-gated, tissue phase mapping (TPM) sequence, using an image-based respiratory navigator calculated from the data itself.

METHODS

The sequence was based on a golden-angle spiral acquisition. Reconstruction of real-time images allowed creation of an image-based navigator. The expiratory spiral interleaves were then retrospectively cardiac-gated using data binning. TPM data were acquired in 20 healthy volunteers and 10 patients with pulmonary hypertension. Longitudinal and radial myocardial velocities were calculated in the left ventricle and right ventricle.

RESULTS

The image-based navigator was shown to correlate well with simultaneously acquired airflow data in 10 volunteers(r=0.93±0.04). The TPM navigated images had a significantly higher subjective image quality and edge sharpness (P<0.0001) than averaged spiral TPM. No significant differences in myocardial velocities were seen between conventional Cartesian TPM with navigator respiratory-gating and the proposed self-navigated TPM technique, in 10 volunteers. Significant differences in the velocities were seen between the volunteers and patients in the left ventricle at systole and end diastole and in the right ventricle at end diastole.

CONCLUSION

The feasibility of measuring myocardial motion using a golden-angle spiral TPM sequence was demonstrated, with an image-based respiratory navigator calculated from the TPM data itself.

MRI temporal acceleration techniques

In recent years, there has been an explosive growth of magnetic resonance imaging (MRI) techniques that allow faster scan speed by exploiting temporal or spatiotemporal redundancy of the images. These techniques improve the performance of dynamic imaging significantly across multiple clinical applications, including cardiac functional examinations, perfusion imaging, blood flow assessment, contrast-enhanced angiography, functional MRI, and interventional imaging, among others. The scan acceleration permits higher spatial resolution, increased temporal resolution, shorter scan duration, or a combination of these benefits. Along with the exciting developments is a dizzying proliferation of acronyms and variations of the techniques. The present review attempts to summarize this rapidly growing topic and presents conceptual frameworks to understand these techniques in terms of their underlying mechanics and connections. Techniques from view sharing, keyhole, k-t, to compressed sensing are covered.

Model-based Acceleration of Parameter mapping (MAP) for saturation prepared radially acquired data

A reconstruction technique called Model-based Acceleration of Parameter mapping (MAP) is presented allowing for quantification of longitudinal relaxation time and proton density from radial single-shot measurements after saturation recovery magnetization preparation. Using a mono-exponential model in image space, an iterative fitting algorithm is used to reconstruct one well resolved and consistent image for each of the projections acquired during the saturation recovery relaxation process. The functionality of the algorithm is examined in numerical simulations, phantom experiments, and in-vivo studies. MAP reconstructions of single-shot acquisitions feature the same image quality and resolution as fully sampled reference images in phantom and in-vivo studies. The longitudinal relaxation times obtained from the MAP reconstructions are in very good agreement with the reference values in numerical simulations as well as phantom and in-vivo measurements. Compared to available contrast manipulation techniques, no averaging of projections acquired at different time points of the relaxation process is required in MAP imaging. The proposed technique offers new ways of extracting quantitative information from single-shot measurements acquired after magnetization preparation. The reconstruction simultaneously yields images with high spatiotemporal resolution fully consistent with the acquired data as well as maps of the effective longitudinal relaxation parameter and the relative proton density.

Two-dimensional radial sodium heart MRI using variable-rate selective excitation and retrospective electrocardiogram gating with golden angle increments

PURPOSE

Two-dimensional projection reconstruction methods provide advantages over three-dimensional techniques because of higher flexibility regarding the resolution and shorter scan time needed. To optimize a two-dimensional radial sequence with respect to signal-to-noise ratio, variable-rate selective excitation and retrospective electrocardiogram gating is investigated.

METHODS

The minimal radiofrequency pulse duration is simulated in dependence of the flip angle and coil parameters using sinc waveforms with two different variable-rate selective excitation approaches and a Fermi pulse. Retrospectively electrocardiogram-gated imaging with Golden Angle incremented projections was implemented to allow for continuous data acquisition enabling the possibility of dynamic electrocardiogram-gated heart imaging.

RESULTS

Especially for abdominal coils with high transmitter voltages required, variable-rate selective excitation strongly reduces the radiofrequency pulse duration and echo time resulting in a signal-to-noise ratio gain up to 15.5% (if the fast relaxation component of sodium is in the order of the radiofrequency pulse duration) compared with standard sinc-shaped radiofrequency pulses. Retrospective electrocardiogram gating shows higher flexibility with regard to the trigger delay enabling the trade-off between heart motion artifacts and signal-to-noise ratio.

CONCLUSION

A two-dimensional radial sequence is optimized for sodium heart imaging regarding signal-to-noise ratio. Different sodium contrasts of the human heart are shown, which can give additional information on heart diseases.

Distributed spirals: a new class of three-dimensional k-space trajectories

This work presents a new class of three-dimensional spiral based-trajectories for sampling magnetic resonance data. The distributed spirals trajectory efficiently traverses a cylinder or sphere or intermediate shape in k-space. The trajectory is shown to be nearly as efficient as a conventional stack of spirals trajectory in terms of scan time and signal-to-noise ratio, while reducing coherent aliasing in all three spatial directions and reducing Gibbs ringing due to the nature of collecting data from a sphere in k-space. The trajectory uses a single two-dimensional spiral waveform with the addition of a single orthogonal waveform which is scaled with each repetition, making it relatively easy to implement. Blurring from off-resonance only occurs in two dimensions due to the temporal nature of the sampling.

Temporal and spectral imaging with micro-CT

PURPOSE

Micro-CT is widely used for small animal imaging in preclinical studies of cardiopulmonary disease, but further development is needed to improve spatial resolution, temporal resolution, and material contrast. We present a technique for visualizing the changing distribution of iodine in the cardiac cycle with dual source micro-CT.

METHODS

The approach entails a retrospectively gated dual energy scan with optimized filters and voltages, and a series of computational operations to reconstruct the data. Projection interpolation and five-dimensional bilateral filtration (three spatial dimensions + time + energy) are used to reduce noise and artifacts associated with retrospective gating. We reconstruct separate volumes corresponding to different cardiac phases and apply a linear transformation to decompose these volumes into components representing concentrations of water and iodine. Since the resulting material images are still compromised by noise, we improve their quality in an iterative process that minimizes the discrepancy between the original acquired projections and the projections predicted by the reconstructed volumes. The values in the voxels of each of the reconstructed volumes represent the coefficients of linear combinations of basis functions over time and energy. We have implemented the reconstruction algorithm on a graphics processing unit (GPU) with CUDA. We tested the utility of the technique in simulations and applied the technique in an in vivo scan of a C57BL∕6 mouse injected with blood pool contrast agent at a dose of 0.01 ml∕g body weight. Postreconstruction, at each cardiac phase in the iodine images, we segmented the left ventricle and computed its volume. Using the maximum and minimum volumes in the left ventricle, we calculated the stroke volume, the ejection fraction, and the cardiac output.

RESULTS

Our proposed method produces five-dimensional volumetric images that distinguish different materials at different points in time, and can be used to segment regions containing iodinated blood and compute measures of cardiac function.

CONCLUSIONS

We believe this combined spectral and temporal imaging technique will be useful for future studies of cardiopulmonary disease in small animals.