ANUSPLIN41用户指南_精品文档Word格式.doc

ANUSPLIN41用户指南_精品文档Word格式.doc

《ANUSPLIN41用户指南_精品文档Word格式.doc》由会员分享，可在线阅读，更多相关《ANUSPLIN41用户指南_精品文档Word格式.doc（18页珍藏版）》请在冰豆网上搜索。

GCV对数据模拟结果的验证由Craven和Wahba（1979）年提出并完成。

Wahba（1990）对薄盘光滑样条技术各种模块进行了全面介绍。

Hutchinson（1991a）对月平均气象要素空间插值的基本理论和应用进行了简要的概述，Hutchinson（1993）和HutchinsonGessler（1994）对运算和相关的统计分析进行了更全面的讨论和与Kriging进行了对比。

Hutchinson（1995，1998ab）对该方法运用到了年和日降水数据的插值上。

同时模拟多个表面非常方便，尤其是对气象数据而言。

ANUSPLIN现在允许任意多的这样的表面和“表面独立变量”，以便独立变量在各表面间做系统地改变。

ANUSPLIN允许系统查询这些表面及其标准误，点文件或者grid文件都行。

ANUSPLIN同样允许对独立变量和依从变量进行转换。

下边将对ANUSPLIN的9个组件进行简要概述。

SPLINA：

适用于<

2000个要素的点文件。

允许一个或多个独立变量的任意数量的（局部）薄盘光滑样条函数。

光滑系数通常由GCV决定。

SPLINB：

与SPLINA大致相同，但是它允许<

10000个要素的点文件，可以用selnot/delnot模块先从大数据集中设置<

2000个knots，同时还可用addnot对knots进行编辑。

AVGCVA：

根据一系列不同的光滑系数，计算出由SPLINA生成的各表面与光滑系数对应的GCV和平均GCV值，计算结果将写成一个文件为进一步的制图服务。

AVGCVB：

为SPLINB服务，功能与AVGCVA类似。

LAPPNT：

进行局部薄盘光滑样条表面估值并计算贝叶斯标准误，并生成一个点文件。

LAPGRD：

功能与LAPPNT类似，只是最后结果将是一个规则的矩形栅格文件。

-----------------------------------

这个程序最好通过一个command-lineshell（即cmd文件）提取.。

在Unix系统，使用终端仿真窗口显示一个shellprompt。

在微软的Windows系统中，启动一个MS-DOSshell，如输入：

Splina<

job.cmd>

job.log

其中，job.cmd是一个输入的命令文件，而job.log则是程序执行后将输出的日志文件。

程序需录入项：

（programinputs）

包括独立样条变量和协变量的个数，每个独立变量的上下限，每个变量的可选的转换试方式以及最低样条次数（决定不同的函数式），要模拟的表面数，以及各表面用以计算光滑系数的方法。

输入输出文件的指定，在SPLINB中还包括knots索引文件名。

位于独立变量界限外的点数据将予以排除。

这个界限有利于在不用切割点文件的情况下指定模拟的表面的范围。

界限值（相当于mask）允许留出一定的边距便于各表面间有重叠部分，对于极大的数据集来说可能需要这样做。

用户指定的界限将会对指定的数据格式和独立变量的样条次数进行了简单的检查。

如果选择的点数过少会出现错误提示。

通过与标准Fortran90的结合，ANUSPLIN4。

1可以为更多的数据和数组动态分配内存。

同时，SPLINA和SPLINB都适用以任意点数来模拟任意数量的表面。

然而，最好还是不要让SPLINA的点数>

2000，SPLINB的点数>

10000.SPLIN适用于设置的knot少于2000的情况。

SPLINA的主要存储空间是点文件所含的点数的平方，处理时间则是点数量的立方。

SPLINB则与knot相关。

需要注意的是，SPLINA目前（4。

1版式本）并不允许点文件中出现坐标重叠的点。

SPLINB中可以出现重坐标的点，但不允许有重坐标的knot。

计算函数略.大致为:

Zi=f（xi）+bTyi+ei（i=1....N）

其中Zi为预估值,xi是一个d维的样条独立变量矩阵,f（）是xi的未知光滑函数,yi是一个p维的协变量矩阵,b是yi的未知p维系数矩阵,ei是各独立变量,均值为0,.........（其中涉及到多个讨厌的参数,不想写了,）上式中T为b的上标.

1）:

当P=0时,即没有协变量时,该函数式就是一个普通的薄盘样条模型,即简单的多变量线性回归.而如果f（xi）为0,即没有独立样条变量时,ANUSPLIN不可用.

2）:

光滑参数P（不是P,应该是密度那个符号）是一个正数.它越趋近于0,则拟合函数越接近精确值,反之,若越趋向无穷,则拟合函数越接近最小二乘多项式（理解得有点怪）,其阶数由d和xi及m的阶数决定.光滑参数通常根据GCV计算拟合表面的最小预测误差决定.（即可自动计算）

---------------

依从变量（因变量）转换

ANUSPLIN目前允许对自变量进行开平方或者求自然对数转换.在用SPLIN或SPLINB进行拟合前,可选用其中任意一种方式对变量进行转换.变量转换可以使其更符合正态分布,前提是变量的值是正数或非负数.该转换会自动传达给拟合表面的系数,因而LAPPNT和LAPGRD可以选择性的输出转换后的表面值或不经转换的表面值.不经转换的表面估值可用依从变量的逆转换（求平方值或者指数）得到.Hutchinson（1998a）指出,在用薄盘光滑样条法进行拟合前,对依从变量日降水数据进行开平方转换,可减少近10%的插值误差.

如果选择输出未经转换的表面估值,LAPPNT和LAPGRD将同时计算其标准差.开平方转换公式的适用性由Hutchinson（1998a）进行了求证...（以下是一大堆公式,略去）

拟合气候表面

最早的拟合程序通常需要至少两个独立样条变量,（即f（xi）中i为2维矩阵）,通常是经度、纬度，以度为单位。

但是在拟合温度或降水量时，可增加第三个独立变量，即海平面以上的高程（海拔）。

在把海拔作为第三个独立样条变量时，可以把它的单位转换为千米，即海拔值除以1000。

稍稍把海拔值缩放一下（此处指转换单位由米到千米）可略微改进拟合结果。

该改动最初由Hutchinson和Bischof（1983）提出，并由Hutchinson（1995，1998b）对此进行了论证.

在某些区域（restrictedareas）,把海拔作为协变量而非独立样条变量时,拟合效果更佳。

在拟合一个温度表面时,作为协变量的海拔的系数将被作为一个先验的温度递减率（Hutchinson,1991a）看待。

其它影响气象变量的因素,如地形效应（topographiceffects）而非海拔,可以作为额外的协变量引入,如果可以求算出合适的参数并且相关数据可用的前提下.Hutchinsonetal.（1984ab,1996a）和Hutchinson（1989am1991ab）对其它的气象插值软件进行了介绍.Hutchinsonetal.（1992,1996b）介绍了以样条法拟合气象表面应用于农业气象分类以及生物多样性评估的其它软件.

在拟合多变量气象表面时,只需要知道样点处的独立变量的值就可以了.因此,气象站点的坐标和海拔信息必须准确.坐标错误的点会在输出的最大残差日志里反映出来,即以降序排列的残差文件中,排在首位的几个极大的残差值对应的站点.最近又出了可检验引入由不同水平分辨率的DEM求算海拔,坡度和坡向进行插值的效用的软件（Hutchinson1995,1998b）.

LAPGRD可用来估算气象要素和标准误,并生成一个规则的栅格图.providedaregulargridofvaluesofeachindependentvariable,additionaltolongitudeandlatitude,issupplied.也就是说通常一个规则/标准的DEM栅格数据（是说DEM数据必须是一个规则的栅格?

并且以什么投影呢?

）是必须的.Hutchinson（1988,1989b,1996,1997）介绍了DEM的获取方法.

SPLINA和SPLINB用户指定项：

双击SPLINA.exe,打开一个基于DOS的界面框，根据屏幕的提示操作，解说如下：

titleoffittedsurfaces：

（最大可输入60个字符）：

surfaces文件的头信息，即对要输出的栅格表面命名

surfacevalueunits（输入0-8间的数字）：

拟合变量的单位，如降水为mm，温度为degree，输入对应的数字即可,

numberofindependentsplinevariables：

（输入非负整数）独立样条变量的个数。

不可以大于指定的数量（当前可输入最大数为10）

numberofindependentcovariates：

（输入非负整数）独立协变量的个数。

受样条自变量个数的限制.

numberofsurfaceindependentsplinevariables：

（输入非负整数）表面独立样条变量的个数。

为各拟合表面输入不同的独立变量时需要填写，若所有表面采用的独立变量都相同,则该项输入0.

numberofsurfaceindependentcovariates：

同上,是各表面的各自的协变量个数.

independentvariablelowerandupperlimits,transformationcode,unitscode,optionalmargins:

（输入2个实数,两个非负整（0-8）,各独立变量的两个实数.）独立变量的上下界限,转换代码,单位代码,可选边界（margins）.（上下界相当于mask,输入的范围可与DEM的范围一致.为左上角,右下角的XY读值,有四个数,以空格隔开）先输入下端界限值,再输入上端界限值,位于上下界范围之外的点,虽然在边界内,也会被忽略掉.一到两个，边界可能会被遗漏掉（datapointsoutsidetheselimits,augmentedbymargins,areignored.oneorbothmarginsmaybeomitted.ifonemarginissupplieditisusedasthecommonloweranduppermargin.ifbothmarginsareomittedthetransformationcodeandunitscodemayalsobeomitted.unitscodeasforsurfacevalueunitscode.）[这段里边,由于对margin的理解不通,整段翻译出错]如果两个