« 2012年12月 | トップページ | 2013年5月 »
WindowsでGMT第9回、前回からの続きです。
連載のまとめはこちら(外部リンク)です。
ブログ内のカテゴリーGMT(地図ツール)からでも一括でみることができます。
前回のgridimageでの色の塗り分けの際、GMTをインストールしたフォルダの中の\share\cptに標準で入っているカラーパレット(GMT_???.cpt)を使用しました。このカラーパレットは範囲があらかじめ決まっているため、範囲外の色は直近の範囲内の色で塗りつぶれていました。たとえばGMT_ocean.cptでは範囲が-8,000m~0mのため、プラスである陸地はすべて0mの色(=白色)で塗りつぶされます。makecptコマンドを使用すると、カラーパレットについて、色調はそのままで、範囲のみを変えることが可能です。
1.makecptでカラーパレットを作って使用する。
前回の最後のバッチファイルの先頭にmakecptコマンドを追加し、途中のパス記載を変更します。
----------------------------------
c:\programs\GMT4\bin\makecpt -Cc:\programs\GMT4\share\cpt\GMT_ocean.cpt -T-10000/8000/100 > GMT_ocean_change.cpt(改行)
c:\programs\GMT4\bin\grdcut c:\programs\GMT4\share\dbase\ETOPO1_Ice_g_gmt4.grd -R120/155/20/50 -Gjapan.grd
c:\programs\GMT4\bin\grdimage japan.grd -JM12c -Bg5a10f5:."Map": -CGMT_ocean_change.cpt -Y5c -P -K>mapbat.ps(改行)
c:\programs\GMT4\bin\psscale -Ba4000g2000f1000:depth(m): -CGMT_ocean_change.cpt -D8c/-1c/9c/0.3ch -O -K>>mapbat.ps(改行)
c:\programs\GMT4\bin\pscoast -JM12c -R120/155/20/50 -W -O>>mapbat.ps(改行)
----------------------------------
1行目でつくったカラーパレットは作業中のフォルダに保存されていますので、gridimage,psscalceでパレットを指定する際はフォルダの指定(d:\data\gmt)を省略できます。
makecptコマンドのオプションは以下の通りです。
-Cオプション:基になるカラーパレットファイルを指定する。
-Tオプション:新しい範囲と段階を記載する(下限値/上限値/段階)。ここでは-10,000mから+8000mまで100段階で変化させる。
2.色調を反転させる。-Iオプション
makecptでIオプションを使うと色調が反転します。
----------------------------------
c:\programs\GMT4\bin\makecpt -Cc:\programs\GMT4\share\cpt\GMT_ocean.cpt -T-10000/8000/100 -I> GMT_ocean_change.cpt(改行)
c:\programs\GMT4\bin\grdcut c:\programs\GMT4\share\dbase\ETOPO1_Ice_g_gmt4.grd -R120/155/20/50 -Gjapan.grd
c:\programs\GMT4\bin\grdimage japan.grd -JM12c -Bg5a10f5:."Map": -CGMT_ocean_change.cpt -Y5c -P -K>mapbat.ps(改行)
c:\programs\GMT4\bin\psscale -Ba4000g2000f1000:depth(m): -CGMT_ocean_change.cpt -D8c/-1c/9c/0.3ch -O -K>>mapbat.ps(改行)
c:\programs\GMT4\bin\pscoast -JM12c -R120/155/20/50 -W -O>>mapbat.ps(改行)
----------------------------------
3.色調をなめらかにする-Zオプション
-Zオプションをつけると色の変化が連続的になります。
----------------------------------
c:\programs\GMT4\bin\makecpt -Cc:\programs\GMT4\share\cpt\GMT_ocean.cpt -T-10000/8000/100 -I -Z> GMT_ocean_change.cpt(改行)
c:\programs\GMT4\bin\grdcut c:\programs\GMT4\share\dbase\ETOPO1_Ice_g_gmt4.grd -R120/155/20/50 -Gjapan.grd
c:\programs\GMT4\bin\grdimage japan.grd -JM12c -Bg5a10f5:."Map": -CGMT_ocean_change.cpt -Y5c -P -K>mapbat.ps(改行)
c:\programs\GMT4\bin\psscale -Ba4000g2000f1000:depth(m): -CGMT_ocean_change.cpt -D8c/-1c/9c/0.3ch -O -K>>mapbat.ps(改行)
c:\programs\GMT4\bin\pscoast -JM12c -R120/155/20/50 -W -O>>mapbat.ps(改行)
----------------------------------
日 | 月 | 火 | 水 | 木 | 金 | 土 |
---|---|---|---|---|---|---|
1 | 2 | 3 | 4 | 5 | 6 | 7 |
8 | 9 | 10 | 11 | 12 | 13 | 14 |
15 | 16 | 17 | 18 | 19 | 20 | 21 |
22 | 23 | 24 | 25 | 26 | 27 | 28 |
29 | 30 | 31 |