来源: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 贵州省土壤类型属性表