Skip to main content

Cardiovascular magnetic resonance myocardial feature tracking using a non-rigid, elastic image registration algorithm: assessment of variability in a real-life clinical setting



Cardiovascular magnetic resonance myocardial feature tracking (CMR-FT) is a promising technique for quantification of myocardial strain from steady-state free precession (SSFP) cine images. We sought to determine the variability of CMR-FT using a non-rigid elastic registration algorithm recently available in a commercial software package (Segment, Medviso) in a real-life clinical setting.


Firstly, we studied the variability in a healthy volunteer who underwent 10 CMR studies over five consecutive days. Secondly, 10 patients were selected from our CMR database yielding normal findings (normal group). Finally, we prospectively studied 10 patients with known or suspected myocardial pathology referred for further investigation to CMR (patient group). In the patient group a second study was performed respecting an interval of 30 min between studies. All studies were manually segmented at the end-diastolic phase by three observers. In all subjects left ventricular (LV) circumferential and radial strain were calculated in the short-axis direction (EccSAX and ErrSAX, respectively) and longitudinal strain in the long-axis direction (EllLAX). The level of CMR experience of the observers was 2 weeks, 6 months and >20 years.


Mean contouring time was 7 ± 1 min, mean FT calculation time 13 ± 2 min. Intra- and inter-observer variability was good to excellent with an coefficient of reproducibility (CR) ranging 1.6% to 11.5%, and 1.7% to 16.0%, respectively and an intraclass correlation coefficient (ICC) ranging 0.89 to 1.00 and 0.74 to 0.99, respectively. Variability considerably increased in the test-retest setting with a CR ranging 4.2% to 29.1% and an ICC ranging 0.66 to 0.95 in the patient group. Variability was not influenced by level of expertise of the observers. Neither did the presence of myocardial pathology at CMR negatively impact variability. However, compared to global myocardial strain, segmental myocardial strain variability increased with a factor 2–3, in particular for the basal and apical short-axis slices.


CMR-FT using non-rigid, elastic registration is a reproducible approach for strain analysis in patients routinely scheduled for CMR, and is not influenced by the level of training. However, further improvement is needed to reliably depict small variations in segmental myocardial strain.


Evaluation of cardiac performance is crucial in daily clinical practice to evaluate disease severity, to assess therapeutic interventions and to predict outcome [1]. The lack of easy-to-determine indices of myocardial contractility imposed the cardiology community to consider ejection fraction (EF) the reference standard of ventricular systolic function despite its inability to evaluate regional contractility and its poor sensitivity to detect subtle alterations in ventricular function [2]. Myocardial strain, defined as the measurement of the fractional change of myocardial dimension in a specific direction, appears a superior parameter to assess both global and regional myocardial function. Furthermore, it has been shown to be an early marker of systolic dysfunction and to precede decline in ejection function [3].

Over the years, cardiovascular magnetic resonance (CMR) has emerged as the reference standard for the evaluation of ventricular morphology and function, including myocardial strain [46]. CMR tagging, introduced about 25 years ago, enables noninvasive quantification of myocardial deformation and strain calculation [7, 8], but due to the need of dedicated tagging sequences and time-consuming analysis has not gained widespread acceptance in the clinical arena [9, 10].

Feature tracking (FT) has recently been introduced for analysis and quantification of myocardial strain based on cine steady-state free-precession (SSFP) images, which are part of a standard study CMR protocol [1113]. Although the feature tracking algorithm was originally developed for echocardiographic imaging studies, it has recently been introduced for analysis and quantification of myocardial strain based on cine steady-state free-precession (SSFP) images [1113]. Important advantages are that this technique does not require additional time consuming sequences (as for CMR tagging) but can be applied to standard cine SSFP exams thereby leaving the potential to retrospectively analyze strain in patients where no specific strain image acquisition sequence was performed. Moreover, it has the advantage that the analyses can be performed in a retrospective manner. The CMR-FT algorithm is based on optical flow technology as commercially introduced by TomTec (Unterschleissheim, Germany). More recently, Circle (cvi, Calgary, Canada) has provided similar CMR-FT software. These algorithms use the left ventricular (LV) boundaries to estimate myocardial strain and strain rates in short (SA) and long-axis (LA) direction. In the present study, we have used a non-rigid, elastic image registration algorithm available from Segment (Medviso, Lund, Sweden) [14, 15]. Although reproducibility values for CMR-FT have been published, all used the optical flow algorithm [1622]. Accordingly, the aim of the present study was to investigate not only intra- and inter-observer variability but also inter-study reproducibility of CMR-FT in a setting reflecting the real clinical practice. Hereto, both subjects with and without pathologic findings at CMR were included and observers with a different level of expertise were involved.


Study population

Part 1/One healthy male volunteer (age 43) underwent 10 consecutive CMR studies over a period of 5 days, i.e., two CMR studies per day. Although per day, the CMR studies were subsequently performed, for the second CMR study the volunteer was re-installed and a new exam was started obtaining new localizers and a new determination of cardiac axes was performed. No Gadolinium-based contrast agent was administered.

Part 2/From the UZ Leuven patient CMR database (University Hospitals Leuven, Belgium) we randomly selected 10 patients (7 male) with suspected cardiomyopathy but normal CMR findings (normal group). Mean age 37 ± 11 years; mean LVEF 57 ± 5%. Each patient underwent a complete CMR study with intravenous administration of Gadolinium-based contrast agent (Gadovist, Bayer).

Part 3/Ten patients (6 male) with known or suspected ischemic or non-ischemic myocardial disease were prospectively enrolled (patient group). Each patient was scanned twice with an interval of 30 min between the two CMR studies after re-installing the patient on the scanner. Gadolinium-based contrast agent (Gadovist, Bayer) was injected once, i.e., in the first exam.

Renal function (eGFR) was checked in all subjects prior to contrast agent administration. Exclusion criteria were standard contra-indications to MRI.

Cardiovascular magnetic resonance acquisition

All studies were performed on a 1.5 T unit (Ingenia; Philips Healthcare, the Netherlands) by using commercially available CMR software, electrocardiographic triggering, and a cardiac-dedicated phase-array coil. For the assessment of LV dimensions and function bSSFP breath-hold cine images were acquired in the following orientations: vertical long axis, horizontal long axis and short axis. Standard parameters were repetition time/echo time 3.6/1.8 ms; sense factor 2, flip angle, 60°; section thickness, 8 mm; matrix, 160 × 256; field of view, 300 mm; pixel size, 1.6 × 1.6 mm; and number of phases, 30, phase percentage, 67%. The vertical long axis was determined on the transverse images by positioning an image plane connecting the middle of the mitral valve with the LV apex. On the vertical long axis, the same anatomical landmark points were used to define the horizontal long axis plane. The cardiac short axis was defined on the horizontal long axis using an image plane perpendicular to the interventricular septum. Care was taken to position the most basal short axis slice at end-diastole exactly through the mitral valve ring. The set of short axis images encompassed the left ventricle entirely. Between slices a gap of 2 mm was used. Per breath hold two short axis cine CMR were acquired. For the Late Gadolinium Enhancement (LGE) studies, a dose of 0.15 mL of gadobutrol (Gadovist, Bayer) per kilogram of body weight, was administered.

Left ventricular volumes and function

The image analysis was performed offline by a single observer with more than 20 years of experience in CMR. LV volumes and function were measured using manual planimetry of the endocardial and epicardial borders from the short-axis set, in accordance with approved protocols using a commercially available software package (ViewForum; Philips Medical Systems, Best, the Netherlands) (see online data material).

Feature tracking

Myocardial strain analysis was performed using Segment Medviso software. This software estimates myocardial strain curves by computing inter-frame deformation fields using a tracking strategy based on non-rigid image registration [14, 15]. Specifically, the deformation field is modeled using a B-spline tensor product transform and is found by deforming one image to the next image in the cine sequence. As such, a dense map that describes the difference between consecutive images for all image positions is obtained. The deformation field is computed through an iterative process guided by an intensity-based similarity metric combined with a regularization term, which enforces spatial smoothness of the recovered deformation field. Note that, instead of myocardial boundaries tracking only, the proposed method use the entire image content (i.e. blood pool, entire myocardium) during the optimization process. Furthermore, in order to increase the robustness of the technique to bad image quality, noisy image and image artifacts, a temporal coherence strategy was proposed. Indeed, the tracking methodology was reformulated to take the temporal information of the entire cardiac sequence into account at once. As a result, a strategy that efficiently penalizes non-smooth deformations while keeping a high enough flexibility to extract the relevant physiological deformation was obtained. Fig. 1 shows an example of the tracking result in a normal patient. Circumferential and radial LV myocardial Lagrangian strain were evaluated on short-axis cine SSFP CMR. Longitudinal Lagrangian strain was derived from both horizontal and vertical long-axis cine SSFP CMR. First, both endo- and epicardial contours were manually drawn at end diastole in the long and short-axis by three readers blinded to patients details and to the results of the other readers. Note that, both experts used the same strategy to perform the delineation. Furthermore, it should be noted that the initial contour defines the relevant components of the dense deformation field for computation of the myocardial strain. Contours were propagated automatically by the software throughout the cardiac cycle generating myocardial strain and strain rate curves. No drift correction was performed. Note that, in case the tracking was suboptimal, the original contour (usually the end-diastolic image) could be manually adjusted and re-propagated without the need to rerun the tracking algorithm completely. Subsequently, a reference axis was manually positioned on the mid-septum in order to establish left ventricle segmentation using the 16-segment approach as defined by the American Heart Association [23]. The level of CMR experience varied between the three readers, expert observer, i.e., >20 years of CMR experience (SCMR level III) and skilled observer i.e., 6 months of CMR experience (EACVI CMR level 2) and a beginner i.e., two weeks of training in CMR analysis (basic knowledge of cardiac anatomy, cardiac imaging planes, and training in cardiac image contouring). To assess the intra-observer variability, the skilled reader re-assessed the CMR studies one week after the first reading, unaware of the initial study results.

Fig. 1
figure 1

Feature tracking in a patient belonging to the normal group in short-axis (upper panels) and horizontal long-axis (lower panels). Manual delineation of the endo- and epicardial contour at the end-diastole (left panels). Next, the software automatically deforms these contours using a dense motion field (yellow arrow) estimated between consecutive frames


Statistical analysis was performed using IBM SPSS Statistics version 20 and Microsoft Excel 2013. All variables with normal distribution are presented as mean ± standard deviation. In the short-axis direction, only slices that were contoured by both observers, or during both readings were considered for calculation of the variability. Next, the position of the papillary muscles was used to identify the basal, mid-ventricular and apical level.

In order to assess inter-study, inter- and intra-observer variability, the intraclass correlation coefficient (ICC), Bland-Altman analysis, coefficient of repeatability (CR) and the coefficient of variation (CV) were computed using the strain value at end-systole (ES). Note that ES timing was defined for each dataset by one expert. The analysis was performed using segmental and global strain values. Segmental values were defined using the AHA 16 segments model [23], while global strain was computed as the mean value of all segments in each level. Note that for the radial and circumferential strain computation only the short-axis slices delineated by all the readers were considered. Long-axis strain was computed using horizontal and vertical long axis acquisitions. For the Bland-Altman analysis, both bias (mean differences) and limits of agreement (LOAs, 1.96 times the standard deviation of the differences) were computed. The CV was defined as the standard deviation of the differences (between multiple readings of the same case) divided by the mean value of all samples [22]. CR is computed as 2.77 multiplied by the standard deviation of the differences between readings, as described in [24]. Regarding the healthy volunteer experiment (10 acquisitions), the CV, bias and LOA was computed using the difference between each acquisition against the first one. A two-tailed paired t-test against zero was used to check for statistically significant differences (p < 0.05) of the observer biases. Moreover, the inter-study, inter- and intra-observer variability was compared through two-tailed F-test’s (p < 0.05). The D’Agostino-Pearson test was used to determine normality of the data. The global strain values in both normal and patient group were compared with LV ejection fraction, using Pearson correlation coefficient.


Demographics of the healthy volunteer, and the subjects in the normal and patient group are summarized in the online addenda (Additional file 1). CMR in the patient group, consisting of a mix of patients with acute and chronic ischemic and non-ischemic myocardial disease, yielded normal findings in two patients (i.e., referral for ventricular extrasystoles and familial history of hypertrophic cardiomyopathy). One long-axis cine study in the patient group was because of motion artefacts not useful for analysis. Time for contouring a complete dataset (horizontal long axis - HLA, vertical long axis – VLA, cine and complete SA cine stack) was 7 ± 1 min, and did not show differences between observers (p = 0.56). Myocardial strain calculation time (Dell XPS 8500 Desktop, Intel i7, 16Gb) for this dataset was 13 ± 2 min.

Mean strain values are shown in Table 1. With regard to variability and myocardial strain, a/the variability was systemically the largest for radial strain, b/variability systemically increased from global to segmental level, and c/variability increased from midventricular to the basal and apical LV level (see Tables 2, 3, 4 and 5 and Fig. 2). For instance, CR for global radial and circumferential strain in the healthy volunteer were 3.6% and 2.6%, 8.8% and 3.7%, and 12.0% and 3.7% for intra-observer, inter-observer and inter-study variability, respectively. Similar values were found in the normal and patient group. Compared to global strain variability, segmental strain variability almost invariably showed two-to threefold higher CR values, e.g. CR for intra-observer variability in the patient group ranged 1.6% to 2.7% for global strains versus 4.7% to 10.1% for regional strain estimation (Table 2). As shown in Fig. 3 and in Table 6, the lowest standard-deviations for variability for radial and circumferential strain were in the midventricular part of the left ventricle but considerably enlarged towards the LV base and LV apex. For segmental longitudinal strain variability, a less consistent pattern was found (Table 7 and Fig. 3).

Table 1 Mean (%) and standard deviation (STD, %) myocardial strain values (at end-systole) in the healthy volunteer, the normal and patient group
Table 2 Intra-observer variability (by the skilled observer)
Table 3 Inter-observer variability in three observers with different level of expertise
Table 4 Inter-study variability in the healthy volunteer
Table 5 Inter-study variability in the patient group
Fig. 2
figure 2

Reproducibility study using the multiple acquisitions of the patient group. Bland-Altman plots for global and segmental ES (a, d) radial, (b, e) circumferential and (c, f) longitudinal strain obtained. Dashed lines represent bias (red) and 95% limits of agreement (b)

Fig. 3
figure 3

Mean strain and variability per slice level and per segment in the healthy volunteer (10 acquisitions). a, d represents the global strain components and b, e the segmental strain. c shows the obtained global radial and circumferential strain curves in a mid-ventricular slice for all the 10 acquisitions. f presents the obtained strain curve for the global longitudinal strain

Table 6 Radial and circumferential strain values obtained for the normal group
Table 7 Longitudinal strain values obtained for the normal group

Inter-study variability contributed the most to the overall variability. Moreover, for inter-study variability a small, but significant bias was found in all groups for radial strain (i.e.,−2.25% to−2.53%) independent of the level of expertise of the observer (Table 5 and Fig. 2). In contrast, intra-observer variability as shown in Table 2 was low. The impact of the level of expertise on variability is shown in Tables 3, 4 and 5. Remarkably, the inter-observer variability between a skilled reader and an expert was not substantially different from the variability between a beginner and an expert (Table 3). A similar pattern was found with regard to differences in inter-study variability between observers in the healthy volunteer and patient group (Tables 4 and 5). This is also reflected in Fig. 4, showing that the differences in strain-time curves were primarily caused by test-retest variability and not by the level of expertise of the observer. As can be appreciated in Figs. 5 and 6, showing a tracking result in a normal subject and a patient with dilated cardiomyopathy, the level of experience has a limited impact of tracking of endo- and epicardial contours over the cardiac cycle.

Fig. 4
figure 4

Strain curve examples from different subjects in the patient group. The continuous line represents the first acquisition, while the dashed line represent the second acquisition. Blue shows the result obtained by the Expert observer, red the result obtained with skilled observer and green the result achieved by the beginner. In order to ease the visualization, the absolute value of the longitudinal strain along the vertical long axis is used

Fig. 5
figure 5

Feature tracking in cardiac short-axis in a subject belonging to the normal group by an expert, a skilled observer, and a beginner. Five time points over the cardiac cycle are shown

Fig. 6
figure 6

Feature tracking in cardiac short-axis in a patient with dilated cardiomyopathy by an expert, a skilled observer, and a beginner. Five time points over the cardiac cycle are shown

Compared to a selected group of patients showing no obvious findings at CMR, the variability, in a group of patients presenting both focal and diffuse myocardial pathology, yielded comparable intra- and inter-observer values, e.g., intra-observer CR for global myocardial strain ranging 1.7% to 2.6% in the normal group, versus 1.6% to 2.7% in the patient group. Finally, as shown in Fig. 7, a moderate correlation was found between global myocardial strain and LV EF ranging from R = 0.70 for longitudinal strain, R = 0.86 for radial strain, to R = 0.91 for circumferential strain.

Fig. 7
figure 7

Correlation between myocardial strain (patient and normal group) and LV ejection fraction. a) radial strain, b) circumferential strain, c) longitudinal strain


In the present study, we tested the variability of a novel CMR-FT software package for myocardial strain calculation using a non-rigid, elastic image registration algorithm in a group of subjects similar to a real-life clinical setting. Moreover, we tested the influence of the level of CMR expertise on study reproducibility. Overall, intra- and inter-observer and inter-study variability in terms of LOA, CV, CR and ICC yielded modest to excellent results, comparable or even superior to values reported in literature using optical flow-based CMR-FT software [1622]. Our study results emphasize CMR-FT is a reproducible technique to study global myocardial strain and can be applied in patients with a wide variety of acute and chronic, ischemic and non-ischemic myocardial disease without significant increase in measurement variability. Moreover, the novel CMR-FT algorithm is not influenced by the level of experience, thus applicable by CMR personnel with a basic training in CMR. Time for manual contouring endo- and epicardial contours (at end-diastole) is brief and as the contours are reliably propagated over the remainder of the cardiac cycle, strain versus time curves can be acquired allowing to appreciate myocardial strain patterns during systole and diastole throughout the left ventricle. Currently, strain computation time is still considerable (13 ± 2 min). However, as strain computation is performed in the background, it does not interfere with image contouring. Moreover, further shortening of analysis time with new software releases would facilitate integration of CMR-FT in daily clinical practice.

Cine SSFP CMR has become the standard for quantification of ventricular volumes and myocardial mass [4]. Although the superior contrast between blood and surrounding myocardium with SSFP enables excellent anatomic depiction of the cardiac chambers, SSFP-based myocardial strain analysis in a routine setting has been hampered by time-consuming manual delineation of cardiac contours over the cardiac cycle and the lack of reliable (semi)-automated analysis packages. The advent of CMR-FT software opened the door towards a more routine use of LV strain analysis [1113]. Currently, the approaches for myocardial motion and deformation estimation are threefold: optical flow-based formulation, non-rigid elastic image registration and model-based [14, 15]. Whereas TomTec (Unterschleissheim, Germany) and Circle (cvi, Calgary, Canada) software use an optical flow-based algorithm, we used a new non-rigid, elastic image registration algorithm employing image warping techniques to estimate LV cardiac motion between subsequent frames. The myocardial deformation field is parametrized using B-splines, but additional conditions are imposed to regularize the obtained deformation field in order to avoid warping solutions that are not desired or physiologically not possible. Although this algorithm was originally developed for three-dimensional myocardial deformation analysis from cardiac ultrasound, this algorithm can be applied on CMR images as well [25, 26].

Although the number of scientific papers using CMR-FT for myocardial strain analysis in normal and patients is rapidly growing, only a few papers by a relatively limited number of research groups have assessed the accuracy and variability of SSFP-based FT analysis comparing this technique to other validated CMR techniques such as SPAMM-tagging and DENSE, or to speckle-tracking echocardiography [1622, 26, 27]. In brief, with regard to variability, these studies have shown that variability of CMR-FT equals or is superior to the other CMR and non-CMR approaches, is not influenced by magnetic field strength nor by diurnal period, while intravenous administration of gadolinium-based contrast agents significantly affects variability [19, 21, 22, 2731]. Invariably, all studies so far have reported a substantial increase in variability at segmental compared to global level, with the lowest variability in the mid segments in short-axis direction [19, 22, 32]. Unfortunately, only a few studies so far has addressed the issue of test-retest variability as well [16, 21, 32]. Compared to reported values in literature, similar or better values for CV and ICC, and smaller LOA were found for intra-observer, inter-observer and inter-study variability in both the normal and patient group in the present study. In particular, variability for radial strain estimation was superior to optical flow-based estimates. For example, intra-observer and inter-observer coefficient of variation for global radial strain in the patient group was 4.1% and 8.6%, respectively versus 8.9% to 20.9% and 12.8% to 19.6% respectively for optical-flow based estimates [20, 33]. Intraclass correlation coefficients were 0.99 and 0.97 versus 0.75 to 0.92 and 0.55 to 0.96, respectively [20, 33]. These findings suggest that a non-rigid registration of the myocardial deformation may be more appropriate for depicting changes in radial direction than the current optical flow-based methods. Indeed, the superior performance found for the radial strain can be explained by the high number of myocardial image samples used during the non-rigid method (i.e. entire image content) and the regularization strategies applied.

Another remarkable finding in the present study is that variability does not significantly increase when applied to a real patient group showing a variety of focal and diffuse myocardial pathology. Moreover, the level of expertise does not have a major impact of myocardial strain variability. Previous work by Schuster et al. has shown that observer dependence is reduced when using CMR-FT rather than visual analysis for interpreting low-dose dobutamine stress CMR in patients with ischemic cardiomyopathy [34]. Although there was a clear difference in CMR experience (two versus seven years of experience) in their study, we opted to include an observer with a basic knowledge in cardiac imaging, opposed to a skilled and an expert observer. Our data show that CMR-FT analyses can be reliably performed by people with a basic training in cardiac imaging (basic knowledge of cardiac anatomy, cardiac imaging planes, and knowledge about the peculiarities of cardiac contouring).

Whereas global strain variability was good to excellent, segmental strain variability was moderate, questioning it use to reliably depict small changes in regional myocardial strain. For example CR for segmental radial strain in the patient groups was as high as 29.1% (Table 1). At this point, the non-rigid image registration algorithm is not superior to the optical-based approach [19, 22, 32]. Apparently, all current FT software encounter problems to reproducibly track in particular the basal and apical short-axis slices. Moreover, as our software necessitates a visual definition of the mid-part of the interventricular septum to define the 16-segment model, and not the insertion points of the right ventricle as anatomic landmarks as done by other vendors, small differences in position between readings, readers or studies may negatively impact segmental strain variability.

A systematic bias in radial strain for inter-study variability in the patient group was observed. As the second CMR study was repeated with a delay of 30 min between studies, wash-out of contrast may have interfered with delineation of the endocardial contour and thus with study variability as recently reported by Kuetting et al. [31].

Although in this study we used CMR-FT for left ventricular strain analysis, the current approach can be also applied for all cardiac chambers including the atria. Since the non-rigid method estimates the motion field between consecutive frames for all image positions instead of specific locations (e.g. myocardial boundaries), no modification is required in the current framework to evaluate strain in the remaining chambers. An example of CMR-FT of the right ventricle and left atrium is shown as additional movie video (see Additional files 2, 3 and 4).

Additional file 2: A movie showing feature tracking applied to the right ventricle (short-axis SSFP cine CMR) in a healthy volunteer. (AVI 1231 kb)

Additional file 3: A movie showing feature tracking applied to the right ventricle (horizontal long-axis SSFP cine CMR) in a healthy volunteer. (AVI 2382 kb)


Several limitations need to be mentioned. The aim of the present study was to report intra-observer, inter-observer and inter-study variability of a novel algorithm for CMR-FT in a group of subjects representing a real-life situation. Awaiting a head-to-head comparison between this algorithm and optical flow-based algorithms, we provided for comparison of our data some variability values as reported in literature. However, as the study groups and approaches may differ, any comparison should be taken with the necessary caution. Moreover, at least two groups have recently published age- and gender related values for CMR-FT based myocardial strain [30, 31]. The current values from the ‘normal’ group cannot be considered as reference values because all subjects had complaints and were therefore referred to CMR. Nevertheless, these strain values compare rather well with the “normal” healthy volunteers as reported by Schuster et al. [19].


We tested a non-rigid, elastic image registration algorithm for CMR-FT-based myocardial strain analysis in a real-life clinical setting. Variability values seem equally or superior than those reported in literature with optical flow-based approaches. Moreover, strain analysis can be reproducibly performed by observers with a different level of CMR expertise in patients presenting a wide range of myocardial diseases. Further improvement, however, is needed to reliably depict subtle changes in segmental strains.



Cardiovascular magnetic resonance


Cardiovascular magnetic resonance – feature tracking


Coefficient of repeatability


Coefficient of variation

EccSAX :

Circumferential strain in short-axis

EllLAX :

Longitudinal strain in long-axis

ErrSAX :

Radial strain in short-axis


Feature tracking


Intraclass correlation coefficient


Limits of agreement


Left ventricle


Steady-state free-precession


  1. Smiseth OA, Torp H, Opdahl A, Haugaa KH, Urheim S. Myocardial strain imaging: how useful is it in clinical decision making? Eur Heart J. 2016;37:1196–207.

    Article  PubMed  Google Scholar 

  2. Kalam K, Otahal P, Marwick TH. Prognostic implications of global LV dysfunction: a systematic review and meta-analysis of global longitudinal strain and ejection fraction. Heart. 2014;100:1673–80.

    Article  PubMed  Google Scholar 

  3. Thavendiranathan P, Poulin F, Lim KD, Plana JC, Woo A, Marwick TH. Use of myocardial strain imaging by echocardiography for the early detection of cardiotoxicity in patients during and after cancer chemotherapy: a systematic review. J Am Coll Cardiol. 2014;63:2751–68.

    Article  PubMed  Google Scholar 

  4. Ichikawa Y, Sakuma H, Kitagawa K, Ishida N, Takeda K, Uemura S, Motoyasu M, Nakano T, Nozaki A. Evaluation of left ventricular volumes and ejection fraction using fast steady-state cine MR imaging: comparison with left ventricular angiography. J Cardiovasc Magn Reson. 2003;5:333–42.

    Article  PubMed  Google Scholar 

  5. Bogaert J, Rademakers FE. Regional nonuniformity of normal adult human left ventricle. Am J Physiol Heart Circ Physiol. 2001;280:H610–620.

    CAS  PubMed  Google Scholar 

  6. Edvardsen T, Gerber BL, Garot J, Bluemke DA, Lima JA, Smiseth OA. Quantitative assessment of intrinsic regional myocardial deformation by Doppler strain rate echocardiography in humans: validation against three-dimensional tagged magnetic resonance imaging. Circulation. 2002;106:50–6.

    Article  PubMed  Google Scholar 

  7. 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.

    Article  CAS  PubMed  Google Scholar 

  8. Axel L, Dougherty L. MR imaging of motion with spatial modulation of magnetization. Radiology. 1989;171:841–5.

    Article  CAS  PubMed  Google Scholar 

  9. Hor KN, Gottliebson WM, Carson C, Wash E, Cnota J, Fleck R, Wansapura J, Klimeczek P, Al-Khalidi HR, Chung ES, Benson DW, Mazur W. Comparison of magnetic resonance feature tracking for strain calculation with harmonic phase imaging analysis. JACC Cardiovasc Imaging. 2010;3:144–51.

    Article  PubMed  Google Scholar 

  10. Ibrahim e-SH. Myocardial tagging by cardiovascular magnetic resonance: evolution of techniques–pulse sequences, analysis algorithms, and applications. J Cardiovasc Magn Reson. 2011;13:36.

    Article  PubMed Central  Google Scholar 

  11. Maret E, Todt T, Brudin L, Nylander E, Swahn E, Ohlsson JL, Engvall JE. Functional measurements based on feature tracking of cine magnetic resonance images identify left ventricular segments with myocardial scar. Cardiovasc Ultrasound. 2009;7:53.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Hor KN, Baumann R, Pedrizzetti G, Tonti G, Gottliebson WM, Taylor M, Benson DW, Mazr W. Magnetic resonance derived myocardial strain assessment using feature tracking. J Vis Exp. 2011;48:e2356.

    Google Scholar 

  13. Schuster A, Hor KN, Kowallick JT, Beerbaum P, Kutty S. Cardiovascular magnetic resonance myocardial feature tracking. Concepts and clinical applications. Circ Cardiovasc Imaging. 2016;9(4):e0004077. doi:10.1161/CIRCIMAGING.115.004077.

    Article  Google Scholar 

  14. Heyde B, Bouchez S, Thieren S, Vandenheuvel M, Jasaityte R, Barbosa D, Claus P, Maes F, Wouters P, D’hooge J. Elastic image registration to quantify 3-D regional myocardial deformation from volumetric ultrasound: experimental validation in an animal model. Ultrasound Med Biol. 2013;39:1688–97.

    Article  PubMed  Google Scholar 

  15. Morais P, Heyde B, Barbosa D, Queirós S, Claus P, D’hooge J. Cardiac motion and deformation estimation from tagged MRI sequences using a temporal coherent image registration framework. Springer, Berlin Heidelberg: International Conference on Functional Imaging and Modeling of the Heart; 2013.

    Book  Google Scholar 

  16. Morton G, Schuster A, Jogiya R, Kutty S, Beerbaum P, Nagel E. Inter-study reproducibility of cardiovascular magnetic resonance myocardial feature tracking. J Cardiovasc Magn Reson. 2012;14:43.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Schuster A, Paul M, Bettencourt N, Morton G, Chiribiri A, Ishida M, Hussain S, Jogiya R, Kutty S, Bigalke B, Perera D, Nagel E. Cardiovascular magnetic resonance myocardial feature tracking for quantitative viability assessment in ischemic cardiomyopathy. Int J Cardiol. 2013;166:413–20.

    Article  PubMed  Google Scholar 

  18. Buss SJ, Breuninger K, Lehrke S, Voss A, Galuschky C, Lossnitzer D, Andre F, Ehlermann P, Franke J, Taeger T, Frankenstein L, Steen H, Meder B, Giannitsis E, Katus HA, Korosoglou G. Assessment of myocardial deformation with cardiac magnetic resonance strain imaging improves risk stratification in patients with dilated cardiomyopathy. Eur Heart J Cardiovasc Imaging. 2015;16:307–15.

    Article  PubMed  Google Scholar 

  19. Schuster A, Morton G, Hussain ST, Jogiya R, Kutty S, Asrress KN, Makowski MR, Bigalke B, Perera D, Beerbaum P, Nagel E. The intra-observer reproducibility of cardiovascular magnetic resonance myocardial feature tracking strain assessment is independent of field strength. Eur J Radiol. 2013;82:296–301.

    Article  PubMed  Google Scholar 

  20. Schuster A, Stahnke V-C, Unterberg-Buchwald C, Kowallick JT, Lamata P, Steinmetz M, Kutty S, Fasshauer M, Staab W, Sohns JM, Bigalke B, Ritter C, Hasenfuss BP, Lotz J. Cardiovascular magnetic resonance feature-tracking assessment of myocardial mechanics: intervendor agreement and considerations regarding reproducibility. Clin Radiology. 2015;70:989–98.

    Article  CAS  Google Scholar 

  21. Singh A, Steadman CD, Khan JN, Horsfield MA, Bekele S, Nazir SA, Kanagala P, Masca NGD, Clarysse P, McCann GP. Intertechnique agreement and interstudy reproducibility of strain and diastolic strain rate at 1.5 and 3 Tesla: a comparison of feature-tracking and tagging in patients with aortic stenosis. J Magn Reson Imaging. 2015;41:1129–37.

    Article  PubMed  Google Scholar 

  22. Augustine D, Lewandowski A, Lazdam M, Rai A, Francis J, Myerson S, Noble A, Becher H, Neubauer S, Petersen SE, Leeson P. Global and regional left ventricular myocardial deformation measures by magnetic resonance feature tracking in healthy volunteers: comparison with tagging and relevance of gender. J Cardiovasc Magn Reson. 2013;15:8.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Cerqueira MD, Weissman NJ, Dilsizian V, Jacobs AK, Kaul S, Laske WK, Pennell DJ, Rumberger JA, Ryan T, Verani MS. American Heart Association Writing Group on Myocardial Segmentation and Registration for Cardiac Imaging. Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart. A statement for healthcare professionals from the Cardiac Imaging Committee of the Council on Clinical Cardiology of the American Heart Association. Circulation. 2002;105:539–42.

    Article  PubMed  Google Scholar 

  24. Vaz S., Falkmer T., Passmore A. E., Parsons R., Andreou P., The Case for Using the Repeatability Coefficent When Calculating Test-Retest Reliability. Plos One. 2013;8(9):e73990.

  25. Tobon-Gomez C, Craene M, McLeod K, Tautz L, Shi W, et al. Benchmarking framework for myocardial tracking and deformation algorithms: An open access database. Med Image Anal. 2013;17:632–48.

    Article  CAS  PubMed  Google Scholar 

  26. Alessandrini M., Heyde B., Queirós S., Cygan S., Zontak M., et al. Detailed Evaluation of Five 3D Speckle Tracking ALgorithms using Synthetic Echocardiographic Recordings. IEEE Transactions on Medical Imaging (in press).

  27. Kowallick JT, Morton G, Lamata P, Jogiya R, Kutty S, Lotz J, Hasenfuss G, Nagel E, Chiribiri A, Schuster A. Inter-study variability of left ventricular torsion and torsion rate quantification using MR myocardial feature tracking. J Magn Reson Imaging. 2016;43:128–37.

    Article  PubMed  Google Scholar 

  28. Schuster A, Kutty S, Padiyath A, Parish V, Gribben P, Danford DA, Makowski MR, Bigalke B, Beerbaum P, Nagel E. Cardiovascular magnetic resonance myocardial feature tracking detects quantitative wall motion during dobutamine stress. J Cardiovasc Magn Reson. 2011;13:58.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Onishi T, Saha SK, Delgado-Montero A, Ludwig DR, Onishi T, Schelbert EB, Schwartzman D, Gorcsan J. Global longitudinal strain and global circumferential strain by speckle-tracking echocardiography and feature-tracking cardiac magnetic resonance imaging: comparison with left ventricular ejection fraction. J Am Soc Echocardiogr. 2015;28:587–96.

    Article  PubMed  Google Scholar 

  30. Aurich M, Keller M, Greiner S, Steen H, aus dem Siepen F, Riffel J, Katus HA, Buss SJ, Mereles D. Left ventricular mechanics assessed by two-dimensional echocardiography and cardiac magnetic resonance imaging: comparison of high-resolution speckle tracking and feature tracking. Eur Heart J Cardiovasc Imaging 2016 Mar 24 (Epub ahead of print).

  31. Kuetting DL, Dabir D, Homsi R, Sprinkart AM, Luetkens J, Schild HH, Thomas DK. The effects of extracellular contrast agent (Gadobutrol) on the precision and reproducibility of cardiovascular magnetic resonance feature tracking. J Cardiovasc Magn Reson. 2016;18(1):30. doi:10.1186/s12968-016-0249-y.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Andre F, Steen H, Matheis P, Westkott M, Breuninger K, Sander Y, Kammerer R, Galuschky C, Giannitis E, Korosoglou G, Katus HA, Buss SJ. Age- and gender-related normal left ventricular deformation assessed by cardiovascular magnetic resonance feature tracking. J Cardiovasc Magn Reson. 2015;17:25.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Taylor RJ, Moody WE, Umar F, Edwards NC, Taylor TJ, Stegemann B, Townend JN, Hor KN, Steeds RP, Mazur W, Leyva F. Myocardial strain measurement with feature-tracking cardiovascular magnetic resonance: normal values. Eur Heart J Cardiovasc Imaging. 2015;16:871–81.

    Article  PubMed  Google Scholar 

  34. Schuster A, Paul M, Bettencourt N, Hussain ST, Morton G, Kutty S, Bigalke B, Chiribiri A, Perera D, Nagel E, Beerbaum P. Myocardial feature tracking reduces observer-dependence in low-dose dobutamine stress cardiovascular magnetic resonance. PLoS One. 2015;10:e0122858.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


Not applicable


Pedro Morias received funding for his PhD scholarship (FCT - Fundacão para a Ciência e a Tecnologia, Portugal, for funding support through the Programa Operacional Capital Humano (POCH) in the scope of the PhD grant SFRH/BD/95438/2013).

Availability of data and materials

Supporting data with regard the healthy volunteer, normal and patient group is added in the form of a word document (Additional file 1. Supporting data is available with regard to CMR-FT applied to the right ventricle and left atrium.

Author’s contributions

PM concept and design, interpretation of data, drafting of manuscript, final approval; AM concept and design, collection of data, drafting of the manuscript, final approval. JAB collection of data, final approval. TD critical revision of manuscript, interpretation of data, final approval. BH interpretation of data, critical revision of manuscript, final approval. JD interpretation of data, critical revision of manuscript, final approval. JB concept and design, collection and interpretation of data, drafting the manuscript, critical revision of the manuscript, final approval, corresponding author.

Author’s information

Not applicable.

Competing interests

As designer of the CMR-FT-algorithm PM, BH and JD receive royalties from Medviso.

Consent for publication

Written informed consent was obtained from the healthy volunteer and the subjects in the patient group for publication of individual details and accompanying in this manuscript. In order not to jeopardize the anonymity of the subjects in the normal and patient group, we removed the age and sex in the online additional file. The consent form is held by the authors and is available by the Editor-in-Chief.

Ethics approval and consent to participate

The study has been approved by the Ethics Committee of the UZ Leuven (S51789). The healthy volunteer and all subjects of the patient group consented in writing to their participation. Because of the retrospective nature of the patient inclusion in the normal group, informed consent in this group was waived.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Jan Bogaert.

Additional files

Additional file 1:

Shows the demographics of the healthy volunteer, and the subjects in the normal and patient group. (DOCX 17 kb)

Additional file 4:

A movie showing feature tracking applied to the left atrium (horizontal long-axis SSFP cine CMR) in a healthy volunteer. (AVI 2111 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Morais, P., Marchi, A., Bogaert, J.A. et al. Cardiovascular magnetic resonance myocardial feature tracking using a non-rigid, elastic image registration algorithm: assessment of variability in a real-life clinical setting. J Cardiovasc Magn Reson 19, 24 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: