Membuat Peta 3D Dengan GMT

Pada postingan kali ini, saya akan berbagi cara membuat peta 3 dimensi dengan GMT. Selain itu saya tampilkan juga cara mengatur parameter-parameter GMT sepeti mengatur ukuran kertas, ukuran font dan jenis font. Serta saya berikan juga cara melakukan resampling data grid. 

Berikut adalah scriptnya :

set F=mentawai.ps
set R=98/100/-2/0

#Setting parameter untuk GMT
gmtset COLOR_NAN 255/255/255 DOTS_PR_INCH 300 ANNOT_FONT_PRIMARY Helvetica-Narrow-Bold ANNOT_FONT_SIZE_PRIMARY 15 ANNOT_FONT_SIZE_SECONDARY 14  HEADER_FONT Helvetica-Narrow-Bold LABEL_FONT Helvetica-Narrow-Bold LABEL_FONT_SIZE 25 HEADER_FONT_SIZE 16 PAPER_MEDIA a4 Y_AXIS_TYPE ver_text PLOT_DEGREE_FORMAT DF

#File warna untuk kontur
echo -6000 255 255 255 0 255 255 255> zero.cpt
echo 0 58 58 58 1000 58 58 58>> zero.cpt

#Membuat file warna untuk elevasi
makecpt -Cglobe -Z > elev.cpt

#Membuat kontur dengan interval 500 meter
grdcontour indo.nc -Jm1.7i -E230/30 -R%R% -C250 -A500 -K > %F%

#Menampilkan batas pantai
pscoast -J -E230/30 -R%R% -B1/1neSW -Gblack -O -K -Df -T97.5/-0.4/1.4i --HEADER_FONT_SIZE=12p>> %F%

#Sampling file topografi menjadi 2 menit disimpan dengan nama menit.nc
grdsample indo.nc -I2m -Gmenit.nc

#Menampilkan peta 3D hasil sampling
grdview menit.nc -Jm1.7i -Jz0.0003i -Czero.cpt -E230/30 -R%R%/-6000/1000 -N-6000/lightgray -Qsm -O -K -B1/1/2000:"Topo (m)":neswZ -Y2i >> %F%

#Membuat Judul Peta
echo 5 4.5 40 0 0 BC S@#iberut@# I@#sland | pstext -R0/10/0/10 -Jx1i -K -O -U/-0.2/-7.1/"Membuat Peta 3D">> %F%

#Menampilkan peta Topografi dan konturnya
grdview indo.nc -Jm1.7i -R%R% -Celev.cpt -Qi  -E230/30 -K -O -Y-2i -X13.8 -N-6000/lightgray -B1/1neSW >> %F%
grdcontour indo.nc -Jm1.7i -E230/30 -R%R% -C500 -A500 -K -O >> %F%
pscoast -J -R -K -O -Dh -Lf98.5/-2.5/17/150+l+jr -W -E230/30>> %F%

#Menampilkan peta 3D topografi
grdview indo.nc -Jm1.7i -R%R%/-6000/1000 -Celev.cpt -Jz -Qi  -E230/30 -O -Y2.i -N-6000/lightgray -B1/1/2000:"Topo (m)":neswZ >> %F%

#Mencetak peta dengan nama mentawai.jpg
ps2raster mentawai.ps -Gc:\programs\gs\gs8.53\bin\gswin32c -Tj -P -Fmentawai

Simpan dengan ekstensi ".bat". Tanda '#' menunjukkan bahwa baris tidak akan dieksekusi oleh GMT.

Perintah grdsample digunakan untuk melakukan resampling data grid menjadi 2 menit busur. Sedangkan untuk mengatur parameter-parameter GMT bisa diatur di gmtset.

Selamat mencoba dan semoga bermanfaat. Salam Orang Indonesia...