Tvorba map pomocí GRASSu

28.05.2006 – 3:06 pm

Tak posledního asi půl roku se bavím klasifikací několika leteckých snímků z
Krkonošského národního parku. V pondělí odjíždíme do Krkonoš, ověřit tyhle
výsledky v terénu. To jsem žádostiv, jak moc jsem se (ne)trefil.

Důvod proč tohle ale píšu je, abych zase jednou ukázal, jak lze i v GRASSu
vytvářet mapové výstupy. No řekněte, není tohle hezká mapa se vším, co tam má
být? Dokonce i s češtinou :-)

Výsledek je poněkud zmenšený originál (na 81%). Originál je PDF soubor na
formát papíru A3.

Poustup, jak takhle mapa vznikla:

  1. Je potřeba připravit si rastrovou mapu. Protože modul pro tvorbu map v
    GRASSu ps.map neumí zobrazovat víc, než jednu rastrovou mapu,
    je potřeba si ji vyrobit předem.

    Na tomto obrázku konkrétně je to mapa s výsledkem klasifikace, s odstíny
    šedi — nejčastěji se používá hillshading, já ale použil jeden z kanálů z
    barevných leteckých snímků. Tyto dvě rastrové mapy musíte složit do jedné:

    r.his h=klasifikace i=snimky.red r=tisk.r b=tisk.b g=tisk.g
    # výsledek se musí složit pomocí r.composite
    r.composite r=tisk.r b=tisk.b g=tisk.g out=tisk
    

    Tím máme vytvořenou mapu tisk, která je vystínovaná a můžeme jí použít v
    ps.map.

  2. Musíte si připravit případné labely, které chcete mít v mapě. Na to je
    modul v.label. Na vrstevnice doporučuji zapřáhnout parametr

    -a, který vytiskne písmenka ,,podél“ vrstevnic.

  3. Připravíte si soubory tisk.hdr, obsahující nadpis mapy
    (,,Klasifikace území Krko…“), a comments.txt, obsahující
    komentář k mapě (,,(c)Ústav geoinformačních technologií…“). Je důležité,
    mít správně vyplněný soubor
    $GISDABASE/$LOCATION_NAME/$MAPSET/cats/klasifikace, aby se
    správně zobrazila legenda k rastrové mapě. Ačkoliv je vytištěna mapa
    tisk, legenda patří k mapě klasifikace.

  4. Nakonec je potřeba všechno pěkně napsat do správného souboru. Člověk
    si chvíli hraje s pozicemi… to je asi nejhorší, ale nakonec to jde.

  5. A teď to nejlepší: po spuštění ps.map, proženete výsledný
    soubor prográmkem ogontify, který převede Post Scriptový
    soubor na kódování iso-8859-2 — tedy české znaky se stanou českými.
    ogontify ale nenastavuje jednu drobnost, proto se to udělá
    nejlépe přes sed, takže:

    ogonkify -eL2 tisk.eps |sed -e 
    s/ISOLatin1Encoding/ISOLatin2Encoding/ > tisk_cesky.eps
    


    Prehledka map
Krkonoš

    A na úplný konec převod do PDFka:

    ps2pdf -sPAPERSIZE=a3 tisk_cesky.eps
    

    A je to!

Napsal jsem skriptík na tvorbu takového atlasu map – rozřeže aktuální region
na oblasti podle zadaného počtu sloupců a řádků a překrytí a vyrobí mapy podle
standardního souboru pro ps.map. Kdyby někdo měl zájem, viz stránka s moduly pro GRASS.

Share Button

Post a Comment