参考链接：https://zhuanlan.zhihu.com/p/477677143
1. 数据简介
土壤质地类型数据来源于世界土壤数据库（HWSD，Harmonized World Soil Database）。该数据是空间分辨率为公里的网格数据，提供了各个格网点的土壤类型、土壤相位、土壤理化性状等信息。
数据库下载地址：https://www.fao.org/soils-portal/soil-survey/soil-maps-and-databases/harmonized-world-soil-database-v12/en/
图片.png
2. 数据下载
FAO SOILS PORTAL下载网址内容如上图，其中全球土壤数据查看器只能查看全球土壤信息，并不能提供下载服务。我们只需要下载后面两个链接即可。
图片.png
下载解压后文件内容如下，包括一个BIL格式的全球栅格数据，和HWSD的数据库文件。
图片.png
3. 数据处理
下载的HWSD.mdb文件可以用Microsoft Access 或 MDB viwer plus软件打开查看并导出xls文件，也可以直接在arcmap中打开。
查看官方手册，可知MU_GLOBAL是提供数据库文件和栅格文件连接的ID。
SU_SYM90和CODE是在FAO90分类系统下的土壤类型代码。
官方手册说明2这里给大家附上数据说明：数据来源于联合国粮农组织（FAO）和维也纳国际应用系统研究所(IIASA)所构建的世界土壤数据库。中国境内数据源为第二次全国土地调查南京土壤所所提供的1：100万土壤数据。该数据可为建模者提供模型输入参数，农业角度可用来研究生态农业分区，粮食安全和气候变化等。数据格式：grid栅格格式，投影为WGS84。采用的土壤分类系统主要为FAO-90。
土壤属性表主要字段包括:
SU_SYM90（FAO90土壤分类系统中土壤名称)；
SU_SYM85(FAO85分类);
T_TEXTURE(顶层土壤质地);
DRAINAGE(19.5);
REF_DEPTH(土壤参考深度);
AWC_CLASS(19.5);
AWC_CLASS(土壤有效水含量);
PHASE1: Real (土壤相位);
PHASE2: String (土壤相位)；
ROOTS: String (到土壤底部存在障碍的深度分类)；
SWR: String (土壤含水量特征)；
ADD_PROP: Real (土壤单元中与农业用途有关的特定土壤类型)；
T_GRAVEL: Real (碎石体积百分比)；
T_SAND: Real (沙含量)；
T_SILT: Real (淤泥含量)；
T_CLAY: Real (粘土含量)；
T_USDA_TEX: Real (USDA土壤质地分类)；
T_REF_BULK: Real (土壤容重)；
T_OC: Real (有机碳含量)；
T_PH_H2O: Real (酸碱度)
T_CEC_CLAY: Real (粘性层土壤的阳离子交换能力)；
T_CEC_SOIL: Real (土壤的阳离子交换能力)
T_BS: Real (基本饱和度)；
T_TEB: Real (交换性盐基)；
T_CACO3: Real (碳酸盐或石灰含量)
T_CASO4: Real (硫酸盐含量)；
T_ESP: Real (可交换钠盐)；
T_ECE: Real (电导率)。
其中以T_开头属性字段表示上层土壤属性（0-30cm），以S_开头属性字段表示下层土壤属性（30-100cm）。
具体属性值代表何意义请参考文件夹下说明文档.pdf及数据库.mdb。
使用arcmap打开HWSD.bil文件
图片.png
由于原始栅格没有空间参考信息，所以需要定义一下投影，选择WGS_1984地理坐标系(可根据需要修改，我下到的好像有)。
图片.png
为栅格创建属性表
打开arctoolbox→data management tools →raster → build raster attribute table
图片.png
连接属性表
右键图层 →joins and relates →join
图片.png
图片.png
点击第二行那个文件夹，点击HWSD.mdb文件，找到HWSD_DATA双击
图片.png
连接字段选择MU_GLOBAL，点击OK就连接好了
单因子栅格图层导出
使用3D analyst tools工具包
3D analyst tools → raster reclass → lookup
图片.png
之后选择相应的字段提取相应的因子图层
图片.png
这样就导出了栅格图层，在导出完成之后，在arcmap中打开该图层右键
网友评论