Dynamic MRI with projection reconstruction and KWIC processing for simultaneous high spatial and temporal resolution.

A method for dynamic imaging in MRI is presented that enables the acquisition of a series of images with both high temporal and high spatial resolution. The technique, which is based on the projection reconstruction (PR) imaging scheme, utilizes distinct data acquisition and reconstruction strategies to achieve this simultaneous capability. First, during acquisition, data are collected in multiple undersampled passes, with the view angles interleaved in such a way that those of subsequent passes bisect the views of earlier ones. During reconstruction, these views are weighted according to a previously described k-space weighted image contrast (KWIC) technique that enables the manipulation of image contrast by selective filtering. Unlike conventional undersampled PR methods, the proposed dynamic KWIC technique does not suffer from low image SNR or image degradation due to streaking artifacts. The effectiveness of dynamic KWIC is demonstrated in both simulations and in vivo, high-resolution, contrast-enhanced imaging of breast lesions.     

Continuous radial data acquisition for dynamic MRI

Since image acquisition times in MRI have been reduced considerably over recent years, several new important application areas of MRI have appeared. In addition to pure static anatomic information, the evolution of a dynamic process may be visualized by a sequence of temporal snapshots of the process acquired within a short time period. This makes applications like interactive or interventional MRI as well as the acquisition of additional functional information feasible. For high temporal resolution, all these applications require a quasi real-time image acquisition during the time the interaction or dynamic process evolves. We present an approach to realtime imaging using a continuous radial acquisition scheme. The intrinsic advantages of radial or projection reconstruction (PR) techniques are used to minimize motion-related image distortions. Modifications of the acquisition scheme as well as dedicated reconstruction techniques are used to further reduce the temporal blurring due to the finite acquisition time of one entire data set in our approach. So far we have used this technique for the visualization of active joint motion.     

Improvement of the diffusion-weighted images acquired with radial trajectories using projection data regeneration

PURPOSE: To reduce the artifacts due to pulsatile motion artifacts in diffusion-weighted imaging (DWI) with radial trajectories and to improve the image quality using projection data regeneration. MATERIALS AND METHODS: The projection data is obtained by a radial spin-echo DWI (rSE-DWI) sequence, from which a temporary image is generated using the inverse Radon transform (IRT). The projection data include some degraded data due to cardiac pulsatile motion. The degraded data are then replaced with data that are regenerated using the Radon transform (RT) of the temporary image. The proposed method for image quality improvement is demonstrated through a computer simulation and in vivo images obtained by rSE-DWI. RESULTS: In general, electrocardiograph (ECG) triggering is used to reduce the degradation of projection data in the DWI with radial trajectories, where the amount of degradation depends on the cardiac phase. Cardiac gating reduces the artifacts resulting from the cardiac pulsatile motion to a certain extent only. The proposed projection data regeneration method successfully improves image quality. CONCLUSION: The regeneration method based on back-projection reconstruction effectively uses the features of the degraded projections having lower signal intensity than the normal projections, resulting in image quality improvement without acquisition of additional data.     

The influence of radial undersampling schemes on compressed sensing reconstruction in breast MRI

Fast imaging applications in magnetic resonance imaging (MRI) frequently involve undersampling of k‐space data to achieve the desired temporal resolution. However, high temporal resolution images generated from undersampled data suffer from aliasing artifacts. In radial k‐space sampling, this manifests as undesirable streaks that obscure image detail. Compressed sensing reconstruction has been shown to reduce such streak artifacts, based on the assumption of image sparsity. Here, compressed sensing is implemented with three different radial sampling schemes (golden‐angle, bit‐reversed, and random sampling), which are compared over a range of spatiotemporal resolutions. The sampling methods are implemented in static scenarios where different undersampling patterns could be compared. Results from point spread function studies, simulations, phantom and in vivo experiments show that the choice of radial sampling pattern influences the quality of the final image reconstructed by the compressed sensing algorithm. While evenly undersampled radial trajectories are best for specific temporal resolutions, golden‐angle radial sampling results in the least overall error when various temporal resolutions are considered. Reduced temporal fluctuations from aliasing artifacts in golden‐angle sampling translates to improved compressed sensing reconstructions overall. Magn Reson Med, 2012. © 2011 Wiley Periodicals, Inc.     

Three-dimensional radial ultrashort echo-time imaging with T2 adapted sampling.

The application of 3D radial sampling of the free-induction decay to proton ultrashort echo-time (UTE) imaging is reported. The effects of T2 decay during signal acquisition on the 3D radial point-spread function are analyzed and compared to 2D radial and 1D sampling. It is found that in addition to the use of ultrashort TE, the proper choice of the acquisition-window duration TAQ is essential for imaging short-T2 components. For 3D radial sampling, a maximal signal-to-noise ratio (SNR) with negligible decay-induced loss in spatial resolution is obtained for an acquisition-window duration of TAQ approximately 0.69 T2. For 2D and 1D sampling, corresponding values are derived as well. Phantom measurements confirm the theoretical findings and demonstrate the impact of different acquisition-window durations on SNR and spatial resolution for a given T2 component. In vivo scans show the potential of 3D UTE imaging with T2-adapted sampling for musculoskeletal imaging using standard MR equipment. The visualization of complex anatomy is demonstrated by extracting curved slices from the isotropically resolved 3D UTE image data.     

A modified projection reconstruction trajectory for reduction of undersampling artifacts

PURPOSE: To reduce undersampling artifacts for a given number of repetitions of the projection reconstruction (PR) sequence by modifying its k-space trajectory to sample more mid-frequencies while reducing the sampling coverage of the peripheral spatial frequencies. MATERIALS AND METHODS: The single k-space spoke measured per repetition in the standard PR was modified so that one complete and two partial spokes were measured per repetition but with decreased k-space extent. The point spread functions (PSFs) and undersampling artifacts of the modified PR were compared with those of the standard PR for various numbers of projections. Phantom and in vivo images were used to assess the relative performance. RESULTS: PSF analysis indicated that the modified PR method provided reduced undersampling artifacts with somewhat reduced spatial resolution. The phantom and in vivo images corroborated this. CONCLUSION: The modified PR trajectory provides reduced undersampling artifact vs. the standard PR, particularly when the number of projections is limited and the artifact level is high.     

Self-calibration method for radial GRAPPA/k-t GRAPPA.

Generalized autocalibrating partially parallel acquisitions (GRAPPA), an important parallel imaging technique, can be easily applied to radial k-space data by segmenting the k-space. The previously reported radial GRAPPA method requires extra calibration data to determine the relative shift operators. In this work it is shown that pseudo-full k-space data can be generated from the partially acquired radial data by filtering in image space followed by inverse gridding. The relative shift operators can then be approximated from the pseudo-full k-space data. The self-calibration method using pseudo-full k-space data can be applied in both k and k-t space. This technique avoids the prescans and hence improves the applicability of radial GRAPPA to image static tissue, and makes k-t GRAPPA applicable to radial trajectory. Experiments show that radial GRAPPA calibrated with pseudo-full calibration data generates results similar to radial GRAPPA calibrated with the true full k-space data for that image. If motion occurs during acquisition, self-calibrated radial GRAPPA protects structural information better than externally calibrated GRAPPA. However, radial GRAPPA calibrated with pseudo-full calibration data suffers from residual streaking artifacts when the reduction factor is high. Radial k-t GRAPPA calibrated with pseudo-full calibration data generates reduced errors compared to the sliding-window method and temporal GRAPPA (TGRAPPA).     

Dynamic radial projection MRI of inhaled hyperpolarized 3He gas.

A radial projection sliding-window sequence has been developed for imaging the rapid flow of (3)He gas in human lungs. The short echo time (TE) of the radial sequence lends itself to fast repetition times, and thus allows a rapid update in the image when it is reconstructed with a sliding window. Oversampling in the radial direction combined with angular undersampling can further reduce the time needed to acquire a complete image data set, without significantly compromising spatial resolution. Controlled flow phantom experiments using hyperpolarized (3)He gas exemplify the temporal resolution of the method. In vivo studies on three healthy volunteers, one patient with chronic obstructive pulmonary disease (COPD), and one patient with hemiparalysis of the right diaphragm demonstrate that it is possible to accurately resolve the passage of gas down the trachea and bronchi and into the peripheral lung.     

Myocardial wall tagging with undersampled projection reconstruction.

Azimuthally undersampled projection reconstruction (PR) acquisition is investigated for use in myocardial wall tagging with MR using grid tags to provide increased temporal and spatial resolution. PR can provide the high-resolution images required for tagging with very few projections, at the expense of artifact. Insight is provided into the PR undersampling artifact, in the context of measuring myocardial motion with tags. For Fourier transform imaging, at least 112 phase-encodings must be collected to image tagging grids spaced 7 pixels apart. PR requires about 80 projections, a 1.4-fold reduction in scan time. Magn Reson Med 45:562-567, 2001. Published 2001 Wiley-Liss, Inc.     

Inversion recovery with embedded self‐calibration (IRES)

With self‐calibrated parallel acquisition, the calibration data used to characterize coil response are acquired within the actual, parallel scan. Although this eliminates the need for a separate calibration scan, it reduces the net acceleration factor of the parallel scan. Furthermore, this reduction gets worse at higher accelerations. A method is described for three‐dimensional inversion recovery gradient‐echo imaging in which calibration is incorporated into the sequence but with no loss of net acceleration. This is done by acquiring the calibration data using very small (≤4°) tip angle acquisitions during the delay interval after acquisition of the accelerated imaging data. The technique is studied at 3 Tesla with simulation, phantom, and in vivo experiments using both image‐space‐based and k‐space‐based parallel reconstruction methods. At nominal acceleration factors of 3 and 4, the newly described inversion recovery with embedded self‐calibration (IRES) method can retain effective acceleration with comparable SNR and contrast to standard self‐calibration. At a net two‐dimensional acceleration factor of 4, IRES can achieve higher SNR than standard self‐calibration having a nominal acceleration factor of 6 but the same acquisition time. Magn Reson Med, 2009. © 2009 Wiley‐Liss, Inc.     

Characterizing radial undersampling artifacts for cardiac applications.

The undersampled radial acquisition has been widely employed for accelerated (by a factor R = N(r)/N(p)) cardiac imaging, but the resulting reduction in image quality has not been well characterized. This investigation presents a method of measuring these artifacts through synthetic undersampling of high SNR images (SNR > or = 30). After validating the method in phantoms, the method was applied to a study of short-axis, long-axis, and coronary MRI imaging in healthy subjects. For 60 projections (60 N(p)), the total artifact is approximately 10% for short and long-axis imaging (R = 2.1) and approximately 15% for coronary MRI (R = 3.7). For 60 N(p), the SD of artifact in the region of the heart is 2% for short- and long-axis imaging (R = 2.1) and 3.5% for coronary MRI (R = 3.7). The artifact content is less in the region of the heart than in the periphery. The artifact is very reproducible among subjects for standard views. A study of coronary MRI at progressively fewer projections (at constant scan time) showed that right coronary MRI images were acceptable if total artifact was <6.5% of image content (N(p) > 120, R = 2.1).     

Sodium MRI using a density‐adapted 3D radial acquisition technique

A density‐adapted three‐dimensional radial projection reconstruction pulse sequence is presented which provides a more efficient k‐space sampling than conventional three‐dimensional projection reconstruction sequences. The gradients of the density‐adapted three‐dimensional radial projection reconstruction pulse sequence are designed such that the averaged sampling density in each spherical shell of k‐space is constant. Due to hardware restrictions, an inner sphere of k‐space is sampled without density adaption. This approach benefits from both the straightforward handling of conventional three‐dimensional projection reconstruction sequence trajectories and an enhanced signal‐to‐noise ratio (SNR) efficiency akin to the commonly used three‐dimensional twisted projection imaging trajectories. Benefits for low SNR applications, when compared to conventional three‐dimensional projection reconstruction sequences, are demonstrated with the example of sodium imaging. In simulations of the point‐spread function, the SNR of small objects is increased by a factor 1.66 for the density‐adapted three‐dimensional radial projection reconstruction pulse sequence sequence. Using analytical and experimental phantoms, it is shown that the density‐adapted three‐dimensional radial projection reconstruction pulse sequence allows higher resolutions and is more robust in the presence of field inhomogeneities. High‐quality in vivo images of the healthy human leg muscle and the healthy human brain are acquired. For equivalent scan times, the SNR is up to a factor of 1.8 higher and anatomic details are better resolved using density‐adapted three‐dimensional radial projection reconstruction pulse sequence. Magn Reson Med, 2009. © 2009 Wiley‐Liss, Inc.     

High-resolution MRI of cardiac function with projection reconstruction and steady-state free precession.

The purpose of this study was to investigate the trabecular structure of the endocardial wall of the living human heart, and the effect of that structure on the measurement of myocardial function using MRI. High-resolution MR images (0.8 x 0.8 x 8 mm voxels) of cardiac function were obtained in five volunteers using a combination of undersampled projection reconstruction (PR) and steady-state free precession (SSFP) contrast in ECG-gated breath-held scans. These images provide movies of cardiac function with new levels of endocardial detail. The trabecular-papillary muscle complex, consisting of a mixture of blood and endocardial structures, is measured to constitute as much as 50% of the myocardial wall in some sectors. Myocardial wall strain measurements derived from tagged MR images show correlation between regions of trabeculae and papillary muscles and regions of high strain, leading to an overestimation of function in the lateral wall.     

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. Magn Reson Med 70:1524–1534, 2013. © 2013 Wiley Periodicals, Inc.     

Lattice permutation for reducing motion artifacts in radial and spiral dynamic imaging.

Radial and spiral trajectories exhibit favorable characteristics for dynamic imaging. Nevertheless, changes in image contents during acquisition lead to inconsistencies in the k-space data, which are manifested as streaks or spiral artifacts, respectively. This work proposes the concept of lattice permutation to reorder the data segments for artifact suppression. This acts to reshuffle the alias pattern along the temporal frequency axis. The proposed approach is well suited to sliding window reconstruction, although more sophisticated methods are also possible. For typical image series where the signal energies are concentrated in the low temporal frequencies, the permutation displaces most of the aliased signals from the low temporal frequencies to the high temporal frequencies, where they are attenuated by sliding window reconstruction, while the signals in the low temporal frequencies are mostly contaminated by aliasing from the much weaker signals in the higher temporal frequencies. This results in considerably reduced artifacts without any increase in scan time. In practice, lattice permutation achieves similar artifact suppression as the bit-reversed order, but with a less stringent restriction on the number of segments. At the same time, it provides a more powerful approach to controlling the alias pattern exactly. Results from real-time cardiac imaging are demonstrated.     

Helical MR: continuously moving table axial imaging with radial acquisitions.

A technique for extended field of view MRI is presented. Similar to helical computed tomography, the method utilizes a continuously moving patient table, a 2D axial slice that remains fixed relative to the MRI magnet, and a radial k-space trajectory. A fully refocused SSFP acquisition enables spatial resolution comparable to current clinical protocols in scan times that are sufficiently short to allow a reasonable breathhold duration. RF transmission and signal reception are performed using the RF body coil and the images are reconstructed in real time. Experimental results are presented that illustrate the technique's ability to resolve small structures in the table-motion direction. Simulation experiments to study the steady-state response of the fully refocused SSFP acquisition during continuous table motion are also presented. Finally, whole body images of healthy volunteers demonstrate the high image quality achieved using the helical MRI approach.