本刊网刊
城市内涝引发道路积水,不仅会降低道路的通行效率,还会对居民的生命财产安全构成威胁
[1‒3]。不同淹没深度和淹没范围的内涝,对应的内涝风险也不尽相同
[4]。如何高效准确地获取城市内涝积水状态及发展态势,对于城市防灾减灾和应急管理至关重要。
2013年,住房和城乡建设部发布了《城市排水防涝设施普查数据采集与管理技术导则(试行)》,各城市陆续开始布设内涝监测设备,基本实现了低洼积水点内涝的监测覆盖
[5‒6]。然而,水位传感器的安装、率定及维护成本较高,限制了其在城市空间上的规模化布设
[7]。以2020年深圳市为例,据统计全市仅有184个内涝专用水位传感器,离内涝积水精细化监测的需求尚存在一定差距。受热岛和雨岛效应的影响,城市暴雨常显示出明显的空间异质性
[8]。特别是短历时的强降雨,其空间异质性尤为显著
[9‒10]。内涝监测点位与内涝实际发生位置在空间上往往存在不匹配的状况,导致内涝积水缺测,尤其是特大暴雨期间。
事实上,在社交媒体动态、视频监控等城市大数据中存储着大量城市内涝积水信息
[11‒13]。内涝发生初期,通过社交媒体发布的内涝受灾信息与内涝事件几乎同步
[14‒15]。此外,城市道路交通摄像头能够完整地记录城市道路的积水过程,且密度远大于内涝积水传感器,成为城市内涝积水的潜在智能监测终端
[16‒17]。以北京市为例,2018年交通摄像头数量约115万个,平均覆盖度约71个/km
2,密度远大于内涝专用传感器
[18]。
随着视觉AI和深度学习的快速发展,利用视频和图像大数据开展内涝积水动态识别成为内涝研究的热点。当前研究主要围绕车辆、行人以及交通参照物来开展内涝积水的深度识别。以汽车为参照物的研究,常利用关键点检测算法识别汽车关键部位,如车轮、排气管、车把手等,并以此来估算积水等级
[11,19‒20]。以行人为参照物的研究,常将目标检测与人体关键点检测相结合,通过识别涉水图像中人脸、手部、身体、脚踝等关键点位来推求内涝积水等级
[21‒23]。以交通标志物为参照的研究,常引入具有一定标准尺寸的交通标志,如交通警示桶
[24]、停车让行标志
[25],来开展内涝积水等级识别。
现有研究多聚焦于利用互联网公开图像或社交媒体数据开展城市内涝积水识别,主要通过不同的深度学习模型或图像增强算法提升检测精度。然而,此类研究的数据来源分散、获取时序不连续,难以在固定空间位置上持续追踪积水变化过程,因而难以支撑内涝积水演变特征的动态反演与定量分析。本文提出了一种基于视觉AI的城市内涝积水特征动态识别方法,依托高清摄像头实时采集的视频监控数据,对同一监测点的内涝积水过程进行时序化识别与反演,重建积水深度随时间变化的动态演变过程,从而刻画内涝发生、发展与维持阶段的关键特征,显著提升了城市内涝监测的时空连续性与工程应用价值。
1 实验算法
1.1 算法流程
本文算法流程按“模型构建—训练策略—性能评估”展开。首先,基于YOLOv8的通用网络架构构建积水等级识别模型;其次,按既定参数迭代训练,采用主干网络冻结等策略提高收敛稳定性;最后,以准确率、平均精度等指标进行性能评估。
1.2 算法介绍
YOLO(you only look once)算法
[26]作为一阶段检测算法的代表,自2015年提出以来经过不断优化迭代,在检测速度和精度上均取得了持续的突破。YOLOv8算法以其实时性强、推理效率高的特点,被广泛应用于各类目标检测场景中。YOLOv8网络由骨干网络(backbone)、特征融合层(neck)、检测头(head)3部分组成,YOLOv8网络结构如
图1所示。
图1中:Input表示输入图像;Conv表示卷积模块,用于提取与变换特征;C2f表示特征提取/融合模块,用于增强特征表达并改善梯度流;SPPF表示快速空间金字塔池化模块,用于引入多尺度上下文信息;Upsample表示上采样操作;Concat表示特征拼接;检测头采用解耦结构,其中Conv2d表示二维卷积层;BboxLoss表示边界框回归损失;ClsLoss表示分类损失,用于指导模型参数优化。
1.3 算法训练
实验在Windows 10操作系统中开展YOLOv8算法训练,深度学习框架采用PyTorch 1.10版本,训练硬件采用NVIDIA RTX 3050 GPU显卡。
YOLOv8模型的主要训练参数设置如下:
1)批次大小(batch size)为16;
2)输入图像尺寸(input size)为640像素×640像素;
3)优化器(optimizer)为AdamW;
4)初始学习率(initial learning rate)为0.001,采用余弦退火策略进行动态学习率调整;
5)最大允许迭代轮数(epochs)为300;
6)backbone冻结策略:前10个epoch阶段冻结主干网络参数,仅训练neck和head部分,以加快早期收敛速度并稳定特征提取;
7)模型检查点(checkpoint)保存频率:每5个epochs保存一次,同时在验证集上评估模型性能。
1.4 算法评估
为有效评估YOLOv8算法在城市内涝积水特征动态识别任务上的表现,本文采用5项关键评估指标,即准确率(Precision,记为P)、召回率(Recall,记为R)、F1分数(F1 Score,记为F1)、精确率‒召回率(P‒R)曲线及平均精度均值(mAP,mean Average Precision,记为mAP)。
1)准确率,衡量模型对阳性类别预测的准确性,定义为真阳性预测与总预测阳性的比率。计算公式为:
式中,TP为真阳性数量,FP为假阳性数量。
2)召回率,评估模型检测数据集中所有相关样本的能力,定义为真阳性与真阳性和假阴性总和的比率。计算公式为:
式中,FN为假阴性数量。
3)F1分数,精确度和召回率的调和平均值,综合反应模型整体检测性能。计算公式为:
4)精确率‒召回率(P‒R)曲线,通过绘制不同置信度阈值下的P和R变化,反应模型在不同检测敏感度下的性能表现。曲线下的面积越大,说明模型整体检测效果越好。
5)平均精度,反应模型在所有类别上的综合检测精度。对单类别而言,平均精度(AP,记为AP)的定义为精确率‒召回率曲线下的面积:
式中,P(R)为准确率随召回率变化的关系。
多分类情况下,取各类别AP的平均值得到mAP:
式中,为第i个分类的平均精度,N为分类总数。
2 实验数据
2.1 样本数据集
本文构建了一个以汽车为参照物的内涝积水样本数据集,共计1 911张,来源包括两部分:互联网公开视频/图片平台收集的典型内涝场景图像1 780张,是以城市内涝、车辆淹水为关键词搜索得到的典型城市内涝事件相关公开图片;城市道路监控摄像头视频的抽帧图像131张,主要来自11场强降雨过程的固定监控视角样本,包含了强降雨期间16场内涝积水过程。两类样本均按“画面清晰、车辆涉水特征可辨识、可标注等级”的标准筛选,剔除重复与模糊帧。总体上,白天/弱光场景样本约占90%,夜间路灯/车灯照明条件下样本约占10%。
将尺寸统一为640像素×640像素、色彩归一化,使训练样本数据集的所有原始图像均完成标准化预处理,以适配模型的数据输入要求。采用YOLOv8默认在线增强设置(Mosaic马赛克拼接、随机仿射、HSV扰动等)。
2.2 样本数据标注
使用LabelImg工具对训练样本数据集中的车辆参照物进行手动标注,通过绘制矩形边界框框定车辆参照物及其积水级别,内涝积水等级标签制作如
图2所示。对1 911张内涝积水图片进行标注后,共生成了5 260个车辆样本实例。每个车辆样本实例包含了标注框左上角与右下角顶点像素坐标、中心点坐标、宽度、高度信息,以及对应的积水级别。
依据积水程度由轻到重,结合车辆行车安全,将积水等级划分为5级,每个级别对应不同的积水等级和车辆通行风险:
第0级(Level 0),无明显积水,车辆可正常行驶;
第1级(Level 1),积水低于轮胎一半,车辆可正常行驶,但存在轻微进水风险;
第2级(Level 2),积水超过轮胎一半,但未完全覆盖保险杠,存在较高进水风险;
第3级(Level 3),积水超过保险杠位置,但未淹没发动机前盖,车辆通行困难;
第4级(Level 4),积水超过发动机前盖,但未完全淹没车辆,车辆无法通行。
2.3 样本标注评估
图3为积水等级标注数据分布分析。
图3(a)中,内涝积水等级标注样本在5个不同积水等级均超过800个,其中:Level 2样本数量最多,占比约35%;其余等级分别在15%~20%之间,整体分布相对稳定。
图3(c)中,
X、
Y为车辆目标标注框中心点在图像中的归一化坐标。由
图3(b)、(c)可以看出,目标中心主要集中于图像中心区域,同时涵盖了从小尺度到大尺度的多类目标,这与图像拍摄的行为习惯一致。
图3(d)中,大多数目标框的尺寸集中在0.1~0.4之间,代表以小中型目标为主,同时也包含部分宽度和高度较大的目标,体现出一定的尺寸延展性,这种尺度分布符合城市内涝视频监控中车辆的特点。
2.4 视频监控数据
1)日间道路监控视频(2024‒04‒27,广东省某市):由道路监控高清摄像头拍摄,时间跨度为37 min(采用倍速存储,获取的视频文件实际时长为1 min55 s),完整地记录了强降雨条件下积水从初现、加深至形成的全过程。视频中一辆灰色面包车因积水过深熄火,成为整个积水过程中的固定参照物,为后续基于车辆的积水等级识别提供了稳定目标。
2)夜间停车场监控视频:由停车场监控摄像头拍摄,时间跨度为66 min(采用倍速存储,获取的视频文件实际时长为55 s),记录了长历时降雨条件下积水由浅到深的演变过程。该场景受路灯/车灯照明、频闪与水面反光影响,积水边缘特征更模糊。画面中多辆静止轿车构成稳定参照系。
3 实验结果
3.1 积水等级识别模型
模型训练设置最大迭代轮数为300,模型在第274个训练轮次到达收敛,自动终止。积水等级识别模型性能评估如
图4所示。
图4中,@0.5表示预测框与真实框的交集面积与并集面积之比,即预测框与真实框的IOU。由
图4(a)、(b)可见:模型在训练初期(第1个训练轮次)损失较高,分别为3.69(训练)和4.12(验证);随着迭代进行,损失持续下降,第274个训练轮次稳定在0.80(训练)和1.08(验证),表现出良好的收敛性和训练稳定性。
图4(c)表明,积水等级识别算法对Leve 4及Level 2识别的表现较好,曲线整体位于图表上方。不同积水等级在IOU大于0.5时检测性能存在显著差异,在Level 4严重积水检测精度最高达0.774,而Level 1轻微积水仅为0.577。
图4(d)显示,不同积水等级在置信度0.5附近的检测性能存在差异,其中,Level 4(严重积水)的F1分数最高,而对Level 1的识别效果较差。这主要是因为Level 1样本的视觉特征与Level 0相似,易造成模型混淆,从而导致了较低的检测性能,说明模型对轻微积水的识别能力相对较弱,但在中低置信度下仍具有较好的稳定性和泛化性。整体而言,模型对道路积水等级的识别具备较好的检测效果。
图5为积水等级识别模型在道路积水等级识别中的部分代表性结果。由
图5(a)可知,积水等级识别模型能够较为准确地识别图像中央横向区域内的车辆,而对边缘位置车辆检测中仍存在漏检现象,在处理小尺度或远距离目标时的识别能力有限。从
图5(b)可以进一步确认,积水等级识别模型能准确识别图像中心区域的道路积水等级,而对边缘区域的车辆目标识别能力有限。
积水等级识别模型总体mAP为0.673,整体性能处于中等水平。其性能差异与
图3的样本分布特征密切相关。一方面,
图3(a)的Level 2样本数量最多,为模型提供了充分的训练数据,因此训练表现最为稳定;Level 4样本虽数量不及Level 2,但因其视觉特征最为显著,模型仍能较好识别。而Level 1样本虽占一定比例,但其视觉特征与Level 0相似,易造成混淆,导致检测性能较低。另一方面,
图3(b)~(d)表明训练数据中车辆样本主要分布在图片中央位置,边缘区域样本数量相对不足,是造成边缘区域小目标的错检和漏检的一个主要原因。考虑到交通监控摄像头多采用中近景拍摄,且车辆主要分布于监控画面的中心区域,因此该模型在基于监控视频的实际应用中受边缘样本偏差影响较小,具有良好的适用性。
为评估本文提出的积水等级识别模型性能水平,进一步与近期基于YOLOv8的相关研究进行对比。Panja等
[27]在雨雾等恶劣天气条件下的车辆检测任务中报告YOLOv8的基准mAP约为0.685;本文在暴雨监控场景下叠加积水等级判别取得的mAP为0.673,整体精度与该基准水平接近。Qiu等
[28]针对公交车积水深度分级任务获得最佳mAP为0.662;相比之下,本文在多车型混合场景下mAP达到0.673,处于同类研究的可比水平。
3.2 实景应用验证
本文的积水过程复原采用“关键帧+线性插值”的技术路线:首先,在监控视频中选取具有明确物理意义的关键帧,基于车辆特征点与参照物几何关系对关键时刻的积水深度进行标定;随后,在相邻关键帧之间采用线性插值,得到连续的积水深度‒时间演变曲线。由于研究区域缺乏同位置水尺/传感器水位观测数据,以车辆参照的图像量测水深为“参考真值”,并在关键帧间隔内选取若干时刻进行人工量测复核;以插值曲线在对应时刻的取值为“模拟值”,两者差异用于表征过程重建误差与趋势非一致性。
本文选取两类具有代表性的真实监控场景开展验证:其一为2024‒04‒27日间道路监控场景(短历时强降雨触发的快速积水过程);其二为夜间停车场监控场景(长历时降雨背景下的缓慢累积过程)。
3.2.1 日间道路监控场景
为复原2024‒04‒27日间内涝积水时序演变过程,本研究选取了4个关键帧,日间场景内涝积水深度识别结果(积水深度‒时间曲线)如
图6所示,关键帧①~④分别对应第一次达到Level 1~Level 4积水等级的时刻,以及两个具有明显物理意义的人工参考点(
图7),涵盖从“首次漫过路缘石(
图7(a)关键帧⑤)”到“积水达到最深(
图7(b)关键帧⑥)”的典型阶段。参照《城市道路工程设计规范》(CJJ 37—2012),此次路缘石外露高度选取10.0 cm。面包车不同积水等级和典型特征点高度由实际测量获得,Level 2~Level 4以及积水最高点尾灯上沿对应深度分别对应28.5、58.0、98.0、113.0 cm。
图6复原了该场景下积水深度随时间的动态演变规律,展示了城市内涝积水过程从积水形成、积水扩张到积水达到巅峰的全过程。由
图6可知,该场积水深度变化表现出显著的阶段性,整体近似呈“S”形曲线。整个内涝积水过程可划分为3个典型阶段。1)积水初期扩张阶段(积水开始至第9 min20 s):随着降雨开始,地表径流迅速汇集并形成路面积水。此刻作为
图6积水过程的0点时刻(视频时间为16:21:00),如
图7(c)关键帧①所示。随后,地表径流开始横向扩张,约4 min后(视频时间为16:25:00)漫过10.0 cm的路缘石,迅速扩散至整个路面(
图7(a)关键帧⑤)。第9 min20 s(视频时间为16:30:20),积水首次达Level 2等级,对应深度为28.5 cm,如
图7(d)关键帧②所示。此阶段内,积水深度在9 min20 s内上涨了28.5 cm,平均上涨速度约3.05 cm/min。2)积水快速加深阶段(第9 min20 s—第17 min20 s):积水进入快速上涨期,积水深度在随后的2 min20 s内从Level 2上涨至Level 3,对应视频中16:32:40、积水深度为58.0 cm,如
图7(e)关键帧③所示。随后在第17 min20 s第一次达到积水等级Level 4,此时对应视频时间为16:38:20、积水深度为98.0 cm,如
图7(f)关键帧④所示。此阶段内,积水深度在8 min内从28.5 cm上涨至98.0 cm,平均上涨速度约8.69 cm/min。3)高水位维持阶段(第17 min20 s以后):随后,积水深度继续缓慢攀升,最终在第31 min43 s(视频时间为16:52:43)达到积水深度最高点,即面包车尾灯上沿高度113.0 cm处(
图7(b)关键帧⑥)。此阶段内,积水上涨速度锐减至1.05 cm/min。在随后几分钟内,水位几乎未见明显回落,形成持续性的高水位状态。这一特征反映出道路排水系统已完全饱和,雨后汇水无法及时排出,呈现“高位滞水—缓慢退水”的典型特征。
图6中,关键帧及人工参考点水深点线性插值得到的过程曲线(实线)与人工量测参考真值点(每隔3 min采点)整体趋势一致,能够刻画积水由缓慢增长到快速加深并进入高水位维持阶段的演变特征。以参考真值为基准,插值曲线在对应时刻的取值作为模拟值,计算误差:RMSE为3.81 cm、MAE 为3.30 cm、最大偏差为6.91 cm;以峰值水深113.0 cm计,3个指标相对误差分别约为3.37%、2.92%与6.12%。
3.2.2 夜间停车场监控场景
为复原2025‒07‒21夜间长历时降雨内涝积水过程,本文选取3个关键帧作为积水深度演变节点,分别对应第一次达到积水等级Level 1~Level 3的时刻,夜间场景内涝积水深度识别结果(积水深度‒时间曲线)如
图8所示。小轿车在积水等级Level 2、Level 3对应的积水深度分别为28.5、39.0 cm。夜间场景积水深度识别关键帧如
图9所示。
由
图8可见,该场积水过程亦呈现分段演变特征,可划分为两个阶段。1)积水缓慢累积阶段(积水开始至第50 min45 s):将首次可辨识积水(Level 1)时刻作为
图8积水过程的0点,视频时间为02:23:05,如
图9(a)关键帧①所示。在第50 min45 s(视频时间为03:13:50)积水首次达到Level 2等级,对应深度28.5 cm,如
图9(b)关键帧②所示。此阶段内,积水深度在50 min45 s内上涨了28.5 cm,平均上涨速度约为0.56 cm/min。2)积水加速增长阶段(第50 min45 s—第56 min47 s):积水深度在随后的6 min2 s内由Level 2上涨至Level 3,对应的视频时间为03:19:52。此阶段水深从28.5 cm快速攀升至39.0 cm,平均上涨速率提高至1.74 cm/min。
在误差验证方面,该场景下视频中每间隔约5 min选取特征时刻开展人工图像量测复核(
图8)。经统计计算得到RMSE为2.34 cm、MAE为1.99 cm、最大偏差为4.61 cm。以该场峰值水深39.0 cm来计算,3个指标相对误差分别为6.00%、5.10%与11.82%。上述结果表明,关键帧标定结合线性插值的方法在夜间低照度场景下仍能以较低误差复原积水深度的时序演变过程。
与日间道路监控场景相比,夜间低照度场景在积水等级判别与过程复原上更具挑战。其一,在积水初始形成阶段,浅水在弱光条件下对比度低,且易受路灯/车灯反射干扰,导致水体边界与地面纹理的区分度显著下降,Level 0与Level 1的临界状态即便人工目视亦难以稳定判定,如
图9(a)关键帧①与
图9(d)关键帧④所示。其二,夜间案例属于长历时缓变过程,单位时间水深增量较小,图像量测误差与水面微波动在相对尺度上更易被放大,从而使相对误差指标较日间快速上涨过程更为敏感(本文夜间RMSE/MAE约为6.00%/5.10%,均高于日间的3.37%/2.92%)。总体而言,夜间场景在浅水分级与缓变过程的敏感性方面更具挑战,但“关键帧标定+线性插值”的路线仍能在低照度条件下为无资料地区提供可用的积水过程。
4 结论与展望
本文针对内涝过程观测缺测问题,提出了一种基于道路监控视频与视觉AI的积水等级识别与过程复原方法。构建以车辆为参照物的积水数据集(1 911张),基于YOLOv8训练得到识别模型,总体mAP为0.673,可实现多等级积水自动识别。实景应用中采用“关键帧标定+线性插值”复原了两场降雨积水过程:日间暴雨场景积水曲线呈典型“S”形特征,积水初期扩张、积水快速加深和高水位维持3个阶段平均上涨速率分别为3.05、8.69和1.05 cm/min;夜间低照度长历时场景呈“缓慢累积—加速增长”特征,关键阶段速率分别为0.56和1.74 cm/min。重建的积水过程误差处于可控范围,表明该方法可依托既有监控网络为缺乏水尺/传感器的区域提供低成本、非接触式的积水过程数据支撑,为内涝风险研判、应急调度与城市雨洪模型参数率定提供数据支撑。
从鲁棒性与适用性看,夜间弱光与反光背景会显著削弱浅水阶段可辨识性,Level 0、Level 1的视觉特征尤其接近,即便人工目视亦难稳定判别;同时,长历时缓变过程单位时间水深增量较小,使量测扰动在相对尺度上更易放大,导致夜间的相对误差较日间快速上涨,场景更敏感。
后续研究将从“数据扩充—反演增强—机理耦合”3方面深化:一是补充夜间与雨雾等恶劣天气样本,形成日间道路暴雨与夜间低照度两类典型应用场景;二是引入路缘石、护栏、交通标线/标志等多参照物联合标定与插值,提升积水过程复原的精度与可靠性;三是尝试与雨洪模型耦合用于率定与精细化模拟,逐步实现监测—模拟—预警的闭环应用。
国家重点研发计划项目(2022YFC3202602)
中央高校基本科研业务费专项项目(226‒2025‒00245)
国家自然科学基金项目(52261160379)