問題タブ [healpy]

量子物理学と大学の数学について私が覚えていることから、すべてのモードlに対して が存在しm = l-1, l, l+1ます。なぜ helpix (私の場合、具体的にはhelppy ) 球面調和ルーチン (たとえば、とhealpy.sphtfunc.alm2mapの配列を同じ長さで返す) を行うのですか?lm

私はヘルピーとパイソンに非常に慣れていませんが、銀河座標で銀河面と赤道面を使用してスカイマップを描きたいと思います。私が抱えている最初の問題は、直線を引くことができないということです:

私はこのドキュメントに従おうとしていました: https://healpy.readthedocs.org/en/1.5.0/generated/healpy.visufunc.projplot.html

また、銀河面と赤道面の座標をどこで取得すればよいかわかりません。

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

パート 1: 私が抱えている主な問題は、healpy mollview を使用して半天プロットを実行しようとすることです。私が使用している現在の回避策は次のとおりです。

これは少し不器用に感じるので、問題にアプローチするためのより良い方法を探しています。

パート 2: 2 番目の問題は、タイトルです。タイトルなしでマップを作成する方法はありますか?

前もって感謝します

宇宙マイクロ波背景放射のヘルシーな地図があります。各マップには 3 つのフィールドがあります。一般的な数のマップの最初のフィールドのパワー スペクトルを計算するために、コード内に次の関数を記述しました。

「files」入力は、文字列としてのマップ ファイル名の配列です。上記のコードは、ipython の 10 個のマップに対して次の出力を返します。

これには、リスト内に 10 個の配列があります。各配列には、特定のマップの最初のフィールド パワー スペクトル データがあります。私ができるようにしたいのは、コード内のこのリスト内の個々の配列を呼び出して、データの個々の配列ごとに関数を実行することです。

配列内のページ配列を見ましたか? Pythonで配列のリスト内の要素にアクセスしますが、期待したほど役に立ちませんでした。

healpy の alm2fits に相当するものは何だろうと思っています。

T、E、および B の 3 つの alm が healpy にあるとします。これらをエクスポートして、fits2alm を介して Fortran コードで読み取ることができるようにするにはどうすればよいですか?

hp.fitsfunc.write_alm('alms.fits',[alm_T,alm_E,alm_B])

リストが気に入らないようで、

hp.fitsfunc.mwrfits(filename= 'alm.fits',data=[alm_T,alm_E,alm_B],colnames=['T','E','B'],keys=None)

それらを間違った順序で保存しているように見えるため、fits2alm で 'Inconsistent l^2+l+m+1 -> l,m mapping' エラーが発生します

ありがとう

A.