Assisted-freehand ultrasound elasticity imaging

Good-quality elasticity imaging requires highly controlled compressions of the breast, which are often challenging to obtain with freehand, even by an experienced radiologist. This paper presents assisted-freehand ultrasound (AFUSON): a fusion of freehand and automated ultrasound systems designed to assisted elasticity imaging acquisition while remaining as flexible as freehand. In the form of a hand-held device, this semi-automatic solution delivers both increased acquisition precision and control. Compared with freehand acquisitions, it reduces out-of-plane motion decorrelation by one-half and lateral motion by one-third, increases within-scan repeatability by 50%, and does so across operators.     

An autocorrelation-based method for improvement of sub-pixel displacement estimation in ultrasound strain imaging

In ultrasound strain and elasticity imaging, an accurate and cost-effective sub-pixel displacement estimator is required because strain/elasticity imaging quality relies on the displacement SNR, which can often be higher if more computational resources are provided. In this paper, we introduce an autocorrelation-based method to cost-effectively improve subpixel displacement estimation quality. To quantitatively evaluate the performance of the autocorrelation method, simulated and tissue-mimicking phantom experiments were performed. The computational cost of the autocorrelation method is also discussed. The results of our study suggest the autocorrelation method can be used for a real-time elasticity imaging system.     

GPU-based real-time small displacement estimation with ultrasound

General purpose computing on graphics processing units (GPUs) has been previously shown to speed up computationally intensive data processing and image reconstruction algorithms for computed tomography (CT), magnetic resonance (MR), and ultrasound images. Although some algorithms in ultrasound have been converted to GPU processing, many investigative ultrasound research systems still use serial processing on a single CPU. One such ultrasound modality is acoustic radiation force impulse (ARFI) imaging, which investigates the mechanical properties of soft tissue. Traditionally, the raw data are processed offline to estimate the displacement of the tissue after the application of radiation force. It is highly advantageous to process the data in real-time to assess their quality and make modifications during a study. In this paper, we present algorithms for efficient GPU parallel processing of two widely used tools in ultrasound: cubic spline interpolation and Loupas' two-dimensional autocorrelator for displacement estimation. It is shown that a commercially available graphics card can be used for these computations, achieving speed increases up to 40x compared with single CPU processing. Thus, we conclude that the GPU-based data processing approach facilitates real-time (i.e., <1 second) display of ARFI data and is a promising approach for ultrasonic research systems.     

High-resolution elasticity imaging for tissue engineering

An elasticity microscope provides high resolution images of tissue elasticity. With this instrument, it may be possible to monitor cell growth and tissue development in tissue engineering. To test this hypothesis, elasticity micrographs were obtained in two model systems commonly used for tissue engineering. In the first, strain images of a tissue-engineered smooth muscle sample clearly identified a several hundred micron thick cell layer from its supporting matrix. Because a one-dimensional mechanical model was appropriate for this system, strain images alone were sufficient to image the elastic properties. In contrast, a second system was investigated in which a simple one-dimensional mechanical model was inadequate. Uncultured collagen microspheres embedded in an otherwise homogeneous gel were imaged with the elasticity microscope. Strain images alone did not clearly depict the elastic properties of the hard spherical cell carriers. However, reconstructed elasticity images could differentiate the hard inclusion from the background gel. These results strongly suggest that the elasticity microscope may be a valuable tool for tissue engineering and other applications requiring the elastic properties of soft tissue at high spatial resolution (75 mum or less).     

Characterization of PVA cryogel for intravascular ultrasound elasticity imaging

The present study characterizes the mechanical properties of polyvinyl alcohol (PVA) cryogel in order to show its utility for intravascular elastography. PVA cryogel becomes harder with an increasing number of freeze-thaw cycles, and Young's modulus and Poisson's ratio are measured for seven samples. Mechanical tests were performed on cylindrical samples with a pressure column and on a hollow cylinder with the calculation of an intravascular elastogram. An image of the Young's modulus was obtained from the elastogram using cylinder geometry properties. Results show the mechanical similitude of PVA cryogel with the biological tissues present in arteries. A good agreement between Young's modulus obtained from pressure column and from elastogram was also observed.     

Lateral displacement estimation using tissue incompressibility

Using the incompressibility property of soft tissue, lateral displacements can be reconstructed from axial strain measurements. Results of simulations and experiments on gelatin-based tissue equivalent phantoms are compared with theoretical displacements, as well as estimates derived from traditional speckle tracking. Incompressibility processing greatly improves the accuracy and signal-to-noise ratio (SNR) of lateral displacement measurements compared with more traditional speckle tracking     

Slow light for deep tissue imaging with ultrasound modulation

Slow light has been extensively studied for applications ranging from optical delay lines to single photon quantum storage. Here, we show that the time delay of slow-light significantly improves the performance of the narrowband spectral filters needed to optically detect ultrasound from deep inside highly scattering tissue. We demonstrate this capability with a 9?cm thick tissue phantom, having 10?cm(-1) reduced scattering coefficient, and achieve an unprecedented background-free signal. Based on the data, we project real time imaging at video rates in even thicker phantoms and possibly deep enough into real tissue for clinical applications like early cancer detection.     

Echo decorrelation from displacement gradients in elasticity and velocity estimation

Several ultrasonic techniques for the estimation of blood velocity, tissue motion and elasticity are based on the estimation of displacement through echo time-delay analysis. A common assumption is that tissue displacement is constant within a short observation time used for time delay estimation (TDE). The precision of TDE is mainly limited by noise sources corrupting the echo signals. In addition to electronic and quantization noise, a substantial source of TDE error is the decorrelation of echo signals because of displacement gradients within the observation time. The authors present a theoretical model that describes the mean changes of the crosscorrelation function as a function of observation time and displacement gradient. The gradient is assumed to be small and uniform within the observation time; the decorrelation introduced by the lateral and elevational displacement components is assumed to be small compared with the decorrelation caused by the axial component. The decorrelation model predicts that the expected value of the crosscorrelation function is a low-pass filtered version of the autocorrelation function (i.e., the crosscorrelation obtained without gradients). The filter is a function of the axial gradient and the observation time. This theoretical finding is corroborated experimentally. Limitations imposed by decorrelation in displacement estimation and potential uses of decorrelation in medical ultrasound are discussed.     

Nonlinear estimation of the lateral displacement using tissue incompressibility

Using the incompressibility property of soft tissue, high quality lateral displacement distributions can be reconstructed from accurate axial displacement measurements and noisy lateral displacement estimates. Previous methods appropriate for small deformations have been extended for high magnitude deformations requiring a nonlinear model. Problems arising in incompressibility processing for large deformations are considered. Applications of nonlinear incompressibility methods to ultrasonic measurements on gel-based, tissue equivalent phantoms are given. Lateral displacement images reconstructed with nonlinear methods are compared to those reconstructed with linear methods for both small and large deformations.     

Time-varying autoregressive spectral estimation for ultrasound attenuation in tissue characterization

In the field of biological tissue characterization, fundamental acoustic attenuation properties have been demonstrated to have diagnostic importance. Attenuation caused by scattering and absorption shifts the instantaneous spectrum to the lower frequencies. Due to the time-dependence of the spectrum, the attenuation phenomenon is a time-variant process. This downward shift may be evaluated either by the maximum energy frequency of the spectrum or by the center frequency. In order to improve, in strongly attenuating media, the results given by the short-time Fourier analysis and the short-time parametric analysis, we propose two approaches adapted to this time-variant process: an adaptive method and a time-varying method. Signals backscattered by an homogeneous medium of scatterers are modeled by a computer algorithm with attenuation values ranging from 1 to 5 dB/cm MHz and a 45 MHz transducer center frequency. Under these conditions, the preliminary results obtained with the proposed time-variant methods, compared with the classical short-time Fourier analysis and the short-time auto-regressive (AR) analysis, are superior in terms of standard deviation (SD) of the attenuation coefficient estimate. This study, based on nonstationary AR spectral estimation, promises encouraging perspectives for in vitro and in vivo applications both in weakly and highly attenuating media.     

Shearwave dispersion ultrasound vibrometry (SDUV) for measuring tissue elasticity and viscosity

Characterization of tissue elasticity (stiffness) and viscosity has important medical applications because these properties are closely related to pathological changes. Quantitative measurement is more suitable than qualitative measurement (i.e., mapping with a relative scale) of tissue viscoelasticity for diagnosis of diffuse diseases where abnormality is not confined to a local region and there is no normal background tissue to provide contrast. Shearwave dispersion ultrasound vibrometry (SDUV) uses shear wave propagation speed measured in tissue at multiple frequencies (typically in the range of hundreds of Hertz) to solve quantitatively for both tissue elasticity and viscosity. A shear wave is stimulated within the tissue by an ultrasound push beam and monitored by a separate ultrasound detect beam. The phase difference of the shear wave between 2 locations along its propagation path is used to calculate shear wave speed within the tissue. In vitro SDUV measurements along and across bovine striated muscle fibers show results of tissue elasticity and viscosity close to literature values. An intermittent pulse sequence is developed to allow one array transducer for both push and detect function. Feasibility of this pulse sequence is demonstrated by in vivo SDUV measurements in swine liver using a dual transducer prototype simulating the operation of a single array transducer.     

人体组织弹性的超声检测和应用前景

组织的软硬或弹性反映了组织本身的特性。对于它们的检测所得信息往往有助于对疾病的诊断。本文介绍了运用超声检测人体组织弹性的原理入检测中信息处理和特征提取的几种方法，包括我们的研究工作并就其应用前景作了展望。     

Supersonic shear imaging: a new technique for soft tissue elasticity mapping

Supersonic shear imaging (SSI) is a new ultrasound-based technique for real-time visualization of soft tissue viscoelastic properties. Using ultrasonic focused beams, it is possible to remotely generate mechanical vibration sources radiating low-frequency, shear waves inside tissues. Relying on this concept, SSI proposes to create such a source and make it move at a supersonic speed. In analogy with the "sonic boom" created by a supersonic aircraft, the resulting shear waves will interfere constructively along a Mach cone, creating two intense plane shear waves. These waves propagate through the medium and are progressively distorted by tissue heterogeneities. An ultrafast scanner prototype is able to both generate this supersonic source and image (5000 frames/s) the propagation of the resulting shear waves. Using inversion algorithms, the shear elasticity of medium can be mapped quantitatively from this propagation movie. The SSI enables tissue elasticity mapping in less than 20 ms, even in strongly viscous medium like breast. Modalities such as shear compounding are implementable by tilting shear waves in different directions and improving the elasticity estimation. Results validating SSI in heterogeneous phantoms are presented. The first in vivo investigations made on healthy volunteers emphasize the potential clinical applicability of SSI for breast cancer detection.     

Intravascular ultrasound tissue harmonic imaging in vivo

Tissue harmonic imaging (THI) has been shown to increase image quality of medical ultrasound in the frequency range from 2 to 10 MHz and might, therefore, also be used to improve image quality in intravascular ultrasound (IVUS). In this study we constructed a prototype IVUS system that could operate in both fundamental frequency and second harmonic imaging modes. This system uses a conventional, continuously rotating, single-element IVUS catheter and was operated in fundamental 20 MHz, fundamental 40 MHz, and harmonic 40 MHz modes (transmit 20 MHz, receive 40 MHz). Hydrophone beam characterization measurements demonstrated the build-up of a second harmonic signal as a function of increasing pressure. Imaging experiments were conducted in both a tissue-mimicking phantom and in an atherosclerotic animal model in vivo. Acquisitions of fundamental 20 and 40 MHz and second harmonic acquisitions resulted in cross sections of the phantom and a rabbit aorta. The harmonic results of the imaging experiments showed the feasibility of intravascular THI with a conventional IVUS catheter both in a phantom and in vivo. The harmonic acquisitions also showed the potential of THI to reduce image artifacts compared to fundamental imaging.     

Acoustic attenuation estimation for soft tissue from ultrasound echo envelope peaks

A theoretical framework of a modified EP (envelope peak) method is developed. In this method, a wideband echo signal from each A line (one echo sequence) is first filtered in parallel by a bank of narrow-bandpass filters using a split-spectrum processing. The attenuation is then estimated from the EPs of each filtered signal using a narrowband technique. The combination of the split-spectrum processing with the narrowband technique enables the accuracy of the attenuation estimation to be well controlled without the precise measurements of the spectral shape and parameters of the transmitted pulses. On the other hand, the precision of the estimation is still determined by the bandwidth of the original echo signal, and is not affected by the split-spectrum processing. As a result, the modified EP method improves the accuracy of the attenuation estimation while retaining the high precision of the original EP method. Results from phantom experiments supported the theoretical analysis.     

Adaptive strain estimation using retrospective processing [medical US elasticity imaging

Because errors in displacement and strain estimates depend on the magnitude of the induced strain, the strain signal-to-noise ratio (SNR) will be a function of the applied deformation. If deformation is applied at the body surface, it is difficult during data acquisition to select a single surface displacement providing the highest strain SNR throughout the image. By applying continuous deformation and capturing data in real-time, the surface displacement providing the highest strain SNR can be selected retrospectively. A method to adaptively optimize strain SNR over the image plane using retrospective processing is presented and demonstrated with experimental results.     

Clutter filters adapted to tissue motion in ultrasound color flow imaging

The quality of ultrasound color flow images is highly dependent on sufficient attenuation of the clutter signals originating from stationary and slowly moving tissue. Without sufficient clutter rejection, the detection of low velocity blood flow will be poor, and the velocity estimates will have a large bias. In some situations, e.g., when imaging the coronary arteries or when the operator moves the probe in search for small vessels, there is considerable movement of tissue. It has been suggested that clutter rejection can be improved by mixing down the signal with an estimate of the mean frequency prior to high pass filtering. In this paper, we compare this algorithm with several other adaptive clutter filtering algorithms using both experimental data and simulations. We found that realistic accelerations of the tissue have a large effect on the clutter rejection. The best results were obtained by mixing down the signal with non-constant phase increments estimated from the signal. This adapted the filter to a possibly accelerated tissue motion and produced a significant improvement in clutter rejection     

Optimizing multicompression approaches to elasticity imaging

Breast lesion visibility in static strain imaging ultimately is noise limited. When correlation and related techniques are applied to estimate local displacements between two echo frames recorded before and after a small deformation, target contrast increases linearly with the amount of deformation applied. However, above some deformation threshold, decorrelation noise increases more than contrast such that lesion visibility is severely reduced. Multicompression methods avoid this problem by accumulating displacements from many small deformations to provide the same net increase in lesion contrast as one large deformation but with minimal decorrelation noise. Unfortunately, multicompression approaches accumulate echo noise (electronic and sampling) with each deformation step as contrast builds so that lesion visibility can be reduced again if the applied deformation increment is too small. This paper uses signal models and analysis techniques to develop multicompression strategies that minimize strain image noise. The analysis predicts that displacement variance is minimal in elastically homogeneous media when the applied strain increment is 0.0035. Predictions are verified experimentally with gelatin phantoms. For in vivo breast imaging, a strain increment as low as 0.0015 is recommended for minimum noise because of the greater elastic heterogeneity of breast tissue.     

Reconstructive elasticity imaging for large deformations

A method is presented to reconstruct the elastic modulus of soft tissue based on ultrasonic displacement and strain images for comparatively large deformations. If the average deformation is too large to be described with a linear elastic model, nonlinear displacement-strain relations must be used and the mechanical equilibrium equations must include high order spatial derivatives of the displacement. Numerical methods were developed to reduce error propagation in reconstruction algorithms, including these higher order derivatives. Problems arising with the methods, as well as results using ultrasound measurements on gel-based, tissue equivalent phantoms, are given. Comparison to reconstructions using a linear elastic model shows that equivalent image quality can be produced with algorithms appropriate for finite amplitude deformations.     

A new method for estimation of velocity vectors

The paper describes a new method for determining the velocity vector of a remotely sensed object using either sound or electromagnetic radiation. The movement of the object is determined from a field with spatial oscillations in both the axial direction of the transducer and in one or two directions transverse to the axial direction. By using a number of pulse emissions, the inter-pulse movement can be estimated and the velocity found from the estimated movement and the time between pulses. The method is based on the principle of using transverse spatial modulation for making the received signal influenced by transverse motion. Such a transverse modulation can be generated by using apodization on individual transducer array elements together with a special focusing scheme. A method for making such a field is presented along with a suitable two-dimensional velocity estimator. An implementation usable in medical ultrasound is described, and simulated results are presented. Simulation results for a flow of 1 m/s in a tube rotated in the image plane at specific angles (0, 15, 35, 55, 75, and 90 degrees) are made and characterized by the estimated mean value, estimated angle, and the standard deviation in the lateral and longitudinal direction. The average performance of the estimates for all angles is: mean velocity 0.99 m/s, longitudinal S.D. 0.015 m/s, and lateral S.D. 0.196 m/s. For flow parallel to the transducer the results are: mean velocity 0.95 m/s, angle 0.10, longitudinal S.D. 0.020 m/s, and lateral S.D. 0.172 m/s.