Surfer软件在台站气象服务中的应用.doc

《Surfer软件在台站气象服务中的应用.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、Surfer软件在台站气象服务中的应用　摘   要 本文以绘制亳州市等雨量图为例，着重介绍利用Surfer软件绘制区域底图和在底图上绘制等雨量线图的方法，展示用VB编程语言调用Surfer7.0的绘图函数，实现Surfer与VB之间的ActiveX自动化连接，从而实现区域等值线图的程序化输出。关键词 Surfer  区域底图  气象要素  等值线　引言GoldenSoftwareSurfer（简称Surfer）是一款绘制三维图和等值线图的软件，最大的特色是具有插值功能和填色绘制。用此软件绘制的区域等雨量线图能直接形象地反映出区域雨

2、量的分布。虽然有很多文献对Surfer软件都进行了较为详尽的介绍，但由于其帮助文件大多是英文版，因此在基层台站的应用还待进一步提高。如何制作区域底图，并将区域自动站雨量数据叠加在底图上，以及怎样利用VB编程语言调用Surfer7.0的绘图函数，实现程序化操作，是本文所要解决的主要问题。　1  资料说明亳州市区域地图和乡镇自动雨量站点的经纬度资料取自Micaps1.0中的地理信息文件。本文以2007年7月2日20时至2007年7月15日20时降水资料（取自安徽省高密度自动观测站网统计数据）为例，介绍亳州市区域雨量图的绘制方法。　2 

3、 Surfer7.0软件绘制区域等雨量线图绘制区域等雨量线图前，必须做好前期基础工作。如生成区域底图文件、站点坐标文件和区域雨量数据文件等。底图文件必须要做两个，一个是闭合区域无县界底图文件，一个是闭合区域分县界底图文件，前者用于Blank等值线，后者作为区域等值线背景图用。2.1  制作区域底图文件绘制区域底图，有两种方法。一是读取区域地图边界上点的经纬度坐标，生成bln文件。另一种方法利用该地区地图，用Surfer软件在map

4、basemap中打开地图文件即可。虽然，第二种方法较为简单，但下载本地区地图并不太容易。这里介绍第一

5、种方法。2.1.1   制作闭合区域无县界底图文件首先用Micaps1.0读出亳州边界经纬度资料，生成闭合无县界的bln文件Bz1.bln。Bz1.bln文件格式如下：Length，flagX1，Y1X2，Y2…….Xn，Yn其中，X1，Y1是闭合曲线初始点的经纬度数据，Xn，Yn是闭合曲线终点的经纬度数据。X是某一点上的纬度数据，Y是某一点上的经度数据。Length是底图曲线中点坐标的个数，当flag＝1时，底图区域内的等高线可以消除掉，当flag＝0时，则可以消除底图区域以外的等高线，本例选flag值为0。注意：一定要是闭合等

6、值线图，所读经纬度数据要首尾相连。2.1.2   制作区域分县界底图生成亳州分县界的bln文件Bz2.bln，方法基本与绘制闭合区域底图相同，但要保证Bz2.bln与Bz1.bln底图外围边界经纬度数据完全相同（可直接粘贴Bz1.bln内数据），在读取经纬度坐标数据时应沿一个方向一次性读取，首尾相连，避免重复。不得已重复时一定要与原来所读的数据相同，否则，底图曲线也会出现重影，影响底图质量。2.1.3   制作站点坐标文件站点坐标文件Bzxz.dat，文件格式：经度1，纬度1，站点名1经度2，纬度2，站点名2…经度n，纬度n，站点

7、名n2.1.4    制作雨量数据文件制作亳州区域自动雨量站雨量数据文件Bzyl.dat，文件格式：经度1，纬度1，雨量1，站点名1经度2，纬度2，雨量2，站点名2…经度n，纬度n，雨量n，站点名n2.2  区域等雨量线图的制作（以Surfer7.0为例）具体步骤：2.2.1   打开Surfer7.0进入主界面（见图1）。　　　　　　　　　　　　　　　　　　　　　2.2.2   将雨量数据文件转化成Surfer7.0所能识别的grd格式文件点击菜单项Grid、Data后，打开区域雨量数据文件Bzyl.dat，出现Seattere

8、dDataTnterpolation对话框（见图2）：　　　　　　其中GridlineGeometry选项X、Y分别表示区域底图的最大外边界经纬度值。本例X选115.54、116.82；Y选32.86、34.08，在GriddingMethodtfad选项框中，选Kriging（插值方法），设置好后，点击“确定”按钮，提示已生成grd格式的文件bzyl.grd，再按“确定”按钮，随后出现Surfer-Report1对话框，关掉对话框不保存文件即可。2.2.3　Blank文件点击菜单项Grid、Blank，打开bzyl.grd，再选

9、中区域无县界闭合图地图文件Bz1.bln，生成out.grd文件。2.2.4　生成等雨量线图从菜单项选Map、contour Map、NewcontourMap，打开out.grd文件，出现ContourMapProperties对话框（见图3），