利用ANUSPLIN进行气象插值步骤详解——以气温为例
1. ANUSPLIN介绍
ANUSPLIN插值软件是澳大利亚国立大学研发的一款专门用于气象插值的一个软件,该方法采用澳大利亚国立大学(Australian National University)首字母和样条函数(splin)结合的方法命名。ANUSPLIN包的目的是提供一种使用薄板平滑样条对有噪声的多变量数据进行透明分析和插值的工具。
本教程以气象站点数据为依据,将高程数据作为协变量,进行空间连续插值。
2.数据的准备及处理
2.1 软件和数据的准备
需要准备的软件和数据有:
- ArcGIS
- SPSS
- Excel
- Notepad++
- ANUSPLIN插值软件(SPLINA.exe & LAPGRD.exe等)
- 全国的DEM数据
- 研究区矢量数据
2.2 数据处理
2.2.1 气象数据的准备和处理
处理流程:①坐标转换 ②气象站展点③选择插值区域④导出插值区域内的站点数据⑤整理数据格式
读者可以通过中国气象数据网获取相应气象数据。
获取的气象数据格式和数据说明如下:
①坐标转换 以2017年12月气温数据为例,我们要将各站点的度和分转换成十进制经纬度坐标,以便空间展示。
将下载的TXT文件用Excel打开进行处理,具体操作如下:
将多余数据删除,仅保留台站号、经度、纬度、海拔和气温数据,然后编辑公式,将十进制经纬度数据替换掉原来的度和分,操作如下:
经度和纬度按照相同的计算公式自动填充后,得到十进制的经纬度信息,将原来数据进行替换,结果如下图:
将其另存为TXT格式,以便在ArcGIS中打开,操作如下:
②在ArcGIS中进行展点 展点过程中的投影信息十分关键,注意不要弄错。
首先加载TXT数据,操作如下:
然后设置图层的地理坐标系和投影坐标系,右击图层后选择属性,在坐标系统中设置投影信息。可以只选择投影坐标系,因为选择投影坐标系之后的地理坐标系会进行自动填充,反之则不能,操作如下: