SWAT| 水文 | SWAT模型（一）：DEM数据和掩膜数据制备

Tips：
        本期向大家分享SWAT模型的DEM数据和掩膜数据的制备方法，后续还会分享SWAT的土地利用数据库、土壤数据库、气象数据库以及天气发生器的制备方法。打算把SWAT做成一个系列，以后再用到的时候能系统的及时回忆起来，也为没有用过SWAT模型的朋友提供一点灵感。

        在开始之前，要注意的是：矢量文件的裁剪、合并，在栅格文件中对应的是按掩膜提取、镶嵌至新栅格，在操作时不要搞错了。操作得到的DEM数据和掩膜数据都必须进行投影，保证投影坐标系一致，这是后续处理最基本的一环。

---

目录

1 研究区矢量范围提取

2 DEM数据制备

2.1 涉及的地市shp裁剪

2.2 DEM数据镶嵌至新栅格

2.3 DEM数据掩膜提取

2.4 根据地理坐标系投影

3 研究区掩膜数据制备

3.1 研究区水库坐标确定和导入

3.2 流域范围描绘

3.3 面转栅格

4 常见问题

---

1 研究区矢量范围提取

将“中国二级流域”矢量文件加载至ArcMap，右键在属性表中选中研究区：澜沧江流域。将选中的要素导出并添加至图层，得到澜沧江流域矢量范围。

①加载“中国二级流域”矢量文件。

②选择研究区。

③将研究区导出并添加至图层。

④得到澜沧江流域shp，提取完成。

2 DEM数据制备

DEM数据选用哥白尼30m分辨率数字高程数据。制备思路：先用上述得到的澜沧江流域shp裁剪全国行政区划地市shp，得到澜沧江流域的涉及地市shp；随后加载对应地市的DEM数据，将其合并；然后将合并后的DEM数据按澜沧江流域shp或涉及地市shp进行掩膜提取，得到澜沧江流域的DEM数据；最后根据其地理坐标系进行投影，DEM数据制备完成。

2.1 涉及的地市shp裁剪

将全国行政区划的地市shp加载至ArcMap，利用上述澜沧江流域shp裁剪地市shp得到澜沧江流域的涉及地市shp。

①加载全国行政区划的地市shp。

②在ArcGIS工具栏中找到“地理处理”——“裁剪”，将地市shp进行裁剪得到涉及地市shp。

③将涉及地市shp添加至ArcMap。

2.2 DEM数据镶嵌至新栅格

加载对应地市的DEM数据，在ArcToolbox中找到“数据管理工具”——“栅格”——“栅格数据集”——“镶嵌至新栅格”，，将对应地市DEM数据进行镶嵌，镶嵌后的结果如下。（由于这一步操作十分简单，所以不放过程直接放结果了）

2.3 DEM数据掩膜提取

将合并后的DEM数据按澜沧江流域shp或涉及地市shp进行掩膜提取，得到澜沧江流域的DEM数据。

①加载涉及地市shp。

②在ArcToolbox中找到“Spatial Analyst 工具”——“提取分析”——“按掩膜提取”，对合并的DEM数据按掩膜提取。

③将按掩膜提取后的澜沧江DEM数据添加至ArcMap。

2.4 根据地理坐标系投影

右键查看澜沧江DEM数据的属性，能够看到其经度范围是93.86°--101.83°，中央经度取97°，地理坐标系是GCS_WGS_1984。采用UTM投影，UTM投影采用6度分带，根据公式计算，带数=（经度整数位/6）的整数部分+31。将97代入得结果为47，所以应该选择WGS_1984_UTM_Zone_47N作为投影坐标系。

①查看澜沧江DEM数据属性。

​

②在ArcToolbox中找到“数据管理工具”——“投影和变换”——“栅格”——“投影栅格”，对澜沧江DEM数据进行UTM投影，同时把输出的像元大小从30调整为90，这样便得到90m分辨率的DEM数据，减少后续计算量。

​

③将投影后的澜沧江DEM数据添加至ArcMap，制备完成。

​

3 研究区掩膜数据制备

制备掩膜数据是为了更精确的缩小研究区范围。比如，我只研究某水库以上的流域，那么我就可以通过掩膜数据描绘范围来实现。制备思路：导入水库坐标，并转为shp格式；创建面要素，绘制研究区流域范围；完成面要素编辑，面转栅格，掩膜数据制备完成。

3.1 研究区水库坐标确定和导入

①在“天地图”平台查得某水库坐标为东经100.77°，北纬22.05°，将坐标输入Excel中，随后另存为转为CSV格式。

​

②将上述CSV文件添加至ArcMap，右键选择“显示XY数据”。在X字段选择经度，Y字段选择纬度。点击“编辑”按钮选择GCS_WGS_1984地理坐标系。（在ArcIS中，X代表纵轴，Y代表横轴，由于输入的数据是地理坐标，所以只需要选择地理坐标系即可，不要投影！）

​

③在ArcMap中看到XY数据成功显示。

​

3.2 流域范围描绘

①在目录栏中新建文件夹“Mask数据制备”用于存储掩膜数据，右键该文件夹创建Shapefile。

​

②在创建选项中，将Shapefile命名为Mask，要素类型选择“面”，点击“编辑”按钮选择投影坐标系。

​

③右键选择“编辑要素”——“开始编辑”。

​

④描绘某水库上游流域范围。

​

⑤描绘的范围是画个大概，一定要比实际的研究区范围大，保证将流域范围都包含进去，将水电站位置也要包含进去。

​

3.3 面转栅格

①在ArcToolbox中找到“转换工具”——“转为栅格”——“栅格”——“面转栅格”，将shp矢量格式转换为grid栅格格式。注意像元大小不要大于DEM数据的分辨率，可以小于或等于90。

​

②将掩膜数据栅格文件添加至ArcMap，制备完成。

​

4 常见问题

两个数据文件坐标系一致却无法叠加显示？
        这是因为错误地直接将不同的地理坐标系进行投影，应该先统一地理坐标系，再去保证投影坐标系一致，不能直接投影。对于错误投影的文件在目录栏右键文件，点击”属性“，选择”XY坐标系“，鼠标放在小地球图表，点击清除，然后在工具箱中先进行”定义投影“保证地理坐标系一致，然后进行”投影“，保证投影坐标系一致，操作完成后两个数据文件成功叠加显示。

        以上是全部内容，欢迎大家评论区留言，批评指正。