とらりもんHOME  Index  Search  Changes  Login

とらりもん - PROJ.4 Diff

  • Added parts are displayed like this.
  • Deleted parts are displayed like this.

参考 PROJ.4 HP  [[http://trac.osgeo.org/proj/]]

!!!!!コマンド

!!!!!cs2cs!!cs2cs (2010/04/14 松吉)

測地系と投影法の変換をしてくれます。

EPSGコードを使ったUTM座標 (zone54) から緯度経度への変換例 (測地系wgs84)。

$ echo "139398.30 3991354.13 0.00" | cs2cs -v +init=epsg:32654 +to +init=epsg:4612

# ---- From Coordinate System ----
#Universal Transverse Mercator (UTM)
# Cyl, Sph
# zone= south
# +init=epsg:32654 +proj=utm +zone=54 +datum=WGS84 +units=m +no_defs
# +ellps=WGS84 +towgs84=0,0,0
# ---- To Coordinate System ----
#Lat/long (Geodetic alias)
#
# +init=epsg:4612 +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
137dE 36dN 0.000



*EPSGコードは様々な測地系や投影法をまとめてID番号をつけたリスト。
*EPSGコードのリストは以下の場所に保存されてます。
$ less /usr/share/proj/epsg

::例えば…

::測地系: wgs84 投影法: 緯度経度 EPSGコード: 4612
::測地系: wgs84 投影法: UTM (zone??, 北半球) EPSGコード: 326??
::測地系: wgs84 投影法: UTM (zone??, 南半球) EPSGコード: 327??

!Python
from pyproj import Proj, transform
x1, y1 = 139398.30, 3991354.13
x2,y2 = transform(Proj(init='epsg:32654') ,Proj(init='epsg:4612'), x1, y1)

(136.99999997910064, 35.99999997210519)