(37) Spectral coherence between gravity and bathymetry grids¶
The next script produces the plot in Figure. We demonstrate how grdfft is used to compute the spectral coherence between two data sets, here multibeam bathymetry and satellite-derived gravity. The grids are detrended and tapered before the Fourier transform is computed; the intermediate plots show the grids being extended and padded to a suitable dimension.
#!/usr/bin/env bash
# GMT EXAMPLE 37
#
# Purpose: Illustrate 2-D FFT and coherence between gravity and bathymetry grids
# GMT modules: basemap, plot, makecpt, grdfft, grdimage, grdinfo
# Unix progs: rm
#
gmt begin ex37
# Testing gmt grdfft coherence calculation with Karen Marks example data
# Prefix of two .nc files
G=grav.V18.par.surf.1km.sq
T=mb.par.surf.1km.sq
gmt set GMT_FFT kiss
gmt grdinfo @$T.nc -Ib > bbox
scl=1.4e-5
sclkm=1.4e-2
gmt makecpt -Crainbow -T-5000/-3000
gmt grdimage @${T}.nc -I+a0+nt1 -R@${T}.nc -Jx${scl}i -C -X1.474i -Y1i
gmt basemap -R-84/75/-78/81 -Jx${sclkm}i -Ba -BWSne+t"Multibeam bathymetry"
gmt makecpt -Crainbow -T-50/25
gmt grdimage @${G}.nc -I+a0+nt1 -R@${G}.nc -Jx${scl}i -C -X3.25i
gmt basemap -R-84/75/-78/81 -Jx${sclkm}i -Ba -BWSne+t"Satellite gravity"
gmt grdfft @${T}.nc @${G}.nc -E+wk -N192/192+d+wtmp > cross.txt
gmt makecpt -Crainbow -T-1500/1500
gmt grdimage ${T}_tmp.nc -I+a0+nt1 -R${T}_tmp.nc -Jx${scl}i -X-3.474i -Y3i
gmt plot bbox -L -W0.5p,-
gmt basemap -R-100/91/-94/97 -Jx${sclkm}i -Ba -BWSne+t"Detrended and extended"
gmt makecpt -Crainbow -T-40/40
gmt grdimage ${G}_tmp.nc -I+a0+nt1 -R${G}_tmp.nc -Jx${scl}i -C -X3.25i
gmt plot bbox -L -W0.5p,-
gmt basemap -R-100/91/-94/97 -Jx${sclkm}i -Ba -BWSne+t"Detrended and extended"
gmt set FONT_TITLE 24p
gmt plot -R2/160/0/1 -JX-6il/2.5i -Bxa2f3g3+u" km" -Byafg0.5+l"Coherency@+2@+" \
-BWsNe+t"Coherency between gravity and bathymetry" -X-3.25i -Y3.3i cross.txt -i0,15 -W0.5p
gmt plot cross.txt -i0,15,16 -Sc0.075i -Gred -W0.25p -Ey
rm -f cross.txt *_tmp.nc ?.cpt bbox
gmt end show