ANUSPLIN插值专题系列—线性插值及研究区站点提取

无论是风里，还是在雨里，我都在这里守候着你～

线性插值及研究区站点选取

       “气象”小伙伴们，小编上期推送的《ANUSPLIN插值专题系列—前期数据准备》么，今天小编就接着来和大家探讨探讨ANUSPLIN插值的研究区站点提取！主要包括以下三个部分：

       （1）线性插值；

       （2）研究区站点的选取；

       （3）研究区站点数据提取。

       涉及的软件包括Excel、Matlab、Arcgis和SPSS。

       通过（二）后我们初步得到了全国气象站点各个气象要素的ANUSPLIN插值的格式，但如果做长时间序列的话，会发现不少站点的记录信息不全，存在着缺失现象，需要通过插补来进行弥补，所幸这样的站点并不多，本文采用线性插值法对缺测站点进行插值。插值完成后，我们可以得到完整的中国各个气象站点的各个气象要素的ANUSPLIN插值的格式，如果研究区的范围并不是全国，则需要从中进行提取，本文提供相应的程序来进行提取，同时该过程中会涉及到研究区站点的选取。

线性插值

       当某个气象站某年的缺测站点数据过多时，不管采用何种方法得出的结果精确度都不高，可选择直接舍弃，但现实中这样的气象站点并不多，这种舍弃不会对最终结果造成大的影响。缺测站点数据较少时，如低于5%，采用线性插值法进行内插，比如第一天和第三天都存在数据，而第二天缺测，则用第一天和第三天的数据求平均后得到第二天的数据。具体代码（https://www.jianshu.com/p/2869eaa5c208），输出的结果中1-4列的表头是缺失的，需要自己手动添加，分别是station，x，y和dem.

研究区站点的选取

       本文以黄土高原为例，进行研究区气象站点的选取。首先加载全国站点和黄土高原的DEM文件，如下图所示：

       右击DEM文件，打开属性，进行符号系统中，将Nodata部分的颜色设置成黄色，如下图（1）所示，出来的结果如下图（2）所示。

（1）

（2）

       然后通过选择工具，将在黄色部分内的站点全部选中，具体如下。

       选择完成后，进入属性表，在图中红色部分选择Export，进入下图右边的窗口，在红色部分选择selected records，输出时保存类型选择dBase Table格式，然后用excel打开其中的DBF文件，将站点另存为excel，共132个点，见右下图。

研究区站点数据提取

       接下来通过匹配站点名称来自动获取研究区站点的气象插值格式，具体代码见(https://www.jianshu.com/p/877219f6e8a0)。本文将时间段设置为1980-2014,1980-2014年黄土高原降水的ANUSPLIN结果则直接出来，如下图所示。出来的结果中1-4列没有表头，需要手动添加，分别是station，x，y，dem。

       通过上述这些操作，即可以得到研究区的ANUSPLIN插值格式的excel数据，后续需要通过SPSS转换为固定ASCII编码的DAT格式及脚本的编写和插值结果的处理，敬请期待下期。

资源仅供学术交流使用，严禁商用！如有侵权，请联系小编微信：ylc107919293

推文期数：2018120

下期预告：重点大学（研究所） 考研招生专业目录（第一批）

责任编辑：姚春宇  尹礼唱  卢雅焱

推文审核：张天舒  梁龙武  骆丹云

总审核：学术无界顾问团

往期回顾：

广告也精彩，小编不易，欢迎点击广告支持小编哦！