
2026/04/25 22:11
HEALPix(ヒールピクス)
RSS: https://news.ycombinator.com/rss
要約▶
Japanese Translation:
HEALPix(Hierarchical Equal Area isoLatitude Pixelisation)は、1997年に Krzysztof M. Górski によって考案された、天球を等積のピクセルに分割するためのアルゴリズムである。ランベル円筒等積投影(赤道座標系)と中断コリニョン投影(極座標系)を組み合わせて 2-球面をユークリッド平面へ写すことで、H=4、K=3 の場合に完全に立方体へと畳み込まれるような柔軟な階層構造を可能にする。1998 年にプレプリントとして初公開され、特に宇宙マイクロ波背景放射マップの分野を含む物理的宇宙論において広く利用されてきた手法は、2006 年 4 月 26 日に国際天文学連合 FITS ワーキンググループによって HPX というキーワード名の下で FITS 世界座標系(WCS)の一部として公式に採択された。HTM(階層三角メッシュ)などの代替手段が存在するにもかかわらず、HEALPix はガイヤ計画の光源同定を含む主要なミッションの基盤技術として引き続き重要な役割を果たしており、C、C++、Fortran、IDL、Java、Python といった言語での公式実装が利用可能で、解像度は最大 0.4 mas(Java の移植版では 0.3 角秒に達する)にまで向上している。Astropy など、healpy および astropy-healpix を介して統合されたフレームワークにおいて、HEALPix は大規模な天文データセットの高解像度解析に必要な計算基盤を提供し続けている。
本文
HEALPix(時に「Healpix」と表記される)とは、2 次元球面に対する階層的かつ等積の緯度準直線画素化(Hierarchical Equal Area isoLatitude Pixelisation of a 2-sphere)の略称であり、2 次元球面の画素化アルゴリズムおよび関連する図法クラスの名称である [1]。この画素化アルゴリズムは、1997 年にデンマークのコペンハーゲンに所在した理論天体物理学センター(Theoretical Astrophysics Center)にて、Krishnótf M. Górski が考案し、1998 年に予習論文として初めて公開された [2][3][4]。
投影法と画素化 [編集]
HEALPix 投影法は、2 次元球面をユークリッド平面へ写像する複数の主要な性質を共有する、一般化された球面投影法の一種である [1]。これらのいずれについても、その後得られた 2 次元平面上の領域に対して、分割(画素化)を適用することができる。特に、上記の投影法のうち H=4 および K=3 の HEALPix 投影法に続き、その結果に対して平面の画素化を施したものを「HEALPix 画素化」と一般的に呼称している [3][4]。この手法は、宇宙マイクロ波背景放射(CMB)のマッピングにおいて物理宇宙論で広く用いられている。本画素化は、球面を平面上の 12 つの正方形ファセット(菱形)へ写像し、その後これらを 2 進的に分割してピクセルとするプロセスと捉えることもできる [5][6][1]。ただし、投影法を用いずにこの構造を導出することも可能である [3][4][7]。関連するソフトウェアパッケージ HEALPix は当該アルゴリズムを実装しており、天文学データファイルの記述に用いられる FITS 標準におけるキーワード「HPX」には、一般化された球面投影法としての HEALPix 投影法が対応している。国際天文連合(IAU)の FITS ワーキンググループは、2006 年 4 月 26 日にこれを公式な FITS 世界座標系(WCS)の一部として承認した [8]。
- 本球面投影法は、赤道近傍領域にはランベルト円筒等積投影(Lambert cylindrical equal-area projection)、極域には準円筒等積かつ断続的なクロニョン投影(interrupted Collignon projection)を組み合わせて構成されている [1]。
- 階層構造のある特定レベルにおける画素はすべて同面積を持つ(H=4、K=3 の HEALPix 投影法においては、正方形を二等分することでこれを達成している)。また、これらの画素の中心点は離散的な数の緯度円上に位置し、各円上では均等間隔で配置されている。このスキームは数学的性質に富み、球面調和関数変換などの特定の計算において効率的である。H=4、K=3 の HEALPix 投影法の案例においては、平面上の画素は正方形であり(これは逆投影操作によって、2 次元球面上で非測地線を持つ四角形へと復元可能)、各頂点に 4 つの画素が集まるが、そのうち 8 つの頂点例外ではそれぞれ 3 つの画素のみが存在する。
- 赤道正交経線と極収束経線との間の遷移緯度は、この投影法を「完全な立方体」へと折りたたむことを可能にするよう選択されている——すなわち、「球体の立方体化(cubing the sphere)」を実現しており、特にこの手法により北極圏が正方形形状となる。
利用法と代替手段 [編集]
H=4 および K=3 投影法に関連する画素化は、宇宙マイクロ波背景放射マッピングの保存および操作において、物理宇宙論で広く採用されている。
- ガイア(Gaia)ミッションでは、天体源の識別のための基礎として HEALPix を用いている [9]。
- 代替となる階層的グリッドとしては、ヒierarchical トリアングル・メッシュ(HTM)がある [10][11]。この場合も階層の特定のレベルにおける画素は類似した大きさを有するものの、完全に同一ではない。本スキームは境界がすべて球面の円弧から構成されるため、複雑な形状の表現に適している。
- その他にも、四角形化した球面立方体(Quadrilateralized Spherical Cube)という別の階層的グリッドが存在する。
- 多くの他の球面グリッドシステムとは異なり [12]、HEALPix は球面を多面体へ写像することに基づいていない。具体的には、標準的な H=4、K=3 の HEALPix 投影法は菱面体十二面体に似た 12 つの四角形面で構成されるが、その頂点配置は異なり、事実上 genus 0 の多面体とは互換性がない [13]。
- H=6 の HEALPIX は、二十面体に基づく別の代替グリッドと類似性を示す [14]。
図法一覧 空間グリッド/測地線グリッド
文献: a, b, c, d Calabretta, Mark R.; Roukema, Boudewijn F. (2007). "Mapping on the HEALPix grid". Monthly Notices of the Royal Astronomical Society 381 (2). Oxford University Press: 865–872. Bibcode:2007MNRAS.381..865C. doi:10.1111/j.1365-2966.2007.12297.x.
^ "HEALPix Background - History". healpix.jpl.nasa.gov. 2019 年 6 月 8 日閲覧。
^ a, b, c, d Górski, Krzysztof M.; Hivon, Éric; Wandelt, Benjamin D. (1999). "Analysis Issues for Large CMB Data Sets". Evolution of Large Scale Structure: From Recombination to Garching. Proceedings of the MAP/ESO Cosmology Conference 'Evolution of Large-Scale Structure'. Netherlands: PrintPartners Ipskamp. p. 37. arXiv:astro-ph/9812350. Bibcode:1999elss.conf...37G.
^ a, b Górski, Krzysztof M.; Wandelt, Benjamin D.; Hansen, Frode K.; Hivon, Éric; Banday, Anthony J. (1999-05-21). "The HEALPix Primer". arXiv:astro-ph/9905275.
^ Roukema, Boudewijn F.; Lew, Bartosz (2004-09-08). "A Solution to the Isolatitude, Equi-area, Hierarchical Pixel-Coordinate System". Public Draft. arXiv:astro-ph/0409533. Bibcode:2004astro.ph..9533R. 2019 年 8 月 4 日にオリジナルアーカイブより削除。2004 年 9 月 8 日閲覧。
^ Roukema, Boudewijn F.; Lew, Bartosz (2004-09-22). "A Solution to the Isolatitude, Equi-area, Hierarchical Pixel-Coordinate System". arXiv:astro-ph/0409533.
^ a, b Górski, Krzysztof M.; Hivon, Éric; Banday, Anthony J.; Hansen, Frode K.; Wandelt, Benjamin D.; Reinecke, M.; Bartelmann, M. (2005). "HEALPix: A Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere". The Astrophysical Journal 622 (2): 759–771. arXiv:astro-ph/0409513. Bibcode:2005ApJ...622..759G. doi:10.1086/427976. S2CID 18743679.
^ Pence, William D. "FITS World Coordinate System (WCS)". High Energy Astrophysics Science Archive Research Center (HEASARC). 2019 年 8 月 4 日にオリジナルアーカイブより削除。2007 年 1 月 9 日閲覧。
^ "Gaia Data Release 1: Datamodel description Documentation release 1.2". gea.esac.esa.int. 2021 年 5 月 31 日閲覧。
^ "SkyServer.org - HTM: Hierarchical Triangular Mesh". SkyServer. 2006 年 6 月 6 日。2007 年 2 月 5 日閲覧。
^ Szalay, Alex; Jim Gray; Gyorgy Fekete; Peter Kunszt; Peter Kukol; Ani Thakar (September 2005). "Indexing the Sphere with the Hierarchical Triangular Mesh". Microsoft Research. arXiv:cs/0701164. Bibcode:2007cs........1164S. 2007 年 2 月 5 日閲覧。
^ Kevin Sahr, Denis White, and A. Jon Kimerling (2003). "Geodesic Discrete Global Grid Systems" (PDF). Cartography and Geographic Information Science 30 (2): 121–134. {{cite journal}}: CS1 メンテナンス – multiple names: authors list
^ オイラーの多面体公式を用いて確認可能である。
^ Tegmark, Max. "Welcome to the icosahedron home page". space.mit.edu.
多言語対応を備えた公式実装(C、C++、Fortran90、IDL、Java、Python)、分解能は最高で 0.4 ミリオング(mas)まで対応可能:
- オリジナルの Fortran コードによる Java ポート版(Nikolay Kuropatkin 氏作成)、分解能は最高で 0.3 角秒まで対応。
- RangeSet を活用するように最適化された Java 実装、高解像度処理に特に優れる。
- healpy: Python ウラッパーパッケージ。
- astropy-healpix: Astropy 用 BSD ライセンス下での HEALPix 実装。
- healpix.cxx: HEALPix 座標変換用の C++ コード。
- JavaScript ポート版: JavaScript で経緯度と HEALPix 座標の相互変換を行う。
- Typescript healpix: JavaScript/TypeScript による HEALPix 実装。
- Healpix.jl: Julia 言語で記述された HEALPix ライブラリ。