*****************************************************************
作业一 (复习流域提取)
---------------------------------------------------
OUDOM村所在河流流域的水系、小流域的提取
数据: oddem05
---------------------------------------------------
******************************************************************
grid
/*填洼
fill oddem05 dem
/*计算水流方向
dir = flowdirection(dem)
/*盆地提取
/*（获得多个流域，且流域出水口均在DEM边界或者内部洼地中心）
/*basn就是用户需要的结果
basn = basin(dir)
/*汇流计算
accum = FLOWACCUMULATION(dir)
/*提取水系（河流），河流编码都为1
stmgrd = con(accum > 4000,1)
/*河段（Link）编码，被赋予唯一ID号
stmlnk = streamlink(stmgrd,dir)
/*提取河段下流出水口
stmlnk_acumax = zonalmax(stmlnk,accum)
outlet = con(stmlnk_acumax == accum,stmlnk)
/*提取小流域 [注意，watershed函数的outlet图层也可以用stmlnk替代]
subwatershed = watershed(dir,outlet)
/*获得河段矢量格式图层（每个河段编码具有唯一ID值）
stmcov = gridline(stmlnk)
/*河流分级编码（stereler, shreve）
stmord = streamorder(stmgrd,dir)
/*获得河流矢量格式图层（多个河流可以具有同样的分级值）
stmcovln = streamline(stmord,dir)
/*退出GRID模块
q
*************************************************************************
作业二
---------------------------------------------------
OUDOM村所在河流流域
小流域参数的提取，包括：
（1）小流域的最大坡长、平均坡度
（2）河道的长度与平均坡度
数据: oddem05以及作业一的结果
---------------------------------------------------
提示：学习grid模块的flowlength用法
*************************************************************************
items basn.vat
list basn.vat
/*在ArcView下显示的时候，File Extensions Spatial Analyst
/*Legend Type:Unique Value
/*Values Field:Value
grid
/*提取编号为15的流域
mask = con(basn == 15, 1)
/*将流域15的小流域加到提出来的15流域上
watershd = con(mask == 1, subwatershed)
/*求坡度
slp =
slope(dem,percentrise)
/*dem为<grid>,百分数/度（这个依量纲定）
/*
odslp = zonalmean(watershd,slp)
/*watershd为感兴趣区的小流域，slp为坡度value图层
stmgrd0 = con(accum >
4000,1,0) /*河道，流经的地区为1，不流经的地区为0
oddir = con(mask == 1,
dir)
/*输出流向
oddir0 = con(stmgrd0 == 0,
oddir)/*除了感兴趣区的流向
/*小流域的ID号=河段ID号，河段是贯穿小流域的，有一个进水口，一个出水口，有方向
odflwlen = flowlength(oddir0)
dflw = zonalmax(watershd,odflwlen)
odlnkslp = zonalmean(stmlnk,slp)
odlnkflw = con(stmlnk > 0,
stmlnk.count)
odinfo = ZONALSTATS(watershd,odflwlen)
q
infodbase odinfo
oddbf.dbf
/*生成.dbf文件