Statistical Analysis of a Three-dimensional Axial Strain and Axial-shear Strain Elastography Algorithm

Pathological phenomena often change the mechanical properties of the tissue. Therefore, estimation of tissue mechanical properties can be of clinical importance. Ultrasound elastography is a well-established strain estimation technique. Until recently, mainly 1D elastography algorithms have been developed. A few 2D algorithms have also been developed in the past. Both of these two types of technique ignore the tissue motion in the elevational direction, which could be a significant source of decorrelation in the RF data. In this thesis, a 3D elastography algorithm that estimates all the three components of tissue displacement is implemented and tested statistically. In this research, displacement fields of mechanical models are simulated. RF signals are then generated based on these displacement fields and used as the input of elastography algorithms. To evaluate the image quality of elastograms, absolute error, SNRe, CNRe and CNRasse are computed. The SNRe, CNRe and CNRasse values are investigated not only under different strain conditions, but also in different frame locations, which forms 3D strain filters. A statistical comparison between image qualities of the 3D technique and 2D technique is also provided. The results of this study show that the 3D elastography algorithm outperforms the 2D elastography algorithm in terms of image quality and robustness, especially under high strain conditions. This is because that the 3D algorithm estimates the elevational displacement, while the 2D technique only estimates the axial and lateral deformation. Since the elevational displacement could be an important source for the decorrelation in the RF data, the 3D technique is more effective and robust compared with the 2D technique

3D estimation of soft biological tissue deformation from radio-frequency ultrasound volume acquisitions

articl

Novel 3D Ultrasound Elastography Techniques for In Vivo Breast Tumor Imaging and Nonlinear Characterization

Breast cancer comprises about 29% of all types of cancer in women worldwide. This type of cancer caused what is equivalent to 14% of all female deaths due to cancer. Nowadays, tissue biopsy is routinely performed, although about 80% of the performed biopsies yield a benign result. Biopsy is considered the most costly part of breast cancer examination and invasive in nature. To reduce unnecessary biopsy procedures and achieve early diagnosis, ultrasound elastography was proposed.;In this research, tissue displacement fields were estimated using ultrasound waves, and used to infer the elastic properties of tissues. Ultrasound radiofrequency data acquired at consecutive increments of tissue compression were used to compute local tissue strains using a cross correlation method. In vitro and in vivo experiments were conducted on different tissue types to demonstrate the ability to construct 2D and 3D elastography that helps distinguish stiff from soft tissues. Based on the constructed strain volumes, a novel nonlinear classification method for human breast tumors is introduced. Multi-compression elastography imaging is elucidated in this study to differentiate malignant from benign tumors, based on their nonlinear mechanical behavior under compression. A pilot study on ten patients was performed in vivo, and classification results were compared with biopsy diagnosis - the gold standard. Various nonlinear parameters based on different models, were evaluated and compared with two commonly used parameters; relative stiffness and relative tumor size. Moreover, different types of strain components were constructed in 3D for strain imaging, including normal axial, first principal, maximum shear and Von Mises strains. Interactive segmentation algorithms were also evaluated and applied on the constructed volumes, to delineate the stiff tissue by showing its isolated 3D shape.;Elastography 3D imaging results were in good agreement with the biopsy outcomes, where the new classification method showed a degree of discrepancy between benign and malignant tumors better than the commonly used parameters. The results show that the nonlinear parameters were found to be statistically significant with p-value \u3c0.05. Moreover, one parameter; power-law exponent, was highly statistically significant having p-value \u3c 0.001. Additionally, volumetric strain images reconstructed using the maximum shear strains provided an enhanced tumor\u27s boundary from the surrounding soft tissues. This edge enhancement improved the overall segmentation performance, and diminished the boundary leakage effect. 3D segmentation provided an additional reliable means to determine the tumor\u27s size by estimating its volume.;In summary, the proposed elastographic techniques can help predetermine the tumor\u27s type, shape and size that are considered key features helping the physician to decide the sort and extent of the treatment. The methods can also be extended to diagnose other types of tumors, such as prostate and cervical tumors. This research is aimed toward the development of a novel \u27virtual biopsy\u27 method that may reduce the number of unnecessary painful biopsies, and diminish the increasingly risk of cancer

Tomographic measurement of all orthogonal components of three-dimensional displacement fields within scattering materials using wavelength scanning interferometry

Experimental mechanics is currently contemplating tremendous opportunities of further advancements thanks to a combination of powerful computational techniques and also fullfield non-contact methods to measure displacement and strain fields in a wide variety of materials. Identification techniques, aimed to evaluate material mechanical properties given known loads and measured displacement or strain fields, are bound to benefit from increased data availability (both in density and dimensionality) and efficient inversion methods such as finite element updating (FEU) and the virtual fields method (VFM). They work at their best when provided with dense and multicomponent experimental displacement (or strain) data, i.e. when all orthogonal components of displacements (or all components of the strain tensor) are known at points closely spaced within the volume of the material under study. Although a very challenging requirement, an increasing number of techniques are emerging to provide such data. In this Thesis, a novel wavelength scanning interferometry (WSI) system that provides three dimensional (3-D) displacement fields inside the volume of semi-transparent scattering materials is proposed. Sequences of two-dimensional interferograms are recorded whilst tuning the frequency of a laser at a constant rate. A new approach based on frequency multiplexing is used to encode the interference signal corresponding to multiple illumination directions at different spectral bands. Different optical paths along each illumination direction ensure that the signals corresponding to each sensitivity vector do not overlap in the frequency domain. All the information required to reconstruct the location and the 3-D displacement vector of scattering points within the material is thus recorded simultaneously in a single wavelength scan. By comparing phase data volumes obtained for two successive scans, all orthogonal components of the three dimensional displacement field introduced between scans (e.g. by means of loading or moving the sample under study) are readily obtained with high displacement sensitivity. The fundamental principle that describes the technique is presented in detail, including the correspondence between interference signal frequency and its associated depth within the sample, depth range, depth resolution, transverse resolution and displacement sensitivity. Data processing of the interference signal includes Fourier transformation, noise reduction, re-registration of data volumes, measurement of the illumination and sensitivity vectors from experimental data using a datum surface, phase difference evaluation, 3-D phase unwrapping and 3-D displacement field evaluation. Experiments consisting of controlled rigid body rotations and translations of a phantom were performed to validate the results. Both in-plane and the out-of-plane displacement components were measured for each voxel in the resulting data volume, showing an excellent agreement with the expected 3-D displacement

Development of novel ultrasound techniques for imaging and elastography. From simulation to real-time implementation

Ultrasound techniques offer many advantages, in terms of ease of realization and patients’ safety. The availability of suitable hardware and software tools is condicio sine qua non for new methods testing. This PhD project addresses medical ultrasound signal processing and seeks to achieve two scientific goals: the first is to contribute to the development of an ultrasound research platform, while the second is introducing and validating, through this platform, non-standard methods. During the thesis, the capabilities of the system were improved by creating advanced software tools, such as acoustic field simulators, and by developing echo-signals elaboration programs. In particular, a novel technique for quasi-static elastography was developed, in-vitro tested and implemented in real-time

Estimation du mouvement de la paroi carotidienne en imagerie ultrasonore par une approche de marquage ultrasonore

This work focuses on the processing of biomedical images. The aim of our study is to estimate the mechanical properties of the carotid artery in vivo using ultrasound imaging, in order to detect cardiovascular diseases at an early stage. Over the last decade, researchers have shown interest in studying artery wall motion, especially the motion of the carotid intima-media complex in order to demonstrate its significance as a marker of Atherosclerosis. However, despite recent progress, motion estimation of the carotid wall is still difficult, particularly in the longitudinal direction (direction parallel to the probe). The development of an innovative method for studying the movement of the carotid artery wall is the main motivation of this thesis. The three main contributions proposed in this work are i) the development, the validation, and the clinical evaluation of a novel method for 2D motion estimation of the carotid wall, ii) the development, the simulation and the experimental validation of the 3D extension of the estimation method proposed, and iii) the experimental evaluation of the 2D proposed method in ultra-fast imaging, for the estimation of the local pulse wave velocity. We propose a motion estimation method combining tagging of the ultrasound images, and a motion estimator based on the phase of the ultrasound images. The ultrasonic tagging is produced by means of transverse oscillations. We present two different approaches to introduce these transverses oscillations, a classic approach using a specific apodization function and a new approach based on filtering. The proposed motion estimator uses the 2D analytical phase of RF images using the Hahn approach. This thesis work shows that, compared with conventional methods, the proposed approach provides more accurate motion estimation in the longitudinal direction, and more generally in directions perpendicular to the beam axis. Also, the experimental evaluation of our method on ultra-fast images sequences from carotid phantom was used to validate our method regarding the estimation of the pulse wave velocity, the Young’s modulus of the vessels wall, and the propagation of a longitudinal movement.Ce travail de thèse est axé sur le domaine du traitement d’images biomédicales. L’objectif de notre étude est l’estimation des paramètres traduisant les propriétés mécaniques de l’artère carotide in vivo en imagerie échographique, dans une optique de détection précoce des pathologies cardiovasculaires. L’étude des comportements dynamiques de l’artère pour le dépistage précoce de l’athérosclérose constitue à ce jour une piste privilégiée. Cependant, malgré les avancées récentes, l’estimation du mouvement de la paroi carotidienne reste toujours difficile, notamment dans la direction longitudinale (direction parallèle au vaisseau). L’élaboration d’une méthode innovante permettant d’étudier le mouvement de la paroi carotidienne constitue la principale motivation de ce travail de thèse. Les trois contributions principales proposées dans ce travail sont i) le développement, la validation, et l’évaluation clinique d’une méthode originale d’estimation de mouvement 2D adaptée au mouvement de la paroi carotidienne, ii) la validation en simulation, et expérimentale de l’extension à la 3D de la méthode d’estimation proposée, et iii) l’évaluation expérimentale de la méthode proposée, en imagerie ultrasonore ultra-rapide, dans le cadre de l’estimation locale de la vitesse de l’onde de pouls. Nous proposons une méthode d’estimation de mouvement combinant un marquage ultrasonore dans la direction latérale, et un estimateur de mouvement basé sur la phase des images ultrasonores. Le marquage ultrasonore est réalisé par l’intermédiaire d’oscillations transverses. Nous proposons deux approches différentes pour introduire ces oscillations transverses, une approche classique utilisant une fonction de pondération spécifique, et une approche originale par filtrage permettant de contrôler de manière optimale leurs formations. L’estimateur de mouvement proposé utilise les phases analytiques des images radiofréquences, extraites par l’approche de Hahn. Ce travail de thèse montre que la méthode proposée permet une estimation de mouvement plus précise dans la direction longitudinale, et plus généralement dans les directions perpendiculaires au faisceau ultrasonore, que celle obtenue avec d’autres méthodes plus traditionnelles. De plus, l’évaluation expérimentale de la méthode sur des séquences d’images ultrasonores ultra-rapides issues de fantômes de carotide, a permis l’estimation locale de la vitesse de propagation de l’onde de pouls, la mise en évidence de la propagation d’un mouvement longitudinal et enfin l’estimation du module de Young des vaisseaux