Astropy：スピーディーフィットファイルからの物理座標

Question

フィッティングfiles （ログインする必要はありません。表示してダウンロードを続行してください）。

ない私の知る限り、関連する（私が間違っているなら、私を修正）：

Find physical coordinates of a pixel in a fits file with python

再現：

from astropy.io import fits 
from astropy.wcs import WCS 
ev=fits.open('sw00032188013xpcw3po_cl.evt.gz') 
w=WCS(ev['events'].header) 
w.wcs_pix2world(290,500,0) 

結果：

[array(291.0), array(501.0)] 

w.to_header（）はOKと思われますが、私はw.wcs_pix2world(290,500,0,ra_dec_order=True)を使用して予告しましたが、奇妙なエラー：

ValueError: WCS does not have longitude type of 'RA', therefore (ra, dec) data can not be returned 

ds9はこのファイルを読み取ってfk5座標を表示しても問題ありません。

事前に感謝します。

PS： 私はあなたがXRTのイベントファイルを使用しているastropy.version.version = '0.4.2' で、

Answer 1

を安定したDebianの8パッケージを使用しています。それはイメージファイルではありません。 DS9はあなたのためにいくつかの解釈を行い、WCSを含むそれから自動的にイメージを作成することができます。 astropy.wcsは、それよりも少ないと解釈し、適切なWCSを使用してイメージを作成することは（簡単に）できません。
（あなたが個々の画像としてDS9の画像を保存することを選ぶことができます。）

代わりに、XRTイベントファイルから画像を作成する方法を見つけるためにSwift XRT manualをお読みください。これは、GTIとバッドピクセルも処理します。
通常、これはHEASoftソフトウェアスイートのusing some toolsで行われます。

また、実際のイメージ（上記の方法で作成されたもの）を含め、自動化されたレベル3の製品available for downloadがあると思います。代わりにそれらを使う方が良いかもしれません。