Analytical method to measure three-dimensional strain patterns in the left ventricle from single slice displacement data

Background: Displacement encoded Cardiovascular MR (CMR) can provide high spatial resolution measurements of three-dimensional (3D) Lagrangian displacement. Spatial gradients of the Lagrangian displacement field are used to measure regional myocardial strain. In general, adjacent parallel slices are needed in order to calculate the spatial gradient in the through-slice direction. This necessitates the acquisition of additional data and prolongs the scan time. The goal of this study is to define an analytic solution that supports the reconstruction of the out-of-plane components of the Lagrangian strain tensor in addition to the in-plane components from a single-slice displacement CMR dataset with high spatio-temporal resolution. The technique assumes incompressibility of the myocardium as a physical constraint. Results: The feasibility of the method is demonstrated in a healthy human subject and the results are compared to those of other studies. The proposed method was validated with simulated data and strain estimates from experimentally measured DENSE data, which were compared to the strain calculation from a conventional two-slice acquisition. Conclusion: This analytical method reduces the need to acquire data from adjacent slices when calculating regional Lagrangian strains and can effectively reduce the long scan time by a factor of two

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

MR imaging of left-ventricular function : novel image acquisition and analysis techniques.

Many cardiac diseases, such as myocardial ischemia, secondary to coronary artery disease, may be identified and localized through the analysis of cardiac deformations. Early efforts for quantifying ventricular wall motion used surgical implantation and tracking of radiopaque markers with X-ray imaging in canine hearts [1]. Such techniques are invasive and affect the regional motion pattern of the ventricular wall during the marker tracking process and, clearly are not feasible clinically. Noninvasive imaging techniques are vital and have been widely applied to the clinic. MRI is a noninvasive imaging technique with the capability to monitor and assess the progression of cardiovascular diseases (CVD) so that effective procedures for the care and treatment of patients can be developed by physicians and researchers. It is 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 imaging sequences. In order to assess cardiac function, there are two categories of indices that are used: global and regional indices. Global indices include ejection fraction, cavity volume, and myocardial mass [2]. They are important indices for cardiac disease diagnosis. However, these global indices are not specific for regional analysis. A quantitative assessment of regional parameters may prove beneficial for the diagnosis of disease and evaluation of severity and the quantification of treatment [3]. Local measures, such as wall deformation and strain in all regions of the heart, can provide objective regional quantification of ventricular wall function and relate to the location and extent of ischemic injury. This dissertation is concerned with the development of novel MR imaging techniques and image postprocessing algorithms to analyze left ventricular deformations. A novel pulse sequence, termed Orthogonal CSPAMM (OCSPAMM), has been proposed which results in the same acquisition time as SPAMM for 2D deformation estimation while keeping the main advantages of CSPAMM [4,5]: i.e., maintaining tag contrast through-out the ECG cycle. Different from CSPAMM, in OCSPAMM the second tagging pulse orientation is rotated 90 degrees relative to the first one so that motion information can be obtained simultaneously in two directions. This reduces the acquisition time by a factor of two as compared to the traditional CSPAMM, in which two separate imaging sequences are applied per acquisition. With the application of OCSPAMM, the effect of tag fading encountered in SPAMM tagging due to Tl relaxation is mitigated and tag deformations can be visualized for the entire cardiac cycle, including diastolic phases. A multilevel B-spline fitting method (MBS) has been proposed which incorporates phase-based displacement information for accurate calculation of 2D motion and strain from tagged MRI [6, 7]. The proposed method combines the advantages of continuity and smoothness of MBS, and makes use of phase information derived from tagged MR images. Compared to previous 2D B-spline-based deformation analysis methods, MBS has the following advantages: 1) It can simultaneously achieve a smooth deformation while accurately approximating the given data set; 2) Computationally, it is very fast; and 3) It can produce more accurate deformation results. Since the tag intersections (intersections between two tag lines) can be extracted accurately and are more or less distributed evenly over the myocardium, MBS has proven effective for 2D cardiac motion tracking. To derive phase-based displacements, 2D HARP and SinMod analysis techniques [8,9] were employed. By producing virtual tags from HARP /SinMod and calculating intersections of virtual tag lines, more data points are obtained. In the reference frame, virtual tag lines are the isoparametric curves of an undeformed 2D B-spline model. In subsequent frames, the locations of intersections of virtual tag lines over the myocardium are updated with phase-based displacement. The advantage of the technique is that in acquiring denser myocardial displacements, it uses both real and virtual tag line intersections. It is fast and more accurate than 2D HARP and SinMod tracking. A novel 3D sine wave modeling (3D SinMod) approach for automatic analysis of 3D cardiac deformations has been proposed [10]. An accelerated 3D complementary spatial modulation of magnetization (CSPAMM) tagging technique [11] was used to acquire complete 3D+t tagged MR data sets of the whole heart (3 dynamic CSPAMM tagged MRI volume with tags in different orientations), in-vivo, in 54 heart beats and within 3 breath-holds. In 3D SinMod, the intensity distribution around each pixel is modeled as a cosine wave front. The principle behind 3D SinMod tracking is that both phase and frequency for each voxel are determined directly from the frequency analysis and the displacement is calculated from the quotient of phase difference and local frequency. The deformation fields clearly demonstrate longitudinal shortening during systole. The contraction of the LV base towards the apex as well as the torsional motion between basal and apical slices is clearly observable from the displacements. 3D SinMod can automatically process the image data to derive measures of motion, deformations, and strains between consecutive pair of tagged volumes in 17 seconds. Therefore, comprehensive 4D imaging and postprocessing for determination of ventricular function is now possible in under 10 minutes. For validation of 3D SinMod, 7 3D+t CSPAMM data sets of healthy subjects have been processed. Comparison of mid-wall contour deformations and circumferential shortening results by 3D SinMod showed good agreement with those by 3D HARP. Tag lines tracked by the proposed technique were also compared with manually delineated ones. The average errors calculated for the systolic phase of the cardiac cycles were in the sub-pixel range

Myocardial strain analysis with high temporal resolution MRI tagging: extended 3D motion tracking in normal and LBBB hearts

Tese de doutoramento em Biofísica, apresentada à Universidade de Lisboa através da Faculdade de Ciências, 200

Comprehensive Cardiovascular magnetic resonance of myocardial mechanics in mice using three-dimensional cine DENSE

<p>Abstract</p> <p>Background</p> <p>Quantitative noninvasive imaging of myocardial mechanics in mice enables studies of the roles of individual genes in cardiac function. We sought to develop comprehensive three-dimensional methods for imaging myocardial mechanics in mice.</p> <p>Methods</p> <p>A 3D cine DENSE pulse sequence was implemented on a 7T small-bore scanner. The sequence used three-point phase cycling for artifact suppression and a stack-of-spirals <it>k</it>-space trajectory for efficient data acquisition. A semi-automatic 2D method was adapted for 3D image segmentation, and automated 3D methods to calculate strain, twist, and torsion were employed. A scan protocol that covered the majority of the left ventricle in a scan time of less than 25 minutes was developed, and seven healthy C57Bl/6 mice were studied.</p> <p>Results</p> <p>Using these methods, multiphase normal and shear strains were measured, as were myocardial twist and torsion. Peak end-systolic values for the normal strains at the mid-ventricular level were 0.29 ± 0.17, -0.13 ± 0.03, and -0.18 ± 0.14 for <it>E_rr</it>, <it>E_cc</it>, and <it>E_ll</it>, respectively. Peak end-systolic values for the shear strains were 0.00 ± 0.08, 0.04 ± 0.12, and 0.03 ± 0.07 for <it>E_rc</it>, <it>E_rl</it>, and <it>E_cl</it>, respectively. The peak end-systolic normalized torsion was 5.6 ± 0.9°.</p> <p>Conclusions</p> <p>Using a 3D cine DENSE sequence tailored for cardiac imaging in mice at 7 T, a comprehensive assessment of 3D myocardial mechanics can be achieved with a scan time of less than 25 minutes and an image analysis time of approximately 1 hour.</p

VALIDATION, OPTIMIZATION, AND IMAGE PROCESSING OF SPIRAL CINE DENSE MAGNETIC RESONANCE IMAGING FOR THE QUANTIFICATION OF LEFT AND RIGHT VENTRICULAR MECHANICS

Recent evidence suggests that cardiac mechanics (e.g. cardiac strains) are better measures of heart function compared to common clinical metrics like ejection fraction. However, commonly-used parameters of cardiac mechanics remain limited to just a few measurements averaged over the whole left ventricle. We hypothesized that recent advances in cardiac magnetic resonance imaging (MRI) could be extended to provide measures of cardiac mechanics throughout the left and right ventricles (LV and RV, respectively). Displacement Encoding with Stimulated Echoes (DENSE) is a cardiac MRI technique that has been validated for measuring LV mechanics at a magnetic field strength of 1.5 T but not at higher field strengths such as 3.0 T. However, it is desirable to perform DENSE at 3.0 T, which would yield a better signal to noise ratio for imaging the thin RV wall. Results in Chapter 2 support the hypothesis that DENSE has similar accuracy at 1.5 and 3.0 T. Compared to standard, clinical cardiac MRI, DENSE requires more expertise to perform and is not as widely used. If accurate mechanics could be measured from standard MRI, the need for DENSE would be reduced. However, results from Chapter 3 support the hypothesis that measured cardiac mechanics from standard MRI do not agree with, and thus cannot be used in place of, measurements from DENSE. Imaging the thin RV wall with its complex contraction pattern requires both three-dimensional (3D) measures of myocardial motion and higher resolution imaging. Results from Chapter 4 support the hypothesis that a lower displacement-encoding frequency can be used to allow for easier processing of 3D DENSE images. Results from Chapter 5 support the hypothesis that images with higher resolution (decreased blurring) can be achieved by using more spiral interleaves during the DENSE image acquisition. Finally, processing DENSE images to yield measures of cardiac mechanics in the LV is relatively simple due to the LV’s mostly cylindrical geometry. Results from Chapter 6 support the hypothesis that a local coordinate system can be adapted to the geometry of the RV to quantify mechanics in an equivalent manner as the LV. In summary, cardiac mechanics can now be quantified throughout the left and right ventricles using DENSE cardiac MRI

Analysis of myocardial contractility with magnetic resonance

Heart failure has considerable morbidity and poor prognosis. An understanding of the underlying mechanics governing myocardial contraction is a prerequisite for interpreting and predicting changes induced by heart disease. Gross changes in contractile behaviour of the myocardium are readily detected with existing techniques. For more subtle changes during early stages of cardiac dysfunction, however, it requires a sensitive method for measuring, as well as a precise criterion for quantifying, normal and impaired myocardial function. Cardiovascular Magnetic Resonance (CMR) imaging is emerging as an important clinical tool because of its safety, versatility, and the high quality images it produces that allow accurate and reproducible quantification of cardiac structure and function. Traditional CMR approaches for measuring contractility rely on tagging of the myocardium with fiducial markers and require a lengthy and often subjective dependant post-processing procedure. The aim of this research is to develop a new technique, which uses velocity as a marker for the visualisation and assessment of myocardial contractility. Two parallel approaches have been investigated for the assessment of myocardial velocity. The first of these is haimonic phase (HARP) imaging. HARP imaging allows direct derivation of myocardial velocity and strain without the need of further user interaction. We investigated the effect of respiration on the accuracy of the derived contractility, and assessed the clinical applicability and potential pitfalls of the technique by analysing results from a group of patients with hypertrophic cardiomyopathy. The second technique we have investigated is the direct measurement of myocardial velocity with phase contrast myocardial velocity mapping. The imaging sequence used employs effective blood saturation for reducing flow induced phase errors within the myocardium. View sharing was used to improve the temporal resolution, which permitted acquisition of 3D velocity information throughout the cardiac cycle in a single breath-hold, enabling a comprehensive assessment of strain rate of the left ventricle. One key factor that affects the derivation of myocardial contractility based on myocardial velocity is the practical inconsistency of the velocity data. A novel iterative optimisation scheme by incorporating the incompressibility constraint was developed for the restoration of myocardial velocity data. The method allowed accurate assessment of both in-plane and through-plan strain rates, as demonstrated with both synthetic and in vivo data acquired from normal subjects and ischaemic patients. To further enhance the clinical potential of the technique and facilitate the visual assessment of contractile abnormality with myocardial velocity mapping, a complementary analysis framework, named Virtual Tagging, has been developed. The method used velocity data in all directions combined with a finite element mesh incorporating geometrical and physical constraints. The Virtual Tagging framewoik allowed velocity measurements to be used for calculating strain distribution within the 3D volume. It also permitted easy visualisation of the displacement of the tissue, akin to traditional CMR tagging. Detailed validation of the technique is provided, which involves both numerical simulation and in vitro phantom experiments. The main contribution of this thesis is in the improvement of the effectiveness and quality of quantitative myocardial contractility analysis from both sequence design and medical image computing perspectives. It is aimed at providing a sensitive means of detecting subtle as well as gross changes in contractile behaviour of the myocardium. The study is expected to provide a clinically viable platform for functional correlation with other functional measures such as myocardial perfusion and diffusion, and to serve as an aid for further understanding of the links between intrinsicOpen acces

QUANTIFICATION OF MYOCARDIAL MECHANICS IN LEFT VENTRICLES UNDER INOTROPIC STIMULATION AND IN HEALTHY RIGHT VENTRICLES USING 3D DENSE CMR

Statistical data from clinical studies indicate that the death rate caused by heart disease has decreased due to an increased use of evidence-based medical therapies. This includes the use of magnetic resonance imaging (MRI), which is one of the most common non-invasive approaches in evidence-based health care research. In the current work, I present 3D Lagrangian strains and torsion in the left ventricle of healthy and isoproterenol-stimulated rats, which were investigated using Displacement ENcoding with Stimulated Echoes (DENSE) cardiac magnetic resonance (CMR) imaging. With the implementation of the 12-segment model, a detailed profile of regional cardiac mechanics was reconstructed for each subject. Statistical analysis revealed that isoproterenol induced a significant change in the strains and torsion in certain regions at the mid-ventricle level. In addition, I investigated right ventricular cardiac mechanics with the methodologies developed for the left ventricle. This included a comparison of different regions within the basal and mid-ventricular regions. Despite no regional variation found in the peak circumferential strain, the peak longitudinal strain exhibited regional variation at the anterior side of the RV due to the differences in biventricular torsion, mechanism of RV free wall contraction, and fiber architecture at RV insertions. Future applications of the experimental work presented here include the construction and validation of biventricular finite element models. Specifically, the strains predicted by the models will be statistically compared with experimental strains. In addition, the results of the present study provide an essential reference of RV baseline evaluated with DENSE MRI, a highly objective technique

Towards automating cine DENSE MRI image analysis : segmentation, tissue tracking and strain computation

Includes bibliographical references (p. 192-206).Over the past two decades, magnetic resonance imaging (MRI) has developed into a powerful imaging tool for the heart. Imaging cardiac morphology is now commonplace in clinical practice, and a plethora of quantitative techniques have also arisen on the research front. Myocardial tagging is an established quantitative cardiac MRI method that involves magnetically tagging the heart with a set of saturated bands, and monitoring the deformation of these bands as the heart contracts