简介面向洪水灾害评估、城市防洪规划及地理信息系统应用人员这份研究报告系统梳理了基于ArcGIS的洪水淹没分析与三维模拟方法。内容先对比基于水位与基于水量两种分析思路选用无源淹没模型并说明其适用于洪水源易确定、地势相对平坦的区域。在ArcGIS实现层面文档详细介绍了从TIN数据预处理生成数字高程模型到利用Spatial Analyst中的栅格计算器提取低于设定水位的淹没区再到统计淹没面积并建立洪水水位与淹没面积关系公式的完整流程结合某建成区案例展示了按五十年、二百年一遇设防水位划分区域天然防洪能力的做法使分析结果更贴合实际应用。随后演示了使用ArcScene对水位渐进抬升过程进行三维模拟直观呈现不同水位下的淹没范围变化。资源为单个doc格式电子文档压缩包仅11KB内容紧凑、可操作性强适合GIS初、中级学习者对照练习。目前已有330人浏览学习有助于快速掌握洪水淹没分析的核心思路与ArcGIS关键工具操作。1. 洪水淹没分析不只是画一条等水位线遇到过不少用ArcGIS做洪水淹没分析的人第一反应是在DEM上按高度值做符号化分级认为设到某个高程以上就是淹没区。这个做法应付汇报图勉强够用但一旦要算淹没面积、与土地利用叠加统计、导出风险图误差会大得没法交代。原因在于符号化只是把连续地形切成了几段颜色既没处理洼地、也没有考虑水体在平面上的连通关系更谈不上有多少水量、从哪里来、往哪里去。基于ArcGIS的洪水淹没分析与三维模拟核心是把“地形淹没”这件事拆成两步先在二维栅格上用确定性的方法算出淹没范围再把范围叠回DEM生成三维体块。二维部分常见有两条路线——无源淹没按水位抬升有源淹没按水流路径生长三维部分则用TIN或栅格拉伸来表达水位面与地形的交切关系。这篇文章就按“数据准备—无源淹没—有源淹没—三维模拟—验证与排错”的顺序把每一步的原理和ArcGIS里的具体操作讲透适合正在做防洪规划、灾损评估或应急预案的技术人员参考。2. 数据准备与无源淹没分析的基本范式2.1 DEM 预处理必须先于一切分析ArcGIS里的任何淹没分析都以数字高程模型DEM为底座栅格质量直接决定后续计算的可信度。我一般要求DEM至少满足三个条件一是投影坐标系单位最好是米避免用经纬度直接算面积导致结果不可用二是浮点型存储整数型高程在处理缓坡时会丢失精度三是已完成填洼或至少知道洼地的数量和分布。具体做法是在ArcToolbox里走Spatial Analyst Tools → Hydrology → Fill填洼阈值默认是无穷大即把所有汇水洼地都抬平。实际项目里不建议一上来就用默认值而要看DEM的噪声水平。如果原始DEM来自5米或10米格网微小的人工沟渠、路基形成的假洼地会被无差别填平导致淹没范围外扩。这时可以先做Focal Statistics取3×3邻域的中值把孤立噪点压掉再用Fill处理真正的凹陷。import arcpy from arcpy.sa import * arcpy.CheckOutExtension(Spatial) dem rD:\flood\dem_5m.tif dem_filled Fill(dem) dem_filled.save(rD:\flood\dem_filled.tif)代码里Fill函数返回的是填洼后的栅格参数ZLimit可选表示填洼的最大深度限制设为None表示全部填平。实际使用时我常用ZLimit来控制填洼深度比如设定为2米只填掉2米以内的洼地保留真实的地形凹陷。这样后续河道和堰塞体周围的淹没分析更接近自然状态。填洼后的栅格需要检查属性表中的最小值和最大值确认没有负值异常——负高程在沿海地区合法在内陆山区多半是原始数据的问题。2.2 无源淹没水平面法下的面积与边界提取无源淹没指假设水位瞬间抬高到某个高度H凡高程低于H的像元都被淹没不关心水源位置和连通路径。这在蓄滞洪区启用、水库溃坝的极端情景估算里很常用优点是计算简单、物理意义明确。实现上就是一条栅格计算器表达式Con(dem_filled H, 1, 0)得到0和1构成的二值栅格1代表淹没。为了后续统计方便我通常把Con的结果再转成整型并做一次RegionGroup把连通的淹没斑块赋予唯一编号这样每个斑块的面积可以用ZonalGeometry或RasterToPolygon转面后按属性表计算。in_dem rD:\flood\dem_filled.tif water_level 35.5 inundated Con(in_dem water_level, 1, 0) inundated.save(rD:\flood\inundated_355.tif) polygons RasterToPolygon(inundated, simplifyTrue, raster_fieldVALUE) polygons.save(rD:\flood\inundated_355.shp)这段代码里Con函数走的是Spatial Analyst的条件判断in_dem 35.5生成布尔栅格满足条件赋1否则0。RasterToPolygon默认会把值为1的区域转成面simplifyTrue表示对边界做抽稀减少面要素的顶点数量。这里要特别说明无源淹没并没有区分“连不连通”如果一个低洼山谷和一个盆地的高程都低于水位它们都会显示为淹没区如果两个区域被山脊隔开也想合并成同一片洪水体那就得换有源淹没的思路。另外水位H的取值通常来自水文频率计算的洪峰水位或某重现期下的控制断面水位没做水动力模拟的项目里H一般用设计洪水位代替。2.3 无源淹没结果的两个常见误读与修正第一个误读是把淹没栅格的像元数量直接当面积。栅格分辨率如果是10米一个像元代表100平方米统计时用GetRasterProperties取COUNT再乘像元面积才对如果转了矢量面属性表里的Shape_Area字段会直接给出平方米数但转面时空白区也就是NoData区域会被跳过边界处的面积略小于栅格真实面积这在极高精度要求下需要注意。第二个误读是忽略高程基准面不一致。同一个研究区里DEM的高程基准若是1985国家高程基准而水位数据来自吴淞基准两者相差约1.8米左右的常数偏移必须先做统一否则1米级别的淹没范围判断会失真。修正做法在转出淹没面后用Tabulate Intersection叠加土地利用或人口分布数据把淹没范围落到承灾体上。landslide rD:\flood\landuse_2020.shp inundated_poly rD:\flood\inundated_355.shp stats TabulateIntersection(inundated_poly, OBJECTID, landslide, LANDUSE, rD:\flood\inund_landuse.dbf)TabulateIntersection参数第一个是区域要素第二个是区域标识字段第三个是被统计图层第四个是统计字段输出是DBF表。表里每一行对应区域要素的每个被统计类别AREA字段给出相交面积。这样可以量化不同地类里被淹的面积为灾损粗估提供直接输入。3. 有源淹没分析的实现路径与关键参数3.1 有源淹没与无源淹没的本质差异有源淹没考虑的是洪水从某个入口持续流入像元被淹不仅因为高程低还因为有一条从源点到达它的路径。这个路径由地形坡度决定水往低处流但不能爬过高山所以有源淹没天然排除了被山脊隔开的低洼盆地。工程上常见做法以种子蔓延为核心从指定的进水口像元出发向四周八邻域扩展当相邻像元高程低于“起始水位加上沿程水头损失”时纳入淹没区否则停止。严格来说这属于简化后的静水淹没模型没有解圣维南方程但它比无源淹没多了一个“水流可达性”约束用来评估溃口周边范围的快速淹没非常合适。ArcGIS里没有现成的“有源淹没”按钮一般两种实现思路。一种是用Cost Distance或Path Distance加Con组合出累积代价面把“到达某个像元需要翻越的最大地形阻力”算出来再与水位比较另一种是直接用Python配合Raster Calculator做区域生长或者借助Watershed工具给定出水口反推汇流区。前一种思路稳定适合批量水位方案比较后一种思路更接近物理直觉但需要写循环。多数情况下我推荐前者下面给出具体参数和表达式。3.2 用 Cost Distance 实现最小阻力可达性分析先基于填洼后的DEM生成每个像元到源点的最小成本距离这里的成本不是路程而是“途经像元时相对源点需要克服的高程差”。具体做法是用Path Distance工具把Cost raster设置为全1栅格把Surface raster设置为原始DEMVertical factor选择Binary这样高程每升高一点成本就大幅增加更简单的办法是先用Minus(dem_filled, source_level)生成相对高差栅格再把正的高差部分作为成本。我常写这样一个步骤先创建源点栅格进水口位置赋0其他像元为NoData。然后用CostDistance得到累计最小成本表面。这个成本表面的含义是从源点到达每个像元沿途需要翻越的高程代价之和。注意这里不是简单的高程差值因为路径可以选择从山脊两侧绕过去。source_point rD:\flood\source.shp # 点要素表示进水口 source_raster RasterToPoint(source_point) # 需配合提取赋值 # 更可靠的做法是先转成栅格 arcpy.env.extent dem_filled.extent arcpy.env.snapRaster dem_filled source_r arcpy.conversion.PointToRaster(source_point, OBJECTID, rD:\flood\source_raster, MOST_FREQUENT, 10) cost_surface CostDistance(source_r, 1, dem_filled, Binary)CostDistance的第一个参数是源栅格第二个是权重栅格即成本栅格第三个是表面栅格用于考虑实际地表距离第四个Vertical factor用Binary表含义是只要上坡成本即乘以一个极大系数而下坡成本为1。这样生成的成本面数值越大说明从源点到达这里越需要爬升。然后与水位做比较假设源点的起始水位为H0沿程水面坡降近似忽略那么凡cost_surface小于某个阈值的区域就是可淹没区。阈值不是水位值本身而是“允许的最大爬升高度”在平原区往往取0.5-1米代表洪水漫过微小田埂的能力山区取2米以上代表快速上涨时水头对局部鞍部的漫越。{{INSERT}}### 3.3 区域生长法的ArcGIS Python实现如果研究区小、源点数量有限更精细的是区域生长。用Python循环处理每次取当前边界像元的八邻域判断相邻像元高程与当前水位的差值来决定是否纳入。这里要注意区域生长如果写成逐像元循环10米分辨率、10000×10000的栅格会非常慢所以我一般把生长过程向量化维护一个候选队列优先处理最低高程的像元这样避免了大面积重复扫描。import numpy as np from scipy import ndimage dem_arr arcpy.RasterToNumPyArray(dem_filled, nodata_to_value-9999) inund_arr np.zeros_like(dem_arr, dtypenp.int8) queue [(r0, c0)] # 源点行列号 water_level 40.8 while queue: r, c queue.pop(0) for dr in [-1, 0, 1]: for dc in [-1, 0, 1]: nr, nc r dr, c dc if 0 nr dem_arr.shape[0] and 0 nc dem_arr.shape[1]: if inund_arr[nr, nc] 0 and dem_arr[nr, nc] water_level: inund_arr[nr, nc] 1 queue.append((nr, nc))这段逻辑是广度优先搜索每纳入一个新像元就把它加入队列继续向外探。循环条件是水位恒定如果模拟随时间上涨的水位就要把water_level改成随迭代次数增加的变量并且每次增长时重新扫描边界像元。不过要注意Python队列处理大范围区域会慢实际工作流建议先用CostDistance粗筛再用区域生长做局部细化。3.4 参数敏感性水位、源点位置和DEM分辨率的影响有源淹没对源点位置非常敏感。同一个DEM源点设在下游河道和设在溃口处结果差异极大源点所在的像元如果是洼地底部洪水会先填满洼地再外溢生长范围受限源点如果选在山脊上则几乎不生长。我一般在源点生成前用Snap Pour Point工具把点对齐到流向累积量最大的像元保证起始水位出现在河流的谷底线。水位参数的敏感性在不同坡度下差别很大——陡峭山区水位抬高5米水平扩展范围可能只有几十米平原地区水位抬高0.5米淹没范围就成片扩张。做方案比选时建议做成参数栅格图水位H按0.5米间隔生成10组淹没面叠加成一张风险图比单一一套结果更有说服力。4. ArcGIS 三维模拟从淹没面到体块表达4.1 三维场景的两种基础表达栅格拉伸与 TIN二维淹没分析得到的是一个面范围三维模拟要表达的是“水位面以下、地形表面以上的空间体”。ArcGIS里常用ArcScene或ArcGIS Pro的Local Scene实现。第一步是构建地形表面小范围用原始DEM直接作为高程源大范围先转TIN再做Terrain优化。第二步是把淹没区做垂直拉伸拉伸方式有两种基于图层属性拉伸和基于栅格值拉伸。如果淹没范围是从Con表达式里来的二值栅格给0值设透明、1值设置基础高度为水位值符号系统里选“按属性拉伸”拉伸字段设为水位常数更好的做法是先用RasterCalculator生成一个常量栅格高程值等于水位的平面再对这个平面做ExtrudeBetween把地形面和平面之间的空间生成实体。我用Python走一遍完整流程核心是把DEM和水位面都转成TIN或带Z值的多面体。这里可以借助GP工具RasterSurface和ExtrudeBetween但ArcGIS原生没有直接生成水体的工具常用做法是把淹没区的矢量面要素增强为三维要素。arcpy.env.workspace rD:\flood\scene.gdb dem_tin arcpy.ddd.RasterTin(rD:\flood\dem_filled.tif, TIN_DEM) water_surface arcpy.ddd.CreateConstantRaster(water_level, FLOAT, 10, extent) # 把水位面转成多面体 water_tin arcpy.ddd.RasterTin(water_surface, TIN_WATER) # 用ExtrudeBetween生成水位面与地形之间的体块 arcpy.ddd.ExtrudeBetween(dem_tin, water_tin, rD:\flood\flood_body, BETWEEN)RasterTin把栅格转成不规则三角网ExtrudeBetween接受两个TIN输出两者之间的闭合体块。BETWEEN模式生成地形与水位面之间的体积ABOVE和BELOW则分别只保留某一侧的体块。生成后的多面体可以直接在ArcScene中设置透明度做成半透明水体叠在影像上。不过ExtrudeBetween在数据量大时容易产生产状物因为两个TIN的三角形边长不一致交切处会出现细长三角形。常见解法是在转TIN时把z_tolerance设为0.5到1米减少三角形数量交切更干净。4.2 场景符号化与飞行漫游导出三维场景的视觉效果取决于两个设置光照和水体透明度。ArcScene里水体图层的符号系统选择Simple Fill把颜色改为浅蓝色透明度调到40%左右地形表面用带山体阴影的立体符号叠加影像底图。淹没体块的底面和顶面如果是分开的需要做一次Union确保体块闭合否则在场景里旋转视角时会发现内部是空的不利于评审汇报。导出动画用ArcGIS Pro的View → Animation功能设定相机的路径和朝向顺时针绕研究区一圈帧数设为300帧左右导出MP4。4.3 体积计算淹没水量和灾损关联三维体块除了看还能量算体积。ArcGIS的Surface Volume工具可以计算水体体积和表面积。输入是表示水位的栅格表面基准面设为ABOVE结果表给出参考高度以下的体积。这个体积就是淹没水量可以与水文站的洪量数据对照反过来验证淹没范围的合理性。也可以用PolygonVolume直接对三维体块要素计算体积。体积计算结果的单位取决于DEM的投影单位如果DEM是米体积就是立方米如果是英尺就得乘0.028317转换。vol_table rD:\flood\volume.dbf arcpy.ddd.SurfaceVolume(rD:\flood\dem_filled.tif, ABOVE, vol_table, 10, 35.5)SurfaceVolume的四个参数依次是输入表面、方向ABOVE表示计算表面以上到参考平面的体积、输出表、参考高程间隔和参考基准面高程。参考高程35.5米即水位输出表里的Volume字段就是淹没水量。与前面无源淹没水面法对比这个体积值更准确因为它在每个像元上动态计算面积和高差而不只是平面面积乘平均水深。5. 结果验证与参数调试的接地气技巧5.1 与历史洪水痕迹对照的三步验证法计算完的淹没范围不能直接采信要和历史数据对照。常见做法是找水文站记录的“历史最高水位”的痕迹线比如墙上的水痕、桥墩上的标尺记录在ArcGIS里把这些痕迹点或线叠加到淹没范围上看吻合度。验证分三步第一步把淹没范围的边界提取出来转成线要素FeatureToLine叠加痕迹线看偏移距离第二步统计痕迹点中有多少落在淹没范围内用Select By Location的INTERSECT淹中率低于70%说明参数偏保守第三步反查错误把漏淹的点导出来看它们的流向和DEM坡向通常原因是水位偏低或者源点位置不对。5.2 水位步长与栅格分辨率的最佳配合做多水位方案时水位步长和DEM分辨率之间有个经验配比。DEM水平分辨率为10米时水位步长选0.5米产生的面积增量往往在平原区突变明显而山区变化平缓若水位步长小于DEM垂直精度比如DEM垂直误差有0.3米那0.1米步长就是虚假精度算出来的面积差异全是噪声。我通常先查DEM属性表里的标准差取标准差的一半作为最小水位步长。研究区是陡峭山区DEM标准差可能有30米那水位步长取1米甚至2米就够了平原地区标准差只有2米水位步长取0.2米才能反映微地形的淹没差异。分辨率方面如果原始DEM是30米不建议直接重采样成10米来追求精度因为插值出来的假细节会误导淹没边界正确做法是保持30米分析但在出图前做一步Aggregate的平滑只改善视觉效果不改变计算值。5.3 三条常见的失败路径和对应的排除方法图表常见问题与检查顺序现象直接原因排查处理手段淹没范围完全空白坐标系不一致或源点落在NoData区用Project Raster统一到DEM坐标系IsNull检查源点像元无源淹没面积远超预期DEM没有填洼内部盆地全部积水做Fill后重跑Con并对比填洼前后DEM的差异三维体块有破洞水位面与地形面三角形尺寸差异过大调整z_tolerance重建TIN或把水位面重采样为DEM两倍像元大小三个问题都出在“数据级别不匹配”。Coordinated Universal Time的错位、NoData填充值不同、像元对齐不一致是最容易忽略的三个点。展开说Con判断时如果水位值落入NoData区结果是不参与计算的而在ArcScene里显示时NoData默认是透明看起来像被掏空。最实用的检查方式是打开栅格的属性表看COUNT和唯一值分布如果唯一值只有0和1且COUNT远小于整个范围的像元数说明有NoData混入用IsNull做掩膜把NoData填成0或极大值再做分析。5.4 一套可复用的批处理模板最后给出一套可以改参数直接跑的批处理骨架把DEM路径、水位列表、输出目录三者分离适合做多情景淹没分析。水位列表可以是多个断面水位也可以按时间增长的演进序列。这个模板我在多个防洪评价项目里用过替换数据路径后基本不用改逻辑。import arcpy from arcpy.sa import * arcpy.CheckOutExtension(Spatial) dem rD:\flood\dem_filled.tif out_dir rD:\flood\scenarios levels [32.5, 33.0, 33.5, 34.0, 34.5, 35.0] for h in levels: out_name finund_{str(h).replace(., _)} inund Con(dem h, 1, 0) inund.save(fr{out_dir}\{out_name}.tif) poly RasterToPolygon(inund, simplifyTrue) poly.save(fr{out_dir}\{out_name}.shp) print(f水位{h}完成范围面积{poly.getArea()})打印的面积实际需要从结果中读取这里展示的是循环写法多情景的关键是水位列表的生成方式可以用range(320, 361, 5)除以10生成等差序列也可以直接从CSV文件读取水文站的预报水位。批量输出的栅格可以用MosaicToNewRaster合并成多层栅格提供后续制图使用。三维模拟方面把每个水位对应的淹没范围批量拉伸并存入同一个场景文件逐层显示就可以制作出随时间上涨的动画效果比单张静态图清晰得多。本文还有配套的精品资源点击获取
简介面向洪水灾害评估、城市防洪规划及地理信息系统应用人员这份研究报告系统梳理了基于ArcGIS的洪水淹没分析与三维模拟方法。内容先对比基于水位与基于水量两种分析思路选用无源淹没模型并说明其适用于洪水源易确定、地势相对平坦的区域。在ArcGIS实现层面文档详细介绍了从TIN数据预处理生成数字高程模型到利用Spatial Analyst中的栅格计算器提取低于设定水位的淹没区再到统计淹没面积并建立洪水水位与淹没面积关系公式的完整流程结合某建成区案例展示了按五十年、二百年一遇设防水位划分区域天然防洪能力的做法使分析结果更贴合实际应用。随后演示了使用ArcScene对水位渐进抬升过程进行三维模拟直观呈现不同水位下的淹没范围变化。资源为单个doc格式电子文档压缩包仅11KB内容紧凑、可操作性强适合GIS初、中级学习者对照练习。目前已有330人浏览学习有助于快速掌握洪水淹没分析的核心思路与ArcGIS关键工具操作。1. 洪水淹没分析不只是画一条等水位线遇到过不少用ArcGIS做洪水淹没分析的人第一反应是在DEM上按高度值做符号化分级认为设到某个高程以上就是淹没区。这个做法应付汇报图勉强够用但一旦要算淹没面积、与土地利用叠加统计、导出风险图误差会大得没法交代。原因在于符号化只是把连续地形切成了几段颜色既没处理洼地、也没有考虑水体在平面上的连通关系更谈不上有多少水量、从哪里来、往哪里去。基于ArcGIS的洪水淹没分析与三维模拟核心是把“地形淹没”这件事拆成两步先在二维栅格上用确定性的方法算出淹没范围再把范围叠回DEM生成三维体块。二维部分常见有两条路线——无源淹没按水位抬升有源淹没按水流路径生长三维部分则用TIN或栅格拉伸来表达水位面与地形的交切关系。这篇文章就按“数据准备—无源淹没—有源淹没—三维模拟—验证与排错”的顺序把每一步的原理和ArcGIS里的具体操作讲透适合正在做防洪规划、灾损评估或应急预案的技术人员参考。2. 数据准备与无源淹没分析的基本范式2.1 DEM 预处理必须先于一切分析ArcGIS里的任何淹没分析都以数字高程模型DEM为底座栅格质量直接决定后续计算的可信度。我一般要求DEM至少满足三个条件一是投影坐标系单位最好是米避免用经纬度直接算面积导致结果不可用二是浮点型存储整数型高程在处理缓坡时会丢失精度三是已完成填洼或至少知道洼地的数量和分布。具体做法是在ArcToolbox里走Spatial Analyst Tools → Hydrology → Fill填洼阈值默认是无穷大即把所有汇水洼地都抬平。实际项目里不建议一上来就用默认值而要看DEM的噪声水平。如果原始DEM来自5米或10米格网微小的人工沟渠、路基形成的假洼地会被无差别填平导致淹没范围外扩。这时可以先做Focal Statistics取3×3邻域的中值把孤立噪点压掉再用Fill处理真正的凹陷。import arcpy from arcpy.sa import * arcpy.CheckOutExtension(Spatial) dem rD:\flood\dem_5m.tif dem_filled Fill(dem) dem_filled.save(rD:\flood\dem_filled.tif)代码里Fill函数返回的是填洼后的栅格参数ZLimit可选表示填洼的最大深度限制设为None表示全部填平。实际使用时我常用ZLimit来控制填洼深度比如设定为2米只填掉2米以内的洼地保留真实的地形凹陷。这样后续河道和堰塞体周围的淹没分析更接近自然状态。填洼后的栅格需要检查属性表中的最小值和最大值确认没有负值异常——负高程在沿海地区合法在内陆山区多半是原始数据的问题。2.2 无源淹没水平面法下的面积与边界提取无源淹没指假设水位瞬间抬高到某个高度H凡高程低于H的像元都被淹没不关心水源位置和连通路径。这在蓄滞洪区启用、水库溃坝的极端情景估算里很常用优点是计算简单、物理意义明确。实现上就是一条栅格计算器表达式Con(dem_filled H, 1, 0)得到0和1构成的二值栅格1代表淹没。为了后续统计方便我通常把Con的结果再转成整型并做一次RegionGroup把连通的淹没斑块赋予唯一编号这样每个斑块的面积可以用ZonalGeometry或RasterToPolygon转面后按属性表计算。in_dem rD:\flood\dem_filled.tif water_level 35.5 inundated Con(in_dem water_level, 1, 0) inundated.save(rD:\flood\inundated_355.tif) polygons RasterToPolygon(inundated, simplifyTrue, raster_fieldVALUE) polygons.save(rD:\flood\inundated_355.shp)这段代码里Con函数走的是Spatial Analyst的条件判断in_dem 35.5生成布尔栅格满足条件赋1否则0。RasterToPolygon默认会把值为1的区域转成面simplifyTrue表示对边界做抽稀减少面要素的顶点数量。这里要特别说明无源淹没并没有区分“连不连通”如果一个低洼山谷和一个盆地的高程都低于水位它们都会显示为淹没区如果两个区域被山脊隔开也想合并成同一片洪水体那就得换有源淹没的思路。另外水位H的取值通常来自水文频率计算的洪峰水位或某重现期下的控制断面水位没做水动力模拟的项目里H一般用设计洪水位代替。2.3 无源淹没结果的两个常见误读与修正第一个误读是把淹没栅格的像元数量直接当面积。栅格分辨率如果是10米一个像元代表100平方米统计时用GetRasterProperties取COUNT再乘像元面积才对如果转了矢量面属性表里的Shape_Area字段会直接给出平方米数但转面时空白区也就是NoData区域会被跳过边界处的面积略小于栅格真实面积这在极高精度要求下需要注意。第二个误读是忽略高程基准面不一致。同一个研究区里DEM的高程基准若是1985国家高程基准而水位数据来自吴淞基准两者相差约1.8米左右的常数偏移必须先做统一否则1米级别的淹没范围判断会失真。修正做法在转出淹没面后用Tabulate Intersection叠加土地利用或人口分布数据把淹没范围落到承灾体上。landslide rD:\flood\landuse_2020.shp inundated_poly rD:\flood\inundated_355.shp stats TabulateIntersection(inundated_poly, OBJECTID, landslide, LANDUSE, rD:\flood\inund_landuse.dbf)TabulateIntersection参数第一个是区域要素第二个是区域标识字段第三个是被统计图层第四个是统计字段输出是DBF表。表里每一行对应区域要素的每个被统计类别AREA字段给出相交面积。这样可以量化不同地类里被淹的面积为灾损粗估提供直接输入。3. 有源淹没分析的实现路径与关键参数3.1 有源淹没与无源淹没的本质差异有源淹没考虑的是洪水从某个入口持续流入像元被淹不仅因为高程低还因为有一条从源点到达它的路径。这个路径由地形坡度决定水往低处流但不能爬过高山所以有源淹没天然排除了被山脊隔开的低洼盆地。工程上常见做法以种子蔓延为核心从指定的进水口像元出发向四周八邻域扩展当相邻像元高程低于“起始水位加上沿程水头损失”时纳入淹没区否则停止。严格来说这属于简化后的静水淹没模型没有解圣维南方程但它比无源淹没多了一个“水流可达性”约束用来评估溃口周边范围的快速淹没非常合适。ArcGIS里没有现成的“有源淹没”按钮一般两种实现思路。一种是用Cost Distance或Path Distance加Con组合出累积代价面把“到达某个像元需要翻越的最大地形阻力”算出来再与水位比较另一种是直接用Python配合Raster Calculator做区域生长或者借助Watershed工具给定出水口反推汇流区。前一种思路稳定适合批量水位方案比较后一种思路更接近物理直觉但需要写循环。多数情况下我推荐前者下面给出具体参数和表达式。3.2 用 Cost Distance 实现最小阻力可达性分析先基于填洼后的DEM生成每个像元到源点的最小成本距离这里的成本不是路程而是“途经像元时相对源点需要克服的高程差”。具体做法是用Path Distance工具把Cost raster设置为全1栅格把Surface raster设置为原始DEMVertical factor选择Binary这样高程每升高一点成本就大幅增加更简单的办法是先用Minus(dem_filled, source_level)生成相对高差栅格再把正的高差部分作为成本。我常写这样一个步骤先创建源点栅格进水口位置赋0其他像元为NoData。然后用CostDistance得到累计最小成本表面。这个成本表面的含义是从源点到达每个像元沿途需要翻越的高程代价之和。注意这里不是简单的高程差值因为路径可以选择从山脊两侧绕过去。source_point rD:\flood\source.shp # 点要素表示进水口 source_raster RasterToPoint(source_point) # 需配合提取赋值 # 更可靠的做法是先转成栅格 arcpy.env.extent dem_filled.extent arcpy.env.snapRaster dem_filled source_r arcpy.conversion.PointToRaster(source_point, OBJECTID, rD:\flood\source_raster, MOST_FREQUENT, 10) cost_surface CostDistance(source_r, 1, dem_filled, Binary)CostDistance的第一个参数是源栅格第二个是权重栅格即成本栅格第三个是表面栅格用于考虑实际地表距离第四个Vertical factor用Binary表含义是只要上坡成本即乘以一个极大系数而下坡成本为1。这样生成的成本面数值越大说明从源点到达这里越需要爬升。然后与水位做比较假设源点的起始水位为H0沿程水面坡降近似忽略那么凡cost_surface小于某个阈值的区域就是可淹没区。阈值不是水位值本身而是“允许的最大爬升高度”在平原区往往取0.5-1米代表洪水漫过微小田埂的能力山区取2米以上代表快速上涨时水头对局部鞍部的漫越。{{INSERT}}### 3.3 区域生长法的ArcGIS Python实现如果研究区小、源点数量有限更精细的是区域生长。用Python循环处理每次取当前边界像元的八邻域判断相邻像元高程与当前水位的差值来决定是否纳入。这里要注意区域生长如果写成逐像元循环10米分辨率、10000×10000的栅格会非常慢所以我一般把生长过程向量化维护一个候选队列优先处理最低高程的像元这样避免了大面积重复扫描。import numpy as np from scipy import ndimage dem_arr arcpy.RasterToNumPyArray(dem_filled, nodata_to_value-9999) inund_arr np.zeros_like(dem_arr, dtypenp.int8) queue [(r0, c0)] # 源点行列号 water_level 40.8 while queue: r, c queue.pop(0) for dr in [-1, 0, 1]: for dc in [-1, 0, 1]: nr, nc r dr, c dc if 0 nr dem_arr.shape[0] and 0 nc dem_arr.shape[1]: if inund_arr[nr, nc] 0 and dem_arr[nr, nc] water_level: inund_arr[nr, nc] 1 queue.append((nr, nc))这段逻辑是广度优先搜索每纳入一个新像元就把它加入队列继续向外探。循环条件是水位恒定如果模拟随时间上涨的水位就要把water_level改成随迭代次数增加的变量并且每次增长时重新扫描边界像元。不过要注意Python队列处理大范围区域会慢实际工作流建议先用CostDistance粗筛再用区域生长做局部细化。3.4 参数敏感性水位、源点位置和DEM分辨率的影响有源淹没对源点位置非常敏感。同一个DEM源点设在下游河道和设在溃口处结果差异极大源点所在的像元如果是洼地底部洪水会先填满洼地再外溢生长范围受限源点如果选在山脊上则几乎不生长。我一般在源点生成前用Snap Pour Point工具把点对齐到流向累积量最大的像元保证起始水位出现在河流的谷底线。水位参数的敏感性在不同坡度下差别很大——陡峭山区水位抬高5米水平扩展范围可能只有几十米平原地区水位抬高0.5米淹没范围就成片扩张。做方案比选时建议做成参数栅格图水位H按0.5米间隔生成10组淹没面叠加成一张风险图比单一一套结果更有说服力。4. ArcGIS 三维模拟从淹没面到体块表达4.1 三维场景的两种基础表达栅格拉伸与 TIN二维淹没分析得到的是一个面范围三维模拟要表达的是“水位面以下、地形表面以上的空间体”。ArcGIS里常用ArcScene或ArcGIS Pro的Local Scene实现。第一步是构建地形表面小范围用原始DEM直接作为高程源大范围先转TIN再做Terrain优化。第二步是把淹没区做垂直拉伸拉伸方式有两种基于图层属性拉伸和基于栅格值拉伸。如果淹没范围是从Con表达式里来的二值栅格给0值设透明、1值设置基础高度为水位值符号系统里选“按属性拉伸”拉伸字段设为水位常数更好的做法是先用RasterCalculator生成一个常量栅格高程值等于水位的平面再对这个平面做ExtrudeBetween把地形面和平面之间的空间生成实体。我用Python走一遍完整流程核心是把DEM和水位面都转成TIN或带Z值的多面体。这里可以借助GP工具RasterSurface和ExtrudeBetween但ArcGIS原生没有直接生成水体的工具常用做法是把淹没区的矢量面要素增强为三维要素。arcpy.env.workspace rD:\flood\scene.gdb dem_tin arcpy.ddd.RasterTin(rD:\flood\dem_filled.tif, TIN_DEM) water_surface arcpy.ddd.CreateConstantRaster(water_level, FLOAT, 10, extent) # 把水位面转成多面体 water_tin arcpy.ddd.RasterTin(water_surface, TIN_WATER) # 用ExtrudeBetween生成水位面与地形之间的体块 arcpy.ddd.ExtrudeBetween(dem_tin, water_tin, rD:\flood\flood_body, BETWEEN)RasterTin把栅格转成不规则三角网ExtrudeBetween接受两个TIN输出两者之间的闭合体块。BETWEEN模式生成地形与水位面之间的体积ABOVE和BELOW则分别只保留某一侧的体块。生成后的多面体可以直接在ArcScene中设置透明度做成半透明水体叠在影像上。不过ExtrudeBetween在数据量大时容易产生产状物因为两个TIN的三角形边长不一致交切处会出现细长三角形。常见解法是在转TIN时把z_tolerance设为0.5到1米减少三角形数量交切更干净。4.2 场景符号化与飞行漫游导出三维场景的视觉效果取决于两个设置光照和水体透明度。ArcScene里水体图层的符号系统选择Simple Fill把颜色改为浅蓝色透明度调到40%左右地形表面用带山体阴影的立体符号叠加影像底图。淹没体块的底面和顶面如果是分开的需要做一次Union确保体块闭合否则在场景里旋转视角时会发现内部是空的不利于评审汇报。导出动画用ArcGIS Pro的View → Animation功能设定相机的路径和朝向顺时针绕研究区一圈帧数设为300帧左右导出MP4。4.3 体积计算淹没水量和灾损关联三维体块除了看还能量算体积。ArcGIS的Surface Volume工具可以计算水体体积和表面积。输入是表示水位的栅格表面基准面设为ABOVE结果表给出参考高度以下的体积。这个体积就是淹没水量可以与水文站的洪量数据对照反过来验证淹没范围的合理性。也可以用PolygonVolume直接对三维体块要素计算体积。体积计算结果的单位取决于DEM的投影单位如果DEM是米体积就是立方米如果是英尺就得乘0.028317转换。vol_table rD:\flood\volume.dbf arcpy.ddd.SurfaceVolume(rD:\flood\dem_filled.tif, ABOVE, vol_table, 10, 35.5)SurfaceVolume的四个参数依次是输入表面、方向ABOVE表示计算表面以上到参考平面的体积、输出表、参考高程间隔和参考基准面高程。参考高程35.5米即水位输出表里的Volume字段就是淹没水量。与前面无源淹没水面法对比这个体积值更准确因为它在每个像元上动态计算面积和高差而不只是平面面积乘平均水深。5. 结果验证与参数调试的接地气技巧5.1 与历史洪水痕迹对照的三步验证法计算完的淹没范围不能直接采信要和历史数据对照。常见做法是找水文站记录的“历史最高水位”的痕迹线比如墙上的水痕、桥墩上的标尺记录在ArcGIS里把这些痕迹点或线叠加到淹没范围上看吻合度。验证分三步第一步把淹没范围的边界提取出来转成线要素FeatureToLine叠加痕迹线看偏移距离第二步统计痕迹点中有多少落在淹没范围内用Select By Location的INTERSECT淹中率低于70%说明参数偏保守第三步反查错误把漏淹的点导出来看它们的流向和DEM坡向通常原因是水位偏低或者源点位置不对。5.2 水位步长与栅格分辨率的最佳配合做多水位方案时水位步长和DEM分辨率之间有个经验配比。DEM水平分辨率为10米时水位步长选0.5米产生的面积增量往往在平原区突变明显而山区变化平缓若水位步长小于DEM垂直精度比如DEM垂直误差有0.3米那0.1米步长就是虚假精度算出来的面积差异全是噪声。我通常先查DEM属性表里的标准差取标准差的一半作为最小水位步长。研究区是陡峭山区DEM标准差可能有30米那水位步长取1米甚至2米就够了平原地区标准差只有2米水位步长取0.2米才能反映微地形的淹没差异。分辨率方面如果原始DEM是30米不建议直接重采样成10米来追求精度因为插值出来的假细节会误导淹没边界正确做法是保持30米分析但在出图前做一步Aggregate的平滑只改善视觉效果不改变计算值。5.3 三条常见的失败路径和对应的排除方法图表常见问题与检查顺序现象直接原因排查处理手段淹没范围完全空白坐标系不一致或源点落在NoData区用Project Raster统一到DEM坐标系IsNull检查源点像元无源淹没面积远超预期DEM没有填洼内部盆地全部积水做Fill后重跑Con并对比填洼前后DEM的差异三维体块有破洞水位面与地形面三角形尺寸差异过大调整z_tolerance重建TIN或把水位面重采样为DEM两倍像元大小三个问题都出在“数据级别不匹配”。Coordinated Universal Time的错位、NoData填充值不同、像元对齐不一致是最容易忽略的三个点。展开说Con判断时如果水位值落入NoData区结果是不参与计算的而在ArcScene里显示时NoData默认是透明看起来像被掏空。最实用的检查方式是打开栅格的属性表看COUNT和唯一值分布如果唯一值只有0和1且COUNT远小于整个范围的像元数说明有NoData混入用IsNull做掩膜把NoData填成0或极大值再做分析。5.4 一套可复用的批处理模板最后给出一套可以改参数直接跑的批处理骨架把DEM路径、水位列表、输出目录三者分离适合做多情景淹没分析。水位列表可以是多个断面水位也可以按时间增长的演进序列。这个模板我在多个防洪评价项目里用过替换数据路径后基本不用改逻辑。import arcpy from arcpy.sa import * arcpy.CheckOutExtension(Spatial) dem rD:\flood\dem_filled.tif out_dir rD:\flood\scenarios levels [32.5, 33.0, 33.5, 34.0, 34.5, 35.0] for h in levels: out_name finund_{str(h).replace(., _)} inund Con(dem h, 1, 0) inund.save(fr{out_dir}\{out_name}.tif) poly RasterToPolygon(inund, simplifyTrue) poly.save(fr{out_dir}\{out_name}.shp) print(f水位{h}完成范围面积{poly.getArea()})打印的面积实际需要从结果中读取这里展示的是循环写法多情景的关键是水位列表的生成方式可以用range(320, 361, 5)除以10生成等差序列也可以直接从CSV文件读取水文站的预报水位。批量输出的栅格可以用MosaicToNewRaster合并成多层栅格提供后续制图使用。三维模拟方面把每个水位对应的淹没范围批量拉伸并存入同一个场景文件逐层显示就可以制作出随时间上涨的动画效果比单张静态图清晰得多。本文还有配套的精品资源点击获取
