grdcut — GMT 6.5.0 documentation
docs.generic-mapping-tools.org
grdgradient — GMT 6.5.0 documentation
docs.generic-mapping-tools.org
grdimage — GMT 6.5.0 documentation
docs.generic-mapping-tools.org
grdcut・grdgradient・grdimage はDEMを使用する時によく使うモジュールです。
DEM3兄弟 とか呼んであげてください。
DEMとは
DEM（数値標高モデル） | 国土地理院
www.gsi.go.jp
GMTではオンラインでDEMを取得することができます。
取得できるものは以下の通り。
Remote Data Sets — GMT 6.5.0 documentation
docs.generic-mapping-tools.org
| Code | Dimensions | Reg | Size | Description |
|---|---|---|---|---|
| 01d | 360 x 180 | g,p | 128 KB | 1 arc degree global relief (SRTM15+V2.1 @ 111 km) |
| 30m | 720 x 360 | g,p | 435 KB | 30 arc minute global relief (SRTM15+V2.1 @ 55 km) |
| 20m | 1080 x 540 | g,p | 918 KB | 20 arc minute global relief (SRTM15+V2.1 @ 37 km) |
| 15m | 1440 x 720 | g,p | 1.6 MB | 15 arc minute global relief (SRTM15+V2.1 @ 28 km) |
| 10m | 2160 x 1080 | g,p | 3.4 MB | 10 arc minute global relief (SRTM15+V2.1 @ 18 km) |
| 06m | 3600 x 1800 | g,p | 8.8 MB | 6 arc minute global relief (SRTM15+V2.1 @ 10 km) |
| 05m | 4320 x 2160 | g,p | 13 MB | 5 arc minute global relief (SRTM15+V2.1 @ 9 km) |
| 04m | 5400 x 2700 | g,p | 19 MB | 4 arc minute global relief (SRTM15+V2.1 @ 7.5 km) |
| 03m | 7200 x 3600 | g,p | 33 MB | 3 arc minute global relief (SRTM15+V2.1 @ 5.6 km) |
| 02m | 10800 x 5400 | g,p | 71 MB | 2 arc minute global relief (SRTM15+V2.1 @ 3.7 km) |
| 01m | 21600 x 10800 | g,p | 258 MB | 1 arc minute global relief (SRTM15+V2.1 @ 1.9 km) |
| 30s | 43200 x 21600 | g,p | 935 MB | 30 arc second global relief (SRTM15+V2.1 @ 1.0 km) |
| 15s | 86400 x 43200 | p | 3.2 GB | 15 arc second global relief (SRTM15+V2.1) |
| 03s | 432000 x 216000 | g | 6.8 GB | 3 arc second global relief (SRTM3S) |
| 01s | 1296000 x 432000 | g | 41 GB | 1 arc second global relief (SRTM1S) |
おすすめは、15s 03s 01s。
かなり広い範囲を描画する場合は30s 01mも検討してみてください。
GMT上では @earthrelief【code】と指定してください。今回は15sを選択します。
DEMを切り出す (grdcut)
gmt grdcut [使用するDEM] -R[切り出し範囲] -G[出力ファイル名]
DEMを使うためには、まず自分が描きたい範囲のDEMを切り出す必要があります。
たとえば、
#!/bin/bash
gmt grdcut @earth_relief_15s -R120/150/20/50 -Gdem.ncこれを行うと、dem.ncというDEMが出力されます。
DEMの勾配計算 (grdgradient)
gmt grdgradient [使用するDEM] -G[出力ファイル名] -A[光を当てる角度] -N[正規化値] -G[出力ファイル名]
かっこいい地図を描く場合は勾配計算を行いましょう。
光を当てる角度は45度。正規化値はe0.8がおすすめです。
#!/bin/bash
gmt grdcut @earth_relief_15s -R120/150/20/50 -Gdem.nc
gmt grdgradient dem.nc -Ggrad.grd -A45 -Ne0.8これを行うと、grad.grdという勾配ファイルが出力されます。
DEMを描画 (grdimage)
gmt grdimage [DEMファイル] -I[勾配ファイル] -C
最後にDEMを地図に描画します。
makecptでDEM用のカラーパレットを作成することを忘れずに。
geo とか globe がおすすめ。
#!/bin/bash
gmt begin hogehoge jpg
gmt basemap -JM12 -R120/150/20/50 -Bafg -BWSNE
gmt coast -Df -W0.25 -Lx11/-1.5+c20+w1000+f
gmt makecpt -Cgeo -T-8000/8000/200 -Z
gmt grdcut @earth_relief_15s -R120/150/20/50 -Gdem.nc
gmt grdgradient dem.nc -Ggrad.grd -A45 -Ne0.8
gmt grdimage dem.nc -Igrad.grd -C
gmt end
rm *.nc *.grdこの結果が、
となる。
コメント