とらりもんHOME  Index  Search  Changes  Login

PROJ.4

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

コマンド

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)
Last modified:2020/03/16 07:27:06
Keyword(s):
References:[とらりもんHOME]