Tips1:土地利用/土壤类型以全国数据分类
鉴于有些步骤我老是忘记处理,我在这里总结如下:
-
1、新建文件夹后,立即将并行设置修改为0,否则会报错
打开Arcgis,点击Geoprocessing——Environments,进入Environment Settings面板,找到Parallel Processing, 输入0,点击OK。 -
2、实测水文站数据记得添加,便于后续径流率定检验
【Outlet and Inlet Definition】中增加点源,数据格式如下:
-
3、土壤数据重分类之前,将新建的土壤类别添加至【SWAT2012.mdb数据库】中的【usersoil】中,土壤数据处理步骤参见此博客
-
4、气象数据导入前,将实测气象站数据导入至【SWAT2012.mdb数据库】中的【WGEN_user表】中
1 土壤类型处理
直接对全国土壤类型数据进行分类,在SWAT模型中输入全国土壤类型,亦可得到研究区域内土壤类型(重分类后),如此可避免SWAT模型反复的输入和处理。但需要提前将研究区域内的土壤重分类数据的物理及化学参数搞清楚,一次操作完成后,后续可直接使用。
具体实现步骤如下:
Step1 从全国土壤数据库中裁剪得到研究区域内土壤类型
将全国土壤类型和研究流域边界图导入,如下:
根据裁剪/掩膜提出研究区域内土壤类型分布图,如下:
确定后,顺利裁剪出区域内土壤类型,如下:
将研究区域内的土壤类型设置为【唯一值】(右键投影后的图像,打开图层属性,选择用唯一值显示),如下:
设置为唯一值的土壤类型如下:
此时,VALUE 值为后面需要计算的各土壤类型唯一编号。
中国土壤数据库结果如下:
世界土壤数据库HWSD结果如下:(84种土壤类型)
Step2 确定研究区域内重分类标准及类型
从以上汇总表可以看出,土壤的种类有几十种,如果不进行重分类,后面的操作过于麻烦,因此考虑根据土壤分组进行重分类。
为什么需要对土壤数据进行重分类?
下载的数据中土壤种类通常很多,数据量大,不便于后续的处理和计算,因此一般先进行重分类,将一些土壤合并,减少种类。
考虑原则如下:同一土壤分组下的所有土壤归为一类,将该土壤分组中百分比最高的这类土作为最终土壤,该组中其他土壤都采用这类土的属性。
首先,从区域土壤数据分布图中属性表复制得到【Rowid VALUE COUNT】三列信息;根据VALUE与MU_GLOBAL一一对应,筛选得到ID和SU_SYM90;从而找到所属土组。
整理得到需要保留的土壤类型,如下:
根据usersoil整理得到所有数据参量,如下:(具体详细步骤参见土壤处理相关博客)
将以上数据添加到SWAT2012.mbd数据库中。
【另】手动重分类步骤如下:
需要用到的工具箱是reclass by table,工具箱路径为:arctoolbox–>spatial analyst tools(空间分析工具)–>reclass(重分类)–>reclass by table(通过表格重分类)。如下:
Step3 ArcSWAT中土壤数据载入
首先,导入全国土壤数据。
其中,input remap table就是你要导入的重分类的表格,最好是csv文件,表格格式如下:
点击【Reclassify】,即可完成土壤数据的重分类。