 Research
 Open Access
 Published:
Undersampled Cine 3D tagging for rapid assessment of cardiac motion
Journal of Cardiovascular Magnetic Resonance volume 14, Article number: 60 (2012)
Abstract
Background
CMR allows investigating cardiac contraction, rotation and torsion noninvasively by the use of tagging sequences. Threedimensional tagging has been proposed to cover the wholeheart but data acquisition requires three consecutive breath holds and hence demands considerable patient cooperation. In this study we have implemented and studied kt undersampled cine 3D tagging in conjunction with kt PCA reconstruction to potentially permit for single breathhold acquisitions.
Methods
The performance of undersampled cine 3D tagging was investigated using computer simulations and invivo measurements in 8 healthy subjects and 5 patients with myocardial infarction. Fully sampled data was obtained and compared to retrospectively and prospectively undersampled acquisitions. Fully sampled data was acquired in three consecutive breath holds. Prospectively undersampled data was obtained within a single breath hold. Based on harmonic phase (HARP) analysis, circumferential shortening, rotation and torsion were compared between fully sampled and undersampled data using BlandAltman and linear regression analysis.
Results
In computer simulations, the error for circumferential shortening was 2.8 ± 2.3% and 2.7 ± 2.1% for undersampling rates of R = 3 and 4 respectively. Errors in ventricular rotation were 2.5 ± 1.9% and 3.0 ± 2.2% for R = 3 and 4. Comparison of results from fully sampled invivo data acquired with prospectively undersampled acquisitions showed a mean difference in circumferential shortening of −0.14 ± 5.18% and 0.71 ± 6.16% for R = 3 and 4. The mean differences in rotation were 0.44 ± 1.8° and 0.73 ± 1.67° for R = 3 and 4, respectively. In patients peak, circumferential shortening was significantly reduced (p < 0.002 for all patients) in regions with late gadolinium enhancement.
Conclusion
Undersampled cine 3D tagging enables significant reduction in scan time of wholeheart tagging and facilitates quantification of shortening, rotation and torsion of the left ventricle without adding significant errors compared to previous 3D tagging approaches.
Background
Patients with myocardial infarction (MI), coronary artery disease (CAD) or impaired signal conduction in the heart such as left bundle branch block (LBBB) suffer from reduced cardiac function. In patients undergoing cardiac resynchronization therapy (CRT) [1] it is of great interest to generate a mechanical map indicating local dyssynchrony [2]. Among the various mechanical parameters, cardiac strain and torsion have been shown valuable clinical parameters in patients with MI [3] or aortic stenosis [4, 5]. Cardiovascular magnetic resonance (CMR) offers a functional imaging modality to assess cardiac motion pattern and synchrony of contraction [6, 7] noninvasively.
Up to date multiple approaches have been used to investigate myocardial motion such as displacement encoding with stimulated echo (DENSE) [8, 9], velocity encoding [10–12], tagging by spatial modulation of magnetization (SPAMM) [13] and complementary spatial modulation of magnetization (CSPAMM) [14]. These methods can either be applied in conjunction with two dimensional (2D) or three dimensional (3D) imaging. In 2D acquisitions, multiple slices are imaged along the left ventricle (LV) and strain maps are calculated slicebyslice. Twodimensional acquisitions require additional techniques in order to compensate for throughplane motion. To this end, slice following [15, 16], acquisition of additional orthogonal slices [17] or the encoding of throughplane displacement [18, 19] (zHARP) have been proposed. While reconstruction of 3D strain patterns from two biplanar acquisitions does require interpolation, slicefollowing techniques provide only a projection of true 3D motion onto a 2D subspace. Using zHARP additional gradients in throughslice direction are applied to estimate throughplane displacement by solving a set of linear equations [18].
To circumvent the need for interpolation and slicefollowing techniques, true 3D tagging has been proposed [20]. By applying modulation of magnetization in all three spatial directions, throughplane motion is intrinsically captured and hence no slicefollowing is needed. Furthermore, 3D acquisition provides full LV coverage and yields an intrinsically higher signaltonoise (SNR) compared to 2D imaging. While cine 2D acquisitions easily fit into a single breath hold, 3D acquisitions require multiple breath holds with long duration [20]. Further approaches attempted to shorten scan time to four breath holds of 21 RR intervals each [21] by only sampling around the harmonic peaks in kspace. These efforts, however, led to a lower spatial resolution in the two phase encoding directions. A different approach was proposed by Zhong et al. [22] using a 3D DENSE acquisition covering the entire left ventricle during free breathing. However, acquisition durations of up to 20 min [23] depending on navigator efficiency were reported and hence feasibility for routine clinical use may be questioned. Rutz et al. [24] introduced an accelerated 3D CSPAMM method only requiring three navigator gated breath holds of 18 RR intervals each. This implementation was tested in a clinical study with patients suffering from MI and LBBB [25] providing 3D maps of synchrony and magnitude of contraction. To address issues of different breath hold levels the data obtained from respiratory navigators was used to correct for potential offsets. In their study the respiratory navigator was placed onto the right hemidiaphragm and thereby respiratory induced displacement of the heart was approximated according to a linear relationship between liver and heart displacement [26]. The accuracy of the linear translation from motion of the lung liver interface to the position of the heart is still being debated. Nehrke et al. [27] reported a strong correlation of the displacement of the right hemidiaphragm and the heart, but also found significant subject dependent variability in the correction coefficients especially comparing inspiration and expiration. Subject dependency of the translation of breathing induced liver motion to bulk motion of the heart has been confirmed by Moghari et al. [28].
In general, data acquisition can be accelerated by undersampling in spatial and temporal dimensions. Among the various approaches, two strategies have gained significant attention. Compressed sensing (CS) [29] employs nonlinear reconstruction methods to recover information from randomly or pseudorandomly undersampled data. Inherently, compressed sensing algorithms require incoherent sampling and hence become applicable if a sufficient number of phaseencodes exists. In case of onedimensional tagging preparation, the number of Cartesian phaseencodes orthogonal to the tag direction can be greatly reduced providing an efficient and simple way for scan time reduction. In consequence the degrees of freedom to generate random sampling patterns become very limited and the application of compressed sensing appears less favourable in this particular application.
The second acceleration strategy involves uniform spatiotemporal undersampling in conjunction with linear reconstruction algorithms. In kt BLAST and kt SENSE [30] low spatial but full temporal resolution training data is used to unfold signal aliasing resulting from data undersampling. The drawback of these methods relates to temporal filtering if undersampling rates increase. To address this issue, Principal Component Analysis (PCA) of the spatialtemporal frequency domain data was introduced and results obtained with kt PCA show improved temporal fidelity [31].
Although kt undersampling has been extensively applied in CMR including cine and realtime imaging [32, 33], perfusion [34–39] and phase contrast imaging [40–42], only few attempts of applying kt undersampling to tagging have been reported [43, 44].
The objective of the present study was to implement and test kt undersampled wholeheart 3D CSPAMM tagging for rapid assessment of cardiac motion. The performance is demonstrated on simulated data, data obtained in healthy subjects and in patients with myocardial infarction.
Methods
kt PCA
To reconstruct undersampled 3D tagging data, kt PCA [31] is used. Similar to kt BLAST [30], acquired data is divided into 1) training data with low spatial resolution (in phase encoding directions k_{y} / k_{z}) but full temporal resolution and 2) kt undersampled data with high spatial resolution and full temporal resolution. The training data p_{ train }(x,t) is Fourier transformed to be represented in the spatialtemporal frequency domain (xf). Using principal component analysis p_{ train }(x,f) is then decomposed into a basis of temporally dependent functions b(f_{ j }) corresponding to the principal components (pc) and spatially dependent weighting coefficients w_{ train }(x) in xpc space according to:
The aliased signal at point (x, f_{j}) resulting from Rfold undersampling can be written as:
with w(x_{ i }) denoting the spatial weighting coefficients of the unaliased image. Hence the aliased signals P_{ alias }(x,f) can be expressed as:
Finally, the unaliased spatial weighting coefficients are obtained by solving:
where Θ represents an estimate of w_{ x } from training data, Ψ denotes noise variance, H conjugate transpose and + the MoorePenrose pseudoinverse.
Computer simulation
All simulations were performed in Matlab (The MathWorks, Natick, MA, USA). Three orthogonal stacks with line tagging modulation in readout direction were generated (Figure 1a). The CSPAMM method [14] was simulated to avoid tag line fading. The model consisted of a contracting left ventricle as well as static tissue representing chest wall and liver. Circumferential shortening and rotation as measured in a healthy subject at basal and apical level was linearly interpolated along the longaxis to create threedimensional motion data. Peak circumferential shortening was 18.8% and 17.8% for base and apex, respectively. Peak rotation was −3.2° and 10.3° for base and apex (Figure 1bd). Longitudinal shortening obtained from the same invivo subject was incorporated. Simulations of undersampled data acquisition were compared to fully sampled simulated data sets with equivalent spatial and temporal resolution. The matrix size was set according to practical values [24, 25] (Table 1). Gaussian noise was added to kspace data before undersampling resulting in a SNR of 25 prior to undersampling. Both undersampled and training data were extracted from the computer model (Figure 2b). Undersampling rates of R = 3, 4, 5 and 8 were simulated. In all simulations, five training profiles were used in k_{y} and k_{z} direction resulting in a total of 25 training profiles.
In order to study regional wall motion abnormalities, myocardial infarctions of different severity were simulated. To this end, radial shortening in the lateral sector was changed from 100% (no infarction) to 0% (completely static) in steps of 10% (Figure 3). The reduction of myocardial motion in the infarcted zone was applied transmurally along the entire long axis of the left ventricle. Reduction of radial shortening directly reduced the circumferential contraction. In order to assure a smooth transition between infarcted and healthy tissue, the infarcted tissue was continuously “attached” to the adjacent healthy tissue, by reducing the motion damping factor continuously over a sector of 40° on both ends of the infarcted region.
Invivo measurements
Eight healthy subjects (7 male, age: 27.5 ± 3.5 years) and five patients with myocardial infarction (5 male, age: 54.8 ± 5.9 years, 1 patient with acute myocardial infarction post percutaneous coronary intervention (PCI), 4 patients with chronic myocardial infarction) were studied. Imaging was performed on a 1.5 T Philips Achieva System (Philips Healthcare, Best, the Netherlands) using a 5 channel cardiac receiver array. Written informed consent was obtained from all subjects and the protocol was approved by the institutional review and ethics boards.
Tagged images were acquired using an ECG triggered multishot EPI sequence with sequence parameters listed in Table 2. As illustrated in the sequence diagram in Figure 2 the tagging preparation was applied within 10 ms after the detection of the Rwave. As the first RF pulse used for tagging preparation was applied selectively in phase encoding direction, signal outside the fieldofview was suppressed after CSPAMM subtraction and hence reduced fieldofview acquisition could be performed [16]. The imaging sequence was repeated to acquire three orthogonally tagged stacks. For each stack 21–24 heart phases were recorded.
Fully sampled data was acquired in three consecutive breath holds. To guarantee similar breath hold position for each tagged stack a pencil beam respiratory navigator placed on the right hemidiaphragm was used. The position of the lung liver interface was monitored at the beginning of each breath hold and displacements within a 5 mm gating window were accepted. The breath hold was repeated if the subject’s breath hold level was not within the given window. Retrospective undersampling was applied with reduction factors of R = 3 and 4. Prospectively undersampled data were acquired according to the schematic shown in Figure 2b in an additional breath hold. To facilitate comparison, all imaging parameters except for undersampling factors were kept identical for fully sampled and undersampled data acquisitions. For undersampling factors R = 3 and 4 the corresponding sets of parameters are given in Table 2. Prior to image reconstruction, the 5channel coil array data were compressed into a single virtual coil data set using the array compression method proposed by Buehrer et al. [45].
Data analysis
Epicardial and endocardial contours were manually drawn for each slice of interest and a midmural contour was calculated. This contour was tracked over time using the peak combination harmonic phase (HARP) [46, 47] algorithm implemented in an inhouse software. Rotational motion and circumferential shortening was calculated from tracked contours as previously described by Ryf et al. [48]. Curves of circumferential shortening and rotation over the cardiac cycle were fitted by 8^{th} order polynomials, and maxima were found by estimating roots of the derivative. From the mean rotation over the contour, torsion was derived by taking the difference in peak rotation at the most apical and the most basal level [49]. Analysis of peak circumferential shortening and rotation was performed in six sectors per slice in 9 slices for computer simulation and 8–9 slices in invivo experiments. Comparisons were performed sectorwise and are reported for the entire LV.
For simulation the tracking results from fully sampled data were used as ground truth reference. For invivo measurements the mean of paired data points was used as reference. In order to facilitate comparison the initial contour was kept the same for the fully sampled reference and the retrospectively undersampled data in simulation and invivo. In the comparison of fully sampled data with prospectively undersampled acquisitions, contours were redrawn to account for changes in breath hold position and/or patient motion inbetween scans.
Analysis of simulated data was done by calculating the relative difference defined as: $erro{r}_{\mathit{rel}}=\leftmax\left(motio{n}_{R=1}\right)max\left(motio{n}_{R>1}\right)\right/\leftmax\left(motio{n}_{R=1}\right)min\left(motio{n}_{R=1}\right)\right$. The fully sampled data was used as reference. Relative differences ± one standard deviation are reported in % for circumferential shortening and rotation. In order to estimate the correlation between undersampled and reference data linear regression was performed correlating circumferential shortening and rotation of reference data and undersampled data. Regression slope, offset and the corresponding 95% confidence interval were estimated as well as the correlation coefficient R^{2} and the standard error of the estimate (SEE). The SEE is given in % of the range of motion and is defined as $SEE=\sqrt{\frac{1}{N2}{\displaystyle {\sum}_{i=1}^{N}{\left({x}_{i}{{x}^{\prime}}_{i}\right)}^{2}}}$ with N being the number of points, x_{ i } denoting measurement data and x_{ i }’ corresponding values derived from linear regression. For invivo imaging, peak circumferential shortening and rotation were compared using Bland Altman analysis. Mean differences and the 95% levels of agreement corresponding to 2 standard deviations (2SD) are reported. Linear regression was performed on peak circumferential shortening and rotation. Regression slopes, offsets, correlation coefficients and the SEE are reported.
Comparing reference data and retrospectively undersampled data, the HARP tracking performance of contour coordinates was studied based on contour vertex definitions in polar coordinates. Similar to the analysis of peak circumferential shortening and rotation, Bland Altman and linear regression analysis were performed.
Patient data are reported using bull’seye plots of peak circumferential shortening, and late gadolinium enhancement (LGE). The centre of the bull’seye plots represents the apex and the outer ring the base of the LV. For LGE images the average myocardial signal intensity per slice and sector normalized to the signal intensity measured in infarcted tissue is shown. Sectors were grouped into two groups: 1) nonviable sectors having more than 50% of the area presenting LGE and 2) viable sectors having less or equal to 50% of the area presenting LGE. Peak circumferential shortening was estimated for all sectors in each group and the median, 50^{th} percentiles and 90^{th} percentiles were estimated and presented in boxplots. Statistical significance of the differences in peak circumferential shortening between both groups were estimated by a twotailed Wilcoxon signedrank test. The results were Bonferroni corrected for repeated testing. A pvalue less than 0.05 was considered statistically significant.
Results
Computer simulation
Normalized differences in peak circumferential shortening and rotation between undersampling factors of R = 3, 4, 5 and 8 and fully sampled data used as reference are shown in Figure 4ab. Errors in circumferential shortening and rotation were below 5% for undersampling factors of R = 3–5 and R = 3, 4, respectively. At R = 8 errors in circumferential shortening and rotation were 2.5 and 3.1 times greater compared to data obtained with R = 4. The differences in torsion between undersampled and fully sampled data were −0.13°, 0.42°, 0.95° and 1.85° for R = 3, 4, 5 and 8, respectively.
Figure 4cj shows regression slope, offset, correlation coefficient and SEE for R = 3, 4, 5 and 8. The 95% confidence intervals are presented as error bars for fitted slopes and offsets.
Figure 3 shows the result of motion tracking in the presence of simulated infarction. Peak circumferential shortening in an equatorial slice is plotted for the reference vs. undersampled data. Linear regression was performed, resulting in a slope of 1.01 (±0.04 95% confidence interval) and an offset of −2.01% of the range of contraction (±4.79% 95% confidence interval).
Invivo measurements
Figure 5 shows short axis slices reconstructed from fully sampled and 3 and 4fold retrospectively and prospectively undersampled data. Frames at 27 ms after detection of the Rwave, at end systole (279 ms) and at mid diastole (559 ms) are shown.
Linear regression analysis performed on fully sampled and retrospectively undersampled data is shown in Figure 6. For each regression the equations of the linear fit, the correlation coefficients R^{2} and the SEE are given. Figure 6ad show the analysis of polar coordinates of tracked points and Figure 6eh demonstrate sectorwise comparison of time curves of circumferential shortening and rotation for R = 3 and 4.
In Figure 7 shows the BlandAltman comparison of fully sampled and retrospectively undersampled invivo data for R = 3 and 4. Dashed lines represent the mean difference (light grey) and the 95% limit of agreement (2SD) (black). In Figure 7ad provides the comparison of radial coordinates (r = radius and ϕ = angle) of tracked contour points. Mean differences ± 2SD of r were 0.0 ± 3.3% and −0.2 ± 4.1% for R = 3 and 4 and −0.2 ± 1.9° and −0.3 ± 2.2° for ϕ. In Figure 7eh BlandAltman plots for peak circumferential shortening and peak rotation are given. The mean differences ± 2SD for circumferential shortening were −0.2 ± 4.1% and −0.1 ± 4.2% for R = 3 and 4. For rotation, mean differences ± 2SD were 0.5 ± 1.8° and 0.7 ± 1.7° for R = 3 and 4. The mean differences of torsion were 0.45 ± 2.22° (p = N.S.) and 0.05 ± 2.24° (p = N.S.) for R = 3 and 4, respectively.
Figure 8 shows the resulting BlandAltman analysis of peak motion (ad) and correlation between full time curves of motion (eh) comparing fully sampled data with data acquired with prospective undersampling (R = 3 and 4). Mean differences ± 2SD of peak circumferential shortening were −0.1 ± 5.2% and −0.7 ± 6.2% for R = 3 and 4. Mean differences of peak rotation were found to be 0.44 ± 1.80° and 0.73 ± 1.67° for R = 3 and 4 while differences in torsion were 0.48 ± 4.20° (p = N.S.) and 0.03 ± 4.48° (p = N.S.) for R = 3 and 4, respectively.
Bull’seye plots of peak circumferential shortening and profiles of circumferential shortening and rotation derived from undersampled (R = 4) and fully sampled data are compared in Figure 9 for one healthy subject. Spatially depend differences are not observed.
Figure 10 compares circumferential shortening obtained from undersampled and fully sampled reference data for the five patients in relation the LGE findings. Differences in peak circumferential shortening in sectors with more than 50% of their area presenting LGE and sectors with less than 50% LGE were statistically significant in all patients. Mean difference ± 2SD in torsion between undersampled and fully sampled reference data in patients was 0.72 ± 2.14°. The SNR in the fully sampled invivo data was 36 ± 12 on average.
Discussion
In this work undersampled 3D CSPAMM tagging in combination with kt PCA has been implemented and validated.
The analysis of simulation results has revealed a maximum applicable acceleration factor of four. Generally, good correlation was found between fully sampled and retrospectively undersampled data over the entire left ventricle. Analysis of relative differences in circumferential shortening showed a slight decrease in error going from the apex towards the base of the heart. As motion was normalized for each slice, this observation is associated with the low magnitude of contraction of apical myocardium.
Different transmural extensions of infarcted tissue were investigated by computer simulation. Motion analysis showed a good correlation between undersampled and fully sampled reference data. Transmural reduction of contraction as low as 20% could be distinguished from healthy fully contracting tissue. Although fully infarcted tissue was simulated as well, sectorwise analysis always showed contraction greater than 6%. This is due to the choice of the position and size of the sectors used for analysis. Sectors with noncontracting tissue contained slightly contracting neighbouring tissue, which leads to a nonzero average circumferential shortening within a sector.
In vivo comparison between fully sampled data and retrospectively undersampled data showed good agreement as the mean difference was less than 1% (% of range of motion). Differences between reference data and prospectively undersampled acquisition were found to be larger in comparison to differences seen relative to retrospectively undersampled data. This finding is related to multiple issues. On one hand, reference data and undersampled data were acquired in two consecutive scans. Despite the use of a gating respiratory navigator, different breath hold levels are possible and slices used for analysis might hence differ in position. On the other hand the fully sampled reference data itself was acquired in three consecutive breath holds, which, despite the use of a respiratory gating window of 5 mm, can lead to stack misalignment within the reference data. The resulting circumferential shortening is dependent on the transmural position of the tracked contour. If stack misalignment occurs a contour that appears to be midmural in one stack can be positioned more epicardially or endocardially in the consecutive two stacks from which the motion orthogonal to the first stack is derived from. Hence a contour point is tracked at three different transmural positions. Therefore the motion profiles obtained from fully sampled data, used as gold standard reference in this work, might have been compromised. Image registration was not performed in this study. The implementation of tailored registration algorithms that can handle orthogonal tagging modulations requires further investigation.
Maps of peak circumferential shortening as well as the corresponding motion profiles were found to agree well for reference and undersampled data. The magnitude of circumferential shortening is comparable to previously reported values [24, 25, 50, 51]. Besides potential motion smoothing, expressed by the positive mean difference for rotation, a slice by slice comparison of prospectively undersampled and fully sampled data might have been biased by an offset in breath hold level for both acquisitions.
In this study, standard 2D single slice tagging data were not available for reference. In order to capture throughplane motion and hence make data comparable to 3D tagging, nonstandard extensions of 2D tagging such as slice following [15] or zHARP [18] are required. Accordingly, assessment of error in the present work was relative to fully sampled 3D CSPAMM data only, which is a clear limitation. Nevertheless, there have been previous records of validating 3D CSPAMM on healthy subjects and patients, which may serve as benchmark [24, 25].
Spatial resolution in each of the three orthogonal stacks in 3D CSPAMM may be considered coarse (3.8 × 5.47.2 × 7.2 mm^{3}). The resolution given in readout direction is, however, directly linked to the tagline spacing when using HARP analysis. Higher temporal resolution (< 30 ms) is desirable as this allows separating data in xpc space further and hence improves reconstruction accuracy in kt PCA. However, this results in more heart phases and hence more RF excitations per cardiac cycle reducing the contrasttonoise ratio of the tagged data.
Five patients with myocardial infarction were examined using undersampled tagging and LGE imaging. Maps of circumferential shortening derived from reference data and prospectively undersampled data agreed well and infarcted regions could be localized. Direct correlation of peak circumferential shortening with the area of delayed enhancement was found to be only moderate. A potential limitation of infarct detection is related to averaging of motion within sectors. A sixsector per slice model [24, 25] results in 16.7% of total myocardial mass captured per sector and per slice. For example, nonviable tissue in border zones of infarction is passively moved and compressed and hence the extent of motion abnormality may be overestimated. This issue may be addressed by increasing the number of sectors per slice and by using multiple circumferential contours covering the entire transmural extent in future work.
Several strategies may be envisioned to shorten the relatively long breath hold durations (~20 sec) of the present protocols. First, sampling of training data can be fully integrated into the acquisition of undersampled data. Such a variable density EPI approach provides shorter overall scan duration at the expense of reduced temporal resolution. Second, the separation of training data and undersampled data allows splitting the data acquisition into two breath holds. Compared to the method proposed by Rutz et al. [24] the three stacks are not acquired within three consecutive breath holds, but could be acquired within one breath hold for high resolution undersampled and a second breath hold for low resolution training data. Since the training data has very low spatial resolution in the phaseencode directions (2.3 × 2.3 cm^{2}), differences in breath hold levels between acquisitions no longer pose a concern.
Conclusion
Undersampled cine 3D tagging in conjunction with kt PCA reconstruction enables significant reduction in scan time of wholeheart tagging and facilitates efficient and accurate quantification of shortening, rotation and torsion of the left ventricle. Using 3fold undersampling the entire cine 3D tagging acquisition could be accommodated in a single breath hold and feasibility in volunteers and patients was demonstrated. Future work is dedicated to shorten breath hold durations further and to apply the method in larger patient cohorts to prove clinical value.
Acknowledgments
The authors acknowledge funding by the Swiss National Science Foundation, grant #CR3213_132671/1.
Abbreviations
 CMR:

Cardiovascular magnetic resonance
 MI:

Myocardial infarction
 CAD:

Coronary artery disease
 LBBB:

Left bundle branch block
 CRT:

Cardiac resynchronization therapy
 LV:

Left ventricle
 PCI:

Percutaneous coronary intervention
 LGE:

Lategadolinium enhancement
 DENSE:

Displacement encoding
 SPAMM:

Spatial modulation of magnetization
 CSPAMM:

Complementary spatial modulation of magnetization
 HARP:

Harmonic phase
 PCA:

Principal component analysis
 SNR:

Signal to noise ratio
 SEE:

Standard error of the estimate
 ref:

Reference
 und:

Undersampling
 csh:

Circumferential shortening.
References
 1.
Bilchick KC, Lardo AC: Cardiac resynchronization therapy: application of imaging to optimize patient selection and assess response. Curr Heart Fail Rep. 2008, 5: 11927. 10.1007/s1189700800202.
 2.
Jehle A, Epstein F, Zhong X, Janiczek R, Tsai WK, Christopher J, Fowler D, Ferguson J, Kramer C, Bilchick K: Cine DENSE MRI for circumferential and radial dyssynchrony in patients referred for cardiac resynchronization therapy. J Cardiovasc Magn Reson. 2009, 11: O9010.1186/1532429X11S1O90.
 3.
Nagel E, Stuber M, Lakatos M, Scheidegger MB, Boesiger P, Hess OM: Cardiac rotation and relaxation after anterolateral myocardial infarction. Coron Artery Dis. 2000, 11: 2617. 10.1097/0001950120000500000009.
 4.
Stuber M, Scheidegger MB, Fischer SE, Nagel E, Steinemann F, Hess OM, Boesiger P: Alterations in the local myocardial motion pattern in patients suffering from pressure overload due to aortic stenosis. Circulation. 1999, 100: 3618. 10.1161/01.CIR.100.4.361.
 5.
Sandstede JJ, Johnson T, Harre K, Beer M, Hofmann S, Pabst T, Kenn W, Voelker W, Neubauer S, Hahn D: Cardiac systolic rotation and contraction before and after valve replacement for aortic stenosis: a myocardial tagging study using MR imaging. AJR Am J Roentgenol. 2002, 178: 9538.
 6.
McVeigh ER: MRI of myocardial function: motion tracking techniques. Magn Reson Imaging. 1996, 14: 13750. 10.1016/0730725X(95)02009I.
 7.
Gotte MJ, Germans T, Russel IK, Zwanenburg JJ, Marcus JT, van Rossum AC, van Veldhuisen DJ: Myocardial strain and torsion quantified by cardiovascular magnetic resonance tissue tagging: studies in normal and impaired left ventricular function. J Am Coll Cardiol. 2006, 48: 200211. 10.1016/j.jacc.2006.07.048.
 8.
Aletras AH, Ding S, Balaban RS, Wen H: DENSE: displacement encoding with stimulated echoes in cardiac functional MRI. J Magn Reson. 1999, 137: 24752. 10.1006/jmre.1998.1676.
 9.
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: 101927. 10.1002/jmri.21317.
 10.
Arai AE, Gaither CC, Epstein FH, Balaban RS, Wolff SD: Myocardial velocity gradient imaging by phase contrast MRI with application to regional function in myocardial ischemia. Magn Reson Med. 1999, 42: 98109. 10.1002/(SICI)15222594(199907)42:1<98::AIDMRM14>3.0.CO;2H.
 11.
Robson MD, Constable RT: Threedimensional strainrate imaging. Magn Reson Med. 1996, 36: 53746. 10.1002/mrm.1910360406.
 12.
Jung B, Zaitsev M, Hennig J, Markl M: Navigator gated high temporal resolution tissue phase mapping of myocardial motion. Magn Reson Med. 2006, 55: 93742. 10.1002/mrm.20808.
 13.
Axel L, Dougherty L: MR Imaging of Motion with Spatial Modulation of Magnetization. Radiology. 1989, 171: 8415.
 14.
Fischer SE, McKinnon GC, Maier SE, Boesiger P: Improved myocardial tagging contrast. Magn Reson Med. 1993, 30: 191200. 10.1002/mrm.1910300207.
 15.
Stuber M, Spiegel MA, Fischer SE, Scheidegger MB, Danias PG, Pedersen EM, Boesiger P: Single breathhold slicefollowing CSPAMM myocardial tagging. MAGMA. 1999, 9: 8591.
 16.
Fischer SE, McKinnon GC, Scheidegger MB, Prins W, Meier D, Boesiger P: True myocardial motion tracking. Magn Reson Med. 1994, 31: 40113. 10.1002/mrm.1910310409.
 17.
Moore CC, O'Dell WG, McVeigh ER, Zerhouni EA: Calculation of threedimensional left ventricular strains from biplanar tagged MR images. J Magn Reson Imaging. 1992, 2: 16575. 10.1002/jmri.1880020209.
 18.
AbdElmoniem KZ, Stuber M, Osman NF, Prince JL: ZHARP: threedimensional motion tracking from a single image plane. Inf Process Med Imaging. 2005, 19: 63951.
 19.
AbdElmoniem KZ, Stuber M, Prince JL: Direct threedimensional myocardial strain tensor quantification and tracking using zHARP. Med Image Anal. 2008, 12: 77886. 10.1016/j.media.2008.03.008.
 20.
Ryf S, Spiegel MA, Gerber M, Boesiger P: Myocardial tagging with 3DCSPAMM. J Magn Reson Imaging. 2002, 16: 3205. 10.1002/jmri.10145.
 21.
Ryf S, Baltes C, Boesiger P: High speed 3D CSPAMM. 2004, 12th Annual Meeting of ISMRM, Kyoto, Japan, 657
 22.
Zhong X, Spottiswoode BS, Meyer CH, Kramer CM, Epstein FH: Imaging threedimensional myocardial mechanics using navigatorgated volumetric spiral cine DENSE MRI. Magn Reson Med. 2010, 64: 108997. 10.1002/mrm.22503.
 23.
Zhong X, Spottiswoode B, Meyer C, Kramer C, Epstein F: Navigatorgated 3D cine DENSE: development and initial evaluation. J Cardiovasc Magn Reson. 2009, 11: O9410.1186/1532429X11S1O94.
 24.
Rutz AK, Ryf S, Plein S, Boesiger P, Kozerke S: Accelerated wholeheart 3D CSPAMM for myocardial motion quantification. Magn Reson Med. 2008, 59: 75563. 10.1002/mrm.21363.
 25.
Rutz AK, Manka R, Kozerke S, Roas S, Boesiger P, Schwitter J: Left ventricular dyssynchrony in patients with left bundle branch block and patients after myocardial infarction: integration of mechanics and viability by cardiac magnetic resonance. Eur Heart J. 2009, 30: 211727. 10.1093/eurheartj/ehp212.
 26.
Wang Y, Riederer SJ, Ehman RL: Respiratory motion of the heart: kinematics and the implications for the spatial resolution in coronary imaging. Magn Reson Med. 1995, 33: 7139. 10.1002/mrm.1910330517.
 27.
Nehrke K, Bornert P, Manke D, Bock JC: Freebreathing cardiac MR imaging: study of implications of respiratory motion–initial results. Radiology. 2001, 220: 8105. 10.1148/radiol.2203010132.
 28.
Moghari MH, Hu P, Kissinger KV, Goddu B, Goepfert L, Ngo L, Manning WJ, Nezafat R: Subjectspecific estimation of respiratory navigator tracking factor for freebreathing cardiovascular MR. Magn Reson Med. 2012, 67: 166572. 10.1002/mrm.23158.
 29.
Lustig M, Donoho D, Pauly JM: Sparse MRI: The application of compressed sensing for rapid MR imaging. Magn Reson Med. 2007, 58: 118295. 10.1002/mrm.21391.
 30.
Tsao J, Boesiger P, Pruessmann KP: kt BLAST and kt SENSE: dynamic MRI with high frame rate exploiting spatiotemporal correlations. Magn Reson Med. 2003, 50: 103142. 10.1002/mrm.10611.
 31.
Pedersen H, Kozerke S, Ringgaard S, Nehrke K, Kim WY: kt PCA: temporally constrained kt BLAST reconstruction using principal component analysis. Magn Reson Med. 2009, 62: 70616. 10.1002/mrm.22052.
 32.
Bo Z, Haldar JP, Brinegar C, ZhiPei L: Low rank matrix recovery for realtime cardiac MRI. In Biomedical Imaging: From Nano to Macro, 2010 IEEE International Symposium. 2010, p. 9969.
 33.
Gamper U, Boesiger P, Kozerke S: Compressed sensing in dynamic MRI. Magn Reson Med. 2008, 59: 36573. 10.1002/mrm.21477.
 34.
Plein S, Ryf S, Schwitter J, Radjenovic A, Boesiger P, Kozerke S: Dynamic contrastenhanced myocardial perfusion MRI accelerated with kt sense. Magn Reson Med. 2007, 58: 77785. 10.1002/mrm.21381.
 35.
Otazo R, Kim D, Axel L, Sodickson DK: Combination of compressed sensing and parallel imaging for highly accelerated firstpass cardiac perfusion MRI. Magn Reson Med. 2010, 64: 76776. 10.1002/mrm.22463.
 36.
Jung B, Honal M, Hennig J, Markl M: ktSpace accelerated myocardial perfusion. J Magn Reson Imaging. 2008, 28: 10805. 10.1002/jmri.21543.
 37.
Adluru G, McGann C, Speier P, Kholmovski EG, Shaaban A, Dibella EV: Acquisition and reconstruction of undersampled radial data for myocardial perfusion magnetic resonance imaging. J Magn Reson Imaging. 2009, 29: 46673. 10.1002/jmri.21585.
 38.
Manka R, Vitanis V, Boesiger P, Flammer AJ, Plein S, Kozerke S: Clinical feasibility of accelerated, high spatial resolution myocardial perfusion imaging. JACC Cardiovasc Imaging. 2010, 3: 7107. 10.1016/j.jcmg.2010.03.009.
 39.
Manka R, Jahnke C, Kozerke S, Vitanis V, Crelier G, Gebker R, Schnackenburg B, Boesiger P, Fleck E, Paetsch I: Dynamic 3dimensional stress cardiac magnetic resonance perfusion imaging: detection of coronary artery disease and volumetry of myocardial hypoenhancement before and after coronary stenting. J Am Coll Cardiol. 2011, 57: 43744. 10.1016/j.jacc.2010.05.067.
 40.
Hansen MS, Baltes C, Tsao J, Kozerke S, Pruessmann KP, Boesiger P, Pedersen EM: Accelerated dynamic Fourier velocity encoding by exploiting velocityspatiotemporal correlations. MAGMA. 2004, 17: 8694. 10.1007/s1033400400628.
 41.
Baltes C, Kozerke S, Hansen MS, Pruessmann KP, Tsao J, Boesiger P: Accelerating cine phasecontrast flow measurements using kt BLAST and kt SENSE. Magn Reson Med. 2005, 54: 14308. 10.1002/mrm.20730.
 42.
Stadlbauer A, van der Riet W, Crelier G, Salomonowitz E: Accelerated timeresolved threedimensional MR velocity mapping of blood flow patterns in the aorta using SENSE and kt BLAST. Eur J Radiol. 2010, 75: e1521.
 43.
Peters DC, Epstein FH, McVeigh ER: Myocardial wall tagging with undersampled projection reconstruction. Magn Reson Med. 2001, 45: 5627. 10.1002/mrm.1076.
 44.
Ryf S, Kozerke S, Tsao J, Boesiger P: 3D CSPAMM tagging accelerated with ktBLAST. 2003, 11th Annual Meeting of ISMRM, Toronto, Canada, p. 1565
 45.
Buehrer M, Pruessmann KP, Boesiger P, Kozerke S: Array compression for MRI with large coil arrays. Magn Reson Med. 2007, 57: 11319. 10.1002/mrm.21237.
 46.
Osman NF, Kerwin WS, McVeigh ER, Prince JL: Cardiac motion tracking using CINE harmonic phase (HARP) magnetic resonance imaging. Magn Reson Med. 1999, 42: 104860. 10.1002/(SICI)15222594(199912)42:6<1048::AIDMRM9>3.0.CO;2M.
 47.
Ryf S, Tsao J, Schwitter J, Stuessi A, Boesiger P: Peakcombination HARP: a method to correct for phase errors in HARP. J Magn Reson Imaging. 2004, 20: 87480. 10.1002/jmri.20186.
 48.
Ryf S, Rutz AK, Boesiger P, Schwitter J: Is postsystolic shortening a reliable indicator of myocardial viability? An MR tagging and lateenhancement study. J Cardiovasc Magn Reson. 2006, 8: 44551. 10.1080/10976640600599494.
 49.
Nagel E, Stuber M, Burkhard B, Fischer SE, Scheidegger MB, Boesiger P, Hess OM: Cardiac rotation and relaxation in patients with aortic valve stenosis. Eur Heart J. 2000, 21: 5829. 10.1053/euhj.1999.1736.
 50.
Kramer CM, Rogers WJ, Theobald TM, Power TP, Petruolo S, Reichek N: Remote noninfarcted region dysfunction soon after first anterior myocardial infarction. A magnetic resonance tagging study. Circulation. 1996, 94: 6606. 10.1161/01.CIR.94.4.660.
 51.
Rosen BD, Gerber BL, Edvardsen T, Castillo E, Amado LC, Nasir K, Kraitchman DL, Osman NF, Bluemke DA, Lima JA: Late systolic onset of regional LV relaxation demonstrated in threedimensional space by MRI tissue tagging. Am J Physiol Heart Circ Physiol. 2004, 287: H17406. 10.1152/ajpheart.00080.2004.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interest.
Authors’ contributions
CS: study design, data acquisition, reconstruction, image analysis, manuscript drafting, statistics, and literature research. RM: study design, image analysis, statistics and revision of manuscript. PB: study design, literature research, and revision of manuscript. SK: Study design, reconstruction, literature research and revision of manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Stoeck, C.T., Manka, R., Boesiger, P. et al. Undersampled Cine 3D tagging for rapid assessment of cardiac motion. J Cardiovasc Magn Reson 14, 60 (2012). https://doi.org/10.1186/1532429X1460
Received:
Accepted:
Published:
Keywords
 Motion analysis
 Myocardial tagging
 kt undersampling
 Principal component analysis
 Image reconstruction