基于Arc／Info的洪水淹没面积的计算方法.docVIP

基于Arc／Info的洪水淹没面积的计算方法 测绘通报2003年第6期 文章编号:0494.0911(20o3)O6.0046.03中图分类号:P208文献标识码:B 基于Arc/lnfo的洪水淹没面积的计算方法 刘小生一,黄玉生 (1.中国矿业大学,江苏徐州221008;2.南方冶金学院,江西赣州341000) TheCalculationoftheDrownAreaBasedonArc/Info LIUXiao-sheng,HUANGYu-sheng 摘要:以Arc/GIS技术为基础,研究了无需编程就可完成对洪水淹没面积的提取及其计算方法. 关键词:Arc/Info;洪水;淹没面积;计算方法 一 ,引言 对于洪水淹没面积的计算,目前国内使用的是 种子蔓延算法.它的核心思想是将给定的种子点作 为一个对象,赋予特定的属性,在某一平面区域上沿 4个(或8个)方向游动扩散,求取满足给定条件,符 合数据采集精度,且具有连通关联分布的点的集合. 利用种子蔓延算法计算淹没区,就是按给定水位条 件,求取满足精度,且要求连通性的点的集合,该集 合给出的连续平面就是所求的淹没范围.这种种子 蔓延算法就是编程实现求算淹没范围的基本思想, 也是国内专家求算淹没区域的基本方法.这种方法 的不足之处在于必须编写程序,并且专家对于程序 代码都采取保密的态度.本文介绍的是利用Ar. cMaD的自带功能,无需编程就可从低于一特定高程 的区域中确定真正的淹没区,并可以直接计算出它 的面积.下面介绍这种方法. 二,TIN的建立过程 由矢量数据建立TIN,首先要做的是对可能淹 没的区域的地形图进行扫描数字化,采集具备高程 值的大量3D点,由河流和道路等采集隔断线 (Breaklines),有湖泊的边界线采集多边形(polv. gon),而后分两步建立TIN:①单击3DAnalyst的 下拉菜单,找到Create/ModifyTIN项的Create TINFromFeatures.在弹出的对话框里,高程来 源(HeightSource)项选择点的高程.这样建成初步 TIN,还不能完整表示地形,必须往TIN添加要素 以使TIN更加完善.②单击3DAnalyst的下拉 菜单,找到Create/ModifyTIN项的AddFeatures toTIN.在弹出的对话框里,选择作为Breaklines 的数据层,并作为软隔断线或硬隔断线;选择p0lv. gon层作为裁切多边形,删除多边形,替换多边形, 或填充多边形四类中的一类. 另外,也可先由等高线来生成TIN,但这样计 算量过大,由于硬件的缘故,我选择由3D点建成 TIN,不但计算较快,对硬件的要求也不会那么高. 三,把ⅡN转化成Raster 在建成TIN后,把TIN转化成Raster便于取 出低于某个高程值的区域,在ArcMap中,有工具可 供把TIN转化成Raster.具体做法是单击3DAn. alyst,找到Convert的TINtoRaster项,生成我 们所需要的Raster,其值表示为高程值. 四,在Raster中取出低于一定高程值的 区域 利用ArcMap中的SpatialAnalyst菜单中的栅 格计算器(RasterCalculator)功能. 下面介绍在栅格计算器中,用来计算出栅格的 函数一Con函数及相关的方法. Con函数:Performsaconditionalif/elseevalua. tion. 语法(Syntax): Setvariable=object.Con(conditionalRaster, trueRasterl,falseRaster]) 参数及其说明: 参数说明 object对象 variable执行IGeoDataset的对象的一个参数. 收稿日期:2003—02—12 基金项目:江西省主要学科学术和技术带头人培养计划资助项目(2o0l2o9008o0) 作者简介:刘小生(1963一),男,江西新余人,教授,博士,主要从事工程测量与GIS的教学和科研工作. 2003年第6期测绘通报47 conditionalRaster必需的.一个[GeoDataset 对象. trueRaster必需的.一个IGeoDataset对象. falseRaster可选择的.表述falseRaster的变量. 介绍: conditionalRaster输入的栅格.其值包含0和 1,1表述为选中状态,0为未选中状态. trueRaster如果conditionalRaster是真,其值 为选中状态.其值被当成导出的值. falseRaster如果conditionalRaster是假,其值 被当成导出的值,如果所有的conditionalRaster的值 为0,conditionalRast