自分で計算せずに地磁気座標を地理座標に変換するにはどうすればよいですか？

地磁気座標のデータがあります。つまり、緯度と経度ですが、地理的な北極ではなく、カナダの磁気北極を参照しています。ご想像のとおり、他の座標系にある他のデータと照合します。

（データはオーロラオーバルの計算として作成されたものです。つまり、それらは地磁気として作成されたものであり、原則として私が相談できる地理的な情報源はありません。）

データはPostGISに保存されます。

これを行う簡単な方法は、地磁気座標系のEPSGコードを見つけて、そのコードでデータにタグを付けることです。その後、すべてが魔法のようにそのときに必要なものに再投影されます。しかし、地磁気または磁気のいずれかを含むEPSGコードをspatialreference.orgで見つけることができません。さまざまなGoogle検索でも何も見つかりません。

必要な精度のレベルについては、異なる北極を持つように変更された緯度/経度座標系（たとえば、WGS84）で十分です。したがって、そのような座標系を定義してPostGISにロードする簡単な方法があれば、それも問題ありません。

明確化

たとえば、地磁気座標から地理座標に変換するのではなく、単にUTM 15NからWGS84に変換したいとします。UTMオブジェクトを変更せずにPostGISに保存し、EPSGコード26915でタグ付けすることができます。その後、後でクエリを実行すると、PostGISに「結果をEPSGコード4326で返してください」と言うことができ、PostGISは自動的にWGS84に変換します。

同様に、地磁気データに何らかの座標系をタグ付けし、reproject()関数を手動で呼び出さなくても、PostGISにバックグラウンドでの再投影を実行させたいと思います。

オープンソースがこの地磁気から地理をサポートしているとは思えません。

しかし、座標が少ない場合は、pls。これを試して
ください。http：//wdc.kugi.kyoto-u.ac.jp/igrf/gggm/index.html

NASAがアルゴリズムを公開しました。http：//idlastro.gsfc.nasa.gov/ftp/pro/astro/mag2geo.proを試すことができ ます

それはクリアvadivelanメイクによって提供されるリンク地磁気座標は球だけ回転しているように（）、N極が現在の磁気通過南物理的及び磁極を通って磁極と、（b）本初子午線通過します。

たとえば、IGRF-95は、1995年の緯度79.30および経度-71.41度での磁北極の位置を使用します。

球面回転は、さまざまな方法で計算できます。通常、3Dの地理中心のデカルト座標に変換し、線形変換（3 x 3行列乗算または四元数乗算による）を適用して、球面座標に変換します。しかし、これをプログラムしたくない場合は、トリックがあります。斜め投影をサポートするGISは、再投影によって回転できます。地磁気座標をあたかも地理座標であるかのように投影します。たとえば、ステレオ投影法の極面（中央子午線が0度）を使用します。次に、同じ投影法の斜めのアスペクトを使用して、起源の緯度が79.30で中央子午線が-71.41である結果を投影解除します。最後に、目的の投影を使用してもう一度投影し、「魔法の再投影」を取得します。