grdcut コマンド メモ †
grd ファイルの範囲切り取りをする．
公式doc †
例 †
topofile="/usr/share/geoware/ttt/data/ttt_topo_5m.i2" topofile_cut="cut.grd" region=120/150/20/50 # 切り取り gmt grdcut $topofile -R$region -G$topofile_cut # 描画で確認 gmt grdimage -Baf -Cetopo1 -JX12/6 -R$topofile $topofile > topo_org.ps # 切り取り前 gmt grdimage -Baf -Cetopo1 -JX12/12 -R$topofile_cut $topofile_cut > topo_cut.ps # 切り取り後
概要 †
基本形は以下の通り．
gmt grdcut [input_file] -R[region] -G[output_file]
-G †
出力ファイル名．surface などの -G オプションと同じ．
-R †
GMT の共通オプション．他の grd ファイルと同一の範囲に切り取りたければこの -R[filename] とする．
