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

Surfer软件在台站气象服务中的应用 　 ??? 要? 本文以绘制亳州市等雨量图为例，着重介绍利用Surfer软件绘制区域底图和在底图上绘制等雨量线图的方法，展示用VB编程语言调用Surfer7.0的绘图函数，实现Surfer与VB之间的ActiveX自动化连接，从而实现区域等值线图的程序化输出。 关键词? Surfer?? 区域底图?? 气象要素?? 等值线 　 Golden Software Surfer（简称Surfer）是一款绘制三维图和等值线图的软件，最大的特色是具有插值功能和填色绘制。用此软件绘制的区域等雨量线图能直接形象地反映出区域雨量的分布。虽然有很多文献对Surfer软件都进行了较为详尽的介绍，但由于其帮助文件大多是英文版，因此在基层台站的应用还待进一步提高。如何制作区域底图，并将区域自动站雨量数据叠加在底图上，以及怎样利用VB编程语言调用Surfer7.0的绘图函数，实现程序化操作，是本文所要解决的主要问题。 　 1?? 亳州市区域地图和乡镇自动雨量站点的经纬度资料取自Micaps 1.0中的地理信息文件。本文以2007年7月2日20时至2007年7月15日20时降水资料（取自安徽省高密度自动观测站网统计数据）为例，介绍亳州市区域雨量图的绘制方法。 　 2?? Surfer 7.0 绘制区域等雨量线图前，必须做好前期基础工作。如生成区域底图文件、站点坐标文件和区域雨量数据文件等。底图文件必须要做两个，一个是闭合区域无县界底图文件，一个是闭合区域分县界底图文件，前者用于Blank等值线，后者作为区域等值线背景图用。 2.1?? 制作区域底图文件 绘制区域底图，有两种方法。一是读取区域地图边界上点的经纬度坐标，生成bln文件。另一种方法利用该地区地图，用Surfer软件在map|base map中打开地图文件即可。虽然，第二种方法较为简单，但下载本地区地图并不太容易。这里介绍第一种方法。 2.1.1??? 制作闭合区域无县界底图文件 首先用Micaps 1.0读出亳州边界经纬度资料，生成闭合无县界的bln文件Bz1.bln。 Bz1.bln文件格式如下： Length，flag X1，Y1 X2，Y2 ……. Xn，Yn 其中，X1，Y1是闭合曲线初始点的经纬度数据，Xn，Yn是闭合曲线终点的经纬度数据。X是某一点上的纬度数据，Y是某一点上的经度数据。Length是底图曲线中点坐标的个数，当flag＝1时，底图区域内的等高线可以消除掉，当flag＝0时，则可以消除底图区域以外的等高线，本例选flag值为0。 注意：一定要是闭合等值线图，所读经纬度数据要首尾相连。 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，站点名n 2.1.4???? 制作雨量数据文件 制作亳州区域自动雨量站雨量数据文件Bzyl.dat，文件格式： 经度1，纬度1，雨量1，站点名1 经度2，纬度2，雨量2，站点名2 … 经度n，纬度n，雨量n，站点名n 2.2?? 区域等雨量线图的制作（以Surfer 7.0为例） 具体步骤： 2.2.1??? 打开Surfer 7.0进入主界面（见图1）。 　 　 　 　 2.2.2??? 将雨量数据文件转化成Surfer 7.0所能识别的grd格式文件 点击菜单项Grid、Data后，打开区域雨量数据文件Bzyl.dat，出现Seattered Data Tnterpolation 对话框（见图2）： 　 　 　 　 　 　 Grid line Geometry选项X、Y分别表示区域底图的最大外边界经纬度值。本例X选115.54、116.82；Y选32.86、34.08，在Gridding Methodtfad选项框中，选Kriging（插值方法），设置好后，点击“确定”按钮，提示已生成grd格式的文件bzyl.grd，再按“确定”按钮，随后出现Surfer-Report1对话框，关掉对话框不保存文件即可。 2.2.3　Blank文件 点击菜单项Grid、Blank，打开bzyl.grd，再选中区域无县界闭合图地图文件Bz1.bln，生成out.grd文件。 2.2.4　生成等雨量线图 从菜单项选Map、contour? M