DEM水文分析_提取水系

#-*-coding:UTF-8-*-# import arcpy from arcpy.sa import * path=unicode(r"D:\PycharmProjects","utf-8") # 創(chuàng)建文件地理數(shù)據(jù) mygdb=arcpy.CreateFileGDB_management(path,"mygdb") # 設置工作空間 arcpy.env.workspace=path+"/mygdb.gdb" # 輸出路徑 outpath=path+"/mygdb.gdb" # 輸入數(shù)據(jù) 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.save(outpath+"/outfill")print u"01填挖方計算完成"# 02計算流向outdirection=arcpy.sa.FlowDirection(outFill,"FORCE") #,flow_direction_type="D8")outdirection.save(outpath+"/outdir")print u"02流向計算完成"# 03計算流量outAcc=arcpy.sa.FlowAccumulation(outdirection)outAcc.save(outpath+"/outAcc")print u"03累積流量計算完成"# 04重分類 非常重要outreclass=arcpy.sa.Reclassify(outAcc,"Value",RemapRange([[0,1000,0],[1001,306993,1]]))outreclass.save(outpath+"/outReclass")print u"04重分類計算完成"# 05柵格轉線backgrVal = "ZERO"dangleTolerance = 50arcpy.RasterToPolyline_conversion(outreclass, "flowstream", backgrVal,dangleTolerance, "SIMPLIFY", "Value")print u"05柵格轉線計算完成" else:print u"Spatial Analysis擴展模塊未啟動"

總結

以上是生活随笔為你收集整理的DEM水文分析_提取水系的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內容還不錯，歡迎將生活随笔推薦給好友。