Доступ к файлу.fits и построение его столбцов

Question

Я пытаюсь получить доступ к файлу.fits и построение двух столбцов (из многих!). Я использовал pyfits для доступа к файлу, и

plt.plotfile('3XMM_DR5.fits', delimiter=' ', cols=(0, 1), names=('x-axis','y-axis'))

но это не работает Есть ли альтернативы? И есть ли способ открыть файл с помощью Python? Для доступа к таблице данных

Answer 1

Согласно документам из matplotlib за plotfile:

Примечание: файл plotf удобен для быстрого вывода данных из плоских файлов; он не предназначен в качестве альтернативы общему графику с pyplot или matplotlib.

Это не очень понятно. Я думаю, что под "плоскими файлами" это просто означает данные CSV или что-то в этом роде - в моем опыте эта функция используется не очень часто, и, конечно, она ничего не знает о файлах FITS, которые редко используются вне астрономии. Вы упомянули в своем посте, что вы что-то сделали с PyFITS, но это нигде не продемонстрировано в вашем вопросе.

Кстати, PyFITS уже несколько лет не рекомендуется, и его функциональность интегрирована в Astropy.

Вы можете открыть таблицу из файла FITS с помощью astropy.Table.read:

from astropy.table import Table
table = Table.read('3XMM_DR5.fits')

затем получите доступ к столбцам с квадратной скобкой, например:

plt.plot(table['whatever the x axis column is named'], table['y axis column name'])