Reconstructions of shear modulus, Poisson's ratio, and density using approximate mean normal stress lambda epsilon alpha alpha as unknown

As a differential diagnosis technique for living soft tissues, we are developing ultrasonic-strain-measurement-based shear modulus reconstruction methods. Previously, we reported three-dimensional (3-D) and 2-D reconstruction methods utilizing a typical Poisson's ratio very close to 0.5 (nearly-incompressible). However, because a decrease in the accuracy of the reconstructed value was confirmed to be due to the difference between the original value and the set value, we proposed 3-D and 2-D methods of reconstructing Poisson's ratio as well. Furthermore, we proposed methods of reconstructing density and dealing with dynamic deformation. However, due to tissue incompressibility, the reconstructions of shear modulus, Poisson's ratio, and density became unstable. In this report, to obtain stable, unique reconstructions, we describe a new reconstruction method using mean normal stress approximated by the product of one of Lame's constants X and volume strain epsilon alpha alpha as an unknown. Regularization is simultaneously applied to the respective distributions to decrease the instability of the reconstructions due to measurement errors of the deformation. This method also enables stable, unique reconstructions of shear modulus and density under the condition that the mean normal stress remains unknown. We also verify the effectiveness of this method through 3-D simulations, while showing erroneous artifacts occurring when 2-D and 1-D reconstructions are performed.     

Ultrafast compound imaging for 2-D motion vector estimation: application to transient elastography

This paper describes a new technique for two-dimensional (2-D) imaging of the motion vector at a very high frame rate with ultrasound. Its potential is experimentally demonstrated for transient elastography. But, beyond this application, it also could be promising for color flow and reflectivity imaging. To date, only axial displacements induced in human tissues by low-frequency vibrators were measured during transient elastography. The proposed technique allows us to follow both axial and lateral displacements during the shear wave propagation and thus should improve Young's modulus image reconstruction. The process is a combination of several ideas well-known in ultrasonic imaging: ultra-fast imaging, multisynthetic aperture beamforming, 1-D speckle tracking, and compound imaging. Classical beamforming in the transmit mode is replaced here by a single plane wave insonification increasing the frame rate by at least a factor of 128. The beamforming is achieved only in the receive mode on two independent subapertures. Comparison of successive frames by a classical 1-D speckle tracking algorithm allows estimation of displacements along two different directions linked to the subapertures beams. The variance of the estimates is finally improved by tilting the emitting plane wave at each insonification, thus allowing reception of successive decorrelated speckle patterns.     

Assessment of elastic parameters of human skin using dynamic elastography

Sonoelastography and transient elastography are two ultrasound-based techniques that facilitate noninvasive characterization of the viscoelastic properties of soft tissues by investigating their response to shear mechanical excitation. Young's modulus is the principle assessment parameter. Because it defines local tissue stiffness, it is of major interest for the medical imaging and cosmetic industries as it could replace subjective palpation by yielding local, quantitative information. In this paper, we describe a new high-resolution device capable of measuring local Young's modulus in very thin layers (1-5 mm) and devoted to the in vivo evaluation of the elastic properties of human skin. It uses an ultrasonic probe (50 MHz) for tracking the displacements induced by a 300 Hz shear wave generated by a ring surrounding the transducer. The displacements are measured using a conventional cross-correlation technique between successive ultrasonic back-scattered echoes. First, this noninvasive technique has been experimentally proven to be accurate for investigating elasticity in different skin-mimicking phantoms. Second, data were acquired in vivo on human forearms. As expected, Young's modulus was found to be higher in the dermis than in the hypodermis and other soft tissues.     

Shear modulus imaging with 2-D transient elastography

In previous work, we have shown that time-resolved 2-D transient elastography is a promising technique for characterizing the elasticity of soft tissues. It involves the measurement of the displacements induced by the propagation of low frequency (LF) pulsed shear waves in biological tissues. In this paper, we present a novel apparatus that contains a LF vibrating device surrounding a linear array of 128 ultrasonic transducers that performs ultrafast ultrasonic imaging (up to 10,000 frames/s) and that is able to follow in real time the propagation of a LF shear wave in the human body. The vibrating device is made of two rods, fixed to electromagnetic vibrators, that produce in the ultrasonic image area a large amplitude shear wave. The geometry has been chosen both to enhance the sensitivity and to create a quasi linear shear wave front in the imaging plane. An inversion algorithm is used to recover the shear modulus map from the spatio-temporal data, and the first experimental results obtained from tissue-equivalent materials are presented.     

Regularization of tissue shear modulus reconstruction using strain variance

An effective setting method (that is, a method using the variances of strain tensor component measurements) is described for the properly spatially varied regularization parameters for our shear modulus reconstruction. At each position, the respective strain variances can be experimentally evaluated using plural field measurements or single field measurement, for example, when using all crosscorrelation- based methods, by using the Ziv-Zakai Lower Bound (ZZLB). The demonstrated regularization by the single field measurement using the cross-spectrum phase gradient method (MCSPGM) in experiments confirms that the use of the axial strain variance estimated by the echo signal-to-noise ratio and correlations (the combined SNRc) effectively stabilizes the 1-D reconstruction on an agar phantom and a human in vivo liver carcinoma during interstitial microwave thermal treatment. The regularization yields a spatially uniform stability in reconstruction.     

Monitoring of thermal therapy based on shear modulus changes: I. shear wave thermometry

The clinical applicability of high-intensity focused ultrasound (HIFU) for noninvasive therapy is today hampered by the lack of robust and real-time monitoring of tissue damage during treatment. The goal of this study is to show that the estimation of local tissue elasticity from shear wave imaging (SWI) can lead to the 2-D mapping of temperature changes during HIFU treatments. This new concept of shear wave thermometry is experimentally implemented here using conventional ultrasonic imaging probes. HIFU treatment and monitoring were, respectively, performed using a confocal setup consisting of a 2.5-MHz single-element transducer focused at 30 mm on ex vivo samples and an 8-MHz ultrasound diagnostic probe. Thermocouple measurements and ultrasound-based thermometry were used as a gold standard technique and were combined with SWI on the same device. The SWI sequences consisted of 2 successive shear waves induced at different lateral positions. Each wave was created using 100-μs pushing beams at 3 depths. The shear wave propagation was acquired at 17,000 frames/s, from which the elasticity map was recovered. HIFU sonications were interleaved with fast imaging acquisitions, allowing a duty cycle of more than 90%. Elasticity and temperature mapping was achieved every 3 s, leading to realtime monitoring of the treatment. Tissue stiffness was found to decrease in the focal zone for temperatures up to 43°C. Ultrasound-based temperature estimation was highly correlated to stiffness variation maps (r2 = 0.91 to 0.97). A reversible calibration phase of the changes of elasticity with temperature can be made locally using sighting shots. This calibration process allows for the derivation of temperature maps from shear wave imaging. Compared with conventional ultrasound-based approaches, shear wave thermometry is found to be much more robust to motion artifacts.     

3-D FDTD simulation of shear waves for evaluation of complex modulus imaging

The Navier equation describing shear wave propagation in 3-D viscoelastic media is solved numerically with a finite differences time domain (FDTD) method. Solutions are formed in terms of transverse scatterer velocity waves and then verified via comparison to measured wave fields in heterogeneous hydrogel phantoms. The numerical algorithm is used as a tool to study the effects on complex shear modulus estimation from wave propagation in heterogeneous viscoelastic media. We used an algebraic Helmholtz inversion (AHI) technique to solve for the complex shear modulus from simulated and experimental velocity data acquired in 2-D and 3-D. Although 3-D velocity estimates are required in general, there are object geometries for which 2-D inversions provide accurate estimations of the material properties. Through simulations and experiments, we explored artifacts generated in elastic and dynamic-viscous shear modulus images related to the shear wavelength and average viscosity.     

On the effects of reflected waves in transient shear wave elastography

In recent years, novel quantitative techniques have been developed to provide noninvasive and quantitative stiffness images based on shear wave propagation. Using radiation force and ultrafast ultrasound imaging, the supersonic shear imaging technique allows one to remotely generate and follow a transient plane shear wave propagating in vivo in real time. The tissue shear modulus, i.e., its stiffness, can then be estimated from the shear wave local velocity. However, because the local shear wave velocity is estimated using a time-of- flight approach, reflected shear waves can cause artifacts in the estimated shear velocity because the incident and reflected waves propagate in opposite directions. Such effects have been reported in the literature as a potential drawback of elastography techniques based on shear wave speed, particularly in the case of high stiffness contrasts, such as in atherosclerotic plaque or stiff lesions. In this letter, we present our implementation of a simple directional filter, previously used for magnetic resonance elastography, which separates the forward- and backward-propagating waves to solve this problem. Such a directional filter could be applied to many elastography techniques based on the local estimation of shear wave speed propagation, such as acoustic radiation force imaging (ARFI), shearwave dispersion ultrasound vibrometry (SDUV), needle-based elastography, harmonic motion imaging, or crawling waves when the local propagation direction is known and high-resolution spatial and temporal data are acquired.     

Time-resolved diffuse optical tomographic imaging for the provision of both anatomical and functional information about biological tissue

We present in vivo images of near-infrared (NIR) diffuse optical tomography (DOT) of human lower legs and forearm to validate the dual functions of a time-resolved (TR) NIR DOT in clinical diagnosis, i.e., to provide anatomical and functional information simultaneously. The NIR DOT system is composed of time-correlated single-photon-counting channels, and the image reconstruction algorithm is based on the modified generalized pulsed spectral technique, which effectively incorporates the TR data with reasonable computation time. The reconstructed scattering images of both the lower legs and the forearm revealed their anatomies, in which the bones were clearly distinguished from the muscles. In the absorption images, some of the blood vessels were observable. In the functional imaging, a subject was requested to do handgripping exercise to stimulate physiological changes in the forearm tissue. The images of oxyhemoglobin, deoxyhemoglobin, and total hemoglobin concentration changes in the forearm were obtained from the differential images of the absorption at three wavelengths between the exercise and the rest states, which were reconstructed with a differential imaging scheme. These images showed increases in both blood volume and oxyhemoglobin concentration in the arteries and simultaneously showed hypoxia in the corresponding muscles. All the results have demonstrated the capability of TR NIR DOT by reconstruction of the absolute images of the scattering and the absorption with a high spatial resolution that finally provided both the anatomical and functional information inside bulky biological tissues.     

Fine elasticity imaging utilizing the iterative RF-echo phase matching method

To non-invasively quantify elasticity of soft tissue, we previously developed the iterative two-dimensional (2-D) rf-echo phase matching method for accurately measuring a 2-D displacement vector field generated in vivo in soft tissue during acquisition of two successive rf-echo data frames. We also developed a stable method for uniquely reconstructing a shear modulus distribution using strains derived from the measurement data. However, as in our measurement method a displacement is determined by using the phase characteristics of the finite local echo data as the index to iteratively search for the corresponding local data, change of the local phase characteristics due to tissue deformation deteriorates the accuracy of the determination. Thus, we improve the previous method such that, in principle, the displacement can be determined using an infinitesimal phase characteristics. That is, we incorporate an effective mechanism into the previous iterative phase matching scheme: the local size is made suitably smaller during the iterative phase matching. The demonstrated ability of measurement and reconstruction in simulation, and experiments on in vitro in pork rib and in vivo in breast tissue, shows this refinement allows not only better spatial resolution of the shear modulus image but also improved accuracy, and indicates that the improved method has a high potential to be applied for various soft tissues.     

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.     

Target detectability in acoustic elastography

The clinically relevant task of visually detecting low contrast targets in noisy strain images estimated from ultrasonic signals is studied. Detectability is measured quantitatively using contrast-to-noise ratio (CNR) analysis. Contrast in strain images is generated by a complex interaction among the soft tissue elasticity shear modulus distribution, target shape and location in the stress field, and external boundary conditions. Although a large strain variation is preferred for enhancing the contrast, this also increases the signal-dependent noise in strain estimates in a nonlinear fashion. Therefore, understanding the tradeoffs between contrast and noise is necessary for improving the diagnostic performance of strain imaging. In this paper, targets with slab, cylindrical, and spherical geometries are studied. Strains in the target and background and the precision of their estimates are described in terms of the corresponding shear modulus values for each geometry. These results are then incorporated into the CNR expression to investigate the changes in target detectability with the variation of shear modulus in the target and the ultrasonic signal parameters (echo signal-to-noise ratio and inverse fractional bandwidth) as well as the signal processing variables (time-bandwidth product and fractional window overlap).     

Electrical Capacitance Volume Tomography

A dynamic volume imaging based on the principle of electrical capacitance tomography (ECT), namely, electrical capacitance volume tomography (ECVT), has been developed in this study. The technique generates, from the measured capacitance, a whole volumetric image of the region enclosed by the geometrically three-dimensional capacitance sensor. This development enables a real-time, 3-D imaging of a moving object or a real-time volume imaging (4-D) to be realized. Moreover, it allows total interrogation of the whole volume within the domain (vessel or conduit) of an arbitrary shape or geometry. The development of the ECVT imaging technique primarily encloses the 3-D capacitance sensor design and the volume image reconstruction technique. The electrical field variation in three-dimensional space forms a basis for volume imaging through different shapes and configurations of ECT sensor electrodes. The image reconstruction scheme is established by implementing the neural-network multicriterion optimization image reconstruction (NN-MOIRT), developed earlier by the authors for the 2-D ECT. The image reconstruction technique is modified by introducing into the algorithm a 3-D sensitivity matrix to replace the 2-D sensitivity matrix in conventional 2-D ECT, and providing additional network constraints including 3-to-2-D image matching function. The additional constraints further enhance the accuracy of the image reconstruction algorithm. The technique has been successfully verified over actual objects in the experimental conditions     

采用稀疏测量数据的有限角度光声层析成像的研究进展

生物光声层析(Photoacoustic Tomography,PAT)成像可以反映生物组织的光吸收分布,定量测量组织的光吸收系数和散射系数,进而分析组织成分,为疾病的早期诊断和治疗提供可靠的依据。由于成像目标特殊的几何结构以及成像装置的机械结构、空间位置和成像时间等的限制,超声探测器只能在有限的角度范围内扫描,采集到稀疏的光声测量数据,导致重建图像中出现伪影和失真。针对有限角度扫描和稀疏测量数据问题,对目前主流的光声图像重建算法进行综述和分析。     

Two-wave propagation imaging to evaluate the structure of cancellous bone

The two-wave phenomenon reflects not only bone mass but also the complex bone structure of cancellous bone. We propose a new simple imaging technique based on the two-wave phenomenon for investigating the anisotropic structure of cancellous bone. A cylindrical specimen of cancellous bone was obtained from a bovine femur. The structure (alignment of trabeculae) of the specimen was obtained from 3-D X-ray micro computed tomography imaging. Using a conventional ultrasonic pulse technique, we rotated the receiver around the specimen to investigate the ultrasonic fields after propagation within the specimen. The ultrasonic propagation image clearly showed the effect of the bone structure. We found that the fast wave showed apparent refraction, whereas the slow wave did not. Fast-wave propagation imaging is thus a simple and convenient technique for easy interpretation of the anisotropic structure. This imaging technique has the potential to become a powerful tool to investigate the structure of trabeculae during in vivo measurements.     

基于超声相控阵的水下船体表面成像方法研究

针对当前水下船体表面成像设备较少、成本较高、成像误差较大的问题,提出一种基于超声相控阵技术的水下船体成像方法。以超声相控阵技术为基础,控制一维线型超声相控阵列换能器发射的聚焦声波对水下船体表面进行相控扫描,声波作用于船体表面后发生散射,回波被换能器接收后利用双曲线交汇法计算出各反射点的点云数据,最终利用点云数据重构出水下船体表面的三维模型。通过搭建水池船模成像实验对方法进行了验证,结果表明,重构图像最大误差小于5mm,具有较高的成像分辨力。     

Nonlinear elasticity imaging: theory and phantom study

In tissue the Young's modulus cannot be assumed constant over a wide deformation range. For example, direct mechanical measurements on human prostate show up to a threefold increase in Young's modulus over a 10% deformation. In conventional elasticity imaging, these effects produce strain-dependent elastic contrast. Ignoring these effects generally leads to suboptimal contrast (stiffer tissues at lower strain are contrasted against softer tissues at higher strain), but measuring the nonlinear behavior results in enhanced tissue differentiation. To demonstrate the methods extracting nonlinear elastic properties, both simulations and measurements were performed on an agar-gelatin phantom. Multiple frames of phase-sensitive ultrasound data are acquired as the phantom is deformed by 12%. All interframe displacement data are brought back to the geometry of the first frame to form a three-dimensional (3-D) data set (depth, lateral, and preload dimensions). Data are fit to a 3-D second order polynomial model for each pixel that adjusts for deformation irregularities. For the phantom geometry and elastic properties considered in this paper, reconstructed frame-to-frame strain images using this model result in improved contrast to noise ratios (CNR) at all preload levels, without any sacrifice in spatial resolution. From the same model, strain hardening at all preload levels can be extracted. This is an independent contrast mechanism. Its maximum CNR occurs at 5.13% preload, and it is a 54% improvement over the best case (preload 10.6%) CNR for frame-to-frame strain reconstruction. Actual phantom measurements confirm the essential features of the simulation. Results show that modeling of the nonlinear elastic behavior has the potential to both increase detectability in elasticity imaging and provide a new independent mechanism for tissue differentiation.     

Application of 1-D transient elastography for the shear modulus assessment of thin-layered soft tissue: comparison with supersonic shear imaging technique

Elasticity estimation of thin-layered soft tissues has gained increasing interest propelled by medical applications like skin, corneal, or arterial wall shear modulus assessment. In this work, the authors propose one-dimensional transient elastography (1DTE) for the shear modulus assessment of thin-layered soft tissue. Experiments on three phantoms with different elasticities and plate thicknesses were performed. First, using 1DTE, the shear wave speed dispersion curve inside the plate was obtained and validated with finite difference simulation. No dispersive effects were observed and the shear wave speed was directly retrieved from time-of-flight measurements. Second, the supersonic shear imaging (SSI) technique (considered to be a gold standard) was performed. For the SSI technique, the propagating wave inside the plate is guided as a Lamb wave. Experimental SSI dispersion curves were compared with finite difference simulation and fitted using a generalized Lamb model to retrieve the plate bulk shear wave speed. Although they are based on totally different mechanical sources and induce completely different diffraction patterns for the shear wave propagation, the 1DTE and SSI techniques resulted in similar shear wave speed estimations. The main advantage of the 1DTE technique is that bulk shear wave speed can be directly retrieved without requiring a dispersion model.     

Theoretical quality assessment of myocardial elastography with in vivo validation

Myocardial elastography (ME), a radio frequency (RF)-based speckle tracking technique with one-dimensional (1-D) cross correlation and novel recorrelation methods in a 2-D search was proposed to estimate and fully image 2-D transmural deformation field and to detect abnormal cardiac function. A theoretical framework was first developed in order to evaluate the performance of 2-D myocardial elastography based on a previously developed 3-D finite-element model of the canine left ventricle. A normal (control) and an ischemic (left-circumflex, LCx) model, which more completely represented myocardial deformation than a kinematic model, were considered. A 2-D convolutional image formation model was first used to generate RF signals for quality assessment of ME in the normal and ischemic cases. A 3-D image formation model was further developed to investigate the effect of the out-of-plane motion on the 2-D, in-plane motion estimation. Both orthogonal, in-plane displacement components (i.e., lateral and axial) between consecutive RF frames were iteratively estimated. All the estimated incremental 2-D displacements from end-diastole (ED) to end-systole (ES) were then accumulated to acquire the cumulative 2-D displacements, which were further used to calculate the cumulative 2-D systolic finite strains. Furthermore, the cumulative systolic radial and circumferential strains, which were angle- and frame-rate independent, were obtained from the 2-D finite-strain components and imaged in full view to detect the ischemic region. We also explored the theoretical understanding of the limitations of our technique for the accurate depiction of disease and validated it in vivo against tagged magnetic resonance imaging (tMRI) in the case of a normal human myocardium in a 2-D short-axis (SA) echocardiographic view. The theoretical framework succeeded in demonstrating that the 2-D myocardial elastography technique was a reliable tool for the complete estimation and depiction of the in-plane myocardial deformation field as well as for accurate identification of pathological mechanical function using established finite-element, left-ventricular canine models. In a preliminary study, the 2-D myocardial elastography was shown capable of imaging myocardial deformation comparable to equivalent tMRI estimates in a clinical setting.