SWAT土壤数据库的建立
我们之前提到土壤数据栅格属性中的数字,不同的数字代表不同的土壤类型,现在我们来对这些土地类型进行处理。找到我们之前下载的文件夹,打开HWSD.mdb文件,文件的打开需要利用office的Access打开。继续对土壤类型进行重分类,继续对SU_SYM90里的土地利用类型,例如LPk为土壤类型的简称。打开HWSD数据库土壤中文名称表格,找到对应的土壤名称的简称,复制下来。土壤数据整理结束后,可以进行重
#需要的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形式的保存名字必须为英文,内容如下:(建模时用)。
更多推荐
所有评论(0)