No Arabic abstract
We present a new algorithm, called Multiresolution Regularized Expectation Maximization (MREM), for the reconstruction of gamma-ray intensity maps from COMPTEL data. The algorithm is based on the iterative Richardson-Lucy scheme to which we added a wavelet thresholding step in order to eliminate image-noise in the reconstruction. The wavelet thresholding explicitly accounts for spatial correlations in the data, and adapts the angular resolution locally, depending on the significance of the signal in the data. We compare the performance of MREM to that of the maximum entropy and the Richardson-Lucy algorithms by means of Monte-Carlo simulations of COMPTEL 1.809 MeV gamma-ray line observations. The simulations demonstrate that the maximum entropy and Richardson-Lucy algorithms provide virtually identical reconstructions which are heavily disturbed by image noise. MREM largely suppresses this noise in the reconstructions, showing only the significant structures that are present in the data. Application of MREM to COMPTEL 1.8 MeV gamma-ray line data results in a 1.809 MeV sky map that is much smoother than the maximum entropy or Richardson-Lucy reconstructions presented previously. The essential features of this map are (1) an asymmetric galactic ridge emission reaching from l=45 deg to l=240 deg, (2) a bright localised emission feature in the Cygnus region around (l,b)=(80 deg,0 deg), (3) two emission spots at l=317 deg and l=332 deg situated in the galactic plane, and (4) an extended emission region around (l,b)=(160 deg, 0 deg). Comparison of the MREM map to the simulated reconstructions demonstrates that the 1.809 MeV emission is confined to the galactic plane.
The instrumental line background of the Compton telescope COMPTEL onboard the Compton Gamma-Ray Observatory is due to the activation and/or decay of many isotopes. The major components of this background can be attributed to eight individual isotopes, namely 2D, 22Na, 24Na, 28Al, 40K, 52Mn, 57Ni, and 208Tl. The identification of instrumental lines with specific isotopes is based on the line energies as well as on the variation of the event rate with time, cosmic-ray intensity, and deposited radiation dose during passages through the South-Atlantic Anomaly. The characteristic variation of the event rate due to a specific isotope depends on its life-time, orbital parameters such as the altitude of the satellite above Earth, and the solar cycle. A detailed understanding of the background contributions from instrumental lines is crucial at MeV energies for measuring the cosmic diffuse gamma-ray background and for observing gamma-ray line emission in the interstellar medium or from supernovae and their remnants. Procedures to determine the event rate from each background isotope are described, and their average activity in spacecraft materials over the first seven years of the mission is estimated.
COMPTEL was the Compton telescope on NASAs Compton Gamma Ray Observatory CGRO launched in April 1991 and which was re-entered in June 2000. COMPTEL covered the energy range 0.75 to 30 MeV, and performed a full-sky survey which is still unique in this range, with no followup mission yet approved. This remains a major uncharted region, and the heritage data from COMPTEL are still our main source of information. Data analysis has continued at MPE however, since the data were never fully analysed during the mission or in the period following, and improvements in analysis techniques and computer power make this possible.
We study inversion of the spherical Radon transform with centers on a sphere (the data acquisition set). Such
We present the time-averaged characteristics of the Crab pulsar in the 0.75-30 MeV energy window using data from the imaging Compton Telescope COMPTEL aboard the Compton Gamma-Ray Observatory (CGRO) collected over its 9 year mission. Exploiting the exceptionally long COMPTEL exposure on the Crab allowed us to derive significantly improved COMPTEL spectra for the Crab nebula and pulsar emissions, and for the first time to accurately determine at low-energy gamma-rays the pulse profile as a function of energy. These timing data, showing the well-known main pulse and second pulse at a phase separation of 0.4 with strong bridge emission, are studied together with data obtained at soft/hard X-ray energies from the ROSAT HRI, BeppoSAX LECS, MECS and PDS, at soft gamma-rays from CGRO BATSE and at high-energy gamma-rays from CGRO EGRET in order to obtain a coherent high-energy picture of the Crab pulsar from 0.1 keV up to 10 GeV.
Investigation of image reconstruction from data collected over a limited angular range in X-ray CT remains a topic of active research because it may yield insight into the development of imaging workflow of practical significance. This reconstruction problem is well-known to be challenging, however, because it is highly ill-conditioned. In the work, we investigate optimization-based image reconstruction from data acquired over a limited-angular range that is considerably smaller than the angular range in short-scan CT. We first formulate the reconstruction problem as a convex optimization program with directional total-variation (TV) constraints applied to the image, and then develop an iterative algorithm, referred to as the directional-TV (DTV) algorithm for image reconstruction through solving the optimization program. We use the DTV algorithm to reconstruct images from data collected over a variety of limited-angular ranges for breast and bar phantoms of clinical- and industrial-application relevance. The study demonstrates that the DTV algorithm accurately recovers the phantoms from data generated over a significantly reduced angular range, and that it considerably diminishes artifacts observed otherwise in reconstructions of existing algorithms. We have also obtained empirical conditions on minimal angular ranges sufficient for numerically accurate image reconstruction with the DTV algorithm.