使用（本地 wcs）PSF 进行 galsim 星系模拟

Question

我想使用Galsim 来模拟一组PSF 卷积的星系。

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

我喜欢使用的 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 已经被探测器倾斜扭曲的事实？

Answer 1

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

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

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