Les-ejk

Tvorba map pomocí GRASSu

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
   

   
   
   

   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.