专业气候数据空间插值软件Anusplin简介

专业气候数据空间插值软件Anusplin简介QQ：2893683302013.1.29主要内容•一、基本原理•二、模块及流程•三、数据格式•四、模型选择•五、Anusplin应用示例一、基本原理Anusplin基于普通薄盘和局部薄盘样条函数插值理论。除了可以引入自变量外，还允许引入协变量。eg:温度和海拔的关系，降水和海岸线的关系等。1、插值方法：2、局部薄盘光滑样条的理论统计模型公式：iZ)(ixfieiiTiieybxfz)(）（Ni,...,1：位于空间i点的因变量：是要估算关于xi的未知光滑函数，xi是独立变量。iy：为p维独立协变量，b为yi的p维系数：随机误差一、基本原理上式中：函数f和系数b通过最小二乘估计来确定：)()(21fJwybxfzmNiiiTii其中Jm(f)是函数f(xi)的粗糙度测度函数，定义为函数f的m阶偏导（称为样条次数，也叫粗糙次数）。ρ是正的光滑参数，在数据保真度与曲面的粗糙度之间起平衡作用，在Anusplin中通常用广义交叉验证GCV的最小化以及最大似然法GML的最小化来确定二、模块及流程•1.Anusplin模块：共8个，分别为：SPLINASPLINBSELNOTADDNOTDELNOTGCVGMLLAPPNTLAPGRD适用于任意个独立变量或多个协变量的薄盘样条函数，站点数小于2000，数据平滑度由GCV或GML决定SPLINASPLINB与SPLINA功能类似，站点数最多可达10000个站点。SELNOT为SPLINB选择初始节点。ADDNOT添加数据节点。说明：删除数据节点DELNOTGCVGML对拟合表面计算GCV或GML误差，用于数据检查或定位。LAPPNT计算预测值或贝叶斯标准误差估计的点文件。LAPGRD生成拟合曲面或贝叶斯标准误差曲面。说明：可生成列表文件(.opt)残差文件最优参数文件表面系数文件协变量误差信息(.res):记录估计与实值的误差(.sur)：记录拟合结果的系数(.lis)(.cov)•2.Anusplin流程二、模块及流程(1)执行SPLINA或SPLINB命令表面系数文件协变量误差信息由第一步生成的插值表面标准误差表面得到•2.Anusplin流程二、模块及流程(2)执行LAPGRD命令三、Anusplin数据格式Anusplin需要将以固定形式存储的气象要素数据进行处理成程序要求的标准格式，生成文本文件。最后存储为.dat文件。（特别注意：如果文件格式不对，在程序执行过程中会产生错误。）包括：经度、纬度、高程、气象要素值。在此过程中要对标准文件进行检查修正处理，剔除重复站点，缺侧值用插值法补齐。注：以上经纬度数据为ALBERS投影，单位为米•Anusplin在日志文件中提供了一系列用于判断误差来源和插值质量的统计参数。包括：四、Anusplin模型选择•有效数量估计Signal(信号自由度)、剩余自由度Error、光滑参数RHO、GCV、期望真实均方误差MSE、最大似然法误差GML、均方残差MSR、方差估计VAR及其平方根。由这些来选择最佳模型。Signal：指示了拟合曲面的复杂程度RHO：平衡了拟合曲面的精确度与平滑度RHO过小和Signal大于观测站点的一半或者RHO过大都预示着拟合过程找不到最优光滑参数，可能数据点过于稀疏或者拟合函数过于复杂，不适用插值。这在Anusplin中会以*符号标出。四、Anusplin模型选择因此：最佳模型判断标准：•（1）GCV或GML最小、•（2）信噪比SNR（信号自由度与剩余自由度之比）最小、•（3）信号自由度小于站点的一半、•（4）模型成功率判断中无*表示。四、Anusplin模型选择•（1）数据准备：将数据准备成如下格式（本例中在spss中整理，也可在其他软件中完成）五、Anusplin应用示例说明：•CODE为站点的label(其要求是字符型,长度小于20即可,图中长度为6,在fortan语言中可用a6表示)X,Y为投影坐标系以米为单位的坐标值,浮点型,在fortan语言中用f表示,即f14.6,说明是带6位小数总长度为14.•ELEV:海拔.单位为米。•Temp代表气象变量温度，本图中只产生一个面（即一个月）也可以同时进行多个月份或多年的数据，只需增加即可。把上图中的类型,宽度,小数三项设置好后,将该数据另存为固定ASCii格式,后缀为了.dat.这样Anusplin的所需数据准备完毕。五、Anusplin应用示例•(2)安装Anusplina,将SPLINA.exe和LAPGRD.exe文件（常用的两个命令）拷入运行的数据路径下。五、Anusplin应用示例•（3）编写运行SPLINA.exe所需要的cmd文件。五、Anusplin应用示例•新建一个文件夹，键入右图显示内容，另存为.cmd格式。•说明：五、Anusplin应用示例•说明：五、Anusplin应用示例•（4）编写运行LAPGRD.exe所需要的cmd文件。•新建一个文件夹，键入右图显示内容，另存为.cmd格式。五、Anusplin应用示例•说明：五、Anusplin应用示例•①开始菜单中,在运行项中输入cmd回车,把cmd运行程序打开•（5）执行命令五、Anusplin应用示例•②用cd命令进入anusplin所在的文件夹，注：anusplin的安装文件和.dat和.cmd文件都要放在同一个文件夹中,路径必须一致.•例:如果放在d:\temp文件夹下,则先输入d:回车,把当前盘符由C盘转到d:•然后输入cdd:\temp命令,以进入该文件夹.五、Anusplin应用示例③输入splina198201tmean.cmd198201tmean.log回车.得到一个198201tmean22.log的文件,记录了该cmd的拟合信息。④输入lapgrd198201tmeangrd.cmd198201tmeanrd.log回车.得到一个198201tmean.grd,为插值结果。五、Anusplin应用示例五、Anusplin应用示例•（6）插值模型选择•总共有spline模型（独立变量、协变量和样条次数多种组合共18种）比照log中提供的信息进行筛选.根据之前所说的判断标准.即无*标示.Signal小于样点数的一半,Signal/Error最小,RTGCV,RTMSR,RTMSE最小,其中RTGCV和RTMSE的差值越大,可间接说明模型的解释率越高。•本示例中以经纬度作为自变量，高程作为协变量，样条次数为3为最优方案。•右图为1982年1月月平均温度插值结果。插值结果：参考文献•钱永兰,吕厚荃,张艳红.基于ANUSPLIN软件的逐日气象要素插值方法应用与评估[J].气象与环境学报,2010,26(2):7-15.•刘志红,TimR.McVicar,LiLingTao等.基于ANUSPLIN的时间序列气象要素空间插值[J].西北农林科技大学学报：自然科学版.2008,36(10):227-234.•刘志红,LiLingtao,TimR.McVicar等.专用气候数据空间插值软件ANUSPLIN及其应用[J].气象,2008,34(2):92-100.•百度“美丽心情”博客：最后一些注意事项：•基本上出错的点都在编写splina.exe和lapgrd.exe的cmd文件上，总结了几点易错处：•（1）编写.dat文件时格式一定要和cmd中编写的对应。•（2）splina的cmd文件中，编写范围时一定注意：记住X的最大值=最小值+cellsize*rows,Y最大值=最小值+cellsize*colums,不满足该公式就会出错。•（3）lapgrd的cmd文件中，编写范围时不能再照搬splina中的范围，要和.sur文件中的范围一致，因为.sur文件对原先范围进行了四舍五入，如果还按照以前的写就会出错。•（4）在编写lapgrd的cmd文件时，若要一次性出来12个面，切不能只写12，而是要写成123456789101112。这样才能出来12个面。•仓促做出的总结非常粗糙，我会在接下来的时间里继续修改增加内容！•Tobecontinued...