Segmentation of 3D radio frequency echocardiography using a spatio-temporal predictor

Three-Dimensional Model-Based Segmentation in Echocardiography Using High Temporal Tissue and Blood Flow Information

Accurate 3-D surface segmentation is a challenging task in echocardiography because of the relatively low image quality. We introduce a new method for 3-D segmentation of the endocardium involving temporal decorrelation of echo signals originating from tissue and blood using radiofrequency (RF) signals acquired in 3-D Doppler mode. Temporal features were extracted in 3-D Doppler mode, where a sequence of RF lines is recorded for each image line. Each set of RF lines is highly correlated because of the high pulse repetition frequency. However, for high blood flow, the RF signals will decorrelate over time in contrast to the endocardium, which will remain relatively highly correlated over time. These decorrelation features permit differentiation between myocardial tissue and blood flow. We describe an implementation of a 3-D segmentation model in which temporal information is used as external constraint. The model was validated in a phantom and in vivo in healthy volunteers (n = 5). The phantom study revealed that the model successfully segmented the artificial blood lumen even for low flow velocity and illustrated the sensitivity of the segmentations to flow rate. In healthy volunteers, high Dice similarity indices indicate that 3-D segmentation of the endocardial border in vivo is feasible.

A new inverse method for estimation of in vivo mechanical properties of the aortic wall

The aortic wall is always loaded in vivo, which makes it challenging to estimate the material parameters of its nonlinear, anisotropic constitutive equation from in vivo image data. Previous approaches largely relied on either computationally expensive finite element models or simplifications of the geometry or material models. In this study, we investigated a new inverse method based on aortic wall stress computation. This approach consists of the following two steps: (1) computing an "almost true" stress field from the in vivo geometries and loading conditions, (2) building an objective function based on the "almost true" stress fields, constitutive equations and deformation relations, and estimating the material parameters by minimizing the objective function. The method was validated through numerical experiments by using the in vivo data from four ascending aortic aneurysm (AsAA) patients. The results demonstrated that the method is computationally efficient. This novel approach may facilitate the personalized biomechanical analysis of aortic tissues in clinical applications, such as in the rupture risk analysis of ascending aortic aneurysms.

Simulating cardiac ultrasound image based on MR diffusion tensor imaging

PURPOSE

Cardiac ultrasound simulation can have important applications in the design of ultrasound systems, understanding the interaction effect between ultrasound and tissue and setting the ground truth for validating quantification methods. Current ultrasound simulation methods fail to simulate the myocardial intensity anisotropies. New simulation methods are needed in order to simulate realistic ultrasound images of the heart.

METHODS

The proposed cardiac ultrasound image simulation method is based on diffusion tensor imaging (DTI) data of the heart. The method utilizes both the cardiac geometry and the fiber orientation information to simulate the anisotropic intensities in B-mode ultrasound images. Before the simulation procedure, the geometry and fiber orientations of the heart are obtained from high-resolution structural MRI and DTI data, respectively. The simulation includes two important steps. First, the backscatter coefficients of the point scatterers inside the myocardium are processed according to the fiber orientations using an anisotropic model. Second, the cardiac ultrasound images are simulated with anisotropic myocardial intensities. The proposed method was also compared with two other nonanisotropic intensity methods using 50 B-mode ultrasound image volumes of five different rat hearts. The simulated images were also compared with the ultrasound images of a diseased rat heart in vivo. A new segmental evaluation method is proposed to validate the simulation results. The average relative errors (AREs) of five parameters, i.e., mean intensity, Rayleigh distribution parameter σ, and first, second, and third quartiles, were utilized as the evaluation metrics. The simulated images were quantitatively compared with real ultrasound images in both ex vivo and in vivo experiments.

RESULTS

The proposed ultrasound image simulation method can realistically simulate cardiac ultrasound images of the heart using high-resolution MR-DTI data. The AREs of their proposed method are 19% for the mean intensity, 17.7% for the scale parameter of Rayleigh distribution, 36.8% for the first quartile of the image intensities, 25.2% for the second quartile, and 19.9% for the third quartile. In contrast, the errors of the other two methods are generally five times more than those of their proposed method.

CONCLUSIONS

The proposed simulation method uses MR-DTI data and realistically generates cardiac ultrasound images with anisotropic intensities inside the myocardium. The ultrasound simulation method could provide a tool for many potential research and clinical applications in cardiac ultrasound imaging.

DTI template-based estimation of cardiac fiber orientations from 3D ultrasound

PURPOSE

Cardiac muscle fibers directly affect the mechanical, physiological, and pathological properties of the heart. Patient-specific quantification of cardiac fiber orientations is an important but difficult problem in cardiac imaging research. In this study, the authors proposed a cardiac fiber orientation estimation method based on three-dimensional (3D) ultrasound images and a cardiac fiber template that was obtained from magnetic resonance diffusion tensor imaging (DTI).

METHODS

A DTI template-based framework was developed to estimate cardiac fiber orientations from 3D ultrasound images using an animal model. It estimated the cardiac fiber orientations of the target heart by deforming the fiber orientations of the template heart, based on the deformation field of the registration between the ultrasound geometry of the target heart and the MRI geometry of the template heart. In the experiments, the animal hearts were imaged by high-frequency ultrasound, T1-weighted MRI, and high-resolution DTI.

RESULTS

The proposed method was evaluated by four different parameters: Dice similarity coefficient (DSC), target errors, acute angle error (AAE), and inclination angle error (IAE). Its ability of estimating cardiac fiber orientations was first validated by a public database. Then, the performance of the proposed method on 3D ultrasound data was evaluated by an acquired database. Their average values were 95.4% ± 2.0% for the DSC of geometric registrations, 21.0° ± 0.76° for AAE, and 19.4° ± 1.2° for IAE of fiber orientation estimations. Furthermore, the feasibility of this framework was also performed on 3D ultrasound images of a beating heart.

CONCLUSIONS

The proposed framework demonstrated the feasibility of using 3D ultrasound imaging to estimate cardiac fiber orientation of in vivo beating hearts and its further improvements could contribute to understanding the dynamic mechanism of the beating heart and has the potential to help diagnosis and therapy of heart disease.

Radial basis functions for combining shape and speckle tracking in 4D echocardiography

Quantitative analysis of left ventricular deformation can provide valuable information about the extent of disease as well as the efficacy of treatment. In this work, we develop an adaptive multi-level compactly supported radial basis approach for deformation analysis in 3D+time echocardiography. Our method combines displacement information from shape tracking of myocardial boundaries (derived from B-mode data) with mid-wall displacements from radio-frequency-based ultrasound speckle tracking. We evaluate our methods on open-chest canines (N=8) and show that our combined approach is better correlated to magnetic resonance tagging-derived strains than either individual method. We also are able to identify regions of myocardial infarction (confirmed by postmortem analysis) using radial strain values obtained with our approach.

Automatic segmentation of right ventricular ultrasound images using sparse matrix transform and a level set

An automatic segmentation framework is proposed to segment the right ventricle (RV) in echocardiographic images. The method can automatically segment both epicardial and endocardial boundaries from a continuous echocardiography series by combining sparse matrix transform, a training model, and a localized region-based level set. First, the sparse matrix transform extracts main motion regions of the myocardium as eigen-images by analyzing the statistical information of the images. Second, an RV training model is registered to the eigen-images in order to locate the position of the RV. Third, the training model is adjusted and then serves as an optimized initialization for the segmentation of each image. Finally, based on the initializations, a localized, region-based level set algorithm is applied to segment both epicardial and endocardial boundaries in each echocardiograph. Three evaluation methods were used to validate the performance of the segmentation framework. The Dice coefficient measures the overall agreement between the manual and automatic segmentation. The absolute distance and the Hausdorff distance between the boundaries from manual and automatic segmentation were used to measure the accuracy of the segmentation. Ultrasound images of human subjects were used for validation. For the epicardial and endocardial boundaries, the Dice coefficients were 90.8 ± 1.7% and 87.3 ± 1.9%, the absolute distances were 2.0 ± 0.42 mm and 1.79 ± 0.45 mm, and the Hausdorff distances were 6.86 ± 1.71 mm and 7.02 ± 1.17 mm, respectively. The automatic segmentation method based on a sparse matrix transform and level set can provide a useful tool for quantitative cardiac imaging.

Automatic Segmentation of Right Ventricle on Ultrasound Images Using Sparse Matrix Transform and Level Set

An automatic framework is proposed to segment right ventricle on ultrasound images. This method can automatically segment both epicardial and endocardial boundaries from a continuous echocardiography series by combining sparse matrix transform (SMT), a training model, and a localized region based level set. First, the sparse matrix transform extracts main motion regions of myocardium as eigenimages by analyzing statistical information of these images. Second, a training model of right ventricle is registered to the extracted eigenimages in order to automatically detect the main location of the right ventricle and the corresponding transform relationship between the training model and the SMT-extracted results in the series. Third, the training model is then adjusted as an adapted initialization for the segmentation of each image in the series. Finally, based on the adapted initializations, a localized region based level set algorithm is applied to segment both epicardial and endocardial boundaries of the right ventricle from the whole series. Experimental results from real subject data validated the performance of the proposed framework in segmenting right ventricle from echocardiography. The mean Dice scores for both epicardial and endocardial boundaries are 89.1%±2.3% and 83.6±7.3%, respectively. The automatic segmentation method based on sparse matrix transform and level set can provide a useful tool for quantitative cardiac imaging.