Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.adass.org/adass/proceedings/adass03/reprints/P6-3.pdf Дата изменения: Sat Aug 28 02:47:27 2004 Дата индексирования: Tue Oct 2 10:15:11 2012 Кодировка: Поисковые слова: п п п п п п п п р п

Astronomical Data Analysis Software and Systems XIII ASP Conference Series, Vol. 314, 2004 F. Ochsenbein, M. Al len, and D. Egret, eds.

Representations of distortions in FITS world co ordinate systems
Mark R. Calabretta Australia Telescope National Facility, PO Box 76, Epping, NSW 1710, Australia Francisco G. Valdes National Optical Astronomy Observatories, PO Box 26732, Tucson, AZ 85719, USA Eric W. Greisen National Radio Astronomy Observatory, PO Box O, Socorro, NM 87801-0387, USA Steven L. Allen UCO/Lick Observatory, University of California, Santa Cruz, CA 95064, USA Abstract. Of particular interest to the fields of astromet troscopy, the fourth in a series of papers defining conventions world coordinate information in FITS headers will consider of representing small systematic errors, or distortions. Here preview of work in progress. ry and specfor encoding the problem we present a

1.

Introduction

Standard methods for specifying world coordinate systems (WCS) in FITS images (Hanisch et al. 2001) have been developed by Greisen & Calabretta (2002) and applied by Calabretta & Greisen (2002) to the problem of celestial coordinates (Papers I & II). The extension to spectral coordinate systems (Greisen et al. 2003, Paper III) is also at an advanced stage. However, these methods implicitly assume ideal astronomical instrumentation and are not easily adapted to describe the complex distortions found in some imaging devices. Examples abound, from classical "plate solutions", to spectrometer wavelength calibration. It may also happen that the distortion is "inherent" to the ob ject of study, for example the oblateness of the Sun and Earth. Irregularity of form is carried to an extreme by various minor bodies of the Solar System. Work in progress, summarised here, extends the current FITS WCS formalism by providing methods to describe the distortions inherent in the image coordinate systems of real astronomical data. It is envisaged that a range of distortion functions will be provided, including N-dimensional polynomial, cubic 551 c Copyright 2004 Astronomical Society of the Pacific. All rights reserved.


552

Calabretta, Valdes, Greisen & Allen
PIXEL COORDINATES optional distortion correction CORRECTED PIXEL COORDINATES linear transformation: translation, rotation, skewness, scale INTERMEDIATE PIXEL COORDINATES optional distortion correction CORRECTED INTERMEDIATE PIXEL COORDINATES scale to physical coordinates INTERMEDIATE WORLD COORDINATES coordinate computation per agreement WORLD COORDINATES CTYPEia CRVALia PVi_ma CDELTia CTYPEia DVi_ma CRPIXjs PCi_ja or CDi_ja CTYPEia DVi_ma

Figure 1. Conversion of pixel coordinates to world coordinates showing optional distortion corrections enclosed in the dashed boxes. spline, B-spline, and table lookup methods; that these will be applicable over multiple image dimensions; and that they may be applied either before or after the standard linear transformation stage of the coordinate calculation. An early draft of Paper IV is available1 for comment. We invite input from the general FITS user community regarding the adequacy of the proposed methods for existing or future applications. 2. Methodology

The new steps to be introduced into the algorithm chain are enclosed in dashed boxes in Figure 1. Key features are: * For each axis a distortion function may be applied to either the pixel coordinates or intermediate world coordinates (but not both). * The distortion function is defined in the pixel-to-world direction. Polynomial, cubic spline, B-spline, and table lookup functions will be provided.
1

http://www.atnf.csiro.au/people/mcalabre/


Distortions in FITS world coordinate systems

553

* The CTYPEia header cards will indicate that a distortion function is to be applied, its type, and whether before or after the linear transformation. * New DVi ma header cards perform several functions: - Define which coordinate axes form the independent variables of the distortion function (axis coupling). - Provide an offset and scale for renormalization of the independent variables of the distortion function. - Encode the parameters required for the distortion function. * DVERRja will record the maximum error of the distortion correction on axis j , and DVERRa will record the maximum error of the combined distortion functions for all axes. * Methods will be provided to define different distortion functions for different regions of the image, e.g. as may be required for arrays of CCD detectors. Of course, the exact details are sub ject to change in response to feedback from the FITS user community. Acknowledgments. The Australia Telescope is funded by the Commonwealth of Australia for operation as a National Facility managed by CSIRO. The National Optical Astronomy Observatory is a facility of the (U.S.) National Science Foundation operated under cooperative agreement by Associated Universities for Research in Astronomy, Inc. The National Radio Astronomy Observatory is a facility of the (U.S.) National Science Foundation operated under cooperative agreement by Associated Universities, Inc. UCO/Lick Observatory is operated by the University of California. References Calabretta, M. R. & Greisen, E. W. 2002, A&A, 395, 1075 (Paper II) Greisen, E. W. & Calabretta, M. R. 2002, A&A, 395, 1059 (Paper I) Greisen, E. W., Valdes, F. G., Calabretta, M. R. & Allen, S. L. 2003, A&A, in preparation, Representations of spectral coordinates in FITS (Paper III) Hanisch, R. J., Farris, A., Greisen, E. W., Pence, W. D., Schlesinger, B. M., Teuben, P. J., Thompson, R. W., & Warnock III, A. 2001, A&A, 376, 359