Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://www.adass.org/adass/proceedings/adass97/jennesst.html
Дата изменения: Fri May 15 22:17:27 1998 Дата индексирования: Tue Oct 2 02:39:53 2012 Кодировка: Поисковые слова: п п п п п п п п |
Next: ISDC Data Access Layer
Up: Data Analysis Applications
Previous: A Software Package for Automatic Reduction of ISOPHOT Calibration Data
Table of Contents -- Index -- PS reprint -- PDF reprint
T. Jenness and J. F. Lightfoot1
Joint Astronomy Centre, 660 N. A`ohoku Place, University Park,
Hilo, HI 96720, USA
1Royal Observatory Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, United Kingdom
The Submillimetre Common-User Bolometer Array (SCUBA) is the recently commissioned instrument on the James Clerk Maxwell Telescope (JCMT) on Mauna Kea. SCUBA is a bolometer array which means that for the first time JCMT has imaging capabilities in the submillimetre.
In order to maximize the number of bolometers in the array, they were packed on a hexagonal (rather than rectangular) grid. This means that the software has to reconstruct the data into an image. Moreover, because the bolometers are larger than half the beam width, the image cannot be fully sampled in a single exposure. In order to fully sample the image, the secondary mirror is ``jiggled'' to several adjacent positions and the software has to keep track of what sky position is associated with each data measurement. Finally, because of how SCUBA is mounted (on the Nasmyth platform), the sky rotates as the telescope tracks; this is also something the software must correct for.
There are three major SCUBA observing modes, and these need some individual treatment from SURF, the Scuba User Reduction Facility, which is run off-line by the user:
The real-time SCUBA observing system first demodulates data in the transputers - it takes data at 128Hz but only supplies data every second. It also provides a `quick-look' display of the data based on some rudimentary data reduction done in the transputers.
The demodulated data is then ``raw data'' as far as SURF, the off-line system, is concerned. SURF aims to provide publication quality data reduction and currently requires user input and feedback.
The `quick-look' display provided by the transputers can also be eavesdropped by a remote observer using a WORF system (Economou et al 1996; Jenness et al 1997).
SURF (like the real-time system) is written in the Starlink software environment and therefore runs on Sun Solaris, Digital Unix and Linux. Computationally intensive routines are written in FORTRAN as ADAM tasks that form part of a monolith. Applications such as data reduction pipelines (calling multiple tasks from the monolith), string handling routines (e.g., generating observing logs) and small data reduction tasks were more suited to modern scripting languages than to FORTRAN and it was decided to use the perl programming language (Wall et al. 1996) for these.
Since the software is written in the Starlink environment, Starlink's N-Dimensional Data Format (NDF) is used for all file I/O. This format provides for error arrays, bad pixel masks, storage of history information and a hierarchical extension mechanism and SURF makes full use of these features. In order to get the most out of using perl it was necessary to write a perl interface to the FORTRAN NDF and HDS libraries. This module (distributed separately from SURF) provides complete access to all of the routines in these libraries allowing for full control of the NDF and HDS structures from perl.
To improve maintainability of the package all source code is kept under revision control (RCS) and subroutine libraries are shared with the on-line observing system.
SURF provides the tasks necessary to convert raw demodulated SCUBA data into a regularly-sampled image:
The different response of each bolometer is removed. Note that the flatfield remains constant with time and does not need to be re-measured every night.
This task simply calculates the extinction due to the atmosphere by using the zenith sky opacity (tau) and calculating the elevation of each bolometer during each 1 second sample.
Occasionally spikes are present in the data (instrumental glitches, cosmic rays etc) and these can be removed manually or by using software.
The submillimetre sky is very sensitive to sky noise caused by fluctuations in the emissivity of the atmosphere passing over the telescope. These variations occur in atmospheric cells that are larger than the array and the resultant noise is seen to be correlated across the array. At present this sky noise is removed by analysing bolometers that are known to be looking at sky and removing this offset from the data.
As described earlier, SCUBA data is under-sampled (for jiggle data), taken on an hexagonal grid and subject to sky rotation. SURF calculates the position of each bolometer during each 1 second sample (in a number of coordinate systems: RA/Dec, Galactic, Nasmyth, Azimuth-Elevation or moving centre (planet)) and regrids these data onto a rectangularly sampled image.
Once the data have been regridded the image can be converted to FITS (a full astrometry header is provided), if necessary, and analysed in the normal way.
The latest version of SURF is freely available from the SURF home page . Additionally the software is distributed by Starlink and available on their CD-ROM.
This software is being developed actively at the Joint Astronomy Centre and we hope to implement the following features in the next few months:
TJ would like to thank Frossie Economou for her helpful discussions during the development of this software.
Economou, F., Bridger, A., Daly, P. N., & Wright, G. S., 1996, in Astronomical Data Analysis Software and Systems V, ASP Conf. Ser., Vol. 101, eds. G. H. Jacoby and J. Barnes (San Francisco, ASP), 384
Economou, F., Bridger, A., Wright, G. S., Rees, N. P., & Jenness, T., 1998, this volume
Emerson, D. T., Klein, U., & Haslam, C. G. T., 1979, ApJ, 76, 92
Jenness, T., Economou, F., & Tilanus, R. P. J., 1997, in Astronomical Data Analysis Software and Systems VI, ASP Conf. Ser., Vol. 125, eds. Gareth Hunt and H. E. Payne (San Francisco, ASP), 401
Richer, J. S., 1992, MNRAS, 254, 165
Wall, L., Christiansen, T., & Schwartz, R. L., 1996, Programming Perl, 2nd edn. (Sebastopol, CA: O'Reilly ) (http://www.perl.org/ )
Next: ISDC Data Access Layer
Up: Data Analysis Applications
Previous: A Software Package for Automatic Reduction of ISOPHOT Calibration Data
Table of Contents -- Index -- PS reprint -- PDF reprint