注釈
クリック here 完全なサンプルコードをダウンロードします
マルチ拡張マッチ(MEF)ファイルのテーブルに格納されているデータへのアクセス¶
FITファイルは、一般に、大量の多次元データおよびテーブルを含むことができる。この例は、ChandraからのHETG-S機器の情報を含むFITSファイルを開く。
この例で使用する astropy.utils.data
マルチ拡張マッチ(MEF)ファイルをダウンロードするには、 astropy.io.fits
ヘッダを調べるには、以下の操作を実行してください astropy.table.Table
これらのデータを探索しています
By: Lia Corrales, Adrian Price-Whelan, and Kelle Cruz
ライセンス:BSD
使用 astropy.utils.data
本例で用いたFITSファイルをダウンロードするためのサブパケット。導入も Table
従 astropy.table
同前の和. astropy.io.fits
from astropy.utils.data import get_pkg_data_filename
from astropy.table import Table
from astropy.io import fits
FITSファイルをダウンロードする
event_filename = get_pkg_data_filename('tutorials/FITS-tables/chandra_events.fits')
FITSファイルの内容に関する情報を表示する.
出て:
Filename: /home/bk/.astropy/cache/download/url/333246bccb141ea3b4e86c49e45bf8d6/contents
No. Name Ver Type Cards Dimensions Format
0 PRIMARY 1 PrimaryHDU 30 ()
1 EVENTS 1 BinTableHDU 890 483964R x 19C [1D, 1I, 1I, 1J, 1I, 1I, 1I, 1I, 1E, 1E, 1E, 1E, 1J, 1J, 1E, 1J, 1I, 1I, 32X]
2 GTI 3 BinTableHDU 28 1R x 2C [1D, 1D]
3 GTI 2 BinTableHDU 28 1R x 2C [1D, 1D]
4 GTI 1 BinTableHDU 28 1R x 2C [1D, 1D]
5 GTI 0 BinTableHDU 28 1R x 2C [1D, 1D]
6 GTI 6 BinTableHDU 28 1R x 2C [1D, 1D]
拡張1、イベントは、チャンドラにヒットしたHETG-S検出器の各X線光子に関する情報を含むテーブルである。
使用 Table
表を読むには、以下の操作を実行してください
events = Table.read(event_filename, hdu=1)
イベントテーブルの列名を印刷する.
print(events.columns)
出て:
<TableColumns names=('time','ccd_id','node_id','expno','chipx','chipy','tdetx','tdety','detx','dety','x','y','pha','pha_ro','energy','pi','fltgrade','grade','status')>
列が単位情報を含む場合,それは関連するものとなる. astropy.units
物体です。
print(events['energy'].unit)
出て:
eV
“エネルギー”列に格納されているデータを印刷します。
print(events['energy'])
出て:
energy
eV
---------
13874.715
2621.1938
12119.018
3253.0364
14214.382
1952.7239
3267.5334
3817.0366
2252.7295
6154.1094
...
4819.8286
12536.866
2599.5652
15535.768
6653.0815
14362.482
14653.954
6652.827
9672.882
1875.9359
Length = 483964 rows
スクリプトの総実行時間: (0分0.310秒)