GMT gyorstalpaló¶

A GMT az eszközök gyűjteménye, melyek lehetővé teszik (x,y) és (x,y,z) adathalmazok kezelését (beleértve a szűrést, a trend illesztést, a rács készítést, vetület váltást, stb.) és Encapsulated PostScript File (EPS) illusztrációk készítését az egyszerű x-y rajzoktól, a szintvonalas térképeken át, a mesterségesen megvilágított felületekig és háromdimenziós perspektivikus nézetekig fekete-fehér, szürke árnyalatos, speciális minták és 24 bites színekkel.

Ebben a rövid GMT bemutatóban digitális terepmodellt (DTM) fogunk használni, hogy árnyalt domborzati térképet hozzunk létre.

Tartalom

Térkép létrehozása
Színpaletta létrehozása
Gradiens meghatározása
Színskála megjelenítése
Kipróbálandó dolgok
Mi a következő?

Térkép létrehozása ¶

Input adathalmazként egy DTM-et fogunk használni a GRASS GIS példa adatok közül, mely kezelhető a GDAL-lal a GDAL-GRASS meghajtóval.

A DTM metaadatok gyors megtekintéséhez nyisson egy terminál ablakot és használjuk a gdalinfo GDAL parancsot:

gdalinfo /home/user/grassdata/nc_basic_spm_grass7/PERMANENT/cellhd/elevation

Láthatjuk, hogy a DTM vetületi koordináta-rendszerben van és az adatok terjedelme 50 és 160 méter közötti.

A DTM-et egy GMT kompatibilis netCDF fájlformátumba konvertálhatjuk a gdal_translate paranccsal így:

gdal_translate -of GMT /home/user/grassdata/nc_basic_spm_grass7/PERMANENT/cellhd/elevation elevation.gmt

A GMT példánkhoz először a az adathalmazt földrajzi koordinátákba alakítjuk át, hosszúság-szélesség (EPSG:4326), a gdalwarp parancsot használva:

gdalwarp elevation.gmt geo_elevation.gmt -t_srs EPSG:4326

Színpaletta létrehozása ¶

A rácsunkhoz egy szín beállítás hozzárendeléséhez egy GMT színpalettát hozhatunk létre a makecpt paranccsal. A színek tartományát és számát a -T kapcsolóval állítjuk be és a -C kapcsolóval rendeljük hozzá a színeket.

gmt makecpt -Chaxby -T50/160/10 -Z > elevation.cpt

Kísérletezhet más színtáblákkal, lásd a man gmt makecpt parancsot a teljes listához.

Gradiens meghatározása ¶

Az árnyékolt domborzat példánkhoz meg kell adnunk hogyan változik a szín intenzitás a rács mentén. Ehhez kiszámoljuk az irány szerinti deriváltat vagy gradienst a rácsunkból a grdgradient parancs segítségével. A grdgradient parancshoz használjuk a -Ne kapcsolót az output normalizálásához a kumulatív Laplace elosztást használva és a -A kapcsolót az irány szerinti deriváltak északtól mért irányához.

gmt grdgradient geo_elevation.gmt -Ne0.8 -A100 -fg -Ggradient.nc

Az output egy netCDF fájl, mely a következőképpen használható mint szín intenzitás a grdimage parancshoz:

gmt grdimage geo_elevation.gmt -Igradient.nc -JM6i -P -Ba -Celevation.cpt -K > GMT_tut_OSGeoLive.ps

Színskála megjelenítése ¶

Végül egy színskálát hozzáadhatunk a psscale paranccsal:

gmt psscale -DjTC+w5i/0.25i+h+o0/-1i -Rgeo_elevation.gmt -J -Celevation.cpt -I0.4 -By+lm -O >> GMT_tut_OSGeoLive.ps

Itt a -D kapcsolót használjuk a színskála pozíciójának kiválasztásához TC: felül középen, a méret 5x0.25 hüvelykben i kifejezve a tájolás h: vizszintes és a -By állítja be az m címkét az y tengelyen.

Az eredmény ps dokumentumnak így kellene kinéznie: