python - 将刻度标签添加到健康的mollview

我正在尝试使用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蒙版吗？

谢谢！

亚历克斯

最佳答案

healpy现在无法执行此操作

healpy对此有一个未解决的问题：https://github.com/healpy/healpy/issues/19。

另一种方法是直接使用matplotlib绘制地图：

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