#需要的HWSD数据库土壤中文名称表格放在资源库里,第一次弄不知道怎么搞,如果找不到,可以私我,发给你。

将裁剪的栅格数据进行重投影,让裁剪的栅格图与研究区的坐标一致。下载的土壤数据类型为地理投影坐标系,通过ArcToolbox Data Management Tools  Projections and Transformations Raster    Project Raster工具实现。

将地理坐标投影成投影坐标,根据需求进行投影。因为我的研究区域在北半球的北纬50N区域,所以我选择的投影坐标为WGS-1984_UTM_Zone_50N。

我们之前提到土壤数据栅格属性中的数字,不同的数字代表不同的土壤类型,现在我们来对这些土地类型进行处理。找到我们之前下载的文件夹,打开HWSD.mdb文件,文件的打开需要利用office的Access打开。

打开后,选择HWSD_DATA,右键导出,导出到Excel。表格名称代表的含义有:

SU SYM90 (FAO90土壤分类系统中土壤名称)

SU SYM85(FAO85分类);TTEXTURE(顶层土壤质地)

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 H20: Real(酸碱度)

T CEC CLAY: Real(粘性层壤的阳离子交换能力)

TCEC 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)。

回到Arcgis中,点击图层鼠标右键打开Open Attribute Table,打开属性表,属性表打开后,点击表选项,选择文件夹后导出。

文件导出成.dbf格式,选择记事本形式打开,打开的界面如下:

将数据复制到新建的Excel文件中。根据我们之前导出的HWSD_DATA的Excel表格筛选我们栅格中需要的土地类型(即裁剪过后我们导出.dbf文件中的数字代码),筛选的列为MU_GLOBAL,这一列储存土壤类型代码,筛选出来后,重新建立一个表格1,如下所示:

表格1建立完成后,我们进一步计算COUNT的占比。COUNT就是该土地利用所占的表格,我们在Arcgis属性表中可以看到这个数据,将其导出粘贴表格1中。将研究区count加起来,用每个土壤利用类型的count除以总数得到占比,得到占比数据后先放一边。

继续对土壤类型进行重分类,继续对SU_SYM90里的土地利用类型,例如LPk为土壤类型的简称。打开HWSD数据库土壤中文名称表格,找到对应的土壤名称的简称,复制下来。

我们将研究区的土壤筛选出来,新建表格2。

结合表格1和表格2,将土壤数据整理成下图形式:

土壤数据整理结束后,可以进行重分类。简而言之,就是把同一类型的土壤归为同一类,归类后的名称以占比最大的部分命名,例如:

了解完重分类后,我们在Arcgis中实现,首先打开ArcToolbox       Spatial Analyst Tools Reclass  Reclassify进行重分类。

点击ok后,土壤重分类完成,重分类成五大类。

将重分类后的value和name设定为.txt形式,txt形式的保存名字必须为英文,内容如下:(建模时用)。

Logo

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

更多推荐