Не могу прочитать файл FITS с GalSim

Question

Настроить

• MacOS 10.12.6
• Python 2.7
• Галсим 1.4.4
• Numpy 1.13.1
• Астропия 1.3.3
• Fitsio

Цель

Моя цель - взять смоделированные широкоформатные изображения HST (12288*12288 пикселей, примерно 600 МБ в формате.fits) и добавить в него функции детектора WFIRST.

проблема

К сожалению, я столкнулся с трудностями при чтении файлов подгонки.
Главный вопрос заключается в том, как превратить подгонку изображений таким образом, чтобы мы могли манипулировать ими с помощью galsim.

я пробовал

Я попробовал метод galsim fits galsim.fits.read ("fine_name.fits") и получил:

galsim.Image(bounds=galsim.BoundsI(xmin=1, xmax=12288, ymin=1, ymax=12288), array=
    array([[ 0., 0., 0., ..., 0., 0., 0.],
[ 0., 0., 0., ..., 0., 0., 0.],
[ 0., 0., 0., ..., 0., 0., 0.],
...,
[ 0., 0., 0., ..., 0., 0., 0.],
[ 0., 0., 0., ..., 0., 0., 0.],
[ 0., 0., 0., ..., 0., 0., 0.]], dtype=float32), wcs=galsim.PixelScale(1.0))

я пытался fitsio.read(“file_name.fits”) но, тем не менее, я не могу свернуть PSF с помощью "galsim.Convolve" или добавить функции детектора.

-Спасибо

Answer 1

Если вы прочтете строку документации для galsim.fits.read, вы увидите, что она возвращает galsim.Image, как вы обнаружили в ходе эксперимента. Чтобы создать GSObject, вам нужно создать экземпляр galsim.InterpolatedImage, который позволит вам свертывать и выполнять другие операции. Вы можете прочитать строку документации для galsim.InterpolatedImage и / или посмотреть galsim demo11.py, чтобы увидеть пример того, как это делается на практике (чтение файла изображения в InterpolatedImage и свертывание чего-то еще).

Функции детектора добавляются к изображениям, а не к объектам GSObject. Пример этого см. В demo13.py от GalSim.

Вопросы базового использования GalSim описаны в Кратком справочном руководстве https://github.com/GalSim-developers/GalSim/blob/releases/1.4/doc/GalSim_Quick_Reference.pdf и демонстрационных материалах в репозитории GalSim.