Point Spread Functions and Deconvolution of Ultrasonic Images

被引:41
作者
Dalitz, Christoph [1 ]
Pohle-Froehlich, Regina [2 ]
Michalk, Thorsten [3 ]
机构
[1] Niederrhein Univ Appl Sci, Inst Pattern Recognit, Math & Comp Sci, Krefeld, Germany
[2] Niederrhein Univ Appl Sci, Inst Pattern Recognit, Krefeld, Germany
[3] GE Sensing & Inspect Technol GmbH, Hurth, Germany
关键词
WAVELET REGULARIZED DECONVOLUTION; SOUND FIELD; PROPAGATION; DISK;
D O I
10.1109/TUFFC.2014.006717
中图分类号
O42 [声学];
学科分类号
070206 ; 082403 ;
摘要
This article investigates the restoration of ultrasonic pulse-echo C-scan images by means of deconvolution with a point spread function (PSF). The deconvolution concept from linear system theory (LST) is linked to the wave equation formulation of the imaging process, and an analytic formula for the PSF of planar transducers is derived. For this analytic expression, different numerical and analytic approximation schemes for evaluating the PSF are presented. By comparing simulated images with measured C-scan images, we demonstrate that the assumptions of LST in combination with our formula for the PSF are a good model for the pulse-echo imaging process. To reconstruct the object from a C-scan image, we compare different deconvolution schemes: the Wiener filter, the ForWaRD algorithm, and the Richardson-Lucy algorithm. The best results are obtained with the Richardson-Lucy algorithm with total variation regularization. For distances greater or equal twice the near field distance, our experiments show that the numerically computed PSF can be replaced with a simple closed analytic term based on a far field approximation.
引用
收藏
页码:531 / 544
页数:14
相关论文
共 39 条
[1]  
Abramowitz M., 1984, Pocketbook of Mathematical Functions
[2]   Recursive least squares adaptive filters for ultrasonic signal deconvolution [J].
Alessandrini, M. ;
De Marchi, L. ;
Speciale, N. .
PROCEEDINGS OF 2008 IEEE INTERNATIONAL SYMPOSIUM ON CIRCUITS AND SYSTEMS, VOLS 1-10, 2008, :2937-2940
[3]  
[Anonymous], 5272 INRI
[4]  
[Anonymous], 10 INT S EL TEL ISET
[5]  
[Anonymous], SIGNAL PROCESSING MA
[6]  
[Anonymous], 2001, Radiography, DOI [10.1053/radi.2001.0341, DOI 10.1053/RADI.2001.0341]
[7]   MEANS FOR COMPUTING THE KIRCHHOFF SURFACE INTEGRAL FOR A DISK RADIATOR AS A SINGLE INTEGRAL WITH FIXED LIMITS [J].
ARCHERHALL, JA ;
BASHTER, AI ;
HAZELWOOD, AJ .
JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA, 1979, 65 (06) :1568-1570
[8]   Acceleration of iterative image restoration algorithms [J].
Biggs, DSC ;
Andrews, M .
APPLIED OPTICS, 1997, 36 (08) :1766-1775
[9]   DIFFRACTION THEORY [J].
BOUWKAMP, CJ .
REPORTS ON PROGRESS IN PHYSICS, 1954, 17 :35-100
[10]   Richardson-Lucy algorithm with total variation regularization for 3D confocal microscope deconvolution [J].
Dey, N ;
Blanc-Feraud, L ;
Zimmer, C ;
Roux, P ;
Kam, Z ;
Olivo-Marin, JC ;
Zerubia, J .
MICROSCOPY RESEARCH AND TECHNIQUE, 2006, 69 (04) :260-266