-1

设置

  • macOS 10.12.6
  • 蟒蛇 2.7
  • 加尔西姆 1.4.4
  • Numpy 1.13.1
  • Astropy 1.3.3
  • 菲西奥

目标

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

问题

不幸的是，我在阅读 fit 文件时遇到了困难。
主要问题是如何以我们可以通过 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 进行卷积或添加检测器功能。

-谢谢

4

1 回答 1

1

如果您阅读 galsim.fits.read 的文档字符串，您会看到它返回一个 galsim.Image，正如您在实验中发现的那样。要创建 GSObject，您需要实例化一个 galsim.InterpolatedImage，这将允许您进行卷积和执行其他操作。您可以阅读 galsim.InterpolatedImage 的文档字符串，和/或查看 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 使用的问题。

于 2017-09-11T15:56:28.753 回答