Question

设置

MacOS 10.12.6
Python 2.7
Galsim 1.4.4
Numpy 1.13.1
Astropy 1.3.3
Fitsio

目标

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

问题

不幸的是，我在阅读拟合文件时遇到了困难主要问题是如何以galsim操纵它们的方式转换图像。

我已尝试

我尝试过galsim fit方法'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”)但是，我仍然无法通过＆＃39; galsim.Convolve＆＃39;来解决PSF问题。或添加探测器功能。

-Thanks

Answer 1

如果您阅读galsim.fits.read的docstring，您会看到它返回一个galsim.Image，就像您在实验中找到的一样。要创建GSObject，您需要实例化一个galsim.InterpolatedImage，它允许您进行卷积和执行其他操作。您可以阅读galsim.InterpolatedImage的docstring，和/或查看galsim demo11.py以查看在实践中完成此操作的示例（将图像文件读入InterpolatedImage并与其他内容进行卷积）。

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

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

无法使用GalSim读取FITS文件

1 个答案: