PARALLEL IMAGE RECONSTRUCTION USING THE MAXIMUM LIKELIHOOD METHOD USING A GRAPHICS PROCESSOR AND THE OpenGL LIBRARY

Abstract

The creation of fast parallel iterative statistical algorithms based on the use of graphics accelerators is an important and urgent task of great scientific and practical importance. An algorithm based on the method of maximizing the mathematical expectation of maximum likelihood (maximum likelihood expectation MLEM) is considered. MLEM is a numerical method for determining maximum likelihood estimates and, since its first application in the field of image reconstruction in 1982, remains one of the most popular statistical methods of image reconstruction, being the foundation for many other approaches. A new version of the MLEM parallel algorithm is proposed, which provides global convergence of the iterative algorithm. To parallelize the algorithm, the texture mapping method is used using the OpenGL graphics library. The parallel algorithm is described in as much detail as possible. Examples of several reconstructions of images of aluminum casting products are given The obtained result can be used for non-destructive testing of various industrial products, including testing of foundry products.