健康:银河和赤道平面图失败

Question

我对健康和巨蟒非常陌生，但我想用银河平面和赤道面在银河坐标下绘制一个天际图。我要面对的第一个问题是，我不能划定一条海峡：

import healpy as hp
import numpy as np
import pylab as pl

hp.mollview(title="Galactic coordinate map")
hp.graticule()

theta = [90., 90.]
phi = [-180., 180.]


hp.projplot(theta, phi, 'r-', coord='G')
pl.show()

我试图跟踪以下文档：https://healpy.readthedocs.org/en/1.5.0/generated/healpy.visufunc.projplot.html

我也不知道在哪里可以得到星系平面和赤道平面的坐标。

Answer 1

在healpy中的bug未被修复之前，就有可能生成足够多的点，从而形成一行：

import healpy as hp
import numpy as np
import pylab as pl

hp.mollview(title="Galactic coordinate map")
hp.graticule()

theta = [0.]*100
phi = np.arange(0, 360, 3.6)


hp.projplot(theta, phi, 'r-', coord='G')
pl.show()