In vivo validation of a novel method for regional myocardial wall motion analysis based on echocardiographic tissue tracking

Optimization-Based Speckle Tracking Algorithm for Left Ventricle Strain Estimation: A Feasibility Study

Speckle tracking echocardiography (STE) is a widespread method for calculating myocardial strains and estimating left ventricle function. Since echocardiographic clips are corrupted by speckle decorrelation noise, resulting in irregular, nonphysiological tissue displacement fields, smoothing is performed on the displacement data, affecting the strain results. Thus, strain results may depend on the specific implementations of 2-D STE, as well as other systems' characteristics of the various vendors. A novel algorithm (called K-SAD) is introduced, which integrates the physiological constraint of smoothness of the displacement field into an optimization process. Simulated B-mode clips, modeling healthy and abnormal cases, were processed by K-SAD. Peak global and subendocardial longitudinal strains, as well as regional strains, were calculated. In addition, 410 healthy subjects were also processed. The results of K-SAD are compared with those of one of the leading commercial product. K-SAD provides global mid-wall strain values, as well as subendocardial and regional strain values, all in good agreement with the ground-truth-simulated phantom data. K-SAD peak global longitudinal systolic strain values for 410 healthy subjects are quite similar for the different regions: - 17.02 ± 4.02%, - 19.00 ± 3.45%, and - 19.72 ± 5.06% at the basal, mid, and apical regions, respectively. Improved performance under noisy conditions was demonstrated by comparing a subgroup of 40 subjects with the best image quality with the remaining 370 cohort: K-SAD provides statistically similar global and regional results for the two cohorts. Our study indicates that the sensitivity of strain values to speckle noise, caused by the post block-matching weighted smoothing, can be significantly reduced and accuracy enhanced by employing an integrated one-stage, physiologically constrained optimization process.

Is strain by Speckle Tracking Echocardiography dependent on user controlled spatial and temporal smoothing? An experimental porcine study

Background

Speckle Tracking Echocardiography (STE) strain analysis relies on both spatial and temporal smoothing. The user is often allowed to adjust these smoothing parameters during analysis. This experimental study investigates how different degrees of user controllable spatial and temporal smoothing affect global and regional STE strain values in recordings obtained from normal and ischemic myocardium.

Methods

In seven anesthetized pigs, left ventricular short- and long-axis B-mode cineloops were recorded before and after left anterior descending coronary artery occlusion. Peak- and postsystolic global STE strain in the radial, circumferential and longitudinal direction as well as corresponding regional strain in the anterior and posterior walls were measured. During post-processing, strain values were obtained with three different degrees of both spatial and temporal smoothing (minimum, factory default and maximum), resulting in nine different combinations.

Results

All parameters for global and regional longitudinal strain were unaffected by adjustments of spatial and temporal smoothing in both normal and ischemic myocardium. Radial and circumferential strain depended on smoothing to a variable extent, radial strain being most affected. However, in both directions the different combinations of smoothing did only result in relatively small changes in the strain values. Overall, the maximal strain difference was found in normal myocardium for peak systolic radial strain of the posterior wall where strain was 22.0 ± 2.2% with minimal spatial and maximal temporal smoothing and 30.9 ± 2.6% with maximal spatial and minimal temporal smoothing (P < 0.05).

Conclusions

Longitudinal strain was unaffected by different degrees of user controlled smoothing. Radial and circumferential strain depended on the degree of smoothing. However, in most cases these changes were small and would not lead to altered conclusions in a clinical setting. Furthermore, smoothing did not affect strain variance. For all strain parameters, variance remained within the corresponding interobserver variance.

Two-dimensional strain imaging of controlled rabbit hearts

Ultrasound strain imaging using 2-D speckle tracking has been proposed to quantitatively assess changes in myocardial contractility caused by ischemia. Its performance must be demonstrated in a controlled model system as a step toward routine clinical application. In this study, a well-controlled 2-D cardiac elasticity imaging technique was developed using two coplanar and orthogonal linear probes simultaneously imaging an isolated retroperfused rabbit heart. Acute ischemia was generated by left anterior descending (LAD) artery ligation. An excitation-contraction decoupler, 2,3-butanedione monoxime, was applied at a 4-mM concentration to reversibly reduce myocardial contractility. Results using a single probe demonstrate that directional changes in the in-plane principal deformation axes can help locate the bulging area as a result of LAD ligation, which matched well with corresponding Evans Blue staining, and strains or strain magnitude, based on principal stretches, can characterize heart muscle contractility. These two findings using asymmetric displacement accuracy (i.e., normal single-probe measurements with good axial but poor lateral estimates) were further validated using symmetric displacement accuracy (i.e., dual-probe measurements using only accurate axial tracking estimates from each). However, the accuracy of 2-D cardiac strain imaging using a single probe depends on the probe's orientation because of the large variance in lateral displacement estimates.

Phase rotation methods in filtering correlation coefficients for ultrasound speckle tracking

In speckle-tracking-based myocardial strain imaging, large interframe/volume peak-systolic strains cause peak hopping artifacts separating the highest correlation coefficient peak from the true peak. A correlation coefficient filter was previously designed to minimize peak hopping artifacts. For large strains, however, the correlation coefficient filter must follow the strain distribution to remove peak hopping effectively. This processing usually means interpolation and high computational load. To reduce the computational burden, a narrow band approximation using phase rotation is developed in this paper to facilitate correlation coefficient filtering. Correlation coefficients are first phase rotated to increase coherence, then filtered. Rotated phase angles are determined by the local strain and spatial position. This form of correlation coefficient filtering enhances true correlation coefficient peaks in large strain applications if decorrelation due to deformation does not completely destroy the coherence among neighboring correlation coefficients. The assumed strain used in the filter can also deviate from the true strain and still be effective. Further improvement in displacement estimation can be expected by combining correlation coefficient filtering with a new Viterbi-based displacement estimator.

Reduced peak-hopping artifacts in ultrasonic strain estimation using the Viterbi algorithm

Internal strain resulting from tissue deformation can be estimated by correlation processing of speckle patterns within complex (i.e., radio frequency) ultrasound images acquired during deformation. At large deformations, the magnitude of the correlation coefficient peak can be significantly lower than unity, so that random speckle correlations will exceed the true peak. This effect is called "peak hopping" and causes significant errors in displacement and deformation estimates. Here we investigate the Viterbi algorithm, a dynamic programming procedure, to overcome peak-hopping artifacts by finding the most likely sequence of hidden states in a sequence of observed events. It is well suited to motion estimation in elasticity-imaging studies because adjacent tissue elements remain adjacent following deformation. Particularly, tissue elements along an ultrasonic beam in one image lie along a 3-D continuous curve in the next image instant. The observed event in this case is the correlation coefficient of a pixel at a certain displacement. Radio-frequency data were generated before and after deformation with an average strain of 6%. Simulations were performed for a homogenous medium and for a medium with a stiffer inclusion. Results show that Viterbi processing of speckle-tracking outputs can significantly reduce peak-hopping artifacts.