世界土壤数据库（HWSD）全球土壤数据下载处理

点击上方“Ai尚研修科研技术平台”选择加"星标"置顶

重磅干货，第一时间送达

热课推荐

TOUGH2系列建模方法与综合实践技术应用高级培训班邀请函 直播时间：3月18日-3月21日

地下水数值模拟技术应用与地下水环评报告编制方法实践培训班

直播时间：3月25日-3月30日

关于举办水文地球化学模拟PHREEQC实践技术应用高级培训班

直播时间：4月23日-4月25日地下水数值模拟GMS实践技术应用培训班预约报名

地下水数值模拟Visual modflow Flex应用培训班

直播时间：4月22日-4月27日

环境土壤物理Hydrus1D2D3D模型实践应用培训班

直播时间：4月16日-4月19日

最新:AQUATOX水环境与生态模型实践技术应用培训班

直播时间：3月19日-3月21日

第二十三期  SWAT模型在水文水资源、面源污染模拟中的实践技术应用及典型案例分析线上直播课

直播时间：3月26日-3月29日

城市内涝水力模型建立方法与风险评估实践技术应用高级培训班3月26日-3月28日

HEC-RAS软件建模实践方法培训班

直播时间：3月27日-4月1日

新《生产建设项目水土保持方案技术审查要点》要求下水土保持方案编制高级实践技术应用培训班通知

直播时间：4月10日-4月11日

EFDC建模及在地表水环境评价、水源地划分、排污口论证等领域中的实践技术应用培训班

直播时间：4月15日-4月18日

“R语言+遥感”的水环境综合评价方法实战培训班

直播时间：4月16日-4月18日

第二期Delft3D水动力-富营养化模型实践高级培训班

直播时间：4月23日-4月26日

一、数据介绍

土壤质地类型数据来源于世界土壤数据库（HWSD，Harmonized World Soil Database）。该数据是空间分辨率为公里的网格数据，提供了各个格网点的土壤类型、土壤相位、土壤理化性状等信息。

国内数据下载地址
http://westdc.westgis.ac.cn/data/611f7d50-b419-4d14-b4dd-4a944b141175
全球数据下载地址：
http://www.fao.org/soils-portal/soil-survey/soil-maps-and-databases/harmonized-world-soil-database-v12/en/

二、数据下载

数据下载说明

FAO SOILS PORTAL下载网址内容如上图，其中全球土壤数据查看器（见下图）仅仅只能查看全球土壤信息，并不能提供下载服务。我们只需要下载后面两个链接即可。

全球土壤数据查看器

下载解压后文件内容如下，包括一个BIL格式的全球栅格数据，和HWSD的数据库文件。

下载结果查看

三、数据处理

本次处理目的得到山东省的土壤类型分类图

打开数据库文件，选择HWSD_DATA右键导出为EXCEL。

数据库文件导出Excel

文件格式最好选择如下，否则在导入ArcGis的时候容易报错。

导出文件格式选择

导出数据查看

查看官方手册，可知MU_GLOBAL是提供数据库文件和栅格文件连接的ID。

官方手册说明1

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。

在ArcGIS中链接HWSD所在的文件夹，直接将栅格文件拖到左边窗口即可打开。

栅格文件导入ArcGIS

由于原始栅格没有空间参考信息，所以需要定义一下投影，选择WGS_1984地理坐标系(可根据需要修改)。

定义投影坐标系

为栅格影像创建属性表，方法如下：
有时候一个栅格数据有属性值，但是“打开属性表”不可用，此时候可以通过ARCGIS中的“创建属性表“工具解决。需要注意的是此栅格数据集必须只有单一波段。不能为具有 32 位浮点像素类型的栅格数据集构建栅格属性表。

单一波段栅格属性创建

如果是浮点型型像素的栅格，可以通过以下两种方法来解决：

1、把栅格转点，”栅格转点“再用”值提取至点“的工具，获得属性表。（或者“采样”工具也可）

栅格转点创建属性表

2、直接重分类，获取属性表。

右键栅格图层连接和关联，连接，选择之前导出的HWSD_DATA.xls文件。

链接数据属性表

在这里我加载了山东的WGS_1984地理坐标的矢量文件。

山东省矢量图

工具箱中裁剪栅格选择山东矢量。然后按照SU_CODE90显示一下，可以看到山东土壤类型情况。

山东土壤类型情况

然后选择自己需要的属性字段栅格转点。在这里我选择了土壤类型的CODE字段。

栅格转点

最后点转栅格，右键栅格图层导出tiff格式文件，导出前确认下有无空间参考，如果没有就再定义一遍。

点转栅格

在ENVI中打开，确认是带有地理坐标的，然后接下来就可以根据自己需要转化投影了。

确认数据有空间参考

四、结果对比

再和HWSD Viewer中分类情况对比下，完全一致，说明处理没有错误。

输出结果与官方对比

整理/编辑：WXW

来源：CSDN，江北20190411；

         新浪博客，CityThoth；

         国家青藏高原科学数据中心

●数据及相关资源共享
●1、【免费领取】中国土地利用数据1980-2015年2、大数据！全国最新地级市和气象站空气质量和60年气候详细数据下载3、最全的全国分省、市、县、乡镇行政区划矢量图（shp格式）免费下载4、干货丨零基础学习R语言 视频+讲义免费下载5、免费资源丨SCI写作全攻略 6、[数据分享]中国地貌矢量数据，撰写报告做图必备数据免费下载7、【免费资源】多个自然资源（大气、土地、水等）面板数据大放送、时间跨度长、地区分布广 快速下载！8、【数据分享】中国各区县PM2.5数据和全国各区县二氧化碳排放数据丨免费下载！

关注“Ai尚研修科研技术平台”公众号，查看更多数据分享，视频教程，学习资料！     

2021年度近期会议

3-4月课程都在这里，点击下边图片了解详情

—END—

关于Ai尚研修