Direct method of image and spectrum recovery in high-energy astronomy

In this paper we propose to use the method of iterative calculation to solve directly the imaging equation for image recovery. The same method can be applied to the recovery of one-dimensional spectra. During the iteration suitable constraints are applied to effect nonlinear control over the process. The use of fitting technique in the treatment of the imaging equation suppresses noise and raises sensitivity. Analyses of both Monte Carlo samples and actual high-energy astronomical data showed that the direct method is more sensitive and having greater resolving power than the traditional method, is capable of revealing simultaneously both diffuse and discrete structures in the target object and can be used on data of lower statistical quality. In principle, the direct method can be used in all types of image and spectrum recovery in other fields.