Michael F. Insana

Phantom materials for elastography

Acoustic and mechanical properties are reported for gelatin materials used to construct tissue-like phantoms for elasticity imaging (elastography). A device and procedure for measuring elastic properties are described. The measured compression forces were comparable to results obtained from finite element analysis when linear elastic media are assumed. Also measured were the stress relaxation, temporal stability, and melting point of the materials. Aldehyde concentration was used to increase the stiffness of the gelatin by controlling the amount of collagen cross-linking. A broad range of tissue-like elastic properties was achieved with these materials, although gels continued to stiffen for several weeks. The precision for elastic modulus measurements ranged from less than 0.1% for 100 kPa samples to 8.9% for soft (<10 kPa), sticky samples.

Describing small‐scale structure in random media using pulse‐echo ultrasound

A method for estimating structural properties of random media is described. The size, number density, and scattering strength of particles are estimated from an analysis of the radio frequency (rf) echo signal power spectrum. Simple correlation functions and the accurate scattering theory of Faran [J.J. Faran, J. Acoust. Soc. Am. 23, 405-418 (1951)], which includes the effects of shear waves, were used separately to model backscatter from spherical particles and thereby describe the structures of the medium. These methods were tested using both glass sphere-in-agar and polystyrene sphere-in-agar scattering media. With the appropriate correlation function, it was possible to measure glass sphere diameters with an accuracy of 20%. It was not possible to accurately estimate the size of polystyrene spheres with the simple spherical and Gaussian correlation models examined because of a significant shear wave contribution. Using the Faran scattering theory for spheres, however, the accuracy for estimating diameters was improved to 10% for both glass and polystyrene scattering media. It was possible to estimate the product of the average scattering particle number density and the average scattering strength per particle, but with lower accuracy than the size estimates. The dependence of the measurement accuracy on the inclusion of shear waves, the wavelength of sound, and medium attenuation are considered, and the implications for describing the structure of biological soft tissues are discussed.

Statistical properties of radio-frequency and envelope-detected signals with applications to medical ultrasound

Both radio-frequency (rf) and envelope-detected signal analyses have lead to successful tissue discrimination in medical ultrasound. The extrapolation from tissue discrimination to a description of the tissue structure requires an analysis of the statistics of complex signals. To that end, first- and second-order statistics of complex random signals are reviewed, and an example is taken from rf signal analysis of the backscattered echoes from diffuse scatterers. In this case the scattering form factor of small scatterers can be easily separated from long-range structure and corrected for the transducer characteristics, thereby yielding an instrument-independent tissue signature. The statistics of the more economical envelope- and square-law-detected signals are derived next and found to be almost identical when normalized autocorrelation functions are used. Of the two nonlinear methods of detection, the square-law or intensity scheme gives rise to statistics that are more transparent to physical insight. Moreover, an analysis of the intensity-correlation structure indicates that the contributions to the total echo signal from the diffuse scatter and from the steady and variable components of coherent scatter can still be separated and used for tissue characterization. However, this analysis is not system independent. Finally, the statistical methods of this paper may be applied directly to envelope signals in nuclear-magnetic-resonance imaging because of the approximate equivalence of second-order statistics for magnitude and intensity.

2-D companding for noise reduction in strain imaging

Companding is a signal preprocessing technique for improving the precision of correlation-based time delay measurements. In strain imaging, companding is applied to warp 2-D or 3-D ultrasonic echo fields to improve coherence between data acquired before and after compression. It minimizes decorrelation errors, which are the dominant source of strain image noise. The word refers to a spatially variable signal scaling that compresses and expands waveforms acquired in an ultrasonic scan plane or volume. Temporal stretching by the applied strain is a single-scale (global), 1-D companding process that has been used successfully to reduce strain noise. This paper describes a two-scale (global and local), 2-D companding technique that is based on a sum-absolute-difference (SAD) algorithm for blood velocity estimation. Several experiments are presented that demonstrate improvements in target visibility for strain imaging. The results show that, if tissue motion can be confined to the scan plane of a linear array transducer, displacement variance can be reduced two orders of magnitude using 2-D local companding relative to temporal stretching.

Analysis Of Ultrasound Image Texture Via Generalized Rician Statistics

Tissue signatures are obtained from the first- and second-order statistics of ultrasonic B-scan texture. Laboratory measurements and early clinical results show that the image may be segmented to discriminate between different normal tissues and to detect abnormal conditions based on a multidimensional feature space. These features describe the intrinsic backscatter properties of the tissues imaged and may be used as the basis of an automatic tissue characterization algorithm.

Fundamental correlation lengths of coherent speckle in medical ultrasonic images

Refinements to previous analyses of the natural correlation lengths within simple images and between images to be compounded are presented. Comparison of theoretical and experimental results show very good agreement for the case of Rayleigh scattering media: the correlation length within a simple image is comparable to the resolution cell size; the correlation length between images to be spatially compounded is comparable to, but smaller than, the transducer on array aperture; and the correlation length between images to be frequency-compounded becomes a frequency comparable to their bandwidth. Complications arising from the presence of specular scattering or due to the presence of just a few scatterers are considered. It is shown that straightforward solutions exist for either of these problems taken by itself. When they occur in combination, calibration techniques may lead to unambiguous identification of the contributions to the scattering from diffuse or incoherent scattering and from specular or coherent scattering, and to estimation of the density of diffuse scatterers.<<ETX>>

Parametric ultrasound imaging from backscatter coefficient measurements: Image formation and interpretation

A broadband method for measuring backscatter coefficients σ b and other acoustic parameters is described. From the σ b measurements, using a commercially-available imaging system, four high-resolution parametric ultrasound images are formed in a C-scan image plane. Scatterer size images are computed from the frequency dependence of σ b and a correlation model function that describes the structure and elastic properties of the medium. Scattering strength images are computed from the absolute magnitude of σ b . Chi-square images are generated to display how well the correlation model represents the interrogated medium. Integrated backscatter coefficient images are formed over the transducer bandwidth. All four images are generated simultaneously and compared with the corresponding B-mode image. Test samples with known physical properties were used to demonstrate experimentally that accurate parametric images are possible if an accurate correlation model is used. Local variations in attenuation, the center frequency and bandwidth of the transducer, and the distribution of scatterer sizes greatly influence the accuracy of estimates and the appearance of the image, thus demonstrating the importance of these factors in parametric image interpretation.

Ultrasound elastography based on multiscale estimations of regularized displacement fields

Elasticity imaging is based on the measurements of local tissue deformation. The approach to ultrasound elasticity imaging presented in this paper relies on the estimation of dense displacement fields by a coarse-to-fine minimization of an energy function that combines constraints of conservation of echo amplitude and displacement field continuity. The multiscale optimization scheme presents several characteristics aimed at improving and accelerating the convergence of the minimization process. This includes the nonregularized initialization at the coarsest resolution and the use of adaptive configuration spaces. Parameters of the energy model and optimization were adjusted using data obtained from a tissue-like phantom material. Elasticity images from normal in vivo breast tissue were subsequently obtained with these parameters. Introducing a smoothness constraint into motion field estimation helped solve ambiguities due to incoherent motion, leading to elastograms less degraded by decorrelation noise than the ones obtained from correlation-based techniques.

Application of autoregressive spectral analysis to cepstral estimation of mean scatterer spacing

The problem of estimation of mean scatterer spacing in an object containing regularly spaced structures is addressed. An autoregressive (AR) spectral estimation method is compared with a conventional fast Fourier transform (FFT)-based approach for this task. Regularly spaced structures produce a periodicity in the power spectrum of ultrasonic backscatter. This periodicity is manifested as a peak in the cepstrum. A phantom was constructed for comparison of the two methods. It contained regularly spaced nylon filaments. It also contained randomly positioned glass spheres that produced incoherent backscatter. In an experiment in which this target was interrogated using broadband ultrasound, the AR spectral estimate offered considerable improvement over the FFT when the analysis gate length was on the order of the structural dimension. Advantages included improved resolution, reduction in bias and variance of scatterer spacing estimates, and greater resistance to ringing artifacts. Data were also acquired from human liver in vivo. AR spectral estimates on human data exhibited a decreased dependence on gate length. These results offer promise for enhanced spatial resolution and accuracy in ultrasonic tissue characterization and nondestructive evaluation of materials.<<ETX>>

Unified approach to the detection and classification of speckle texture in diagnostic ultrasound

Second order statistics have been derived for the speckle in diagnostic ultrasound that arises from diffuse (incoherent) scattering in the presence of distributed and organized specular (coherent) scattering. They serve as the basis for a three-dimensional feature space in which tissue textures can be classified. The covariance matrix of the measurements in this space is a generalization of the speckle spot number or sampling concept that arises in the study of signal or lesion detectability.