Saturday 11 May 2013

Membuat Peta Topografi Dengan GMT

Pembaca yang budiman..:D
Mungkin anda bosan melihat tampilan peta yang cuman terdiri dari dua warna, yaitu warna daratan dan lautan. Kali ini akan saya bagi cara membuat peta topografi, sehingga tampilan peta anda menjadi cantik dan lebih enak dipandang mata.. :D

Sebelumnya download dulu data topografi dan bathymetri disini. Data tersebut didapat dari GEBCO (General Bathymetric Chart of the Ocean). Untuk mendownload data dari GEBCO anda tinggal mendaftarkan alamat email anda. Selain itu anda dapat pula menggunakan data topografi dan bathymetri dari sumber lain asalkan dalam format netCDF.


Langkah berikutnya, buatlah file color palette/file warna sebagai 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

lalu simpan dengan nama "elev.cpt". Kemudian tulislah baris-baris perintah berikut :

psbasemap -JM25c -R95/141/-13/9 -K -X2 -Y5.5 -Ba10g5WSne:: --HEADER_OFFSET=0.01i> map.ps
grdimage indo.nc -R -JM -Celev.cpt -K -O -Ba5g5wsne:: >> map.ps
echo 95.2 -7 > legenda
echo 101 -7 >> legenda
echo 101 -12 >> legenda
echo 95.2 -12 >> legenda
echo 95.2 -7 >> legenda
psxy -R -JM -W5.0 -M -O -K legenda -Gwhite >> map.ps
psscale -Celev.cpt -D0.3/0.1/2.5/0.4 -B4000:meter: -O -K --LABEL_FONT_SIZE=13p --LABEL_OFFSET=0.1p -Y1.85 >> map.ps

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


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


4 comments:

  1. bang Yos,
    habis data bathimetry ny d donlod,
    diapain bang?

    ReplyDelete
    Replies
    1. habis di donlod ya diplot om.. :D
      taruh aja satu folder ama script ".bat" nya

      Delete
  2. mas, abis data topografi didownload diapain ya?

    ReplyDelete
  3. bang cara input DEM di GMT gimna?

    ReplyDelete