Anusplin使用方法（以降水数据为例）

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

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 dataerror;但此时只是读入了站点的一半;然后删除所有生成文件,将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个表面时都可行。 ……

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说教育文库Anusplin使用方法（以降水数据为例）在线全文阅读。