Scholarly article on topic 'Accelerated CMR using zonal, parallel and prior knowledge driven imaging methods'

Accelerated CMR using zonal, parallel and prior knowledge driven imaging methods Academic research paper on "Medical engineering"

Share paper
OECD Field of science

Academic research paper on topic "Accelerated CMR using zonal, parallel and prior knowledge driven imaging methods"

Journal of Cardiovascular Magnetic

ReSOnanCe BioMed Central

Open Access


Accelerated CMR using zonal, parallel and prior knowledge driven imaging methods

Sebastian Kozerke*1 and Sven Plein2

Address: institute for Biomedical Engineering, University of Zurich and Swiss Federal Institute of Technology, Zurich, Switzerland and 2Academic Unit of Cardiovascular Medicine, University of Leeds, Leeds, UK

Email: Sebastian Kozerke* -; Sven Plein - * Corresponding author

Published: 5 June 2008 Received: 19 May 2008

, , , „. r, , ,„,,„,,,,-,„ „„„w Accepted: 5 June 2008

Journal of Cardiovascular Magnetic Resonance 2008, 10:29 doi:l0.ll86/l 532-429X-10- v '

This article is available from: http://www.jcmr-online.eom/content/l0/l/29 © 2008 Kozerke and Plein; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Accelerated imaging is highly relevant for many CMR applications as competing constraints with respect to spatiotemporal resolution and tolerable scan times are frequently posed. Three approaches, all involving data undersampling to increase scan efficiencies, are discussed in this review. Zonal imaging can be considered a niche but nevertheless has found application in coronary imaging and CMR flow measurements. Current work on parallel-transmit systems is expected to revive the interest in zonal imaging techniques. The second and main approach to speeding up CMR sequences has been parallel imaging. A wide range of CMR applications has benefited from parallel imaging with reduction factors of two to three routinely applied for functional assessment, perfusion, viability and coronary imaging. Large coil arrays, as are becoming increasingly available, are expected to support reduction factors greater than three to four in particular in combination with 3D imaging protocols. Despite these prospects, theoretical work has indicated fundamental limits of coil encoding at clinically available magnetic field strengths. In that respect, alternative approaches exploiting prior knowledge about the object being imaged as such or jointly with parallel imaging have attracted considerable attention. Five to eight-fold scan accelerations in cine and dynamic CMR applications have been reported and image quality has been found to be favorable relative to using parallel imaging alone.

With all acceleration techniques, careful consideration of the limits and the trade-off between acceleration and occurrence of artifacts that may arise if these limits are breached is required. In parallel imaging the spatially varying noise has to be considered when measuring contrast- and signal-to-noise ratios. Also, temporal fidelity in images reconstructed with prior knowledge driven methods has to be studied carefully.


The use of Magnetic Resonance (MR) methods in the cardiovascular realm continues to gain impact. A number of clinical indications have already been identified with Class I rating, including the assessment of global ventricu-

lar function and mass and the detection of acute and chronic myocardial infarction and myocardial scar [1]. A range of further applications have been indicated [2]. For example, recent results from a multi-vendor, multi-center clinical trial showed high diagnostic performance of myo-

cardial perfusion MR imaging [3] compared with nuclear scintigraphy, suggesting that this method will have an important future role in clinical practice.

The technological development in cardiovascular MR (CMR) has been enormous with considerable innovations refining existing or even enabling new applications. A key driving factor for this innovation process has been the challenge to image in the presence of considerable object motion as it is found in the heart. Satisfying competing demands on total scan duration, image quality and spatiotemporal resolution has been a core motivation. With many CMR applications being performed during a breath-hold, improving spatiotemporal resolution for a given, tolerable scan duration is of great importance. Improved spatiotemporal resolution allows discerning finer details of object structure and dynamics. Also, on a fundamental level, spatial and temporal resolutions determine image quality as the degree of partial volume artifacts, image fidelity of dynamic features and motion induced image distortions are affected.

In general, the relative discomfort due to lengthy breath-holds, or in case of free-breathing acquisitions the overall scan duration, are considered drawbacks of CMR over other imaging modalities. In addition, the likelihood of image degradation due to irregular cardiac and respiratory motion tends to amplify with longer measurement times. Therefore, there has been high demand to improve the efficiency of data collection, thereby permitting improved spatiotemporal resolutions or reduced scan times or combinations thereof.

The relatively low imaging speed of MR is inherently linked to the sequential spatial encoding procedure by using time-varying magnetic field gradients. This principle leaves little freedom to reduce scan time for a given resolution and field-of-view in a conventional fashion due to constraints on the tolerable eddy currents induced in the body. A remedy to this problem is to reduce the number of encoding steps. This, inevitably, results in either reduced resolution if the density of the sampled points remains unchanged or in a reduced field-of-view if the density of sampled points is decreased. Since spatial resolution is not intended to be compromised, reduced scan time is typically traded for a reduced field-of-view. Because of the inverse relation between sampling density and field-of-view, one may refer to undersampling or reduced field-of-view methods collectively.

In some applications in which the object-of-interest is considerably smaller than the total field-of-view, scan time reduction by reducing the field-of-view and thus allowing foldover artifacts purposely is a straightforward approach and acceptable as long as the aliased portions

do not fold on the object-of-interest itself. However, this strategy requires some experience in controlling foldover artifacts in particular with double-oblique scanning and only allows for very moderate reductions in scan time.

The concept of using a reduced field-of-view to shorten scan time can be refined by modifying the imaging pulse sequence such that signals from outside a desired reduced region-of-interest are suppressed or not excited [4]. This principle, which will be referred to as zonal imaging hereafter has been successfully used in a range of CMR applications [5-10].

On general grounds, however, suppression of magnetization cannot be ensured or achieved at all with any pulse sequence. Also, the amount of field-of-view reduction and thus scan time reduction is typically very limited given the extent of the object being investigated. To this end, treatment of foldover artifacts from a reduced field-of-view acquisition is required. In mathematical terms, an under-determined problem is posed in such that a full field-of-view image is desired to be reconstructed from data acquired at a reduced field-of-view. In order to make this mathematical problem tractable, complementary spatial encoding functions are required. These are provided if the object signals are received with multiple, independent coils as proposed in first concepts in the late 1980s [11,12]. First practical implementations of parallel imaging using a k-space formalism [13] as well as the mathematical foundations from an image-domain perspective [14] followed ten years later. Today SMASH [13], its successor GRAPPA [15] and SENSE [14] are considered true milestones in MR technology with SENSE and GRAPPA being widely available on commercial MR systems. Instrumental to parallel imaging is the availability of dedicated receive coil arrays [16,17]. State-of-the-art systems permit the connection of up to 32-element coil arrays [18,19] with even more channels to become available in the near future. The quest for a larger number of independent channels is directly linked to expected improvements in parallel imaging performance with increasing numbers of channels. In addition, the extended coverage of many-element coil arrays can simplify patient setup as the actual positioning of the array on the patient becomes less critical. It has been shown that many-element coil arrays can be efficiently tailored to a specific target anatomy by deriving optimal virtual coil configurations [20].

Besides using parallel imaging principles, interesting opportunities for speeding up data acquisition arise in dynamic imaging. In dynamic imaging of the heart, the object-of-interest is embedded in static or only moderately dynamic structures such as the chest wall and the liver. In addition, individual, neighboring time frames of the heart are often very similar suggesting that considera-

ble information redundancy is present in the data. A range of methods have been proposed attempting to exploit this information redundancy on its own or in conjunction with parallel imaging principles. The assumption of information redundancy being present can be regarded as prior knowledge and, accordingly, these methods are referred to as prior knowledge driven. Several methods have been proposed exploiting prior knowledge without and in combination with parallel imaging principles. Among those methods are UNFOLD [21], k-t BLAST/k-t SENSE [22], Compressed Sensing [23] and FOCUSS [24].

It is the purpose of the article to review the concepts and applications of CMR techniques that exploit zonal, parallel and prior knowledge driven imaging methods to speed up the imaging procedure. To give structure, the paper is subdivided into a concepts section addressing zonal imaging, parallel imaging and prior knowledge driven methods separately followed by an application section presenting preclinical or clinical examples of accelerated CMR imaging including a discussion of relative advantages and limitations of each approach.

Concepts and implementation

In the following a brief review of the underlying principles of undersampling techniques is given. For simplicity and for the fact that most clinically relevant CMR pulse sequences employ rectilinear sampling, only the Cartesian imaging case is considered herein. Accordingly, imaging scan time is assumed to be directly proportional to the number of phase-encode steps. If resolution is kept constant, scan time reduction by reducing the number of phase-encode steps straightforwardly implies a reduction of the field-of-view with resulting foldover artifacts. This is illustrated in Figure 1 for a short-axis slice of the heart reconstructed from 2-fold undersampled data.

Zonal imaging

In zonal or inner volume imaging, foldover artifacts from field-of-view reduction are removed by suppressing any signal from areas outside the reduced field-of-view as shown schematically in Figure 2. This signal suppression can be achieved by choosing perpendicular selection directions of two successive selective radio-frequency pulses in a spin-echo experiment [10]. Alternatively, two-dimensionally selective excitation pulses [5,6,8] or pre-pulses [9] as used with tagging [25] or strain-encoded acquisitions [26] can be used. These approaches are applicable to non spin-echo imaging modes. In general, the net speed gain from zonal imaging is inversely proportional to the extent of the object of interest relative to the full field-of-view. As a consequence, speed-up factors much greater than two cannot be achieved for common anatomy unless individual vessel cross-sections are to be imaged [27]. There are a few advantages of zonal imaging

over other undersampling techniques. First, phased-array receive coils are not necessarily needed and image reconstruction involves the conventional Fourier transforms with its ease and speed. Second, suppressing unwanted signal contributions can be advantageous when using CMR motion tracking and correction techniques as signals from static regions do not interfere. Limitations of multi-dimensionally selective radio-frequency pulse are their relatively long duration and sensitivity to gradient imperfections and off-resonances. Parts of those problems can be alleviated by reducing the pulse duration in combination with the transmit analogue of parallel imaging [28], an area which sees significant research efforts in view of tackling eminent problems at very high magnetic field strengths.

Parallel imaging

The principle of parallel imaging is illustrated in Figure 3. In this exemplary illustration, two-fold undersampling results in two signals overlapping in each point of the image, if reconstructed conventionally. For example, the image intensity seen at point P as indicated in Figure 3 is the superposition of the original signal at that point and the signal folded in from position P' located at position P plus the field-of-view divided by two (= the reduction factor). In order to unfold the signal received in point P additional information is required since two unknowns exist with only their superposition being known. This additional information is derived from the difference in signal reception in the receive coils placed around the object. In the example, two receive coils having different sensitivities with respect to point P and point P' are required. By calculating the difference in sensitivity of coil #1 and coil #2 for all image points, all folded image pixels can be reassigned successively.

The ability to reassign folded signal intensities depends on how many signals actually fold on top of each other. It also depends on the orthogonality of the coil sensitivity functions and thus the coil array configuration and the wave length of the radio-frequency fields involved. The latter is determined by the static magnetic field strength. In the range of clinically available magnets for CMR ranging from 0.5 to 3.0 Tesla, wave length effects play only a minor role. Accordingly, the performance of parallel imaging is primarily determined by the speed-up or reduction factor and the number of coil elements and their configuration with respect to the object being imaged.

Signal-to-noise considerations

In general, and this applies to all undersampling methods, reducing the number of sampled data points reduces the signal-to-noise ratio or, if assuming unit signal response, increases the noise in the reconstructed images by the square-root of the reduction factor. Beyond this


Relation between sampling density and field-of-view. Reduction of sampling density results in signals components folding into the desired, reduced field-of-view.

obvious penalty, noise can further be amplified in parallel imaging due to nonorthogonality of the coil encoding functions as shown in Figure 4 for increasing reduction factors and different imaging modes. This additional noise penalty is reflected by the geometry or g-factor [14] and leads to spatially varying noise across the image. In theoretical work it has been shown that the g-factor relation is governed by fundamental electrodynamics [29,30]. In the 2D Cartesian imaging situation, undersampling can only be applied along one phase-encode direction. The

critical reduction factor is defined as the cut-off beyond which image degradation due to noise increases exponentially. According to Figure 4, the critical factor for 2D imaging is around 3 to 4. It is emphasized again, that the g-factor does not reflect the inherent noise amplification due to scan time reduction; instead it gives the additional noise penalty specific to parallel imaging. In the 3D imaging case, in which undersampling along two dimensions becomes possible [31] the critical limit can be extended to 16, with great prospects for highly accelerated volumetric


Zonal imaging principle. Signals from outside the desired field-of-view are suppressed using appropriate pulse sequence modifications. Accordingly, a foldover-free image is obtained.

3D imaging [19]. Obviously, a prerequisite of applying these large reduction factors is sufficient base signal-to-noise, since 16-fold scan time reduction alone accounts already for 4-fold decrease in the signal-to-noise ratio.


On today's systems, SENSE [14] and GRAPPA [15] are most widely available. The two methods differ in terms of the domain in which the reconstruction problem is solved and with respect to how coil sensitivities are derived. As the image and k-space domains are linked linearly, both SENSE and GRAPPA reconstructions can be formulated

using a common mathematical framework [32]. Using this framework it can be shown that GRAPPA is an approximation of SENSE. A real difference between SENSE and GRAPPA implementations refers to the derivation of coil sensitivities. While in most SENSE implementations coil sensitivities are obtained from a separate low-resolution volumetric reference scan, GRAPPA acquires a few additional profiles encoding the full field-of-view during actual data acquisition for coil calibration. The latter is often referred to as auto-calibration [33]. Auto-calibration has the advantage of ensuring consistency between coil sensitivities and undersampled data if motion during or

Figure 3

Parallel imaging principle. A reduced field-of-view acquisition results in image foldover (for example, point P contains signal from two points (P and P') in the original). The folded image pixels are unfolded using differences in coil sensitivity of the two antenna elements (coil #1 and #2).

in-between successive scans occurs. The downside of this approach is a reduction of the net acceleration factor of the actual acquisition which is smaller than the nominal undersampling factor as reported frequently in the literature.

In dynamic imaging application the undersampling pattern can be shifted as a function of time, such that coil sensitivities can be estimated from a time-average image.

Thereby additional calibration data are not required. The use of interleaved sampling for calibration purposes in combination with SENSE reconstruction is known as TSENSE [34]. In UNFOLD-SENSE [35], the time-interleaving and with it the displacement of static and dynamic signal distributions in the field-of-view along the temporal frequency axis is exploited in combination with a temporal filter to reduce residual artifact from parallel image reconstruction.

Figure 4

Noise amplification in parallel imaging. The additional noise penalty as expressed by the geometry factor is shown. For 2D imaging the critical reduction factor is between 3-4. In 3D imaging with comparable dimensions and reduction factors along the two phase-encode axes the critical reduction factor is about 16.

Prior knowledge driven imaging

Prior knowledge driven methods are used to exploit information redundancy in dynamic CMR image series. The amount of redundant information present in a cine series of the heart can be visualized using a simple experiment. If one subtracts the time-average over image points from each individual image frame of the series, only the differences between the average and each individual frame

remain (Figure 5). Those difference maps represent the image in a much sparser fashion with only very few significant signals. Based on this observation, undersampling and reconstruction strategies have been developed to accelerate the acquisition of objects presenting with large degrees of image correlation in space and time such as the heart embedded in the chest.

Several strategies have emerged using image correlations as such [21,22,24,36,37] or in combination with parallel imaging [22,24,35,38]. Widely available on current systems are k-t BLAST and k-t SENSE [22]. These methods also use a time-interleaved, shifted undersampling similar to the TSENSE and UNFOLD-SENSE methods. Besides sequentially shifted undersampling, optimal undersam-pling pattern can be derived [39]. Like in TSENSE, auto-calibration for parallel imaging from a time-average image is applicable.

In k-t BLAST and k-t SENSE an image estimate is obtained in a so-called training stage. In this training stage, images are acquired at very low spatial but at full temporal resolution. These data may either be obtained in a separate scan or in an interleaved fashion during the actual acquisition [40,41]. During the undersampling stage, sparse sampling is applied along k-space and time (k-t imaging). To obtain a high-resolution, unaliased image series, the prior knowledge as derived from the training data is used to unfold the image as outlined in Figure 6. In parlance of image reconstruction, the training data determines a filter which suppresses folding artifacts as a result of undersam-pling. While this reconstruction filter is the only handle to

Cine series (a)

Difference (a-b)

Time-average (b)

Figure 5

Information redundancy. In this example the time-average of all cine image frames is subtracted from each original image frame to illustrate the sparseness of information contained in the cine series. It is seen that the difference images contains only very few image pixels with significant values.


k-t BLAST principle. Based on low-resolution training data, an estimate of the expected signal intensities of the object is obtained and used to unfold the aliased components resulting from undersampled data acquisition.

resolve foldover in k-t BLAST, coil sensitivities are used in addition in k-t SENSE.

In general, k-t SENSE tends to outperform its frame-by-frame SENSE counterpart for comparable reduction factors (Figure 7), even if state-of-the-art 32-channel coil arrays are employed. It is, however, to be noted that k-t SENSE is only applicable to cine or dynamic applications, while parallel imaging applies to all kinds of acquisition strategies.

Limitations of k-t BLAST and k-t SENSE can include insufficient un-mixing of different temporal frequency compo-

nents folded on top of each other which can lead to temporal blurring and flickering in the data. To alleviate the effect of temporal blurring at high reduction factors, training data plug-in has been proposed [42]. Thereby training data are substituted into the final k-space reconstruction ensuring data consistency. A comparison of blood velocities obtained in the ascending aorta with up to 8x net reduction using k-t SENSE with training data plug-in relative to SENSE/TSENSE reconstructions is shown in Figure 8. Superior noise suppression with k-t methods at high reduction factors is traded for pronounced temporal filtering as exemplified for 10x net reduction. Temporal filtering seen at at high reduction fac-

8xfr-f SENSE 8x SENSE

12x/c-f SENSE 12x SENSE


Comparison of k-t SENSE and SENSE at high acceleration using a 32-element coil array. Short axis views obtained from 2D imaging at nominal reduction factors of 8 and 12 are compared (net reduction factors: 5.6 and 7.3). Spatiotemporal plots were taken along the dashed line indicated.

tors can be reduced by using an increased number of training data profiles [43] which, however, offsets the initial speed gain. As an alternative, parallel imaging can be employed to boost the training data resolution without time penalty. To this end, however, careful considerations as to the minimum acceptable SNR in the training data are required. The low resolution in the training data also demands dedicated SENSE reconstruction methods to avoid residual folding artifacts [44].

Signal-to-noise considerations

As for parallel imaging, image noise is non-uniform across the image with k-t methods. In parallel imaging, this non-uniformity is directly accessible from g-factor maps and

depends on the reduction factor and coil array configuration. In k-t methods, the noise also depends on object dynamics. For example, in static or slowly moving structures, the noise variance will be much smaller compared to areas with highly dynamic features. This noise non-uniformity found with parallel imaging and k-t methods has implications for practical signal-to-noise ratio measurements. To this end methods have been proposed [45,46]. However, difference methods [45] may render inaccurate results in particular in cardiac imaging, because truly identical object states are difficult to capture with successive images. Instead of using a difference method it is proposed to repeat the sequence with all radio-frequency pulse switched off or obtain noise only data during a short


Comparison of TSENSE and k-t SENSE reconstructions of data decimated to different reduction factors from a fully sampled phase-contrast flow measurement acquired with a 32-channel coil array. Aortic velocity data as a function of cardiac phase from TSENSE and k-t SENSE with net reduction factors of 2, 3, 4 and 6 are compared relative to values obtained from the fully sampled reference (left, middle). Training data plug-in is used in all k-t SENSE reconstructions. Significant temporal filtering is seen for net reduction beyond factor 8 with k-t SENSE (right).

prescan [46]. This does not require respiratory control or breathholding as spatial information is not encoded and would provide together with the signal data true signal-to-noise ratio assessment. Such a scheme is currently offered by at least one major vendor.

Accelerations achievable with parallel and prior knowledge driven CMR

A number of parallel imaging and prior knowledge driven methods as well as hybrid methods of both concepts have been proposed and applied to CMR. In Table 1 examples of these approaches are given with references for their practical CMR application as well as nominal and net reduction factors reported. While nominal reduction factors specify the amount of undersampling without taking into account data sampled to populate auto-calibration or training data, net reduction factors denote effective reduction factors. Unfortunately, net reduction factors are not always reported in the literature and actual scanner implementations may differ from what is described in publications. Accordingly, for Table 1, assumptions had to be made to estimate net reduction factors. For methods with external calibration, nominal and net reduction factors are assumed to be identical.

Critical reduction factors beyond which the diagnostic value of reconstructed images is compromised have been indicated and are largely in agreement with theoretical considerations and experimental g-factor measurements for parallel imaging methods [29,30]. In 2D imaging applications using SENSE, TSENSE, GRAPPA, TGRAPPA this critical factor amounts to about four. In 3D imaging with large coil arrays, greater reduction factors have been

achieved [47,48] but the critical value with respect to diagnostic image quality achievable in CMR applications remains to be determined. Using prior knowledge driven methods in combination with parallel imaging, net reduction factors up to 16 have been realized [24]. The diagnostic value of images acquired at such high reduction factors remains, however, to be shown. CMR studies using k-t BLAST and k-t SENSE have indicated temporal filtering effects for net reduction factors greater than eight compromising the true temporal resolution in dynamic data.

Clinical application

Parallel imaging and k-t methods that are available on clinical scanners, such as SENSE, GRAPPA, TSENSE, k-t SENSE and k-t BLAST, provide operators with flexible and important tools to optimize a range of CMR protocols. The increased scanning efficiency afforded by these methods can be used to shorten imaging times especially in patients with poor breathhold capacity with resulting improvements in image quality. Alternatively, spatial and temporal resolution as well as coverage can be increased.

Cine imaging of cardiac function

Cine imaging is an essential component of virtually every CMR study. Measurements of regional and global ventricular function are typically derived from a stack of multiple 2D cines. Depending on the clinical application, further cine series in customized imaging planes are often performed. With current balanced steady-state free-precession pulse sequences the acquisition time for each 2D slice is on the order of 10 seconds to achieve acceptable spatial and temporal resolution [49]. Accordingly, cine imaging accounts for a large proportion of the scan time

Table 1:

Parallel Imaging Prior Knowledge Hybrids Prior Knowledge CMR Rnom/Rnet References

Driven Imaging Applications

SENSE cine, rt, dce, lge, 1.5-8.0 [19,55,64,66-

qflow, cmra 68,73,78,82-84]

SMASH - cine, cmra 2.0 [16,85-87]

GRAPPA - cine, cmra 2.0/1.6 [15,53,88]

UNFOLD constant x-f cine, rt, dce 2.0 [21,89-92]


k-t BLAST Measured x-f cine, rt, dce, qflow 5.0-8.0/4.0- 5.7 [22,39,40,42,51,57,6

support 9,70,93]

CS sparse x-f support cine, qflow 2.0-7.0 [23,24,94]

TCR temporal constraint dce 4.0-5.0 [37]

TSENSE temporal cine, rt, dce, dense 2.0-7.0 [34,52,65,75,91,95-

correlation1) 98]

TGRAPPA temporal cine 2.0-4.0 [99]


k-t GRAPPA temporal cine, cmra 3.0-7.0/2.0- 5.2 [38,100 ]


k-t SENSE Measured x-f cine, rt, dce, qflow 5.0-8.0/4.0- 5.7 [22,39,41,42]


FOCUSS measured, sparse x- cine 4.0-16.0 [24]

f support

UNFOLD-SENSE constant x-f cine, dce 2.0-3.0/1.8- 2.7 [35]


'icoil sensitivity weights from temporal average of data sampled on a sheared grid (cmp. DC signal in x-f support)

2)auto-calibration signals (acs) from temporal average of data sampled on a sheared grid

3)with or without acs lines from temporal average of data sampled on a sheared grid x: spatial position, f: temporal frequency

cine: cardiac-gated cine imaging; rt: non-gated real-time imaging; dce: dynamic contrast-enhanced imaging; lge: late gadolinium enhancement; qflow: quantitative flow imaging; cmra: coronary magnetic resonance angiography; dense: displacement-encoded imaging

Rnom: nominal reduction factor; Rnet: net reduction factor taking into account calibration or training data (only reduction factors used in 2D CMR applications are reported here; in case net reduction factors are not available from literature a matrix of 200 in phase-encode direction and 11 auto-calibration or training profiles acquired interleaved are assumed).

in most CMR studies. Therefore, parallel imaging methods with two to three-fold acceleration are often used in clinical cine protocols (Figure 9) reducing the scan time per slice to about 4-5 seconds. As expected from theory, beyond three-fold acceleration, SNR limitations increasingly affect the clinical utility of the acquired data. Several studies have shown that cine imaging with SENSE and other parallel imaging methods yields comparable measurements of left and right ventricular volumes and function as conventional cine imaging [50].

Further acceleration can be achieved with k-t methods reducing acquisition times for a 2D cine scan to 2-3 seconds per slice. As seen in Figure 9, image quality is better preserved at higher acceleration than with parallel imaging alone. Initial studies have suggested that 2D cine imaging with k-t BLAST [51] and TSENSE [52] provide similar measurements of cardiac function as non-accelerated acquisition as long as the net acceleration factor does not exceed four to five. For larger reduction factors, underestimation of ejection fraction has been noted [52].

Parallel imaging methods also allow for so called "realtime" cine imaging, in which multi-phase data sets are acquired with a continuous dynamic acquisition [53]. While such methods have not yet found a role for the assessment of cardiac function and dimensions, potential applications of 2D "real-time" cine imaging include for example the study of the motion of the interventricular septum in suspected constrictive pericardial disease [54] (Figure 10).

A different use of increased scan efficiency is time-resolved 3D imaging. In particular in combination with many-element coil arrays, parallel imaging can provide sufficient acceleration to image the entire heart within a single breathhold at spatial and temporal resolutions comparable to conventional 2D cine methods [55]. Prior knowledge driven techniques can also be employed to facilitate single-breathhold 3D cine imaging [40,51]. The clinical role of 3D cine imaging is yet to be determined, but its potential benefits are easier planning, faster data acquisition and avoidance of misregistration that can affect 2D


a^h éh^hí j

2x/c-fSENSE 3xBÍSENSE 4xfr-fSENSE 6x/f-fSENSE


Comparison of clinical 2D cine imaging in the short axis orientation acquired with SENSE and k-t SENSE using a 6-element coil array. Each image depicts the end-diastolic frame of separate breath-hold cines acquired with 1* to 3* SENSE and 2* to 6* k-t SENSE. 2* SENSE provides excellent image quality. With 3* SENSE, the SNR reduction and image artifacts become apparent (arrows). k-t SENSE up to 4-fold net acceleration results in good image quality. At 6* k-t SENSE, slight blurring of the endocardial wall and papillary muscles is visible (arrow).


Real-time imaging with 2* SENSE in the short axis orientation in a patient with recent mitral valve replacement surgery. Images captured during deep expiration (left) and inspiration (right) show diastolic inversion of the interven-tricular septum in inspiration consistent with constrictive pericardial disease.

imaging with multiple breathholds. For volumetric measurements, initial studies have suggested a small systematic underestimation in particular of the end-diastolic volume in 3D compared with 2D cine imaging, which may be related to the effects of misregistration or to the encoding process [51,56].

Three-dimensional cine imaging has also been successfully applied to dobutamine-stress MRI, where cardiac function has to be assessed repeatedly and rapidly at different stages of pharmacological stress. Segmental wall motion analysis from k-t BLAST 3D cine imaging was similar to conventional cine SSFP imaging at rest and during dobutamine stress [57].

Real-time imaging

Besides its use in interventional applications [58-60], advanced acquisition schemes combined with rapid reconstruction also permit the acquisition of CMR data in real-time [61-63]. Real-time methods are of great importance for imaging patients with irregular sinus rhythm and poor breathhold capability. There has been a number of studies indicating the benefits of improved spatiotemporal resolution by using parallel imaging methods for assessing functional indices [53,64,65], diagnosing wall motion abnormalities [66] or quantifying blood flow velocities [67] using real-time methods.

Velocity-encoded imaging

Accelerated quantitative flow imaging has been of interest in a range of applications. Given the small diameter of blood vessel, zonal imaging techniques have been used to facilitate real-time flow quantification during hyperemic response [7,27]. Flow imaging accelerated with parallel imaging techniques has been used for shunt quantification in children [68]. The speed gain with k-t BLAST has also been demonstrated to reduce very long scan time associated with velocity spectrum mapping or Fourier velocity imaging [69]. Application for velocity quantification in stenotic valves has been described [70]. Quantitative flow data obtained with k-t SENSE have been shown to agree well with the reference up to reduction factors of eight [42]. Beyond net reduction factors of 8 temporal filtering of blood velocity curves is noticeable in accordance to the data presented in Figure 8.

Coronary and anatomical imaging

The reliable non-invasive visualization of the coronary arterial system remains one of the most rewarding prospects of tomographic imaging. Although coronary computed tomography currently provides superior image quality, CMR coronary angiography remains an important alternative because it does not require the use of contrast agents and does not expose patients to ionizing radiation. Current clinical applications for coronary MRI

are the assessment of anomalous coronary arteries and other congenital pathologies [1]. Imaging of the coronary venous system for the planning of electrophysiological studies is emerging as a potential future application.

Accelerated data acquisition is essential to permit the acquisition of a 3D imaging slab that covers the entire heart with the required spatial resolution and within an acceptable time frame. So called "whole-heart" CMR with typically 90 to 120 slices and a spatial resolution of approximately 1 x 1 x 1 mm3 is usually acquired with the help of respiratory navigators [71]. Parallel imaging is used in most whole-heart protocols in order to keep imaging times below 10 minutes. Further acceleration is feasible, especially with larger coil arrays [19]. As with all applications of parallel imaging, SNR and image quality are offset against the advantage of shorter acquisition time as demonstrated in Figure 11.

The same "whole-heart" imaging methods developed for coronary imaging can be employed for anatomical assessment of other cardiac structures. The high spatial resolution of the acquired 3D data sets provides a useful basis for assessment of congenital heart disease and other morphological studies [72].

Dynamic imaging of myocardial perfusion

The principal challenges of dynamic first pass myocardial perfusion imaging are to acquire data with high temporal resolution to suppress motion artifacts, to provide sufficient spatial coverage to delineate ischemia in all coronary segments and to deliver sufficient spatial resolution to minimize imaging artifacts such as endocardial dark rims. In order to satisfy all of the requirements, gradient echo perfusion CMR has been combined with echo planar imaging and SENSE [73,74] or TSENSE [75]. With these methods, three to four imaging sections can typically be acquired at every heart beat in an acquisition time of around 100 ms per slice and an in-plane spatial resolution of 2-3 mm [73,74]. Recently, prior-knowledge driven myocardial perfusion imaging using k-t BLAST and k-t SENSE has been proposed. The additional acceleration has been invested in either a shorter acquisition time per slice [57] or improved spatial resolution [41]. In particular the latter approach appears to be promising, as imaging artifacts such as endocardial dark rims, were reduced in line with the spatial resolution and the ability to assess transmural perfusion distribution was enhanced (Figure 12).

Imaging of myocardial viability

Late gadolinium-enhanced (LGE) CMR for the assessment of myocardial viability has become a cornerstone in the diagnosis of ischemic and non-ischemic cardiomyopathies [76]. The method provides unparalleled delineation


Free-breathing, whole-heart coronary MRA with SENSE using a 32-element coil array. The path of the right coronary artery has been tracked. Up to 4-fold acceleration SNR losses are acceptable. At 8-fold acceleration noise enhancement compromises diagnostic quality.

of fibrotic or destructive processes in the myocardium with high contrast to normal tissue. However, each section of a 2D stack of late contrast-enhanced data has an acquisition time of 10-12 seconds [76], so that coverage of the whole heart can require up to 10 minutes of scan time. Because of the typically high contrast between normal myocardium and scar, LGE is well suited for acceleration with parallel imaging methods with acceleration factors of 2 and 3 in 2D and 3D imaging protocols [77]. As can be seen in Figure 13, image quality and diagnostic

yield are comparable between non-accelerated and 2x SENSE imaging. Image-based motion correction techniques in combination with single-shot SENSE have been developed to improve signal- and contrast-to-noise ratios in LGE by averaging [78].


Accelerated CMR using parallel imaging and prior knowledge driven methods are increasingly utilized to improve scanning efficiency and image quality of 2D and 3D appli-

5x k-t SENSE


Adenosine stress myocardial perfusion imaging with an acquired in-plane spatial resolution of 1.3 x 1.3 mm2 facilitated by 5x k-t SENSE. Four short axis slices are shown covering the heart from the apex (left) to the base (right). A large infero-lateral perfusion defect and anteroseptal ischemia in the apex is seen. Due to the very high spatial resolution, endocardial dark rim artifacts are minimal and the spatial extent of the perfusion defects is clearly delineated.



Comparison of 2D late-contrast enhanced MR imaging without and with 2* SENSE acceleration and using a 6-element coil array. Both images provide similar image quality with clear delineation of an antero-septal scar. The SNR decrease is seen with 2* SENSE, best appreciated in the lateral wall.

cations. Many clinical 2D protocols make already use of these methods with net accelerations of up to four-fold for parallel imaging and up to five-fold for prior knowledge driven methods.

Speed-up techniques continue to develop. The wider availability of 32-channel MR systems will foster studies into the practical limits. Clearly, 3D imaging techniques will be at the focus according to theoretical predictions of noise amplification in parallel imaging. It is foreseen that the range of practically applicable reduction factors will extend to three to four in 2D imaging while large 3D volumetric imaging will be accelerated by up to factors of eight to ten [48]. Beyond those limits, the resulting signal-to-noise ratio will inherently set the bound. With the advent of many-element coil arrays careful considerations of coil calibration are required as individual coil elements are small with considerable changes in sensitivity as a function of space. It is likely that the resolution of calibration data needs to be improved leading to additional time to be spent thereby offsetting the initial speed-gain if auto-calibration is used. When acquiring separate calibration data as with SENSE implementations, the effect of spatial misregistration due to motion in-between the calibration and actual data is amplified prompting for some precautions or registration steps.

On the other side, a very practical advantage of using large coil arrays with their inherent coverage is the relative insensitivity to misplacement of the coil on the patient's chest. This simplifies patient setup and avoids coil readjustment for different target territories.

The field of prior knowledge driven methods currently enjoys vivid research activities. Several modifications to the original k-t BLAST/k-t SENSE scheme have been proposed addressing some of the current shortcomings [24,79]. Another very interesting and potential approach uses the Compressed Sensing framework [80,81]. In Compressed Sensing data are undersampled randomly leading to incoherent aliasing. Using iterative reconstruction methods, the aliasing can be resolved. It has been shown that such an approach is well suited for speeding up object data that are particularly sparse in the image domain itself as it is case with angiography or using a suitable sparsify-ing transformation [23].

Ultimately, speed-up techniques are required to provide the most appropriate balance between noise amplification and image fidelity for a given reduction factor. It is the purpose of ongoing and future research to determine the faithful range of reduction factors in combination with a suitable encoding and reconstruction strategy for

all required CMR applications. It has to be a clear aim to provide unifying assessment of image quality and fidelity for all the different acceleration techniques available to allow comparisons and robust predictions.

In conclusion, parallel imaging and prior knowledge driven methods have become very valuable tools to streamline and simplify CMR exams. Their impact is expected to develop rapidly over the next years with the wider availability of multi-channel MR systems. Several methods are available which differ in their performance depending on application requiring careful considerations when setting up imaging protocols.

Competing interests

The authors declare that they have no competing interests. Acknowledgements

The authors would like to acknowledge contributions and support from: Profs. Peter Boesiger and Klaas P. Pruessmann, Martin Buehrer and Viton Vitanis, Institute for Biomedical Engineering, University and ETH Zurich, Switzerland; Jeffrey Tsao, Novartis Institutes of Biomedical Research Inc., Cambridge, MA, USA; and Philips Healthcare, Best, The Netherlands.


1. Pennell DJ, Sechtem UP, Higgins CB, Manning WJ, Pohost GM, Rade-makers FE, van Rossum AC, Shaw LJ, Yucel EK: Clinical indications for cardiovascular magnetic resonance (CMR): Consensus panel report. Journal of Cardiovascular Magnetic Resonance 2004, 6(4):727-765.

2. Gershlick AH, de Belder M, Chambers J, Hackett D, Keal R, Kelion A, Neubauer S, Pennell DJ, Rothman M, Signy M, Wilde P: Role of non-invasive imaging in the management of coronary artery disease: an assessment of likely change over the next ! 0 years. A report from the British Cardiovascular Society Working Group. Heart 2007, 93(4):423-43l.

3. Schwitter J, Wacker C, Wilke N, Ai-Saadi N, Berlin C, Hoebel N, Simor T: Magnetic resonance imaging for myocardial perfusion assessment in coronary artery disease trial (MR-IMPACT II): A phase III multicenter, multivendor trial comparing perfusion cardiac magnetic resonance versus single photon emission computed tomography for the detection of coronary artery disease. Circulation 2006, !!4(!8):806-806.

4. Feinberg DA, Hoenninger JC, Crooks LE, Kaufman L, Watts JC, Arakawa M: Inner Volume Mr Imaging - Technical Concepts and Their Application. Radiology 1985, !56(3):743-747.

5. Yang GZ, Gatehouse PD, Keegan J, Mohiaddin RH, Firmin DN: Three-dimensional coronary MR angiography using zonal echo planar imaging. Magnetic Resonance in Medicine 1998, 39(5):833-842.

6. Yang GZ, Burger P, Gatehouse PD, Firmin DN: Locally focused 3D coronary imaging using volume-selective RF excitation. Magnetic Resonance in Medicine 1999, 4!(!): 171-178.

7. Mohiaddin RH, Gatehouse PD, Moon JCC, Youssuffidin M, Yang GZ, Firmin DN, Pennell DJ: Assessment of reactive hyperaemia using real time zonal echo-planar flow imaging. Journal of Cardiovascular Magnetic Resonance 2002, 4(2):283-287.

8. Oelhafen M, Schwitter J, Kozerke S, Luechinger R, Boesiger P: Assessing arterial blood flow and vessel area variations using real-time zonal phase-contrast MRI. Journal of Magnetic Resonance Imaging 2006, 23(3):422-429.

9. Pan L, Stuber M, Kraitchman DL, Fritzges DL, Gilson WD, Osman NF: Real-time imaging of regional myocardial function using Fast-SENC. Magnetic Resonance in Medicine 2006, 55(2):386-395.

10. Gamper U, Boesiger P, Kozerke S: Diffusion imaging of the in vivo heart using spin echoes-considerations on bulk motion sensitivity. Magnetic Resonance in Medicine 2007, 57(2):33l-337.

11. Carlson JW: An Algorithm for Nmr Imaging Reconstruction Based on Multiple Rf Receiver Coils. Journal of Magnetic Resonance 1987, 74(2): 376-380.

12. Hutchinson M, Raff U: Fast Mri Data Acquisition Using Multiple Detectors. Magnetic Resonance in Medicine 1988, 6(1 ):87-9l.

13. Sodickson DK, Manning WJ: Simultaneous acquisition of spatial harmonics (SMASH): fast imaging with radiofrequency coil arrays. Magn Reson Med 1997, 38(4):59l-603.

14. Pruessmann KP, Weiger M, Scheidegger MB, Boesiger P: SENSE: sensitivity encoding for fast MRI. Magn Reson Med l999, 42(5):952-962.

15. Griswold MA, Jakob PM, Heidemann RM, Nittka M, Jellus V, Wang J, Kiefer B, Haase A: Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn Reson Med 2002, 47(6): 1202-1210.

16. Griswold MA, Jakob PM, Edelman RR, Sodickson DK: A multicoil array designed for cardiac SMASH imaging. Magma 2000, 10(2): 105-113.

17. Weiger M, Pruessmann KP, Osterbauer R, Bornert P, Boesiger P, Jez-zard P: Sensitivity-encoded single-shot spiral imaging for reduced susceptibility artifacts in BOLD fMRI. Magn Reson Med 2002, 48(5):860-866.

18. Hardy CJ, Cline HE, Giaquinto RO, Niendorf T, Grant AK, Sodickson DK: 32-Element receiver-coil array for cardiac imaging. Magnetic Resonance in Medicine 2006, 55(5): 1142-1 149.

19. Niendorf T, Hardy CJ, Giaquinto RO, Gross P, Cline HE, Zhu YD, Kenwood G, Cohen S, Grant AK, Joshi S, Rofsky NM, Sodickson DK: Toward single breath-hold whole-heart coverage coronary MRA using highly accelerated parallel imaging with a 32-channel MR system. Magnetic Resonance in Medicine 2006, 56(1): 167-176.

20. Buehrer M, Pruessmann KP, Boesiger P, Kozerke S: Array compression for MRI with large coil arrays. Magnetic Resonance in Medicine 2007, 57(6): 1131-1 139.

21. Madore B, Glover GH, Pelc NJ: Unaliasing by fourier-encoding the overlaps using the temporal dimension (UNFOLD), applied to cardiac imaging and fMRI. Magn Reson Med 1999, 42(5):8l3-828.

22. Tsao J, Boesiger P, Pruessmann KP: k-t BLAST and k-t SENSE: dynamic MRI with high frame rate exploiting spatiotemporal correlations. Magn Reson Med 2003, 50(5): 1031-1042.

23. Lustig M, Donoho D, Pauly J: Sparse MRI: The application of compressed sensing for rapid MR imaging. Magn Reson Med 2007, 58(6): 1182-1195.

24. Jung H, Ye JC, Kim EY: Improved k-t BLAST and k-t SENSE using FOCUSS. Physics in Medicine and Biology 2007, 52(1 l):320l-3226.

25. Fischer SE, Mckinnon GC, Scheidegger MB, Prins W, Meier D, Boesiger P: True Myocardial Motion Tracking. Magnetic Resonance in Medicine 1994, 31 (4):40l-4l3.

26. Osman NF: Detecting stiff masses using strain-encoded (SENC) imaging. Magnetic Resonance in Medicine 2003, 49(3):605-608.

27. Schwitter J, Oelhafen M, Wyss BM, Kozerke S, Amann-Vesti B, Luscher TF, Boesiger P: 2D-spatially-selective real-time magnetic resonance imaging for the assessment of microvascular function and its relation to the cardiovascular risk profile. Journal of Cardiovascular Magnetic Resonance 2006, 8(5):759-769.

28. Katscher U, Bornert P, Leussler C, Brink JS van den: Transmit SENSE. Magnetic Resonance in Medicine 2003, 49(1): 144-l 50.

29. Ohliger MA, Grant AK, Sodickson DK: Ultimate intrinsic signal-to-noise ratio for parallel MRI: electromagnetic field considerations. Magn Reson Med 2003, 50(5): l0l8-1030.

30. Wiesinger F, Boesiger P, Pruessmann KP: Electrodynamics and ultimate SNR in parallel MR imaging. Magnetic Resonance in Medicine 2004, 52(2):376-390.

3 l. Weiger M, Pruessmann KP, Boesiger P: 2D SENSE for faster 3D MRI. Magma 2002, 14(1): l0-l9.

32. Pruessmann KP: Encoding and reconstruction in parallel MRI. Nmr in Biomedicine 2006, l9(3):288-299.

33. Jakob PM, Griswold MA, Edelman RR, Sodickson DK: AUTOSMASH: a self-calibrating technique for SMASH imaging. SiMultaneous Acquisition of Spatial Harmonics. Magma l998, 7(l):42-54.

34. Kellman P, Epstein FH, McVeigh ER: Adaptive sensitivity encoding incorporating temporal filtering (TSENSE). Magn Reson Med 2001, 45(5):846-852.

35. Madore B: UNFOLD-SENSE: A parallel MRI method with self-calibration and artifact suppression. Magnetic Resonance in Medicine 2004, 52(2):310-320.

36. Tsao J, Behnia B, Webb AG: Unifying linear prior-information-driven methods for accelerated image acquisition. Magn Reson Med 2001, 46(4):652-660.

37. Adluru G, Awate SP, Tasdizen T, Whitaker RT, DiBella EVR: Temporally constrained reconstruction of dynamic cardiac perfusion MRI. Magnetic Resonance in Medicine 2007, 57(6): 1027-1036.

38. Huang F, Akao J, Vijayakumar S, Duensing GR, Limkeman M: k-t GRAPPA: A k-space implementation for dynamic MRI with high reduction factor. Magnetic Resonance in Medicine 2005, 54(5): 1 172-1184.

39. 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. Magnetic Resonance in Medicine 2005, 53(6): 1372-1382.

40. Kozerke S, Tsao J, Razavi R, Boesiger P: Accelerating cardiac Cine 3D Imaging using k-t BLAST. Magnetic Resonance in Medicine 2004, 52(1): 19-26.

41. Plein S, Ryf S, Schwitter J, Radjenovic A, Boesiger P, Kozerke S: Dynamic contrast-enhanced myocardial perfusion MRI accelerated with k-t SENSE. Magnetic Resonance in Medicine 2007, 58(4):777-785.

42. Baltes C, Kozerke S, Hansen MS, Pruessmann KP, Tsao J, Boesiger P: Accelerating cine phase-contrast flow measurements using k-t BLAST and k-t SENSE. Magnetic Resonance in Medicine 2005, 54(6): 1430-1438.

43. Hansen MS, Kozerke S, Pruessmann KP, Boesiger P, Pedersen EM, Tsao J: On the influence of training data quality in k-t BLAST reconstruction. Magnetic Resonance in Medicine 2004, 52(5): 1 175-1183.

44. Sanchez-Gonzalez J, Tsao J, Dydak U, Desco M, Boesiger P, Pruessmann KP: Minimum-norm reconstruction for sensitivity-encoded magnetic resonance spectroscopic imaging. Magnetic Resonance in Medicine 2006, 55(2):287-295.

45. Reeder SB, Wintersperger BJ, Dietrich O, Lanz T, Greiser A, Reiser MF, Glazer GM, Schoenberg SO: Practical approaches to the evaluation of signal-to-noise ratio performance with parallel imaging: Application with cardiac imaging and a 32-channel cardiac coil. Magnetic Resonance in Medicine 2005, 54(3):748-754.

46. Kellman P, McVeigh ER: Image reconstruction in SNR units: A general method for SNR measurement. Magnetic Resonance in Medicine 2005, 54(6): 1439-1447.

47. Zhu YD, Hardy CJ, Sodickson DK, Giaquinto RO, Dumoulin CL, Kenwood G, Niendorf T, Lejay H, McKenzie CA, Ohliger MA, Rofsky NM: Highly parallel volumetric imaging with a 32-element RF coil array. Magnetic Resonance in Medicine 2004, 52(4):869-877.

48. Niendorf T, Sodickson DK: Parallel imaging in cardiovascular MRI: methods and applications. Nmr in Biomedicine 2006, I9(3):325-34I.

49. Alfakih K, Plein S, Thiele H, Jones T, Ridgway JP, Sivananthan MU:

Normal human left and right ventricular dimensions for MRI as assessed by turbo gradient echo and steady-state free precession imaging sequences. Journal of Magnetic Resonance Imaging 2003, I7(3):323-329.

50. Kacere RD, Pereyra M, Nemeth MA, Muthupillai R, Flamm SD: Quantitative assessment of left ventricular function: Steady-state free precession MR imaging with or without sensitivity encoding. Radiology 2005, 235(3): 1031-1035.

51. Jahnke C, Nagel E, Gebker R, Bornstedt A, Schnackenburg B, Kozerke S, Fleck E, Paetsch I: Four-dimensional single breathhold magnetic resonance imaging using kt-BLAST enables reliable assessment of left- and right-ventricular volumes and mass. Journal of Magnetic Resonance Imaging 2007, 25(4):737-742.

52. Wintersperger BJ, Reeder SB, Nikolaou K, Dietrich O, Huber A, Greiser A, Lanz T, Reiser MF, Schoenberg SO: Cardiac CINE MR Imaging with a 32-channel cardiac coil and parallel imaging: Impact of acceleration factors on image quality and volumetric accuracy. Journal of Magnetic Resonance Imaging 2006, 23(2):222-227.

53. Wintersperger BJ, Nikolaou K, Dietrich O, Rieber J, Nittka M, Reiser MF, Schoenberg SO: Single breath-hold real-time cine MR

imaging: improved temporal resolution using generalized autocalibrating partially parallel acquisition (GRAPPA) algorithm. Eur Radiol 2003, 13(8): 1931-1936.

54. Francone M, Dymarkowski S, Kalantzi M, Bogaert J: Real-time cine MRI of ventricular septal motion: A novel approach to assess ventricular coupling. Journal of Magnetic Resonance Imaging 2005, 21 (3):305-309.

55. Niendorf T, Saranathan M, Lingamneni A, Pedrosa I, Spencer M, Cline H, Foo TKF, Rofsky NM: Short breath-hold, volumetric coronary MR angiography employing steady-state free precession in conjunction with parallel imaging. Magnetic Resonance in Medicine 2005, 53(4):885-894.

56. Greil GF, Germann S, Kozerke S, Baltes C, Tsao J, Urschitz MS, Seeger A, Tangcharoen T, Bialkowsky A, Miller S, Sieverding L: Assessment of left ventricular volumes and mass with fast 3D cine steady-state free precession k-t space broad-use linear acquisition speed-up technique (k-t BLAST). Journal of Magnetic Resonance Imaging 2008, 27(3):5I0-515.

57. Jahnke C, Paetsch I, Gebker R, Bornstedt A, Fleck E, Nagel E: Accelerated 4D dobutamine stress MR imaging with k-t BLAST: Feasibility and diagnostic performance. Radiology 2006, 24I(3):718-728.

58. Guttman MA, Lederman RJ, Sorger JM, McVeigh ER: Real-time volume rendered MRI for interventional guidance. Journal of Cardiovascular Magnetic Resonance 2002, 4(4):43 I-U43I.

59. Lederman RJ, Guttman MA, Peters DC, Thompson RB, Sorger JM, Dick AJ, Raman VK, McVeigh ER: Catheter-based endomyocardial injection with real-time magnetic resonance imaging. Circulation 2002, 105(11): 1282-1284.

60. McVeigh ER, Guttman MA, Lederman RJ, Li M, Kocaturk O, Hunt T, Kozlov S, Horvath KA: Real-time interactive MRI-guided cardiac surgery: Aortic valve replacement using a direct apical approach. Magnetic Resonance in Medicine 2006, 56(5):958-964.

61. Setser RM, Fischer SE, Lorenz CH: Quantification of left ventricular function with magnetic resonance images acquired in real time. Journal of Magnetic Resonance Imaging 2000, I2(3):430-438.

62. Klein C, Nagel E, Schalla S, Lehmkuhl H, Schnackenburg B, Bornstedt A, Fleck E: Quantification of flow velocity and volume with real time cardiovascular magnetic resonance. Journal of the American College of Cardiology 2000, 35(2):454a-454a.

63. Schalla S, Nagel E, Klein C, Paetsch I, Bonrstedt A, Schnackenburg B, Lehmkuhl H, Fleck E: Real-time high-dose dobutamine stress magnetic resonance imaging for the detection of left ventricular wall motion abnormalities in patients with coronary artery disease. European Heart Journal 2000, 21:570-570.

64. Weiger M, Pruessmann KP, Boesiger P: Cardiac real-time imaging using SENSE. Magnetic Resonance in Medicine 2000, 43(2): 177-184.

65. Yamamuro M, Tadamura E, Kanao S, Okayama S, Okamoto J, Urayama S, Kimura T, Komeda M, Kita T, Togashi K: Cardiac functional analysis by free-breath real-time cine CMR with a spatiotemporal filtering method, TSENSE: Comparison with breath-hold cine CMR. Journal of Cardiovascular Magnetic Resonance 2006, 8(6):80I-807.

66. Strach K, Meyer C, Schild H, Sommer T: Cardiac stress MR imaging with dobutamine. European Radiology 2006, 16(1 2):2728-2738.

67. Nezafat R, Kellman P, Derbyshire JA, McVeigh ER: Real-time blood flow imaging using autocalibrated spiral sensitivity encoding. Magnetic Resonance in Medicine 2005, 54(6): I557-I56I.

68. Beerbaum P, Korperich H, Gieseke J, Barth P, Peuster M, Meyer H: Rapid left-to-right shunt quantification in children by phase-contrast magnetic resonance imaging combined with sensitivity encoding. Circulation 2003, 108(11): I355-I36I.

69. Hansen MS, Baltes C, Tsao J, Kozerke S, Pruessmann KP, Boesiger P, Pedersen EM: Accelerated dynamic Fourier velocity encoding by exploiting velocity-spatio-temporal correlations. Magnetic Resonance Materials in Physics Biology and Medicine 2004, I7(2):86-94.

70. Baltes C, Hansen M, Tsao J, Kozerke S, Rezavi R, Pedersen E, Boesiger P: Determination of Peak Velocity in Stenotic Areas: Echocardiography versus k-t SENSE Accelerated MR Fourier Velocity Encoding. Radiology 2007:I9.

71. Weber OM, Martin AJ, Higgins CB: Whole-heart steady-state free precession coronary artery magnetic resonance angiography. Magn Reson Med 2003, 50(6): I223-1228.

72. Sorensen TS, Korperich H, Greil GF, Eichhorn J, Barth P, Meyer H, Pedersen EM, Beerbaum P: Operator-independent isotropic three-dimensional magnetic resonance Imaging for morphology in congenital heart disease - A validation study. Circulation 2004, 1 10(2): l63-l69.

73. Plein S, Radjenovic A, Ridgway JP, Barmby D, Greenwood JP, Ball SG, Sivananthan MU: Coronary artery disease: Myocardial perfusion MR imaging with sensitivity encoding versus conventional angiography. Radiology 2005, 235(2):423-430.

74. Wolff SD, Schwitter J, Coulden R, Friedrich MG, Bluemke DA, Bied-erman RW, Martin ET, Lansky AJ, Kashanian F, Foo TKF, Licato PE, Comeau CR: Myocardial first-pass perfusion magnetic resonance imaging - A multicenter dose-ranging study. Circulation 2004, 1 l0(6):732-737.

75. Weber S, Kronfeld A, Kunz RP, Fiebich M, Horstick G, Kreitner KF, Schreiber WG: Comparison of three accelerated pulse sequences for semiquantitative myocardial perfusion Imaging using sensitivity encoding incorporating temporal filtering (TSENSE). Journal of Magnetic Resonance Imaging 2007, 26(3):569-579.

76. Simonetti OP, Kim RJ, Fieno DS, Hillenbrand HB, Wu E, Bundy JM, Finn JP, Judd RM: An improved MR imaging technique for the visualization of myocardial infarction. Radiology 2001, 218(1 ):2l5-223.

77. Dewey M, Laule M, Taupitz M, Kaufels N, Hamm B, Kivelitz D: Myo-cardial viability: Assessment with three-dimensional MR imaging in pigs and patients. Radiology 2006, 239(3):703-709.

78. Kellman P, Larson AC, Hsu LY, Chung YC, Simonetti OP, McVeigh ER, Arai AE: Motion-corrected free-breathing delayed enhancement Imaging of myocardial infarction. Magnetic Resonance in Medicine 2005, 53(1): 194-200.

79. Xu D, King KF, Liang ZP: Improving k-t SENSE by adaptive reg-ularization. Magnetic Resonance in Medicine 2007, 57(5):9l8-930.

80. Candes EJ, Romberg J, Tao T: Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information. Ieee Transactions on Information Theory 2006, 52(2):489-509.

81. Donoho DL: Compressed sensing. Ieee Transactions on Information Theory 2006, 52(4): l289-l306.

82. Pruessmann KP, Weiger M, Boesiger P: Sensitivity encoded cardiac MRI. J Cardiovasc Magn Reson 200l, 3(1): l-9.

83. Kostler H, Sandstede JJW, Lipke C, Landschutz W, Beer M, Hahn D: Auto-SENSE perfusion imaging of the whole human heart. Journal of Magnetic Resonance Imaging 2003, l8(6):702-708.

84. Gharib AM, Herzka DA, Ustun AO, Desai MY, Locklin J, Pettigrew RI, Stuber M: Coronary MR angiography at 3T during diastole and systole. Journal of Magnetic Resonance Imaging 2007, 26(4):92l -926.

85. Jakob PM, Griswold MA, Edelman RR, Manning WJ, Sodickson DK: Accelerated cardiac imaging using the SMASH technique. J Cardiovasc Magn Reson 1999, l(2): l53-l57.

86. Sodickson DK: Tailored SMASH image reconstructions for robust in vivo parallel MR imaging. Magn Reson Med 2000, 44(2):243-25l.

87. McKenzie CA, Yeh EN, Sodickson DK: Improved spatial harmonic selection for SMASH image reconstructions. Magn Reson Med 200l, 46(4):83l-836.

88. Hunold P, Maderwald S, Ladd ME, Jellus V, Barkhausen J: Parallel acquisition techniques in cardiac cine magnetic resonance imaging using TrueFISP sequences: Comparison of image quality and artifacts. Journal of Magnetic Resonance Imaging 2004, 20(3):506-5 ll.

89. Kellman P, Sorger JM, Epstein FH, McVeigh ER: Low-latency temporal filter design for real-time MRI using UNFOLD. Magn Reson Med 2000, 44(6):933-939.

90. Di Bella EV, Wu YJ, Alexander AL, Parker DL, Green D, McGann CJ: Comparison of temporal filtering methods for dynamic contrast MRI myocardial perfusion studies. Magn Reson Med 2003, 49(5):895-902.

91. Guttman MA, Kellman P, Dick AJ, Lederman RJ, McVeigh ER: Realtime accelerated interactive MRI with adaptive TSENSE and UNFOLD. Magn Reson Med 2003, 50(2):3l 5-32l.

92. Ablitt NA, Gatehouse PD, Firmin DN, Yang GZ: Respiratory reordered UNFOLD perfusion imaging. Journal of Magnetic Resonance Imaging 2004, 20(5):8l7-825.

93. Sigfridsson A, Wigstrom L, Kvitting JPE, Knutsson H: k-t(2) BLAST: Exploiting spatiotemporal structure in simultaneously car-

diac and respiratory time-resolved volumetric Imaging. Magnetic Resonance in Medicine 2007, 58(5):922-930.

94. Gamper U, Boesiger P, Kozerke S: Compressed sensing in dynamic MRI. Magnetic Resonance in Medicine 2008, 59(2):365-373.

95. Kellman P, Derbyshire JA, Agyeman KO, McVeigh ER, Arai AE: Extended coverage first-pass perfusion imaging using slice-interleaved TSENSE. Magn Reson Med 2004, 5l(l):200-204.

96. Wintersperger BJ, Bauner K, Reeder SB, Friedrich D, Dietrich O, Sprung KC, Picciolo M, Nikolaou K, Reiser MF, Schoenberg SO: Cardiac steady-state free precession CINE magnetic resonance imaging at 3.0 tesla - Impact of parallel imaging acceleration on volumetric accuracy and signal parameters. Investigative Radiology 2006, 4l (2): l4l-l47.

97. Wintersperger BJ, Sincleair S, Runge VM, Dietrich O, Huber A, Reiser MF, Schoenberg SO: Dual breath-hold magnetic resonance cine evaluation of global and regional cardiac function. European Radiology 2007, l7(l):73-80.

98. Kim D, Kellman P: Improved cine displacement-encoded MRI using balanced steady-state free precession and time-adaptive sensitivity encoding parallel imaging at 3T. Nmr in Biomed-icine 2007, 20(6):59l-60l.

99. Breuer FA, Kellman P, Griswold MA, Jakob PM: Dynamic autocali-brated parallel imaging using temporal GRAPPA (TGRAPPA). Magnetic Resonance in Medicine 2005, 53(4):98l-985.

100. Lai P, Huang F, Larson AC, Li DB: Fast four-dimensional coronary MR angiography with k-t GRAPPA. Journal of Magnetic Resonance Imaging 2008, 27(3):659-665.

Publish with BioMedCentral and every scientist can read your work free of charge

"BioMed Central will be the most significant development for disseminating the results of biomedical research in our lifetime." Sir Paul Nurse, Cancer Research UK

Your research papers will be:

• available free of charge to the entire biomedical community

• peer reviewed and published immediately upon acceptance

• cited in PubMed and archived on PubMed Central

• yours — you keep the copyright

Submit your manuscript here: i } BioMedcentral ^