Performance analysis of a parallel algorithm for restoring large-scale CT images

详细信息    查看全文

• 作者：Stanislav Harizanov^a ; ^{sharizanov@parallel.bas.bg" class="auth_mail" title="E-mail the corresponding author} ; ; Ivan Lirkov^a ; ^{ivan@parallel.bas.bg" class="auth_mail" title="E-mail the corresponding author} ; ; Krassimir Georgiev^a ; ^{georgiev@parallel.bas.bg" class="auth_mail" title="E-mail the corresponding author} ; ; Marcin Paprzycki^b ; ^d ; ^{marcin.paprzycki@ibspan.waw.pl" class="auth_mail" title="E-mail the corresponding author} ; ; Maria Ganzha^b ; ^c ; ^{maria.ganzha@ibspan.waw.pl" class="auth_mail" title="E-mail the corresponding author} ;
• 关键词：68U10 ; 94A08 ; 68W10 ; 65Y20 ; 65D18 ; 65F50
• 刊名：Journal of Computational and Applied Mathematics
• 出版年：2017
• 出版时间：15 January 2017
• 年：2017
• 卷：310
• 期：Complete
• 页码：104-114
• 全文大小：455 K

文摘

In multiple areas of image processing, such as Computed Tomography, in which data acquisition is based on counting particles that hit a detector surface, Poisson noise occurs. Using variance-stabilizing transformations, the Poisson noise can be approximated by a Gaussian one, for which classical denoising filters can be used. This paper presents an experimental performance study of a parallel implementation of the Poissonian image restoration algorithm, introduced in Harizanov et al. (2013). Hybrid parallelization based on MPI and OpenMP standards is investigated. The convergence rate of the algorithm heavily depends on both the image size and the choice of input parameters (ρ,σ)$(\rho ,\sigma )$, thus maximizing its parallel efficiency is vital for real-life applications. The implementation is tested for high-resolution radiographic images, on Linux clusters with Intel processors and on an IBM supercomputer.