GMT(The Generic Mapping Tools)

2012年6月29日 (金)

GMTの学習(SESSION FOUR)

■最後のセッションです
●4.1    Cpt files

■cptファイルの作成
$ makecpt  -Crainbow  -T-20/60/10  >  disc.cpt
$ makecpt  -Crainbow  -T-20/60/10  -Z  >  cont.cpt

■psscaleコマンドを用いてカラーバー(カラースケール)を描く
$ psbasemap  -R0/8.5/0/11  -Jx1i  -P  -B0  -K  >  bar.ps
$ psscale  -D3i/3i/4i/0.5ih  -Cdisc.cpt  -B:discrete:  -O  -K  >>  bar.ps
$ psscale  -D3i/5i/4i/0.5ih  -Ccont.cpt  -B:continuous:  -O  -K  >>  bar.ps
$ psscale  -D3i/7i/4i/0.5ih  -Cdisc.cpt  -B:discrete:  -I0.5  -O  -K  >>  bar.ps
$ psscale  -D3i/9i/4i/0.5ih  -Ccont.cpt  -B:continuous:  -I0.5  -O  >>  bar.ps
【出力】

出力ファイル

●4.1.1    Exercises
1.  Redo the makecpt exercise using the master table hot and redo the bar plot.

$ makecpt  -Chot  -T-20/60/10  >  disc_h.cpt                                  
$ makecpt  -Chot  -T-20/60/10  -Z  >  cont_h.cpt 
                            
$ psbasemap  -R0/8.5/0/11  -Jx1i  -P  -B0  -K  >  bar_b.ps                     
$ psscale  -D3i/3i/4i/0.5ih  -Cdisc_h.cpt  -B:discrete:  -O  -K  >>  bar_b.ps 
$ psscale  -D3i/5i/4i/0.5ih  -Ccont_h.cpt  -B:continuous:  -O  -K  >>  bar_b.ps
$ psscale  -D3i/7i/4i/0.5ih  -Cdisc_h.cpt  -B:discrete:  -I0.5  -O  -K  >>  bar_b.ps
$ psscale  -D3i/9i/4i/0.5ih  -Ccont_h.cpt  -B:continuous:  -I0.5  -O  >>  bar_b.ps
※-I0.5 -Hのパラメータ(陰影効果をつける)の値は-1.0~1.0の範囲か.
【出力】

出力ファイル

2.  Try specifying -B10g5.
$ psbasemap  -R0/8.5/0/11  -Jx1i  -P  -B0  -K  >  bar_c.ps                           
$ psscale  -D3i/3i/4i/0.5ih  -Cdisc_h.cpt  -B10g5:discrete:  -O  -K  >>  bar_c.ps
$ psscale  -D3i/5i/4i/0.5ih  -Ccont_h.cpt  -B10g5:continuous:  -O  -K  >>  bar_c.ps 
$ psscale  -D3i/7i/4i/0.5ih  -Cdisc_h.cpt  -B10g5:discrete:  -I0.5  -O  -K  >>  bar_c.ps
$ psscale  -D3i/9i/4i/0.5ih  -Ccont_h.cpt  -B10g5:continuous:  -I0.5  -O  >>  bar_c.ps
【出力】

出力ファイル

●4.3    Color images
【グリッドファイル us.nc の入手】
・GMTサイトのtutorialディレクトリにあると思います.
例えば,here
ftp://ftp.scc.u-tokai.ac.jp/pub/gmt/windows/GMT_tut.zip

「us.nc」をダウンロード (←ここにも置かせていただきました)

$ grdraster  9  -R-108/-103/35/40  -Gus.nc
は関連ファイルが不明のため未実行です.


######################
●us.ncに相当するグリッドファイル(my_us_nc.grd)を以下のサイトを参考に作成しました(感謝).

※参照サイト
GTOPO30から作成する
http://shimizus.hustle.ne.jp/wiki/wiki.cgi?page=GTOPO30%A4%AB%A4%E9%BA%EE%C0%AE%A4%B9%A4%EB
GTOPO30 Global 30 Arc Second Elevation Data Set
http://www1.gsi.go.jp/geowww/globalmap-gsi/gtopo30/gtopo30.html

$ xyz2grd W140N40.DEM -Dm/m/m/1/0/=/= -GW140N40.grd -R-140/-100/-10/40  -I30c -N-9999 -V -ZTLhw -F

C:\GMT4\bin\xyz2grd.exe: nx = 4800  ny = 6000
C:\GMT4\bin\xyz2grd.exe: Working on file W140N40.DEM


$ grdinfo W140N40.grd


W140N40.grd: Title: W140N40.grd
W140N40.grd: Command: C:\GMT4\bin\xyz2grd.exe -V -R-140/-100/-10/40 -I30c W140N40.DEM -Dm/m/m/1/0/=/= -GW140N40.grd -N-9999 -ZTLhw -F
W140N40.grd: Remark:
W140N40.grd: Pixel node registration used
W140N40.grd: Grid file format: nf (# 18) GMT netCDF format (float)  (COARDS-compliant) [DEFAULT]
W140N40.grd: x_min: -140 x_max: -100 x_inc: 0.00833333333333 name: m nx: 4800
W140N40.grd: y_min: -10 y_max: 40 y_inc: 0.00833333333333 name: m ny: 6000
W140N40.grd: z_min: -79 z_max: 4328 name: m
W140N40.grd: scale_factor: 1 add_offset: 0


$ grdcut W140N40.grd -G
my_us_nc.grd -R-108/-103/35/40
$ grdimage 
my_us_nc.grd  -JM6i  -P  -B2  -Ctopo_43.cpt  -V  -K > my_topo_43.ps

C:\GMT4\bin\grdimage.exe: Allocates memory and read data file
C:\GMT4\bin\grdimage.exe: project grid files
C:\GMT4\bin\grdimage.exe: Evaluate pixel colors
C:\GMT4\bin\grdimage.exe: Creating PostScript image 24-bit [colorimage]
GMT_grd_is_global: no!
GMT_grd_is_global: no!
GMT_grd_is_global: no!
C:\GMT4\bin\grdimage.exe: Grid projection from size 600x600 to 600x600
GMT_grd_is_global: no!
GMT_boundcond_param_prep determined edgeinfo: gn = 0, gs = 0, nxp = 0, nyp = 0
pslib: Too many colors to make colormap - using 24-bit direct color instead.
pslib: LZW compressed 1080000 to 586543 bytes


$ psscale  -D3i/8.5i/5i/0.25ih  -Ctopo_43.cpt  -I0.4  -B/:m:  -O  >>  my_topo_43.ps


※my_topo_43.ps は以下に示すtopo_43.psと同じ図となります.
######################

$ makecpt  -Crainbow  -T1000/5000/500  -Z  >  topo_43.cpt

$ grdimage  us.nc  -JM6i  -P  -B2  -Ctopo_43.cpt  -V  -K  >  topo_43.ps
C:\GMT4\bin\grdimage.exe: Allocates memory and read data file
C:\GMT4\bin\grdimage.exe: project grid files
C:\GMT4\bin\grdimage.exe: Evaluate pixel colors
C:\GMT4\bin\grdimage.exe: Creating PostScript image 24-bit [colorimage]

GMT_grd_is_global: no!
GMT_grd_is_global: no!
GMT_grd_is_global: no!

C:\GMT4\bin\grdimage.exe: Grid projection from size 600x600 to 600x600
GMT_grd_is_global: no!
GMT_boundcond_param_prep determined edgeinfo: gn = 0, gs = 0, nxp = 0, nyp = 0
pslib: Too many colors to make colormap - using 24-bit direct color instead.
pslib: LZW compressed 1080000 to 586543 bytes


※GMT_grd_is_global: no! 内容不明ですが・・・
$ psscale  -D3i/8.5i/5i/0.25ih  -Ctopo_43.cpt  -I0.4  -B/:m:  -O  >>  topo_43.ps


出力ファイル

※リオ・グランド国立森林公園 (コロラド州 37.857507,-106.75415)周辺のようです.

$ grdgradient  us.nc  -Ne0.8  -A100  -M  -Gus_i.nc
$ grdimage  us.nc  -Ius_i.nc  -JM6i  -P  -B2  -Ctopo_43.cpt  -K  >  topo_43_b.ps
$ psscale  -D3i/8.5i/5i/0.25ih  -Ctopo_43.cpt  -I0.4  -B/:m:  -O  >>  topo_43_b.ps

※GMT_grd_is_global: no! 内容不明ですが・・・
$ grdimage  us.nc  -Ius_i.nc  -JM6i  -P  -B2  -Ctopo_43.cpt  -V -K  >  topo_43_b.ps
とすると,GMT_grd_is_global: no! が表示される.
【出力】

出力ファイル

●4.3.1    Exercises
1.  Force a gray-shade image.
$ makecpt  -Cgray  -T1000/5000/500  -Z  >  topo_gray.cpt
$ grdimage  us.nc  -Ius_i.nc  -JM6i  -P  -B2  -Ctopo_gray.cpt  -K  >  topo_gray.ps
$ psscale  -D3i/8.5i/5i/0.25ih  -Ctopo_gray.cpt  -I0.4  -B/:m:  -O  >>  topo_gray.ps
出力ファイル
※GMT_grd_is_global: no!
2.  Rerun grdgradient with -N1.
$ grdgradient  us.nc  -N1  -A100  -M  -Gus_i_b.nc                              
$ grdimage  us.nc  -Ius_i_b.nc  -JM6i  -P  -B2  -Ctopo_43.cpt  -K  >  topo_43_c.ps
$ psscale  -D3i/8.5i/5i/0.25ih  -Ctopo_43.cpt  -I0.4  -B/:m:  -O  >>  topo_43_c.ps
出力ファイル
※GMT_grd_is_global: no!
●4.4    Multi-dimensional maps
otemp.anal1deg.ncの入手】
ftp://ftp.cdc.noaa.gov/Datasets/nodc.woa98/temperat/seasonal/otemp.anal1deg.nc

$ ncdump -h otemp.anal1deg.nc
netcdf otemp.anal1deg {
dimensions:
        lon = 360 ;
        lat = 180 ;
        level = 33 ;
        time = UNLIMITED ; // (4 currently)
variables:
        float lat(lat) ;
                lat:units = "degrees_north" ;
                lat:long_name = "Latitude" ;
                lat:actual_range = 89.5f, -89.5f ;
                lat:standard_name = "latitude" ;
                lat:axis = "Y" ;
        float lon(lon) ;
                lon:units = "degrees_east" ;
                lon:long_name = "Longitude" ;
                lon:actual_range = 0.5f, 359.5f ;
                lon:standard_name = "longitude" ;
                lon:axis = "X" ;
        float level(level) ;
                level:units = "meters" ;
                level:positive = "down" ;
                level:long_name = "Level" ;
                level:actual_range = 0.f, 5500.f ;
                level:axis = "Z" ;
        double time(time) ;
                time:units = "days since 1-1-1 00:00:0.0" ;
                time:long_name = "Time" ;
                time:actual_range = 0., 273. ;
                time:delta_t = "0000-03-00 00:00:00" ;
                time:avg_period = "0097-00-00 00:00:00" ;
                time:prev_avg_period = "0000-01-00 00:00:00" ;
                time:ltm_range = 693597., 729026. ;
                time:standard_name = "time" ;
                time:axis = "T" ;
        float otemp(time, level, lat, lon) ;
                otemp:long_name = "Ocean Temperature, analyzed mean, 1-deg grid, Seasonal" ;
                otemp:valid_range = -10.f, 90.f ;
                otemp:actual_range = -2.1f, 32.7039f ;
                otemp:units = "deg C" ;
                otemp:add_offset = 0.f ;
                otemp:scale_factor = 1.f ;
                otemp:missing_value = -9.96921e+36f ;
                otemp:var_desc = "Ocean Temperature, analyzed mean" ;
                otemp:dataset = "NODC World Ocean Atlas 1998" ;
                otemp:level_desc = "Multiple Levels" ;
                otemp:statistic = "Analyzed Mean" ;
                otemp:parent_stat = "Mean" ;

// global attributes:
                :title = "NODC World Ocean Atlas 1998" ;
                :history = "created February 2000 by Hoop" ;
                :platform = "Marine Analyses" ;
                :Conventions = "COARDS" ;
                :references = "http://www.esrl.noaa.gov/psd/data/data.nodc.woa98.html" ;
}

$ makecpt -Cno_green -T-2/30/2 > otemp.cpt
$ grdimage -Rg -JW180/9i "otemp.anal1deg.nc?otemp[2,0]" -Cotemp.cpt -Bg30 > otemp.ps

出力ファイル

●4.4.1    Exercises
1.  Plot the temperatures for Spring at 5000 m depth. (Hint: use ncdump -v level to figure out what level number that is).
$ ncdump -v level  otemp.anal1deg.nc
netcdf otemp.anal1deg {
dimensions:
        lon = 360 ;
        lat = 180 ;
        level = 33 ;
        time = UNLIMITED ; // (4 currently)
variables:
        float lat(lat) ;
                lat:units = "degrees_north" ;
                lat:long_name = "Latitude" ;
                lat:actual_range = 89.5f, -89.5f ;
                lat:standard_name = "latitude" ;
                lat:axis = "Y" ;
        float lon(lon) ;
                lon:units = "degrees_east" ;
                lon:long_name = "Longitude" ;
                lon:actual_range = 0.5f, 359.5f ;
                lon:standard_name = "longitude" ;
                lon:axis = "X" ;
        float level(level) ;
                level:units = "meters" ;
                level:positive = "down" ;
                level:long_name = "Level" ;
                level:actual_range = 0.f, 5500.f ;
                level:axis = "Z" ;
        double time(time) ;
                time:units = "days since 1-1-1 00:00:0.0" ;
                time:long_name = "Time" ;
                time:actual_range = 0., 273. ;
                time:delta_t = "0000-03-00 00:00:00" ;
                time:avg_period = "0097-00-00 00:00:00" ;
                time:prev_avg_period = "0000-01-00 00:00:00" ;
                time:ltm_range = 693597., 729026. ;
                time:standard_name = "time" ;
                time:axis = "T" ;
        float otemp(time, level, lat, lon) ;
                otemp:long_name = "Ocean Temperature, analyzed mean, 1-deg grid, Seasonal" ;
                otemp:valid_range = -10.f, 90.f ;
                otemp:actual_range = -2.1f, 32.7039f ;
                otemp:units = "deg C" ;
                otemp:add_offset = 0.f ;
                otemp:scale_factor = 1.f ;
                otemp:missing_value = -9.96921e+36f ;
                otemp:var_desc = "Ocean Temperature, analyzed mean" ;
                otemp:dataset = "NODC World Ocean Atlas 1998" ;
                otemp:level_desc = "Multiple Levels" ;
                otemp:statistic = "Analyzed Mean" ;
                otemp:parent_stat = "Mean" ;

// global attributes:
                :title = "NODC World Ocean Atlas 1998" ;
                :history = "created February 2000 by Hoop" ;
                :platform = "Marine Analyses" ;
                :Conventions = "COARDS" ;
                :references = "http://www.esrl.noaa.gov/psd/data/data.nodc.woa98.html" ;
data:

level = 0, 10, 20, 30, 50, 75, 100, 125, 150, 200, 250, 300, 400, 500, 600,
    700, 800, 900, 1000, 1100, 1200, 1300, 1400, 1500, 1750, 2000, 2500,
    3000, 3500, 4000, 4500,
5000, 5500 ;
}


5000mは配列要素番号[32]とした.配列添え字は0から数えた.
Springは[1]と推測したが・・・

$ makecpt -Cno_green -T-2/30/2 > otemp.cpt
$ grdimage -Rg -JW180/9i "otemp.anal1deg.nc?otemp[1,32]" -Q -K -Cotemp.cpt -Bg30 > otemp_s5000.ps
$ psscale  -D4.5i/-0.25i/5i/0.25ih  -Cotemp.cpt  -I0.4  -B/:degC:  -O  -K >>  otemp_s5000.ps
$ pscoast  -Rg -JW180/9i -Gchocolate  -O >>  otemp_s5000.ps      
                        
出力ファイル

2.  Include a color scale at the bottom of the plot.
$ makecpt -Cno_green -T-2/30/2 > otemp.cpt
$ grdimage -Rg -JW180/9i "otemp.anal1deg.nc?otemp[2,0]" -Cotemp.cpt -Bg30 -K > otemp_b.ps
$ psscale  -D4.5i/-0.25i/5i/0.25ih  -Cotemp.cpt  -I0.4  -B/:C:  -O  >>   otemp_b.ps

出力ファイル

●4.5.1    Mesh-plot
$ grd2cpt  bermuda.nc  -Cocean  -V >  bermuda.cpt   
C:\GMT4\bin\grd2cpt.exe:  Mean and S.D. of data are -4651.91174489 497.003975175
C:\GMT4\bin\grd2cpt.exe: z = -5475 and CDF(z) = 0
C:\GMT4\bin\grd2cpt.exe: z = -5288.84718928 and CDF(z) = 0.0253818508535
C:\GMT4\bin\grd2cpt.exe: z = -5070.20023048 and CDF(z) = 0.0500898472597
C:\GMT4\bin\grd2cpt.exe: z = -4912.54062947 and CDF(z) = 0.200808625337
C:\GMT4\bin\grd2cpt.exe: z = -4777.827702 and CDF(z) = 0.454851752022
C:\GMT4\bin\grd2cpt.exe: z = -4651.91174489 and CDF(z) = 0.570530098832
C:\GMT4\bin\grd2cpt.exe: z = -4525.99578778 and CDF(z) = 0.789083557951
C:\GMT4\bin\grd2cpt.exe: z = -4391.28286031 and CDF(z) = 0.908580413297
C:\GMT4\bin\grd2cpt.exe: z = -4233.6232593 and CDF(z) = 0.940925426774
C:\GMT4\bin\grd2cpt.exe: z = -4014.97630051 and CDF(z) = 0.959568733154
C:\GMT4\bin\grd2cpt.exe: z = -89 and CDF(z) = 1
C:\GMT4\bin\grd2cpt.exe: Reading colortable C:/GMT4/share/cpt/GMT_ocean.cpt


$ cat bermuda.cpt
#       cpt file created by: C:\GMT4\bin\grd2cpt.exe -V bermuda.nc -Cocean
#COLOR_MODEL = RGB
#
-5475   0       2       10      -5288.84718928  0       2       10
-5288.84718928  0       6       30      -5070.20023048  0       6       30
-5070.20023048  0       10      50      -4912.54062947  0       10      50
-4912.54062947  0       66      110     -4777.827702    0       66      110
-4777.827702    0       122     170     -4651.91174489  0       122     170
-4651.91174489  34      169     194     -4525.99578778  34      169     194
-4525.99578778  103     207     181     -4391.28286031  103     207     181
-4391.28286031  172     245     168     -4233.6232593   172     245     168
-4233.6232593   203     249     202     -4014.97630051  203     249     202
-4014.97630051  234     253     237     -89     234     253     237
B       0       0       0
F       255     255     255
N       128     128     128


$ grdview  bermuda.nc  -JM5i  -P  -JZ2i  -E135/30  -B2  -Cbermuda.cpt  >  bermuda_map.ps
出力ファイル
ワイヤフレームの表示になる.

 $ grdview  bermuda.nc  -JM5i  -P  -JZ2i  -E135/30  -B2  -Cbermuda.cpt  -Qs -X1.5c -V>  bermuda_map_b.ps
C:\GMT4\bin\grdview.exe: Processing shape file
C:\GMT4\bin\grdview.exe: Find contours
C:\GMT4\bin\grdview.exe: Trace and bin contours...
C:\GMT4\bin\grdview.exe: Now tracing contour interval      -89
C:\GMT4\bin\grdview.exe: Start creating PostScript plot
C:\GMT4\bin\grdview.exe: Done!
GMT_grd_is_global: no!
GMT_grd_is_global: no!
GMT_grd_is_global: no!
GMT_boundcond_param_prep determined edgeinfo: gn = 0, gs = 0, nxp = 0, nyp = 0
GMT_grd_is_global: no!

出力ファイル
3-D表面塗り,

-X1.5c 描画位置調整

$ grdview  bermuda.nc  -JM5i  -P  -JZ2i  -E135/30  -B2  -Cbermuda.cpt  -Qsm -X1.5c -W0.25p,. -V>  bermuda_map_c.ps
C:\GMT4\bin\grdview.exe: Processing shape file
C:\GMT4\bin\grdview.exe: Find contours
C:\GMT4\bin\grdview.exe: Trace and bin contours...
C:\GMT4\bin\grdview.exe: Now tracing contour interval      -89
C:\GMT4\bin\grdview.exe: Start creating PostScript plot
C:\GMT4\bin\grdview.exe: Done!
GMT_grd_is_global: no!
GMT_grd_is_global: no!
GMT_grd_is_global: no!
GMT_boundcond_param_prep

出力ファイル
3-D表面塗り,メッシュプロット(点線)

※3-D塗り分け,メッシュプロット,メッシュプロットの線種などのオプションの設定
-Qm/cyan

-Qs
-Qsm

●Exercises
1.  Select another vantage point and vertical height.

$ grdview  bermuda.nc  -JM5i  -P  -JZ2i  -E135/30  -B2  -Cbermuda.cpt  >  bermuda_map.ps
または
$ grdview  bermuda.nc  -JM5i  -P  -JZ2i -E135/30  -B2  -Cbermuda.cpt  -Qsm -X1.5c -W0.25p,. -V>  bermuda_map_c.ps

※赤文字部を変えてみよう.
視点の位置(-E135/30)および垂直スケール(-JZ2i)を変更してみましょう.

●4.5.2    Color-coded view
$ grdview  us.nc  -JM6i  -E135/35  -Qi50  -Ius_i.nc  -Ctopo_43.cpt  -V  -B2  -JZ0.5i  >  view.ps
出力ファイル
※GMT_grd_is_global: no! 気になるが内容不明

$ grdview  us.nc  -JM6i  -E135/35  -Qi200  -Ius_i.nc  -Ctopo_43.cpt  -V  -B2  -JZ0.5i  >  view_200dpi.ps
【出力ファイル】サイズが大きいため(3994KB) 掲載なし
※GMT_grd_is_global: no!

●Exercises
1.  Choose another vantage point and scaling.
$ grdview  us.nc  -JM6i  -E135/15  -Qi50  -Ius_i.nc  -Ctopo_43.cpt  -V  -B2  -JZ0.5i  >  view15.ps
出力ファイル
※-JZ0.5i も変えてみよう.
※GMT_grd_is_global: no!

2.  Redo grdgradient with another illumination direction and replot.
【4.3  Color images】で実行した以下のパラメータを変えてみる.未確認
$ grdgradient  us.nc  -Ne0.8  -A100  -M  -Gus_i.nc

3.  Select a higher dpi, e.g., 200.
$ grdview  us.nc  -JM6i  -E135/15  -Qi200  -Ius_i.nc  -Ctopo_43.cpt  -V  -B2  -JZ0.5i  >  view15_200dpi.ps
出力ファイル
※GMT_grd_is_global: no!

※grdview,grdimageをus.ncに対して実行するとGMT_grd_is_global: no!が表示される・・・
内容や対応わからず.
※psファイルをpdfファイルに変換
$ ps2pdf sample.ps sample.pdf

■終わりに
理解できていない箇所ありますが,学習履歴とします.
使い込んでいないため,誤りが多々あると思います.ご容赦ください.
利用できそうなコマンドをみつけて作図に役立てたいと思っています.
緑色の文字の箇所は未確認・棚上げ中などの項目です.
地震や原発事故による土壌,海域,山林,大気など環境の汚染などの
可視化に役立つソフトと思います.できるところから使ってみましょう.
【謝辞】
多くのGMT関係サイトを参考にさせて頂きました.
サイト掲載の皆様に感謝いたします.




2012年6月27日 (水)

GMTの学習(SESSION THREE)

■PC環境によって変更が必要になることがあると思いますので注意しましょう.

●3.1    Contouring gridded data sets
1) etopo5.i2
2) grdraster.info
の2つのファイルをC:\GMT4\share\dbase におく.

【入手先】
1) http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO5/TOPO/ETOPO5/
ETOPO5.DATをetopo5.i2とファイル名を変更するようです.
※このデータでは実行できなかった.理由はわかりません(grdraster.infoが不適だったか).
2) http://www.volcanoloco.org/nico/Mac_OS_X/Entries/2008/11/30_Generic_Mapping_Tool%3A_I_-_Installation.html
※参考
3) http://www.geophysik.uni-muenchen.de/~oeser/LV/GMT/
※使用

ここでは,3) からの二つのファイルをDLして C:\GMT4\share\dbase におく.

※環境によりうまく実行できない場合もあるかもしれません.

■ bermuda.nc の作成  (ファイル名をbermuda.grdとしてもよいでしょう)
以上の準備のもとで,次を実行するとグリッドファイルbermuda.ncが生成される.
$ grdraster  1  -R-66/-60/30/35  -Gbermuda.nc  -V

-V 実行の進捗状況を画面に表示
※バーミューダ海域

■ bermuda.nc グリッドファイルの内容についての情報表示

$ grdinfo  bermuda.nc
bermuda.nc: Title: ETOPO5 global topography
bermuda.nc: Command: C:\GMT4\bin\grdraster.exe -V -R-66/-60/30/35 1 -Gbermuda.nc
bermuda.nc: Remark: C:/GMT4/share/dbase/etopo5.i2
bermuda.nc: Gridline node registration used
bermuda.nc: Grid file format: nf (# 18) GMT netCDF format (float)  (COARDS-compl    iant) [DEFAULT]
bermuda.nc: x_min: -66 x_max: -60 x_inc: 0.0833333333333 name: longitude [degree    s_east] nx: 73
bermuda.nc: y_min: 30 y_max: 35 y_inc: 0.0833333333333 name: latitude [degrees_n    orth] ny: 61
bermuda.nc: z_min: -5475 z_max: -89 name: m
bermuda.nc: scale_factor: 1 add_offset: 0


※bermuda.ncまたはbermuda.grdはGMTのサイトから入手できると思います.
直接入手したものをここにおきます.

■グリッドデータから等値線図を描く
$ grdcontour  bermuda.nc  -JM6i  -C250  -A1000  -P  -B2  > map.ps
-C250  コンター間隔 (m)
-A1000 等値線に描く数値の間隔(m)
【出力】

出力ファイル

●3.1.1    Exercises 解答例
1.  Add smoothing with -S4.
コンターをスムージング
$ grdcontour  bermuda.nc  -JM6i  -C250  -A1000  -P  -B2 -S4 > map.ps

2.  Try tick all highs and lows with -T.
$ grdcontour  bermuda.nc  -JM6i  -C250  -A1000  -P  -B2 -T > map.ps

3.  Skip small features with -Q10.
$ grdcontour  bermuda.nc  -JM6i  -C250  -A1000  -P  -B2 -Q10 > map.ps

4.  Override region using -R-70/-60/25/35.
$ grdraster  1  -R-70/-60/25/35  -Gover_bermuda.grd  -V
$ grdcontour  over_bermuda.grd  -JM6i  -C250  -A1000  -P  -B2  > map.ps
出力ファイル

5.  Try another region that clips our data domain.
$ grdraster  1  -R-75/-55/20/40  -Gano_bermuda.grd  -V
$ grdcontour  ano_bermuda.grd  -JM6i  -C250  -A1000  -P  -B2  > map.ps
または
$ grdraster  1  -R-80/-50/15/45  -Gano_bermuda_b.grd  -V
$ grdcontour  ano_bermuda_b.grd  -JM6i  -C250  -A1000  -P  -B5  > map.ps
出力ファイル