#!/usr/bin/python
#?-*-?coding:?utf-8?-*-
#?author:zhoulong_GISER
#?-*-?coding:?utf-8?-*-
#?blog:https://blog.csdn.net/qq_33356563
‘‘‘
Arcpy?水文分析
步驟：
1．洼地填平，（填洼工具）Fill
2．水流方向計算。Flow?Direction
3．水流積聚計算。Flow?Accumulation
4.使用spatial?analyst中的柵格計算器，提取河網柵格（二值化）（計算大小根據實際而定，這里設置為200）
5．生成河網矢量，（屬性表選擇grid_code?=?1）

‘‘‘
import?arcpy
from?arcpy?import?env
from?arcpy.sa?import?*
def?main（dem_pathout_path）:
????#?1．洼地填平，（填洼工具）Fill
????outFill?=?Fill（dem_path）
????outFill.save（out_path+“/fill_dem“）
????print（‘完成填洼。。。。。。‘）
????#?2．水流方向計算。Flow?Direction
????outFlowDirection?=?FlowDirection（outFill?“NORMAL“）
????outFlowDirection.save（out_path+“/outflowdir01“）
????print（‘完成水流方向計算。。。。‘）
????#?3．水流積聚計算。FlowAccumulation
????outFlowAccumulation?=?FlowAccumulation（outFlowDirection）
????outFlowAccumulation.save（out_path+“/outflowacc01“）
????print（‘完成水流積聚計算。。。。‘）
????#?4.?使用spati