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!

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