来源:GIS前沿

一、数据介绍

土壤属性表主要字段包括(图1):详细描述请参考Harmonized World Soil Database (version 1.1).pdf文件,其中以T开头的土壤属性表示土壤上层的属性(0-30cm),以S开头的土壤属性代表土壤下层的土壤属性(30-100cm)。

图1土壤属性表内容示意图

二、数据下载

世界土壤数据库(HWSD)的中国土壤数据集(v1.1)从国家青藏高原科学数据中心(http://data.tpdc.ac.cn/zh-hans/)下载获得,无需注册。其中数据共有2个版本,在官方网站输入关键词“土壤”即可搜索得到。

数据下载解压后可以得到栅格数据、土壤属性表和数据使用说明等文件(图3),注意使用该数据时,请遵循相关引用规定,合理合法运用(非商业用途)。

图3数据内容

*本教程仅供参考,不足之处欢迎斧正。

三、数据处理

3.1 土壤属性表的转换

打开数据库文件,右键选择HWSD_DATA导出为Excel文件(图4),文件格式要选择97-2003类型,否则ArcGIS容易出错。

图4 导出Excel

导出的数据在Excel中打开找到“MU_GLOBAL”,该字段是与栅格图中ID字段一一对应的(图5)。

图5土壤属性表

3.2 加载数据

打开ArcMap,加载数据并查看数据是否有坐标,本实验下载的数据是自带了WGS_1984坐标系的(图6)。如果没有坐标则利用【定义投影】工具投影到对应坐标系即可。

图6加载栅格数据

3.3 栅格数据构建属性表

在工具箱中打开【数据管理工具】——【栅格】——【栅格属性】,找到【构建栅格属性表】,为栅格创建属性表,为后续“链接”属性字段奠定基础(图7)。值得注意的是,该栅格数据仅有单一波段,若是具有32位浮点像素类型的栅格数据则不能构建属性表。

图7构建栅格属性表

3.4 数据链接

在图层“HWSD_China_Geo”右键,选择【连接和关联】,设置连接的对应字段(图8)。从数据说明中可以查到,字段“MU_GLOBAL”与栅格图中Value字段是一一对应的。

图8 数据连接

3.5 数据导出

将数据与土壤属性表连接后,一定要单击右键导出栅格数据(图9)。导出数据时设置像元大小为1km×1km。

图9导出连接后的栅格数据

3.6 栅格裁剪

根据研究需要,利用研究区矢量边界对栅格进行裁剪,以贵州省为例。打开【空间分析工具】——【提取分析】——【按掩膜提取】,输入贵州省边界数据和栅格数据,最终得到结果(图10)。

图10贵州省土壤类型

该栅格数据的属性表中对应了贵州省的土壤类型属性的相关参数,打开属性表可以查看,也可以在【符号系统】中根据不同的需要进行制图或表达(图11)。

图11贵州省土壤类型属性表