Membuat Kontur Dengan GMT

Pada postingan kali ini, saya akan berbagi cara membuat garis-garis kontur dengan GMT. Kali ini kita akan menggunakan perintah psbasemap sebagai dasar, grdimage untuk menampilkan peta topografi, grdcontour untuk membuat kontur serta ps2raster untuk menghasilkan gambar akhir. Sebelumnya akan saya jelaskan dulu alasan penggunaan perintah ps2raster. Jika kita membuat peta dengan data raster, kita tidak bisa langsung melakukan convert untuk menghasilkan file JPEG atau file image lainnya. Untuk menghasilkan file image kita harus menggunakan perintah ps2raster.

Sebagai tambahan untuk ps2raster, pastikan Ghostscript dan GSview anda instal di C:\programs bukan di C:\Program Files. Sebab pada Windows, Ghostscript tidak ditambahkan pada path. Sehingga nantinya kita harus menulis path ke Ghotscript tersebut. Dan untuk memudahkan, Ghostscript dan GSview installah di C:\programs. Jika anda terlanjur menginstal di C:\Program Files, uninstallah dan instal ulang di C:\programs.

Langkah selanjutnya, pastikan anda memiliki file topografi raster grid 2D. Anda dapat mendownloadnya disini. Kemudian buatlah file warna seperti berikut :

  -10000    153    0    255    -9500    153    0    255
-9500    153    0    255    -9000    153    0    255
-9000    153    0    255    -8500    153    0    255
-8500    136    17    255    -8000    136    17    255
-8000    119    34    255    -7500    119    34    255
-7500    102    51    255    -7000    102    51    255
-7000    85    68    255    -6500    85    68    255
-6500    68    85    255    -6000    68    85    255
-6000    51    102    255    -5500    51    102    255
-5500    34    119    255    -5000    34    119    255
-5000    17    136    255    -4500    17    136    255
-4500    0    153    255    -4000    0    153    255
-4000    27    164    255    -3500    27    164    255
-3500    54    175    255    -3000    54    175    255
-3000    81    186    255    -2500    81    186    255
-2500    108    197    255    -2000    108    197    255
-2000    134    208    255    -1500    134    208    255
-1500    161    219    255    -1000    161    219    255
-1000    188    230    255    -500    188    230    255
-500    215    241    255    -200    215    241    255
-200    241    252    255    0    241    252    255
  0      0       97      71      50      16      122     47
  50     16      122     47      500     232     215     125
  500    232     215     125     1200    161     67      0      
  1200   161     67      0       1700    130     30      30
  1700   130     30      30      2800    110     110     110
  2800   110     110     110     4000    255     255     255
  4000   255     255     255     6000    255     255     255

Simpan dengan nama "elev.cpt". Kemudian buatlah baris-baris perintah berikut :

psbasemap -JM18c -R100/105/-5/0 -K -Ba2g2WSne:: > map.ps
grdimage indo.nc -R -JM -Celev.cpt -K -O -Ba2g2wsne:: >> map.ps
grdcontour indo.nc -R -JM -O -C300 -A600+f14 >> map.ps
ps2raster map.ps -Gc:\programs\gs\gs8.53\bin\gswin32c -Tj -P -Fpeta

Simpan dengan ekstensi ".bat". Setelah dieksekusi, hasilnya akan seperti berikut :

Anda dapat mengatur interval garis kontur dan interval labelnya di grdcontour. -C300 artinya interval garis konturnya setiap 300 meter. -A600+f14 artinya interval labelnya setiap 600 meter dengan ukuran font 14. Sedangkan pada ps2raster, -Tj dan -Fpeta artinya hasil akhir berupa file JPEG dengan nama "peta.JPEG". Yang harus anda perhatikan adalah perintah terakhir sebelum ps2raster tidak boleh ada -K.

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