自从GrADS2.0.a8版本开始，GrADS引入了对shp图形的支持，关于此格式在这里不多说，
于是今晚就简单测试了一下最简单画图和查询命令（后续还将测试输出shp图形的命令）
测试数据采用的是160站降水的站点资料，已经写为grads可用的二进制格式，
1、画图
整个过程相当简单，就是一个命令的使用而已（之前我以为很复杂一直没有信心去测试）：
draw shp shapefile <n>
<m>
其中的
shapefile就是shp文件的完整文件名（如果改路径已经加入grads环境变量的GADDIR中则可以直接输入文件名）；
n：可选参数，表示需要画的元素的序号
m：可选参数，和n一起用表示画从n到m范围的元素。
下面是一个gs文件示例：
‘reinit‘
‘open E:\projects\rain.ctl‘
‘open E:\projects\china_mask.ctl‘
*‘set mpdset
hires‘
‘setmpdrawoff‘
‘setlon 70 140‘
‘setlat 15 55‘
‘d
maskout(oacres(mask.2,rain),mask.2(t=1))‘
‘draw shp e:\projects\china.shp‘
‘printim E:\projects\test.png white‘
;
需要注意的有：
1、在画图时需要先关闭默认的地图输出；
2、如果mask文件的范围输出时选择要和你数据的ctl一致（时间维除外），否则会出错；
3、设置好需要显示的经纬度范围；
4、先d后draw shp。
输出的图像如下：
test.png (20.27 KB, 下载次数: 14)
下载附件 保存到相册
2011-9-6 20:48 上传
注：地图的颜色等控制方法感兴趣的可以参考grads官网，后面有空再翻译为中文给大家，建议看英文原版，然后来这里测试讨论啦
2、查询
查询的命令有两种：q shp shpfile和q dbf shpfile，有兴趣的可以测试一下啦