問題タブ [healpy]

私は健全な全天地図を持っており、指定されたピクセルを中心としたいくつかの空のパッチを選択する必要があります。hp.query_discは、指定されたディスク内のピクセルを見つけます。これは素晴らしいことです。しかし、異なる空のパッチで加算/減算できるようにしたいので、異なるディスクを同じサイズにすることができないという事実が問題を引き起こします。

私も使用しました

hp.gnomview(return_projected_map=True)

私の空のパッチを見て、返された各 2D 配列が同じサイズであることがわかりました。ただし、おそらく多くのサブマップで作業を行うため、毎回マップを描画する必要はありません。二次元配列のみ取得することはできますか？もしそうなら、plt.imshowそれらの配列を加算/減算した後に使用しても問題ありませんか?

私が望むものを達成するためのよりエレガントな方法があるべきだと思います-私が見落とした既存のルーチンはありますか? ありがとう！

AKARI Far Infrared Surveyor データベース (公開) からの HEALPix 全天マップを持っています。を使用してマップを「滑らかに」しようとしましhealpyたが、結果は非常に奇妙に見えます。より良い方法はありますか？ただし、私の質問は、全天 HEALPix マップ (つまり、IRAS、Planck、WISE、WMAP) に関するものです。

私の目的は、この「あかり」マップの実効点広がり関数を 1 度の角度分解能に「滑らかにする」ことです (元のデータの PSF は約 1 分角です)。これは、遠赤外線「あかり」マップを低解像度のマイクロ波マップ (特に異常なマイクロ波前景のもの) と比較できるようにするためです。

以下の例では、劣化したバージョンのマップを使用しているため、Github にアップロードするのに十分な大きさです。これは、ピクセルが約 3.42 分角であることを意味します。通常、PSF スムージングの前に、ピクセル スケールをそれほど劣化させることはありませんが、これは単なる例です。

私は、healpy.sphtfunct.smoothing ルーチン ( https://healpy.readthedocs.org/en/latest/generated/healpy.sphtfunc.smoothing.html#healpy.sphtfunc.smoothingsmoothing ) を試しました。球面調和関数にマップし、ガウスと畳み込み、空間マップに変換します。

ipython ノートブックと低解像度の .FITS HEALpix マップを github リポジトリに保存しました。

https://github.com/aaroncnb/healpy_smoothing_test

healpy(パッケージをインストールする必要があります)

ノートブックでコードを実行すると、私が抱えている問題を簡単に視覚化できます。マップを平滑化した後、ピクセルが円形のガシアンで平滑化されるのではなく、ボックス平均化が繰り返されたかのように、いくつかの奇妙な「アーティファクト」があります。プロフィール。私が期待しているのは、入力マップのぼやけたバージョンです。

平滑化が行われる前に、球面調和関数への変換に関する基本的な何かが欠けていると思います。

HEALPixマップで、この種の全天スムージングを以前に試みた人はいますか?

別のオプションは、マップを標準の長方形配列に変換してから、平滑化を行うことだと思います。しかし、私は HEALPix 形式を離れずに問題を解決することに興味があります。

頂点によって定義される幾何学的領域 (球面三角形、球面多角形など) に属するピクセル インデックスを見つけるように設計された Healpix IDL ルーチンが多数あります。query_*ルーチンがあります（例query_triangle）。こちらのドキュメントを参照してください。

http://healpix.jpl.nasa.gov/html/idlnode45.htm

Healpy プログラムでこれらのピクセル インデックスを使用したいと考えています。また

query_*(A) IDLルーチンからのピクセル インデックスの出力リストをdata_file.save形式で保存できました。.sav次に、さまざまなモジュールを使用して、このピクセル インデックスのファイルを Python にインポートできます。たとえば、 http ://www.astropython.org/packages/idlsave94/

(B) どうにかして IDL をまったく使用しない方がはるかに便利です。query_*Healpy にはピクセル関連の関数がいくつかありますが、healpy だけを使用して IDL ルーチンを「変換」する方法はないようです。

query_polygonヘルピーを使う方法はありますか？これを行うことは可能ですか？

Healpy.mollview で他のカラーマップを使用しようとしています このコードで成功しました

しかし、予期しない青い背景が表示され、それを白に設定する方法がありません

コードのデバッグを試みましたが、AltAz 座標を .csv ファイルに保存しようとすると、最終的には壊れることに気付きました。これは、numpy 配列ではなく、SkyCoord オブジェクトであるためです。赤道座標の大きなテーブルを AltAz に変換する簡単な方法や、コードをファイルに保存する方法を誰かが提案できますか?