Anusplin使用方法(以降水数据为例).pdfVIP

Anusplin 使用方法 一、解压安装包 可见anusplin 核心模块如下： 我们主要用到的是splina 和lapgrd 程序，具体模块介绍和流程图可查询用户手册 二、准备数据格式 2．1 气象台站数据 原始数据格式是以站点名、日期、气温、降水等变量为列，不便插值处理。利用spss 转 化为气象变量固定的文件，每个气象文件中均以站点和日期为列，变量唯一。 原始数据 整理数据 注：数据不含表头，从第 一行起就必须与下面的 splina cmd 文件中的 format 吻合 操作步骤： ① Spss 打开txt 数据文件→菜单中点击“数据”→重组→将个案重组为变量，其中站点名 作为标识变量，年月日作为索引变量，点击下一步直至完成，即可将数据重组。 ② Spss 打开txt 站点信息文件（包含站点经纬度和高程）→在步骤①中的数据文件界面， 菜单中点击“数据”→合并文件、添加变量→选择刚刚打开的站点信息文件→勾选匹配 关键变量的个案，将站点名作为关键变量，活动集为基于变量的表，确定。即可完成按 照站点名进行站点信息的匹配和添加。 ③ 将已包含了经纬高程信息的站点数据另存为固定的ASCII 格式.dat,注意不是制表符分 隔的.dat。可供splina 读取。 2.2 DEM 数据 下载 STRM90m 分辨率的原始高程数据(未经投影) ，在ArcGIS 中对DEM 文件进行镶 嵌、投影（根据研究情况与其他数据统一投影即可）、重采样，最终利用toolbox 中Conversion 工具→From Raster to ASCII ，转化为txt 文本文件，可供lapgrd 读取。 三、编写命令行 3.1 Splina cmd 在写字板中编辑，另存为.cmd 后缀即可，以下是文件和说明： 输出表面名 要插数据的单位（由于是降水数据已处理为m，选择1） 独立变量个数3 （经度、纬度、高程） 协变量个数0 表面样条变量个数 表面协变量个数 经度的左右极限、是否转换（0-不转换）、单位 （1-m ） 纬度的左右极限、是否转换（0-不转换）、单位 （1-m ） 高程的上下限、是否转换 （1-转换方式为x/a ）、单位 （1-m） 转换系数a= 1000 独立变量转换 （2-将所有表面值调整为原来的平方根） 样条次数（2 次样条） 输出表面个数，365 个 0- 为每个表面采用同样权重 1- 优化参数，通常为1 平滑参数选择方法（1-GCV 法）