The non-Fourier image reconstruction method for the STIX instrument

Abstract

AbstractIn this work we aimed to develop the image reconstruction algorithm without any analytical simplifications and restrictions. In our method we abandon Fourier’s approach to image reconstruction, and instead use the number of counts recorded in each detector pixel, and then reconstruct each image using a classical Richardson-Lucy algorithm. Among similar works performed in the past, our approach is based, for the first time, on the real geometry of STIX. We made a preliminary analysis of expected differences in STIX imaging which may occur due to usage of slightly different geometries. The other difference is that we use single-pixel-response maps. Namely, knowing the instrument geometry we are able to calculate the detector response for point sources covering entire the solar disc. Next, we iteratively combine them with varying weights until the best match between reconstructed and observed detector responses is achieved. Preliminary tests revealed that the developed algorithm reproduces high quality images. The algorithm is moderately fast, but the result comparable to CLEAN algorithm is obtained within 20-50 iteration steps which takes less than 2 seconds on typical portable computer configuration. The location, size and intensity of reconstructed sources are very close to simulated ones. Therefore the algorithm is very well suited for the detailed photometry of the solar HXR sources. Moreover, its simplicity allows to improve photon transmission calculation in case of any grids uncertainties measured after the launch.