Real-time cardiovascular magnetic resonance at 1.5 T using balanced SSFP and 40 ms resolution
- Dirk Voit†1,
- Shuo Zhang†1, 4,
- Christina Unterberg-Buchwald2, 3, 4,
- Jan M Sohns3, 4,
- Joachim Lotz3, 4 and
- Jens Frahm1, 4Email author
© Voit et al.; licensee BioMed Central Ltd. 2013
Received: 10 May 2013
Accepted: 4 September 2013
Published: 12 September 2013
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.
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.
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.
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.
KeywordsReal-time CMR Cardiovascular magnetic resonance SSFP Arrhythmias Abnormal wall motion
State-of-the-art cardiovascular magnetic resonance (CMR) relies on ECG-synchronized cine acquisitions with balanced SSFP contrast, typically at a magnetic field strength of 1.5 T [1–3]. On the other hand, recent work at 3 T reports significant potential for real-time CMR of cardiac function  and cardiovascular flow [5, 6] at high spatiotemporal resolution usingT1-weighted radial FLASH with pronounced undersampling and image reconstruction by regularized nonlinear inversion (NLINV) [7, 8]. This is because real-time SSFP CMR at 3 T is frequently affected by off-resonance “banding” artifacts due to magnetic field inhomogeneities and further limited by SAR regulations when trying to achieve sufficiently high flip angles for repetition times as short as 3 ms. In order to expand real-time CMR applications to a wider clinical community, this work combines the aforementioned acquisition and reconstruction advances with SSFP CMR at 1.5 T. Real-time CMR not only improves patient compliance because of free breathing and eventually shorter examination times, it also offers extended diagnostic opportunities by providing functional information about individual cardiac cycles and access to immediate physiologic responses to stress and exercise.
The proposed method employs a highly undersampled radial gradient-echo CMR technique with fully balanced gradients  in conjunction with serial image reconstruction by NLINV . Extending previous studies by others [10–12], the method is directly applicable without the need for any calibration scan or averaging of multiple k-space data over time to obtain coil sensitivity maps for image reconstruction. This is because the underlying principle relies on an advanced parallel imaging method  that jointly estimates the image and all coil sensitivities from the undersampled dataset of a single frame in a self-consistent manner and thereby optimally exploits the available limited information.
The purpose of this work was twofold: First, to demonstrate the technical feasibility and achievable image quality of real-time radial SSFP CMR with NLINV reconstruction using studies of healthy subjects on two different 1.5 T CMR systems, and second, to evaluate the future clinical potential using preliminary applications to patients with cardiac arrhythmias and abnormal wall motion.
Six participants with no known illness were recruited among the students of the local university. All subjects gave written informed consent before each CMR examination. A total of 27 consecutive patients (19 male, 8 female, age range 16–82 years, mean 54.0 ± 19.8 years) scheduled for cardiac diagnostics (mainly ischemia, function and scar) underwent additional real-time CMR. Respective movies along all relevant anatomic orientations were recorded as a final adjunct to the standard examination protocol lasting for a total duration of about 3 min. Informed consent was obtained from all patients. The local ethics board Ethikkommission Universitätsmedizin Göttingen approved the study.
Undersampled radial gradient-echo CMR with balanced gradients
All acquisitions were based on a gradient-echo CMR sequence with radial encoding and fully balanced gradients as previously described . The design results in a zero gradient moment (zero phase) for each repetition interval TR and a symmetric echo at TE = TR/2. The actual implementations employed a highly undersampled radial encoding scheme (i.e., 13 to 15 spokes per frame) where 5 successive acquisitions comprised complementary sets of spokes. The spokes of each frame were equally distributed over a full 360° circle in order to homogeneously sample k space. NLINV reconstructions  benefit from complementary sets of spokes by temporal regularization to the preceding frame, which constrains the range of possible solutions to the inverse problem. Moreover, the concept forces residual streakings to “flicker” from frame to frame which facilitates their removal by a post-processing temporal median filter.
Image reconstruction by NLINV was performed offline. Online monitoring at the same frame rate was accomplished by a sliding-window technique that relied on gridding and inverse FFT of a composite dataset of 5 consecutive radial acquisitions (i.e., 65 to 75 spokes) followed by subsequent shifts of the composite dataset by 1/5 . Offline reconstructions of serial images from the highly undersampled single-frame acquisitions employed a parallelized algorithm of the NLINV method . It was implemented on a computer equipped with 8 graphical processing units (GPU) and a specially designed library of newly developed basis functions for the efficient computing on multiple GPUs . The current version achieves a reconstruction speed of about 20 frames per second.
Real-time SSFP CMR
The real-time SSFP CMR sequence with radial undersampling was implemented on two commercial 1.5 T CMR systems (Siemens Healthcare, Erlangen, Germany). They were equipped with different gradient systems and used for studying healthy subjects (1.5 T Avanto) and patients (1.5 T TIM Symphony), respectively. In either case, radiofrequency excitation was accomplished with a body coil, while signal reception employed a multi-channel coil arrangement consisting of an anterior 8-coil array and two or three four-coil arrays from the spine coil (automatically selected depending on image position and orientation).
Acquisition parameters for real-time SSFP CMR
Maximum gradient/mT m-1
Slew rate/mT m-1 ms-1
1.8 × 1.8 × 6.0
1.8 × 1.8 × 8.0
256 × 256
256 × 256
Nyquist undersampling factor
Real-time CMR movies with a typical duration of 15 s (i.e., 12 to 15 heartbeats) covered the heart in sequential scans along anatomically defined orientations including three short-axis views (base, midcardial, apical) as well as a two-chamber, three-chamber, and four-chamber view.
Cine SSFP CMR
Patient studies with conventional ECG-synchronized cine SSFP CMR (Cartesian encoding) and breath-holding conditions involved multiple views at 1.5 to 1.8 mm resolution and 6 or 7 mm section thickness. The flip angle was 60°, the repetition time 40 ms and the bandwidth 930 Hz pixel-1. The FOVs of 380 × 309 mm2 to 492 × 400 mm2 were generally larger than necessary for real-time radial SSFP CMR to avoid aliasing. The magnetic field homogeneity (product shim) was neither optimized for cine nor real-time SSFP CMR.
Image quality was assessed by visual inspection of real-time and cine SSFP CMR movies by two blinded observers with 3 years (JMS) and 8 years (CU) experience of CMR evaluation, respectively. Scans were evaluated according to an image quality score ranging from 0 = no diagnostic quality to 1 = reduced diagnostic quality, 2 = many artifacts, 3 = some artifacts, and 4 = optimal diagnostic quality. Artifacts that were considered to affect the image quality primarily referred to SSFP banding artifacts. Cine CMR showed additional problems due to involuntary motions and real-time CMR presented with residual streakings.
Functional analyses were done by a team of experienced radiologists and cardiologists. A quantitative analysis of cardiac mass, chamber volumes and ejection fraction was achieved with the use of a standard software package (Medis Qmass, Leiden, The Netherlands) applied to each individual cardiac cycle separately. For real-time CMR the analyses were restricted to 5 consecutive heartbeats, while cine breath-hold CMR typically averaged data of 10–12 cardiac cycles. Comparative evaluations for real-time and cine SSFP CMR are reported for 19 out of 27 patients where a sufficient number of similar views were available. Nevertheless, residual differences in section alignment are unavoidable due to different positions during free breathing (real-time CMR) and breath-holding (cine CMR). Patients with normal chamber anatomy were analyzed using long-axis views (four-chamber and vertical long axis), while short-axis views served for the evaluation of patients with abnormal wall motion.
Image quality scores for real-time and cine SSFP CMR of patients
Number of views#
3.33 ± 0.77
3.70 ± 0.58
3.25 ± 0.67
3.66 ± 0.69
3.20 ± 0.72
3.73 ± 0.54
3.21 ± 0.88
3.71 ± 0.67
3.67 ± 0.51
3.85 ± 0.43
3.50 ± 0.64
3.81 ± 0.48
3.10 ± 0.98
3.56 ± 0.64
Figure 3 depicts diastolic and systolic frames selected from three short-axis movies of a patient with abnormal motion of the inferior and infero-lateral segments of the left-ventricular myocardial wall due to known myocardial infarction. The lack or severe reduction of wall thickening and contraction is well depicted in all three sections though less pronounced in the basal plane and more clearly in the two lower sections closer to the apex. These motion abnormalities are even more easily detected in corresponding real-time CMR movies [see Additional files 3 and 4] in a short-axis and three-chamber view, respectively.
Evaluations of real-time and cine SSFP CMR (19 patients)
Real-time CMR free breathing
Cine CMR breath-holding
35.5 ± 10.5
35.1 ± 10.5
Enddiastolic left-ventricular volume/ml
49.5 ± 13.9
52.1 ± 18.2
Endsystolic left-ventricular volume/ml
26.7 ± 13.9
24.4 ± 11.9
22.8 ± 5.2
28.5 ± 15.7
49 ± 12#
55 ± 14
Cardiovascular imaging probably represents the most relevant field of application for real-time MRI. This is because examinations of myocardial anatomy, function, and blood flow during free breathing and independent of an ECG-synchronized data acquisition promise improved patient compliance, extended diagnostic capabilities, and temporal (economic) efficiency. The present work demonstrates that such applications may not only be performed at 3 T using T1-weighted FLASH [4–8], but also at a lower field strength of 1.5 T and with the use of the hitherto more “conventional” SSFP contrast. Real-time CMR acquisitions based on undersampled radial SSFP sequences may be implemented on existing CMR scanners without any hardware modification. On the other hand, online visualization of the NLINV reconstructions requires a parallelized version of the algorithm and a GPU-based computer bypassing the conventional reconstruction pipeline (see below).
Previous radial SSFP approaches to CMR predominantly aimed at a more efficient (i.e., accelerated) 3D coverage of the heart, e.g. see [16–19]. Conversely, radial SSFP CMR has only rarely been applied for cross-sectional real-time CMR [20–22]. In comparison to the present results, however, such studies were restricted to healthy volunteers and offered lower spatial resolution of ≥ 2.5 mm as well as much longer acquisition times of 100–300 ms, while high display rates in terms of frames per second could only be achieved by view sharing and sliding-window reconstruction. More recent attempts that provide similar high frames rates as achieved here are k-t BLAST and SENSE approaches with optimized sampling schemes , through-time radial GRAPPA techniques  or combinations based on k-t SPARSE-SENSE . These techniques enhance the degree of undersampling of conventional parallel imaging by exploiting temporal information, but require either the acquisition of lengthy calibration scans prior to each real-time acquisition or extensive data averaging over time to obtain the coil sensitivity maps needed for image reconstruction. In contrast, this work describes the implementation and first application of a robust real-time SSFP CMR method which is truly self-consistent in the sense that no pre-scan or any other supplementary information is required for serial image reconstruction of individual frames.
In a technical sense, the current real-time images eventually suffered from a local and transient appearance of residual streaking artifacts which are due to imperfect estimations from highly undersampled data. It turns out that such problems mostly originate from the site of a pronounced SSFP artifact which, for example, may be due to air in the abdominal regions below the heart. This is because the SSFP problems usually involve signals with both high contrast and high spatial frequencies, which are particularly sensitive to extreme undersampling. This may also apply to small anatomical structures with bright intensities. Future applications will be able to further reduce residual streakings by taking advantage of receiver coils with a larger number of independent elements.
In general, none of the aforementioned problems resulted in a reduced diagnostic assessment of the left or right ventricle. In absolute terms, the achieved image quality of the real-time CMR movies was slightly lower than that of cine recordings – provided the latter were not affected by respiratory motions or irregular heartbeats. If such problems can be excluded, cine acquisitions that extend over several seconds must of course yield images with better SNR and spatial resolution than obtainable within 30 to 40 ms. Nevertheless, the quality of the real-time short-axis views at 1.8 mm nominal resolution turned out to be almost comparable to that of the cine studies at 1.5 mm resolution (see also Additional files 1, 2, 3 and 4). This is because the ability to readily identify wall motions and blood-myocardium borders may be more relevant for a clinical assessment than the mere spatial resolution. Thus, the contrast and temporal accuracy of the real-time images may lead to a similar diagnostic utility as obtained for cine images at slightly higher spatial resolution. Moreover, 1.5 T CMR systems with multi-element coils and faster gradient systems are expected to further improve the technical quality of real-time images.
Functional evaluations of the preliminary patient data resulted in quantitatively comparable results for real-time and cine CMR with the exception of a 10% lower ejection fraction (averaged across patients). While such differences are difficult to interpret in the absence of an in vivo gold standard, they may reflect systematic differences in image acquisition and physiological condition when trying to compare ECG-synchronized cine CMR with breath-holding, which merges data from 10–12 heartbeats, with real-time CMR of individual heartbeats during free breathing. The present findings suggest that real-time SSFP CMR will generally reach diagnostic quality, while being superior for patients with irregular heartbeats.
At this stage, practical limitations of the proposed method are mainly due to the high computational demand of the NLINV reconstruction. In order to fully exploit the real-time CMR capacity, i.e. to overcome the online sliding-window visualization, the NLINV algorithm requires a multi-GPU computer. Such a system is routinely working in our 3 T CMR system where it has been integrated into the software architecture as an “invisible” bypass computer without the need for any user interference . However, it is not yet commercially available and so far has not been implemented into the 1.5 T CMR system used for patient studies. Another challenge is the need for adaptations of the existing post-processing software for parametric evaluations. In order to retain and extract information about pathophysiologic variations, such programs should be able to efficiently and reliably deal with several hundreds of images from multiple heartbeats rather than analyze 20–30 images from a single cardiac cycle. This particularly refers to a fully automated segmentation of the myocardium without the need for manual corrections.
This work describes a robust acquisition and reconstruction technique for real-time CMR at 1.5 T using conventional SSFP contrast. The NLINV method is directly applicable and fully self-consistent using only data from a single frame without the need for any calibration scan or other computation of averaged coil sensitivity maps. The acquisition technique operates on commercial CMR systems without hardware modification, while rapid image reconstructions require a dedicated computer equipped with GPUs. Preliminary applications to patients with arrhythmias and abnormal wall motion promise considerable clinical potential. Real-time vs cine CMR validation studies of well-defined cohorts of patients are now warranted.
Broad-use linear acquisition speed-up technique
Fast Fourier transform
Fast low-angle shot
Field of view
Graphical processing unit
Generalized autocalibrating partially parallel acquisition
Cardiovascular magnetic resonance
Regularized nonlinear inversion
Specific absorption rate
Steady-state free precession
- TE :
- TR :
We are grateful to Prof. Dr. Frank Wacker and Dr. Marcel Gutberlet from the Department of Radiology at the Medizinische Hochschule Hannover, Hannover, Germany, to let us use their 1.5 T Avanto CMR system. We thank Martin Faßhauer for extensive analyses of serial real-time CMR acquisitions. We thank for financial support (S.Z.) by the DZHK (German Centre for Cardiovascular Research) and BMBF (German Ministry of Education and Research).
- Earls JP, Ho VB, Foo TK, Castillo E, Flamm SD: Cardiac MRI: Recent progress and continued challenges. J Magn Reson Imag. 2002, 16: 111-27. 10.1002/jmri.10154.View ArticleGoogle Scholar
- Schär M, Kozerke S, Fischer SE, Boesiger P: Cardiac SSFP imaging at 3 Tesla. Magn Reson Med. 2004, 51: 799-806. 10.1002/mrm.20024.View ArticlePubMedGoogle Scholar
- Ishida M, Kato S, Sakuma H: Cardiac MRI in ischemic heart disease. Circ J. 2009, 73: 1577-88. 10.1253/circj.CJ-09-0524.View ArticlePubMedGoogle Scholar
- Zhang S, Uecker M, Voit D, Merboldt KD, Frahm J: Real-time cardiovascular magnetic resonance at high temporal resolution: radial FLASH with nonlinear inverse reconstruction. J Cardiovasc Magn Reson. 2010, 12: 39-10.1186/1532-429X-12-39.PubMed CentralView ArticlePubMedGoogle Scholar
- Joseph AA, Merboldt KD, Voit D, Zhang S, Uecker M, Lotz J, Frahm J: Real-time phase-contrast MRI of cardiovascular blood flow using undersampled radial fast low-angle shot and nonlinear inverse reconstruction. NMR Biomed. 2012, 25: 917-24. 10.1002/nbm.1812.View ArticlePubMedGoogle Scholar
- Joseph AA, Kowallick JT, Merboldt KD, Voit D, Schaetz S, Zhang S, Sohns JM, Lotz J, Frahm J: Real-time flow MRI of the aorta at a resolution of 40 ms. J Magn Reson Imag. 2013, in pressGoogle Scholar
- Uecker M, Zhang S, Voit D, Karaus A, Merboldt KD, Frahm J: Real-time MRI at a resolution of 20 ms. NMR Biomed. 2010, 23: 986-94. 10.1002/nbm.1585.View ArticlePubMedGoogle Scholar
- Uecker M, Zhang S, Voit D, Merboldt KD, Frahm J: Real-time MRI – Recent advances using radial FLASH. Imag Med. 2012, 4: 461-76. 10.2217/iim.12.32.View ArticleGoogle Scholar
- Zhang S, Block KT, Frahm J: Magnetic resonance imaging in real time: advances using radial FLASH. J Magn Reson Imag. 2010, 31: 101-9. 10.1002/jmri.21987.View ArticleGoogle Scholar
- Tsao J, Kozerke S, Boesiger P, Pruessmann KP: Optimizing spatiotemporal sampling for k-t BLAST and k-t SENSE: application to high-resolution real-time cardiac steady-state free precession. Magn Reson Med. 2005, 53: 1372-82. 10.1002/mrm.20483.View ArticlePubMedGoogle Scholar
- Seiberlich N, Ehses P, Duerk J, Gilkeson R, Griswold M: Improved radial GRAPPA calibration for real-time free-breathing cardiac imaging. Magn Reson Med. 2011, 65: 492-505. 10.1002/mrm.22618.PubMed CentralView ArticlePubMedGoogle Scholar
- Feng L, Srichai MB, Lim RP, Harrison A, King W, Adluru G, Dibella EVR, Sodickson DK, Otazo R, Kim D: Highly accelerated real-time cardiac cine MRI using k–t SPARSE-SENSE. Magn Reson Med. 2013, 70: 64-74. 10.1002/mrm.24440.PubMed CentralView ArticlePubMedGoogle Scholar
- Uecker M, Hohage T, Block KT, Frahm J: Image reconstruction by regularized nonlinear inversion – Joint estimation of coil sensitivities and image content. Magn Reson Med. 2008, 60: 674-82. 10.1002/mrm.21691.View ArticlePubMedGoogle Scholar
- Uecker M, Zhang S, Frahm J: Nonlinear inverse reconstruction for real-time MRI of the human heart using undersampled radial FLASH. Magn Reson Med. 2010, 63: 1456-62. 10.1002/mrm.22453.View ArticlePubMedGoogle Scholar
- Schätz S, Uecker M: A multi-GPU programming library for real-time applications. In: algorithms and architectures for parallel processing (Springer). Lect Notes Comp Sci. 2012, 7439: 114-28. 10.1007/978-3-642-33078-0_9.View ArticleGoogle Scholar
- Peters DC, Ennis DB, Rohatgi P, Syed MA, McVeigh ER, Arai AE: 3D breath-held cardiac function with projection reconstruction in steady state free precession validated using 2D cine MRI. J Magn Reson Imag. 2004, 20: 411-6. 10.1002/jmri.20145.View ArticleGoogle Scholar
- Arunachalam A, Samsonov A, Block WF: Self-calibrated GRAPPA method for 2D and 3D radial data. Magn Reson Med. 2007, 57: 931-8. 10.1002/mrm.21223.View ArticlePubMedGoogle Scholar
- Liu J, Wieben O, Jung Y, Samsonov AA, Reeder SB, Block WF: Single breathhold cardiac CINE imaging with multi-echo three-dimensional hybrid radial SSFP acquisition. J Magn Reson Imag. 2010, 32: 434-40. 10.1002/jmri.22269.View ArticleGoogle Scholar
- Liu J, Spincemaille P, Codella NCF, Nguyen TD, Prince MR, Wang Y: Respiratory and cardiac self-gated free-breathing cardiac CINE imaging with multiecho 3D hybrid radial SSFP acquisition. Magn Reson Med. 2010, 63: 1230-7. 10.1002/mrm.22306.PubMed CentralView ArticlePubMedGoogle Scholar
- Kühl HP, Spuentrup E, Wall A, Franke A, Schröder J, Heussen N, Hanrath P, Günther RW, Buecker A: Assessment of myocardial function with interactive non–breath-hold real-time MR imaging: comparison with echocardiography and breath-hold cine MR imaging. Radiology. 2004, 231: 198-207. 10.1148/radiol.2311021237.View ArticlePubMedGoogle Scholar
- Winkelmann S, Schaeffter T, Koehler T, Eggers H, Doessel O: An optimal radial profile order based on the golden ratio for time-resolved MRI. IEEE Trans Med Imag. 2007, 26: 68-76.View ArticleGoogle Scholar
- Bauer RW, Radtke I, Block KT, Larson MC, Kerl JM, Hammerstingl R, Graf TG, Vogl TJ, Zhang S: True real-time cardiac MRI in free breathing without ECG synchronization using a novel sequence with radial k-space sampling and balanced SSFP contrast mode. Int J Cardiovasc Imag. 2013, 29: 1059-67. 10.1007/s10554-013-0183-0.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.