Извлечение третьих осей, кубов данных кубиков

Question

У меня есть файл данных куба fitsfile, состоящий из интенсивности на регулярной (512x512) 2d сетке. Однако есть и третьи оси; каждое значение на этой оси соответствует частоте. На каждой частоте выполняется такое 2d-изображение.

Я пытался найти эти значения частоты в файле fits, но до сих пор мне не повезло. Я знаю, что есть 140 различных частотных каналов, я знаю единицы и т. Д., Глядя на заголовок файла fits. Проблема заключается в нахождении значений этих 140 каналов.

Программирование в PYTHON, кстати.

Answer 1

Обычно это указывается так же, как и для оси «n». Контрольное значение «CRVALn», опорная позиция «CRPIXn», приращение для каждого пикселя/канала «CDELTn» и тип координат, заданных в этих ключевых словах «CTYPEn». Так с, что вы можете просто отработать значения каждого «канала» в кубе, как

(channel_number - CRPIXn) * CRDELTn + CRVALn 

Так как вы сказали, что вы решили ее, возможно, вы могли бы разместить свое решение в пользу других?

---

Редактировать/Дополнение:

Более подробной информации о том, как это выражается в FITS файлов: http://fits.gsfc.nasa.gov/fits_wcs.html

особенно бумаги "Представление спектральных координат в FITS" Грейзен и др. (2005) (link на бумаге по НАСА ADS)