3

我正在尝试使用 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 掩码吗?

谢谢!

亚历克斯

4

1 回答 1

2

现在没有办法做到这一点healpy

关于这个有一个未解决的问题healpyhttps ://github.com/healpy/healpy/issues/19 。

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

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

于 2014-06-23T16:18:22.633 回答