python - 将刻度标签添加到 healpy mollview

Question

我正在尝试使用 healpy 绘制 HEALPix 掩码，例如:

    import healpy as hp
    import matplotlib
    from pylab import *

    # Read in mask...
    # ...

    # Plot mask
    fig = figure(figsize=(12,8))
    ax = hp.projaxes.HpxMollweideAxes(fig,[0.1,0.1,0.8,0.8],rot=(180.0,0.0,0.0),coord=["C"])
    fig.add_axes(ax)
    ax.projmap(mask,nest=False)
    hp.visufunc.graticule()
    show()

但是，我想在轴上显示 RA/Dec (phi/theta) 刻度标签。有人知道怎么做这个吗？

哦，还有，有人有 python 代码可以在锤子投影图上绘制 HEALPix 掩码吗？

谢谢!

亚历克斯

Answer 1

现在无法使用healpy

healpy 上有一个 Unresolved 问题:https://github.com/healpy/healpy/issues/19 .

另一种方法是直接使用 matplotlib 绘制 map :

https://github.com/zonca/paperplots/blob/master/python/scripts/PlanckFig_map.py