Arcpy实现dem河网提取

#!/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