流域提取步骤.docx

流域提取步骤——以沱沱河为例

WatershedwithDEMdata—AcaseofTuoTuoRiver

1、找到流域坐标范围及水文站（出水口）坐标；（查文献或水文站资料）

沱沱河流域范围：33°40~34°57N，91°00~93°00E

沱沱河水文站坐标：34°135.725N，92°2814.798E

2、去地理空间数据云→用流域坐标范围下载DEM数据；

选择SRTMDEMUTM（投影好的）数据库

90m精度的为5°×5°一个网格，30m精度的为3°×3°一个网格，区域范围一个网格以内的可以直接用下载数据，两个及以上网格的需要在ArcGIS中做影像拼接。

“筛选”中填写坐标范围的中心经纬度，按90m精度下载数据，沱沱河所处范围在30°~35°N，90°~95°E，为一个网格内，故中心经纬度填写92.5°和32.5°，下载筛选出的DEM数据。

3、解压下载的DEM数据到文件夹，文件名自取，比如“WStth”；

4、打开ArcMap，点“+”，加载下载好的DEM数据

5、在右侧Catalog中右击WStth文件夹，新建shapefile，新建一个图层裁取流域范围

要素类型：面Polygon，Name：Shapefile

点“编辑Edit”，“Import导入”底图所在坐标

点“应用”→“确定”→OK

6、选择GotoXY，在新建的图层上用具体流域坐标范围标出四个流域顶点；范围比流域坐标范围稍大，如：

90°30E33°00N90°30E35°00N93°30E35°00N93°30E33°00N

Createfeatures，添加标注点

点“编辑器Edit”开始编辑此图层

画出流域的大概区域，然后选择“保存编辑SaveEdits→停止编辑StopEditing”

去ArcToolbox下的

SpatialAnalystTools空间分析工具→Extraction提取→Extractbymask按掩膜提取流域图

Inputraster

提取完毕后，内容列表中出现裁剪后的流域的图层extent

此时可将不需要的部分去掉√，或者直接删除，得到：

7、用水文站（出水口）具体坐标添加站点到流域图中；

纬度→Y坐标（带小数点的）

经度→X坐标（带小数点的）

XY坐标粘入新的Excel，格式如下

用ArcGIS导入坐标点

右击内容列表下Sheet1$，选择显示XY数据displayXYdata，然后选择坐标参考，将两图层坐标统一。

选择Edit→Select

→GeographicCoordinateSystems

→World→WGS1984→应用→确认→OK

出现水文站的点

右击“Sheet1$Events”，选择Data→ExportData到当前文件夹，saveastype：shapefile

生成水文站点station后，可将不需要的删除

双击DataManagementTools→“ProjectionsandTransformations”→“Feature”→“Project”，将水文站点重投影到平面上：

出现了投影后的水文站点station_utm，station的点可以勾除

8、双击“Hydrology”→“Fill”，Inputsurfaceraster输入裁剪好的流域区域数据，填充区域中表面洼地。

9、双击“Hydrology”→“FlowDirection”，分析地表径流方向，输入数据为填充好的流域DEM数据；

9、用水流方向分析流域中的所有河流

双击“Hydrology”→“FlowAccumulation”，根据径流方向分析得到流域内的所有河流：

10、生成水文站范围的流域图。

先放大看水文站点是否在河道上，如不在则提取不出流域图，需先Edit站点位置到河道上。

双击“Hydrology”→“Watershed”，根据水流方向和水文站点投影坐标提取流域

流域面积即为COUNT×精度2×10-6km2=2177399×902×10-6km2=17637km2。