color shaded relief - 陰影図をつくる（光をあてて立体的にみせる）

---

  使うコマンド　　makecpt, grdgradient, grdimage, colorbar, coast   basemap



bash スクリプト例
#  plot bathymetry color fill map

# parameter setting
region=130/145/25/35                    # map region east/west/south/north
proj=M10                                # map projection and scale
ticks=a5f1g5                            # boundary tick info
frame=WSne+tcolor_shade                     # boundary frame info
climit=-8000/0/1000                  # color table min/max/interval
light=0/90                             # shade direction
grdfile=JTOPO1_30.grd                   # input bathymetry grid file
intfile=JTOPO1_30.int                   # intencity (shade) file
#
#
# plot
gmt begin colshade
     gmt makecpt -Chaxby -T$climit -Z
     gmt grdgradient $grdfile -A$light -G$intfile -Ne0.6
     gmt grdimage $grdfile -R$region -J$proj -I$intfile
     gmt colorbar -DJMR -B1000a2000
     gmt coast -Di -Ggray -Wthin,black
     gmt basemap -B$ticks -B$frame
gmt end show
#


Tips

• 地形の傾斜から特定の方向から光をあてた場合の陰影を計算し、grdimageで-I指定すると浮き彫り式の地形図が得られる。直感的に地 形を理解することができ、鳥瞰図と違って、死角がなく、距離や方向を計測しやすいのところが利点。モノクロでも比較的立体感を出せるので、用途は広い。
• ある方向の地形要素を強調したい（検出したい）場合は$lightで指定している方向を１方向とし、ーNオプションの平滑化パラメタを調整して特定走向の断層地形などを強調することができる
• 全体として見栄えのよい図をつくる場合は、このスクリプトのように２方向（90°ずれの必要はない）から光をあてるとよい。
• 卓越する方向性が２つありじっくり解析したい場合は、それぞれの方向の光をあてた同じ仕様の図を別々に2枚つくり、Illustratorで重ねて透明機能を使う方法もある