ELTE Meteorológiai Tanszék => Oktatás => GMT oktatóanyag

GMT Example 14


A script megtekinthető és a postscript fájl letölthető a GMT oldalán az  almenünél vagy itt: (35.tar - UNIX tar formátumban, illetve alternatív módon 35.zip - WinZip tömörítéssel)

A GMT jónéhány rácskészítő algoritmust támogat. Az alábbiakban a GMT 14-es példájához talál segítséget. A parancsok és egyes kapcsolóik magyarázata csak egyszer található meg. Ha nem találja valamely parancs vagy paraméter magyarázatát akkor az azért van, mert már korábban sor került erre.



A script:

gmtset GRID_PEN_PRIMARY thinnest,-

A gmtset parancs segítségével a GMT alap beállítását módosíthajuk.

Az alapparancs utáni paraméterek jelentése:

paraméter
jelentés
GRID_PEN_PRIMARY thinnest,-

Vonalvastagság beállítása, itt a legvékonyabb.


psxy table_5.11 -R0/7/0/7 -JX3.06i/3.15i -B2f1WSNe -Sc0.05i -Gblack -P -K -Y6.45i >! example_14.ps

A psxy parancs (x,y) párokat olvas ki fájlból vagy valamilyen standard inputból és PostScript kódot generál belőlük és az így kapott helyeken vonalakat, poligonokat vagy szimbólumokat rajzol a térképre, ábrára.

Az alapparancs utáni paraméterek jelentése:

paraméter
jelentés
table_5.11 Fájl neve, amiből az (x,y) párok kiolvasásra kerülnek. Ha nincs megadva fájlnév akkor a standard inputból történik az olvasás.
-R0/7/0/7 A -R tartomány beállítására szolgáló paraméter (xmin/xmax/ymin/ymax), földrajzi területek esetén a határok rendben: nyugat, kelet, dél és észak
-JX3.06i/3.15i A -J paraméter a térképvetület beállítására szolgál. Itt a -JX egy  nem földrajzi vetületet jelöl, amin az x tengely 3,06 inch, az y tengely 3,15 inch.
-B2f1WSNe A B paraméterrel a térképet határoló vonalak jelölését lehet beállítani. Itt azt állítjuk be, hogy kettesével legyen egy jelölés és csak az égtájak szerint vagyis bal oldalon 'W', felül 'N', alul 'S', az 'e' pedig azt, hogy a jobb oldalon ne legyen jelölés
-Sc0.05i A -S szimbólumokat rajzol, itt az -Sc0.05i azt jelenti, hogy a szimbólum kör és az átmérője 0,05 inch
-Gblack A -G a szimbólum szinét állítja be, a 'black' egyértelmű: fekete
-P Portrait mód, vagyis a lapméret legyen álló.
-K Nagyon fontos beállítás: azt, jelöli, hogy a készülő PostScript fájl nem zárható le, még további információ kerül az állományba (jelen esetben a example_14.ps-be).


blockmean table_5.11 -R0/7/0/7 -I1 >! mean.xyz


psbasemap -R0.5/7.5/0.5/7.5 -J -O -K -B0g1 -X3.25i >> example_14.ps

A psbasemap parancs egy alaptérképet készít, amire később tetszés szerint helyezhetünk további tartalmat (szimbólumokat, országhatárokat, jelmagyarázatot, feliratokat, stb.).

Az alapparancs utáni paraméterek jelentése:

paraméter
jelentés
-R0.5/7.5/0.5/7.5 A -R tartomány beállítására szolgáló paraméter (xmin/xmax/ymin/ymax), földrajzi területek esetén a határok rendben: nyugat, kelet, dél és észak
-J A -J a térképvetületet adja meg.
-O A létező koordinátarendszert használjuk, nem definiálunk új térképet. Alapértelmezés szerint új térkép készül.
-B0g1 A B paraméterrel a térképet határoló vonalak jelölését lehet beállítani.

psxy -R0/7/0/7 -J -B2f1eSNw mean.xyz -Ss0.05i -Gblack -O -K >> example_14.ps

A psxy parancs magyarázata feljebb.

Az alapparancs utáni paraméterek jelentése:

paraméter
jelentés
mean.xyz Fájl, amiből az adatok olvasása történik
-Ss0.05i A -S szimbólumokat rajzol, itt az -Ss0.05i azt jelenti, hogy a szimbólum négyzet, amely átlója 0,05 inch.


surface mean.xyz -R -I1 -Gdata.grd

A surface parancs fájlból olvas számhármasokat (x,y,z) és bináris grid fájlt készít, amiben az értékek rácspontokra vannak átszámolva.

Az alapparancs utáni paraméterek jelentése:

paraméter
jelentés
mean.xyz xyz fájl
-I1 Rácsköz/távolság beállítása. Itt 1 fok.
-Gdata.grd A  a kimenő fájl nevét  adja meg. A kimeneti fájl egy 2_D.grd bináris fájl

grdcontour data.grd -J -B2f1WSne -C25 -A50 -Glct/cb -S4 -O -K -X-3.25i -Y-3.55i >> example_14.ps

A grdcontour parancs a kontúrok rajzolásában segít.

Az alapparancs utáni paraméterek jelentése:

paraméter
jelentés
data.grd A grd fájl, amiben az adatok vannak.
-C25 Kontúr, vonal felrajzolása. 25-ös léptékben rajzol egy kontúrvonalat.
-A50 Vonalak jelölése. 50-es léptékben a kontúrvonalat megvastagítja.
-G3i/10 A -G  segítségével a kontúrvonalak jelölőszámainak helyét lehet beállítani
-S4 A kontúrvonalakat lehet vele simítani

psxy -R -J mean.xyz -Ss0.05i -Gblack -O -K >> example_14.ps

Az alapparancs utáni paraméterek jelentései már fentebb tisztázva vannak



grdtrend data.grd -N10 -Ttrend.grd

A grdtrend parancs egy grd fájl adataira polinomiális trendet illeszt (legkisebb négyzetek módszere).

Az alapparancs utáni paraméterek jelentése:

paraméter
jelentés
data.grd A .grd fájl
-N10 A -N meghatározza a paraméterek számát, hogy a modell hány paraméterrel számol.
-Ttrend.grd A -T beleírja a trendet a megadott grd fájlba.

project -C0/0 -E7/7 -G0.1 > track

A project parancs metszetek készítéséhez nyújt segítséget.

Az alapparancs utáni paraméterek jelentése:

paraméter
jelentés
-C0/0 A kezdőpont beállítása. Itt (0,0).
-E7/7 A végpont beállítása. Itt (7,7)
-G0.1 A metszet mentén az adatok távolsága egymástól.


grdcontour trend.grd -J -B2f1wSne -C25 -A50 -Glct/cb -S4 -O -K -X3.25i >> example_14.ps

A grdcontour parancs a kontúrok rajzolásában segít.

Az alapparancs utáni paraméterek jelentése:

paraméter
jelentés
-Glct/cb A -G  segítségével a kontúrvonalak jelölőszámainak helyét lehet beállítani.


psxy -R -J track -Wthick,. -O -K >> example_14.ps

Az alapparancs utáni paraméterek jelentése:

paraméter
jelentés
track Fájlnév
-Wthick,. A felrajzolt szimbólum méretét lehet beállítani.


grdtrack track -Gdata.grd | cut -f3,4 >! data.d
grdtrack track -Gtrend.grd | cut -f3,4 >! trend.d

A grdtrack parancssal az x,y mezőből egy vonal mentén lehet mintát venni.

Az alapparancs utáni paraméterek jelentése:


paraméter
jelentés
track Az x,y fájl neve, ez adja meg az x,y koordinátákat a mintához.
-Gdata.grd
-Gtrend.grd
A -G  a grd fájl neve, ami az adatokat tartalmazza az x,y pontokban
cut -f3,4 A bemenet (standard input) vagy a paraméterben megadott fájlok minden sorának egy adott részét vágja ki. A '-f3,4' mező kivágása, itt 3.,4. mező - sorrend nem számít


psxy `minmax data.d trend.d -I0.5/25` -JX6.3i/1.4i data.d -Wthick -O -K -X-3.25i -Y-1.9i \ -B1/50WSne >> example_14.ps

A minmax parancs kikeresi a megadott fájlok adatainak szélsőértékeit, a minimumot és a maximumot.

Az alapparancs utáni paraméterek jelentése:


paraméter
jelentés
-I0.5/25 Jelentést ad a min/max-ról a megadott tömbökben


psxy -R -J trend.d -Wthinner,- -O -U"Example 14 in Cookbook" >> example_14.ps

Az alapparancs utáni paraméterek jelentése:

paraméter
jelentés
-U"Example 14 in Cookbook" A -U  kapcsoló kirak egy GMT logót az aktuális dátummal és idővel.


Tovább