Tips：
本期向大家分享SWAT模型的土地利用数据库的制备方法。在开始之前，要注意的是：我的土地利用数据是CNLUCC分类，需要和SWAT土地分类对照表进行比对，重新分类。最终要完成的任务有3项，分别是：
（1）重分类后的土地利用类型数据（栅格数据文件）
（2）土地利用类型索引表（txt文件）
（3）作物生长参数数据库（SWAT自带的参数库“SWAT2012“中的”crop“表）
三者关系如下：索引表起到桥梁的作用，SWAT模型通过索引表来连接栅格数据与crop作物数据库
目录
1 研究区土地利用数据提取
将“中国土地利用类型”栅格文件加载至ArcMap，通过上一期内容得到的澜沧江DEM数据进行掩膜提取，得到澜沧江土地利用类型数据。
①加载“中国土地利用”栅格数据文件。
②在ArcToolbox中找到“Spatial Analyst 工具”——“提取分析”——“按掩膜提取”，利用澜沧江DEM数据lcrdem_prj对中国土地利用数据Id2020按掩膜提取。
③将按掩膜提取后的澜沧江土地利用数据LCRlu添加至ArcMap。
④对澜沧江土地利用数据进行投影，直接选择与澜沧江DEM数据相同的投影坐标系。
⑤将投影后的澜沧江土地利用数据lcrlu_prj添加至ArcMap。
2 研究区土地利用数据重分类
2.1 CNLUCC与SWAT对照重分类
下面有两个长表格，分别是：
（1）中国多时期土地利用遥感监测数据分类表（简称：CNLUCC分类表）
（2）二调土地利用与SWAT土地分类对照表（简称：SWAT对照表）
上述澜沧江土地利用数据是按照CNLUCC分类表进行分类的，我们需要对照SWAT对照表进行土地利用类型重分类，具体的操作可以在Excel中进行，下面会详细介绍。
（1）CNLUCC分类表：
（2）SWAT对照表
将澜沧江土地利用数据属性表中的VALUE列数据导出，导入Excel表格中，VALUE值即为CNLUCC分类表中的编号，然后对照SWAT对照表在Excel表格第二列填入SWAT分类代码，第三列填入新的重分类序号，这个重分类序号没有严格要求，只要自己能区分开就好。比如：第一个VALUE值是11，在CNLUCC分类表中是“水田”，“水田”在SWAT对照表中的代码是AGRL，所以Excel表第二列填AGRL，由于这是我分的第一个类，所以我把AGRL对应的重分类序号定为1，在Excel表第三列填1，往后只要属于AGRL的分类，都填1。重复上述过程直至全部完成重分类。
2.2 ArcGIS对土地利用栅格数据重分类
完成上述操作后，返回ArcMap界面。
①在ArcToolbox中找到“Spatial Analyst 工具”——“重分类”——“重分类”，按照上述Excel表格的重分类结果，输入VALUE值对应的新值。
②将重分类后的澜沧江土地利用数据Re_lcrlu添加至ArcMap，至此土地利用重分类完成，完成了第一项任务。
（补充：像这种分类比较少的情况，用这个方法比较合适。但是如果分类很多，那么为了减少工作量，就需要用VLOOKUP函数和连接表的方法来实现重分类，下一期土壤数据库的制备会做介绍。）
3 研究区土地利用索引表制作
①打开先前的Excel表格，后续根据此表制作土地利用索引表 。
②删除第一列，在将剩余两列数据互换，并修改列名为“VALUE”、“NAME” 。
③在Excel工具栏中选择“数据”——“删除重复值”，这样能够保证每个VALUE都是唯一值。
④在Excel工具栏中选择“文件”，另存为逗号分隔的CSV格式，再将文件名后缀改为txt。
——>修改后缀名为txt——>
⑤打开索引表.txt文件，为列名添加英文双引号（eg：VALUE——>"VALUE"）， 土地利用索引表制作完成，完成了第二项任务。
4 crop表的位置查找
crop表是在下载ArcSWAT时自带的，一个SWAT工程包含3个数据库，分别是：工程库（工程名.mdb）、栅格库(RasterStore.mdb)、参数库(SWAT2012.mdb)。crop表在参数库中可以找到，参数库的默认路径是:
电脑下载ArcSWAT的路径\ArcSwat\Databases\SWAT2012.mdb
我们打开SWAT2012.mdb后找到crop表，会发现其中“CPNM”列即为土地利用索引表对应的“NAME”列，SWAT在运行时会先读取澜沧江土地利用栅格数据的VALUE值， 再读取土地利用索引表里对应VALUE值的NAME，最后在crop表中将NAME和CPNM进行匹配，从而获取研究区对应的土地利用参数，所以索引表十分重要。
以上是全部内容，欢迎大家评论区留言，批评指正。