.. include:: references.txt .. _astropy-wcs: *************************************** World Coordinate System (`astropy.wcs`) *************************************** Introduction ============ World Coordinate Systems (WCSs) describe the geometric transformations between one set of coordinates and another. A common application is to map the pixels in an image onto the celestial sphere. Another common application is to map pixels to wavelength in a spectrum. `astropy.wcs` contains utilities for managing World Coordinate System (WCS) transformations defined in several elaborate `FITS WCS standard`_ conventions. These transformations work both forward (from pixel to world) and backward (from world to pixel). For historical reasons and to support legacy software, `astropy.wcs` maintains two separate application interfaces. The ``High-Level API`` should be used by most applications. It abstracts out the underlying object and works transparently with other packages which support the `Common Python Interface for WCS `_, allowing for a more flexible approach to the problem and avoiding the `limitations of the FITS WCS standard `_. The ``Low Level API`` is the original `astropy.wcs` API and originally developed as ``pywcs``. It ties applications to the `astropy.wcs` package and limits the transformations to the three distinct types supported by it: - Core WCS, as defined in the `FITS WCS standard`_, based on Mark Calabretta's `wcslib`_. (Also includes ``TPV`` and ``TPD`` distortion, but not ``SIP``). - Simple Imaging Polynomial (`SIP`_) convention. (See :doc:`note about SIP in headers `.) - Table lookup distortions as defined in the FITS WCS `distortion paper`_. .. _pixel_conventions: Pixel Conventions and Definitions --------------------------------- Both APIs assume that integer pixel values fall at the center of pixels (as assumed in the `FITS WCS standard`_, see Section 2.1.4 of `Greisen et al., 2002, A&A 446, 747 `_). However, there’s a difference in what is considered to be the first pixel. The ``High Level API`` follows the Python and C convention that the first pixel is the 0-th one, i.e. the first pixel spans pixel values -0.5 to + 0.5. The ``Low Level API`` takes an additional ``origin`` argument with values of 0 or 1 indicating whether the input arrays are 0- or 1-based. The Low-level interface assumes Cartesian order (x, y) of the input coordinates, however the Common Interface for World Coordinate System accepts both conventions. The order of the pixel coordinates ((x, y) vs (row, column)) in the Common API depends on the method or property used, and this can normally be determined from the property or method name. Properties and methods containing “pixel” assume (x, y) ordering, while properties and methods containing “array” assume (row, column) ordering. A Simple Example ================ One example of the use of the high-level WCS API is to use the `~astropy.wcs.wcs.WCS.pixel_to_world` to yield the simplest WCS with default values, converting from pixel to world coordinates:: >>> from astropy.io import fits >>> from astropy.wcs import WCS >>> from astropy.utils.data import get_pkg_data_filename >>> fn = get_pkg_data_filename('data/j94f05bgq_flt.fits', package='astropy.wcs.tests') >>> f = fits.open(fn) >>> w = WCS(f[1].header) >>> sky = w.pixel_to_world(30, 40) >>> print(sky) # doctest: +FLOAT_CMP >>> f.close() Similarly, another use of the high-level API is to use the `~astropy.wcs.wcs.WCS.world_to_pixel` to yield another simple WCS, while converting from world to pixel coordinates:: >>> from astropy.io import fits >>> from astropy.wcs import WCS >>> from astropy.utils.data import get_pkg_data_filename >>> fn = get_pkg_data_filename('data/j94f05bgq_flt.fits', package='astropy.wcs.tests') >>> f = fits.open(fn) >>> w = WCS(f[1].header) >>> x, y = w.world_to_pixel(sky) >>> print(x, y) # doctest: +FLOAT_CMP 30.00000214673885 39.999999958235094 >>> f.close() Using `astropy.wcs` =================== .. toctree:: :maxdepth: 2 Shared Python Interface for World Coordinate Systems Legacy Interface Supported Projections Examples creating a WCS programmatically ======================================== .. toctree:: :maxdepth: 2 Example of Imaging WCS Example of Cube WCS Loading From a FITS File WCS Tools ========= .. toctree:: :maxdepth: 1 wcstools.rst Relax Constants =============== .. toctree:: :maxdepth: 1 relax Other Information ================= .. toctree:: :maxdepth: 1 history validation .. note that if this section gets too long, it should be moved to a separate doc page - see the top of performance.inc.rst for the instructions on how to do that .. include:: performance.inc.rst .. _wcs-reference-api: Reference/API ============= .. toctree:: :maxdepth: 1 reference_api See Also ======== - `wcslib`_ Acknowledgments and Licenses ============================ `wcslib`_ is licenced under the `GNU Lesser General Public License `_.