-
数据介绍
-
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。
-
官网(https://gaez.fao.org/pages/hwsd)下载数据,包含数据查看器、土壤数据(数据库和栅格数据)、帮助文档
-
数据库——主要两个数据是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)为例,选择输出路径,导出为栅格即可。