使用 healpy 和 projplot 的 Aitoff 预测

Question

我使用 healpy pix2ang 命令有一个 theta 和 phi 的范围，然后转换为 RA, Decl.::

ra = np.rad2deg(phi)
dec = np.rad2deg(0.5 * np.pi - theta)

我只想将这些投影到例如Aitoff 类型的投影，但对于我的生活来说，我无法弄清楚如何通过 :: https://healpy.readthedocs.io/en/latest/generated/healpy.visufunc.projplot.html

projplot(ra, dec, 'bo')

并没有真正做任何事情。

Answer 1

hp.projplot 用于向 现有绘图 添加线条。如果您只是对在不同的投影上绘制线条感兴趣，我建议您查看 matplotlib's projections.

对于 healpy，请在下面找到一个快速示例。

import healpy as hp
import numpy as np

nside = 64
npix = hp.nside2npix(nside)
arr = np.random.randn(npix)

# Draw a circle
r = np.full(100, 20.)
phi = np.linspace(0., 2*np.pi, 100)
x = np.cos(phi)*r
y = np.sin(phi)*r

# Plot the map and the circle
hp.mollview(arr)
hp.projplot(x, y, c='r', lonlat=True)

使用 healpy 和 projplot 的 Aitoff 预测

Aitoff projections using healpy and projplot

python

plot

matplotlib

healpy