SWAT水文模型建立及应用:土地利用数据的准备
- 1 简介
- 2 土地利用数据的下载
- 1.1 数据下载方式
- 1.2 数据下载
- 2 土地利用数据的准备
- 2.1 矢量转栅格
- 2.2 土地利用类型的重分类
- 2.3 土地利用分布图投影调整
- 2.4 土地利用类型索引表建立
- 参考
SWAT水文模型建立及应用第一期主要介绍了DEM数据的利用及子流域的划分,具体可参见另一博客-【SWAT水文模型】SWAT水文模型建立及应用第一期:DEM数据利用/子流域的划分。
本博客主要介绍土地利用数据的准备,分为数据下载和数据处理两部分。
1 简介
土地利用是SWAT模型建立的重要数据之一,影响着降水在陆面的成流过程,对模拟结果有重要影响。
SWAT模型需要的土地利用数据包括土地利用分布图(矢量或栅格)及土地利用类型索引表。土地利用分布图的属性数据中必须含有说明图层中土地利用类型的字段,并且每种类型与SWAT landcover/plant数据库中的某条记录相对应。
SWAT模型在加载土地利用类型时可以选择矢量图也可以选择栅格图。
土地利用类型索引表是连接土地利用类型栅格图Value值与SWAT landcover/plant数据库中已有分类的纽带。在建立模型时,用户可以在ArcSWAT界面输入各种土地利用类型与数据记录的对应关系,也可以导入事先准备好的土地利用类型索引表将两者进行关联。
关键步骤 :
- 土地利用分布图矢量转栅格
- 土地利用类型重分类 (可选)
- 土地利用分布图投影调整
- 建立土地利用类型索引表
2 土地利用数据的下载
1.1 数据下载方式
下载地址如下:
-
方式1:中科院1km土地利用数据
空间分辨率:1km
分类等级:含二级分类
-
方式2:清华大学高精度土地利用数据
1.2 数据下载
下载后的LUCC数据如下:
在资源管理器中把【info】和叫【lucc2010(年份)】两个文件夹同时选中,然后拖到arcgis界面,土地利用类型在GIS中显示如下:
LUCC分类体系如下(部分):
2 土地利用数据的准备
2.1 矢量转栅格
SAWT模型在加载土地利用类型时可以选择矢量图也可以选择栅格图,为了方便我们常选用提前处理好的栅格图层。
处理步骤如下:
(1)打开ArcToolbox/Conversion Tool/To Raster/Feature to Raster,打开如下图所示对话框:
(2)输入的矢量图层的属性数据中必须要有一个数值型的能区分各土地利用类型的字段Field,定义输出栅格字段。点击OK。
2.2 土地利用类型的重分类
为什么要进行土地利用类型的重分类?
SWAT模型中建议土地利用类型不超过10种,如果拿到的数据类型超过10种,则需要进行重分类,减少模型最后生成的HRU的数量。
下载的土地利用数据类型通常种类过多,不利于后面SWAT的计算与处理,在不影响研究的情况下,可以进行适当的合并,当然,如果研究关注各种土地利用的影响,则另当别论。
具体实现步骤如下:
(1)ArcToolbox/Spatial Analyst Tools/Reclass/Reclassify,出现如下图所示对话框,输入待重分类的栅格图层,选择区分各土地利用类型的字段,在New value下进行重编号,对想要合并的类型重编同样的新号,即可以合并为一类。点击OK。
2.3 土地利用分布图投影调整
2.4 土地利用类型索引表建立
在本次建模过程中,土地利用类型为五类,其土地利用类型索引表如下图所示,其中“VALUE”即图层中代表了各类型的字段编号,“SNAME”即在数据库中的名称简写,该表的作用就是将研究区的土地利用类型与SWAT数据库中的类型进行关联,该表以.txt格式存储。