HWSD2.0数据处理

• 数据介绍

• D1：0-20cm、D2：20-40cm、D3：40-60cm、D4：60-80cm、D5：80-100cm、D6：100-150cm、D7：150-200cm

• HWSD2_SMU是基础图层，提供了土壤测绘单元中占比（Share）最大的土壤类型的第一层（D1）的土壤理化特征。

• HWSD2_LAYERS是所有土壤类型的D1-D7层的土壤理化特征。HWSD2_SMU_ID是提供数据库文件和栅格文件连接的唯一ID。

1. 官网（https://gaez.fao.org/pages/hwsd）下载数据，包含数据查看器、土壤数据（数据库和栅格数据）、帮助文档

2. 数据库——主要两个数据是HWSD2_SMU和HWSD2_LAYERS，对应字段的介绍可看对应的._METADATA。

单因子tiff提取

• 先根据HWSD2_SMU确定占比最大的土壤类型，然后再在HWSD2_LAYERS筛选出占比最大的土壤类型的7层土壤理化特征 。

  • 先将.mdb中的HWSD2_SMU导出为.xlsx

  • HWSD2_LAYERS导出为txt（逗号分割）（由于数据量比较大，建议导出为txt，xlsx也可以，但是别导出xls，会缺失数据）

  • 利用R语言进行筛选。R+RStudio（4.3.2）

  

• 将筛选好的HWSD2_LAYERS_Max.xlsx另存为xls(记得加上标题)，便于与ArcGIS中的栅格连接。（D1-D7.xls）

  • arcgis——加载研究区矢量边界，HWSD2.bil，对边界做缓冲区（1000m），按掩膜提取（以缓冲区边界进行提取土壤数据）

  • 导入裁剪后的栅格文件，右击选择连接和关联→连接（若是灰色或者没有属性表，可以先在属性→符合系统，显示唯一值，获取属性表），分别连接D1-D7层，分别进行不同土壤层的因子提取。

  • 这里以D6层为例：选择连接D6表，连接字段选择HWSD2_SMU_ID，将D6层的数据导入。

  • 右键导出栅格数据（D1-D7）

单因子栅格导出---选择3D Analyst工具→栅格重分类→查找表工具，选择需要的字段，这里以有机碳（OC）为例，选择输出路径，导出为栅格即可。