0

この質問には 2 つの部分があります。

パート 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: 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 に答える