我使用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')
真的什么也没做。
最佳答案
hp.projplot
用于将线添加到现有图。如果您只想在其他投影上绘制线,建议您参考matplotlib's projections。
要恢复健康,请在下面找到一个简单的示例。
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)
关于python - 使用healpy和projplot的Aitoff投影,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/57560428/