Abstract-The point spread function (PSF), namely the response of an ultrasound system to a point source, is a powerful measure of the quality of an imaging system. The lack of an analytical formulation inhibits many applications ranging from apodization optimization, array-design, and deconvolution algorithms. We propose to fill this gap through a general PSF derivation that is flexible with respect to the type of transmission (synthetic aperture, plane-wave, diverging-wave etc.), while faithfully capturing the spatially-variant blurring of the Tissue Reflectivity Function as caused by Delay-And-Sum reconstruction. We validate the derived PSF against simulation using Field II, and show that accounting for PSF spatial-variability in sparsebased deconvolution improves reconstruction.Index Terms-Point-Spread-Function, deconvolution, image enhancement, ultrasonic image simulation, apodization design.
I. BACKGROUND AND MOTIVATIONIn ultrasound (US) imaging, the finite bandwidth and aperture of transducer elements limits image resolution. This limitation on the resolving capability of the tissue reflectivity function (TRF) can be modelled explicitly by re-casting the radio-frequency (RF) image as the results of an operator between the point-spread function (PSF) and the TRF. In this sense, the PSF, defined as the response of the imaging method in presence of a single scatterer, contains the blurring due to the instrument, and is a powerful tool in assessing the equipment performance in terms of imaging quality.Deconvolution methods use RF images to retrieve the TRF by accounting for the effect of PSF and prior knowledge of the image type. Many authors thus discuss deconvolution in conjunction with PSF assessment [1]- [3]. Two approaches to account for PSF blurring can be distinguished: deterministic Most deconvolution methods, blind or not, assume a spatially-invariant PSF model, primarily for computational purposes. In the blind-deconvolution context, [3], [5] argue that tissue-dependent attenuation and dispersive effects require the PSF to be estimated during the TRF computation. To avoid too complex an optimization, the PSF is usually assumed spatially-invariant across the imaging domain. In non-blind deconvolution [4], [6], a deterministic model of the PSF is obtained by means of simulation, such as Field II [7] or numerical approximation [1]. Again, the convenience of a spatiallyinvariant PSF avoids time-consuming repeated simulation.In both cases, a spatial invariance assumption on the PSF hence leads to important computational gain, since the operator linking the TRF and RF image becomes a convolution, which can efficiently be implemented in the Fourier domain. This assumption can however have profound consequences on the quality of the recovered TRF images: in practice, the PSF can indeed vary quite dramatically across the imaging domain, leading to non-uniform recovery performances.A few attempts have been made to account for spatiallyvarying PSFs [3]. However, they usually come at the cost of some simplifying...