Knihovna PROJ.4

• Knihovna proj se stará ve většině Open Source GIS projektů o práci se souřadnými systémy.
• Součástí její distribuce jsou i konzolové programy sloužící k převodu souřadnic či rastrových nebo vektorových dat z jednoho systému

na druhý

• http://proj.maptools.org

Výpis podporovaných mapových projekcí:

GRASS> cs2cs -lp

Podporované parametry konkrétní projekce:

GRASS> cs2cs -l=utm
GRASS> cs2cs -l=krovak

Podporované elipsoidy:

GRASS> cs2cs -le

Podporované datumy:

GRASS> cs2cs -ld

Převod S-JTSK -> WGS84:

GRASS> echo "-734542.9904476  -1045586.42733531" |cs2cs +proj=krovak
+ellps=bessel +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +to +init=epsg:4326

    14d32'26.049"E  50d4'28.608"N 45.158

GRASS> echo "-734542.9904476  -1045586.42733531" |cs2cs +init=epsg:2065 +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +to +init=epsg:4326

    14d32'26.049"E  50d4'28.608"N 45.158

GRASS> echo "-734542.9904476  -1045586.42733531" |cs2cs -f "\%f" +init=epsg:2065 +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +to +init=epsg:4326

    14.540569       50.074613 45.158022

Převod souřadnic v souboru:

GRASS> cat soubor.txt|cs2cs ...

Kódy EPSG

• European Petroleum Survey Group

• http://epsg.org

• Soubor s definicemi souřadných systémů /usr/share/proj/epsg nebo /usr/local/share/proj:

  # Unknown datum based upon the Airy 1830 ellipsoid
  <4001> +proj=longlat +ellps=airy +nodefs  <>
  # Unknown datum based upon the Airy Modified 1849 ellipsoid
  <4002> +proj=longlat +a=6377340.189 +b=6356034.447938534 +nodefs  <>
  # Unknown datum based upon the Australian National Spheroid
  <4003> +proj=longlat +ellps=austSA +nodefs  <>
  ...
  # WGS 84
  <4326> +proj=longlat +ellps=WGS84 +datum=WGS84 +nodefs  <>

• S-JTSK mají v názvu tři definice (2007):

  # S-JTSK
  <4156> +proj=longlat +ellps=bessel +nodefs  <>
  # S-JTSK (Ferro)
  <4818> +proj=longlat +ellps=bessel +pm=ferro +nodefs  <>
  # S-JTSK (Ferro) / Krovak
  <2065> +proj=krovak +lat0=49.5 +lon0=42.5 +alpha=30.28813972222222 +k=0.9999 +x0=0 +y0=0 +ellps=bessel +pm=ferro +units=m +nodefs  <>

• Žádný EPSG kód s +proj=krovak neobsahuje transformační parametry na elipsoid WGS84. Evropské datumy lze nalézt na http://crs.ifag.de nebo v souboru $GISBASE/etc/datumtransform.table

Vlastní S-JTSK jako EPSG

Vložte na začátek souboru /usr/share/proj/epsg řádky:

# S-JTSK
<102067> +proj=krovak +lat_0=49.5 +lon_0=24.83333333333333 +alpha=30.28813975277->778 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel +units=m +towgs84=498.17,136.89,510.->08,6.007,4.343,3.831,3.38 no_defs <>

Nyní můžete jednoduše převádět:

echo "-734542.9904476  -1045586.42733531" |cs2cs +init=epsg:102067 +to +init=epsg:4326

Nebo transformovat rastrové a vektorové soubory, nebo zakládat GRASS Location z kódu EPSG.