1

我想Galsim用来模拟一组由 PSFs 卷积的星系。

星系是我使用 galsim 中的类创建的简单双 seric 轮廓Sersic（然后通过剪切引入一些椭圆度）。

我喜欢使用的PSF 是从 Zemax或 CodeV 光线追踪模拟计算的像素化图像。它们是在对应于检测器表面的网格上计算的。该表面相对于主射线倾斜，因此这些 PSF 已经包含探测器倾斜的 WCS失真。

我想计算探测器表面上的 PSF 卷积星系。一种可能的方法是：

  • 加载 psf

psf = galsim.InterpolatedImage(galsim.Image(psf))

  • 然后由星系卷积：

gal = galsim.Convolve(psf, gal)

  • 然后在检测器表面上绘制（在 local_wcs 中倾斜）：

gal_image = gal.drawImage(wcs=local_wcs)

我认为我犯了一个错误，因为 PSF 蜜蜂受到了两次扭曲的影响（Zemax PSF 中的原始倾斜 + drwImage 方法的 local_wcs）。我的担心正确吗？

我是否应该将 local_wcs 扭曲应用于原始未卷积星系（通过应用与 local_wcs 对应的剪切），然后通过 psf 进行卷积并将其绘制在未扭曲的 wcs 上？这是否正确考虑到我的 PSF 已经被探测器倾斜扭曲的事实？

4

1 回答 1

1

我遇到了一个可能的解决方案，那就是在加载 PSF 时指定 local_wcs ：

psf = galsim.InterpolatedImage(galsim.Image(psf), wcs=local_wcs)

我们会正确修复吗？galsim 会注意到它只需要扭曲星系并通过 PSF 卷积（而不是不扭曲 PSF，用 PSF 卷积星系并扭曲卷积星系）？

于 2018-03-16T13:47:41.513 回答