GrADS中shp地图使用初探

mofangbao

测试本文的方法前请确保你的GrADS版本是高于2.0.a8的，如果不是，可以在这里下载最新版的GrADS:
http://www.06climate.com/view/1501.html
========================================
    自从GrADS2.0.a8版本开始，GrADS引入了对shp图形的支持，关于此格式在这里不多说，射月楼主在有一篇帖子中也有提到，同时google一下介绍很多。

    于是今晚就简单测试了一下最简单画图和查询命令（后续还将测试输出shp图形的命令）
    测试数据采用的是160站降水的站点资料，已经写为grads可用的二进制格式，同时预先根据@传说中的谁的一篇帖子（摸我过去看看），使用MeteoInfo生成了mask文件。

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。

输出的图像如下：


注：地图的颜色等控制方法感兴趣的可以参考grads官网，后面有空再翻译为中文给大家，建议看英文原版，然后来这里测试讨论啦

2、查询

查询的命令有两种：q shp shpfile和q dbf shpfile，有兴趣的可以测试一下啦

artmunich

看着就很有技术水准，支持清风，研究后有不懂的再请教啦~

topmad

刚刚从实验室回来，就看到清风的大作了,很好很强大，支持

Mouse

学习下先，回头测试一下

athlon8048

很给力！支持清风！

传说中的谁

先看再学习，学会了给钱

Mouse

athlon8048 发表于 2011-9-6 22:25

                               
登录/注册后可看大图

很给力！支持清风！

@athlon8048这是Grads版第1000个帖子……

喜上眉梢

恩，初步了解，用到时过来细究！感谢老大的孜孜贡献！

zsmnantong

清风不愧是老大，学习一下，顶你！

飞雪残剑

清风，顶，高手就是高手！