【python】WGS84から平面直角座標系への変換 - Qiita

はじめに 測量では成果を緯度経度でなく、平面直角座標系という２次元の平面の座標値で表現します。その相互変換を行う方法のメモです。 日本では現在、19の直角座標系が用いられています。(「分かりやすい平面直角座標系」from 国土地理院) よくある(?)誤解: 慣性航法で使う座標系とは違う！ 座標系の変換というと簡単に聞こえますが、これは WGS84の楕円体の極座標(緯度,経度,楕円体高)を地球中心の座標系(ECEF)の直交座標系の座標値XYZに変換 ECEFの直交座標系の値を、該当する楕円体における接平面に対応するように回転させる という話とは『異なり』ます。ガウス・クリューゲル投影法という計算方法を使うようなのですが、私は数式が複雑で着いていけませんでした。 当初、自分は上記のような単純な回転計算と誤解していて、周りに「そんなの簡単な計算だよ」などと話していて、思い返すと恥ずかしいです。