1

FITSファイルヘッダーから情報を取得し、その情報をASCIIテーブルに転送する方法を学びたいです。たとえば、これは私が情報を取得する方法です。

import pyfits
a = pyfits.open('data.fits')
header = a[0].header # Which should return something like this (It is  BinHDUlist)
SIMPLE  =                T / conforms to FITS standards                    
                               / institution responsible for creating this file 
TELESCOP= 'Kepler  '           / telescope                                      
INSTRUME= 'Kepler Photometer'  / detector type                                  
OBJECT  = 'KIC 8631743'        / string version of KEPLERID                     
RA_OBJ  =           294.466516 / [deg] right ascension                          
DEC_OBJ =            44.751131 / [deg] declination

RA_OBJとDEC_OBJを含むASCIIテーブルを作成するにはどうすればよいですか？

編集：ヘッダーから2つの列（RAとDEC）を含む.datファイルを作成したいと思います。これが私が試していることの例です：

import asciitable
import asciidata
import pyfits
import numpy as np

# Here I have taken all the fits files in my current directory and did the following:
# ls > z.txt so that all the fits files are in one place.

a = asciidata.open('z.txt')
i = 0 #There are 371 fits files in z.txt
while i<=370:
    b = pyfits.open(a[0][i])
    h = b[0].header
    RA = np.array([h['RA_OBJ']])
    DEC = np.array(h['DEC_OBJ']])
    asciitable.write({'RA': RA, 'DEC': DEC}, 'coordinates.dat', names=['RA', 'DEC'])
    i = i+1

次のようなものを含む.datファイルを作成したいと思います。

RA    DEC
###   ###
...   ...
...   ...
...   ...

代わりに、私のコードは前のファイルのキーを上書きするだけです。何か案は？

4

1 に答える 1

1

pyfitsのドキュメントをもっと注意深く読むことでメリットが得られると思います。ヘッダー属性はpyfits.header.Header、辞書のようなオブジェクトであるオブジェクトです。したがって、次のようなことができます。

import pyfits

keys = ['SIMPLE', 'TELESCOP', 'INSTRUME', 'OBJECTS', 'RA_OBJ', 'DEV_OBJ']

hdulist = pyfits.open("data.fits")
header = hdulist[0].header
for k in keys:
    print k, "=", header[k]

より凝った出力を追加したり、結果の文字列を変数に入れたり、欠落しているキーをチェックしたりすることができます。

編集：

asciitableこれをとと組み合わせる方法は次のnumpyとおりです。

import asciitable
import numpy as np

keys = ['RA', 'DEC']
data = {}

# Initialize "data" with empty lists for each key
for k in keys:
    data[k] = []

# Collect all data in the "data" dictionary
for i in range(0, 50):
    data['RA'].append(np.array(i))
    data['DEC'].append(np.array(i+1))

asciitable.write(data, "coords.dat", names=keys)
于 2012-07-17T16:57:32.120 に答える