Detailed Evaluation of Five 3D Speckle Tracking Algorithms Using Synthetic Echocardiographic Recordings

Robust unsupervised texture segmentation for motion analysis in ultrasound images

PURPOSE

Ultrasound imaging has emerged as a promising cost-effective and portable non-irradiant modality for the diagnosis and follow-up of diseases. Motion analysis can be performed by segmenting anatomical structures of interest before tracking them over time. However, doing so in a robust way is challenging as ultrasound images often display a low contrast and blurry boundaries.

METHODS

In this paper, a robust descriptor inspired from the fractal dimension is presented to locally characterize the gray-level variations of an image. This descriptor is an adaptive grid pattern whose scale locally varies as the gray-level variations of the image. Robust features are then located based on the gray-level variations, which are more likely to be consistently tracked over time despite the presence of noise.

RESULTS

The method was validated on three datasets: segmentation of the left ventricle on simulated echocardiography (Dice coefficient, DC), accuracy of diaphragm motion tracking for healthy subjects (mean sum of distances, MSD) and for a scoliosis patient (root mean square error, RMSE). Results show that the method segments the left ventricle accurately ( DC = 0.84 ) and robustly tracks the diaphragm motion for healthy subjects ( MSD = 1.10 mm) and for the scoliosis patient ( RMSE = 1.22 mm).

CONCLUSIONS

This method has the potential to segment structures of interest according to their texture in an unsupervised fashion, as well as to help analyze the deformation of tissues. Possible applications are not limited to US image. The same principle could also be applied to other medical imaging modalities such as MRI or CT scans.

Coordinate-Independent 3-D Ultrasound Principal Stretch and Direction Imaging

OBJECTIVE

In clinical ultrasound, current 2-D strain imaging faces challenges in quantifying three orthogonal normal strain components. This requires separate image acquisitions based on the pixel-dependent cardiac coordinate system, leading to additional computations and estimation discrepancies due to probe orientation. Most systems lack shear strain information, as displaying all components is challenging to interpret.

METHODS

This paper presents a 3-D high-spatial-resolution, coordinate-independent strain imaging approach based on principal stretch and axis estimation. All strain components are transformed into three principal stretches along three normal principal axes, enabling direct visualization of the primary deformation. We devised an efficient 3-D speckle tracking method with tilt filtering, incorporating randomized searching in a two-pass tracking framework and rotating the phase of the 3-D correlation function for robust filtering. The proposed speckle tracking approach significantly improves estimates of displacement gradients related to the axial displacement component. Non-axial displacement gradient estimates are enhanced using a correlation-weighted least-squares method constrained by tissue incompressibility.

RESULTS

Simulated and in vivo canine cardiac datasets were evaluated to estimate Lagrangian strains from end-diastole to end-systole. The proposed speckle tracking method improves displacement estimation by a factor of 4.3 to 10.5 over conventional 1-pass processing. Maximum principal axis/direction imaging enables better detection of local disease regions than conventional strain imaging.

CONCLUSION

Coordinate-independent tracking can identify myocardial abnormalities with high accuracy.

SIGNIFICANCE

This study offers enhanced accuracy and robustness in strain imaging compared to current methods.

2-Dimensional Echocardiographic Global Longitudinal Strain With Artificial Intelligence Using Open Data From a UK-Wide Collaborative

BACKGROUND

Global longitudinal strain (GLS) is reported to be more reproducible and prognostic than ejection fraction. Automated, transparent methods may increase trust and uptake.

OBJECTIVES

The authors developed open machine-learning-based GLS methodology and validate it using multiexpert consensus from the Unity UK Echocardiography AI Collaborative.

METHODS

We trained a multi-image neural network (Unity-GLS) to identify annulus, apex, and endocardial curve on 6,819 apical 4-, 2-, and 3-chamber images. The external validation dataset comprised those 3 views from 100 echocardiograms. End-systolic and -diastolic frames were each labelled by 11 experts to form consensus tracings and points. They also ordered the echocardiograms by visual grading of longitudinal function. One expert calculated global strain using 2 proprietary packages.

RESULTS

The median GLS, averaged across the 11 individual experts, was -16.1 (IQR: -19.3 to -12.5). Using each case's expert consensus measurement as the reference standard, individual expert measurements had a median absolute error of 2.00 GLS units. In comparison, the errors of the machine methods were: Unity-GLS 1.3, proprietary A 2.5, proprietary B 2.2. The correlations with the expert consensus values were for individual experts 0.85, Unity-GLS 0.91, proprietary A 0.73, proprietary B 0.79. Using the multiexpert visual ranking as the reference, individual expert strain measurements found a median rank correlation of 0.72, Unity-GLS 0.77, proprietary A 0.70, and proprietary B 0.74.

CONCLUSIONS

Our open-source approach to calculating GLS agrees with experts' consensus as strongly as the individual expert measurements and proprietary machine solutions. The training data, code, and trained networks are freely available online.

Using averaged models from 4D ultrasound strain imaging allows to significantly differentiate local wall strains in calcified regions of abdominal aortic aneurysms

Abdominal aortic aneurysms are a degenerative disease of the aorta associated with high mortality. To date, in vivo information to characterize the individual elastic properties of the aneurysm wall in terms of rupture risk is lacking. We have used time-resolved 3D ultrasound strain imaging to calculate spatially resolved in-plane strain distributions characterized by mean and local maximum strains, as well as indices of local variations in strains. Likewise, we here present a method to generate averaged models from multiple segmentations. Strains were then calculated for single segmentations and averaged models. After registration with aneurysm geometries based on CT-A imaging, local strains were divided into two groups with and without calcifications and compared. Geometry comparison from both imaging modalities showed good agreement with a root mean squared error of 1.22 ± 0.15 mm and Hausdorff Distance of 5.45 ± 1.56 mm (mean ± sd, respectively). Using averaged models, circumferential strains in areas with calcifications were 23.2 ± 11.7% (mean ± sd) smaller and significantly distinguishable at the 5% level from areas without calcifications. For single segmentations, this was possible only in 50% of cases. The areas without calcifications showed greater heterogeneity, larger maximum strains, and smaller strain ratios when computed by use of the averaged models. Using these averaged models, reliable conclusions can be made about the local elastic properties of individual aneurysm (and long-term observations of their change), rather than just group comparisons. This is an important prerequisite for clinical application and provides qualitatively new information about the change of an abdominal aortic aneurysm in the course of disease progression compared to the diameter criterion.

Fast strain mapping in abdominal aortic aneurysm wall reveals heterogeneous patterns

Abdominal aortic aneurysm patients are regularly monitored to assess aneurysm development and risk of rupture. A preventive surgical procedure is recommended when the maximum aortic antero-posterior diameter, periodically assessed on two-dimensional abdominal ultrasound scans, reaches 5.5 mm. Although the maximum diameter criterion has limited ability to predict aneurysm rupture, no clinically relevant tool that could complement the current guidelines has emerged so far. In vivo cyclic strains in the aneurysm wall are related to the wall response to blood pressure pulse, and therefore, they can be linked to wall mechanical properties, which in turn contribute to determining the risk of rupture. This work aimed to enable biomechanical estimations in the aneurysm wall by providing a fast and semi-automatic method to post-process dynamic clinical ultrasound sequences and by mapping the cross-sectional strains on the B-mode image. Specifically, the Sparse Demons algorithm was employed to track the wall motion throughout multiple cardiac cycles. Then, the cyclic strains were mapped by means of radial basis function interpolation and differentiation. We applied our method to two-dimensional sequences from eight patients. The automatic part of the analysis took under 1.5 min per cardiac cycle. The tracking method was validated against simulated ultrasound sequences, and a maximum root mean square error of 0.22 mm was found. The strain was calculated both with our method and with the established finite-element method, and a very good agreement was found, with mean differences of one order of magnitude smaller than the image spatial resolution. Most patients exhibited a strain pattern that suggests interaction with the spine. To conclude, our method is a promising tool for investigating abdominal aortic aneurysm wall biomechanics as it can provide a fast and accurate measurement of the cyclic wall strains from clinical ultrasound sequences.

Imaging the subcellular viscoelastic properties of mouse oocytes

In recent years, cellular biomechanical properties have been investigated as an alternative to morphological assessments for oocyte selection in reproductive science. Despite the high relevance of cell viscoelasticity characterization, the reconstruction of spatially distributed viscoelastic parameter images in such materials remains a major challenge. Here, a framework for mapping viscoelasticity at the subcellular scale is proposed and applied to live mouse oocytes. The strategy relies on the principles of optical microelastography for imaging in combination with the overlapping subzone nonlinear inversion technique for complex-valued shear modulus reconstruction. The three-dimensional nature of the viscoelasticity equations was accommodated by applying an oocyte geometry-based 3D mechanical motion model to the measured wave field. Five domains-nucleolus, nucleus, cytoplasm, perivitelline space, and zona pellucida-could be visually differentiated in both oocyte storage and loss modulus maps, and statistically significant differences were observed between most of these domains in either property reconstruction. The method proposed herein presents excellent potential for biomechanical-based monitoring of oocyte health and complex transformations across lifespan. It also shows appreciable latitude for generalization to cells of arbitrary shape using conventional microscopy equipment.

Estimating Echocardiographic Myocardial Strain of Left Ventricle with Deep Learning

The global longitudinal strain of the myocardial tissue has been shown to be a better indicator of cardiac pathologies in the subclinical stage than other indices, such as the ejection fraction. This article presents a new deep learning approach for strain estimation in 2D echocardiograms. The proposed method improves the performance of the state of the art without losing stability with noisy echocardiograms and achieved an average end point error of 0.14 ± 0.17 pixels in the estimation of the optical flow in the myocardium and an error of 1.34 ± 2.34 % in the estimation of the global longitudinal strain indicator when evaluated in a synthetic echocardiographic dataset. Further research will validate the proposed method by a clinical in-vivo dataset. Clinical relevance- This paper presents a method to estimate the global longitudinal strain index in noisy echocardiograms, which promises to be a better indicator of cardiac pathologies in the subclinical stage than other indices such as the ejection fraction.

Segmentation Enhanced Elastic Image Registration for 2D Speckle Tracking Echocardiography-Performance Study In Silico

Although the two dimensional Speckle Tracking Echocardiography has gained a strong position among medical diagnostic techniques in cardiology, it still requires further developments to improve its repeatability and reliability. Few works have attempted to incorporate the left ventricle segmentation results in the process of displacements and strain estimation to improve its performance. We proposed the use of mask information as an additional penalty in the elastic image registration based displacements estimation. This approach was studied using a short axis view synthetic echocardiographic data, segmented using an active contour method. The obtained masks were distorted to a different degree, using different methods to assess the influence of the segmentation quality on the displacements and strain estimation process. The results of displacements and circumferential strain estimations show, that even though the method is dependent on the mask quality, the potential loss in accuracy due to the poor segmentation quality is much lower than the potential accuracy gain in cases where the segmentation performs well.

Speckle Tracking Accuracy Enhancement by Temporal Super-Resolution of Three-Dimensional Echocardiography Images

Background:

Speckle tracking has always been a challenging issue in echocardiography images due to the lowcontrast and noisy nature of ultrasonic imaging modality. While in ultrasound imaging, framerate is limited by image size and sound speed in tissue, speckle tracking results get worse inthree-dimensional imaging due to its lower frame rate. Therefore, numerous techniques have beenreported to overcome this limitation and enhance tracking accuracy.

Methods:

In this work, we have proposedto increase the frame rate temporally for a sequence of three-dimensional (3D) echocardiographyframes to make tracking more accurate. To increase the number of frames, cubic B-spline is usedto interpolate between intensity variation time curves extracted from every single voxel in theimage during the cardiac cycle. We have shown that the frame rate increase will result in trackingaccuracy improvement.

Results:

To prove the efficiency of the proposed method, numerical evaluation metricsfor tracking are reported to make a comparison between high temporal resolution sequences andlow temporal resolution sequences. Anatomical affine optical flow is selected as the state-of-the-artspeckle tracking method, and a 3D echocardiography dataset is used to evaluate the proposedmethod.

Conclusion:

Results show that it is beneficial for speckle tracking to perform on temporally condensedframes rather than ordinary clinical 3D echocardiography images. Normalized mean enhancementvalues for mean absolute error, Hausdorff distance, and Dice index for all cases and all frames are0.44 ± 0.09, 0.42± 0.09, and 0.36 ± 0.06, respectively.

Non-rigid image registration using a modified fuzzy feature-based inference system for 3D cardiac motion estimation

BACKGROUND AND OBJECTIVE

Non-rigid image registration is a well-established method for estimating cardiac motion on 3D echocardiographic images. However, such images have relatively poor spatio-temporal resolution making registration challenging. Some of the main challenges are extracting features relevant to the registration problem and defining a suitable geometrical transformation to be applied. The latter can be tackled using a fuzzy inference system considering its potential in transformation modeling. From this point of view, feature-based image registration can be considered an identification problem in which the transformation parameters are computed through an optimization process. This study, thus, aims to estimate cardiac motion on 3D echocardiographic images based on feature-based non-rigid image registration through sets of modified fuzzy rules.

METHODS

The 3D volume features were extracted with the popular scale-invariant feature transform (SIFT) descriptors in 3D space. Sets of fuzzy rules were generated according to the extracted features to register every two consecutive frames. Finally, some supplementary rules modified the registration rule for estimating cardiac motion.

RESULTS

Applying the fuzzy feature-based inference system on the STRAUS synthetic database showed the proposed method to be competitive with other well-established registration algorithms in terms of tracking error and accuracy of strain estimates. The proposed algorithm yielded a tracking error of 1 mm and a relative circumferential strain error of 0.82±4.69%. In addition, the potential of the proposed algorithm for clinical applications was confirmed by evaluating its performance on an in-vivo database called CETUS.

CONCLUSION

This paper proposes a novel registration method based on fuzzy logic which was shown to enable tracking complex cardiac deformations in 3D echocardiographic images with high accuracy.

Temporal Super-resolution of Ultrasound Imaging Using Matrix Completion

The temporal super-resolution of the dynamic ultrasound imaging, a means to observe rapid heart movements, is considered an important subject in medical diagnosis of cardiac conditions. Here, a new technique based on the acquisition scheme using the matrix completion (MC) theory is offered for the temporal super-resolution of the two-dimensional (2D) and three-dimensional (3D) ultrasound imaging. MC mentions the problem of completing a low-rank matrix when only a subset of its elements can be observed. Here, the lower scan lines are acquired. Whereby, the proposed method uses temporal and spatial information of the radio frequency (RF) image sequences for the reconstruction of skipped RF lines. This is performed using the construction of the MC images and then reconstruction of them by the MC theory. The results of the proposed method are compared with the compressive sensing (CS) reconstruction methods. The qualitative and quantitative evaluations of 2D and 3D data demonstrate that in the proposed method, which uses the spatial and temporal relation of RF images and the MC theory, the reconstruction is more accurate, and the reconstruction error is lower. The computational complexity of this method is very low. It also does not require hardware adjustments. Therefore, it can be easily implemented in current ultrasound-imaging devices with the frame-rate enhancement. For instance, the frame rate up to two times the original sequence is feasible using the proposed methods, while root mean square error is decreased by about 35% and 30% for 2D and 3D data, respectively, compared with the CS reconstruction method.

Variability in the Assessment of Myocardial Strain Patterns: Implications for Adequate Interpretation

Variability in global and regional peak strain has been thoroughly studied, but variability in spatiotemporal myocardial strain patterns has not been studied as well. This study reports on such variability and its implications for adequate disease interpretation. Forty in-training operators, distributed on 20 workstations, analyzed six cases with representative deformation patterns with commercial speckle tracking. Inter-operator differences were quantified through the variability in myocardial delineations, spatiotemporal longitudinal strain patterns and peak longitudinal strain. Intra-operator differences were assessed similarly using 10 repeated measurements from a single clinician expert. Delineations varied mainly along the lateral wall and at the valve level. Peak longitudinal strain variability was low to moderate. The spatiotemporal strain patterns were consistent despite high variability at the apex and near the valve. The results indicate that relevant pattern assessment is possible despite heterogeneous experience with speckle tracking and that careful interpretation of pattern abnormalities should be recommended before a more systematic quantitative analysis.

Machine Learning Approaches for Myocardial Motion and Deformation Analysis

Information about myocardial motion and deformation is key to differentiate normal and abnormal conditions. With the advent of approaches relying on data rather than pre-conceived models, machine learning could either improve the robustness of motion quantification or reveal patterns of motion and deformation (rather than single parameters) that differentiate pathologies. We review machine learning strategies for extracting motion-related descriptors and analyzing such features among populations, keeping in mind constraints specific to the cardiac application.

Intervendor Discordance of Fetal and Neonatal Myocardial Tissue Doppler and Speckle-Tracking Measurements

BACKGROUND

Fetal and neonatal studies report a wide range of cardiac parameters derived by pulsed-wave Doppler tissue imaging (DTI) and two-dimensional speckle-tracking echocardiographic (STE) imaging. The use of different ultrasound systems and their vendor-specific software compromises the ability to compare echocardiographic findings among various studies. The aim of this study was to evaluate intervendor reproducibility as well as intra- and interobserver repeatability of DTI and STE measurements in normal-term fetuses and neonates.

METHODS

A prospective study was conducted of term fetuses (n = 196) from uncomplicated pregnancies assessed days before the onset of labor and a few hours after birth. Fetal and neonatal DTI and STE parameters were obtained and analyzed using vendor-specific software on three ultrasound systems: Toshiba Aplio MX versus GE Vivid E9 and GE Vivid E9 versus Philips EPIQ. A reproducibility study in fetuses and neonates (n = 118) was performed by systematic scanning with head-to-head comparison.

RESULTS

DTI reproducibility showed moderate to good correlation, with good agreement for fetuses and neonates on Toshiba versus GE (intraclass correlation coefficient [ICC] = 0.4-0.8). Correlation of DTI measurements on GE versus Philips was poor to moderate for fetuses (ICC = 0.1-0.6) and moderate to good for neonates (ICC = 0.5-0.8), with wider limits of agreement. Fetal and neonatal STE parameters revealed very poor correlation (ICC = 0.1-0.3) and agreement among ultrasound vendors. Intra- and interobserver repeatability demonstrated good to excellent correlation of all fetal and neonatal DTI and STE measurements, with good agreement irrespective of the ultrasound platform used.

CONCLUSIONS

These findings demonstrate reliable assessment of fetal and neonatal DTI and STE measurements when performed on the same ultrasound platform, whereas ultrasound machines and software from different vendors give significantly divergent estimates of DTI and STE parameters in fetuses and neonates. These intervendor discrepancies have significant clinical and research implications and should be considered when interpreting and comparing study findings, establishing reference standards, or performing systematic reviews.

MITT: Medical Image Tracking Toolbox

Over the years, medical image tracking has gained considerable attention from both medical and research communities due to its widespread utility in a multitude of clinical applications, from functional assessment during diagnosis and therapy planning to structure tracking or image fusion during image-guided interventions. Despite the ever-increasing number of image tracking methods available, most still consist of independent implementations with specific target applications, lacking the versatility to deal with distinct end-goals without the need for methodological tailoring and/or exhaustive tuning of numerous parameters. With this in mind, we have developed the medical image tracking toolbox (MITT)-a software package designed to ease customization of image tracking solutions in the medical field. While its workflow principles make it suitable to work with 2-D or 3-D image sequences, its modules offer versatility to set up computationally efficient tracking solutions, even for users with limited programming skills. MITT is implemented in both C/C++ and MATLAB, including several variants of an object-based image tracking algorithm and allowing to track multiple types of objects (i.e., contours, multi-contours, surfaces, and multi-surfaces) with several customization features. In this paper, the toolbox is presented, its features discussed, and illustrative examples of its usage in the cardiology field provided, demonstrating its versatility, simplicity, and time efficiency.

Quantitative Assessment of the Effect of the Out-of-Plane Movement of the Homogenous Ellipsoidal Model of the Left Ventricle on the Deformation Measures Estimated Using 2-D Speckle Tracking-An In-Silico Study

Effect of the out-of-plane (OOP) movement amplitude on estimates of global displacements (radial, circumferential) and strains (radial , circumferential ) was studied in an ellipsoidal model of the left ventricle using finite-element modeling (FEM), synthetic ultrasonic data, and short-axis view. This effect was assessed using median of the absolute relative error (RE) of the global parameters. FEM provided node displacements for synthetic ultrasonic data and reference data generation. Displacements were estimated using block-matching (BM) and B-spline (BS) methods. FEM-derived data analysis, free from errors resulting from speckle tracking, indicated that the tissue motion introduced REs of global strain estimates below 4.5%. The effect of the OOP motion amplitude on strain estimates was strain specific and depended on the displacement estimation method. In the case of , the increase of the OOP amplitude resulted in quasi-linear increase of the RE from approximately 10% to 15%. The modulus of the end-systolic (ES) errors of the estimates almost linearly increased with increasing OOP amplitude approximately from 10% to 16%. REs of the estimate were close to 80% and 40%, respectively, in the case of the BM and BS methods, and increased with increasing OOP amplitude. The modulus of the ES errors of the estimates in the case of the BS method was about -40% and showed low sensitivity to the OOP amplitude; in the BM case, these errors varied approximately from -70% to -58% for OOP amplitude from 0 to 15 mm.

Realistic Vendor-Specific Synthetic Ultrasound Data for Quality Assurance of 2-D Speckle Tracking Echocardiography: Simulation Pipeline and Open Access Database

Two-dimensional (2-D) echocardiography is the modality of choice in the clinic for the diagnosis of cardiac disease. Hereto, speckle tracking (ST) packages complement visual assessment by the cardiologist by providing quantitative diagnostic markers of global and regional cardiac function (e.g., displacement, strain, and strain-rate). Yet, the reported high vendor-dependence between the outputs of different ST packages raises clinical concern and hampers the widespread dissemination of the ST technology. In part, this is due to the lack of a solid commonly accepted quality assurance pipeline for ST packages. Recently, we have developed a framework to benchmark ST algorithms for 3-D echocardiography by using realistic simulated volumetric echocardiographic recordings. Yet, 3-D echocardiography remains an emerging technology, whereas the compelling clinical concern is, so far, directed to the standardization of 2-D ST only. Therefore, by building upon our previous work, we present in this paper a pipeline to generate realistic synthetic sequences for 2-D ST algorithms. Hereto, the synthetic cardiac motion is obtained from a complex electromechanical heart model, whereas realistic vendor-specific texture is obtained by sampling a real clinical ultrasound recording. By modifying the parameters in our pipeline, we generated an open-access library of 105 synthetic sequences encompassing: 1) healthy and ischemic motion patterns; 2) the most common apical probe orientations; and 3) vendor-specific image quality from seven different systems. Ground truth deformation is also provided to allow performance analysis. The application of the provided data set is also demonstrated in the benchmarking of a recent academic ST algorithm.

A Framework for the Generation of Realistic Synthetic Cardiac Ultrasound and Magnetic Resonance Imaging Sequences From the Same Virtual Patients

The use of synthetic sequences is one of the most promising tools for advanced in silico evaluation of the quantification of cardiac deformation and strain through 3-D ultrasound (US) and magnetic resonance (MR) imaging. In this paper, we propose the first simulation framework which allows the generation of realistic 3-D synthetic cardiac US and MR (both cine and tagging) image sequences from the same virtual patient. A state-of-the-art electromechanical (E/M) model was exploited for simulating groundtruth cardiac motion fields ranging from healthy to various pathological cases, including both ventricular dyssynchrony and myocardial ischemia. The E/M groundtruth along with template MR/US images and physical simulators were combined in a unified framework for generating synthetic data. We efficiently merged several warping strategies to keep the full control of myocardial deformations while preserving realistic image texture. In total, we generated 18 virtual patients, each with synthetic 3-D US, cine MR, and tagged MR sequences. The simulated images were evaluated both qualitatively by showing realistic textures and quantitatively by observing myocardial intensity distributions similar to real data. In particular, the US simulation showed a smoother myocardium/background interface than the state-of-the-art. We also assessed the mechanical properties. The pathological subjects were discriminated from the healthy ones by both global indexes (ejection fraction and the global circumferential strain) and regional strain curves. The synthetic database is comprehensive in terms of both pathology and modality, and has a level of realism sufficient for validation purposes. All the 90 sequences are made publicly available to the research community via an open-access database.

Efficient Two-Pass 3-D Speckle Tracking for Ultrasound Imaging

Speckle tracking based on block matching is the most common method for multi-dimensional motion estimation in ultrasound elasticity imaging. Extension of two-dimensional (2-D) methods to three dimensions (3-D) has been problematic because of the large computational load of 3-D tracking, as well as performance issues related to the low frame (volume) rates of 3-D images. To address both of these problems, we have developed an efficient two-pass tracking method suited to cardiac elasticity imaging. PatchMatch, originally developed for image editing, has been adapted for ultrasound to provide first-pass displacement estimates. Second-pass estimation uses conventional block matching within a much smaller search region. 3-D displacements are then obtained using correlation filtering previously shown to be effective against speckle decorrelation. Both simulated and in vivo canine cardiac results demonstrate that the proposed two-pass method reduces computational cost compared to conventional 3-D exhaustive search by a factor of 10. Moreover, it outperforms one-pass tracking by a factor of about 3 in terms of root-mean-square error relative to available ground-truth displacements.

Motion Estimation in Echocardiography Using Sparse Representation and Dictionary Learning

This paper introduces a new method for cardiac motion estimation in 2-D ultrasound images. The motion estimation problem is formulated as an energy minimization, whose data fidelity term is built using the assumption that the images are corrupted by multiplicative Rayleigh noise. In addition to a classical spatial smoothness constraint, the proposed method exploits the sparse properties of the cardiac motion to regularize the solution via an appropriate dictionary learning step. The proposed method is evaluated on one data set with available ground-truth, including four sequences of highly realistic simulations. The approach is also validated on both healthy and pathological sequences of in vivo data. We evaluate the method in terms of motion estimation accuracy and strain errors and compare the performance with state-of-the-art algorithms. The results show that the proposed method gives competitive results for the considered data. Furthermore, the in vivo strain analysis demonstrates that meaningful clinical interpretation can be obtained from the estimated motion vectors.This paper introduces a new method for cardiac motion estimation in 2-D ultrasound images. The motion estimation problem is formulated as an energy minimization, whose data fidelity term is built using the assumption that the images are corrupted by multiplicative Rayleigh noise. In addition to a classical spatial smoothness constraint, the proposed method exploits the sparse properties of the cardiac motion to regularize the solution via an appropriate dictionary learning step. The proposed method is evaluated on one data set with available ground-truth, including four sequences of highly realistic simulations. The approach is also validated on both healthy and pathological sequences of in vivo data. We evaluate the method in terms of motion estimation accuracy and strain errors and compare the performance with state-of-the-art algorithms. The results show that the proposed method gives competitive results for the considered data. Furthermore, the in vivo strain analysis demonstrates that meaningful clinical interpretation can be obtained from the estimated motion vectors.

Fully automatic left ventricular myocardial strain estimation in 2D short-axis tagged magnetic resonance imaging

Cardiovascular diseases are among the leading causes of death and frequently result in local myocardial dysfunction. Among the numerous imaging modalities available to detect these dysfunctional regions, cardiac deformation imaging through tagged magnetic resonance imaging (t-MRI) has been an attractive approach. Nevertheless, fully automatic analysis of these data sets is still challenging. In this work, we present a fully automatic framework to estimate left ventricular myocardial deformation from t-MRI. This strategy performs automatic myocardial segmentation based on B-spline explicit active surfaces, which are initialized using an annular model. A non-rigid image-registration technique is then used to assess myocardial deformation. Three experiments were set up to validate the proposed framework using a clinical database of 75 patients. First, automatic segmentation accuracy was evaluated by comparing against manual delineations at one specific cardiac phase. The proposed solution showed an average perpendicular distance error of 2.35  ±  1.21 mm and 2.27  ±  1.02 mm for the endo- and epicardium, respectively. Second, starting from either manual or automatic segmentation, myocardial tracking was performed and the resulting strain curves were compared. It is shown that the automatic segmentation adds negligible differences during the strain-estimation stage, corroborating its accuracy. Finally, segmental strain was compared with scar tissue extent determined by delay-enhanced MRI. The results proved that both strain components were able to distinguish between normal and infarct regions. Overall, the proposed framework was shown to be accurate, robust, and attractive for clinical practice, as it overcomes several limitations of a manual analysis.

Homeostatic Left Heart integration and disintegration links atrio-ventricular covariation's dyshomeostasis in Hypertrophic Cardiomyopathy

Left ventricle and left atrium are and have been practically always analyzed separately in common clinically and non-clinically oriented cardiovascular investigations. Both classic and speckle tracking echocardiographic data contributed to the knowledge about deformational impairments occurring in systo-diastolic differences. Recently new trajectory based approaches allowed a greater awareness about the entire left ventricle or left atrium revolution and on their deficiencies that take place in presence of hypertrophic cardiomyopathy. However, surprisingly, the concomitant function of the two left heart chambers has not been analyzed for their geometrical/mechanical relationship. For the first time we study here, by acquiring left ventricle and left atrial geometries on the same heartbeat, the trajectory attributes of the entire left heart treated as a whole shape and the shape covariation of its two subunits. We contrasted healthy subjects with patients affected by hypertrophic cardiomyopathy. We found impaired left heart trajectory mainly in terms of orientation and size. More importantly, we found profound differences in the direction of morphological covariation of left ventricle and left atrium. These findings open to new perspectives in pathophysiological evaluation of different diseases by allowing the appreciation of concomitant functioning of both left heart whole geometry and of its two chambers.

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

BACKGROUND

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.

METHODS

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.

RESULTS

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.

CONCLUSIONS

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.