HWSD全球土壤数据下载处理

数据介绍

      土壤质地类型数据来源于世界土壤数据库（该数据由联合国粮农组织维也纳国际应用系统研宄所、荷兰、中国科学院南京土壤研宄所、欧洲委员会联合研宄中心共同构建）。该数据是空间分辨率为公里的网格数据，提供了各个格网点的土壤类型、土壤相位、土壤理化性状等信息。
      之前国内寒旱所提供中国区域的数据下载，不过最近那个网站进不去了，链接还是放到这里，万一哪天又好了呢。

国内数据下载地址
http://westdc.westgis.ac.cn/data/611f7d50-b419-4d14-b4dd-4a944b141175
全球数据下载地址：
http://www.fao.org/soils-portal/soil-survey/soil-maps-and-databases/harmonized-world-soil-database-v12/en/

数据下载

FAO SOILS PORTAL下载网址内容如上图，其中全球土壤数据查看器（见下图）仅仅只能查看全球土壤信息，并不能提供下载服务。我们只需要下载后面两个链接即可。

下载解压后文件内容如下，包括一个BIL格式的全球栅格数据，和HWSD的数据库文件。

数据处理

本次处理目的得到山东省的土壤类型分类图

打开数据库文件，选择HWSD_DATA右键导出为EXCEL。

文件格式最好选择如下，否则在导入ArcGis的时候容易报错。

查看官方手册，可知MU_GLOBAL是提供数据库文件和栅格文件连接的ID

SU_SYM90和CODE是在FAO90分类系统下的土壤类型代码。

在ArcGis中链接HWSD所在的文件夹，直接将栅格文件拖到左边窗口即可打开。

由于原始栅格没有空间参考信息，所以需要定义一下投影，选择WGS_1984地理坐标系。

为栅格影像创建属性表，创建方法可见下面博客。
ARCGIS为栅格数据创建属性表
右键栅格图层连接和关联，连接，选择之前导出的HWSD_DATA.xls文件。

在这里我加载了山东的WGS_1984地理坐标的矢量文件

工具箱中裁剪栅格选择山东矢量。然后按照SU_CODE90显示一下，可以看到山东土壤类型情况

然后栅格转点，选择自己需要的属性字段。在这里我选择了土壤类型的CODE字段

最后点转栅格，右键栅格图层导出tiff格式文件，导出前确认下有无空间参考，如果没有就再定义一遍。

在ENVI中打开，确认下是带有地理坐标的，然后接下来就可以根据自己需要转化投影了。

再和Viewer中分类情况对比下，完全一致，说明处理没有错误