GMTサンプルスクリプトページ@東京海洋大学

HOME | sample.1 | sample.2 | sample.3 | sample.4 | sample.5 | sample.6 | sample.7 | sample.8 | sample.9 | sample.11 |

8. ラスター形式の図を枠に当てはめる
SeaDASで作ったクロロフィルa画像にGMTの枠をあてはめる。
青い部分はあとから重ね合わせたもの(pscoast)。
入力ファイル サンプルシェル 出力ファイル
2001spring-OC4.ras 以下に示す SeaWiFS_chl.ps
#!/bin/csh -f
gmtset BASEMAP_TYPE plain DEGREE_FORMAT 3
gmtset ANOT_FONT 5 ANOT_FONT_SIZE 12 PAGE_COLOR 255/255/255

set image=2001spring-OC4.ras
set ofile=SeaWiFS_chl.ps

set jj=0.4
set rr=175/204/50/65
set xpos=0
set ypos=0 # 2001

psimage $image -W11.6c/11.3c -K -P -X4 -Y15> $ofile      # Sun raster形式の図を読み込む -Wオプションで図の大きさを調整。はっきりいってチカラワザ。
gmtset BASEMAP_FRAME_RGB 0/0/255 FRAME_PEN 2
pscoast -Jm$jj -Df -Ba4f2/a2f1/WESN -R$rr -O -G0/0/255 -X$xpos -Y$ypos>> $ofile

exit