Anusplin气象插值总结（气温和降水）

        

• 气温插值

1. 数据准备：

需要气象站点数据，其包括站点名称,经度，纬度，高程，插值变量（降水，气温。注意其单位需要转换，mm，摄氏度），经纬度需转换为十进制。转换：F=C*9/5+32,C=(F-32)*5/9,毫米=英寸*25.4.整理如下：

将上述Excel导入arcgis.

具体步骤：A、数据框设置为常用投影：WGS_84_Albers.

B、导入Excel，file-add xy date-x,y（longitude，latitude）-descriptiol- WGS_84-导出shp-clip出研究区的气象站点

C、打开shp属性表-add field-x，y-calculate geometry-分别选择对应x，y—导出txt

D、将txt导入Excel重新整理成如下并保存为csv格式：

E、打开spss，加载csv文件，本地编码，删除数据视图第一行，改变变量视图如下：

设置好之后将该数据另存为"固定ASCii格式"及本地编码，后缀为了.dat。这时SPSS会输出这五个变量的数据类型，分别为A6，F14.4，F14.4，F8.1，F8.2。数据类型在.cmd文件中将会用到。

1.  Dem准备：

将dem转为ASCii格式：toolbox-conversion tools-from raster-raster to ASCII.保存是TXT文件，留着备用。

   

1.  编写splina .cmd 

TMIN_JJA_82S_MG       #是对该项文件命名,可自定义,但两个cmd文件的命名不可重

5                     #代表t的单位是度

2                     #代表两个自变量(X,Y)

1                     #代表一个样条协变量，即高程(降水不需要)

0                     #代表无表面变量

0                     #无表现协变量

89200 3129300 0 1    #代表每一个自变量的上下限,有无转换,单位是什么.本例中该项为拟合区域的X的最大最小坐标,即最左边界和最右边界的X读值.0代表不转换,1代表单位为米.不要最大值最小值坐标弄成度为单位，要以米为单位。

3963450 6829600 0 1    #与上同.代表的是Y值上下限,即最高最低处值.

(记住X的最大=Xllcorner+cellsize*nocols,Y最大值=Yllcorner +cellsize*nrows,不满足该公式就会出错,其数值在dem转为ASCII的TXT文件中有)

-145 2974 1 1   #代表ELEV的上下限,1 1代表用ELEV进行X/A的转换,即把单位由米转换成千米,后边的1代表原单位是米.

1000            #代表上述的A值等于1000.0

0               #独立变量转换参数,0 为不转换

2               #样条次数（样条次数2较好，3标准差较大）

1               #输出表面个数

0               #相对误差方差个数,0 为每个表面采用统一权重

1               #优化参数指标,通常为1

1                 #平滑参数选择方法1-GCV ,2-MSE ,3-固定值,4-GML

TMIN_JJA_82_M.dat    #需要插值的数据文件，spss输出对应的dat文件

200                  #样本总数,比实际样本稍多

6                      #站标字符数

(A6,F14.4,F14.4,F8.1,F8.2) #数据格式，spss输出结果

TMIN_JJA_82_M.res      #输出残差文件

TMIN_JJA_82_M.opt    #光滑参数文件

TMIN_JJA_82_M.sur   #表面文件,用于下一步LAPGRD的输入数据

TMIN_JJA_82_M.lis   #列表文件,观测数据、预测数据、误差数据

TMIN_JJA_82_M.cov  #拟合表面系数的误差协方差文件,下一步LAPGRD 的输入

    

注意：编写时要空多行。无理由

1. 运行splina  cmd

开始菜单中,在"运行"项中输入cmd 回车,把cmd运行程序打开；用cd命令进入anusplin所在的文件夹.切记,anusplin的安装文件和.dat和.cmd文件都要放在同一个文件夹中,路径必须一致.例:如果放在L:\data\anusplin文件夹下,则先输入L:回车,把当前盘符由C盘转到L:

然后输入cd data\anusplin命令,以进入该文件夹；输入splina<t.cmd>t.log 回车.将获得一个t.log的文件,里边记录了该cmd的拟合信息.此处的t.cmd即上述所说的cmd文件；比照t.log中提供的信息进行筛选.即无*标示.Signal小于样点数的一半,Signal/Error最小,RTGCV,RTMSR,RTMSE最小,其中RTGCV和RTMSE的差值越大,可间接说明模型的解释率越高.

     

5、   写lapgrd  cmd

  

TMIN_JJA_82_M.sur  #输入表面文件名

1                  #表面个数

1                  #转换恢复

TMIN_JJA_82_M.cov  #输入误差表面文件名

2                #误差类型1-模型标准误差2-预测标准误差3- 95 %模型置信区间4- 95 %预测置信区间

空格

1                  #默认为1

1                  #第一栅格变量指标,通常为1

89280 3129280 1000  #X范围,1000表示输出分辨率，注意这里的范围尽量和t.cmd的范围保持一致， t.cmd稍微大一点点。但此试验中需和dem一样大，不然出不来。无解

2                   #第二栅格变量指标

3963511 6829511 1000   #Y范围

0                  #掩膜方式,0-不提供MASK 文件

2                #独立协变量数据格式,2-ARC/ INFO 格式,输入、输出栅格大小一致

amurdem_big.txt       #协变量文件(DEM 文件)

2                  #ARC/ GRID 格式

-9999              #输出规定(指示空值数据)

TMIN_JJA_82_M.grd   #输出文件名

空格

2                 #输入误差表面格式,ARCINFO 格式

-9999             #无数据指示

TMIN_JJA_82_Mcovlap.grd        #输出误差文件名

此处也需要多空行！

6、     运行lapgrd cmd

开始菜单中,在"运行"项中输入cmd 回车,把cmd运行程序打开；用cd命令进入anusplin所在的文件夹.切记,anusplin的安装文件和.dat和.cmd文件都要放在同一个文件夹中,路径必须一致.例:如果放在L:\data\anusplin文件夹下,则先输入L:回车,把当前盘符由C盘转到L:

然后输入cd data\anusplin命令,以进入该文件夹；输入lapgrd<grid.cmd>grid.log 回车.将获得一个grid.log的文件,里边记录了该cmd的拟合信息.同时，如果没有错误，将会产生两个.grd格式的文件：temp.grd和temp_cov.grd

7、将.grd文件导入argis，处理，输出tif。

• 降水插值

1. 数据准备、dem准备，与气温相同。
2. splina cmd  编写

prcp                          #文件名,可自定义,但两个cmd文件的命名不可重复

7                           #代表t的单位是mm

3                           #代表三个自变量(X，Y，高程)

0                           #代表协变量个数

0                           #代表表面样条变量个数

0                           #代表表面协变量个数

671617 3079617 0 1           #和气温一样

3774469 6388469 0 1          #和气温一样

-210 2784 1 1                #高程上下限，位（1-米m）

1000                        #代表上述的a为1000，同样是cellsize的大小

2                           #独立变量转换参数，0为不转换

2                           #样条次数（2次样条）

1                           #输出表面个数

0                           #相对误差方差个数,0 为每个表面采用统一权重

1                           #优化参数指标,通常为1

1                           #平滑参数选择方法1-GCV ,2-MSE ,3-固定值,4-GML

prcp201507.dat                #输入需要插值的数据文件.dat

110                         #站点个数,比实际个数稍多

6                           #站点标识符个数，6个字符

(a6,2f14.4,f8.1,f9.2)            #数据格式，与spss数据相通

prcp.res                     #输出残差文件

prcp.opt                     #光滑参数文件

prcp.sur                      #表面文件,用于下一步LAPGRD的输入数据

prcp.lis                       #列表文件,观测数据、预测数据、误差数据

prcp.cov                       #拟合表面系数的误差协方差文件

                    也需空格

1. 运行和气温一样

1. 编写lapgrd   cmd文件

prcp.sur…………………………………………………输入表面文件的名称（上个cmd生成）

0…………………………………………………拟合表面个数（0-输出所有表面，如果输出的多，需要依次排列，如1 2 3 4 5）（可删除此行）

1…………………………………………………计算表面类型(1-拟合表面值)

1…………………………………………………表面值转换（1-应用表面值转换，即之前的平方根）

prcp.cov…………………………………………………输入协方差表面文件名

2…………………………………………………误差类型1-模型标准误差2-预测标准误差3- 95 %模型置信区间4- 95 %预测置信区间

…………………………………………………空行，最大标准误差不填

1…………………………………………………单元中心值1,

1…………………………………………………变量索引，第一个变量为经度，写1

671617 3079617 1000…………………………………………………X范围,1000表示输出分辨率，注意这里的范围尽量和t.cmd的范围保持一致，t.cmd稍微大一点点。

2…………………………………………………变量索引，第二个变量为纬度，写2

3774469 6388469 1000…………………………………………………Y范围,1000表示输出分辨率，注意这里的范围尽量和t.cmd的范围保持一致，t.cmd稍微大一点点。

0…………………………………………………掩膜方式（0-无掩膜文件）

2…………………………………………………独立协变量数据格式,2-ARC/ INFO 格式,输入、输出栅格大小一致

amurdem_big.txt…………………………………………………输入dem文件

2…………………………………………………输出栅格文件格式（2- ARC/ GRID 格式）

-9999…………………………………………………dem中定义无效值

prcp.grd…………………………………………………输出文件名

…………………………………………………空行，默认输出二进制文件格式，占内存少速度快

2…………………………………………………输出表面误差栅格文件格式

-9999…………………………………………………误差栅格特殊值

prcp ­_covlap.grd…………………………………………………误差栅格命名

1.  运行和后续与气温操作一样