应用Win-Surfer软件绘制降水等值线图-生态环境学报.doc

PAGE 362 土壤与环境 第11卷第4期（2002年11月） 文雅等：应用Win-Surfer软件绘制降水等值线图 PAGE 361  土壤与环境 2002, 11(4): 360~362 Soil and Environmental Sciences E-mail: ses@ 基金项目：广东省“百项工程”项目（99B05801G） 作者简介：文 雅（1969－），女，助理研究员，主要研究方向为土地资源与土壤信息系统。E-mail: wenya@ 收稿日期：2002-07-02 应用Win-Surfer软件绘制降水等值线图 文 雅，郭治兴 广东省生态环境与土壤研究所，广东省农业环境综合治理重点实验室，广东 广州 510650 摘要：空间插值是各种地学、生态学研究中经常遇到的问题。文章以中国1950年至1998年160个测站年均降水量数据为例，介绍采用空间插值绘图软件Win-Surfer绘制等值线图的基本方法与步骤，并对软件中几种常用的空间插值方法的优缺点进行简要的对比分析。 关键词：空间插值；等值线；Win-surfer软件 中图分类号：P426.61+3 文献标识码：A 文章编号：1008-181X（2002）04-0360-03 在进行各种生态环境问题的研究过程中，常常会涉及到绘制等值线图的问题，如绘制降水等值线、土壤养分空间分布、生物量区域变化等，其目的是将离散的台站数据或野外观测数据进行空间插值处理，进而绘出等值线，为下一步的研究提供直观的图形支持。过去这样的工作常常由手工完成，费时费力，随着计算机软硬件水平的不断提高，已经能够很方便地借助计算机完成这项工作。这里以中国1950年至1998年160个测站年均降水量为例，介绍一个小型实用的空间插值绘图软件Win-Surfer，应用它能很方便地完成样点的空间插值，并绘出等值线图。 1 原始数据的格式及输入 Win-Surfer自己带有打开和编辑数据文件的工作表格（work-sheet），类似Microsoft Execl中的工作表，可以通过它输入和编辑数据文件，但其功能有限，对于文本格式的数据文件有时会读成错误的排列方式或出现乱码，通过Microsoft Execl转换成扩展名为xls或ASCII格式扩展名为dat的数据文件后，Surfer都能正常打开。本例中数据文件名为chnrain.dat，它包含的数据项有样点的位置坐标（经纬度）及属性项（年均降水量），属性项可以有多项，但样点的位置坐标是唯一的。下面我们就可以对原始数据进行空间插值处理。 2 离散数据的空间插值 对离散的原始样点数据进行空间插值是绘制等值线的前提条件。选择Surfer主菜单Grid中Data项，出现选项对话框，选定chnrain.dat文件后，出现离散数据插值对话框，在这一对话框中有几个选项区，在数据栏（data columns）中将X选为column B（经度值），Y选为column C（纬度值），Z选为column D（降水量），选好后可按data info可以了解数据的基本统计特征，包括最大值、最小值、均值、方差等；然后选择插值处理方法，虽然在应用Surfer进行空间插值时不需要知道插值方法的具体内容，但了解几种常用的空间插值方法的原理及其优缺点，对于我们针对不同空间数据选择适当的插值方法及其相关参数是十分重要的[1, 2]。 2.1 距离倒数加权法（Inverse Distance Weighted Method） 该方法以数据点与网格点的距离平方倒数作为权值，距离越大权值越小，反之则越大。在具体应用时，既可以把它当作精确估值，又可当作平滑估值。 距离倒数加权法其优点在于算法简单，易于实现。在进行插值时其结果出在用于插值数据的最大和最小值之间。缺点是：距离加权法没有考虑数据场在空间的分布，往往会因为采样点的分布不均而使得估值结果产生偏差。另外，由于插值结果肯定介于估值点的最大和最小值之间，因此会因为采样点的不均而对空间场产生错误的估计。 2.2 最小曲率法（Minimum Curvature） 该法构造出具有最小曲率的曲面，使其穿过空间场的每一点，并尽可能使曲面变得光滑。最小曲率法是平滑