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

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 ( ź , ź ) , 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.

[1]  Rohit Chandra,et al.  Parallel programming in openMP , 2000 .

[2]  J. Demmel,et al.  Sun Microsystems , 1996 .

[3]  Nelly Pustelnik,et al.  Epigraphical Projection and Proximal Tools for Solving Constrained Convex Optimization Problems: Part I , 2012, ArXiv.

[4]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[5]  Antonin Chambolle,et al.  A First-Order Primal-Dual Algorithm for Convex Problems with Applications to Imaging , 2011, Journal of Mathematical Imaging and Vision.

[6]  Barbara Chapman,et al.  Using OpenMP - portable shared memory parallel programming , 2007, Scientific and engineering computation.

[7]  R. A. Leibler,et al.  On Information and Sufficiency , 1951 .

[8]  F. J. Anscombe,et al.  THE TRANSFORMATION OF POISSON, BINOMIAL AND NEGATIVE-BINOMIAL DATA , 1948 .

[9]  Maria Ganzha,et al.  Performance analysis of a scalable algorithm for 3D linear transforms , 2014, 2014 Federated Conference on Computer Science and Information Systems.

[10]  Gabriele Steidl,et al.  Homogeneous Penalizers and Constraints in Convex Image Restoration , 2013, Journal of Mathematical Imaging and Vision.

[11]  Barbara Chapman,et al.  Using OpenMP: Portable Shared Memory Parallel Programming (Scientific and Engineering Computation) , 2007 .

[12]  Alessandro Foi,et al.  Optimal Inversion of the Anscombe Transformation in Low-Count Poisson Image Denoising , 2011, IEEE Transactions on Image Processing.

[13]  Gabriele Steidl,et al.  Epigraphical Projection for Solving Least Squares Anscombe Transformed Constrained Optimization Problems , 2013, SSVM.

[14]  W. Walker,et al.  Mpi: a Standard Message Passing Interface 1 Mpi: a Standard Message Passing Interface , 1996 .

[15]  Fionn Murtagh,et al.  Image restoration with noise suppression using a multiresolution support. , 1995 .