3

拡張子が 1 つの Fits ファイルがあり、データがそれぞれ 100 要素の 2 つの列のテーブルで構成されているとします。

data = pyfits.open('path2myfile')[1].data
head = pyfits.open('path2myfile')[1].header
print data['field1'] # print an array with 100 elements
print data['field2'] # print another array with 100 elements

ここで、テーブルに新しい列を追加したいと思います。たとえば、100 要素の別の配列である data['field3'] としましょう。

正確にはどうすればよいですか？

4

2 に答える 2

1

Iguananaut が示したように、答えはここにあります: http://pyfits.readthedocs.org/en/latest/users_guide/users_table.html#merging-tables しかし、この質問を回答済みとしてマークするだけです:

cols = [] 
cols.append(
    pyfits.Column(name='field3', format='D', array= arrayContainingTheData)
    )
orig_cols = data.columns
new_cols = pyfits.ColDefs(cols)
hdu = pyfits.BinTableHDU.from_columns(orig_cols + new_cols)
hdu.writeto('newtable.fits')
于 2014-08-09T09:50:13.713 に答える