日志
（转载）绘制地形图并标出研究区域 http://bbs.06climate.com/forum.php?mod=viewthre ...
绘制地形图并标出研究区域
http://bbs.06climate.com/forum.php?mod=viewthread&tid=54483&fromuid=53836
(出处: 气象家园)
grads 画图，非常好的工作
第一贴出图像。
第二贴出脚本。其中用红框将（108-119E 19-25N），（116-121E 26-29N）和（112-122E 29-33N）三个研究区域
标出（见上图）。
'reinit'
'set vpage 0.0 11.0 0.0 8.5'
'set parea 1.0 10.5 1.5 8.0'
'set grads off'
'set grid off'
'open topo10.ctl'
'set dfile 1'
'set x 600 781'
'set y 630 781'
'set vrange 15 40'
'set vrange2 100 130'
'set ylint 5'
'set xlint 5'
'set ylopts 1 4 0.18'
'set xlopts 1 4 0.18'
'set gxout shaded'
'set rgb 20 24 116 205'
'set rgb 21 51 204 255'
'set rgb 22 0 255 255'
'set rgb 23 0 255 100'
'set rgb 24 0 255 0'
'set rgb 25 90 149 44'
'set rgb 26 155 207 52'
'set rgb 27 173 223 80'
'set rgb 28 197 239 116'
'set rgb 29 229 243 167'
'set rgb 30 249 242 191'
'set rgb 31 243 234 169'
'set rgb 32 231 216 133'
'set rgb 33 218 197 93'
'set rgb 34 210 174 56'
'set rgb 35 203 153 28'
'set rgb 36 167 116 37'
'set rgb 37 131 84 40'
'set rgb 38 100 64 30'
'set rgb 39 82 54 25'
'set rgb 40 122 100 79'
'set clevs 0 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500'
'set ccols 20 21 22 23 24 25 26 27 28 29 30 33 35 36 37 38 39 40'
'd tp'
'cbarn 1 0 5.7 0.35'
*添加经纬度标注
'set string 1 c 5 90'
'set strsiz 0.24'
'draw string 1.5 4.9 `1Latitude'
'set string 1 c 5 0'
'set strsiz 0.24'
'draw string 6.0 0.9 `1Longitude'
*将大海填成蓝色（20号颜色）
'set mproj latlon'
'set mpdset cnworld'
'set map 1 1 4'
'draw map'
'basemap O 20 1 M'
*用红框标出三个研究区域
'q w2xy 108 19'
x1=subwrd(result,3)
y1=subwrd(result,6)
'q w2xy 119 25'
x2=subwrd(result,3)
y2=subwrd(result,6)
'set line 2 1 6'
'draw rec 'x1' 'y1' 'x2' 'y2''
'q w2xy 116 26'
x3=subwrd(result,3)
y3=subwrd(result,6)
'q w2xy 121 29'
x4=subwrd(result,3)
y4=subwrd(result,6)
'set line 2 1 6'
'draw rec 'x3' 'y3' 'x4' 'y4''
'q w2xy 112 29'x5=subwrd(result,3)
y5=subwrd(result,6)
'q w2xy 122 33'
x6=subwrd(result,3)
y6=subwrd(result,6)
'set line 2 1 6'
'draw rec 'x5' 'y5' 'x6' 'y6''
'printim topo10.jpg x1000 y800 white'
'print -R topo10.eps'
'disable print'
;
第三给出附件。给出用到的topo10.ctl，topo10.bin，cnworld和basemap.gs文件，并给出原始脚本topo10.gs。
topo10.gs (1.83 KB, 下载次数: 133) topo10.ctl (330 Bytes, 下载次数: 173) topo10.bin (8.9 MB, 下载次数: 1019) cnworld (188.44 KB, 下载次数: 92) basemap.gs (6.91 KB, 下载次数: 105)