0

这个问题有两个部分:

第 1 部分:我遇到的主要问题是尝试使用 healpy mollview 绘制半边天图。我正在使用的当前解决方法是:

import numpy as np
import healpy as hp
import matplotlib.pyplot as plt

w, h = 6, 4
fig = plt.figure(1, figsize=(w,h))

map = np.arange(hp.nside2npix(64))
hp.mollview(map, fig=1)
hp.graticule(verbose=False)

for ax in fig.get_axes():
    if type(ax) is hp.projaxes.HpxMollweideAxes:
        ax.set_ylim(-1, 0.005)
fig.set_size_inches(w, h/2.8, forward=True)

plt.show()

这感觉有点笨拙,所以我正在寻找一种更好的方法来解决这个问题。

第 2 部分:第二个问题是标题。有没有办法制作完全没有标题的地图?

提前致谢

4

1 回答 1

0

这是一段代码,它绘制了地图的上半部分,没有标题,也没有彩色地图。

#!/usr/bin/env python
import healpy as hp
import pylab as pl
import numpy as np

NSIDE=128
m=np.arange(hp.nside2npix(NSIDE))

hp.mollview(m, title="", cbar=False)
pl.xlim([0,1])

pl.show()

希望能帮助到你。

于 2014-12-17T14:04:20.253 回答