grads 读取shp

自从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'
' set mpdraw off '
' set lon 70 140'
' set lat 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,有兴趣的可以测试一下啦

你可能感兴趣的:(HP)