***************
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ů :file:`/usr/share/proj/epsg` nebo :file:`/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 :file:`$GISBASE/etc/datumtransform.table`

.. _102067-label:
========================
Vlastní S-JTSK jako EPSG
========================
Vložte na začátek souboru :file:`/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.