File: wcs_utils.rst

package info (click to toggle)
specutils 2.3.0-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 12,880 kB
  • sloc: python: 13,735; makefile: 111
file content (39 lines) | stat: -rw-r--r-- 1,850 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
=============
WCS Utilities
=============

The `~specutils.utils.wcs_utils` module has functions for converting spectral values
between air and vacuum, as well as a function for calculating the refractive index
of air, which is used in the air to vacuum conversions. There are multiple methods
available for the refractive index calculation, most of which agree with each other
to a high level of precision. However, the Greisen (2006) equation seems to have a small
offset from the others of about 0.07 Angstrom (at the wavelengths shown in the plot below)
and thus in specutils 1.17.0 we have changed from using that method as the default to
using the Morton (2000) equation by default, which is consistent with the IAU standard.

The source of the discrepancy between Greisen (2006) and the other methods seems to be
that Greisen (2006) assumes an air temperature of 0C, versus an apparent assumed air temperature of 15C
for the other methods (thanks to Jon Holtzman for investigating this). Users should keep this
in mind when choosing which method to use for air to vacuum conversions.

.. image:: air_to_vac_offset.png
   :alt: Results of air to vac conversion at optical wavelengths.

The downside of all methods but Greisen (2006) is that they have mathematical singularities
in the far UV, and thus are only valid at wavelengths longer than 200 nm. The specutils
conversion functions will raise an error if these methods are used for wavelengths shorter
than this limit.

.. image:: refractive_index_singularities.png
   :alt: Demonstration of singularities in the equations for refractive index.

For additional context and discussion, see https://github.com/astropy/specutils/issues/1162.

Reference/API
-------------
.. automodapi:: specutils.utils.wcs_utils
    :no-heading:

    :skip: GWCS
    :skip: Tabular1D
    :skip: SpectralGWCS