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

Anusplin 使用方法
一、 解压安装包
可见anusplin 核心模块如下：
我们主要用到的是splina 和lapgrd 程序，具体模块介绍和流程图可查询用户手册
二、准备数据格式
2．1 气象台站数据
原始数据格式是以站点名、日期、气温、降水等变量为列，不便插值处理。

利用spss 转化为气象变量固定的文件，每个气象文件中均以站点和日期为列，变量唯一。

操作步骤：
①
Spss 打开txt 数据文件→菜单中点击“数据”→重组→将个案重组为变量，其中站点名作为标识变量，年月日作为索引变量，点击下一步直至完成，即可将数据重组。

原始数据
整理数据
注：数据不含表头，从第
一行起就必须与下面的
splina cmd 文件中的
format 吻合
② 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法）
输入数据文件.dat
站点个数，设置得大一些
站点标识符个数，5个字符
输入数据格式，与spss中数据宽度和小数位一致。

注意：开始时须在365f5.1后加/，并在原始数据文件
的高程值后面加回车，才能读入，否则提示input data
error；但此时只是读入了站点的一半；然后删除所有
生成文件，将cmd文件中的/去掉再次运行，方可全
部读入。

……并不知道原理是什么，试验了好久才发
现这个规律
5个输出文件名的定义
3.2 lapgrd cmd
…… ……输入表面文件的名称
拟合表面个数（0-输出所有表面，如果输出前三个需写123） 计算表面类型（1-拟合表面值）
表面值转换（1-应用表面值转换，即之前的平方根） 输入协方差文件的名称
2-计算预测标准误差
空行，最大标准误差不填
1用单元中心值
变量索引，第一个变量为经度，写1
左右范围与分辨率（小于splina 生成的sur 范围） 变量索引，第二个变量为经度，写2
上下范围与分辨率（小于splina 生成的sur 范围） 0- 无掩模文件
2- Arcinfo 独立可变网格格式
输入dem 文件
输出栅格文件格式2-arcinfo
Dem 中无效值定义
输出grd 命名
空行-推荐默认输出二进制格式，占内存少速度快。

如设置成(365f5.1)批量处理后在arcgis 中打开无响应（可能只是个人的个别情况），只输出1个表面时都可行。

输出表面误差栅格文件的格式 2-arcinfo
误差栅格特殊值
误差栅格命名
空行-推荐默认输出二进制格式，占内存少速度快。

如设置成(365f5.1)批量处理后在arcgis 中打开无响应（可能只是个人的个别情况），只输出1个表面时都可行。

……。