一、数据介绍土壤质地类型数据来源于世界土壤数据库(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/

二、数据下载085de399d5a0e622679774d713324ec2.png数据下载说明FAO SOILS PORTAL下载网址内容如上图,其中全球土壤数据查看器(见下图)仅仅只能查看全球土壤信息,并不能提供下载服务。我们只需要下载后面两个链接即可。d71bd9f2e2d725b03d3c35da7ad5ccfe.png全球土壤数据查看器下载解压后文件内容如下,包括一个BIL格式的全球栅格数据,和HWSD的数据库文件。a070f79f9e02df113d13fdd2d6d4ce64.png下载结果查看三、数据处理

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

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

26d78b2e538abd128fc9ddc3951df94f.png

数据库文件导出Excel

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

edb2d41ef4300c439780f537c71e4bd0.png

导出文件格式选择

e22742dfc42535ef4ccf6cabef09fff0.png导出数据查看查看官方手册,可知MU_GLOBAL是提供数据库文件和栅格文件连接的ID。f16a2c178198245d810b2854a53c5541.png官方手册说明1SU_SYM90和CODE是在FAO90分类系统下的土壤类型代码。376329773fbd91f85c9617fb1a8196c3.png官方手册说明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所在的文件夹,直接将栅格文件拖到左边窗口即可打开。3e9da93a5108651fece10863fcf70991.pngc7eb1b4c1b4b41a347890850a376b241.png栅格文件导入ArcGIS由于原始栅格没有空间参考信息,所以需要定义一下投影,选择WGS_1984地理坐标系(可根据需要修改)。60c3d7d706848eea3c69a05066080ea3.png定义投影坐标系

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

90d38ec43a19591123ace971c4fdc4c0.png

单一波段栅格属性创建

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

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

6cc75b6ee712a5c04178797fe0071386.png

栅格转点创建属性表

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

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

c681b83e86ddbfa9517d38a8a3b50f01.png

链接数据属性表

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

b087c9cb44fcdfc3ec322c69e94f1086.png

山东省矢量图

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

c4e75b8ae0bdf5a46a593c7e05984838.png

山东土壤类型情况

然后选择自己需要的属性字段栅格转点。在这里我选择了土壤类型的CODE字段。aefea1c23669bab08e43d4017e98cbe1.png栅格转点最后点转栅格,右键栅格图层导出tiff格式文件,导出前确认下有无空间参考,如果没有就再定义一遍。331a6c8118c1b6b8588f280c856f2977.png9036a23e97da28e1a6a846cf7464e022.png点转栅格在ENVI中打开,确认是带有地理坐标的,然后接下来就可以根据自己需要转化投影了。441ad531c519f40d08be189de5948bca.png确认数据有空间参考四、结果对比再和HWSD Viewer中分类情况对比下,完全一致,说明处理没有错误。b93ab2ad1025bc56951cd4c3d28f996e.png8400708de0e558c7366a00e4ccd56e3c.png输出结果与官方对比

来源:CSDN,江北20190411;

         新浪博客,CityThoth;

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

如有侵权请联系删除

e4d0b41d5de62a14faa02012da948d37.png

298f71463656ea936943d542bd6828ff.png

0b0c743120dfc1303a6fc896a3109fe9.png

0b0c743120dfc1303a6fc896a3109fe9.png

往期精选

09 October 2020

●收藏|这里有一份关于自然资源负债表编制的干货帖●关注|《第三次全国国土调查耕地资源质量分类工作方案》印发●关注|广东省全面铺开国家公园建设,打造粤港澳大湾区生态屏障●聚焦|草木蕃成为“红树林保护与恢复国家创新联盟”首批成员单位●草木蕃十周年系列活动之趣味运动会视频回顾●关注|《市级国土空间总体规划编制指南(试行)》公开征求意见●新发现|海南热带雨林国家公园体制试点区发现19个新物种●建军93周年|草木蕃党支部八一建军节专题党课●十周年系列活动|草木蕃首届趣味运动会圆满落幕!●开课啦!|“菁英院”三期新员工培训正式开班授课●添新彩|广东大亚湾省级自然保护区首现尖笔帽螺“潮”

35fca045081fcd2c877fb69f8e11e4aa.png

Logo

腾讯云面向开发者汇聚海量精品云计算使用和开发经验,营造开放的云计算技术生态圈。

更多推荐