healpy で 3 つの numpy 配列をプロットしようとしています。そのうちの 2 つは角度位置のシータとファイに対応し、もう 1 つは温度です。この以前の回答からの情報を使用しました: Plotting a numpy array in healpy to do it、次のように実行しました:
NSIDE = 512
m_sm = numpy.arange(healpy.nside2npix(NSIDE) )
m_sm = m_sm*0.
indx = healpy.pixelfunc.ang2pix(NSIDE, theta, phi, lonlat=False)
m_sm[indx] = Temp
theta
、phi
、は私Temp
の配列です (それぞれに 100 万個近くの要素がありますが、この行との非互換性はないようですm_sm[indx] = Temp
)。次に、次のように投影を実行しmollview
ました。
gmap=healpy.smoothing(m_sm, sigma=numpy.radians(0.3),verbose=False)
cmap = cm.get_cmap('RdBu_r')
cmap.set_under('w')
healpy.mollview(gmap, title=ur"$T(\mathbf{x})$ for Reconstruction 1",
unit=ur"$T(\mathbf{x})$", nest=False, min=gmap.min(),
max=gmap.max(), remove_dip=False, cmap = cmap,
coord=['C','G'])
しかし、すべてのマップを通過する白い線の奇妙なパターンで、次のマップを取得しました。
でもNSIDE = 256
同じパターンが得られました。立方体グリッドで作業しているためかどうかはわかりません。グリッドの各セルについて、ボックスの中心に対するデカルト位置から角度位置を計算します。私に何ができる?または、私のコードの何が問題になっていますか? 前もって感謝します。