PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
J Magn Reson Imaging. Author manuscript; available in PMC 2013 June 1.
Published in final edited form as:
PMCID: PMC3354015
NIHMSID: NIHMS349369

DIfferential Subsampling with Cartesian Ordering (DISCO): a High Spatio-temporal Resolution Dixon Imaging Sequence for Multiphasic Contrast Enhanced Abdominal Imaging

Abstract

Purpose

To develop and evaluate a multiphasic contrast-enhanced MRI method called DIfferential Sub-sampling with Cartesian Ordering (DISCO) for abdominal imaging.

Materials and Methods

A three-dimensional, variable density pseudo-random k-space segmentation scheme was developed and combined with a Dixon-based fat-water separation algorithm to generate high temporal resolution images with robust fat suppression and without compromise in spatial resolution or coverage. With IRB approval and informed consent, 11 consecutive patients referred for abdominal MRI at 3T were imaged with both DISCO and a routine clinical 3D SPGR-Dixon (LAVA FLEX) sequence. All images were graded by two radiologists using quality of fat suppression, severity of artifacts, and overall image quality as scoring criteria. For assessment of arterial phase capture efficiency, the number of temporal phases with angiographic phase and hepatic arterial phase was recorded.

Results

There were no significant differences in quality of fat suppression, artifact severity or overall image quality between DISCO and LAVA FLEX images (p > 0.05, Wilcoxon signed rank test). The angiographic and arterial phases were captured in all 11 patients scanned using the DISCO acquisition (mean number of phases were 2 and 3 respectively).

Conclusion

DISCO effectively captures the fast dynamics of abdominal pathology such as hyperenhancing hepatic lesions with a high spatio-temporal resolution. Typically, 1.1×1.5×3 mm spatial resolution over 60 slices was achieved with a temporal resolution of 4–5 seconds.

Keywords: Dynamic contrast enhanced MRI, Tumor imaging, Time resolved imaging, Dixon fat-water separation

INTRODUCTION

Multi-phasic contrast enhanced MRI (CEMRI) is commonly used for detection and characterization of primary and metastatic lesions in the abdomen and pelvis. However, compromises are often made between temporal and spatial resolution. A high spatial resolution is required to characterize lesion morphology and relationship to adjacent structures whereas a high temporal resolution is required to accurately characterize and/or quantify contrast uptake as well as to capture the arterial phase of enhancement. In hepatic imaging in particular, the late arterial phase is critical, particularly for tumors such as hepatocellular carcinoma (HCC) and hypervascular metastases. These lesions often exhibit brief arterial phase enhancement and potentially become iso-intense with the liver on subsequent portal venous and equilibrium phase images. For contrast-enhanced abdominal imaging applications in a clinical setting, a T1-weighted 3D RF spoiled gradient recalled echo (SPGR) sequence with fat suppression is still the most commonly used imaging method, referred to as LAVA/VIBE/THRIVE by three of the major manufacturers (General Electric/Siemens/Philips). The temporal resolution is typically on the order of 15–20s for adequate spatial resolution and coverage.

Various methods have been proposed to address the spatio-temporal resolution tradeoff in MRI. Some of the earlier keyhole-based methods (1,2) suffered from increased blurring and ghosting artifacts in regions where rapid dynamic changes occur, often the regions of high clinical interest (2). Generalized series modeling methods like RIGR or TRIGR (34) can help overcome this problem by using more sophisticated algorithms to extrapolate the missing data but require solution of large systems of equations, introducing instability and noise amplification. A number of k-t segmentation based time resolved imaging schemes have been proposed and demonstrated but have mostly been confined to high spatio-temporal resolution MR angiography (MRA). Time Resolved Imaging of Contrast KineticS (TRICKS) (5) and Cartesian Acquisition with Projection Reconstruction like sampling (CAPR) (6) segment elliptically ordered 3D k-space (7) into annular or radial regions and time interleave the acquisition of central and peripheral k-space regions to achieve temporal acceleration. Highly constrained back projection (HYPR) (8) imaging employs radial undersampling and a composite image to constrain the reconstruction of the subsampled data and eliminate streaking artifacts. Undersampled projection reconstruction imaging with k-space weighted image combination (KWIC) has been demonstrated in breast imaging by Dougherty et al. (9). Parallel imaging techniques (1012) have also shown considerable potential in improving spatio-temporal resolution, often in conjunction with these other methods.

In addition to high spatio-temporal resolution, uniform fat suppression can be critical for lesion conspicuity as well as for avoiding inadvertent water excitation. This is even more relevant at higher field strengths where fat suppression can be challenging due to more heterogeneous B1. Conventional fat suppression techniques apply inversion or saturation pulses centered on the lipid resonance followed by acquisition of multiple lines of k-space in order to restrict scan time increase. Recently, variants of the in-phase/out-of-phase method originally proposed by Dixon (13)- the 3-point Dixon based IDEAL (Iterative Decomposition of water and fat with Echo Asymmetry and Least-squares estimation) (14) and the 2-point based MEDAL (Multi-Echo with 2-Point Dixon Reconstruction for Decomposition of Aqua/Lipid) (15)- have been shown to be robust for fat-water separation, in the presence of B0 inhomogeneity. The elimination of magnetization preparation pulses provides these Dixon-based methods complete flexibility in the order of k-space acquisition, and in particular enables combination with elliptical centric or other novel non-sequential k-space ordering approaches. However, the temporal resolution of these sequences for abdominal imaging is on the order of 15–20s (1617), making them suboptimal for visualization of fast enhancing lesions.

Recently, TRICKS (5) was applied to liver DCEMRI in conjunction with a 2-point Dixon based fat suppression method (18). In this acquisition approach called META, k-space was divided into elliptical regions and the central k-space regions were acquired more frequently than the peripheral regions. At each reconstructed time frame, missing k-space data were estimated using linear interpolation from previous and succeeding time frames. One drawback of this method was the interpolation across successive breath-holds, which could lead to ghosting artifacts due to inconsistent breath-holding. This could partly be addressed using higher channel-count receive array coils, which would enable higher parallel imaging acceleration and better view sharing schemes.

In this study, we develop and demonstrate clinical feasibility of a new high spatio-temporal resolution DCEMRI technique, referred to as DISCO (DIfferential Sub-sampling with Cartesian Ordering), that combines a dual-echo SPGR sequence with pseudo-random variable density k-space segmentation and a view sharing reconstruction. DISCO samples an elliptically ordered central k-space region every time, and sub-samples the outer regions with pseudo-random segmentation such that the aliasing artifacts from subsampling are rendered somewhat incoherent. Robust fat suppression is achieved using a two-point Dixon fat-water reconstruction algorithm, which is completely compatible with the irregular k-space ordering scheme as described above.

MATERIALS AND METHODS

Pulse sequence and reconstruction

A variable density k-space undersampling strategy was developed to generate a pseudo-random distribution of k-space points. All points in ky-kz space are first sorted in the order of increasing kr = √(ky2 + kz2). This elliptically ordered ky-kz space is segmented into N annular regions and each region i is subsampled by a factor of Si (typically Si >= i) with the central region being fully sampled and the outer regions being progressively subsampled by retaining every Sith sample (again in order of increasing kr within the region). For each region, Si different sampling patterns can be generated by staggering the starting sample from 1 through Si. Note that all k-space points are still confined to a Cartesian grid to retain the advantages of Cartesian image reconstruction. We call this method DIfferential Sub-sampling with Cartesian Ordering or DISCO.

The temporal footprint (i.e. total time extent of samples used to reconstruct a temporal phase) was minimized by nearest-neighbor view sharing of the undersampled views. View sharing was restricted to nearest neighbors and to k-space segments acquired within the same breath-hold, minimizing temporal blurring and motion mis-registration. For N=2 and Si=i, the region schedule is AB1AB2AB3AB1… where Xi refers to the under-sampled subset of the fully sampled region X (Fig 1b). Since each sub-region is elliptically ordered in increasing kr, points are traversed smoothly and eddy current artifacts that can accrue due to true stochastic sampling are mitigated while still imparting a degree of “randomness”, when viewed in ky-kz space.

Figure 1
The DISCO k-space segmentation scheme (a), acquisition schedule (b) and the pulse sequence diagram of the proposed DISCO scheme (c). The central red region is acquired every time and the outer pseudorandom subsampled regions (blue, green, yellow) are ...

The DISCO k-space segmentation scheme was incorporated into a dual-echo bipolar-readout 3D SPGR sequence depicted in Figure 1c. A 2-point Dixon reconstruction with a region-growing based fat-water separation method proposed in (15) was used for fat-water separation. The use of high receiver bandwidths (±167 kHz) enabled the acquisition of opposed-phase and in-phase echoes at approximately 1.2 ms and 2.4 ms (optimal for 3T) in a single TR, greatly reducing scan times and enabling reliable fat-water separation. Self-calibrated hybrid space parallel imaging (ARC) enabled the use of two dimensional acceleration (along ky and kz) (19). Following nearest neighbor view-sharing and ARC parallel imaging reconstruction, full k-space in-phase and opposed-phase complex images were available at each reconstructed time point. These were then passed through the Dixon reconstruction algorithm with a phase correction as described in (15) to yield fat-only and water-only image volumes at each reconstructed phase. All reconstruction algorithms were implemented on-line on the vendor-supplied (GE Healthcare) 8-processor (2.4 GHz) cluster to minimize image reconstruction latency. Each temporal phase took approximately 60s to reconstruct after the acquisition of the requisite k-space data.

Experiments

Simulations and phantom experiments

MATLAB (MathWorks, Natick MA) simulations were used to compute the Point Spread Function (PSF) of sub-sampled TRICKS, CAPR and DISCO trajectories. For each trajectory, a central k-space region (A) and one peripheral region- a ring in the case of TRICKS, a set of spokes in the case of CAPR, a set of pseudo-randomly distributed points confined to an annulus in the case of DISCO- was used to generate the PSF. The size of the central k-space region and the total number of points was kept constant for all three trajectories. Using a 120×120 matrix, the size of the central A region was 1600 and the peripheral B region 2800 for all three cases.

In order to evaluate the impact of PSF on image quality, a doped resolution phantom was used to compare the effect of motion on the three trajectories. Halfway through a 25s scan, the phantom was moved by 5 mm to simulate motion that could, for example, arise from loss of breath-holding. The scan parameters used were identical to the human subject scans detailed below with the exception that parallel imaging was not used. For both the phantom experiments and PSF simulations above, a single A region (fraction 0.16) and 3 B regions with equal number of k-space points (fraction 0.28) were used to span full k-space, with the B regions being elliptical annuli for TRICKS, radial sectors for CAPR and pseudorandom regions for DISCO.

Human subjects

Human experiments were conducted following IRB approval and informed consent. In order to evaluate DISCO’s motion robustness and reduced temporal footprint, healthy volunteers were asked to hold their breath at slightly different positions in the first and second end-inspiration breath-holds and two sets of images reconstructed- (i) view sharing within a breath-hold and (ii) across breath-holds. Imaging was performed on a GE 3T MR750 system (GE Healthcare, Waukesha, WI) with the following parameters: 12° flip, ±167 kHz bandwidth, TR/TE1/TE2 4.1/1.2/2.4 ms, 320×224 matrix, 30–35 cm FOV, 3 mm slice thickness, 60 slices, ARC 2×2 outer acceleration. Using a net acceleration factor of 3.4 (to account for ARC calibration), the total ky-kz points is ~ 3000, assuming an elliptical k-space region. This yields an A region size of 0.16*3000=480 and 3 B regions of size of 0.28*3000=840. The total points for each phase is hence 1320 points yielding ~5.2s temporal resolution.

After validation experiments in volunteers, eleven consecutive patients referred for abdominal MRI were recruited into the study. A 32-channel torso array coil optimized for high acceleration factors was used with the upper 19 elements enabled for the data acquisition and reconstruction of the abdominal scans. Following localizer scans and axial fat-suppressed fast spin echo T2 weighted imaging, DISCO was run to acquire a full k-space breath-held data set, which also served as a pre-contrast reference phase. For contrast enhanced imaging, gadoxetate was used in 3 patients and gadobenate in 8 patients. After injection of a single dose of gadolinium contrast (0.2 mL/kg for gadobenate, 0.1 mL/kg for gadoxetate) at 2 cc/sec followed by 20 mL saline chaser at the same rate, the remaining k-space data were acquired as shown in Figure 1a. This second 24–28s breath-hold typically yielded 6–7 temporal phases with a temporal resolution of 4–5s. The timing for this dynamic breath-hold was fixed at 15s following contrast injection to minimize any possibility of early arrival of contrast. For image quality comparisons, a 3D SPGR with 2-point Dixon reconstruction (3D LAVA FLEX) scan was acquired in an immediately subsequent breath-hold with the same sequence parameters including parallel imaging acceleration factors. The temporal resolution of this sequence was ~14s. All breath-hold image acquisitions including the scout scans were performed during end-inspiration to maximize patient comfort.

Data analysis

In all patients, the 3D LAVA FLEX scan acquired immediately after DISCO was used for comparing overall image quality, severity of artifacts, and quality of fat suppression. The last DISCO phase was used for the comparisons to use a temporal phase closest to that of the LAVA FLEX scan. The number of angiographic and arterial temporal phases that were captured by the DISCO acquisition were also recorded. The criteria for designating a temporal phase as an angiographic phase were both a bright hepatic artery and no portal vein enhancement. The criteria for a hepatic arterial phase were some portal venous enhancement that was less intense than that of the hepatic artery, as well as no enhancement of hepatic veins. The images were randomized and graded for quality of liver delineation on a scale of 0–4 by consensus review of two experienced board-certified radiologists. The categories for image scoring were overall image quality, artifact severity and quality of fat suppression. Table 1 summarizes the criteria and the scoring scale. A non-parametric Wilcoxon signed rank test was used to compare qualitative ordinal image scores.

Table 1
Scoring criteria for image comparisons between DISCO and LAVA FLEX

RESULTS

Figure 2 compares subsampling point spread functions obtained from simulations of TRICKS, CAPR and DISCO k-space segmentation schemes (a–c) for the same total number of sampled points. Note the significant dispersal of ghosting energy in DISCO (f) compared to that of TRICKS (d) and CAPR (e) due to the pseudo-random nature of sampling in DISCO suggesting that artifacts arising due to view sharing and potential signal discontinuities between regions being shared may be much more dispersed in DISCO compared to TRICKS or CAPR. Comparable sections from a 3D acquisition on a resolution phantom, where the phantom was displaced by 5 mm midway through the scan are shown in (g–i) for the three different k-space segmentation schemes. Note that the artifacts (arrows) are stronger and more coherent in TRICKS and CAPR (g–h) compared to DISCO (i), agreeing with the point spread function simulations.

Figure 2
The ky-kz segmentation scheme for TRICKS (a), CAPR (b) and DISCO (c) and the corresponding subsampled point spread functions (d–f). Note that the point spread function for DISCO (f) shows significantly reduced coherent artifacts around the main ...

Comparable slices from the first temporal phase of the second breath-hold of a DISCO scan on a healthy human subject with view sharing spanning across two end inspiration breath-holds (a) and view sharing restricted to within the second end inspiration breath-hold (b) are shown in Figure 3. Note that the ghosting artifacts in (a) (white arrow) due to mis-registration of the first and second breath-holds, are completely eliminated in (b), allowing the use of end-inspiration breath-holding in DISCO.

Figure 3
Comparable slices from the first temporal phase of the second breath-hold of a DISCO scan on a healthy human subject with view sharing spanning across two end inspiration breath-holds (a) and view sharing restricted to within each breath-hold (b). Note ...

Figure 4 shows box-whisker plots of radiologist image score distribution using overall image quality, artifact severity and degree of fat suppression as scoring criteria of the proposed DISCO method and conventional 3D LAVA FLEX sequence. There were no significant differences between the two techniques (p > 0.05) for any of the scoring criteria. At least one angiographic phase and one arterial phase were captured in each of the 11 cases (mean number of angiographic and arterial phases identified- 1.91 ± 0.8 and 3.27 ± 1.5 respectively) attesting to the high temporal resolution of the proposed DISCO technique.

Figure 4
Distribution of scores from radiologist rankings in overall image quality, artifact severity and quality of fat suppression between DISCO and conventional LAVA FLEX acquisitions. There were no significant differences found amongst the two methods for ...

Representative results with eight temporal phases (~4s temporal resolution) obtained on a 67 year old male with two foci of likely hepatocellular carcinoma are shown in Figure 5. Non-contrast-enhanced fully sampled data was obtained in one breath-hold (a), followed by a second breath-hold of seven DISCO dynamic contrast enhancement phases (b–h) showing first arteries (black arrow), tumors (non-dashed and dashed white arrow), portal veins (solid blue arrow), and hepatic veins (thin long black arrow) enhancing with a single-dose gadoxetate injection. Note that the tumors show hyper-enhancement on only two temporal phases, mostly at the time the portal vein has its first blush of contrast (thin dashed black arrow in e). Also note that by this time, contrast has already washed out from the hepatic arteries.

Figure 5
67 year old male with two foci of likely hepatocellular carcinoma. Non-contrast-enhanced fully sampled data obtained in one breath-hold (a), followed by a second breath-hold of seven DISCO contrast enhancement phases (b–h) with ~4s temporal resolution ...

Another multiphasic DISCO acquisition from a 53 year old female with a likely hepatocellular carcinoma is shown in Figure 6. Non-contrast-enhanced fully sampled data was acquired in one breath-hold (a) followed by a second breath-hold of six DISCO dynamic contrast enhancement phases with ~4 s temporal resolution (b–g) showing an enhancing tumor (white arrow) and an adjacent transient hepatic intensity difference (dashed white arrow) with a single-dose gadobenate injection. Note regions of hyper-enhancement of the tumor are rapidly reaching same signal intensity as background liver (black arrows in g). The fully sampled LAVA FLEX data obtained in a third breath-hold a few seconds later (h) shows comparable image quality.

Figure 6
53 year old female with a likely hepatocellular carcinoma. Non-contrast-enhanced data acquired in one breath-hold (a) followed by a second breath-hold of six DISCO contrast enhancement phases (b–g) with ~4s temporal resolution showing enhancing ...

DISCUSSION

Contrast enhanced body MRI applications, especially hepatic and renal imaging, require high spatio-temporal resolution to reliably capture the angiographic and arterial phases with robust fat suppression and minimal motion blurring/artifacts. The proposed DISCO pulse sequence combining a variable density pseudo-random k-space segmentation scheme, view sharing, a two-point Dixon-based fat-water separation, and 2D parallel imaging using ARC provided excellent fat suppression at 3T where conventional fat saturation techniques are often suboptimal. It also provided a high spatio-temporal resolution for excellent visualization of structures with rapid contrast enhancement and washout and visualization of the angiographic and arterial phases in all 11 patient cases scanned. Approximately 6–7 phases were reconstructed from data acquired in a 28–30s end-inspiration breath-hold, resulting in a temporal resolution of ~4s with ~1×1.5×3 mm spatial resolution. The use of a Dixon-based fat-water separation scheme was critical to enable the use of a pseudo random k-space segmentation scheme since conventional fat suppression schemes require sequential or linearly segmented k-space for optimal performance. The image quality of the DISCO images were comparable to and statistically no worse than the 3D LAVA FLEX despite higher frame rates that captured arterial and angiographic phases in all 11 cases and enabled clear visualization of hypervascular tumors.

One advantage of our k-space segmentation and view sharing scheme over the previously proposed META technique (18) is its immunity to breath-hold mismatching. In (18), k-space data was linearly interpolated across multiple breath-holds necessitating the use of end-expiratory breath-holding, which is challenging for many patients. While end-expiratory breath-holding is more consistent, it is still possible to have mismatches between successive breath-holds, resulting in motion artifacts. DISCO eliminates this by restricting the view sharing to within a breath-hold and thereby reducing its temporal footprint compared to META. This can be seen in Figure 3 where the across breath-hold reconstruction resulted in artifacts compared to the within breath-hold reconstruction. Furthermore, the pseudorandom nature of outer k-space in the DISCO k-space segmentation scheme disperses motion artifacts in much the same way as undersampling artifacts are dispersed in radial imaging as illustrated in the PSF simulations of Figure 2. The use of higher channel count surface coils could enable META acquisition with higher acceleration than was reported in (18) but the limitations of a linear interpolation based reconstruction method would still require that some of the k-space data be shared between breath holds or reduce the total number of dynamic phases acquired compared to nearest-neighbor viewsharing.

The DISCO k-space segmentation bears similarity to the recently proposed TWIST scheme, (20) which has been applied to peripheral MRA (21), renal imaging (22) and more recently to breast DCEMRI (23). DISCO uses a variable density undersampling scheme i.e. k-space is progressively undersampled away from the center and is capable of achieving higher subsampling rates than the spiral in- spiral out scheme of TWIST. Unlike TWIST, where the k-space radius kr of the acquired k-space points oscillates as a function of time within each temporal phase, it monotonically increases from the center outwards in DISCO, perhaps making it better suited for capture of contrast kinetics, akin to true elliptical centric k-space ordering. Furthermore, the work of Herrmann et al. (23) used a mask subtraction scheme in the breast to eliminate fat, which is suboptimal for breath-holding applications. Even in breast imaging with relatively minimal motion, the authors had to resort to a retrospective image registration scheme. The use of Dixon fat-water separation in conjunction with DISCO makes it ideal for applications like breast and liver imaging, where fat suppression is especially challenging due to B0/B1 heterogeneity. Recently, Dixon fat-water separation has been shown to be superior to conventional fat suppression in reducing motion artifacts in breast imaging (24).

This preliminary study had some limitations. We studied a small number of patients (n=11) to assess the feasibility and image quality of the DISCO technique. While the last DISCO phase which was temporally close to the LAVA FLEX scan was used for image comparisons, in some cases, they were in different stages of enhancement (late arterial vs. portal venous for example). However, the comparisons were more to evaluate potential artifacts that could be introduced by the DISCO viewsharing scheme rather than features specific to stage of enhancement. A comparison of the arterial phase capture efficiency of the conventional method (3D LAVA FLEX for example) using, say, a cohort of n successive patients imaged was not possible due to large variations in the protocol parameters of the conventional sequence stemming from the tailoring of imaging parameters by the technologist to individual patient size and predicted breath-holding capacity. Using DISCO, the angiographic and arterial phases were captured on each of the 11 cases, spanning an average of 2 and 3 temporal phases respectively. While there were no artifacts introduced due to view sharing compared to the conventional 3D LAVA FLEX scans, some degree of temporal blurring could be expected with view sharing especially in fast enhancing lesions. The inherent pseudo-random sampling of the DISCO trajectory (which results in a relatively incoherent PSF) is compatible with a compressed sensing reconstruction. Future developments would investigate the use of compressed sensing reconstruction techniques to reconstruct a single phase without view sharing. A larger clinical study comparing the DISCO technique to conventional LAVA FLEX as well with alternative reconstruction strategies is currently underway.

Optimizing arterial phase timing for contrast-enhanced MRI can be challenging in hepatic MRI, particularly for lesions such as hepatocellular carcinoma and neuroendocrine metastases which are sometimes visualized briefly during the arterial phase (25). This is particularly relevant when using a contrast agent with a fast vascular clearance such as gadoxetate as illustrated in Figure 5, where the rapid contrast wash-in and wash-out has been effectively captured by the DISCO technique. In addition, the angiographic phase can be of considerable value in planning chemo-embolization procedures for treatment planning in order to visualize collateral arteries to lesions. Our feasibility study demonstrates some of the advantages of the DISCO technique compared to conventional 3D LAVA FLEX or SPGR imaging including significantly improved efficiency of arterial phase capture, and comparable overall image quality and artifacts. Beyond hepatic imaging, other potential applications of this technique include MR urography, and dynamic contrast enhanced imaging of the breast and prostate.

In conclusion, the DISCO technique allows acquisition of 3D volumes with an effective temporal resolution of 4–5 seconds while preserving spatial resolution, fat suppression and image quality and its feasibility has been demonstrated in a clinical setting. This technique may help to alleviate some of the variability in the appearance of arterial phase images from examination to examination and thereby allow more confident assessment of interval growth or stability of small enhancing lesions.

Acknowledgments

Grant support: NIH R01-EB009055 and P41-RR009784

The authors would like to thank Dr. Pauline W. Worters for helpful discussions and assistance with the figures.

References

1. van Vaals JJ, Brummer ME, Dixon WT, et al. "Keyhole" method for accelerating imaging of contrast agent uptake. J Magn Reson Imaging. 1993;3:671–675. [PubMed]
2. Plewes DB, Bishop J, Soutar I, Cohen E. Errors in quantitative dynamic three-dimensional keyhole MR imaging of the breast. J Magn Reson Imaging. 1995;5:361–364. [PubMed]
3. Liang ZP, Lauterbur PC. An efficient method for dynamic magnetic resonance imaging. IEEE Trans Med Imaging. 1994;13:677–686. [PubMed]
4. Hanson JM, Liang ZP, Wiener EC, Lauterbur PC. Fast dynamic imaging using two reference images. Magn Reson Med. 1996;36:172–175. [PubMed]
5. Korosec FR, Frayne R, Grist TM, Mistretta CA. Time-resolved contrast-enhanced 3D MR angiography. Magn Reson Med. 1996;36:345–351. [PubMed]
6. Madhuranthakam AJ, Hu HH, Barger AV, Haider CR, Kruger DG, Glockner JF, Huston J, 3rd, Riederer SJ. Undersampled elliptical centric view-order for improved spatial resolution in contrast-enhanced MR angiography. Magn Reson Med. 2006;55:50–58. [PubMed]
7. Wilman AH, Riederer SJ. Performance of an elliptical centric view order for signal enhancement and motion artifact suppression in breath-hold three-dimensional gradient echo imaging. Magn Reson Med. 1997;38:793–802. [PubMed]
8. Mistretta CA, Wieben O, Velikina J, Block W, Perry J, Wu Y, Johnson K, Wu Y. Highly constrained backprojection for time-resolved MRI. Magn Reson Med. 2006;55:30–40. [PMC free article] [PubMed]
9. Dougherty L, Isaac G, Rosen MA, Nunes LW, Moate PJ, Boston RC, Schnall MD, Song HK. High frame-rate simultaneous bilateral breast DCE-MRI. Magn Reson Med. 2007;57:220–225. [PubMed]
10. Pruessmann KP, Weiger M, Scheidegger MB, Boesiger P. SENSE: sensitivity encoding for fast MRI. Magn Reson Med. 1999;42:952–962. [PubMed]
11. Sodickson DK, Manning WJ. Simultaneous acquisition of spatial harmonics (SMASH): fast imaging with radiofrequency coil arrays. Magn Reson Med. 1997;38:591–603. [PubMed]
12. Vogt FM, Antoch G, Hunold P, et al. Parallel acquisition techniques for accelerated volumetric interpolated breath-hold examination magnetic resonance imaging of the upper abdomen: assessment of image quality and lesion conspicuity. J Magn Reson Imaging. 2005;21:376–382. [PubMed]
13. Dixon WT. Simple proton spectroscopic imaging. Radiology. 1984;153:189–194. [PubMed]
14. Reeder SB, Markl M, Yu H, Hellinger JC, Herfkens RJ, Pelc NJ. Cardiac CINE imaging with IDEAL water-fat separation and steady-state free precession. J Magn Reson Imaging. 2005;22:44–52. [PubMed]
15. Ma J. Breath-hold water and fat imaging using a dual-echo two-point Dixon technique with an efficient and robust phase-correction algorithm. Magn Reson Med. 2004;52:415–419. [PubMed]
16. Rosenkrantz AB, Mannelli Lorenzo, Kim S, Babb JS. Gadolinium-Enhanced Liver Magnetic Resonance Imaging Using a 2-Point Dixon Fat-Water Separation Technique: Impact Upon Image Quality and Lesion Detection. J Comp. Assist. Tomo. 2011;35:96–101. [PubMed]
17. Low RN, Bayram E, Panchal NJ, Estkowski L. High-resolution double arterial phase hepatic MRI using adaptive 2D centric view ordering: initial clinical experience. Am J Roentgenol. 2010;194:947–956. [PubMed]
18. Saranathan M, Rettmann D, Bayram E, Lee C, Glockner J. Multiecho time-resolved acquisition (META): a high spatiotemporal resolution Dixon imaging sequence for dynamic contrast-enhanced MRI. J Magn Reson Imaging. 2009;29:1406–1413. [PubMed]
19. Beatty PJ, Brau AC, Chang S, et al. A method for autocalibrating 2-D accelerated volumetric parallel imaging with clinically practical reconstruction times. Proceedings of the 15th Annual Meeting of ISMRM; Berlin, Germany. 2007. p. 1749. In:
20. Vogt FM, Eggebrecht H, Laub G, Kroeker R, Schmidt M, Barkhausen J, Ladd SC. High spatial and temporal resolution MRA (TWIST) in acute aortic dissection. Proceedings of the 15th Annual Meeting of ISMRM; Berlin, Germany. 2007. p. 92. In:
21. Lim RP, Shapiro M, Wang EY, Law M, Babb JS, Rueff LE, Jacob JS, Kim S, Carson RH, Mulholland TP, Laub G, Hecht EM. 3D time-resolved MR angiography (MRA) of the carotid arteries with time-resolved imaging with stochastic trajectories: comparison with 3D contrast-enhanced Bolus-Chase MRA and 3D time-of-flight MRA. Am J Neuroradiol. 2008;29:1847–1854. [PubMed]
22. Song T, Laine AF, Chen Q, et al. Optimal k-space sampling for dynamic contrast-enhanced MRI with an application to MR renography. Magn Reson Med. 2009;61:1242–1248. [PMC free article] [PubMed]
23. Herrmann KH, Baltzer PA, Dietzel M, Krumbein I, Geppert C, Kaiser WA, Reichenbach JR. Resolving arterial phase and temporal enhancement characteristics in DCE MRM at high spatial resolution with TWIST acquisition. J Magn Reson Imaging. 2011;34:973–982. [PubMed]
24. Shi K, Zhou K, Niu X, Luo B, Han H, Bao S, Ma J. Investigation of motion artifacts associated with fat saturation technique in 3D flash imaging. Med Phys. 2011;38:4556–4562. [PubMed]
25. Ito K, Fujita T, Shimizu A, et al. Multiarterial phase dynamic MRI of small early enhancing hepatic lesions in cirrhosis or chronic hepatitis: differentiating between hypervascular hepatocellular carcinomas and pseudolesions. Am J Roentgenol. 2004;183:699–705. [PubMed]