上文中生成的等值线村子一条曲线,通过放大gs中坐标范围可以发现,这是默认的世界底图的一部分。可以通过在gs文件中添加'set map draw off' 使底图不显示。但我要绘制的是营口地区的降雨量等值线图,所以需要营口市底图。而Grads中并没有市级底图,所以需要
自己动手绘制。
一种方式是通过Grads命令画线生成边界,需要准备边界数据文件。我的是通过在surfer中有shapefile底图得来,数据如下,供598个点:
把底图坐标点读入的gs文件如下:
最后在等值线gs文件中加入如下'run D:/YK/Board.gs'即可。生成效果图如下:
第二种方式是通过某些软件直接生成等值线,如借助MeteoInfo软件(在动力论坛上),使用说明如下:
输出GrADS地图数据示例(生成全球各国行政界限地图):
1、打开MeteoInfo,添加'country1.shp'图层。
2、点击'Tools'菜单下的'Output Map Data'子菜单,打开地图数据输出对话框。
3、'Map Layer'选'Map_country1.shp',选择所有图元,'Output Format'选择'GrADS Map File'。
4、点击'Output'按钮,在另存为对话框中设定输出文件名即可。
第三种方式是通过MapInfo得到MIF格式文件,最后写fortran程序生成(本人对MapInfo不熟,没亲自做过),在动力论坛上也有此贴。