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