前文ArcGIS分析城市在大暴雨中有洪涝隐患的区域(上)_Atom数据的博客-CSDN博客通过DEM数据识别了城市在大暴雨中可能遭受洪涝的区域——蓝点BlueSpots,但是,有蓝点并不代表就一定会被淹。一般来说,蓝点深的地方可以储蓄更多的降雨,相比蓝点浅的地方更不易被淹,位于可快速充满的蓝点内的建筑物要比位于填充速度较慢的建筑物面临的洪水风险更高。因此,为了更好地评估洪水风险,我们需要计算在暴雨中填充每个蓝点所需要的降雨量。
在水文概念里,景观中每个汇(蓝点)都有集水区或本地流域(其中的降水只流向某个特定蓝点的区域)。通过计算蓝点的体积及其流域的面积,可以求出将蓝点填充满所需的降雨量:即体积除以面积。这里我们考虑的是比较理想的状态,在大暴雨中可以忽略植被拦截、蒸发、土壤渗透等,同时也忽略流至下游蓝点的过量径流、蓝点内建筑物的基础高程等。
具体分析在ArcGIS分析城市在大暴雨中有洪涝隐患的区域(上)_Atom数据的博客-CSDN博客的基础上,增加了蓝点体积、流域面积等计算。
(1)(2)(3)(4)步不变,得到初始蓝点图层PreliminaryBS。
(5)使用”面转栅格“工具,像元大小2,将建筑物buildings转换为建筑物栅格BldgRaster,再使用”栅格计算器“,得到不包括建筑物的蓝点图层Bluespots Excluding Buildings。
(6)使用“区域分组”工具,将无建筑物蓝点图层按区域合并,并赋予唯一编号,得到BSRegions。
(7)计算蓝点容量。
首先使用”以表格显示分区统计“工具,统计类型选”总和“,计算未填洼深度,得到UnfilledDepth表格;再次使用”以表格显示分区统计“工具,计算已填洼深度,得到FilledDepth表格。
然后使用”连接字段“工具,将UnfilledDepth表格按照VALUE连接进BSRegions属性表,传输字段选”SUM“;再次使用”连接字段“工具,将FilledDepth表格以同样的方式连接进BSRegions属性表。
使用”添加字段“工具,给BSRegions添加”Volume“字段,字段类型选”浮点型(单精度)”;使用“计算字段”工具,用FilledDepth的SUM减去UnfilledDepth的SUM,再乘以栅格面积2*2,得到蓝点容量。
(8)计算流域面积。
先使用“流向”工具,计算出流向FlowDireciton;再使用“集水区”工具,得到分水岭Watersheds。
使用”构建栅格属性表“工具,为Watersheds创建属性表;使用”添加字段“工具,为Watersheds添加”WatershedArea“字段,字段类型选”浮点型(单精度)“;再使用”计算字段“工具,用”Count“字段乘以栅格面积2*2,得到流域面积WatershedArea。
(9)计算蓝点填充满所需降雨量。
先使用”连接字段“工具,把Watersheds的属性表按照”Value"连接到BSRegions的属性表里,传输字段"WatershedArea“;再使用”添加字段“工具,给BSRegions添加”FillUp“字段,字段类型选”浮点型(单精度)“;再使用”计算字段“工具,用Volume*1000/WatershedArea,得到FillUp值。
(10)使用“栅格转面“工具,将BSRegions转为面要素BSPoly;再使用”融合“工具,将BSPoly融合为BlueSpots,融合字段选择gridcode,得到Bluespot Polygons with IDs;再使用“连接字段”工具,把BSRegions的Volume和FillUp字段连接到Bluespot Polygons with IDs,得到具有容量和降雨量的Bluespots。至此,我们已经求出了蓝点及将蓝点填充满所需的降雨量。(注:在实际分析中,还需考虑下水道容量)
图中颜色越红的区域,越容易受洪涝影响。