防风固沙功能中多年雪深数据的处理方法

数据来源

中国雪深长时间序列数据集（1978-2012）来源于寒区旱区科学数据中心（http://westdc.westgis.ac.cn/data/f957004d-0756-4796-b4c5-c02265be64e1），
该数据集包括1978年10月24日到2012年12月31日逐日的中国范围的积雪厚度分布数据，其空间分辨率为25km，采用GCS_WGS_1984地理坐标系统。数据集格式为txt文本格式，每日雪深数据为一个txt，一年365/366个txt。

txt转raster

为了对雪深数据进行计算，需要先将txt转换成raster。格式转换可以采用arcgis中conversion tools→to raster→ASCII to raster实现转换，小量数据还可以右键工具batch进行批量转换，数据转换过程中需注意转换数据不能存储在gdb数据库中，也不能存储过深，否则数据处理将报错。具体示例如下：

对于如此大量的数据需调用GIS中的pyhton功能，具体如下：

然后等待两小时左右，数据将全部转换完成。再将转换出来的raster按照中国边界进行剪裁。为了确定全国各地各月份雪深大于5mm的天数，需要对雪深数据进行重新分类，将0-5mm赋值为0，>5mm赋值为1，此过程需注意分类的最高限可赋一个大于全国的最大雪深的数，比如说50cm，且输出栅格需要注意输出到数据库gdb中，而且命名最好有字母，否则将报错。再将全国每月雪盖大于5mm的天数统计求和，在此数据的table中新加一列，属性为float，通过栅格计算器用value值除以当月天数得到雪盖因子，使用weighted sum工具将因子数据提取出来，最终多年1-12月雪盖因子。