用(局部wcs) PSF模拟galsim星系

Question

我想用Galsim来模拟一组星系，这些星系由PSFs聚合而成。

星系是我用galsim的Sersic类创建的简单的双螺旋轮廓(然后剪切来介绍一些椭圆度)。

我喜欢使用的PSF是从泽马克斯、或CodeV射线追踪模拟计算出来的像素图像。它们是在与探测器表面相对应的网格上计算的。这个表面是倾斜的wrt，所以，这些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两次受扭转影响(原始倾斜在泽马克斯PSF + local_wcs的drwImage方法)。我的担心正确吗？

我是否应该将local_wcs变形应用于原始的非卷积星系(通过应用对应于local_wcs的剪切)，然后由psf收缩，并将其绘制到一个非扭曲的wcs上？这是否正确地考虑到我的PSF已经被探测器倾斜扭曲了？

Answer 1

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

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

这会不会是个正确的解决办法？盖尔西姆是否会注意到，它只需要扭曲星系，并由PSF (而不是不扭曲的PSF，由PSF收缩星系，扭曲旋转的星系)？