2

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

thetaphi、は私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同じパターンが得られました。立方体グリッドで作業しているためかどうかはわかりません。グリッドの各セルについて、ボックスの中心に対するデカルト位置から角度位置を計算します。私に何ができる?または、私のコードの何が問題になっていますか? 前もって感謝します。

4

0 に答える 0