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!

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