ANUSPLIN详细使用说明.docx

《ANUSPLIN详细使用说明.docx》由会员分享，可在线阅读，更多相关《ANUSPLIN详细使用说明.docx（6页珍藏版）》请在优知文库上搜索。

1、ANUSP1.lN运用说明anusplin是供应好用的转换分析和对多变量数据采纳薄盘光滑样条插值进行插值的工具。它供应了完整的统计分析、数据诊断以及空间分布标准误。同样也支持多种数据输入和表面查询功能。薄盘光滑样条表面拟合法最早由Wahba出(1979),然后由Bates和Wahba(1982),EkIen(1984),HUIChinSon(1984)和deHoog(1985)对其进行了改进以适用于大数据集。Batesetal.(1987)将其扩展为局部样条法，这样就可以把参数线性亚模型(或协变量)添加插值中，而不像以前只能考虑独立样条变量(即自变量)。这为这些因变量供应由这些自因变量确定的参

2、数form供应了一种极好的方式。在没有独立样条变量的状况下(当前不允许)，将进行简洁的多变量线性回来O薄盘光滑样条插值事实上可以被看作广义的标准多变量线性回来，但是参数模型由一个适用的光滑非参数函数所替代。拟合函数的光滑度，或者与之相对的困难度通常会依据数据拟合表面的最小预料误差通过GCV(广义交叉验证)自动计算。GCV数据模拟结果的验证由Craven和Wahba(1979)年提出并完成。Wahba(1990)对薄盘光滑样条技术各种模块进行了全面介绍。HutChinSon(199Ia)对月平均气象要素空间插值的基本理论和应用进行T简要的概述，Huichinson(1993)和HUtChinSO

3、nGeSSler(1994)对运算和相关行了对比。Hutchinson(1995,1998ab)对该方法运用到了年和日降水数据的插值上。同时模拟多个表面特别便利，尤其是对气象数据而言。ANUSP1.lN在允许随意多的这样的表面和“表面独立变量”，以便独立变量在各表面间做系统地变更。anusplin允许系统查询这些表面及其标准误，点文件或者grid文件都行。anusplin同样允许对独立变量和依从变量进行转换。华中农业商校国家2“工程商校华中农业高校国家211工程高校下边将对ANUSP1.lN个组件进行简要概述。SP1.INA：适用于2000个要素的点文件。允许一个或多个独立变量的随意数量的（局

4、部）薄盘光滑样条函数。光滑系数通常由GCV确定SP1.INB:与SP1.INA大致相同,但是它允许10000个要素的点文件，可以用selnot/delnot模块先从大数据集中设置2000个knots,同时还可用addnot对knots进行编辑。AVGCVA:依据一系列不同的光滑系数，计算出由SP1.INA生成的各表面与光滑系数对应的GCV和平均GCV值，计算结果将写成一个文件为进一步的制图服务AVGCVB:为SP1.lNB服务，功能与AVGCVA类似。1.APPNT：进行局部薄盘光滑样条表面估值并计算贝叶斯标准误,并生成一个点文件。1.APGRD：功能与1.APPNT类似，只是最终结果将是一个

5、规则的矩形栅格文件。这个程序般好通过一个Command-IineShell（即cmd文件）提取.。在Unix系统，运用终端仿真窗口显示一个Shellprompto在微软的WindOWS系统中，启动一个MS-DOSShe11,如输入:Splinajob.log其中，jobd是一个输入的吩咐文件，而job.log则是程序执行后将输出的日志文件。程序需录入项：(programinputs)包括独立样条变量和协变量的个数，每个独立变量的上下限，每个变量的可选的转换试方式以及最低样条次数(确定不同的函数式)，要模拟的表面数，以及各表面用以计算光滑系数的方法。输入输出文件的指定，在SP1.INB中还包括k

6、nois索引文件名。位于独立变量界限外的点数据将予以解除。这个界限有利于在不用切割点文件的状况下指定模拟的表面的范囹。界限值(相当于mask)允许留出肯定的边更便于各表面间有重我部分，对于极大的数据集来说可能须要这样做。用户指定的界限将会对指定的数据格式和独立变量的样条次数进行了简洁的检查。假如选择的点数过少会出现错误提示。通过与标准FOnran9()的结合，ANUSP1.IN4.1可以为更多的数据和数组动态安排内存。同时，SP1.lNA和SPuNB都适用以随意点数来模拟随意数量的表面。然而，最好还是不要让SP1.INA的点数2000,SP1.INB的点数1000().SP1.IN适用于设置的

7、knot于2000的状况。SP1.INA的主要存储空间是点文件所含的点数的平方，处理时间则是点数量的立方。SP1.lNB则与knot须要留意的是，SP1.INA目前(4.1版式本)并不允许点文件中出现坐标重叠的点。SP1.lNB中可以出现重坐标的点，但不允许有重坐标的knot。计算函数略.大致为:Zi=f（xi）+bli+ei（i=其中Zi为预估值,xi是一个d维的样条独立变量矩阵,f（）是Xi的未知光滑函数,yi是一个p维的协变量矩阵,b是yi的未知P维系数矩阵,ei是各独立变量,均值为0（其中涉及到多个厌烦的参数,不想写了,）上式中T为b的上标1）当P=O时,即没有协变量时,该函数式就是一

8、个一般的薄盘样条模型,即简洁的多变量线性回来.而假如f（xi）为0,即没有独立样条变量时,anusplin不行用2）:光滑参数P（不是P,应当是密度那个符号）是一个正数.它越趋近于0,则拟合函数越接近精确值,反之,若越趋向无穷,则拟合函数越接近最小二乘多项式（理解得有点怪）,其阶数由d和Xi的阶数确定.光滑参数通常依据GCV计算拟合表面的最小预料误差确定.（即可自动计算）依从变量（因变量）转换anusplin目前允许对自变量进行开平方或者求自然对数转换.在用Splinsplinb进行拟合前,可选用其中随意一种方式对变量进行转换.变量转换可以使其更符合正态分布,前提是变量的值是正数或非负数.该转

9、换会自动传达给拟合表面的系数，因而1.APPNT和1.APGRD可以选择性的输出转换后的表面值或不经转换的表面值.不经转换的表面估值可用依从变量的逆转换（求平方值或者指数）得到.Hutchinson（l998a）指出,在用薄盘光滑样条法进行拟合前,对依从变量日降水数据进行开平方转换,可削减近10%的插值误差.假如选择输出未经转换的表面估值.1.APPNT和1.APGRD将同时计算其标准差.开平方转换公式的适用性由HUtChinson(1998a)进行了求证.(以下是一大堆公式,略去)拟合气候表面最早的拟合程序通常须要至少两个独立样条变量,(即f(xi)中i维矩阵),通常是经度、纬度，以度为单位

10、。但是在拟合温度或降水量时，可增加第三个独立变量，即海平面以上的高程(海拔)。在把海拔作为第三个独立样条变量时，可以把它的单位转换为千米，即海拔值除以1000。稍梢把海拔值缩放一下(此处指转换单位由米到千米)可略微改进拟合结果。该改动最初由HUtChinSonBiSChOf(1983)提出,并由Hutchinson(1995,1998b)对此进行了论证.在某些区i(Festrictedareas),把海拔作为协变量而非独立样条变量时,拟合效果更佳。在拟合一个温度表面时,作为协变量的海拔的系数将被作为一个先验的温度递减率(HUlChinSOn.199Ia)看待。其它影响气象变量的因素,如地形效应

11、(IOPOgraPhiCeffeCtS)而非海拔,可以作为额外的协变量引入，假如可以求算出合适的参数并且相关数据可用的前提下.Hutchinsonetal.(1984ab.1996a)和Hutchinson(l989am199lab)对其它的气象插值软件进行了介绍.HutchinsoneiaM1992.1996b)介绍了以样条法拟合气象表面应用于农业气象分类以及生物多样性评估的其它软件.在拟合多变量气象表面时,只须要知道样点处的独立变量的值就可以了.因此,气象站点的坐标和海拔信息必需精确.坐标错误的点会在输出的最大残差日志里反映出来,即以降序排列的残差文件中,排在首位的几个极大的残差值对应的站点.最近又出了可检验引入由不同水平辨别率的DEM求算海拔,坡度和坡向进行插值的效用的软件(Hutchinson1995,1998b).1.APGRD可用来估算气象要素和标准误,并生成一个规则的栅格图.rovidedregulargrideachindeendentvariable,additionallatitude,issulied也就是说通常一个规则/标准的DEM栅格数据(是说DEM数据必需是一个规则的栅格？并且以什么投影呢？)是必需的.14111心由5011(1988,1989人1996,1997)介绍了DEM的获得方法