无法使用 GalSim 读取 FITS 文件

Question

设置

MacOS 10.12.6
Python 2.7
Galsim 1.4.4
麻木 1.13.1
Astropy 1.3.3
Fitsio

目标

我的目标是拍摄模拟的宽视场 HST 图像（12288*12288 像素，.fits 格式大约 600MB）并向其中添加 WFIRST 探测器功能。

问题

不幸的是，我在阅读 fits 文件时遇到了困难。
主要问题是如何以我们可以通过 galsim 操作它们的方式来调整图像。

我试过了

我尝试了 galsim 拟合方法‘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”)，但仍然无法通过 'galsim.Convolve' 对 PSF 进行卷积或添加检测器功能。

-谢谢

Answer 1

如果您阅读 galsim.fits.read 的文档字符串，您会发现它 returns 是 galsim.Image，正如您在实验中发现的那样。要创建一个 GSObject，您需要实例化一个 galsim.InterpolatedImage，这将允许您进行卷积和执行其他操作。您可以阅读 galsim.InterpolatedImage 的文档字符串，and/or 查看 galsim demo11.py 以查看在实践中完成的示例（将图像文件读入 InterpolatedImage 并与其他内容进行卷积） .

检测器功能 do 被添加到图像，而不是 GSObjects。有关此示例，请参阅 GalSim 的 demo13.py.

快速参考指南 https://github.com/GalSim-developers/GalSim/blob/releases/1.4/doc/GalSim_Quick_Reference.pdf 和 GalSim 存储库中的演示涵盖了基本 GalSim 使用问题。

无法使用 GalSim 读取 FITS 文件

Can't read FITS file with GalSim

fits

galsim