【SWAT水文模型】SWAT水文模型建立及应用第二期：土地利用数据的准备

SWAT水文模型建立及应用第一期主要介绍了DEM数据的利用及子流域的划分，具体可参见另一博客-【SWAT水文模型】SWAT水文模型建立及应用第一期：DEM数据利用/子流域的划分。
本博客主要介绍土地利用数据的准备，分为数据下载和数据处理两部分。

1 简介

土地利用是SWAT模型建立的重要数据之一，影响着降水在陆面的成流过程，对模拟结果有重要影响。
SWAT模型需要的土地利用数据包括土地利用分布图（矢量或栅格）及土地利用类型索引表。

• 土地利用分布图的属性数据中必须含有说明图层中土地利用类型的字段，并且每种类型与SWAT landcover/plant数据库中的某条记录相对应。SWAT模型在加载土地利用类型时可以选择矢量图也可以选择栅格图。
• 土地利用类型索引表是连接土地利用类型栅格图Value值与SWAT landcover/plant数据库中已有分类的纽带。在建立模型时，用户可以在ArcSWAT界面输入各种土地利用类型与数据记录的对应关系，也可以导入事先准备好的土地利用类型索引表将两者进行关联。

关键步骤 ：

• 土地利用分布图矢量转栅格
• 土地利用类型重分类 （可选）
• 土地利用分布图投影调整
• 建立土地利用类型索引表

2 土地利用数据的下载

2.1 数据下载方式

2.1.1 中科院1km土地利用数据

下载地址如下：中科院1km土地利用数据
空间分辨率：1km
分类等级：含二级分类
在数据集（库）目录中选择土地利用遥感监测数据，点击进入，可以看到数据详情介绍，包含有数据的具体分类系统表。

2.1.2 清华大学高精度土地利用数据

下载地址如下：清华大学高精度土地利用数据

2.2 数据下载

下载后的全国LUCC数据如下：

在资源管理器中把【info】和叫【lucc2010（年份）】两个文件夹同时选中，然后拖到arcgis界面，土地利用类型在GIS中显示如下：

LUCC分类体系如下（部分）：

根据研究区域对全球LUCC数据进行掩膜处理（Extract by Mask），得到研究区域内LUCC图，如下：

3 土地利用数据的准备

3.1 矢量转栅格

SAWT模型在加载土地利用类型时可以选择矢量图也可以选择栅格图，为了方便我们常选用提前处理好的栅格图层。
具体处理步骤如下：
（1）打开ArcToolbox/Conversion Tool/To Raster/Feature to Raster，打开如下图所示对话框：

（2）输入的矢量图层的属性数据中必须要有一个数值型的能区分各土地利用类型的字段Field，定义输出栅格字段。点击OK。

3.2 土地利用类型的重分类

为什么要进行土地利用类型的重分类？
SWAT模型中建议土地利用类型不超过10种，如果拿到的数据类型超过10种，则需要进行重分类，减少模型最后生成的HRU的数量。
下载的土地利用数据类型通常种类过多，不利于后面SWAT的计算与处理，在不影响研究的情况下，可以进行适当的合并，当然，如果研究关注各种土地利用的影响，则另当别论。

考虑使用一级类型，即将11和12两种类型合并为一种，21、22、23、24等合并为第二种，依次类推。具体实现步骤如下：
（1）ArcToolbox/Spatial Analyst Tools/Reclass/Reclassify，出现如下图所示对话框，输入待重分类的栅格图层，选择区分各土地利用类型的字段，在New value下进行重编号，对想要合并的类型重编同样的新号，即可以合并为一类。点击OK。

合并后的土地类型如下图所示：

3.3 土地利用分布图投影调整

在SWAT模型中，DEM、土地利用数据和土壤数据需要投影到统一坐标系。

3.4 土地利用类型索引表建立

为什么需要构建索引表？
下载的数据集与SWAT模型中自带的数据库类型不完全一致，因此需要根据重分类后的土地利用情况与SWAT数据库中的类型进行一一对应起来，即建立索引表。

在本次建模过程中，土地利用类型为六类，分别是：

• 耕地（1）
• 林地（2）
• 草地（3）
• 水域（4）
• 城乡、工矿和居民用地（5）
• 未利用土地（6）

土地利用类型索引表如下图所示，其中“VALUE”即图层中代表了各类型的字段编号，“SNAME”即在数据库中的名称简写，该表的作用就是将研究区的土地利用类型与SWAT数据库中的类型进行关联，该表以.txt格式存储。

至此，土地利用数据准备完毕。

参考

1、CSDN博客-SWAT模型教程—土地利用、土壤数据、气象数据的处理