GEBCO-2014 30 arc-second(約500m)の海底地形グリッドデータをダウンロードし,GMT-5を用いて図化や解析処理を行っていきます.

データのダウンロード

GEBCO_2014 Gridデータはこちらからダウンロードできます.ダウンロード範囲を指定するボックスがありますが,デフォルトでは全範囲となっているようです.GEBCO_2014 Grid (30 arc-second interval)をクリックし,扱いやすそうな2D netCDFを選択し,Add data to basketをクリックします.
View basketをクリックし,データのリクエストを確認します.ここで,まだユーザー登録をしていない場合は登録を済ませる必要があります.
確認できたら,Checkout your requestをクリックし,Licensingで「I accept the above terms and conditions.」をクリックし,Submit licenseをクリックします(毎回必要になるのかわかりません).
My accountのRequest statusをクリックし,Ready to downloadとなっていたら,Downloadできます.zipのファイルサイズは1,124MBでした.適当なディレクトリに解凍しておきます.
データフォーマットの情報はこちらにありますが,GMTでは直接読み込むことができます.

鳥瞰図を描く Create a Bird’s-eye view

以下のスクリプトで鳥瞰図を描きます.

#!/bin/bash
# bash script for viewing GEBCO_2014 30 arc-second grid using GMT-5.2.1 or later
#   Windows Subsystem for Linux (WSL)の利用可
#   WSLではGMTはWindows版とLinux版の両方が利用できる
#     $gmtをgmt.exe (Windows)またはgmt (Linux)と定義し,$gmtを使う
#     WindowsのCドライブの参照の方法を定義する
#       Windows: c:  Linux: /mnt/c
#   Windows 10のJupyter Notebook (Python)からsubprocessで実行可能
# Author Jun Sasaki @UTokyo
# Coded on Nov. 25, 2017, Revised on Nov. 25, 2017
# USAGE: $ bash gmt_gebco2014_30.sh [windows] [ngg]
#   最初の引数がwindowsまたはlinuxの場合は特別扱い
#   ngg: Using exsisting grdgradient file
#   Input: GEBCO_2014_2D.nc (original GEBCO_2014 30 arc-second grid data)
#   Output: GEBCO_2014_2D.eps GEBCO_2014_2D_png
##################################
gmt_system=linux  # set defalut: =windows or =linux
### check the 1st argument
if [ "$1" = "windows" ] ; then
  gmt_system=windows; shift
elif [ "$1" = "linux" ] ; then
  gmt_system=linux; shift
fi
if [ "${gmt_system}" = "windows" ] ; then
  dat_dir="c:" # Set dir containing data
  gmt="gmt.exe"
  echo -n Windows GMT-; $(echo $gmt --version)
else
  dat_dir="/mnt/c" # Set dir containing data
  gmt="gmt"
  echo -n Linux GMT-; $(echo $gmt --version)
fi
### Set parameters
rxy=125/150/25/50/-20000/5000  ### region to be drawn
proj=X15c  # x-y map projection, X: linear, 15c: 15 cm width
vmin=-10000; vmax=10000; inc=500
palette=globe
cpt=my_globe.cpt  # color palette
### Input data dir and file name
fdata=GEBCO_2014_2D  # GEBCO_2014_2D.nc without extension
indir=../GEBCO/
outdir=./output/
ps=${outdir}${fdata}.ps # output file name
fdata=${indir}${fdata}.nc # GEBCO_2014_2D.nc
gg=grad30.grd
### LongitudeとLatitudeのaxis labelは+lで設定しても現れない
# http://thread.gmane.org/gmane.comp.gis.gmt.user/16895
Bf=wESnZ
Bx=a5f2.5
By=a5f2.5
Bz=a5000f1000+lElevation
Bs=a5000f2500g500+lElevation
### Set parameters of GMT for renewing gmt.conf
$gmt gmtset MAP_FRAME_TYPE          = plain
$gmt gmtset FONT_ANNOT_PRIMARY      = 14p,Helvetica,black
$gmt gmtset FORMAT_GEO_MAP          = ddd:mm:ssF  # F means
### Make a color pallete
$gmt makecpt -C${palette} -D -T${vmin}/${vmax}/${inc} >$cpt
if [ "$1" = "ngg" ] ; then
  echo "Using existing grdgradient file"
else
  $gmt grdgradient ${fdata} -G${gg}  -A45 -Ne0.4
fi
$gmt grdview $fdata -R${rxy} -JZ8c -Y6c -JM15c -C${cpt} -I${gg} -p170/40 -B"${Bf}" -Bx"${Bx}" -By"${By}" -Bz"${Bz}" -N-20000+glightgray -Qi300 -V -P -K > $ps
$gmt psscale -DjBC+w10c/0.3c+h+o0/-2c -R${rxy} -JM15c -JZ8c -C${cpt} -Bx"${Bs}" -By+lm -V -P -O >> $ps
$gmt psconvert ${ps} -A -Tg -E300 # PNG
$gmt psconvert ${ps} -A -Te       # EPS
rm $ps

以下のような図ができます.