緯度経度を平面座標に変換する方法（北米）

北米において、緯度経度を平面座標に変換する方法を知りたいです。
緯度経度が二組与えられた時の、二直線の交点を求めたいと思っています。

数十メートルと短い距離なので、平面座標に置き換えることを考え、
直角平面座標を調べてpythonでの実装例も見つけましたが、
これは日本国内にしか使えないようです。
https://sw1227.hatenablog.com/entry/2018/11/30/200702

北米の緯度経度を平面座標に変換するための式と、できれば実装例を教えていただけませんか。

hoshi-takanori

2020/04/07 07:40

平面直角座標系の境界線ってかなり恣意的に設定されているように見えるので、よその国に勝手に適用したら顰蹙なのでは。 https://www.gsi.go.jp/sokuchikijun/jpc.html

行動規範の内容に同意します

回答2件

ベストアンサー

座標変換できるpyprojがあります。

pyproj

transform関数で座標変換できます。

緯度経度の投影法は何で(WGS84,EPSG4326?)、どの投影法の平面座標に変換したいのかを、確認する必要があります。

投稿2020/04/16 12:14

hasami

総合スコア1277

Satomi07

2020/04/19 04:14

下記サイトで使い方と実装例をみつけることができました、ありがとうございます。 https://ikatakos.com/pot/programming/python/packages/pyproj 英語のもので調べるのはGIS初心者の私には難しかったですが、下記を参考に、北米のコードを割り出すことができました。 https://club.informatix.co.jp/?p=1225

行動規範の内容に同意します

平面直角座標系の大まかなイメージとしては、ある地点を原点とし、その近辺で地球が平坦であるとみなしたものであると考えれば良いでしょう。もちろん地球の楕円体を歪みなく平面に変換することは不可能ですが、なるべく「ずれ」が無いように工夫してあります。

ただ、日本全土のスケールになってくると地球の丸みは無視できないものになるため、「ずれ」を抑えるためには原点を複数点用意する必要が生じます。実際、日本国内では平成十四年国土交通省告示第九号によって19の平面直角座標系が定められており、それぞれの座標系原点の経緯度と適用範囲が定められています。平面直角座標はこの適用範囲それぞれの内部においてのみ比較可能なものであることに留意しましょう。*2

同じことを北米で行えばいいだけでは？面積から考えて基準点が19箇所では足りないでしょうが。
というかこれ別にPython関係ないですよね…？

投稿2020/04/04 05:12