DEM水文分析 提取水系

    #-*-coding:UTF-8-*-#
    import arcpy
    from arcpy.sa import *
    path=unicode(r"D:\PycharmProjects","utf-8")
    # 创建文件地理数据
    mygdb=arcpy.CreateFileGDB_management(path,"mygdb")
    # 设置工作空间
    arcpy.env.workspace=path+"/mygdb.gdb"
    # 输出路径
    outpath=path+"/mygdb.gdb"
    # 输入数据 DEM
    input_dem=unicode(r"D:\PycharmProjects\riverResult.gdb\elevation","utf-8")
    
    #检查SA扩展模块
    if arcpy.CheckExtension("Spatial") == "Available":
    arcpy.CheckOutExtension("Spatial")
    print "OK"
    
    # 01 填挖
    outFill = arcpy.sa.Fill(input_dem)
    outFill.sav