Membuat Peta Anomali Gravitasi dan Magnetik Dengan GMT

Pada kesempatan kali ini, saya akan berbagi cara membuat peta anomali gravitasi dan magnetik. Sebelumnya downloadlah data gravitasi dan magnet dahulu. Ada banyak model data gravitasi dan magnetik yang bisa anda download. Contohnya seperti di TOPEX, EarthByte dan NOAA untuk data gravitasi, dan GEOMAG untuk data magnet. Pada contoh kali ini saya akan menggunakan data gravitasi dari Earthbyte dan data model magnet EMAG3 dari GEOMAG.

Sebelum anda buat petanya, pertama harus anda pahami format data yang akan anda download. Pada data gravitasi dari EarthByte, format datanya adalah netCDF dengan satuan mGal. Data tersebut dikali 10 oleh EarthByte, sehingga kita harus membaginya dengan 10 kembali agar didapat nilai aslinya. Untuk membagi data tersebut kita gunakan perintah "grdmath".

Sedangkan pada data model magnetik EMAG3 dari GEOMAG format datanya adalah ASCII. Sehingga anda harus melakukan gridding untuk menampilkannya di peta. Sebelum melakukan gridding, data tersebut harus difilter dengan perintah "blockmean", "blockmedian", atau "blockmode". Setelah itu baru dilakukan gridding dengan perintah "surface".

Setelah anda download datanya, buatlah script berikut untuk peta anomali gravitasi :

set F=gravity.ps
set R=94/108/-6/6

#Membagi grav_down.15.2.nc dengan 10 agar didapat nilai sebenarnya
grdmath grav_down.15.2.nc 10 DIV = grav_down.nc

#Membuat file warna
grd2cpt grav_down.nc -Crainbow -Z > grav.cpt

#Plotting data
grdview grav_down.nc -JM18 -R%R% -Cgrav.cpt -Qi -K -Ba3f3g3 > %F%
pscoast -JM -R -B -O -W4 -K -Dh >> %F%
echo 5 6.7 25 0 0 BC GRAVITY ANOMALY MAP | pstext -R0/10/0/10 -Jx1i -K -O -U/0/-1.7/"Membuat Peta Anomali Gravitasi">> %F%
psscale -Cgrav.cpt -D0.7/0.1/12/0.4 -B50:mGal: -O -X20 -Y7 >> %F%

#Print peta
ps2raster %F% -Gc:\programs\gs\gs8.53\bin\gswin32c -Tj -P -Fgravity

Simpan dengan ekstensi ".bat".

Setelah dieksekusi, hasilnya akan seperti berikut :

Sedangkan untuk peta anomali magnetik buatlah script berikut :

set F=magnetic.ps
set R=94/108/-6/6

#Membuat file warna
echo -310 64 0 128 -200 0 0 117 > mag.cpt
echo -200 0 0 117 -150 0 0 157 >> mag.cpt
echo -150 0 0 157 -100 0 0 221 >> mag.cpt
echo -100 0 0 221 -50 202 202 255 >> mag.cpt
echo -50 202 202 255 0 255 255 0 >> mag.cpt
echo 0 255 255 0 50 255 140 85 >> mag.cpt
echo 50 255 140 85 100 230 0 230 >> mag.cpt
echo 100 230 0 230 150 187 0 187 >> mag.cpt
echo 150 187 0 187 303 136 0 136 >> mag.cpt

#Filtering data

blockmean WDMAM_NGDC_V1.1.xyz -R%R% -I1.5m > mag.xyz

#Gridding data

surface mag.xyz -R%R% -I1.5m -GWDMAM_NGDC_V1.1.grd

#Plotting data
grdview WDMAM_NGDC_V1.1.grd -JM18 -R%R% -Cmag.cpt -Qi -K -Ba3f3g3 > %F%
pscoast -JM -R -B -O -W4 -K -Dh >> %F%
echo 5 6.7 25 0 0 BC MAGNETIC ANOMALY MAP | pstext -R0/10/0/10 -Jx1i -K -O -U/0/-1.7/"Membuat Peta Anomali Magnetik" >> %F%
psscale -Cmag.cpt -D0.7/0.1/12/0.4 -B50:nT: -O -X20 -Y7 >> %F%

#Print peta
ps2raster %F% -Gc:\programs\gs\gs8.53\bin\gswin32c -Tj -P -Fmagnetic

Simpan dengan ekstensi ".bat".

Setelah dieksekusi hasilnya akan seperti berikut :

Sekian untuk kali ini. Selamat mencoba dan semoga bermanfaat. Salam Orang Indonesia...