Global DEM - SRTM數據的整備(以北京地區為例)?srtm/ver2

1、將srtm數據轉為Arcinfo DEM

批處理 bat.aml：

&run srtmgrid.aml N39E114.hgt N39E114

&run srtmgrid.aml N39E115.hgt N39E115

&run srtmgrid.aml N39E116.hgt N39E116

&run srtmgrid.aml N39E117.hgt N39E117

&run srtmgrid.aml N40E114.hgt N40E114

&run srtmgrid.aml N40E115.hgt N40E115

&run srtmgrid.aml N40E116.hgt N40E116

&run srtmgrid.aml N40E117.hgt N40E117

&run srtmgrid.aml N41E114.hgt N41E114

&run srtmgrid.aml N41E115.hgt N41E115

&run srtmgrid.aml N41E116.hgt N41E116

&run srtmgrid.aml N41E117.hgt N41E117

&run srtmgrid.aml N42E114.hgt N42E114

&run srtmgrid.aml N42E115.hgt N42E115

&run srtmgrid.aml N42E116.hgt N42E116

&run srtmgrid.aml N42E117.hgt N42E117

Arc: &run bat.aml

Arc: grid

Grid：bjdem =

mosaic(n39e114,n39e115,n39e116,n39e117,n40e114,n40e115,n40e116,n40e117,n41e114,n41e115,n41e116,n41e117,n42e114,n42e115,n42e116,n42e117)

Grid: dem = con(isnull(bjdem),0,bjdem)?/*

填平小黑點:bjdem上有很多小黑點，雷達的盲區，無值區

Grid: fill dem demf?/*

填洼

Grid: water = con(slope(demf) == 0 and demf ==

143,1,0)?/* 提取出密云水庫

#去掉水庫周圍散點：轉成矢量?編輯?再轉回來

Grid: watercov = gridpoly(water)

/* arcview：在extensions 勾選 spatial analyst?選擇三塊大區域

theme: convert to grid?彈出的對話框中選擇 same as

water?grid-code

得到waterbody

Grid: demm = con(isnull(waterbody),demf)

Grid: dir = flowdirection(demm)

Grid: accum = flowaccumulation(dir)

Grid: stm = con(accum > 2000,1)

Grid: stmlnk = streamlink(stm,dir)

Grid: accummax = zonalmax(stmlnk,accum)

Grid: outlet = con(accummax == accum,stmlnk)

Grid: gridclip outlet outletp box 116.76 40.43 117.07

40.58?/*

水庫大致范圍可在arcview中查看

Grid: wtsd1 =

watershed(dir,outletp)?/* 兩種方法

Grid: wtsd2 = watershed(dir,waterbody)

Grid: dirp = con(wtsd1 > 0,dir)

Grid: accump = con(wtsd1 > 0,accum)

Grid: wtsdp = watershed(dirp,stmlnk)

Grid: wtsdpp = con(wtsd1 > 0,wtsdp)

bjdem

dem

demf

watercov

waterbody

demm

wtsd1

wtsd2

wtsdp

wtsdpp

總結

以上是生活随笔為你收集整理的gis怎么通过水库划分子流域_ArcGIS提取小流域——以北京为例，提取出密云水库...的全部內容，希望文章能夠幫你解決所遇到的問題。

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