GrADS和地理信息系统的结合初探

mofangbao

    很多人都说GrADS跟不上时代了，什么什么才是现在的主流，不管怎么说，GrADS既然能红火这么多年，必然有其优越之处，而且，在最新的GrADS版本中增加的一些功能也紧跟时代的步伐，比如这张图片：

                               
登录/注册后可看大图

趁着端午节放假之前，我今天就抛砖引玉的介绍下GrADS中kml图的生成，这个其实很简单，甚至没有多说的必要。

（请先下载最新版的GrADS或者OpenGrADS，这里下载：GrADS2.0.1.oga.1）

上面这张图的地理地图是ArcGis的在线地图，上面叠加的是一层GrADS绘制的kml数据，kml数据是谷歌定义的一种扩展的xml标记语言，具体的介绍可以这里查看：https://developers.google.com/kml/documentation/

生成KML数据是GrADS2.0.a5版本之后开始支持的，当时可以支持polyline(contour)图，到了2.0.0版之后增加了polygon(shaded)图。

要生成KML文件，先使用

'set gxout kml'
将输出文件类型设置为 kml，注意，不要再出现 set gxout shaded或者contour之类的。
然后使用
'set kml <-type> fnameroot'
来设置是线图还是填色图，fnameroot是生成文件的存放路径。
比如下面这段gs代码就是用来生成上面那幅图的：

1. 'reinit'
   
2. 'open E:\program\example\kml\china.ctl'
   
3. 'open e:\program\example\sta2grd\sta.ctl'
   
4. 'set lon 70 140'
   
5. 'set lat 15 55'
   
6. 'define  rgrid=maskout(oacres(mask,rain.2),mask)'
   
7. 'set gxout kml'
   
8. 'set kml -poly E:\program\example\kml\china.kml'
   
9. 'd rgrid'
   
10. ;
    
11. 
    

'reinit'<br /> 'open E:\program\example\kml\china.ctl'<br /> 'open e:\program\example\sta2grd\sta.ctl'<br /> 'set lon 70 140'<br /> 'set lat 15 55'<br /> 'define  rgrid=maskout(oacres(mask,rain.2),mask)'<br /> 'set gxout kml'<br /> 'set kml -poly E:\program\example\kml\china.kml' <br /> 'd rgrid'<br /> ;<br /> <br />

-poly 选项表示生成多边形，改成 -ln则表示生成线段。（其他选项请点这里看GrADS官网介绍）
运行之后会生成 china.kml文件，下面你就可以拿着这个KML文件去谷歌地球里面试试啦！
如果你没有安装谷歌地球，那也不要紧，我做了一个测试的页面，你在这里打开你的kml文件就可以得到上面的图了。

测试页面地址：http://w.06climate.com




这个数据浏览器现在只是让大家测试一下你的kml！有好的开发建议也可以告诉我。如果使用这个请按照页面提示安装silverlight控件。

最后送大家一份测试数据，里面包含了我生成kml所用到的文件已经生成好的kml文件，你也可以直接先试试。

kml.rar (9.14 KB, 下载次数: 146)

2012-6-20 18:32 上传

点击文件名下载附件

{:soso__10189826289261162925_4:}祝大家端午节快乐！

wanghp

楼主好给你啊，祝楼主端午快乐，顶起

善人/jw

顶起顶起~~~

duan2008

不错 版主又法标了

njzqxt

真是太弓虽大了。

wlzhongouc

图挺好看的  我想任何一种绘图软件用到极致都能绘出很漂亮的图  能精通一种软件就很好

luoziwuhui

学习一下版主的大作！！！！！！！！

夏天爱裤衩

哇哦，终于发帖了，很好很强大

qxtlyf

好漂亮，grads真是与时俱进，祝楼主端午快乐，顶起

做个霸气的木头

感谢楼主，端午节快乐