Three-dimensional iterative reconstruction of pulsed radiation sources using spherical harmonic decomposition

Neutron and x-ray imaging are essential ways to diagnose a pulsed radiation source. The three-dimensional (3D) intensity distribution reconstructed from two-dimensional (2D) radiation images can significantly promote research regarding the generation and variation mechanisms of pulsed radiation sources. Only a few (≤5) projected images at one moment are available due to the difficulty in building imaging systems for high-radiation-intensity and short-pulsed sources. The reconstruction of a 3D source with a minimal number of 2D images is an ill-posed problem that leads to severe structural distortions and artifacts of the image reconstructed by conventional algorithms. In this paper, we present an iterative method to reconstruct a 3D source using spherical harmonic decomposition. Our algorithm improves the representation ability of spherical harmonic decomposition for 3D sources by enlarging the order of the expansion, which is limited in current analytical reconstruction algorithms. Prior knowledge of the source can be included to obtain a reasonable solution. Numerical simulations demonstrate that the reconstructed image quality of the iterative algorithm is better than that of the analytical algorithm. The iterative method can suppress the effect of noise in the integral projection image and has better robustness and adaptability than the analytical method.