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

Plotting Basemap Dengan CPT Asimetri

Basemap Nusa Tenggara Timur
Basemap Nusa Tenggara Timur

Kita akan melihat kembali script GMT pada posting sebelumnya (plotting peta jaringan YS) dan akan menggunakan script tersebut untuk membuat basemap dengan color palette table asimetri. Warna batimetri akan kita ambil dari abyss.cpt sedangkan warna topografi akan kita ambilkan dari arctic.cpt dan memofidikasi file ini dengan menghilangkan warna-warna untuk kedalaman laut.

Sebelum membuat tabel warna (cpt) mari kita lihat informasi grid file yang akan kita pakai (srtm_ntt.grd) yang merupakan hasil pemotongan dari grid global srtm 15 plus. Untuk keperluan tersebut kita memanggil perintah gmt grdinfo

user@kkomputer ~/GMT_Project $ gmt grdinfo srtm_ntt.grd
srtm_ntt.grd: Title: Produced by grdcut
srtm_ntt.grd: Command: grdcut srtm15-plus/topo15.grd -R118.5/127.5/-11.5/-7.5 -Gsrtm_ntt.grd
srtm_ntt.grd: Remark: Derived from the full resolution shorelinen
srtm_ntt.grd: Gridline node registration used [Geographic grid]
srtm_ntt.grd: Grid file format: nf = GMT netCDF format (32-bit float), COARDS, CF-1.5
srtm_ntt.grd: x_min: 118.5 x_max: 127.5 x_inc: 0.00416666666667 name: longitude [degrees_east] nx: 2161
srtm_ntt.grd: y_min: -11.5 y_max: -7.5 y_inc: 0.00416666666667 name: latitude [degrees_north] ny: 961
srtm_ntt.grd: z_min: -7079.01904297 z_max: 2847 name: z
srtm_ntt.grd: scale_factor: 1 add_offset: 0
srtm_ntt.grd: format: netCDF-4 chunk_size: 136,138 shuffle: on deflation_level: 3

Sekarang kita tahu nilai z minimal pada file srtm_ntt.grd adalah -7079 meter dan nilai z maksimumnya adalah 2847 meter. Menggunakan keterangan ini kita bisa membuat color palette table asimetri yang kita inginkan secara terpisah dengan memberikan perintah berikut ini:

$ gmt makecpt -Cabyss -T-7500/0/250 -Z -N > mynew.cpt
$ gmt makecpt -Cmod_arctic.cpt -T0/3000/500 -Z >> mynew.cpt

Atau kita bisa memasukkannya menjadi bagian dari keseluruhan script plotting basemap berikut ini:

#!/usr/bin/bash

# Plot Basemap Menggunakan CPT asimetri
# GMT 5.2.1
# oleh: Hendro Nugroho

gmt set MAP_FRAME_TYPE plain

# Data grid topografi: SRTM15 plus
ntt='srtm_ntt.grd'
ntti='srtm_ntt_i.grd'
out='basemap_ntt.ps'
area='-R118.5/127.5/-11.5/-7.5'
# abyss.cpt dam modofikasi arctic.cpt [download dari cpt-city]
cpt='asimetri.cpt'

gmt makecpt -Cabyss -T-7500/0/250 -Z -N > $cpt
gmt makecpt -Cmod_arctic.cpt -T0/3000/500 -Z >> $cpt

# plot grd plus iluminasi
gmt grdimage $ntt -JM8.5i -Xc $area -I$ntti -C$cpt -K > $out

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

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

Peta hasil plotting menggunakan GMT ini bisa kita lihat dengan gs

$ gs basemap_ntt.ps

Sekarang Anda bisa mencoba sendiri membuat peta di atas menggunakan tabel warna yang sama atau menggunakan warna kesukaan Anda.
Mungkin Anda akan memadukan abyss.cpt dengan dem2.cpt? Silakan mencoba!

Plotting Peta Jaringan YS

stasiun-YS

Jaringan YS adalah jaringan seismic recorder sementara dalam rangka penelitian di daerah Busur Banda. Data titik koordinat dan nama stasiun ini akan saya ambil untuk memaparkan pengeplotan titik, jaringan tringulasi, dan plot nama stasiun menggunakan aplikasi Generic Mapping Tools. Background topografi yang akan dipakai kali ini diambil dari SRTM 15 plus dari NOAA. Master color palette table yang akan saya pakai adalah gray.cpt.

Berikut ini adalah script untuk plotting peta jaringan YS tersebut:

#!/usr/bin/bash

# Plot stasiun seismometer network YS
# Sumber: http://www.fdsn.org/networks/detail/YS_2014/
# Ekstraksi informasi: lon, lat, nama_stasiun
# Data grid topografi: SRTM15 plus
# GMT 5.2.1
# oleh: Hendro Nugroho

data='/home/eagle1/Data/srtm15-plus/'
srtm15_ori='topo15.grd'
ntt='srtm_ntt.grd'
ntti='srtm_ntt_i.grd'
out='net_YS.ps'
area='-R118.5/127.5/-11.5/-7.5'
cpt='ntt_gray.cpt'

# potong srtm15+
gmt grdcut $data$srtm15_ori $area -G$ntt

# membuat cpt abu-abu
gmt grd2cpt $ntt -Cgray > $cpt

# Membuat gradient dari file grid NTT
gmt grdgradient $ntt -A345 -Ne0.6 -G$ntti

# plot grd plus iluminasi
gmt grdimage $ntt -JM8.5i -Xc $area -I$ntti -C$cpt -K > $out 

# buat jaringan triangulasi
# Sumber: http://ds.iris.edu/mda/YS?timewindow=2014-2017
# Ekstraksi informasi:
# Lon Lat elev station
# 124.411500  -8.219400 69  ALRB
# Nama file: stasiun_YS.gmt; hasil triangulasi: ys.net.gmt
gmt triangulate stasiun_YS.gmt -M > ys.net.gmt

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

# siapkan clipper
gmt pscoast -R -JM -Xc -Di -Gc -O -K >> $out

# plot grd plus iluminasi
gmt grdimage $ntt -JM -Xc $area -I$ntti -C$cpt -O -K >> $out 

# undo clipping
gmt pscoast -R -JM -B -Xc -O -K -Q >> $out
#
# peta garis pantai
gmt pscoast -R -JM -Xc -B -Di -Wthinnest -O -K >> $out
# plot triangulasi
gmt psxy ys.net.gmt -R -JM -Xc -B -W -O -K >> $out
# plot stasiun
gmt psxy stasiun_YS.gmt -R -JM -Xc -Sc0.06i -Gred -O -K >> $out
gmt psxy stasiun_YS.gmt -R -JM -Xc -Sc0.04i -Gwhite -O -K >> $out
# plot nama stasiun dari file nama_stasiun.gmt
# Isi file adalah ekstraksi informasi yg terdiri dari
# Lon Lat Justify Station
# CB = Center Bottom
# 124.411500 -8.219400 CB ALRB
gmt pstext nama_stasiun.gmt -R -JM -F+f6p,Helvetica,red+j -D0/0.05i -N -Gwhite -Wthinnest -C0.02i -O >> $out

Silakan Anda coba!