Fast non-blind deconvolution based on 2D point spread function database for real-time ultrasound imaging

In the ultrasound medical imaging system, blurring which occurs after passing through ultrasound scanner system, represents Point Spread Function (PSF) that describes the response of the ultrasound imaging system to a point source distribution. So, de-blurring can be achieved by de-convolving the images with an estimated of PSF. However, it is hard to attain an accurate estimation of PSF due to the unknown properties of the tissues of the human body through the ultrasound signal propagates. In addition to, the complexity is very high in order to estimate point spread function and de-convolve the ultrasound image with estimated PSF for real-time implementation of ultrasound imaging. Therefore, conventional methods of ultrasound image restoration are based on a simple 1D PSF estimation [8] that axial direction only by restoring the performance improvement is not in the direction of Lateral. And, in case of 2D PSF estimation, PSF estimation and restoration of the high complexity is not being widely used. In this paper, we proposed new method for selection of the 2D PSF (estimated PSF of the average speed sound and depth) simultaneously with performing fast non-blind 2D de-convolution in the ultrasound imaging system. Our algorithm works on the beam-formed uncompressed radio-frequency data, with pre-measured and estimated 2D PSFs database from actual probe used. In the 2d PSF database, there are pre-measured and estimated 2D PSFs that classified the each different depth (about 5 different depths) and speed of sound (about 1450 or 1540m/s). Using a minimum variance and simple Weiner filter method, we present a novel way to select the optimal 2D PSF in pre-measured and estimated 2D PSFs database that acquired from the actual transducer being used. For de-convolution part with the chosen PSF, we focused on the low complexity issue. So, we are using the Weiner Filter and fast de-convolution technique using hyper-Laplacian priors [11], [12] which is several orders of magnitude faster than existing techniques that use hyper-Laplacian priors. Then, in order to prevent discontinuities between the differently restored each depth image regions, we use the piecewise linear interpolation on overlapping regions. We have tested our algorithm with vera-sonic system and commercial ultrasound scanner (Philips C4-2), in known speed of sound phantoms and unknown speeds in vivo scans. We have applied a non-blind de-convolution with 2D PSFs database for ultrasound imaging system. Using the real PSF from actual transducer being used, our algorithm produces a better restoration of ultrasound image than de-convolution by simulated PSF, and has low complexity for real-time ultrasound imaging. This method is robust and easy to implement. This method may be a realistic candidate for real-time implementation.

[1]  T. Taxt Comparison of cepstrum-based methods for radial blind deconvolution of ultrasound images , 1997, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[2]  T. Taxt,et al.  Two-dimensional noise-robust blind deconvolution of ultrasound images , 2001, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[3]  Aggelos K. Katsaggelos,et al.  Digital image restoration , 2012, IEEE Signal Process. Mag..

[4]  R. Jirik,et al.  Two-dimensional blind Bayesian deconvolution of medical ultrasound images , 2008, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[5]  Rob Fergus,et al.  Fast Image Deconvolution using Hyper-Laplacian Priors , 2009, NIPS.

[6]  Rob Fergus,et al.  Blind deconvolution using a normalized sparsity measure , 2011, CVPR 2011.

[7]  J. Arendt Paper presented at the 10th Nordic-Baltic Conference on Biomedical Imaging: Field: A Program for Simulating Ultrasound Systems , 1996 .

[8]  Graham M. Treece,et al.  Estimation of average speed of sound using deconvolution of medical ultrasound data. , 2010, Ultrasound in medicine & biology.

[9]  Cishen Zhang,et al.  A blind deconvolution approach to ultrasound imaging , 2012, IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control.

[10]  Murali Subbarao,et al.  Focused image recovery from two defocused images recorded with different camera settings , 1995, IEEE Transactions on Image Processing.

[11]  João C Machado,et al.  Ultrasonic wave speed measurement using the time-delay profile of rf-backscattered signals: simulation and experimental results. , 2002, The Journal of the Acoustical Society of America.

[12]  James Ng,et al.  Sensitivity to point-spread function parameters in medical ultrasound image deconvolution. , 2009, Ultrasonics.

[13]  Oleg V. Michailovich,et al.  Blind Deconvolution of Medical Ultrasound Images: A Parametric Inverse Filtering Approach , 2007, IEEE Transactions on Image Processing.

[14]  Bernard Heit,et al.  Monocular depth perception by evaluation of the blur in defocused images , 1994, Proceedings of 1st International Conference on Image Processing.

[15]  J A Jensen,et al.  Deconvolution of Ultrasound Images , 1992, Ultrasonic imaging.

[16]  Torfinn Taxt,et al.  High-resolution ultrasonic imaging using fast two-dimensional homomorphic filtering. , 2006, IEEE transactions on ultrasonics, ferroelectrics, and frequency control.

[17]  J-F Synnevåg,et al.  A low-complexity data-dependent beamformer , 2011, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.