基于Anusplin制作降雨数据插值操作中可能遇到的问题

1.暂停正在运行的cmd ，输入“Ctrl+C”键。

2.RDATA - NO DATA POINTS WITHIN SPECIFIED LIMITS

        调整高程范围（可略大一些）。

3.GETLIM - INCOMPATIBLE GRID SPACING

        不兼容的网格间距，在CMD文件中修改成正确的分辨率 。

4.Cmd单位可米、度，本文为度（°）。（需注意单位不同需调整坐标系、降雨数据对应X、Y计算几何、dem范围、两个cmd文件中的X、Y范围与对应单位等）

5.原始dem数据合并后需转为1000m分辨率，因为后续操作需把dem转为txt文件，数据太大notepad++打不开。

6.当站点较少（需大于7个），得到的插值文件为空时，可修改lapgrd.cmd中的标准差，可删去标准差范围，空一行（代表标准差不指定），但这样得到的数据精度回受影响，有条件可尝试增加站点数据。

7.安装Anusplin软件包时，需将bin文件中splina.exe、lapgrd.exe放到后续实验数据所在文件中，其中，可能需要修改电脑的环境变量，将splina.exe、lapgrd.exe所在文件增加到环境中。

8.站点高程需从dem中提取，需与dem中匹配才能顺利运行、插值。（笔者站点高程原为爬虫获取，但站点高程与dem中香味位置高程并不完全一致，有误差，执行两个cmd文件时就会报错，无法运行）

9.两个cmd文件中多余的空格与注释需删除再运行。

10.SPSS软件导出dat文件时，注意编码格式，默认为UTF-88编码，需改为本地编码格式。

11.splina.cmd、lapgrd.cmd文件中都有经纬度范围，splina.cmd的经纬度范围可略大于dem的经纬度范围，、lapgrd.cmd的经纬度范围必须与dem文件一致。

12.splina.cmd、lapgrd.cmd文件末都需空六行或以上。

13.当站点超于dem范围时，需新增dem重新合并直至包括全部站点，或将多余站点移除。

14.插值完后降水存在负值时，通过arcgis中的栅格计算器来进行异常值处理。此时的负值是由于引进了协变量DEM造成的，可以认为出现的负值的地方降水为0的，在栅格计算器中输入Con("**.tif"<0,0,"**.tif")，**.tif为插值后的降雨数据。 也可通过matlab进行批量处理

[a,R]=geotiffread('F:\项目\dem.tif');%先导入投影信息

info=geotiffinfo('F:\项目\dem.tif');

e=dir('*.tif'); for i=1:size(e,1)

      data=importdata(e(i).name);

      data(data<0)=0;

      geotiffwrite(e(i).name,data,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);

end

参考：基于ANUSPLIN的批量气象插值-从数据处理到最终结果（6） - 简书

15.注意坐标系的选择，当选择站点坐标以米为单位时，dem与降雨数据均需设置投影坐标，且投影坐标需一致。

        可定义 Albers 投影信息

  1）打开 [地图处理-ArcToolBox-数据管理工具-投影和变换-创建空间参考] 工具，利用

<空间参考> 属性（或其他方法打开 <空间参考> 属性），新建一个投影坐标系。

  2）投影坐标系参数设置。<名称> 项为新建坐标文件的名称，这里命名为

<Asia_North_Albers_WGS84_LCR> ， <投影名称> 选择 Albers，<Central_Meridian> 设置为 105.0，<Standard_Parallel_1> 设置为 25.0，<Standard_Parallel_2> 设置为 47.0，线性单位 < 名称> 选择 Meter，<每单位米数> 设置 1，<地理坐标系> 选择 WGS 1984，其余参数为默认。

参考：气象 | 基于ANUSPLIN的气象数据插值过程_anusplin插值_洛的地理研学的博客-CSDN博客

        

本文参考：

利用ANUSPLIN进行气温插值_卍晗:-*◎的博客-CSDN博客

http://t.csdn.cn/cLQWM

气象 | 基于ANUSPLIN的气象数据插值过程_anusplin插值_洛的地理研学的博客-CSDN博客

基于ANUSPLIN的批量气象插值-从数据处理到最终结果（6） - 简书