Comparison of temporal and spectral scattering methods using acoustically large breast models derived from magnetic resonance images.

Accurate and efficient modeling of ultrasound propagation through realistic tissue models is important to many aspects of clinical ultrasound imaging. Simplified problems with known solutions are often used to study and validate numerical methods. Greater confidence in a time-domain k-space method and a frequency-domain fast multipole method is established in this paper by analyzing results for realistic models of the human breast.

Aberration compensation of an ultrasound imaging instrument with a reduced number of channels.

Focusing and imaging qualities of an ultrasound imaging system that uses aberration correction were experimentally investigated as functions of the number of parallel channels. Front-end electronics that consolidate signals from multiple physical elements can be used to lower hardware and computational costs by reducing the number of parallel channels. However, the signals from sparse arrays of synthetic elements yield poorer aberration estimates.

Reduced-Rank Approximations to the Far-Field Transform in the Gridded Fast Multipole Method.

The fast multipole method (FMM) has been shown to have a reduced computational dependence on the size of finest-level groups of elements when the elements are positioned on a regular grid and FFT convolution is used to represent neighboring interactions. However, transformations between plane-wave expansions used for FMM interactions and pressure distributions used for neighboring interactions remain significant contributors to the cost of FMM computations when finest-level groups are large. The transformation operators, which are forward and inverse Fourier transforms with the wave space confined to the unit sphere, are smooth and well approximated using reduced-rank decompositions that further reduce the computational dependence of the FMM on finest-level group size.

Acoustic scattering by arbitrary distributions of disjoint, homogeneous cylinders or spheres.

A T-matrix formulation is presented to compute acoustic scattering from arbitrary, disjoint distributions of cylinders or spheres, each with arbitrary, uniform acoustic properties. The generalized approach exploits the similarities in these scattering problems to present a single system of equations that is easily specialized to cylindrical or spherical scatterers. By employing field expansions based on orthogonal harmonic functions, continuity of pressure and normal particle velocity are directly enforced at each scatterer using diagonal, analytic expressions to eliminate the need for integral equations.

The Fast Multipole Method and Fourier Convolution for the Solution of Acoustic Scattering on Regular Volumetric Grids.

The fast multipole method (FMM) is applied to the solution of large-scale, three-dimensional acoustic scattering problems involving inhomogeneous objects defined on a regular grid. The grid arrangement is especially well suited to applications in which the scattering geometry is not known a priori and is reconstructed on a regular grid using iterative inverse scattering algorithms or other imaging techniques. The regular structure of unknown scattering elements facilitates a dramatic reduction in the amount of storage and computation required for the FMM, both of which scale linearly with the number of scattering elements.

Large-scale propagation of ultrasound in a 3-D breast model based on high-resolution MRI data.

A 40 x 35 x 25-mm(3) specimen of human breast consisting mostly of fat and connective tissue was imaged using a 3-T magnetic resonance scanner. The resolutions in the image plane and in the orthogonal direction were 130 microm and 150 microm, respectively. Initial processing to prepare the data for segmentation consisted of contrast inversion, interpolation, and noise reduction.

A mesh-free approach to acoustic scattering from multiple spheres nested inside a large sphere by using diagonal translation operators.

A multiple-scattering approach is presented to compute the solution of the Helmholtz equation when a number of spherical scatterers are nested in the interior of an acoustically large enclosing sphere. The solution is represented in terms of partial-wave expansions, and a linear system of equations is derived to enforce continuity of pressure and normal particle velocity across all material interfaces. This approach yields high-order accuracy and avoids some of the difficulties encountered when using integral equations that apply to surfaces of arbitrary shape.

A model of distributed phase aberration for deblurring phase estimated from scattering.

Correction of aberration in ultrasound imaging uses the response of a point reflector or its equivalent to characterize the aberration. Because a point reflector is usually unavailable, its equivalent is obtained using statistical methods, such as processing reflections from multiple focal regions in a random medium. However, the validity of methods that use reflections from multiple points is limited to isoplanatic patches for which the aberration is essentially the same.

A k-space method for acoustic propagation using coupled first-order equations in three dimensions.

A previously described two-dimensional k-space method for large-scale calculation of acoustic wave propagation in tissues is extended to three dimensions. The three-dimensional method contains all of the two-dimensional method features that allow accurate and stable calculation of propagation. These features are spectral calculation of spatial derivatives, temporal correction that produces exact propagation in a homogeneous medium, staggered spatial and temporal grids, and a perfectly matched boundary layer.

Scattering calculation and image reconstruction using elevation-focused beams.

Pressure scattered by cylindrical and spherical objects with elevation-focused illumination and reception has been analytically calculated, and corresponding cross sections have been reconstructed with a two-dimensional algorithm. Elevation focusing was used to elucidate constraints on quantitative imaging of three-dimensional objects with two-dimensional algorithms. Focused illumination and reception are represented by angular spectra of plane waves that were efficiently computed using a Fourier interpolation method to maintain the same angles for all temporal frequencies.

Born iterative reconstruction using perturbed-phase field estimates.

A method of image reconstruction from scattering measurements for use in ultrasonic imaging is presented. The method employs distorted-wave Born iteration but does not require using a forward-problem solver or solving large systems of equations. These calculations are avoided by limiting intermediate estimates of medium variations to smooth functions in which the propagated fields can be approximated by phase perturbations derived from variations in a geometric path along rays.

An eigenfunction method for reconstruction of large-scale and high-contrast objects.

A multiple-frequency inverse scattering method that uses eigenfunctions of a scattering operator is extended to image large-scale and high-contrast objects. The extension uses an estimate of the scattering object to form the difference between the scattering by the object and the scattering by the estimate of the object. The scattering potential defined by this difference is expanded in a basis of products of acoustic fields.

A ring transducer system for medical ultrasound research.

An ultrasonic ring transducer system has been developed for experimental studies of scattering and imaging. The transducer consists of 2048 rectangular elements with a 2.5-MHz center frequency, a 67% -6 dB bandwidth, and a 0.

Reduction of variance in spectral estimates for correction of ultrasonic aberration.

A variance reduction factor is defined to describe the rate of convergence and accuracy of spectra estimated from overlapping ultrasonic scattering volumes when the scattering is from a spatially uncorrelated medium. Assuming that the individual volumes are localized by a spherically symmetric Gaussian window and that centers of the volumes are located on orbits of an icosahedral rotation group, the factor is minimized by adjusting the weight and radius of each orbit. Conditions necessary for the application of the variance reduction method, particularly for statistical estimation of aberration, are examined.

Statistical estimation of ultrasonic propagation path parameters for aberration correction.

Parameters in a linear filter model for ultrasonic propagation are found using statistical estimation. The model uses an inhomogeneous-medium Green's function that is decomposed into a homogeneous-transmission term and a path-dependent aberration term. Power and cross-power spectra of random-medium scattering are estimated over the frequency band of the transmit-receive system by using closely situated scattering volumes.

Comparisons of lesion detectability in ultrasound images acquired using time-shift compensation and spatial compounding.

The effects of aberration, time-shift compensation, and spatial compounding on the discrimination of positive-contrast lesions in ultrasound b-scan images are investigated using a two-dimensional (2-D) array system and tissue-mimicking phantoms. Images were acquired within an 8.8 x 12-mm2 field of view centered on one of four statistically similar 4-mm diameter spherical lesions.

Spectral estimation for characterization of acoustic aberration.

Spectral estimation based on acoustic backscatter from a motionless stochastic medium is described for characterization of aberration in ultrasonic imaging. The underlying assumptions for the estimation are: The correlation length of the medium is short compared to the length of the transmitted acoustic pulse, an isoplanatic region of sufficient size exists around the focal point, and the backscatter can be modeled as an ergodic stochastic process. The motivation for this work is ultrasonic imaging with aberration correction.

Simulation of ultrasonic focus aberration and correction through human tissue.

Ultrasonic focusing in two dimensions has been investigated by calculating the propagation of ultrasonic pulses through cross-sectional models of human abdominal wall and breast. Propagation calculations used a full-wave k-space method that accounts for spatial variations in density, sound speed, and frequency-dependent absorption and includes perfectly matched layer absorbing boundary conditions. To obtain a distorted receive wavefront, propagation from a point source through the tissue path was computed.

Spatial coherence analysis applied to aberration correction using a two-dimensional array system.

Complex degree of coherence functions are computed using synthetic and measured ultrasound data to demonstrate noteworthy aspects of coherence analysis in the context of aberration correction. Coherence functions calculated from synthetic data illustrate the importance of proper normalization of the constituent cross-correlation integrals when weak elements and receiver directivity are significant factors. The synthetic data also show that a spike can occur at the zero-lag position of the coherence function when the signal-to-noise ratio is reduced by element directivity near the edges of a large aperture.

Examples of design curves for multirow arrays used with time-shift compensation.

The effect of 1.75-D array element height on time-shift compensation is analyzed using a two-dimensional array system and a distributed aberration phantom. The analysis uses hydrophone measurements of transmit beams.

Estimation and compensation of ultrasonic wavefront distortion using a blind system identification method.

A common random input filter model is described for estimation and correction of wavefront aberration in ultrasonic b-scan imaging. In the model, aberration between the focus and the transducer elements is represented by the response of a linear filter bank to a common random signal. The response of each filter in the bank is found using a two-level extension of an existing subspace method for blind system identification.

A k-space method for coupled first-order acoustic propagation equations.

A k-space method for large-scale simulation of ultrasonic pulse propagation is presented. The present method, which solves the coupled first-order differential equations for wave propagation in inhomogeneous media, is derived in a simple form analogous to previous finite-difference methods with staggered spatial and temporal grids. Like k-space methods based on second-order wave equations, the present method is exact for homogeneous media, unconditionally stable for "slow" [c(r) < or = c0] media, and highly accurate for general weakly scattering media.