python - 使用 astropy 执行坐标投影

Question

我在 equatorial coordinates 中有一组点我需要将它们投影到平面上，即:天顶角或 azimuthal projection .

astropy 显然能够执行这种类型的投影，among many others .

问题是我不知道如何将这些模块应用到我的数据中，而且文档也很少。

对于此示例，假设可以通过以下方式生成以度为单位的赤道 ra, dec 坐标:

ra = np.random.uniform(130., 135., 1000)
dec = np.random.uniform(-55., -57., 1000)

Answer 1

这些投影与 wcs(世界坐标系)对象和 FITS format 一起使用.主要步骤是:

import numpy as np
from astropy import wcs

#create data
ra = np.random.uniform(130., 135., 1000)
dec = np.random.uniform(-55., -57., 1000)

#create wcs object
w = wcs.WCS(naxis=2)
# set coordinate and projection types for each axis
w.wcs.ctype = ["RA---ARC", "DEC--ARC"]

#do the projection
pixel_coords = w.wcs_world2pix(ra, dec, 1) 

注意wcs项还有几个额外的配置参数可能很重要，比如crpixj，轴j的图像引用点，和pvi_m，中间世界坐标轴i参数值。

有关一些用法示例，请参阅 astropy.wcs documentation和 projection unit tests .