Menggambar Kontur Batimetri

Kontur Batimetri Hitam-Putih
Kontur Batimetri Hitam-Putih

File grid srtm_ntt.grd bisa kita interpolasi untuk menggambarkan kontur batimetri. SRTM 15 plus merupakan gabungan data batimetri dan topografi dari berbagai sumber akuisisi data, baik data dari satelit maupun data survei kelautan.

Generic Mapping Tools akan membantu kita melakukan plotting kontur batimetri  dari kedalaman -7500m sampai -500m dengan kontur interval 500m. Anotasi diberikan setiap 2000m. Daratan digambarkan memakai warna abu-abu.  Script singkat untuk melakukan plotting peta ada di bawah ini:

#!/usr/bin/bash

# Plot Basemap Dan Kontur Batimetri
# GMT 5.2.1
# oleh: Hendro Nugroho

gmt set MAP_FRAME_TYPE plain

# Data grid topografi: SRTM15 plus
ntt='srtm_ntt.grd'
out='kontur_batimetri_ntt.ps'
area='-R118.5/127.5/-11.5/-7.5'
k_int=500	# interval kontur
a_int=2000	# interval anotasi kontur
k_limit=-7500/-500

# frame peta
gmt psbasemap $area -JM8.5i -Xc -Bxa2f1 -Bya1f1 -K > $out

# peta garis pantai
gmt pscoast -R -JM -Xc -B -Di -Ggray -Wthinnest -O -K >> $out

# plot kontur batimetri
gmt grdcontour $ntt -R -JM -Xc -C$k_int -A$a_int -L$k_limit -O >> $out

Peta hasil pengeplotan bisa dilihat pada featured image di atas.

Dengan sedikit tambahan warna, hasil akhir peta akan lebih menarik. Level warna kontur bisa kita berikan dengan terlebih dahulu membuat tabel warna. Lalu kita plot file grid dengan perintah grdimage.

Mari kita sisipkan kode berikut ini di bawah deklarasi variabel terakhir.

cpt='kontur.cpt'
#membuat tabel warna level kontur
gmt makecpt -Chaxby -T-8000/0/1000 -Z > $cpt

# plot grd dengan tabel warna baru
# -Yc ditambahkan untuk memberi tempat skala warna
# yang akan diplot di bawah peta.
gmt grdimage $ntt -JM8.5i -Xc -Yc $area -C$cpt -K > $out

Perubahan kode script serta tambahan perintah untuk menggambar skala warna di bawah peta ada di bawah ini:

# ==> tambahan kode -Yc
# frame peta
gmt psbasemap $area -JM8.5i -Xc -Yc -Bxa2f1 -Bya1f1 -O -K >> $out

# peta garis pantai
gmt pscoast -R -JM -Xc -Yc -B -Di -Ggray -Wthinnest -O -K >> $out

# plot kontur batimetri
gmt grdcontour $ntt -R -JM -Xc -Yc -C$k_int -A$a_int -L$k_limit -O -K >> $out

# plot skala warna
gmt psscale -D5/-1/10/0.5h -Xc -C$cpt2 -Bf1000a2000 -O >> $out

Inilah hasil akhir peta dengan kontur batimetri berwarna:

peta_kontur_batimetri_berwarna

Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s