基于DEM对洪水预报模型的改进初探.pdfVIP

汇流 图l基于DEM的水文信息提取流程图 chart lnfonllationdenVedn_omDEM Flg．I同owofllydmlo酆cal 1．1DEM预处理 由于DEM水平分辨率和垂向分辨率的限制，以及DEM生成过程中的系统误差，不可 避免存在洼地(sink)，造成水流方向不正确，从而给水流线的跟踪和流域边界线的确定带 来困难，以及其它特征信息的不准确。因此，必须进行洼地填充(Fill)。但是，不是所有的 洼地都是由于数据误差造成的，有些洼地是地表形态的真实反映。因此，在进行洼地填充之 li 前，必须计算洼地深度，设置合理的填充阈值(也即Zlimit)。zmjt其数值表示了小于该深 度的洼地才‘会被填平，深度超过此数值的洼地将被视为合理存在狮不被处理。洼地填充过程 是一个反复的过程，直到旧的洼地被填平，新的洼地不再产生才‘能终止。zljmit确定步骤如 下：(1)Flow direction：输入原始的DEM，计算流向。(2)sink：通过流向识别出所有的洼地。 (3)watershed．找出洼地的影响区域。(4)zonalfill：计算出每个洼地区域出口的墩低高程， 也即每个洼地贡献区域的最高高程。也可用spatialanalysttools：[具中的MapA19ebra工具 集直接输入命令来计算。(5)zonalsta石stics：计算出每个洼地贡献区域的最低高程。同理也 可用MapAlgebra来完成。(6)Rastercalculator：计算洼地深度。洼地深度等于最高高程减 去最低高程。 1．2水流方向(Fl洲direction) 要在栅格DEM上提取流域信息，栅格间流向的判别是基础。目前，确定流向的算法主 FlowDirection， Flow Direction，MFD) 要有单流向算法(simple sFD)和多流向算法(Mulciple 两种，因单流l门算法简单方便f『11：『得到广泛的应用。 单流向法假定一个栅格中的水流只从一个方向流出栅格，然后根据栅格高程判断水流方 向。关于单流向的确定，也出现了很多方法，目前应用最为广泛的是D8法，该法假设单个 网格中的水流只有八种可能的流向，即流入与之相邻的八个网格当中，用最陡坡度法来确定 水流方向，I：i口在3*3的窗口上，计算中心格网与各相邻格网间的距离权落差，取距离权落差 最大的格网为中心格网的流出格网，该方向即为中心格网的流向。被处理格网单元同相邻8 个格网单元之间坡降的算法为：s=△z，D，△z为两个格网单元之间的高差，D为两个格网 单元中心之问的距离。在流向矩阵中用8个数字l，2，4，8，16，32，64，128分别代表 方向东，东南，南，西南，西，西北，北，东北8个流向。 1．3汇流累积量(FI佣accumuIation) 汇流累积量表示区域地形每点的流水累积量，假设以规则格网表示的数字高程模型每点 处有一个币．位的水量，按照自然水流从高处流往低处的自然规律，根据区域地形的水流方向 数据计算每点处所流过的水量数值，便得到了该区域的汇流累积量。因此，一个栅格的汇流 累积量数值越大，表示水流能够流入其中的栅格数目越多，其汇流能力也就越强，该栅格所 代表的地形特征就有可能是河谷；反之，汇流累积量数值为零的地方则可能代表流域的分水 岭。 1．4栅格河网生成(singIeoutputmapaIgebra) 河网的生成是通过设置汇流能力阈值来确定的，大于该阈值的栅格被看作是沟谷线上的 点，连接各沟谷线上的点就形成了河流网络。利用c0N语句进行有条件的查询可得到栅格 Ja 河网。例如：con(flowaccmucionJ000，1)其中阈值大于looO，其属性值赋为1，其余设置为 no calculacor也可完成此 data。或者利用arctoolbox中的spatialanalyst分析模块下的raster 功能。至此，栅格河网就提取出来了，然后可以