伝統的なインターフェース

Asterpy.wcsインタフェース

♪the Low Level API あるいは…。 Legacy Interface オリジナルですか。 astropy.wcs 原料の薬。3つのタイプの変換をサポートしています

  • コアWCS、例えば FITS WCS standard マーク·カラブレタの wcslib それがそうです。(他にも含む) TPV そして TPD 歪んでいますがそうではありません SIP ）。

  • 単純結像多項式 (SIP) 大会です。(参照されたい) note about SIP in headers ()

  • FITS WCSで定義されているテーブルは歪みを探す. distortion paper それがそうです。

これらの変換の各々は単独で使用することができ、標準パイプで一緒に使用することもできる。すべてのメソッドはスカラーと配列入力をサポートしている.全ての方法には追加の位置パラメータが必要です origin 入力されたデータ。2つの可能な値があります 0 -Numpy配列のようなゼロからの座標の場合、または 1 −FITS規格のような1に基づく座標、またはDS 9からの座標。

基本ワークフローは、使用するWCSオブジェクトを作成することです Header および/または HDUList オブジェクトを呼び出し、以下の方法の1つを呼び出す：

>>> from astropy import wcs
>>> from astropy.io import fits
>>> 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)
>>> wcsobj = wcs.WCS(f[1].header)

代替的に、FITSファイルが使用を推奨しない機能または非標準機能を使用する場合、呼び出しが必要となる場合がある。 fix 方法です。

以下の変換方法の1つを用いる.

  1. 画素と世界座標の間にすべての歪みを使用します：

  • all_pix2world ：画素から世界座標へのすべての3つの変換(コアWCS、SIP、およびルックアップテーブル歪み)を連続的に実行する。どれを使うべきかわからなければ、これを使います。

    >>> lon, lat = wcsobj.all_pix2world(30, 40, 0)
>>> print(lon, lat)  
5.528442425094046 -72.05207808966726
  • all_world2pix ：全3つの操作を実行する

    世界座標から画素座標への変換(コアWCS，SIP，テーブル探索歪み)は，必要に応じて反復方法を用いる.

    >>> x, y = wcsobj.all_world2pix(lon, lat, 0)
>>> print(x, y) # 
30.00000214673885 39.999999958235094

  1. パフォーマンス SIP 変換のみ：

    • sip_pix2foc ：画素からへの変換

      焦点面座標使用 SIP 多項式係数。

      >>> xsip, ysip = wcsobj.sip_pix2foc(30, 40, 0)
>>> print(xsip, ysip)  
-1985.8600487630586 -984.4223711273145
    • sip_foc2pix ：焦点からの変換

      平面から画素への座標使用 SIP 多項式係数。この方法は、ヘッダに逆SIP歪みが指定されている場合にのみ有効であることに留意されたい。

  2. パフォーマンス distortion paper 変換のみ：

    • p4_pix2foc ：画素からへの変換

      FITS WCSで記述されたルックアップテーブル歪み方法を用いた焦点面座標 distortion paper それがそうです。

    • det2im ：検出器からの変換

      座標から画像座標まで。一般に狭列補正に用いられる。

コアwcslb API

コアwcslb APIは、WCSファイルI、II、III、IVで定義されているFITS WCS規格をサポートしています。コアライブラリ内の関数を使用すると、歪みが適用されないことに注意されたい。

  1. ピクセルから世界座標まで：

    • wcs_pix2world ：コアWCSのみを実行する

      画素から世界座標への変換。

      >>> lon, lat = wcsobj.wcs_pix2world(30, 40, 0)
>>> print(lon, lat)  
5.527103615238458 -72.0522441352217

  2. 世界座標から画素座標まで：

    • wcs_world2pix ：コアWCS変換を実行する

      世界座標から画素座標まで.

      >>> x, y = wcsobj.wcs_world2pix(lon, lat, 0)
>>> print(x, y)  
30.000000000223267 40.0000000003696