Proceedings of the Workshop
"The Magellanic Clouds and Other Dwarf Galaxies"
of the Bonn/Bochum-Graduiertenkolleg

A detailed view of the Magellanic Clouds in the FIR

M. Braun1, R. Assendorp1, Tj.R. Bontekoe2,3, D.J.M. Kester4, and G. Richter1

1Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 14482 Potsdam-Babelsberg, Germany
2Sterrenkundig Instituut Anton Pannekoek, Kruislaan 403, 1098 SJ Amsterdam, NL
3Bontekoe Data Consultancy, Herengracht 47, 2312 LC Leiden, NL
4Space Research Organization Netherlands, P.O.Box 800, 9700 AV Groningen, NL

Received 15th March 1998
Abstract. High-resolution 12, 25, 60, and 100 µm IRAS images of the Magellanic Clouds have been produced using a maximum entropy method. These images which are available for public use have the highest dynamic range and the highest resolution currently possible for IRAS data. In this contribution we describe the basics of the data processing. Additionally, the increased image quality is demonstrated by some few examples.

1. Introduction

The Magellanic Clouds have been observed by the Infrared Astronomical Satellite (IRAS, IRAS Explanatory Supplement (IES), 1988) during the all-sky survey and several Pointed Observations (POs). In order to study global properties and extended regions it is necessary to make images from the data. However, creating images from IRAS data is a complicated matter due to the scan nature of the data. The best results for the Magellanic Clouds to date have been achieved by Schwering and Israel (1991) who processed Pointed Observation data into images. The dynamic range and spatial resolution of these images is high, higher than e.g. that of the Infrared Sky Survey Atlas (ISSA, 1994).

In this contribution we present high resolution images of the Magellanic Clouds produced by a maximum entropy algorithm. The spatial resolution of these images is the best currently possible. At present we are in the process of scientific validation of the images in order to determine the reliability of positions and fluxes of the numerous point sources and extended structures that can be recognized from the images. To achieve this, the images are compared with data taken at other wavelengths, ranging from radio to X-ray.

The images are available for public use. Due to their size we decided not to distribute them via Internet. Requests for copies on tape or other storage mediums should be directed to the authors (e.g. MBraun@aip.de).

Readers interested in the second part of the talk dealing with `the possible interaction of high-velocity clouds with the disk of the LMC' shall be refered to Braun (1996) where this topic is discussed in detail.

2. Data processing

IRAS has surveyed the entire sky at four wavelengths (12, 25, 60, and 100 µm) with a regular pattern of long scans running from pole to pole in the ecliptic coordinate system. IRAS images, such as in the ISSA, are constructed by combining the scans. They have a linear spatial resolution of about 5', corresponding to the largest dimension of the IRAS detectors in the focal plane (IES). The imaging quality of the optics, however, allows a much better resolution: about 25" - 100" depending on the wavelength band. The sampling rate during the scans was such that the sky was almost critically sampled for each wavelength, but only in the in-scan direction. High-resolution information in the cross-scan direction is available because the IRAS observing strategy ensured sufficient redundant coverage.

It is possible to reconstruct reliable images at a resolution of nearly the diffraction limit of the IRAS telescope (Bontekoe et al. 1994). This is achieved by application of the maximum entropy package MemSys5 (Gull and Skilling 1991) which largely removes the blur of the IRAS detector response functions. The resulting spatial resolution and the dynamic range in the high-resolution image depends on the number of data samples locally in a map, and is usually not uniform over the map.

An essential ingredient for the production of high-resolution images is the consistency from scan to scan, which have different background levels. In order to bring the scan on a common baseline a first order fit to the background is subtracted. This, and other pre-processing steps are performed using the IRAS data processing routines described by Assendorp et al. (1995), followed by a series of data self-calibration steps during the maximum entropy image reconstruction process.

The scan coverage of the Magellanic Clouds is much larger than average because they are near one of the poles of the orbit of the IRAS spacecraft. The number of scans that have overlap e.g. with the LMC is roughly 600; the number of data points at 60 µm is about 6.3 million. Though the large amount of IRAS data ensures a high dynamic range and spatial resolution, the disadvantage is the sheer computation time which surpasses the capabilities of most high-end workstations. The images presented here are processed on the Convex SPP 1200 computer of the Astrophysikalisches Institut Potsdam. This is a machine with 8 cpu's with 120 MHz clock cycle and 2 GB workspace memory. Even on this machine the processing of the four images of the LMC took more than two year (single) CPU time.

3. Examples

Fig. 1 shows in comparison two 60 µm maps of the Large Magellanic Cloud. The first one was constructed by combining the IRAS survey scans, the second one was computed using the MemSys5 algorithm. The angular resolution is increased by a factor of about 3, while the large scale structures are also conserved (diffuse emission).

Examples of images produced by the different techniques can be compared in more detail in Fig. 2 where the 60 µm emission of the south-western region of the Small Magellanic Cloud is presented. In addition to images as similarly shown in Fig. 1, the upper right panel gives the outcome of the application of MemSys5 to a combined set of IRAS survey and Pointed Observation data. A further increase in resolution is obvious. However, some artifacts are easily detactable close to bright point sources. They can be addressed to problems with the definition of the detector response functions at the lower limit. The lower right panel was constructed by scan combination from both the survey and the PO data. But in this case, only the data provided by the small detectors - 2 of the 16 60 µm detectors had smaller sizes of 1.3'×1.5' instead of 4.8'×1.5' - were used. The construction of this image was only possible due to the extremely high coverage of this area by IRAS POs. Though the signal-to-noise ratio is low in this map, it is a convincing proof of the reliability of the high resolution images provided by the MemSys5 algorithm.

[Click here to see Fig. 1 and 2!]

References


Links (back/forward) to:

First version: 09thAugust,1998
Last update: 25thSeptember,1998

Jochen M. Braun   &   Tom Richtler
 (E-Mail: jbraun|richtler@astro.uni-bonn.de)