数据介绍
土壤质地类型数据来源于世界土壤数据库(该数据由联合国粮农组织维也纳国际应用系统研宄所、荷兰、中国科学院南京土壤研宄所、欧洲委员会联合研宄中心共同构建)。该数据是空间分辨率为公里的网格数据,提供了各个格网点的土壤类型、土壤相位、土壤理化性状等信息。
之前国内寒旱所提供中国区域的数据下载,不过最近那个网站进不去了,链接还是放到这里,万一哪天又好了呢。
国内数据下载地址
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的数据库文件。
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/f39ab3d8908a5ea54e8abcbb292371e4.png)
数据处理
本次处理目的得到山东省的土壤类型分类图
打开数据库文件,选择HWSD_DATA右键导出为EXCEL。
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/4428043bd334fcd8d41df688aeb57738.png)
文件格式最好选择如下,否则在导入ArcGis的时候容易报错。
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/3d1877ea12a379dce781f2a395877cb9.png)
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/5cca27e0ff79d14b5aed64f1c107ab1d.png)
查看官方手册,可知MU_GLOBAL是提供数据库文件和栅格文件连接的ID
SU_SYM90和CODE是在FAO90分类系统下的土壤类型代码。
在ArcGis中链接HWSD所在的文件夹,直接将栅格文件拖到左边窗口即可打开。
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/9a8374312d370045dbb3b32e3996b606.png)
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/4b59044afc31082ec5041c834c8588d9.png)
由于原始栅格没有空间参考信息,所以需要定义一下投影,选择WGS_1984地理坐标系。
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/e479b2fc250c6b6df78af773c78607cc.png)
为栅格影像创建属性表,创建方法可见下面博客。
ARCGIS为栅格数据创建属性表
右键栅格图层连接和关联,连接,选择之前导出的HWSD_DATA.xls文件。
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/c7cfadfe2dc84068d2bdf6407dd0cb89.png)
在这里我加载了山东的WGS_1984地理坐标的矢量文件
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/274f1a56289a7bf3a478c5b9fe3e161f.png)
工具箱中裁剪栅格选择山东矢量。然后按照SU_CODE90显示一下,可以看到山东土壤类型情况
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/85aa4b08a7370d015dcf863cf7074044.png)
然后栅格转点,选择自己需要的属性字段。在这里我选择了土壤类型的CODE字段
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/6d78c249f342a3f01455ecadbceea76d.png)
最后点转栅格,右键栅格图层导出tiff格式文件,导出前确认下有无空间参考,如果没有就再定义一遍。
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/3113597c418424b39d3a682cf45403a5.png)
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/acf9486591662b92811bb2169a921177.png)
在ENVI中打开,确认下是带有地理坐标的,然后接下来就可以根据自己需要转化投影了。
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/882731dba92d518719e5fac4d30f4f4e.png)
再和Viewer中分类情况对比下,完全一致,说明处理没有错误
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/c2beed93e7d2d15a72a22b65f72f01e5.png)
![图片替换文本](https://i-blog.csdnimg.cn/blog_migrate/f0cced8188fa77daeb813e8d061fc15a.png)