Skip to main content

Acceleration of tissue phase mapping by k-t BLAST: a detailed analysis of the influence of k-t-BLAST for the quantification of myocardial motion at 3T



The assessment of myocardial motion with tissue phase mapping (TPM) provides high spatiotemporal resolution and quantitative motion information in three directions. Today, whole volume coverage of the heart by TPM encoding at high spatial and temporal resolution is limited by long data acquisition times. Therefore, a significant increase in imaging speed without deterioration of the quantitative motion information is required. For this purpose, the k-t BLAST acceleration technique was combined with TPM black-blood functional imaging of the heart. Different k-t factors were evaluated with respect to their impact on the quantitative assessment of cardiac motion.


It is demonstrated that a k-t BLAST factor of two can be used with a marginal, but statistically significant deterioration of the quantitative motion data. Further increasing the k-t acceleration causes substantial alteration of the peak velocities and the motion pattern, but the temporal behavior of the contraction is well maintained up to an acceleration factor of six.


The application of k-t BLAST for the acceleration of TPM appears feasible. A reduction of the acquisition time of almost 45% could be achieved without substantial loss of quantitative motion information.


Quantification of myocardial mechanics is supposed to provide an improved understanding of cardiac motion as well as to enable a more detailed assessment of certain myocardial diseases such as cardiac insufficiency. A major limitation in quantification of cardiac function is the long measurement time required for three-dimensional (3D) velocity encoded imaging. However, in diagnosis and staging of certain cardiac diseases and for therapy selection, 3D functional information of the myocardial motion appears mandatory. Especially for the selection of patients eligible for cardiac resynchronization therapy (CRT), quantification of the 3D-cardiac motion appears paramount to reduce or completely avoid non-responders, which represent 30% of treated patients using current selection criteria [1].

Four main approaches have been introduced for the assessment of myocardial mechanics including tagging [24], displacement encoding with stimulated echoes (DENSE) [58], strain encoding (SENC) [9] and tissue phase mapping (TPM) [1014], which has also been introduced as phase contrast velocity encoded imaging [15, 16] of tissue.

In the tagging technique, lines or a grid are mapped onto the myocardium by either spatial modulation techniques [2, 3] or a DANTE pulse train in the presence of a frequency-encoding gradient [17]. Direct analysis of the tag-deformation over the cardiac cycle provides access to the inter-voxel strain and velocity of the myocardium, but is limited by the spatial resolution of the tag pattern. This can partly be solved by applying dedicated post-processing techniques such as the harmonic phase approach (HARP) [18].

The DENSE approach directly encodes displacements over long time intervals at high spatial density [5]. Due to the long displacement encoding intervals, data acquisition is very slow.

In the SENC technique, an intra-voxel tag-pattern is used for the assessment of the intra-voxel strain, which enables the assessment of the stiffness of the myocardium. The application of the SENC technique as the sole technique for the assessment of the cardiac function is limited by the lack of information on the inter-voxel strain and myocardial velocities.

In TPM, the myocardial velocity is directly encoded by the application of bipolar gradients causing the spins to acquire a phase that is directly proportional to their velocity. Since the direction of the velocity encoding gradients can be chosen freely, TPM enables the quantitative assessment of the 3D flow vector. Wide application of TPM is still limited by the long acquisition times, which preclude large volume coverage at sufficient spatial resolution and may introduce image deterioration due to varying respiratory or irregular cardiac motion [19].

For acceleration of the image acquisition, several methods have been introduced. Local imaging techniques aim at reducing the field-of-view (FOV) to a confined area containing the heart [1921]. Its sensitivity to patient motion and the required complicated planning of the anatomy have limited their clinical application. More promising techniques employ correlations in k-space or image space like sensitivity encoding (SENSE) [22], generalized autocalibrating partially parallel acquisitions (GRAPPA) [23] and partial Fourier methods [24].

View sharing exploits temporal correlations by reusing some of the same k-space data in order to reconstruct additional images [2528]. With view sharing, a decrease of the acquisition time of 37.5% could be obtained without significant deterioration of the velocity mapping data [28]. Temporal correlations are also exploited in the UNFOLD approach (unaliasing by Fourier-encoding the overlaps using temporal dimension) [29, 30], which avoids aliasing resulting from undersampling by shifting the sampling function in time, such that Fourier transformation through time can resolve these overlaps.

More recently dedicated techniques like k-t BLAST and k-t SENSE exploiting both correlations in k-space and in time by sparse sampling have been introduced [3134]. The resulting aliased images in the reciprocal spatio-temporal frequency domain are resolved using the information of prior acquired low resolution data leading to aliasing free images. It seems feasible to use k-t BLAST for accelerated image data acquisition in applications involving quasiperiodic motion such as the heart.

K-t BLAST has been applied to various applications in the medical field [3539]. Especially it has proven its applicability to velocity encoded imaging for flow quantification [32, 4042].

The objective of this contribution is to investigate the potential of k-t BLAST for the acceleration of TPM image acquisition at 3 T. The k-t BLAST technique was evaluated applying different acceleration factors for the assessment of quantitative myocardial motion and compared to the non accelerated technique.


Data acquisition

20 adult volunteers (7 females, 13 males, age 29 ± 11 years) were enrolled in the study [kt-group]. All volunteers enrolled underwent one non-k-t BLAST TPM data acquisition and 6 k-t BLAST accelerated sequences with k-t BLAST factors R ranging from 2 to 7. The sequence order was randomized to reduce the influence of physiological variations on the myocardial motion between the acquired sequences. To assess of the reproducibility of the approach, the non k-t BLAST protocol (no-kt) was repeated twice in 20 additional volunteers (fourteen males, six females, age 28 ± 6 years) [reference group]. The study protocols were approved by the local ethics committee and informed written consent was obtained from all volunteers prior to the MRI examination.

All MRI scans were performed on a 3 T whole body MR scanner (Achieva 3.0 T, Philips, Best, The Netherlands) with a 32 [2 × 4 × 4] channel phased array cardiac coil. A vector ECG was applied for cardiac triggering. Breath-hold cine cardiac two-chamber and four-chamber views were acquired to define the short axis image geometry, which was used in all subsequent acquisitions.

The TPM acquisition was performed applying a respiratory navigated segmented and velocity encoded cardiac triggered gradient echo sequence. The acquisition parameters were as follows: TE/TR = 4.7 ms/7.1 ms, flip angle α = 15°, FOV = 340 × 340 mm2, acquisition matrix (M*P) = 172*168, slice thickness = 8 mm, in plane resolution: 2 × 2 mm2, 3 k-lines per segment and one startup echo, VENC = 30 cm/s in all 3 encoding directions. The acquisition window was 90% of the RR interval. The duration of the saturation module consisting of saturation pulses and spoiler gradients was 12 ms. The phase interval without k-t BLAST acceleration was 40.4 ms. The navigator was 15.5 ms long, the navigator feedback time 5 ms. To solely investigate the impact of the k-t BLAST, data acquisition was not combined with further acceleration techniques like view sharing or parallel imaging. The TPM encoding was performed in all 3 spatial orientations in consecutive heart beats. A conventional pencil beam navigator through the dome of the right hemi-diaphragm was applied at each start of the cardiac cycle for respiratory gating [43, 44].

To avoid flow related artifacts in the phase images caused by the strong blood flow in the ventricle and to improve the delineation of myocardium and blood, two saturation pulses were incorporated superior and inferior of the imaged slice to generate black-blood contrast [13]. To avoid idle times due to high SAR demands of the sequence at 3 T, the saturation pulses were applied alternating and the maximal B1-amplitude of the RF pulses was optimized to 8 μT [45].

For the non-accelerated technique and a heart rate of 60 beats per minute, the maximum number of heart phases was 21 and the nominal scan duration was 225 s. In this study, the number of acquired cardiac phases was a multiple integer of the k-t BLAST factor to ensure similar undersampling factors in all regions of k-space. The respective nominal scan times and maximum number of acquired phases for different acceleration factors are listed in table 1 exemplary for a heart rate of 60 beats per minute. For all k-t BLAST measurements, the number of ky phase encoding steps for the training scan was set to eleven. Please note, that the acceleration factor only impacts the acquisition stage, whereas the length of the training stage remains constant. Hence the relative gain in acquisition speed decreases with increasing k-t factors.

Table 1 Relationship between k-t BLAST factor, maximum number of heart phases and scan time.

Data analysis

The TPM MR images were processed by in house developed MATLAB programs (MATLAB R2008a; Mathworks, Natick, Mass). The segmentation of the myocardium was performed automatically, relying on active contour techniques by incorporating a shape model. After the segmentation of the first phase, the information was propagated through the entire sequence by tracking profile intensities [46, 47]. Background phase error correction was performed using a linear fit to the phase of static tissue as suggested earlier [48]. The radial (towards the center of the blood pool) and longitudinal (towards the apex of the heart) velocity curves were calculated from the acquired three-directional velocity vector. Prior to the analysis the velocity data were interpolated over time by cubic splines to provide a continuous velocity profile thus enabling the comparison of sequences with a different number of sampled heart phases. Physiologically, the accumulated phase over the entire heart cycle must result to zero. To compensate for non-linear phase error contributions, in a subsequent correction step the resulting velocity curves were shifted accordingly to meet the physiological conditions.

To assess the impact of the k-t factor on the quantitative velocity information the following parameters were derived:

  • The systolic and diastolic peak velocities vp,sys and vp,dias were determined and the resulting velocity difference Δv = vp,sys - vp,dias was calculated for each sequence. Bland-Altman analysis was performed for the velocity differences Δv and the differences between Δv with and without kt-BLAST is denoted as Δv Diff. The peak factor PF was defined as the ratio of the velocity differences with and without k-t acceleration (kt-group) and accordingly as the ratio between the two reproducibility measurements (reference group): PF (seq. 1, seq.2) = Δvseq.2/Δvseq.1. This parameter was determined to evaluate whether k-t acceleration has an impact on clinical main features of the motion curve. Ideally, PF should be one. Due to the temporal smoothing of the k-t BLAST algorithm, it is expected that especially sharp peaks will be abraded. The PF for the radial and longitudinal velocity curves were referred as PFr and PFl.

  • The normalized root mean square deviation nRMSD between the velocity curves with and without k-t acceleration (kt-group) and accordingly between the curves obtained by the reproducibility measurements (reference group) was calculated. The normalization was performed by dividing the root mean square deviation by Δvnokt. The radial and longitudinal nRMSDs were denoted as nRMSDr and nRMSDl.

  • The correlation coefficients were determined for both groups to evaluate the statistical dependency between the velocity curves. Let u be the velocity of the first sequence to compare, and w be the velocity of the second sequence. Than ui and wi are the spline interpolated data at different time steps (step size: 0.01 ms), n is the number of time steps, ū and w ¯ are the mean values of ui and wi and σu and σw the corresponding standard deviation. The correlation coefficient c is than calculated as: c = Σ i = 1 n [ ( u i u ¯ σ u ) ( w i w ¯ σ w ) ] . Ideally the correlation coefficient should be one. The correlation coefficient between the radial velocity curves was denoted as cr, the correlation coefficient between the longitudinal velocity curves as cl.

  • The times to the peak diastolic velocity of the radial and longitudinal velocity tr,dias and tl,dias was identified for each acquisition technique in order to determine, whether the temporal behavior of the myocardium is preserved for different acceleration factors. The temporal behavior is of special interest in various cardiac diseases such as cardiac asynchrony. Bland-Altman analysis was performed for tr,dias and tl,dias and the mean difference times Δtr,dias (seq.1, seq.2) = tr,dias,seq.2- tr,dias,seq.1 and Δtl,dias (seq.1, seq.2) = tl,dias,seq.2- tl,dias,seq.1 and their standard deviations over all volunteers were compared.

For the evaluation of significances an unpaired two-tailed student's t-test was performed. Values below 0.05 were considered to be significant. The variance of both reference and kt-group was considered to be the same.


Figure 1 exemplarily shows anatomical and velocity encoded images of the myocardium acquired without and with k-t BLAST acceleration applying acceleration factors of two, four and six. All data were acquired approximately 50 ms after the R-Wave. For low k-t BLAST acceleration factors, the anatomical and phase encoded images do not show obvious deterioration. With increasing k-t BLAST factors, however, an increasing blurring can be observed in the anatomical images. In the TPM data, a reduction of velocity amplitude with increasing k-t BLAST factors can be noticed.

Figure 1
figure 1

Anatomical and phase contrast images with and without k-t BLAST. Anatomical (left) images of the myocardium and respective velocity encoded phase images in AP, FH and RL directions for the conventional sequence and the k-t BLAST accelerated sequences (acceleration factors 2, 4 and 6, heart phase with maximum longitudinal velocity (approx. 50 ms after the R-Wave)

Figure 2 shows the radial (a, c) and the longitudinal (b, d) velocity curves over time exemplarily for one volunteer out of the kt-group and one volunteer out of the reference group. The data clearly reveal a decrease of the measured peak velocities with increasing k-t BLAST acceleration factors, which is substantially higher than the reproducibility of the TPM-technique.

Figure 2
figure 2

Velocity curves for different k-t accelerating factors and for the reproducibility study. In Figure 2 a) and b) the radial and longitudinal velocity curves of a healthy volunteer over the time after the R-peak are displayed for different k-t accelerations (kt-group). In c) and d) two velocity curves without k-t acceleration are displayed (reference group).

The PF analysis provides more quantitative information of the amplitudinal decrease of the velocity curves. Table 2 provides the peak factors PFr and PFl for the different k-t factors for the radial and longitudinal velocity of the myocardium. PFr and PFl are both very sensitive to k-t acceleration. A deviation of less than 10% from the expected value of one can only be obtained for a k-t factor of two. Even for this k-t factor, the standard deviation of PFr is doubled compared to the reproducibility data, whereas the standard deviation of PFl is in the same range. For all investigated k-t BLAST factors, there is a significant (p-value ≤ 0.05) underestimation of the peak-factor.

Table 2 Mean radial and longitudinal peak-factors PFr and PFl and their standard deviations

The results of the Bland-Altman analysis of the velocity ranges Δvr and Δvl for the kt-group and reference group are listed in table 3. Significant increased mean differences ΔvrDiff and ΔvlDiff were obtained for increasing acceleration factors.

Table 3 Bland-Altman results for the velocity range differences ΔvrDiff and ΔvlDiff and their standard deviations

The mean nRMSDr and nRMSDl values and their standard deviations are listed in table 4. The values reveal an increasing nRMSD value with increasing k-t acceleration. The small nRMSD values for the repeated no-kt measurements (nRMSDr: 0.04 ± 0.02; nRMSDl: 0.04 ± 0.02) indicate good reproducibility for the TPM technique. All nRMSDr and nRMSDl values are significantly higher (p-value ≤ 0.05) than the reproducibility of the technique. However, for a k-t BLAST factor of two and four the nRMSD values can be kept below 10%.

Table 4 Normalized mean root mean square deviations nRMSDr and nRMSDl and their standard deviations

The mean radial and longitudinal correlation coefficients cr and cl and their standard deviations are listed in table 5. A significant decrease (p-value < 0.05) of cr and cl was obtained with increasing k-t BLAST factors. All correlation coefficients are above 0.8.

Table 5 Mean radial and longitudinal correlation coefficients cr and cl and their standard deviations

The results of the Bland-Altman analysis for the time to the diastolic peak velocity are shown in table 6. For all acceleration factors R Δtr,dias and Δtl,dias are below 12 ms. For Δtr,dias no significant differences could be obtained between the reproducibility measurement and the k-t accelerated measurements. Only for Δtl,dias significant differences were obtained for R equal to four, five and seven.

Table 6 Results of the Bland-Altman analysis for Δtr,dias and Δtl,dias and their standard deviations


The application of k-t BLAST to accelerate black-blood TPM imaging for the quantification of the myocardial velocity appears feasible. With increasing k-t accelerating factors R an increasing deterioration of the velocity curve and a decrease of the peak velocity is observed. This effect can likely be explained by the inherent slight temporal smoothing of the k-t BLAST algorithm [32, 40, 42]. With R = 2, the impact on the velocity can be kept very small (≤ 10% deviation from the no-kt measurement for the peak-factors, normalized RMSDs and correlation coefficient and less than 5 ms temporal deviation between the time to the diastolic peak). Nevertheless, the standard deviation of the radial peak-factor was doubled, which might have implications for clinical studies, where the overlap between velocity ranges of healthy volunteers and patients might be increased. Higher acceleration factors show substantial degradation of the peak factors and motion pattern. There is only minimal influence of k-t BLAST on the time to the minimum, indicating the applicability of k-t BLAST for the assessment of temporal behavior.

A similar reduction of the peak velocities was observed in prior work published on the application of k-t BLAST to quantitative flow measurements [32, 40, 42]. However, in this study a substantial reduction was also observed at lower acceleration factors, which might be attributed to the more complex motion pattern in the case of the myocardium.

In this study, the number of acquired cardiac phases was a multiple integer of the k-t BLAST factor to ensure same undersampling factors in all regions of k-space. Depending on the heart rate and the specific acceleration factor, this might introduce worse temporal sampling, which can also cause lower peak factors [11]. Further improvement may be obtained by an application specific adaptation of the number of cardiac phases. Since only a k-t BLAST factor of two appears reasonable for quantification in tissue phase mapping, choosing multiple integers of the acceleration factor results in the loss of one heart phase at maximum.

Since the possible acceleration by k-t BLAST appears limited to 2, its combination with other acceleration techniques must be considered for further increasing imaging speed. Most promising here may be the combination with parallel imaging techniques (k-t SENSE, t-Grappa [31, 49]) or the combination with view sharing.


In summary, a k-t BLAST factor of two can be applied with statistically significant but not substantial loss of motion information of the myocardium, enabling a 45% decrease in scan duration or a 1.8 fold increase in volume coverage. Higher accelerating factors show substantial degradation of the motion pattern and should therefore be avoided.

The k-t BLAST sequence has the potential to enable 3D tissue phase mapped myocardial imaging with reasonable image acquisition times. The possible combination with parallel imaging techniques offers a way to further reduce the overall scan time.



Cardiac Resynchronization Therapy


Delay alternating with nutation for tailored excitation


Displacement encoding with stimulated echoes




Generalized autocalibrating partially parallel acquisitions


Harmonic phase approach


Normalized root mean square deviation


Peak factor


k-t Blast acceleration factor


Specific absorption rate




Sensitivity encoding


Tissue Phase Mapping


Unaliasing by Fourier-encoding the overlaps using temporal dimension


  1. Ypenburg C, Westenberg JJ, Bleeker GB, VAN de Veire N, Marsan NA, Henneman MM, van der Wall EE, Schalij MJ, Abraham TP, Barold SS, Bax JJ: Noninvasive imaging in cardiac resynchronization therapy--part 1: selection of patients. Pacing Clin Electrophysiol. 2008, 31: 1475-99. 10.1111/j.1540-8159.2008.01212.x.

    Article  PubMed  Google Scholar 

  2. Fischer SE, McKinnon GC, Scheidegger MB, Prins W, Meier D, Boesiger P: True myocardial motion tracking. Magn Reson Med. 1994, 31: 401-13. 10.1002/mrm.1910310409.

    Article  CAS  PubMed  Google Scholar 

  3. Fischer SE, McKinnon GC, Maier SE, Boesiger P: Improved myocardial tagging contrast. Magn Reson Med. 1993, 30: 191-200. 10.1002/mrm.1910300207.

    Article  CAS  PubMed  Google Scholar 

  4. Axel L, Dougherty L: Heart wall motion: improved method of spatial modulation of magnetization for MR imaging. Radiology. 1989, 172: 349-50.

    Article  CAS  PubMed  Google Scholar 

  5. Aletras AH, Ding S, Balaban RS, Wen H: DENSE: Displacement Encoding with Stimulated Echoes in Cardiac Functional MRI. J Magn Reson. 1999, 137: 247-252. 10.1006/jmre.1998.1676.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Spottiswoode BS, Zhong X, Lorenz CH, Mayosi BM, Meintjes EM, Epstein FH: 3D myocardial tissue tracking with slice followed cine DENSE MRI. J Magn Reson Imaging. 2008, 27: 1019-27. 10.1002/jmri.21317.

    Article  PubMed  Google Scholar 

  7. Aletras AH, Balaban RS, Wen H: High-resolution strain analysis of the human heart with fast-DENSE. J Magn Reson. 1999, 140: 41-57. 10.1006/jmre.1999.1821.

    Article  CAS  PubMed  Google Scholar 

  8. Aletras AH, Wen H: Mixed echo train acquisition displacement encoding with stimulated echoes: an optimized DENSE method for in vivo functional imaging of the human heart. Magn Reson Med. 2001, 46: 523-34. 10.1002/mrm.1223.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Osman NF: Detecting stiff masses using strain-encoded (SENC) imaging. Magn Reson Med. 2003, 49: 605-8. 10.1002/mrm.10376.

    Article  PubMed  Google Scholar 

  10. Jung B, Markl M, Föll D, Hennig J: Investigating myocardial motion by MRI using tissue phase mapping. Eur J Cardiothorac Surg. 2006, 29 (Suppl 1): S150-7. 10.1016/j.ejcts.2006.02.066.

    Article  PubMed  Google Scholar 

  11. Jung B, Zaitsev M, Hennig J, Markl M: Navigator gated high temporal resolution tissue phase mapping of myocardial motion. Magn Reson Med. 2006, 55: 937-42. 10.1002/mrm.20808.

    Article  PubMed  Google Scholar 

  12. Jung B, Föll D, Böttler P, Petersen S, Hennig J, Markl M: Detailed analysis of myocardial motion in volunteers and patients using high-temporal-resolution MR tissue phase mapping. J Magn Reson Imaging. 2006, 24: 1033-9. 10.1002/jmri.20703.

    Article  PubMed  Google Scholar 

  13. Delfino JG, Johnson KR, Eisner RL, Eder S, Leon AR, Oshinski JN: Three-directional myocardial phase-contrast tissue velocity MR imaging with navigator-echo gating: in vivo and in vitro study. Radiology. 2008, 246: 917-25. 10.1148/radiol.2463062155.

    Article  PubMed  Google Scholar 

  14. Petersen SE, Jung BA, Wiesmann F, Selvanayagam JB, Francis JM, Hennig J, Neubauer S, Robson MD: Myocardial tissue phase mapping with cine phase-contrast mr imaging: regional wall motion analysis in healthy volunteers. Radiology. 2006, 238: 816-26. 10.1148/radiol.2383041992.

    Article  PubMed  Google Scholar 

  15. Pelc NJ, Herfkens RJ, Shimakawa A, Enzmann DR: Phase contrast cine magnetic resonance imaging. Magn Reson Q. 1991, 7: 229-54.

    CAS  PubMed  Google Scholar 

  16. Thompson RB, McVeigh ER: High temporal resolution phase contrast MRI with multiecho acquisitions. Magn Reson Med. 2002, 47: 499-512. 10.1002/mrm.10079.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Mosher TJ, Smith MB: A DANTE tagging sequence for the evaluation of translational sample motion. Magn Reson Med. 1990, 15: 334-9. 10.1002/mrm.1910150215.

    Article  CAS  PubMed  Google Scholar 

  18. Osman NF, Kerwin WS, McVeigh ER, Prince JL: Cardiac motion tracking using CINE harmonic phase (HARP) magnetic resonance imaging. Magn Reson Med. 1999, 42: 1048-60. 10.1002/(SICI)1522-2594(199912)42:6<1048::AID-MRM9>3.0.CO;2-M.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Kozerke S, Plein S: Accelerated CMR using zonal, parallel and prior knowledge driven imaging methods. J Cardiovasc Magn Reson. 2008, 10: 29-10.1186/1532-429X-10-29.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Feinberg DA, Hoenninger JC, Crooks LE, Kaufman L, Watts JC, Arakawa M: Inner volume MR imaging: technical concepts and their application. Radiology. 1985, 156: 743-7.

    Article  CAS  PubMed  Google Scholar 

  21. Weiß S, Rasche V: Projection-reconstruction reduced FOV imaging. Magn Reson Imag. 1999, 17: 517-525.

    Article  Google Scholar 

  22. Pruessmann KP, Weiger M, Scheidegger MB, Boesiger P: SENSE: sensitivity encoding for fast MRI. Magn Reson Med. 1999, 42: 952-62. 10.1002/(SICI)1522-2594(199911)42:5<952::AID-MRM16>3.0.CO;2-S.

    Article  CAS  PubMed  Google Scholar 

  23. 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: 1202-1210. 10.1002/mrm.10171.

    Article  PubMed  Google Scholar 

  24. McGibney G, Smith MR, Nichols ST, Crawley A: Quantitative evaluation of several partial Fourier reconstruction algorithms used in MRI. Magn Reson Med. 1993, 30: 51-9. 10.1002/mrm.1910300109.

    Article  CAS  PubMed  Google Scholar 

  25. Foo TK, Bernstein MA, Aisen AM, Hernandez RJ, Collick BD, Bernstein T: Improved ejection fraction and flow velocity estimates with use of view sharing and uniform repetition time excitation with fast cardiac techniques. Radiology. 1995, 195: 471-8.

    Article  CAS  PubMed  Google Scholar 

  26. Doyle M, Walsh EG, Blackwell GG, Pohost GM: Block regional interpolation scheme for k-space (BRISK): a rapid cardiac imaging technique. Magn Reson Med. 1995, 33: 163-70. 10.1002/mrm.1910330204.

    Article  CAS  PubMed  Google Scholar 

  27. Parrish T, Hu X: Continuous update with random encoding (CURE): a new strategy for dynamic imaging. Magn Reson Med. 1995, 33: 326-36. 10.1002/mrm.1910330307.

    Article  CAS  PubMed  Google Scholar 

  28. Markl M, Hennig J: Phase contrast MRI with improved temporal resolution by view-sharing: k-space related velocity mapping properties. Magn Reson Imag. 2001, 19: 669-676. 10.1016/S0730-725X(01)00386-1.

    Article  CAS  Google Scholar 

  29. 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: 813-28. 10.1002/(SICI)1522-2594(199911)42:5<813::AID-MRM1>3.0.CO;2-S.

    Article  CAS  PubMed  Google Scholar 

  30. Tsao J: On the UNFOLD method. Magn Reson Med. 2002, 47: 202-7. 10.1002/mrm.10024.

    Article  PubMed  Google Scholar 

  31. 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: 1031-42. 10.1002/mrm.10611.

    Article  PubMed  Google Scholar 

  32. 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. Magn Reson Med. 2005, 54: 1430-8. 10.1002/mrm.20730.

    Article  PubMed  Google Scholar 

  33. Maredia N, Kozerke S, Larghat A, Abidin N, Greenwood JP, Boesiger P, Plein S: Measurement of left ventricular dimensions with contrast-enhanced three-dimensional cine imaging facilitated by k-t SENSE. J Cardiovasc Magn Reson. 2008, 10: 27-10.1186/1532-429X-10-27.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Plein S, Ryf S, Schwitter J, Radjenovic A, Boesiger P, Kozerke S: Dynamic contrast-enhanced myocardial perfusion MRI accelerated with k-t sense. Magn Reson Med. 2007, 58: 777-85. 10.1002/mrm.21381.

    Article  PubMed  Google Scholar 

  35. Kozerke S, Tsao J, Razavi R, Boesiger P: Accelerating cardiac cine 3D imaging using k-t BLAST. Magn Reson Med. 2004, 52: 19-26. 10.1002/mrm.20145.

    Article  PubMed  Google Scholar 

  36. Hamdan A, Kelle S, Schnackenburg B, Wellnhofer E, Fleck E, Nagel E: Single-breathhold four-dimensional assessment of left ventricular volumes and function using k-t BLAST after application of extracellular contrast agent at 3 Tesla. J Magn Reson Imaging. 2008, 27: 1028-36. 10.1002/jmri.21332.

    Article  PubMed  Google Scholar 

  37. 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. J Magn Reson Imaging. 2007, 25: 737-42. 10.1002/jmri.20877.

    Article  PubMed  Google Scholar 

  38. Gebker R, Jahnke C, Paetsch I, Schnackenburg B, Kozerke S, Bornstedt A, Fleck E, Nagel E: MR myocardial perfusion imaging with k-space and time broad-use linear acquisition speed-up technique: feasibility study. Radiology. 2007, 245: 863-71. 10.1148/radiol.2453061701.

    Article  PubMed  Google Scholar 

  39. 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, 241: 718-28. 10.1148/radiol.2413051522.

    Article  PubMed  Google Scholar 

  40. 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. MAGMA. 2004, 17: 86-94. 10.1007/s10334-004-0062-8.

    Article  PubMed  Google Scholar 

  41. Marshall I: Feasibility of k-t BLAST technique for measuring "seven-dimensional" fluid flow. J Magn Reson Imaging. 2006, 23: 189-96. 10.1002/jmri.20494.

    Article  PubMed  Google Scholar 

  42. Stadlbauer A, van der Riet W, Globits S, Crelier G, Salomonowitz E: Accelerated phase-contrast MR imaging: comparison of k-t BLAST with SENSE and Doppler ultrasound for velocity and flow measurements in the aorta. J Magn Reson Imaging. 2009, 29: 817-24. 10.1002/jmri.21706.

    Article  PubMed  Google Scholar 

  43. Wang Y, Rossman PJ, Grimm RC, Riederer SJ, Ehman RL: Navigator-echo-based real-time respiratory gating and triggering for reduction of respiration effects in three-dimensional coronary MR angiography. Radiology. 1996, 198: 55-60.

    Article  CAS  PubMed  Google Scholar 

  44. Bellenger NG, Gatehouse PD, Rajappan K, Keegan J, Firmin DN, Pennell DJ: Left ventricular quantification in heart failure by cardiovascular MR using prospective respiratory navigator gating: comparison with breath-hold acquisition. J Magn Reson Imaging. 2000, 11: 411-417. 10.1002/(SICI)1522-2586(200004)11:4<411::AID-JMRI9>3.0.CO;2-B.

    Article  CAS  PubMed  Google Scholar 

  45. Lutz A, Bornstedt A, Hombach V, Manzke R, Rasche V: Low SAR black-blood CINE imaging at 3T. Proc. ESMRM B: 1-3 October 2009; Antalya. 2009, Antalya, 547-

    Google Scholar 

  46. Fradkin M, Ciofolo C, Mory B, Hautvast G, Breeuwer M: Comprehensive segmentation of cine cardiac MR images. Med Image Comput Comput Assist Interv. 2008, 11: 178-185.

    PubMed  Google Scholar 

  47. Hautvast G, Lobregt S, Breeuwer M, Gerritsen F: Automatic contour propagation in cine cardiac magnetic resonance images. IEEE Trans Med Imaging. 2006, 25: 1472-1482. 10.1109/TMI.2006.882124.

    Article  PubMed  Google Scholar 

  48. Walker PG, Cranney Gregory B, Scheidegger Markus B, Waseleski Gena, Pohost Gerald M, Yoganathan Ajit P: Semiautomated method for noise reduction and background phase error correction in MR phase velocity data. J Magn Reson Imaging. 1993, 3: 521-530. 10.1002/jmri.1880030315.

    Article  CAS  PubMed  Google Scholar 

  49. Huang F, Akao J, Vijayakumar S, Duensing GR, Limkeman M: k-t GRAPPA: a k-space implementation for dynamic MRI with high reduction factor. Magn Reson Med. 2005, 54: 1172-84. 10.1002/mrm.20641.

    Article  PubMed  Google Scholar 

Download references

Author information

Authors and Affiliations


Corresponding author

Correspondence to Anja Lutz.

Additional information

Competing interests

VR and AL have a research grant with Philips Healthcare. PE is employed by Philips Healthcare. RM is employed by Philips Research.

Authors' contributions

AL developed the sequence protocol, performed the data acquisition, the analysis and interpretation of data and drafted the manuscript. AB was involved in developing the sequence and in interpreting the data. PE was involved in the analysis of the data. RM and GUN were involved in the interpretation of the data. VR and AB made substantial contributions to conception and design and revised the manuscript critically for important intellectual content. All authors have given final approval of the version to be published.

Authors’ original submitted files for images

Rights and permissions

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Lutz, A., Bornstedt, A., Manzke, R. et al. Acceleration of tissue phase mapping by k-t BLAST: a detailed analysis of the influence of k-t-BLAST for the quantification of myocardial motion at 3T. J Cardiovasc Magn Reson 13, 5 (2011).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: