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

3.4. Fájlban tárolt adatok ábrázolása térképen (jelöléssel) II.

A script és a hozzá tartozó fájlok itt (34.tar - UNIX tar formátumban, illetve alternatív módo 34.zip - WinZip tömörítéssel) letölthetőek.

 - Vetület és térképméret beállítása
 
 psbasemap "--PAPER_MEDIA=Custom_640x455" "--ANNOT_FONT_SIZE=10p" -R15.5/45.5/23/49r -X50p -Y50p \ -JX18c/12c -Ba1/a1 -P -V -K  > tmp.ps

A psbasemap mögötti paraméterek már korábban tisztázva lettek.
 
- Interpoláció

A blockmedian és surface parancsok már korábban tisztázva.

blockmedian -R -I5m  -V  synop.txt > syn_grd.txt
surface syn_grd.txt -R -I5m -V  -Gsyn.grd

A színskála elemeinek beállítása.
 

 c_min="-10."
 c_max="20."
 c_step="1."
 makecpt -Crainbow -T${c_min}/${c_max}/${c_step} > haho.cpt

- A szinezett mező megrajzolása

 grdview syn.grd -R -JX -Chaho.cpt -Qs  -K -O -V  >> tmp.ps

- Skála felrajzolása

 psscale "--ANNOT_FONT_SIZE=10p" -D19c/6c/12c/0.4c -B${c_step}g${c_step} -Chaho.cpt -E -K -O \
>> tmp.ps


- Térképkontúrok, folyók stb. megrajzolása
 
 pscoast -R -JX  -I1/1p/0/0/0  -I2/1p/0/0/0  -N1p/0/0/0  -W1p/0/0/0 -Dh -S46/153/255 -V -K -O \
>> tmp.ps

 
- A megfigyelések felrajzolása a térképre, itt szükség van kis programozási ismeretekre

exec < synop.txt
while read acl
do
  lat=`echo $acl | awk '{print $2}'`
  lon=`echo $acl | awk '{print $1}'`
  tval=`echo $acl | awk '{print $3}'`

  echo $lon $lat  | psxy  -R -JX -P -St0.2c -G40/0/255 -K -O >> tmp.ps
  echo $lon $lat 10 0 0 1 $tval | pstext -R  -JX -G0/0/255 -P -K -O >> tmp.ps        
done 

 echo $lon $lat  10 0 0 1 $tval | pstext -R  -JS -G0/0/255 -P -O  >> tmp.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:

parancs
jelentés
-St0.2c A -S szimbólumokat rajzol, itt az -St0.2c azt jelenti, hogy a szimbólum háromszög, amel köré írható kör átmérője 0,2 cm.
-G40/0/255 A jelek, szimbólumok színének beállítása.

A pstext paranccsal szövegeket rajzolhatunk fel a térképre.

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

parancs
jelentés
-G0/0/255 A szöveg szinének beállítása, jelen esetben kék.
.
- Képpé alakítás

 convert +antialias tmp.ps tmp.png

A kapott ábra:




A scriptet készítette: Kertész Sándor


Tovább