書式設定,書式設定,第 2,第 3,第 4,第 5,*,書式設定,書式設定,第 2,第 3,第 4,第 5,*,专业气候数据空间插值软件Anusplin简介,主要内容,一、根本原理,二、模块及流程,三、数据格式,四、模型选择,五、Anusplin应用例如,一、根本原理,Anusplin基于,普通薄盘,和,局部薄盘样条函数,插值理论。,除了可以引入,自变量,外，还允许引入,协变量,。,eg:温度和海拔的关系，降水和海岸线的关系等。,1、插值方法：,2、局部薄盘光滑样条的理论统计模型公式：,：位于空间i点的因变量,：是要估算关于x,i,的未知光滑函数，x,i,是独立变量。,：为p维独立协变量，b为y,i,的p维系数,：随机误差,一、根本原理,上式中：函数f和系数b通过,最小二乘估计,来确定：,其中Jm(f)是函数f(xi)的粗糙度测度函数，定义为函数f的m阶偏导称为样条次数，也叫粗糙次数。,是正的光滑参数，在数据保真度与曲面的粗糙度之间起平衡作用，在Anusplin中通常用广义交叉验证GCV的最小化以及最大似然法GML的最小化来确定,二、模块及流程,1.Anusplin模块：共8个，分别为：,SPLINA,SPLINB,SELNOT,ADDNOT,DELNOT,GCVGML,LAPPNT,LAPGRD,适用于任意个独立变量或多个协变量的薄盘样条函数，站点数小于2000，数据平滑度由GCV或GML决定,SPLINA,SPLINB,与SPLINA功能类似，站点数最多可达10000个站点。,SELNOT,为SPLINB选择初始节点。,ADDNOT,添加数据节点。,说明：,删除数据节点,DELNOT,GCVGML,对拟合外表计算GCV或GML误差，用于数据检查或定位。,LAPPNT,计算预测值或贝叶斯标准误差估计的点文件。,LAPGRD,生成拟合曲面或贝叶斯标准误差曲面。,说明：,可生成,列表文件,(.opt),残差文件,最优参数文件,外表系数文件,协变量误差信息,(.res):记录估计与实值的误差,(.sur)：记录拟合结果的系数,(.lis),(.cov),2.Anusplin流程,二、模块及流程,(1)执行SPLINA或SPLINB命令,另外，最后还可生成log日志文件，给出统计结果，并可以检查运行时产生的错误。,外表系数文件,协变量误差信息,由第一步生成的,插值外表,标准误差外表,得到,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模型选择,因此：最正确模型判断标准：,1GCV或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应用例如,新建一个文件夹，键入右图显示内容，另存为 d格式。,说明：,五、Anusplin应用例如,说明：,五、Anusplin应用例如,4编写运行LAPGRD.exe所需要的cmd文件。,新建一个文件夹，键入右图显示内容，另存为 d格式。,五、Anusplin应用例如,说明：,五、Anusplin应用例如,开始菜单中,在运行项中输入cmd 回车,把cmd运行程序翻开,5执行命令,五、Anusplin应用例如,用cd命令进入anusplin所在的文件夹，注：anusplin的安装文件和.dat和 d文件都要放在同一个文件夹中,路径必须一致.,例:如果放在d:temp文件夹下,那么先输入d:回车,把当前盘符由C盘转到d:,然后输入cd d:temp命令,以进入该文件夹.,五、Anusplin应用例如,输入splina198201tmean.log 回车.得到一个198201tmean22.log的文件,记录了该cmd的拟合信息。,输入lapgrd198201tmeanrd.log 回车.得到一个198201tmean.grd,为插值结果。,五、Anusplin应用例如,五、Anusplin应用例如,6插值模型选择,总共有spline模型独立变量、协变量和样条次数多种组合共18种,比照log中提供的信息进行筛选.,根据之前所说的,判断标准.即无*标示.Signal小于样点数的一半,Signal/Error最小,RTGCV,RTMSR,RTMSE最小,其中RTGCV和RTMSE的差值越大,可间接说明模型的解释率越高,。,本例如中以经纬度作为自变量，高程作为协变量，样条次数为3为最优方案。,右图为1982年1月月平均温度插值结果。,插值结果：,参考文献,钱永兰,吕厚荃,张艳红.基于ANUSPLIN软件的逐日气象要素插值方法应用与评估J.气象与环境学报,2021,26(2):7-15.,刘志红,Tim R.McVicar,Li LingTao等.基于ANUSPLIN的时间序列气象要素空间插值J.西北农林科技大学学报：自然科学版.2021,36(10):227-234.,刘志红,Li Lingtao,Tim R.McVicar等.专用气候数据空间插值软件ANUSPLIN及其应用J.气象,2021,34(2):92-100.,最后一些本卷须知：,根本上出错的点都在编写splina.exe和lapgrd.exe的cmd文件上，总结了几点易错处：,1编写.dat文件时格式一定要和cmd中编写的对应。,2splina的cmd文件中，编写范围时一定注意：记住X的最大值=最小值+cellsize*rows,Y最大值=最小值+cellsize*colums,不满足该公式就会出错。,3 lapgrd的cmd文件中，编写范围时不能再照搬splina中的范围，要和.sur文件中的范围一致，因为.sur文件对原先范围进行了四舍五入，如果还按照以前的写就会出错。,4在编写lapgrd的cmd文件时，假设要一次性出来12个面，切不能只写12，而是要写成1 2 3 4 5 6 7 8 9 10 11 12。这样才能出来12个面。,仓促做出的总结非常粗糙，我会在接下来的时间里继续修改增加内容！,To be continued.,