Menampilkan Relief Bathymetri Pada Peta

Setelah sekian lama vakum, kali ini saya akan menjelaskan cara untuk menampilkan relief bathymetri pada peta. Pada contoh ini, saya menggabungkan antara data usia kerak samudera global dengan relief bathymetri. Perhatikan pada peta dibawah. Wilayah di sebelah kiri bujur 80 derajat hanya menampilkan data usia kerak samudera, sedangkan data di sebelah kanan bujur 80 derajat menampilkan data usia kerak samudera yang digabung dengan relief bathymetri. Data usia kerak samudera global didapat dari EarthByte dan data bathymetri dari GEBCO. Saya menggunakan GMT versi 5.

 

Untuk menggabungkan kontur usia lempeng samudera dan relief bathymetri, hal yang perlu diperhatikan adalah  node registration dan resolusi data grid. Node registration dan resolusi data harus sama agar bisa digabungkan. Ada dua jenis node registration pada data grid, yaitu gridline node registration dan pixel node registration.

1. Gridline registration
Node terletak pada perpotongan grid-line. Data-point nya adalah rata-rata pada dimensi cell yang terletak di tengah tiap-tiap node.

2. Pixel registration
Node terletak di tengah grid cell. Data-point nya adalah nilai rata-rata dari masing-masing cell.

Gridline- (kiri) dan pixel-registration (kanan)

Untuk mengetahui jenis node-registration dan resolusi data, gunakan perintah grdinfo. Ketikkan:

grdinfo age.3.2.my.nc

kemudian

grdinfo indo.nc 

Pada contoh ini, data usia kerak samudera memiliki gridline registration sedangkan data bathymetri memiliki pixel-registration. Resolusi data usia kerak samudera adalah 0.0333333333 atau 2 menit, sedangkan data bathymetri resolusinya 0.00833333333 atau 30 detik. Untuk mengatasinya digunakan perintah grdsample. Untuk membuat file iluminasi atau pencahayaan untuk menampilkan relief bathymetri nya digunakan perintah grdgradient.

Pertama-tama buatlah file pencahayaannya. Ketikkan perintah berikut:

grdgradient indo.nc -A0 -Ne0.6 -Giluminasi.nc=nb/a

grdsample iluminasi.nc -T -I2m -Ggrid_iluminasi.nc

-T artinya merubah node-registration, dari pixel ke gridline dan sebaliknya

-I2m artinya merubah resolusi data menjadi 2 menit

Selanjutnya buatlah batch file dan ketik perintah berikut:

set F=gradasi.ps
grd2cpt age.3.2.my.nc -Crainbow -I -Fr -L0/280 -Z > age.cpt
psbasemap -JM20 -R50/110/-20/15 -K -Ba20g20WSne -Y5 > %F%
grdimage age.3.2.my.nc -R -JM -Cage.cpt -K -O -Ba20g20wsne >> %F%
grdimage age.3.2.my.nc -R -JM -Cage.cpt -K -O -Ba20g20wsne -Igrid_iluminasi.nc >> %F%
pscoast -JM -R -K -O -Dl -Glightgray -Wthinnest >> %F%
psxy -R -JM -Wthin,red -Sf0.1i/0.02i+l+t -Gblack -O -K trench.gmt >> %F%
psxy -R -JM -Wthin -O -K transform.gmt  >> %F%
psxy -R -JM -Wthin -O -K ridge.gmt  >> %F%
psxy -R -JM -Wthinner,-  -O -K GSFML_SF_FZ_KM.gmt >> %F%
psxy -R -JM -Wthinnest,-  -O -K GSFML_SF_ER_KM.gmt >> %F%
psscale -D10/-1/20/0.4h -Cage.cpt -B20:"Age of Oceanic Lithosphere [m.y.]": -O >> %F%
ps2raster %F% -TG -P -A

Simpan dan jalankan. Maka anda akan mendapatkan peta seperti gambar diatas.

Selamat mencoba dan semoga bermanfaat.