我一直在疯狂搜索文档,但没有找到答案。 我正在使用Python生成FITS图像,并需要为图像分配WCS坐标。 我知道有很多方法可以通过将点源与已知目录进行匹配来实现此目的,但在这种情况下,我正在生成一个尘埃地图,因此点源匹配不起作用(据我所知)。 因此,该图像是一个形状为(240,240)的2D Numpy数组。 它被编写成下面这样(x和y坐标分配有些奇怪,但它可以正常工作):
H, xedges, yedges = np.histogram2d(glat, glon, bins=[ybins, xbins], weights=Av)
count, x, y = np.histogram2d(glat, glon, bins=[ybins, xbins])
H/=count
hdu = pyfits.PrimaryHDU(H)
hdu.writeto(filename)
>>> print H.shape
(240,240)
这些都可以单独工作得很好。对于分配银河系坐标,似乎你只需要做以下操作:
glon_coords = np.linspace(np.amin(glon), np.amax(glon), 240)
glat_coords = np.linspace(np.amin(glat), np.amax(glat), 240)
但我不理解FITS图像如何存储这些坐标,所以我不知道该如何编写它们。我也尝试在SAO DS9中分配它们,但没有成功。我只需要一个简单直接的方法来将这些坐标分配给图像。
感谢您提供的任何帮助。