2014/08/29

Projで座標変換

ある1点の緯度経度からUTMに変換する方法.

proj.4 に含まれている cs2cs で出来る

% echo "140.0 36.0" | cs2cs +init=EPSG:4326 +to +init=EPSG:32654 -f %.4f
409870.9506     3984410.7881 0.0000

入力データの座標系(測地系+投影図法)は,EPSG:4326でWGS84の緯度経度座標.
出力データの座標系は,EPSG:32654でWGS84のUTM 54N.

EPSGコードは,proj.4 パッケージの epsg ファイルで確認できる

% grep -B 1 32654 /opt/local/share/proj/epsg
# WGS 84 / UTM zone 54N
<32654> +proj=utm +zone=54 +datum=WGS84 +units=m +no_defs  <>

ちなみに,ここ(Google Maps)

参考:
http://ryuiki.agbi.tsukuba.ac.jp/~nishida/lecture/GIS/proj.html

余談:
Macの場合,Ports に入っているのでついでによく使うgdalも入れちゃうのが良い.

% port install gdal

で,projも依存関係で入る.