GEBCO_2014 Grid (30 arc-second interval) with GMT-5

Table of Content

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

以下のスクリプトで鳥瞰図を描きます.
[cc]
#!/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
[/cc]
以下のような図ができます.

コメントを残す

メールアドレスが公開されることはありません。 * が付いている欄は必須項目です

このサイトはスパムを低減するために Akismet を使っています。コメントデータの処理方法の詳細はこちらをご覧ください