galsim galaxy simulation with（local wcs）PSF

Question

我想用Galsim来模拟一组由PSF卷入的星系。

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

Answer 1

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

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

我们会正确修复吗？ galsim会注意到它只需要扭曲星系并通过PSF进行卷积（并且不会使PSF失真，使PSF对星系进行卷积并扭曲旋转的星系）？