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

《专业气候数据空间插值软件Anusplin简介.ppt》由会员分享，可在线阅读，更多相关《专业气候数据空间插值软件Anusplin简介.ppt（31页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、专业气候数据空间插值软专业气候数据空间插值软件件AnusplinAnusplin简介简介主要内容主要内容一、基本原理二、模块及流程三、数据格式四、模型选择五、Anusplin应用示例一、基本原理一、基本原理Anusplin基于普通薄盘普通薄盘和局部薄盘样条函数局部薄盘样条函数插值理论。除了可以引入自变量自变量外，还允许引入协变协变量量。eg:温度和海拔的关系，降水和海岸线的关系等。1 1、插值方法：、插值方法：2、局部薄盘光滑样条的理论统计模型公式：、局部薄盘光滑样条的理论统计模型公式：位于空间：位于空间i i点的因变量点的因变量：是要估算关于：是要估算关于x xi i的未知光滑函数，的未知光

2、滑函数，x xi i是独立变量。是独立变量。：为：为p p维独立协变量，维独立协变量，b b为为y yi i的的p p维系数维系数：随机误差：随机误差一、基本原理一、基本原理上式中：函数f和系数b通过最小二乘估计最小二乘估计来确定：其中Jm(f)是函数f(xi)的粗糙度测度函数，定义为函数f的m阶偏导（称为样条次数，也叫粗糙次数）。是正的光滑参数，在数据保真度与曲面的粗糙度之间起平衡作用，在Anusplin中通常用广义交叉验证GCV的最小化以及最大似然法GML的最小化来确定二、模块及流程二、模块及流程1.Anusplin模块：共8个，分别为：SPLINASPLINBSELNOTADDNOTDE

3、LNOTGCVGMLLAPPNTLAPGRD适用于任意个独立变量或多个协变量的薄盘样条函数，站点数小于2000，数据平滑度由GCV或GML决定SPLINASPLINB与SPLINA功能类似，站点数最多可达10000个站点。SELNOT为SPLINB选择初始节点。ADDNOT添加数据节点。说明：说明：删除数据节点DELNOTGCVGML对拟合表面计算GCV或GML误差，用于数据检查或定位。LAPPNT计算预测值或贝叶斯标准误差估计的点文件。LAPGRD生成拟合曲面或贝叶斯标准误差曲面。说明：说明：可生成列表文件列表文件(.opt)残差文件残差文件最优参数文件最优参数文件表面系数文件表面系数文件协

4、变量误差信息协变量误差信息(.res):记录估计与实值的误差(.sur)：记录拟合结果的系数(.lis)(.cov)2.Anusplin流程二、模块及流程二、模块及流程(1)执行SPLINA或SPLINB命令另外，最后还可生成log日志文件，给出统计结果，并可以检查运行时产生的错误。表面系数文件表面系数文件协变量误差信息协变量误差信息由第一步生成的插值表面插值表面标准误差表面标准误差表面得到得到2.Anusplin流程二、模块及流程二、模块及流程(2)执行LAPGRD命令三、三、Anusplin数据格式数据格式Anusplin需要将以固定形式存储的气象要素数据以固定形式存储的气象要素数据进行处

5、理成程序要求的标准格式，生成文本文件。最后存储为最后存储为.dat文件。文件。（特别注意：如果文件格式不对，在程序执行过程中会产生错误。）包括：经度、纬度、高程、气象要素值经度、纬度、高程、气象要素值。在此过程中要对标准文件进行检查修正处理，剔除重复站点，缺侧值用插值法补齐。注：以上经纬度数据为ALBERS投影，单位为米Anusplin在日志文件中提供了一系列用于判断误差来源和插值质量的统计参数。包括：四、四、Anusplin模型选择模型选择有效数量估计Signal(信号自由度)、剩余自由度Error、光滑参数RHO、GCV、期望真实均方误差MSE、最大似然法误差GML、均方残差MSR、方差估

6、计VAR及其平方根及其平方根。由这些来选择最佳模型。Signal：指示了拟合曲面的复杂程度：指示了拟合曲面的复杂程度RHO：平衡了拟合曲面的精确度与平滑度：平衡了拟合曲面的精确度与平滑度RHO过小和过小和Signal大于观测站点的一半大于观测站点的一半或者RHO过大过大都预示着拟合过程找不到最优光滑参数，可能数据点过于稀疏或者拟合函数过于复杂，不适用插值。这在Anusplin中会以*符号标出。四、四、Anusplin模型选择模型选择因此：最佳模型判断标准：（1）GCV或或GML最小、最小、（2）信噪比）信噪比SNR（信号自由度与剩余自由（信号自由度与剩余自由度之比）最小、度之比）最小、（3）信

7、号自由度小于站点的一半、）信号自由度小于站点的一半、（4）模型成功率判断中无）模型成功率判断中无*表示。表示。四、四、Anusplin模型选择模型选择（1）数据准备：将数据准备成如下格式（本例中在spss中整理，也可在其他软件中完成）五、五、Anusplin应用示例应用示例说明：CODE为站点的label(其要求是字符型,长度小于20即可,图中长度为6,在fortan语言中可用a6表示)X,Y为投影坐标系以米为单位的坐标值,浮点型,在fortan 语言中用f表示,即f14.6,说明是带6位小数总长度为14.ELEV:海拔.单位为米。Temp代表气象变量温度，本图中只产生一个面（即一个月）也可以

8、同时进行多个月份或多年的数据，只需增加即可。把上图中的类型,宽度,小数三项设置好后,将该数据另存为固定ASCii格式,后缀为了.dat.这样Anusplin的所需数据准备完毕。五、五、Anusplin应用示例应用示例(2)安装Anusplina,将SPLINA.exe和LAPGRD.exe文件（常用的两个命令）拷入运行的数据路径下。五、五、Anusplin应用示例应用示例（3）编写运行SPLINA.exe所需要的cmd文件。五、五、Anusplin应用示例应用示例新建一个文件夹，键入右图显示内容，另存为.cmd格式。说明：五、五、Anusplin应用示例应用示例说明：五、五、Anusplin应

9、用示例应用示例（4）编写运行LAPGRD.exe所需要的cmd文件。新建一个文件夹，键入右图显示内容，另存为.cmd格式。五、五、Anusplin应用示例应用示例说明：五、五、Anusplin应用示例应用示例开始菜单中,在运行项中输入cmd 回车,把cmd运行程序打开（5）执行命令五、五、Anusplin应用示例应用示例用cd命令进入anusplin所在的文件夹，注：anusplin的安装文件和.dat和.cmd文件都要放在同一个文件夹中,路径必须一致.例:如果放在d:temp文件夹下,则先输入d:回车,把当前盘符由C盘转到d:然后输入cd d:temp命令,以进入该文件夹.五、五、Anusp

10、lin应用示例应用示例输入splina198201tmean.log 回车.得到一个198201tmean22.log的文件,记录了该cmd的拟合信息。输入lapgrd198201tmeanrd.log 回车.得到一个198201tmean.grd,为插值结果。五、五、Anusplin应用示例应用示例五、五、Anusplin应用示例应用示例（6）插值模型选择总共有spline模型（独立变量、协变量和样条次数多种组合共18种）比照log中提供的信息进行筛选.根据之前所说的判断标准.即无*标示.Signal小于样点数的一半,Signal/Error最小,RTGCV,RTMSR,RTMSE最小,其中

11、RTGCV和RTMSE的差值越大,可间接说明模型的解释率越高。本示例中以经纬度作为自变量，高程作为协变量，样条次数为3为最优方案。右图为1982年1月月平均温度插值结果。插值结果：参考文献钱永兰,吕厚荃,张艳红.基于ANUSPLIN软件的逐日气象要素插值方法应用与评估J.气象与环境学报,2010,26(2):7-15.刘志红,Tim R.McVicar,Li LingTao等.基于ANUSPLIN的时间序列气象要素空间插值J.西北农林科技大学学报：自然科学版.2008,36(10):227-234.刘志红,Li Lingtao,Tim R.McVicar等.专用气候数据空间插值软件ANUSPLIN及其应用J.气象,2008,34(2):92-100.百度“美丽心情”博客：http:/ 2 3 4 5 6 7 8 9 10 11 12。这样才能出来12个面。仓促做出的总结非常粗糙，我会在接下来的时间里继续修改增加内容！To be continued.结束语结束语谢谢大家聆听！谢谢大家聆听！31