Projで座標変換

ある１点の緯度経度から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も依存関係で入る．