Open Access

Acceleration of tissue phase mapping with sensitivity encoding at 3T

  • Anja Lutz1Email author,
  • Axel Bornstedt1,
  • Robert Manzke2,
  • Patrick Etyngier3,
  • G Ulrich Nienhaus4,
  • Wolfgang Rottbauer1 and
  • Volker Rasche1
Journal of Cardiovascular Magnetic Resonance201113:59

DOI: 10.1186/1532-429X-13-59

Received: 15 March 2011

Accepted: 12 October 2011

Published: 12 October 2011

Abstract

Background

The objective of this study was to investigate the impact of sensitivity encoding on the quantitative assessment of cardiac motion in black blood cine tissue phase mapping (TPM) sequences. Up to now whole volume coverage of the heart is still limited by the long acquisition times. Therefore, a significant increase in imaging speed without deterioration of quantitative motion information is indispensable.

Methods

20 volunteers were enrolled in this study. Each volunteer underwent myocardial short-axis TPM scans with different SENSE acceleration factors. The influence of SENSE acceleration on the measured motion curves was investigated.

Results

It is demonstrated that all TPM sequences with SENSE acceleration have only minimum influence on the motion curves. Even with a SENSE factor of four, the decrease in the amplitude of the motion curve was less than 3%. No significant difference was observed for the global correlation coefficient and deviation between the motion curves obtained by the reproducibility and the SENSE accelerated measurements.

Conclusions

It is feasible to accelerate myocardial TPM measurements with SENSE factors up to 4 without losing substantial information of the motion pattern.

Background

Quantification of myocardial motion provides insight in the myocardial mechanics and enables a more detailed assessment of motion abnormalities in several cardiac diseases such as cardiac insufficiency. Among other imaging modalities like echocardiographic tissue Doppler imaging (TDI) [1], cardiovascular Magnetic Resonance (CMR) has proven its value for noninvasive assessment of global and regional cardiac function [25]. Myocardial motion can be quantified with CMR using a variety of techniques including tagging [6, 7], strain encoding (SENC) [8], displacement encoding with stimulated echoes (DENSE) [9], and phase contrast velocity encoding (tissue phase mapping, TPM) [2, 1014].

TPM directly encodes the myocardial velocity by the application of bipolar gradients and enables the quantitative assessment of three-directional motion. The spatial resolution of TPM is only limited by the voxel size. Despite the obvious advantages of TPM, its wide application is still limited by the long acquisition times, which precludes large volume coverage at sufficient spatial resolution and may cause image degradation due to irregular respiratory and cardiac motion [15]. For wide clinical applications an acceleration of acquisition time without information loss is mandatory.

Several generic methods accelerating the image acquisition have been introduced. Local imaging techniques reduce the field-of-view (FOV) to a confined area in order to reduce the scan time [1517]. Clinical applications of these techniques are still limited due to their sensitivity to patient motion and the required complicated planning of the anatomy. Accelerating techniques like view sharing [18, 19] and k-t BLAST [20] exploit temporal correlations. View sharing has shown to measure an accurate velocity profile with a gain of acquisition speed in the order of 37.5% [19]. Accelerating techniques like k-t BLAST allow higher accelerating factors. A disadvantage of the k-t BLAST algorithm is a reduction of peak velocity values due to temporal smoothing, which was observed in quantitative phase contrast angiography (PCA) [2123] as well in myocardial tissue phase mapping [24].

Parallel imaging techniques exploiting coil sensitivities for unaliasing like sensitivity encoding (SENSE) [25] and generalized auto-calibrating partially parallel acquisitions (GRAPPA) [26] enable high acceleration factors. No temporal smoothing is expected. The aim of this work was to quantify the feasibility of using SENSE for tissue phase mapping of the left ventricle. An analysis of the influence of SENSE on TPM based flow quantification was performed and investigated for different accelerating factors. A case example is added to show the possibility of the combination of SENSE with TPM to measure the left-ventricular myocardial motion in 3D.

Methods

Volunteers

A total of 20 volunteers (8 females, 12 males, age 30 ± 10 years) were enrolled in this study. One additional volunteer was examnied to show the possibility of the use of SENSE to enable 3D TPM imaging (male,24 years). The study protocol was approved by the local ethics committee. All volunteers provided written informed consent prior to the CMR examination.

Data acquisition

Image acquisition was performed on a 3T whole body MR scanner (Achieva 3.0T, Philips, Best, The Netherlands) with a 32 [2 × 4 × 4] channel phased array cardiac coil.

A coil-sensitivity reference scan was acquired for obtaining the coil sensitivity profiles as prerequisite for the subsequent SENSE reconstructions. Breath-hold cine cardiac two- and four chamber views were acquired to define the equatorial short-axis image orientation.

The TPM acquisition was performed for one equatorial short axis applying a black blood prepared [14, 27], respiratory navigated, segmented and velocity encoded cardiac triggered gradient echo sequence. The acquisition parameters are listed in table 1.
Table 1

Acquisition protocol for the 2D and 3D data acquisitions

Acquisition parameter

Value

 

sequence type

velocity encoded segmented gradient echo

velocity encoded segmented gradient echo

TR/TE

7 ms/4.7 ms

7.1 ms/4.9 ms

FOV

340 × 340 mm2

380 × 380 mm2

flip angle

15°

15°

k-lines per segment

3 + 1 startup echo

3

acquisition matrix

172 × 168

128 × 124

spatial resolution

2 × 2 × 8 mm3

3 × 3 × 3 mm3

  

(with over-contiguous slices)

black blood mode

alternating

alternating

number of slices

1

21

cardiac phases at 60 bpm

21

23

phase interval

40 ms

37.3 ms

VENC

30 cm/s

20 cm/s

navigator duration + evaluation

15.5 ms + 5 ms

15.5 ms + 5 ms

number of chunks

 

3 (to enable black blood suppression)

RR interval

90%

90%

SENSE acceleration factor R

1,2,3,4

4

nominal scan duration [min]

3:45, 1:57, 1:17, 0:57 (corresponding to R)

15:30

Isotropic velocity encoding of 30 cm/s was performed in a Hadamard fashion by a four-point velocity vector method [28, 29]. To improve the temporal resolution, the different flow encoding directions were encoded in subsequent heart beats [14, 24].

Black blood preparation was performed by two presaturation slabs of 40 mm thickness with 8 mm distance on either side of the imaged slice in order to reduce flow artifacts [14, 27, 30]. The duration of the saturation module consisting of the saturation pulses and spoiler gradients was 12 ms. The saturation slabs were applied alternating in subsequent cardiac phases to avoid SAR limitations at 3T [31]. Thus the effective distance between two saturation pulses at the same position was 80 ms. Further SAR optimization was achieved by limiting the B1-amplitude to 8 μT in all experiments. Further reducing the B1-amplitude prolonged the presaturation pulse part of the sequence, whereas higher B1-amplitudes resulted in too high SAR demands.

For respiratory motion compensation, conventional navigator-gating and tracking with an acceptance window of 10 mm was performed applying a pencil beam navigator through the dome of the right hemi-diaphragm [32]. The navigator was applied at each start of the cardiac cycle [33]. The navigator duration was 15.5 ms, the navigator evaluation time was 5 ms. Prospective gating was performed since it reveals superior image quality than retrospective gating [34]. Cardiac triggering was performed using a vector electrocardiogram.

In each volunteer 5 different sequences were performed: two scans without SENSE acceleration for the assessment of the reproducibility and 3 scans with SENSE acceleration factors R of 2,3 and 4. The direction of SENSE was chosen dependent on the orientation of the heart either in RL (8 volunteers) or SI direction (12 volunteers).

For a pulse rate of 60 beats per minute 21 heart phases could be acquired. The nominal imaging times were 3 : 45 minutes without SENSE acceleration, 1 : 57 minutes with a SENSE factor of 2, 1 : 17 minutes with a SENSE factor of 3 and 57 seconds for a SENSE factor of 4, so that the nominal acceleration factors are given by 1.92, 2.92 and 3:95 for a SENSE factor of R = 2, 3 or 4 respectively.

One volunteer underwent a black blood prepared, respiratory gated, cardiac phase resolved 3D anatomical and velocity encoded acquisition. For maximal acceleration a SENSE acceleration factor of R = 4 was chosen. Data were acquired on a 3T whole body MR scanner (Achieva 3.0T, Philips, Best, The Netherlands) with a 32 [2 × 4 ×4] channel phased array cardiac coil. The acquisition protocol is listed in table 1.

Data analysis

The TPM MR images were analyzed by an in-house developed MATLAB software (R2008a; Mathworks, Natick, Mass). For each scan 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 [35, 36].

Before quantification of the resulting myocardial velocities, a background phase error correction was performed using a linear fit to the phase of static tissue [37].

The radial (towards the center of the blood pool) and longitudinal (towards the apex of the heart) velocity curves were calculated. Velocity-time curves for either direction of motion were generated using the average myocardial velocity of the respective slice. Prior to data analysis the velocity data acquired over time were interpolated by cubic splines to provide a continuous velocity profile. Physiologically, the accumulated phase over the entire heart cycle must result to zero. To reduce remaining phase errors after the linear background phase error correction described above, in a subsequent correction step the resulting velocity curves were shifted accordingly to meet the physiological conditions.

The systolic and diastolic peak velocities vp,sys and vp,dias were determined and the velocity range Δ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 SENSE acceleration was denoted as Δv Diff. The peak factor PF seq .  1 ,  seq .  2 = Δ v seq .2 Δ v seq .1 was calculated. The quality of the resulting velocity curves was additionally quantified by the normalized root mean square deviation (nRMSD(seq. 1, seq. 2)) between the velocity curves obtained by the SENSE accelerated and non-accelerated technique and the correlation coefficient s c(seq:1, seq.2). Thereby, the correlation coefficient is given by
c = i = 1 N [ ( v i,seq . 1 v ¯ seq .1 σ seq .1 ) ( v i,seq . 2 v ¯ seq .2 σ seq .2 ) ] ,
(1)

where vi,seq.1 and vi,seq.2 are the spline interpolated data at different time steps (step size: 0.01 ms).

For the regional analysis of velocity information, the velocity range, peak factor an systolic and diastolic peak velocities are determined for six segments of the investigated equatorial slice from anteroseptal to anterior.

The times to the minima radial and longitudinal velocity, tr,dias and tl,dias were determined for each acquisition technique. Bland-Altman analysis was performed for tr,dias and tl,dias and the mean time differences Δ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 assessment of the inherently reduced signal-to-noise ratios (SNR) in parallel imaging, SNR maps of the anatomical and velocity encoded data were calculated. Local SNR values were derived from the mean value of a 3 × 3 matrix divided by its standard deviation. SNR values were calculated for the whole myocardium and the respective regional 6 segments.

For the evaluation of significances a paired two-tailed student's t-test was performed. P-values below 0.05 were considered to be significant. The 3D anatomical + 3D velocity encoded + time acquisition was analyzed regarding the radial and longitudinal motion at different position of the left ventricle.

Results

The scan protocol could be finished in all volunteers. The navigator gating efficiencies of this study were 68.8% ± 14.2% for no SENSE (1), 72.6% ± 12.0% for no SENSE (2), 71.5% ± 12.9% for SENSE with R = 2, 70.0% ± 15.3% for SENSE with R = 3 and 71.9% ± 13.6% for SENSE with R = 4 and did not show any significant differences for all acceleration factors.

Figure 1 shows anatomical and velocity encoded images of a myocardial short-axis scan acquired without SENSE and with SENSE acceleration factors R = 2 - 4 exemplary for one volunteer. Visually, only small intensity differences can be appreciated in the anatomical images, whereas a decrease of the velocity magnitude could only be observed for R = 4 in the velocity encoded images.
Figure 1

Anatomical and phase contrast images of the left ventricle with different SENSE acceleration factors. Anatomical and phase contrast images of a short-axis TPM measurement of the left ventricle without SENSE acceleration and with SENSE acceleration factors 2 - 4. The velocity images are divided into 3 columns for each direction of motion (right-left, anterior-posterior and superior-inferior). Visually only small intensity differences can be revealed in the anatomical images, whereas a decrease of the velocity magnitude was only observed for R = 4 in the velocity encoded images.

Figures 2a) and 2b) show the respective radial and longitudinal velocity curves over time for the investigated sequences. Visually, only for an acceleration factor of four a small decrease in the peak amplitude was observed. Table 2 provides the radial and longitudinal peak velocities separately for systole and diastole and the velocity ranges. Significant differences between the accelerated and non-accelerated values were only obtained for the longitudinal diastolic peak velocity for an acceleration factor of R = 4. Table 3 provides the peak factors PFr and PFl for the different SENSE acceleration factors for the radial and longitudinal velocity of the myocardium. A small but significant (p-value < 0.05) reduction of the PF was obtained in the case of the radial motion for a SENSE accelerating factor R of 2 and 4. The lowest peak-factor value obtained with SENSE was 3% lower than the optimal value of 1.
Figure 2

Radial and longitudinal velocity curves over time exemplary for one volunteer. Radial (a) and longitudinal (b) velocity curves over time for the different sequences no SENSE (1), no SENSE (2), SENSE 2, SENSE 3 and SENSE 4 exemplary for one volunteer.

Table 2

Radial and longitudinal systolic and diastolic peak velocities and velocity ranges for the sequences no SENSE (1), no SENSE (2) and SENSE with R = 2 - 4

seq.

vr,sys

vr,dias

Δvr

vl,sys

vl,dias

Δvl

no SENSE (1)

2.90 ± 0.55

-5.03 ± 0.53

7.93 ± 0.83

6.08 ± 2.24

-8.58 ± 2.58

14.66 ± 4.28

no SENSE (2)

3.02 ± 0.63

-5, 15 ± 0.72

8.17 ± 1.00

6.18 ± 2.36

-8.46 ± 2.69

14.63 ± 4.44

SENSE 2

2.83 ± 0.56

-4.99 ± 0.68

7.82 ± 0.98

6.14 ± 2.26

-8.76 ± 2.72

14.91 ± 4.46

SENSE 3

2.94 ± 0.65

-4.89 ± 0.66

7.83 ± 1.11

6.04 ± 2.21

-8.38 ± 2.61

14.41 ± 4.22

SENSE 4

2.93 ± 0.57

-4.78 ± 0.58

7.71 ± 0.96

6.07 ± 2.21

-8.15 ± 2.29

14.22 ± 4.06

Table 3

PF values of the sequence no SENSE (1) (seq. 1) and the sequences no SENSE (2) and SENSE with R = 2 - 4. (seq. 2)

seq. 1

seq. 2

PFr

ΔPFr

PFl

ΔPFl

no SENSE (1)

no SENSE (2)

1.03

0.09

1.00

0.10

no SENSE (1)

SENSE 2

0.99

0.08

1.02

0.10

no SENSE (1)

SENSE 3

0.99

0.07

0.99

0.12

no SENSE (1)

SENSE 4

0.97

0.07

0.98

0.12

The results of the Bland-Altman analysis of the velocity range Δv are shown in table 4. There were no significant differences between Δvr Diff and Δvl Diff obtained by the reproducibility measurement and the SENSE accelerated measurements.
Table 4

Results of the Bland-Altman analysis of the velocity ranges v r and v l

R(seq. 1)

R(seq. 2)

Δvr Diff

Δvl Diff

1

1

-0.20 ± 0.71

0.02 ± 1.27

1

2

0.16 ± 0.62

-0.25 ± 1.19

1

3

0.14 ± 0.59

0.25 ± 1.54

1

4

0.27 ± 0.56

0.43 ± 1.62

Figures 3 and 4 display the radial and longitudinal velocities over time for each segment exemplary for one volunteer. The segmental analysis reveals some deviations between the different acquisitions. However, all deviations appear to be within the reproducibility of the investigated technique.
Figure 3

Segmental radial velocity curves over time exemplary for one volunteer. Regional radial velocity curves over time for the different sequences no SENSE (1), no SENSE (2), SENSE 2, SENSE 3 and SENSE 4 exemplary for one volunteer.

Figure 4

Segmental longitudinal velocity curves over time exemplary for one volunteer. Longitudinal radial velocity curves over time for the different sequences no SENSE (1), no SENSE (2), SENSE 2, SENSE 3 and SENSE 4 exemplary for one volunteer.

Figures 5 displays the radial and longitudinal velocity ranges and peak factors. Significant deviations were found for the longitudinal (R = 4, anteroseptal) and radial (R = 2, anterolateral) velocity range and peak factor (R = 4, PF r and PF l , anterolateral), and radial (R = 2,4, anterolateral) and longitudinal (R = 4, anteroseptal, anterolateral and anterior) diastolic peak velocities. Thus, for most acceleration factors and segments, the longitudinal and radial velocity ranges and peak factors are not altered.
Figure 5

Longitudinal and radial velocity ranges and peak factors averaged over all volunteers. The longitudinal (a,c) and radial (b,d) velocity ranges (on the top) and peak factors (at the bottom) are displaced for each segment.

Figure 6 displays the radial and longitudinal systolic and diastolic peak velocities for the different segments. The absolute value of peak velocities was only significantly reduced for vr,dias in the anterolateral segment for acceleration factors R = 2, 4 and for vl,dias in the anteroseptal, anterolateral and anterior segment for an acceleration factor of R = 4. The mean nRMSDr and nRMSDl values and their standard deviations are listed in table 5. A significant difference between the obtained mean nRMSD values was not observed.
Figure 6

Longitudinal and radial peak velocities averaged over all volunteers. The longitudinal (a) and radial (b) peak systolic and diastolic velocities are displaced for each segment.

Table 5

nRMSD r and nRMSD l values between the sequence no SENSE (1) (seq. 1) and the sequences no SENSE (2) and SENSE with accelerating factors from 2 - 4 (seq. 2)

R(seq. 1)

R(seq. 2)

nRMSDr[%]

nRMSDl[%]

1

1

5.1 ± 2.4

4.6 ± 1.9

1

2

4.9 ± 1.4

4.4 ± 1.7

1

3

5.4 ± 2.4

4.4 ± 2.4

1

4

5.2 ± 2.3

4.7 ± 1.6

Table 6 provides the correlation coefficients cr and cl between the reproducibility and accelerated scans. No significant difference between these values could be obtained for the different acceleration techniques in the global analysis over all segments. The segmental analysis only reveals a significant difference of the longitudinal correlation coefficient for the anterior segment with an acceleration factor of 4.
Table 6

Correlation coefficients cr and cl between no SENSE (1) (seq. 1) and no SENSE (2), SENSE 2, SENSE 3 or SENSE 4 (seq. 2)

R(seq. 1)

R(seq. 2)

cr

cl

1

1

0.98 ± 0.02

0.98 ± 0.02

1

2

0.98 ± 0.01

0.98 ± 0.01

1

3

0.98 ± 0.02

0.98 ± 0.02

1

4

0.98 ± 0.03

0.98 ± 0.01

The time to the diastolic peak radial and longitudinal velocities was calculated for all volunteers. The results of the Bland-Altman analysis for the mean time differences Δtr,dias and Δtl,dias over all volunteers and their standard deviations are provided in table 7. Significant differences between the relative time differences obtained by the reproducibility measurement and the time differences between the accelerated and non-accelerated SENSE TPM measurement were not observed.
Table 7

Results of the Bland-Altman analysis for Δtr,dias and Δtl,dias

R(seq. 1)

R(seq. 2)

Δtr,dias[ms]

Δtl,dias[ms]

1

1

1.00 ± 9.26

1.72 ± 16.60

1

2

0.48 ± 9.24

1.73 ± 13.12

1

3

2.50 ± 7.70

-0.89 ± 5.45

1

4

2.41 ± 7.51

-1.82 ± 13.60

Figure 7 displays the SNR maps for the anatomical and velocity encoded images. A steady decrease in the SNR with increasing acceleration can be appreciated. Table 8 provides the mean SNR values of the entire myocardium and for all equatorial AHA heart segments. The decrease of SNR was statistically significant. Application of the investigated technique to the three-dimensional coverage of the left ventricle is provided in Figure 8. The continous change of the radial and longitudinal velocity patterns from the apex to the base can be clearly appreciated.
Figure 7

Anatomical and phase contrast SNR maps of the left ventricle with different SENSE acceleration factors. Anatomical and phase contrast SNR maps of the left ventricle without SENSE acceleration and with SENSE acceleration factors 2 - 4. The velocity images are divided into 3 columns for each direction of motion (anterior-posterior, superior-inferior and right-left). Visually the SNR decreases with increasing acceleration factor.

Table 8

Results of the SNR analysis for the anatomical and velocity encoded images

selected segments

no SENSE (1)

no SENSE (2)

SENSE 2

SENSE 3

SENSE 4

anatomical information

all

32.36 ± 4.81

32.21 ± 5.00

27.64 ± 3.48

26.81 ± 4.43

23.91 ± 2.90

anteroseptal

39.23 ± 7.43

38.61 ± 6.77

33.02 ± 4.55

31.98 ± 6.00

28.73 ± 4.68

inferoseptal

37.22 ± 8.72

36.25 ± 8.70

31.21 ± 5.69

31.28 ± 7.08

27.65 ± 5.98

inferior

26.93 ± 5.52

27.06 ± 5.25

24.54 ± 4.19

23.81 ± 5.51

21.25 ± 4.66

inferolateral

31.28 ± 5.74

32.03 ± 7.25

27.73 ± 4.55

25.09 ± 5.81

22.92 ± 4.31

anterolateral

28.68 ± 5.13

28.58 ± 6.60

23.36 ± 4.65

22.67 ± 4.23

20.38 ± 4.81

anterior

31.14 ± 5.58

31.01 ± 5.63

26.09 ± 5.11

25.84 ± 5.47

22.41 ± 4.53

velocity information in AP direction

all

6.73 ± 1.49

6.57 ± 1.28

5.23 ± 1.06

4.76 ± 1.09

4.30 ± 0.75

anteroseptal

5.94 ± 1.54

6.01 ± 1.42

4.96 ± 1.64

4.69 ± 1.38

3.95 ± 1.04

inferoseptal

5.93 ± 1.40

5.66 ± 1.32

4.67 ± 1.22

4.15 ± 0.97

3.92 ± 0.85

inferior

7.21 ± 1.86

6.80 ± 1.56

5.54 ± 1.32

4.79 ± 1.24

4.45 ± 1.26

inferolateral

7.44 ± 2.60

7.24 ± 2.08

5.47 ± 1.89

4.63 ± 1.07

4.32 ± 1.16

anterolateral

6.99 ± 1.80

7.15 ± 1.96

5.50 ± 1.46

5.27 ± 1.67

4.75 ± 1.19

anterior

6.84 ± 2.23

6.55 ± 1.65

5.21 ± 1.46

5.03 ± 1.79

4.42 ± 0.87

velocity information in FH direction

all

7.03 ± 1.61

7.01 ± 1.61

5.80 ± 1.35

5.12 ± 1.36

4.41 ± 0.92

anteroseptal

7.23 ± 1.44

7.09 ± 1.79

5.77 ± 1.62

5.23 ± 1.39

4.56 ± 1.11

inferoseptal

6.41 ± 1.79

6.24 ± 2.01

5.49 ± 1.58

4.68 ± 1.40

4.23 ± 1.18

inferior

7.21 ± 2.06

7.35 ± 1.77

5.90 ± 1.66

5.29 ± 1.32

4.57 ± 1.09

inferolateral

7.40 ± 2.61

7.68 ± 2.54

6.29 ± 1.87

5.47 ± 2.13

4.81 ± 1.49

anterolateral

6.94 ± 3.09

6.66 ± 2.84

5.42 ± 1.76

4.77 ± 2.06

4.09 ± 1.31

anterior

6.94 ± 1.72

6.98 ± 1.87

5.85 ± 1.95

5.18 ± 1.07

4.17 ± 0.99

velocity information in RL direction

all

7.67 ± 1.84

7.54 ± 1.51

5.89 ± 1.42

5.51 ± 1.18

4.94 ± 1.10

anteroseptal

7.47 ± 1.96

7.56 ± 1.77

5.97 ± 1.24

5.91 ± 1.25

4.95 ± 1.28

inferoseptal

6.86 ± 1.90

6.64 ± 2.13

5.64 ± 1.74

5.00 ± 1.37

4.66 ± 1.13

inferior

7.95 ± 2.17

7.59 ± 1.68

6.31 ± 2.18

5.64 ± 1.72

5.29 ± 1.40

inferolateral

8.64 ± 2.29

8.86 ± 1.87

6.60 ± 1.84

5.93 ± 1.64

5.55 ± 1.48

anterolateral

7.58 ± 2.62

7.38 ± 2.38

5.43 ± 1.73

5.11 ± 1.38

4.66 ± 1.58

anterior

7.44 ± 2.44

7.18 ± 1.98

5.32 ± 1.65

5.39 ± 2.03

4.47 ± 1.50

Figure 8

Radial and longitudinal velocity curves in 3D exemplary for one volunteer. Radial (a) and longitudinal (b) velocities over time are displaced from apex to basis with the 3D anatomical + 3D velocity + time data acquisition with a SENSE acceleration factor of 4. A strong variation appears dependent on the position of the left ventricle.

Radial motion: In systole, all slice first move toward the center of the heart. The highest radial velocities occur in equatorial regions. In systole, apical slices perform a single outward motion, whereas for basal segments two peaks of outward motion can be observed.

Longitudinal motion: At the beginning of systole, apical, equatorial and basal slices move towards the apex of the heart. This motion is stronger for basal than for apical slices. At the end of systole, the velocity of the basal slices approach a plateau, whereas most apical slices reach small negative values. During diastole, all slices move back towards the basis. Again, this motion is more pronounced in the basis than in the apex. Finally, a short motion in opposite direction occurs in all regions.

Discussion

The application of SENSE to accelerate the acquisition of TPM data for the quantification of myocardial velocities appears feasible. All investigated SENSE accelerating factors can be applied without fundamental loss of motion information of the myocardium. No significant differences were obtained for the nRMSD, the global correlation coefficient c and the temporal differences Δt dias between the reproducibility and accelerated measurements and the absolute velocity differences Δv for both the radial and longitudinal velocity curves. Significant differences were obtained for the peak-factor. This differences of the peak factor can be caused by the increased noise in the SENSE accelerated data, which is an inherently property of SENSE. Nevertheless, these differences were not clinical relevant (less than 3%). Delfino et al. found a mean decrease in diastolic peak velocity values in heart failure patients between 68% (longitudinal) and 13% (circumferential) in comparison to healthy volunteers [14]. Therefore, the use of SENSE for TPM imaging appears promising. In most segments, the analysis of the velocity ranges, peak factors and peak velocities obtains no significant deviations.

The major limitation of SENSE in tissue phase mapping is the inherent decrease of SNR. The presented data indicate that the reduction in SNR does not cause clinical relevant deviations of the velocity information. Deviations of the velocity range and peak velocities appear to be within the reproducibility of the underlying technique and hence no significant limitations are expected from using SENSE. Whether the intrinsically lower SNR values at lower field strength will pose limitations remain to be investigated. Since in this study SENSE acceleration was restricted to the phase-encoding direction, the results should be representative also for coil arrays with a lower number of elements (e.g. a 6- or 12 channel coil) as long the acceleration factors does not exceed the number of receive coils in phase encoding direction and the resulting SNR is not limited by the performance of the coil.

This high performance of SENSE to measure accurate velocity profiles was also found in previous studies Investigating PC-CMR combined with SENSE in vessels [38, 39]. GRAPPA, which uses correlations ink-space instead of image-space, has also shown high performance in PC-CMR acquisitions in vessels [40]. It can be assumed, that the accuracy of the measurement of myocardial velocity profiles with GRAPPA is similar to the accuracy of SENSE.

One method to even further accelerate TPM data acquisition would be the combination of SENSE with view sharing, thus exploiting correlation in both, k-space and time. Since the gain of acquisition speed in view sharing is based on less data acquisition in the outer part of k-space, no further velocity degradation would be expected.

In this study a temporal resolution of 40 ms was used. Higher temporal resolutions might increase the accuracy in the peak velocity quantification and can for example be obtained with view sharing [19], which also decreases scan time. Since the aim of this study was to investigate SENSE independent on other accelerating techniques, in this study view sharing was not applied.

A combination of k-t BLAST and SENSE for velocity encoded imaging was performed for blood flow measurements in the ascending aorta [22]. With an high acceleration factor of 8 a slow temporal low-pass filtering was obtained, whereas the flow measurement with an acceleration factor of 5 agreed to the reference flow measurement of this study. Further, it has been demonstrated, that there is only minor deviation of the myocardial motion for the k-t GRAPPA approach for an accelerating factor of 6 [41]. Therefore k-t SENSE might be used to further accelerate the TPM data acquisition.

A further method to accelerate TPM data acquisition with even higher acceleration factors would be the combination of SENSE with a temporal constraint k-t BLAST approach using principal component analysis (k-t PCA) [42]. First perfusion experiments in vivo show a higher performance of k-t PCA compared with k-t SENSE for high acceleration factors [42]. Future velocity encoded k-t PCA experiments of the human myocardium are necessary to validate the applicability of this technique.

In this study, the SENSE approach was only applied for accelerating the data acquisition. Alternatively, the gained acquisition time might be used to increase the spatial resolution or to encode more imaging slices. The combination of TPM and SENSE may provide wide application of myocardial motion measurements with sufficient spatial and temporal resolution in research as well as in clinical routine. Additionally it has The potential to enable 3D anatomical + 3D velocity encoded whole heart cine imaging, which is shown exemplary in one healthy volunteer. Here the strong variation of motion pattern from the apex to basis can be assessed. The results are in conformity with earlier studies investigating a few slices in 2D with 3D velocity encoding [43, 44]. This 3D anatomical + 3D velocity encoded whole heart cine imaging may in future be used to classify motion abnormalities in diverse cardiac diseases such as cardiac insufficiency or asynchrony.

Conclusions

In summary, a combination of TPM and SENSE for the acceleration of phase-encoded data acquisition is possible with only negligible motion information loss for all investigated acceleration factors. This enables a decrease of scan duration of 75% or increase of volume coverage of 4.

The SENSE sequence has the potential to enable TPM measurements in clinical routine and 3D whole heart TPM measurements in reasonable image acquisition times. The possible combination of SENSE with imaging techniques using correlations in k-space and time like k-t BLAST may offer a way to further reduce the overall scan time.

Declarations

Authors’ Affiliations

(1)
Department of Internal Medicine II, University Hospital of Ulm
(2)
Philips Research North America
(3)
Medisys Research Lab, Philips Healthcare
(4)
Karlsruhe Institute of Technology

References

  1. l De Boeck BW, Cramer MM, Oh JK, van der Aa RPLM, Jaarsma W: Spectral pulsed tissue Doppler imaging in diastole: a tool to increase our insight in and assessment of diastolic relaxation of the left ventricle. Am Heart J. 2003, 146 (3): 411-419. 10.1016/S0002-8703(03)00322-3.View ArticleGoogle Scholar
  2. Jung B, Markl M, Foll D, Hennig J: Investigating myocardial motion by MRI using tissue phase mapping. Eur J Cardiothorac Surg. 2006, 29 (Suppl-1): S150-157.View ArticlePubMedGoogle Scholar
  3. Mohiaddin RH, Longmore DB: Functional aspects of cardiovascular nuclear magnetic resonance imaging. Techniques and application. Circulation. 1993, 88: 264-281.View ArticlePubMedGoogle Scholar
  4. Reichek N: Magnetic resonance imaging for assessment of myocardial function. Magn Reson Q. 1991, 7 (4): 255-274.PubMedGoogle Scholar
  5. Masood S, Yang GZ, Pennell DJ, Firmin DN: Investigating intrinsic myocardial mechanics: the role of MR tagging, velocity phase mapping, and diffusion imaging. J Magn Reson Imaging. 2000, 12 (6): 873-883. 10.1002/1522-2586(200012)12:6<873::AID-JMRI10>3.0.CO;2-J.View ArticlePubMedGoogle Scholar
  6. Zerhouni EA, Parish DM, Rogers WJ, Yang A, Shapiro EP: Human heart: tagging with MR imaging-a method for noninvasive assessment of myocardial motion. Radiology. 1988, 169: 59-63.View ArticlePubMedGoogle Scholar
  7. Axel L, Dougherty L: Heart wall motion: improved method of spatial modulation of magnetization for MR imaging. Radiology. 1989, 172 (2): 349-350.View ArticlePubMedGoogle Scholar
  8. Osman NF: Detecting stiff masses using strain-encoded (SENC) imaging. Magn Reson Med. 2003, 49 (3): 605-608. 10.1002/mrm.10376.View ArticlePubMedGoogle Scholar
  9. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  10. Pelc NJ, Herfkens RJ, Shimakawa A, Enzmann DR: Phase contrast cine magnetic resonance imaging. Magn Reson Q. 1991, 7 (4): 229-254.PubMedGoogle Scholar
  11. Pelc LR, Sayre J, Yun K, Castro LJ, Herfkens RJ, Miller DC, Pelc NJ: Evaluation of myocardial motion tracking with cine-phase contrast magnetic resonance imaging. Investigative Radiology. 1994, 29 (12): 1038-1042. 10.1097/00004424-199412000-00005.View ArticlePubMedGoogle Scholar
  12. Hennig J, Schneider B, Peschl S, Markl M, Krause T, Laubenberger J: Analysis of myocardial motion based on velocity measurements with a black blood prepared segmented gradient-echo sequence: methodology and applications to normal volunteers and patients. J Magn Reson Imaging. 1998, 8 (4): 868-877. 10.1002/jmri.1880080417.View ArticlePubMedGoogle Scholar
  13. 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 (3): 816-826. 10.1148/radiol.2383041992.View ArticlePubMedGoogle Scholar
  14. 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 (3): 917-925. 10.1148/radiol.2463062155.View ArticlePubMedGoogle Scholar
  15. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  16. Feinberg DA, Hoenninger JC, Crooks LE, Kaufman L, Watts JC, Arakawa M: Inner volume MR imaging: technical concepts and their application. Radiology. 1985, 156 (3): 743-747.View ArticlePubMedGoogle Scholar
  17. Weiss S, Rasche V: Projection-reconstruction reduced [correction of reduces] FOV imaging. Magn Reson Imaging. 1999, 17 (4): 517-525. 10.1016/S0730-725X(98)00212-4.View ArticlePubMedGoogle Scholar
  18. 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 (2): 471-478.View ArticlePubMedGoogle Scholar
  19. Markl M, Hennig J: Phase contrast MRI with improved temporal resolution by view sharing: k-space related velocity mapping properties. Magn Reson Imaging. 2001, 19 (5): 669-676. 10.1016/S0730-725X(01)00386-1.View ArticlePubMedGoogle Scholar
  20. 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. 10.1002/mrm.10611.View ArticlePubMedGoogle Scholar
  21. 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 (2): 86-94. 10.1007/s10334-004-0062-8.View ArticlePubMedGoogle Scholar
  22. 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 (6): 1430-1438. 10.1002/mrm.20730.View ArticlePubMedGoogle Scholar
  23. 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 (4): 817-824. 10.1002/jmri.21706.View ArticlePubMedGoogle Scholar
  24. Lutz A, Bornstedt A, Manzke R, Etyngier P, Nienhaus GU, Rasche V: 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. 2011, 13: 5-10.1186/1532-429X-13-5.PubMed CentralView ArticlePubMedGoogle Scholar
  25. Pruessmann KP, Weiger M, Scheidegger MB, Boesiger P: SENSE: sensitivity encoding for fast MRI. Magn Reson Med. 1999, 42 (5): 952-962. 10.1002/(SICI)1522-2594(199911)42:5<952::AID-MRM16>3.0.CO;2-S.View ArticlePubMedGoogle Scholar
  26. 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. 10.1002/mrm.10171.View ArticlePubMedGoogle Scholar
  27. Felmlee J, Ehman R: Spatial presaturation: a method for suppressing flow artifacts and improving depiction of vascular anatomy in MR imaging. Radiology. 1987, 164 (2): 559-564.View ArticlePubMedGoogle Scholar
  28. Pelc NJ, Bernstein MA, Shimakawa A, Glover GH: Encoding strategies for three-direction phase-contrast MR imaging of flow. J Magn Reson Imaging. 1991, 1 (4): 405-413. 10.1002/jmri.1880010404.View ArticlePubMedGoogle Scholar
  29. Souza SP, Szumowski J, Dumoulin CL, Plewes DP, Glover G: SIMA: simultaneous multislice acquisition of MR images by Hadamard-encoded excitation. J Comput Assist Tomogr. 1988, 12 (6): 1026-1030. 10.1097/00004728-198811000-00021.View ArticlePubMedGoogle Scholar
  30. Nayak KS, Rivas PA, Pauly JM, Scott GC, Kerr AB, Hu BS, Nishimura DG: Real-time black-blood MRI using spatial presaturation. J Magn Reson Imaging. 2001, 13 (5): 807-12. 10.1002/jmri.1112.View ArticlePubMedGoogle Scholar
  31. Lutz A, Bornstedt A, Manzke R, Nienhaus GU, Etyngier P, Rasche V: SAR reduced black-blood cine TPM for increased temporal resolution at 3T. Magma (New York, N.Y.). 2011, 24 (3): 127-135. 10.1007/s10334-010-0242-7.Google Scholar
  32. 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.View ArticlePubMedGoogle Scholar
  33. 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 (4): 411-417. 10.1002/(SICI)1522-2586(200004)11:4<411::AID-JMRI9>3.0.CO;2-B.View ArticlePubMedGoogle Scholar
  34. Du Y, McVeigh E, Bluemke D, Silber H, Foo T: A comparison of prospective and retrospective respiratory navigator gating in 3D MR coronary angiography. Int J Cardiovasc Imaging. 2001, 17 (4): 287-94. 10.1023/A:1011675130904.PubMed CentralView ArticlePubMedGoogle Scholar
  35. Fradkin M, Ciofolo C, Mory B, Hautvast G, Breeuwer M: Comprehensive Segmentation of Cine Cardiac MR Images. MICCAI (1). 2008, 178-185.Google Scholar
  36. Hautvast G, Lobregt S, Breeuwer M, Gerritsen F: Automatic contour propagation in cine cardiac magnetic resonance images. IEEE Trans Med Imaging. 2006, 25 (11): 1472-1482.View ArticlePubMedGoogle Scholar
  37. Walker PG, Cranney GB, Scheidegger MB, Waseleski G, Pohost GM, Yoganathan AP: Semiautomated method for noise reduction and background phase error correction in MR phase velocity data. J Magn Reson Imaging. 1993, 3 (3): 521-530. 10.1002/jmri.1880030315.View ArticlePubMedGoogle Scholar
  38. Thunberg P, Karlsson M, Wigström L: Accuracy and reproducibility in phase contrast imaging using SENSE. Magn Reson Med. 2003, 50 (5): 1061-1068. 10.1002/mrm.10634.View ArticlePubMedGoogle Scholar
  39. Beerbaum P, Körperich 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 (SENSE). Circulation. 2003, 108 (11): 1355-1361. 10.1161/01.CIR.0000087603.97036.C2.View ArticlePubMedGoogle Scholar
  40. Lew C, Alley MT, Spielman DM, Bammer R, Chan FP: Breathheld autocalibrated phase-contrast imaging. J Magn Reson Imaging. 2010, 31 (4): 1004-1014. 10.1002/jmri.22127.PubMed CentralView ArticlePubMedGoogle Scholar
  41. Jung B, Honal M, Ullmann P, Hennig J, Markl M: Highly k-t-space-accelerated phase-contrast MRI. Magn Reson Med. 2008, 60 (5): 1169-1177. 10.1002/mrm.21764.View ArticlePubMedGoogle Scholar
  42. Pedersen H, Kozerke S, Ringgaard S, Nehrke K, Kim WY: k-t PCA: temporally constrained k-t BLAST reconstruction using principal component analysis. Magn Reson Med. 2009, 62 (3): 706-716. 10.1002/mrm.22052.View ArticlePubMedGoogle Scholar
  43. Codreanu I, Robson M, Golding S, Jung B, Clarke K, Holloway C: Longitudinally and circumferentially directed movements of the left ventricle studied by cardiovascular magnetic resonance phase contrast velocity mapping. J Cardiovasc Magn Reson. 2010, 12: 48-10.1186/1532-429X-12-48.PubMed CentralView ArticlePubMedGoogle Scholar
  44. 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 (5): 1033-1039. 10.1002/jmri.20703.View ArticlePubMedGoogle Scholar

Copyright

© Lutz et al; licensee BioMed Central Ltd. 2011

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.