Kamis, 12 Juli 2012

GMT(Generic Mapping Tools) untuk pemula

GMT (Generic Mapping tools) adalah software open source (alias gratis) khusus menangani masalah pemetaan. Mula-mulai GMT hanya bisa diinstal di linux, akan tetapi lambat laut dikembangkan oleh Paul Wessel dan Walter Smith sehingga bisa diinstal di windows.

GMT memiliki beberapa keuntungan dan kekurangan, misalnya :
(+) dapat membuat peta dalam berbagai format
(- ) tidak adanya GUI sehingga harus lihai dalam menguasai scrips, otomatis akan semakin rumit scripnya jika ingin membuat peta yang kompleks.

Sebenarnya, inti dari GMT hanya 3 hal, yaitu :
1) notepad
2) GSView 49
3) Saya menyebutnya sbg progam run

Steps-steps memulai GMT :
1) Progam-->rum--> cmd-->pilih direktori yang akan kita gunakan untuk menyimpan file kita, misalnya :

C: D:
D: cd "anggav"
D\anggav>

2) Membuat script di notepad lalu simpan dalam format .(dot) bat

3) Running script yang udah kita buat tadi (misal namanya peta.ps) di progam run (D\anggav>peta.ps )

4) Buka hasilnya di GSView 49

Oh ya, ini salah satu contoh scrip peta sederhana :

set R=90/105/-5/6 (ini maksudnya membuat region long 90-105 dan lat -5s/d6)
set O=peta.ps
(nama outputnya peta.ps)

gmtset ANNOT_FONT_SIZE 11
gmtset LABEL_FONT_SIZE 11

psbasemap -R%R% -JM15C -Bg2a1:.PETA_INDONESIA: --HEADER_FONT_SIZE=14 -P -K -Y17 > %O%



pscoast -R -JM -Bg2a1 -Ggreen -Slightblue -P -Dh -W -O -K >> %O%
psxy stasiun.dat -R -JM -St0.5c -Gyellow -P -O -K >> %O%
pstext stasiun.dat -R -JM -P -O -K>> %O%
psmeca moment2.dat -R -JM -P -O -Sm1c -T0 -C -M -Wthin -K >> %O%
gawk "{print $1,$2}" moment2.dat | psxy -R -JM -Sc0.3c -Gred -Wthin -P -O -K>> %O%
psxy trench.gmt -R -JM -P -O -Gblack -Sf1/0.3lt -K >> %O%
psxy transform.gmt -R -JM -P -O -Gblack -Sf2/0.4rs -K>> %O%

echo 94 2 10 0 1 LT A > point.dat
echo 97 4 10 0 1 LT A' >> point.dat
gawk "{print $1,$2}" point.dat | psxy -R -JM -WS/red -P -O -K >> %O%
pstext point.dat -R -JM -P -O -K >>%O%


gawk "{print $1,$2,$3}" moment2.dat | project -C94/2 -E97/4 -Qxyzpqrs W-100/100 > cross.dat

psbasemap -JX15/-10 -R0/400/0/50 -B100:"Distance (km)":/30NWes:"Depth (km)": -P -O -Y-15 -K >> %O%

gawk "{print $4,$3}" cross.dat | psxy -JX -R -Sc0.5c -Gred -P -O >> %O%

0 komentar:

Posting Komentar

next page