Combined tag tracking and strain reconstruction from tagged cardiac MR images without user-defined myocardial contours

Computational Analysis of Cardiac Contractile Function

PURPOSE OF REVIEW

Heart failure results in the high incidence and mortality all over the world. Mechanical properties of myocardium are critical determinants of cardiac function, with regional variations in myocardial contractility demonstrated within infarcted ventricles. Quantitative assessment of cardiac contractile function is therefore critical to identify myocardial infarction for the early diagnosis and therapeutic intervention.

RECENT FINDINGS

Current advancement of cardiac functional assessments is in pace with the development of imaging techniques. The methods tailored to advanced imaging have been widely used in cardiac magnetic resonance, echocardiography, and optical microscopy. In this review, we introduce fundamental concepts and applications of representative methods for each imaging modality used in both fundamental research and clinical investigations. All these methods have been designed or developed to quantify time-dependent 2-dimensional (2D) or 3D cardiac mechanics, holding great potential to unravel global or regional myocardial deformation and contractile function from end-systole to end-diastole. Computational methods to assess cardiac contractile function provide a quantitative insight into the analysis of myocardial mechanics during cardiac development, injury, and remodeling.

Artifacts reduction in strain maps of tagged magnetic resonance imaging using harmonic phase

Tagged Magnetic Resonance Imaging (MRI) is a noninvasive technique for examining myocardial function and deformation. Tagged MRI can also be used in quasi-static MR elastography to acquire strain maps of other biological soft tissues. Harmonic phase (HARP) provides automatic and rapid analysis of tagged MR images for the quantification and visualization of myocardial strain. We propose a new artifact reduction method in strain maps. Image intensity of the DC component is estimated and subtracted from spatial modulation of magnetization (SPAMM) tagged MR images. DC peak interference in harmonic phase extraction is greatly reduced after DC component subtraction. The proposed method is validated using both simulated and MR acquired tagged images. Strain maps are obtained with better accuracy and smoothness after DC component subtraction.

An image-based gait simulation study of tarsal kinematics in women with hallux valgus

BACKGROUND

Although not well understood, foot kinematics are changed with hallux valgus.

OBJECTIVE

The purpose of this study was to examine tarsal kinematics in women with hallux valgus deformity.

DESIGN

A prospective, cross-sectional design was used.

METHODS

Twenty women with (n=10) and without (n=10) deformity participated. Data were acquired with the use of a magnetic resonance scanner. Participants were posed standing to simulate gait, with images reconstructed into virtual bone datasets. Measures taken described foot posture (hallux angle, intermetatarsal angle, arch angle). With the use of additional computer processes, the image sequence was then registered across gait conditions to compute relative tarsal position angles, first-ray angles, and helical axis parameters decomposed into X, Y, and Z components. An analysis of variance model compared kinematics between groups and across conditions. Multiple regression analysis assessed the relationship of arch angle, navicular position, and inclination of the first-ray axis.

RESULTS

Both the hallux and intermetatarsal angles were larger with deformity; arch angle was not different between groups. The calcaneus was everted by ≥6.6 degrees, and the first ray adducted (F=44.17) by ≥9.3 degrees across conditions with deformity. There was an interaction (F=5.06) for the first-ray axis. Follow-up comparisons detected increased inclination of the first-ray axis over middle stance compared with late stance in the group with deformity.

LIMITATIONS

Gait was simulated, kinetics were not measured, and sample size was small.

CONCLUSIONS

There were group differences. Eversion of the calcaneus and adduction of the first ray were increased, and the first-ray axis was inclined 24 degrees over middle stance in women with deformity compared with 6 degrees in control participants. Results may identify risk factors of hallux valgus and inform nonoperative treatment (orthoses, exercise) strategies.

Accurate two-dimensional cardiac strain calculation using adaptive windowed Fourier transform and Gabor wavelet transform

PURPOSE

Cardiac strain calculated from tagged magnetic resonance (MR) images provides clinicians information about abnormalities of heart-wall motion in patients. It is important to develop an accurate method to determine the cardiac strain efficiently. An adaptive windowed harmonic phase (AWHARP) method is proposed for cardiac strain calculation.

MATERIALS AND METHODS

AWHARP is based on adaptive windowed Fourier transform (AWFT) and 2D Gabor wavelet transform (2D-GWT). The AWFT provides a spatially varying representation of the signal spectra, which allows the harmonic phase (HARP) image to be extracted with high accuracy. Instantaneous spatial frequencies are calculated using 2D-GWT, and the widths of the adaptive windows are then determined according to the instantaneous spatial frequencies for multi-resolution analysis of phase extraction. The proposed method was studied using simulated images and patients' MR images. Both single tagged images (SPAMM) and subtracted tagged images (CSPAMM) were generated using our simulation method, and their results calculated using AWHARP and HARP methods were compared. Normal and pathological tagged MR images were also processed to evaluate the performance of our method.

RESULTS

Our experimental results show that the accuracies of phase and strain images calculated using the AWHARP method are higher than that calculated using the HARP method especially for large tag line deformation. The improvement in accuracies can be up to 3.2 strain (E1) and 17.3 calculation from MR images reveals that the cardiac strain in the end-systolic state is significantly reduced for patients with hypertrophic cardiomyopathy (HCM) compared to that of healthy subjects.

CONCLUSION

The proposed AWHARP is an accurate and efficient method for cardiac strain estimation from MR images. This new algorithm can help clinicians to detect left ventricle dysfunctions and myocardial diseases with accurate cardiac strain analysis.

Cardiac motion and deformation recovery from MRI: a review

Magnetic resonance imaging (MRI) is a highly advanced and sophisticated imaging modality for cardiac motion tracking and analysis, capable of providing 3D analysis of global and regional cardiac function with great accuracy and reproducibility. In the past few years, numerous efforts have been devoted to cardiac motion recovery and deformation analysis from MR image sequences. Many approaches have been proposed for tracking cardiac motion and for computing deformation parameters and mechanical properties of the heart from a variety of cardiac MR imaging techniques. In this paper, an updated and critical review of cardiac motion tracking methods including major references and those proposed in the past ten years is provided. The MR imaging and analysis techniques surveyed are based on cine MRI, tagged MRI, phase contrast MRI, DENSE, and SENC. This paper can serve as a tutorial for new researchers entering the field.

Motion tracking and strain map computation for quasi-static magnetic resonance elastography

This paper presents a new imaging method for quasi-static magnetic resonance elastography (MRE). Tagged magnetic resonance (MR) imaging of human lower leg was acquired with probe indentation using a MR-compatible actuation system. Indentation force was recorded for soft tissue elasticity reconstruction. Motion tracking and strain map of human lower leg are calculated using a harmonic phase (HARP)-based method. Simulated tagged MR images were constructed and analyzed to validate the HARP-based method. Our results show that the proposed imaging method can be used to generate accurate motion distribution and strain maps of the targeted soft tissue.

Myocardial tagging by cardiovascular magnetic resonance: evolution of techniques--pulse sequences, analysis algorithms, and applications

Cardiovascular magnetic resonance (CMR) tagging has been established as an essential technique for measuring regional myocardial function. It allows quantification of local intramyocardial motion measures, e.g. strain and strain rate. The invention of CMR tagging came in the late eighties, where the technique allowed for the first time for visualizing transmural myocardial movement without having to implant physical markers. This new idea opened the door for a series of developments and improvements that continue up to the present time. Different tagging techniques are currently available that are more extensive, improved, and sophisticated than they were twenty years ago. Each of these techniques has different versions for improved resolution, signal-to-noise ratio (SNR), scan time, anatomical coverage, three-dimensional capability, and image quality. The tagging techniques covered in this article can be broadly divided into two main categories: 1) Basic techniques, which include magnetization saturation, spatial modulation of magnetization (SPAMM), delay alternating with nutations for tailored excitation (DANTE), and complementary SPAMM (CSPAMM); and 2) Advanced techniques, which include harmonic phase (HARP), displacement encoding with stimulated echoes (DENSE), and strain encoding (SENC). Although most of these techniques were developed by separate groups and evolved from different backgrounds, they are in fact closely related to each other, and they can be interpreted from more than one perspective. Some of these techniques even followed parallel paths of developments, as illustrated in the article. As each technique has its own advantages, some efforts have been made to combine different techniques together for improved image quality or composite information acquisition. In this review, different developments in pulse sequences and related image processing techniques are described along with the necessities that led to their invention, which makes this article easy to read and the covered techniques easy to follow. Major studies that applied CMR tagging for studying myocardial mechanics are also summarized. Finally, the current article includes a plethora of ideas and techniques with over 300 references that motivate the reader to think about the future of CMR tagging.

3D left ventricular strain from unwrapped harmonic phase measurements

PURPOSE

To validate a method for measuring 3D left ventricular (LV) strain from phase-unwrapped harmonic phase (HARP) images derived from tagged cardiac magnetic resonance imaging (MRI).

MATERIALS AND METHODS

A set of 40 human subjects were imaged with tagged MRI. In each study the HARP phase was computed and unwrapped in each short-axis and long-axis image. Inconsistencies in unwrapped phase were resolved using branch cuts manually placed with a graphical user interface. 3D strain maps were computed for all imaged timeframes in each study. The strain from unwrapped phase (SUP) and displacements were compared to those estimated by a feature-based (FB) technique and a HARP technique.

RESULTS

3D strain was computed in each timeframe through systole and mid-diastole in approximately 30 minutes per study. The standard deviation of the difference between strains measured by the FB and the SUP methods was less than 5% of the average of the strains from the two methods. The correlation between peak circumferential strain measured using the SUP and HARP techniques was over 83%.

CONCLUSION

The SUP technique can reconstruct full 3D strain maps from tagged MR images through the cardiac cycle in a reasonable amount of time and user interaction compared to other 3D analysis methods.

Deformation analysis of 3D tagged cardiac images using an optical flow method

BACKGROUND

This study proposes and validates a method of measuring 3D strain in myocardium using a 3D Cardiovascular Magnetic Resonance (CMR) tissue-tagging sequence and a 3D optical flow method (OFM).

METHODS

Initially, a 3D tag MR sequence was developed and the parameters of the sequence and 3D OFM were optimized using phantom images with simulated deformation. This method then was validated in-vivo and utilized to quantify normal sheep left ventricular functions.

RESULTS

Optimizing imaging and OFM parameters in the phantom study produced sub-pixel root-mean square error (RMS) between the estimated and known displacements in the x (RMSx = 0.62 pixels (0.43 mm)), y (RMSy = 0.64 pixels (0.45 mm)) and z (RMSz = 0.68 pixels (1 mm)) direction, respectively. In-vivo validation demonstrated excellent correlation between the displacement measured by manually tracking tag intersections and that generated by 3D OFM (R >or= 0.98). Technique performance was maintained even with 20% Gaussian noise added to the phantom images. Furthermore, 3D tracking of 3D cardiac motions resulted in a 51% decrease in in-plane tracking error as compared to 2D tracking. The in-vivo function studies showed that maximum wall thickening was greatest in the lateral wall, and increased from both apex and base towards the mid-ventricular region. Regional deformation patterns are in agreement with previous studies on LV function.

CONCLUSION

A novel method was developed to measure 3D LV wall deformation rapidly with high in-plane and through-plane resolution from one 3D cine acquisition.

A normalized framework for the design of feature spaces assessing the left ventricular function

A through description of the left ventricle functionality requires combining complementary regional scores. A main limitation is the lack of multiparametric normality models oriented to the assessment of regional wall motion abnormalities (RWMA). This paper covers two main topics involved in RWMA assessment. We propose a general framework allowing the fusion and comparison across subjects of different regional scores. Our framework is used to explore which combination of regional scores (including 2-D motion and strains) is better suited for RWMA detection. Our statistical analysis indicates that for a proper (within interobserver variability) identification of RWMA, models should consider motion and extreme strains.

Direct three-dimensional myocardial strain tensor quantification and tracking using zHARP

Images of myocardial strain can be used to diagnose heart disease, plan and monitor treatment, and to learn about cardiac structure and function. Three-dimensional (3D) strain is typically quantified using many magnetic resonance (MR) images obtained in two or three orthogonal planes. Problems with this approach include long scan times, image misregistration, and through-plane motion. This article presents a novel method for calculating cardiac 3D strain using a stack of two or more images acquired in only one orientation. The zHARP pulse sequence encodes in-plane motion using MR tagging and out-of-plane motion using phase encoding, and has been previously shown to be capable of computing 3D displacement within a single image plane. Here, data from two adjacent image planes are combined to yield a 3D strain tensor at each pixel; stacks of zHARP images can be used to derive stacked arrays of 3D strain tensors without imaging multiple orientations and without numerical interpolation. The performance and accuracy of the method is demonstrated in vitro on a phantom and in vivo in four healthy adult human subjects.

Left ventricular motion reconstruction with a prolate spheroidal B-spline model

Tagged cardiac magnetic resonance (MR) imaging can non-invasively image deformation of the left ventricular (LV) wall. Three-dimensional (3D) analysis of tag data requires fitting a deformation model to tag lines in the image data. In this paper, we present a 3D myocardial displacement and strain reconstruction method based on a B-spline deformation model defined in prolate spheroidal coordinates, which more closely matches the shape of the LV wall than existing Cartesian or cylindrical coordinate models. The prolate spheroidal B-spline (PSB) deformation model also enforces smoothness across and can compute strain at the apex. The PSB reconstruction algorithm was evaluated on a previously published data set to allow head-to-head comparison of the PSB model with existing LV deformation reconstruction methods. We conclude that the PSB method can accurately reconstruct deformation and strain in the LV wall from tagged MR images and has several advantages relative to existing techniques.