pythonで緯度経度を平面直角座標に変換したい

Question

緯度経度→平面直角座標に変換するプログラムを作成したいです。

当方、プログラミングの知識が乏しく自力で書けなかったため、知恵をお借りさせていただければと思います。

使用言語はpython3.5です。

何卒、よろしくお願いします。

Comments

[matobaa]

平面直角座標に変換する式をご存じであれば教えてください

[sekaikan_ozaki]

http://vldb.gsi.go.jp/sokuchi/surveycalc/surveycalc/algorithm/bl2xy/bl2xy.htm　こちらが計算式となっています。

Answer 1

pyproj を使用して座標変換するのが簡単かと思います。

http://jswhit.github.io/pyproj/

import pyproj

# 東京タワーの位置
lat,lon = 35.658581, 139.745433

# 変換元の座標系：緯度経度
p1 = pyproj.Proj(init='epsg:4612')
# 変換先の座標系：平面直角座標系 ９系
p2 = pyproj.Proj(init='epsg:2451')

# 変換
y,x = pyproj.transform(p1, p2, lon, lat)

# 結果表示
print(x,y)
# -37875.0849492072 -7958.643508325489

Comments

[sekaikan_ozaki]

error: command 'cl.exe' failed: No such file or directoryというエラーメッセージが出てしまい、インストールすることができませんでした。
pip以外で、pyprojをインストールする方法はありますでしょうか？