Restoration of medical ultrasound images using two-dimensional homomorphic deconvolution

Describes how two-dimensional (2D) homomorphic deconvolution can be used to improve the lateral and radial resolution of medical ultrasound images recorded by a sector scanner. The recorded radio frequency ultrasound image in polar coordinates is considered as a 2D sequence of angle and depth convolved with a 2D space invariant point-spread function (PSF). Each polar coordinate sequence is transformed into the 2D complex cepstrum domain using the fast Fourier transform for Cartesian coordinates. The low-angle and low-depth portion of this sequence is taken as an estimate of the complex cepstrum representation of the PSF. It is transformed back to the Fourier frequency domain and is used to compute the deconvolved angle and depth sequence by 2D Wiener filtering. Two-dimensional homomorphic deconvolution produced substantial improvement in the resolution of B-mode images of a tissue-mimicking phantom in vitro and of several human tissues in vivo. It was better than lateral or radial homomorphic deconvolution alone, and better than 2D Wiener filtering with a PSF recorded in vitro.<>