专利名称：基于视频图像的河流流量在线测验实现方法
技术领域：
本发明涉及一种河流流量测验方法，尤其涉及一种基于视频图像的河流流量在线测验实现方法，属于非接触式明渠测流技术领域。
背景技术：
随着工农业生产和生活用水量的增加，水资源供需矛盾日益尖锐，导致河道断流现象及水污染灾害频发，已对沿河地区的经济、居民生活和生态环境等产生诸多不利影响。 同时，污水流量的长期监测也是一项艰巨的任务。由于排污渠道环境恶劣，其湿度高、杂物多，并可能具有腐蚀性的气体和液体，会给设备和操作人员带来危险。因此，要求流量测验仪器能够以最小的维护投入进行长期的自动测量。对于平原灌区，在整个渠系的排水过程中，若精确了解到某时某处的流量，可为整个渠系的计划用水和水量调配等工作提供量化参考，长期观测也可为水利规划设计和科学研究提供数据积累，因此如何迅速、准确地测得河流流量显得至关重要。目前用于河流流量在线测验的方法可分为以下两大类(1)建筑物测流法，采用堰槽量水设备(如巴歇尔槽)，让水流通过束缩的过水断面形成缓流到急流的过渡即临界流，保证稳定的水位-流量关系，然后采用非接触式水位计(如超声波水位计)测出水位，根据率定曲线计算出动态流量。该方法易于实现自动化， 但需对已有河道进行改造，投资大；并且，率定曲线是一种全局的基于模型的曲线，不能够捕获局部的水力影响，对于非均勻流及紊流通常不能给出精确的结果。因此只适用于少数小型的人工渠道。(2)流速-面积法，以测速垂线为界将河流断面划分为几个子区域，根据实测流速及水深计算各子区域的断面面积、平均流速及流量，然后将各个子区域的流量求和得到总的断面流量。根据传感器相对水流的位置关系，流速测量方法又可分为以下三类a)接触式。其中基于转子式流速仪的方法被认为是河流流量测验的标准方法。但由于流速仪测得的是点流速，历时长，需配合横跨河道的桥梁或缆道系统才能实现自动测验，因此采用铅鱼流速仪所测流速值的空间和时间分辨率都较差。此外，对于污水渠道，这种流速仪必须放置于污水内，因此不能可靠地用于长期监测。b)微侵入式。典型的是目前较为先进的声学多普勒流速剖面仪法(ADCP)。座底式或水平式ADCP可在短时内直接测得断面子区域中各水层的平均流速，实现河流流量在线测验，但存在盲区且其对含沙量较为敏感，通常只适用于形状规则的人工渠道；此外，由于换能器必须浸没于水中，易受高速水流及漂浮物的冲击并受河床冲淤影响，导致仪器运行的不可靠甚至损毁。c)非接触式。比较成熟的是电波流速仪法，它测量的是位于波束与水面交界处的一维表面流速，空间分辨率有限，并且当风或湍流在水面产生细小的波动时才能产生足够强的反射回波，因此不适用于水面平静的低流速河流。可见，现有的流量测验手段在原理上存在先天不足，受现场条件下流态和环境复杂度的影响以及技术和经济因素的限制，使目前河流流量的实时监测中存在着许多不精确和不确定性问题。尤其是在洪水、决提、大坝合龙等极端高流速环境下，以及针对流速极低而过水面积巨大的大尺度河流流量测验，目前还没有行之有效的测验手段。
发明内容
本发明针对现有河流测流技术存在的不足，而提出一种基于视频图像的河流流量在线测验实现方法。在该方法中设置有如下硬件系统a.用于实现水面流动可视化的水流示踪系统；b.用于拍摄和存储水面实时视频流的视频流采集系统；c.用于给图像正射变换提供地面控制点及其坐标信息的摄像机标定系统；d.用于处理数字图像的图像处理系统；e.用于计算河流流量的流量估计系统；该方法分为三个阶段(一 )测站部署阶段，包括如下步骤(1)选择测流河段；(2)获取测流断面水下地形；(3)部署现场硬件系统；(4)勘测地面控制点在世界坐标系下的平面坐标和高程坐标；(5)计算用于图像畸变校正的DLT变换系数；( 二)流场测量阶段，包括如下步骤(1)启动水流示踪系统投放示踪粒子；(2)启动视频流采集系统拍摄包含示踪粒子的水面视频流，然后按固定的时间间隔取帧并数字化处理；(3)由图像处理系统处理数字图像包括图像畸变校正、图像增强、流速矢量提取、流速场后处理，最后得到修正后的水面瞬时流速场；(三)流量估计阶段，该阶段是先勘测测流断面水位值，然后流量估计系统采用流速-面积法计算出河流流量。技术效果1.时空分辨率高。仪器测得的是水面瞬时二维流速矢量场，包含了水面流速分布的细节信息，而非传统仪器测量的点流速或线流速；本方法能在数秒内完成视频图像采集、 流速场获取和流量估计，并能在常规和极端水文过程期间对难以到达的测点实施密集测量。2.测量值可靠。当缺乏水面天然示踪物时，通过播撒人工示踪粒子的方式实现流动可视化，相比传统浮标法测流，粒子的跟随性更好，测量值更可靠。3.抗毁性强。仪器不接触水体，不会因为高速水流的冲击和漂浮物的缠绕而损毁， 维护成本低，可进行远程控制和长期监测。4.操作安全。仪器可布设于岸边、桥上或空中自治运行，工作人员可通过计算机网络远程操控，无需涉水测量，极大保障了测流人员的生命安全。
5.原理直观。相比声学和雷达测流技术，基于视频图像的方法更加直观明了 ；数字图像易于理解、处理、存储和传输，因此还可用于系统自身的故障诊断以及现场的工情监测。6.成本低廉，配置灵活。仪器采用了市面上成熟且通用的硬件产品，如CCTV摄像机、LED补光灯、图像采集卡和PC机等，系统易于获取和搭建，可根据具体应用灵活配置，并且随着新兴成像技术和图像处理技术的发展及应用，仪器的性能可得到进一步提升。7.可测量极低流速。非接触式的测流方式不会对水流造成干扰，适用于浅水、宽断面、低流速的河流。8.可测量极高流速。视频图像具有较高的分辨率和帧速率，覆盖了从几平方米到几百平方米的水面，可捕捉高速流动的水面目标。
图1是本发明方法涉及的系统架构示意图，图中标号名称1、水流示踪系统；1-1、 缆道；1-2、物料斗；1-3、示踪粒子；2、水面照明系统；3、视频流采集系统中的摄像机；4、水面；5、河岸；6、摄像机标定系统中的地面控制点；7表示时间平均水面流速；8表示水面插值流速；9表示深度平均流速；10、流量估计(测流)断面；11、勘测的水下地形。图2是本发明方法的流程图。图3是基于直接线性变换的图像畸变校正算法示意图。图4是基于时域中值滤波的背景差分图像增强算法示意图。图5是基于相关跟踪的流速矢量提取算法示意图。图6是基于空域中值滤波的流速场后处理算法示意图。
具体实施例方式下面对本发明方法作进一步说明。本发明方法中涉及的硬件系统架构如图1所示，主要是由水流示踪系统、水面照明系统、视频流采集系统、摄像机标定系统、图像处理系统和流量估计系统六个部分组成， 其中前四个部署在测验现场，后两个为监控中心计算机运行的处理和计算系统。下面对上述六个系统进行分别介绍1.水流示踪系统，该系统由示踪粒子和粒子投放装置组成。理想条件下，标识流动的粒子应当具备以下两个条件一是要能够精确跟随局部流体的运动，惯性及浮力是反应其是否符合流动特性的主要因素；二是要足够大使其能够被视频流采集系统的成像器件检测到。尽管在河流水面上可能存在一些天然示踪物，如轻质漂浮物、泡沫、冰块、悬浮质及细小的波纹，但它们的出现是随机且不稳定的，在低流速条件下的跟随性和可见性也参差不齐。本方法采用投放人工粒子的方式以在天然示踪物不足时进行补充，可获得比较理想的效果。具体选用一种廉价的用作包装填充材料的缓冲粒(又称泡泡粒)产品作为人工示踪粒子，它是由99%的玉米淀粉制成，是一种生物可降解的环保材料。粒子投放于水中数分钟后就会自行分解，不会对水体和生态环境造成污染；吸水后其密度与待测水流的密度相近， 能满足跟随性的要求；其部分浸没于水中可避免风力干扰引起的复杂度；并且白色的粒子能够与深色的水体形成对比，便于检测。粒子投放装置包括悬索式缆道和物料斗，物料斗悬挂于悬索上，其经缆道控制在水面上方做往复运动，内盛的示踪粒子从底部漏出播撒于水2.水面照明系统，该系统由光源及光强传感器组成。光源采用高亮卤素灯或LED 补光灯，功率和发射角度应当和照射距离及镜头焦距相匹配，这样才能在视场内获得较为均勻的光照；光强传感器用于感知环境光的强弱，以实现控制照明系统自动开启或关闭。3.视频流采集系统，该系统由测验现场的CCTV摄像机及监控中心的视频服务器组成。摄像机安装于云台上的防水盒内，可对其旋转角和俯仰角进行调节。摄像机拍摄PAL 制式(720X576像素，25帧/秒)或NTSC制式(640 X 480像素，30帧/秒)的模拟视频流， 通过同轴电缆上传至监控中心的视频服务器。视频服务器通过视频采集卡将模拟视频流进行数字化处理，一路经过压缩编码(MPEG-4或H. 264格式)转换成数字视频流存储于硬盘阵列上，可由本地/远程的客户端访问，或通过视频解码盒解码转换为视频信号在测控屏上显示；另一路不经压缩直接存储为原始数字视频帧数据(RAW或BMP格式)，由图像服务器调用。4.摄像机标定系统，该系统由至少六个地面控制点(GCP)组成，用于后期图像的几何校正。地面控制点是采用对测距激光具有强反射的铝合金材料制成的形状规则的十字形或正方形靶面。5.图像处理系统，该系统运行于监控中心的图像服务器上。它首先从视频服务器上读取原始的数字视频帧数据，然后对图像进行畸变校正和图像增强的预处理，接下来采用模式识别的方法获得水面示踪物在无畸变连续视频帧中的位置，再将测得的位移除以连续帧的时间间隔得到每个分析区域的流速矢量，对视场中包含的各分析区域重复上述过程得到原始的瞬时流速场，最后利用邻域信息对流速矢量进行校正得到修正的瞬时流速场。6.流量估计系统，该系统运行于监控中心的图像服务器上。它是利用水面流速场以及勘测的断面水下地形和水位值，基于流速-面积法估计出河流流量。图像服务器将得到的流速场及流量结果存储在本地硬盘上共客户端调用，并将信息叠加在水面的视频帧上通过VGA接口显示于测控屏。本发明方法的流程如图2所示，总体分为测站部署阶段、流场测量阶段和流量估计阶段三个阶段，下面对每一阶段分别说明。(一 )测站部署阶段，包括如下步骤内容1、选择测流河段。具体实施时按照以下条件选取测流河段(a)河槽的底坡、断面形状、糙率等因素比较稳定，槽段内的断面规则、水流平顺、流速分布对称均勻；(b)河段内无巨大块石、丛生水草阻水，无巨大漩涡、乱流等现象，当天然河道达不到以上要求时，必须进行人工整治使其符合测流条件；(c)宜选择河面较窄、通视条件好、横断面较单一、受风浪影响较小的河段，这样有利于示踪粒子的投放，并能减小图像远场分辨率不足对测流精度的影响。有条件的测点还可利用桥梁、堰闸布置测验。2、勘测断面水下地形。人工渠道的断面结构通常是已知的，而天然河道的水下地形可通过专用的声学仪器(如多波束声纳或声学多普勒流速剖面仪)直接测得。河流的水下地形勘测也可在测流时或测流后进行，前提是在该时间段内河床和边坡不会发生显著的变化。
3、部署现场仪器系统。水流示踪系统部署于紧邻测流断面的上游；水面照明系统部署于河流测流断面的两岸，光源可与视频流采集系统中的摄像机布置于同一位置，也可根据实际情况分散部署； 视频流采集系统中的CCTV摄像机架设于岸边制高点(可以是高层建筑、电线塔或临时搭建的高杆)上或横跨断面的桥梁上，调整云台的角度使其光轴尽量垂直于水流方向，且成像平面尽量平行于水面；由于在水面上放置控制点较为困难，在布设摄像机标定系统时，可将地面控制点平行放置于待测河段两岸河面最高水位以上的某一水平面处，并且尽量均勻而分散地覆盖完整的摄像机视场区域。4、勘测地面控制点坐标。采用全站仪、差分全球定位系统(DGPQ或激光测距仪测量各地面控制点(GCP)在世界坐标系Ow-XwYJw下的平面坐标(Ζ^, ^)以及高程坐标Zg，i = 1，L，N，这里N为GCP的数量。5、计算变换系数。拍摄一幅视频帧，提取出GCP的图像坐标，并基于直接线性变换的方法求解出用于图像畸变校正的DLT变换系数，供镜头畸变校正和图像正射校正使用。至此完成测站的前期准备工作。( 二)流场测量阶段，包括如下步骤内容1、控制照明并投放粒子。启动水面照明系统，光强传感器根据环境光强判断是否开启补光光源。启动粒子投放装置，往复运动的物料斗将内盛的示踪粒子均勻地播撒于水面，粒子群顺流而下进入待测区域。2、拍摄视频流并取帧。摄像机拍摄包含示踪粒子的水面视频流，由视频服务器取帧并数字化为8位灰度图像。3、处理数字图像。数字视频图像的处理流程如下(1)图像畸变校正。CCTV摄像机的光学系统一般为由多片透镜组成的光学镜头， 按理想化的小孔成像原理建模时存在着模型误差，引起包括径向畸变、切向畸变和薄棱镜畸变在内的二维图像的非线性畸变。此外，为了使图像包含较大的水流区域(通常是一个完整的测验断面)，摄像机通常处于一个倾斜于流动平面的视角下进行拍摄，由此引入了图像的透视畸变。本方法采用了一种基于直接线性变换的图像畸变校正算法，首先通过建立理想针孔摄像机模型、平面间变高单应关系以及镜头非线性畸变模型，建立镜头畸变模型下的图像坐标和2D世界坐标间的映射关系，接着采用最小二乘法求解变换系数，包括DLT 正变换系数、非线性畸变系数和DLT逆变换系数，最后进行镜头畸变校正和图像正射校正以获得真实世界中的度量信息，并用双线性内插法对整数坐标位置上的像素灰度进行插值，以便在图像上进行显示。(2)图像增强。通过观察可发现，水面背景像素的分布通常是均衡的，即单模态的，并且示踪粒子相对于整个水面场景是一种运动连续的小目标。但相比传统实验室中使用的粒子图像测速技术，现场环境下的测量条件更为复杂由于河流的开放性，不可避免地会在一定程度上受到外界环境的影响，如强风降雨在水面产生的扰动(动态)、阳光在水面产生的镜面反射(准静态)以及树木或建筑的倒影(静态)等。如果在图像预处理阶段不能合理地消除这些不利的水面模式，将会对后续的流速矢量提取环节造成影响，因此需要采用图像增强算法对水面灰度图像进行预处理。本方法结合水面图像的特点，设计了一种基于时域中值滤波的背景差分图像增强算法，背景差分可充分利用图像中变化部分的信息，使相关曲面的峰值锐化，而时域中值滤波背景估计的速度较快，保证所估计的背景像素可取到各帧输入图像中稳定时间最长的像素值。因此该算法能较好地抑制多种不利的水面模式，突出前景的运动粒子，从而达到图像增强的效果。(3)流速矢量提取。流速矢量提取算法是本方法的关键。在粒子密度中等、模式细节丰富且分布较为均勻的情况下，河流水面示踪粒子及其它水面模式的相对位置在相邻帧间隔内不会发生显著变化。由此，本方法依据模式识别的原理进行目标跟踪，具体采用了一种基于相关跟踪的流速提取算法，它将归一化的互相关系数作为模板图像与匹配图像灰度分布的相似性度量，通过对相关曲面中的峰值位置进行搜索及亚像素定位，确定第一帧图像中的分析区域相对于第二帧图像中对应窗口的最大似然位移。由于该算法利用图像的统计特征而不是粒子的几何特征进行匹配，单个示踪粒子的尺寸并不重要，因此对图像质量要求不高，可在低信噪比的条件下正常工作，能适应变化较复杂的流体运动目标。(4)流速场后处理。受示踪粒子浓度低和分布不均的影响以及外界环境的干扰，视场中的某些待测区域可能会暂时性地缺乏明显的示踪物，导致无法估计出这些区域(标记为NULL的分析区域)的瞬时表面流速场(矢量缺失)，或导致出现伪矢量(幅值或方向的偏差较大)。因此，需要采用一些准则以及全局或局部的滤波过程来减少错误矢量的个数， 本方法采用空域中值滤波的方法修正伪矢量，并对缺失矢量进行插补。至此获得修正后的水面瞬时流速场。(三)流量估计阶段。河流流量是单位时间内通过测量断面的水流体积，即流量=断面平均流速X过水断面面积。因此，本方法利用河流水下地形、水面流速场、水下纵向流速分布及水位值，基于流速-面积法估计河流流量。流速-面积法的原理是先以测速垂线作为划分界线，将河流断面划分为若干个子断面，再根据实测流速及水深计算各子断面的过水面积、平均流速及流量，然后将各个子断面的流量求和得到总的断面流量。从本质上看它依赖于河流的断面流速分布，在本方法中，水面流速场反映了水面横向流速分布的规律，而建立水面流速与深度平均流速间关系的流速系数则反映了纵向流速分布的规律。流量估计方法流程如下1、勘测河流测流断面水位。在确认河流水下地形参数有效(即在测流期间断面的水下地形未发生显著变化)的基础上，通过超声波水位计或水尺勘测实时水位值1。2、计算时间平均流速场。对于一段长度为N帧的视频序列共有N-I个图像对及其对应的瞬时流速场。对分析区域i的瞬时流速矢量Vi在时域进行信息融合，获得时间平均流速矢量Ui, i = 1，2，L，M，M为时间平均流速矢量总数，简单的融合方法可采用算术平均
Σ Y-υ = Af-i
1 ~ N-I然后，将所选勘测断面中分析区域i的时间平均流速矢量Ui投影到流动主方向X， 得到标量的时间平均水面流速U”
3、计算水面插值流速。将测流断面划分为以n+1条测速垂线为边界的η个子断面， 则子断面Sv的边界为垂线ν和ν+1。由于测速垂线ν不一定对应于分析区域i，因此对垂线相邻分析区域的时间平均水面流速Ui进行线性插值，求得测速垂线的水面插值流速U' v。4、计算深度平均流速。用水面插值流速U' v计算垂线的深度平均流速Z7V。假设每条测速垂线处纵向流速剖面的形状都相同，则深度平均流速Ζ7ν=α /，其中α为流速系数，它反映了纵向流速分布的规律，其取值依赖于剖面形状和河床的特性，在对数律的纵向流速分布和典型河床糙率的条件下，通常取α =0.85。5、计算子断面流量。根据经典的流速-面积法计算子断面Sv的流量Dv=Av(I)-Usv其中是子断面的平均流速，假定其大小为子断面两条边界垂线深度平均流速的均值，即是水位值为1时子断面Sv的过水面积，
Av (/) = ^(K+hv+l)wv ； hv和hv+1是子断面Sv两条边界垂线处的水深，由水下地形及水位值1
共同给出；wv是子断面Sv的宽度，即子断面Sv两条边界垂线之间的水平距离。6、计算断面流量。将各个子断面的流量求和得到河流的总断面流量至此完成河流流量估计。下面结合附图对本方法中涉及的部分算法展开说明。基于直接线性变换的图像畸变校正算法如图3所示。图像畸变校正实质上是建立 3D世界坐标系中的点在2D世界坐标系中的坐标与它在镜头畸变模型下的图像坐标之间的关系。校正原理建立在含畸变的针孔摄像机模型基础上，校正过程包括镜头畸变校正和图像正射校正，具体内容如下I、理想针孔摄像机模型的建立CCTV摄像机的成像系统建立在理想针孔摄像机模型的基础上，模型和图像畸变校正算法共包含5种坐标系反映真实物理空间的3D世界坐标系Ow-XwYwZw、反映摄像机方位的摄像机坐标系0。-xyz、反映(XD/CM0S感光器件所在平面的像平面坐标系O1-X1Y1、反映采集的视频帧在计算机上显示的图像坐标系Os-XsYs，以及正射校正后水面对应的2D世界坐标系 Ow-XwYw。其中，摄像机坐标系Oc-Xyz的坐标原点为透镜的光学中心0C，ζ轴与光轴重合；像平面坐标系O1-X1Y1的中心位于光轴与像平面的交点O1, X1, Y1平行于χ、y轴；图像坐标系 Os-XsYs的坐标原点&与2D世界坐标系Ow-XwYw的坐标原点Ow重合。对于3D世界坐标系中的点Pw0(w，Yff, Zw)，其在摄像机坐标系中的坐标为Pc(x，y， ζ)，在理想针孔模型下的像平面坐标为P1 (X1, Y1)，在镜头畸变模型下的像平面坐标为P1 (X1, Y1),在理想针孔模型下的图像坐标为Ps(xs，Ys)，在镜头畸变模型下的图像坐标为Ps(xs， ys)，在理想针孔模型下的2D世界坐标为Pw (Xw，Yw)，在镜头畸变模型下的2D世界坐标为 Pw (xw，yw) οII、平面间变高单应关系的建立
在理想针孔摄像机模型中，摄像机中心、图像点和目标点间存在共线的约束，图像和世界坐标系间的坐标空间变换可由共线方程推导得到如下线性关系
χ _ h-^w + h^w + h^w + h ν s j γ ι 7 γ ι 7 ν ι ι
ι9Λ W "‘" 1IO1IV """'11 乙『Tl
Y _ h-^w + h^w + Z7Zfj7 + Is
1S _
h-^w ho^w + hi^w +1
该变换即为直接线性变换(DLT)的正变换表达式，Ii称作DLT正变换系数，i = 1， 2，K，11。通常情况下，可以把水面近似看作一个平面，Zff = 0的直接线性变换描述了它和图像平面间的单应关系。由于在水面上放置控制点较为困难，在布设标定系统时可将GCP平行放置于河面最高水位以上高程坐标为Zg的某一平面处，Zg可表示为Z^ =CilX^+ dX +d3(2)其中，Clpd2分别为高程相对Zg、if方向变化的系数，d3(t)是t时刻GCP所在平面与水平面间的高度差，水位值由自记式水位计测得。假设河流表面的倾斜很小，则可将上式简化为Zg =式( )。将其代入直接线性变换公式中可得Xg = 炉(0 + 4 ⑶
h-^-w jrKJw +/^3(,) + 1YG = + ！6巧 + + K
h-^-w+/^3(,) + 1可见，结合水位变化进行修正的直接线性变换的方程组共有11个待求的未知量。III、镜头非线性畸变模型的建立为消除或减轻非线性畸变，在像平面坐标系中建立以下非线性畸变模型X1=X1+ R1X1 (Xj2 +yj) + Ii2X1 (Xj2 + yj)2 + px (3xf +yj) + 2p2xIyI + S1 (Xj2 + yj) (4)Y1 =yI+ Uj1 (Xj2 + 2) + ^y1 (Xj2 + yj f + IplXJ1 + p2 ^x21 +yj) + S2 (Xj2 + yj)式中1^1 为径向畸变系数，Pi、P2为切向畸变系数，S1^s2为薄棱镜畸变系数。根据实际情况可对该畸变模型进行简化，例如可只需要考虑径向畸变。因为引入过多的非线性参数有可能会引起解的不稳定性。此外，畸变像平面坐标(Xl，Y1)与畸变图像坐标(xs，ys)的关系可表示为X1=(Xs-Cx)S-Jdx (5)Y1 = (ys-Cy) dy式中CX、Cy为镜头光轴与像平面的交点在图像坐标系Xs、Ys轴上的坐标值，Sx为垂直比例因子，dx为C⑶在χ方向相邻像敏单元中心距，dy为C⑶在y方向相邻像敏单元中心距。同理，理想像平面坐标(X1, Y1)与理想图像坐标(Xs，Ys)的关系可表示为X1=(Xs-Cx)S-Jdx(Q)Y1 = (Ys-Cy) dy
11
对应的逆关系式为Xs =X1SxCi-J+Cx (7)Ys=Y1Ci-^Cy设a。=Sxd-J，则最终待求的系数为(Cx，Cy, a0, k2，Pl, p2, S1, s2)。IV、变换系数求解1.求解DLT正变换系数在测站部署阶段勘测N ^ 6个共面的GCP，将N个平面坐标值(Zf, 及其高程坐
标值Zg =式代入式(1)直接线性变换方程组的矩阵形式WL=I(8)式中
权利要求
1.一种基于视频图像的河流流量在线测验实现方法，其特征在于该方法中设置有如下硬件系统a.用于实现水面流动可视化的水流示踪系统；b.用于拍摄和存储水面实时视频流的视频流采集系统；c.用于给图像正射变换提供地面控制点及其坐标信息的摄像机标定系统；d.用于处理数字图像的图像处理系统；e.用于计算河流流量的流量估计系统；该方法分为三个阶段(一)测站部署阶段，包括如下步骤(1)选择测流河段；(2)获取测流断面水下地形；(3)部署现场硬件系统；(4)勘测地面控制点在世界坐标系下的平面坐标和高程坐标；(5)计算用于图像畸变校正的DLT变换系数；(二)流场测量阶段，包括如下步骤(1)启动水流示踪系统投放示踪粒子；(2)启动视频流采集系统拍摄包含示踪粒子的水面视频流，然后按固定的时间间隔取帧并数字化处理；(3)由图像处理系统处理数字图像包括图像畸变校正、图像增强、流速矢量提取、流速场后处理，最后得到修正后的水面瞬时流速场；(三)流量估计阶段，该阶段是先勘测测流断面水位值，然后流量估计系统采用流速-面积法计算出河流流量。
2.根据权利要求1所述的基于视频图像的河流流量在线测验实现方法，其特征在于 该方法中还设置有由光源和光强传感器组成的水面照明系统，所述水面照明系统部署于测流断面的两岸。
3.根据权利要求1所述的基于视频图像的河流流量在线测验实现方法，其特征在于 所述阶段(一)的步骤(3)的内容如下水流示踪系统部署于紧邻测流断面的上游，架设视频流采集系统中的CCTV摄像机，使其光轴垂直于水流方向，将摄像机标定系统中的地面控制点平行设置于测流河段两岸最高水位以上一水平面处，地面控制点覆盖完整的摄像机视场区域。
4.根据权利要求1或3所述的基于视频图像的河流流量在线测验实现方法，其特征在于所述摄像机标定系统是由至少六个地面控制点组成。
5.根据权利要求1所述的基于视频图像的河流流量在线测验实现方法，其特征在于 所述阶段(一)的步骤(5)的方法为利用视频流采集系统拍摄一幅视频帧，提取出地面控制点的图像坐标，基于直接线性变换的方法求解出用于图像畸变校正的DLT变换系数。
6.根据权利要求1所述的基于视频图像的河流流量在线测验实现方法，其特征在于 所述采用流速-面积法计算出河流流量的步骤如下(1)计算时间平均流速场对第i个分析区域的瞬时流速矢量Vi在时域进行信息融合，得到时间平均流速矢量Ui, 将Ui投影到流动主方向X，得到标量的时间平均水面流速Ui, i = 1,2, ... , M，M为时间平均流速矢量总数；(2)计算水面插值流速根据流速-面积法，将测流断面划分为以n+1条测速垂线为边界的η个子断面，对测速垂线相邻分析区域的时间平均水面流速仏进行线性插值，得到测速垂线的水面插值流速 U' ν；—(3)计算深度平均流速Z7v=af/，其中α为流速系数；(4)计算子断面流量计算子断面Sv的流量Dv = <(/)_ZL，其中&是子断面的平均流速，AJl)是水位值为1 时子断面Sv的过水面积；(5)将各个子断面的流量求和，得到所测河流的总断面流量，即求得河流流量。
7.根据权利要求1、3、4或5所述的基于视频图像的河流流量在线测验实现方法，其特征在于所述地面控制点是采用铝合金材料制成的十字形或正方形靶面。
全文摘要
本发明公开了一种基于视频图像的河流流量在线测验实现方法，属于非接触式明渠测流技术领域。该方法中涉及水流示踪系统、水面照明系统、视频流采集系统、摄像机标定系统、图像处理系统和流量估计系统六个部分，水流示踪系统在测流河段投放示踪粒子，摄像机拍摄水面视频流，然后按固定时间间隔取帧并数字化处理，图像处理系统根据摄像机标定系统提供的地面控制点坐标信息对数字图像进行处理，最后流量估计系统采用流速-面积法计算出河流流量。本发明在降低河流测流复杂度的同时能保证较高的测量精度，可实现极端环境下河流流量的动态监测。
文档编号G01F1/704GK102564508SQ20111041946
公开日2012年7月11日 申请日期2011年12月14日 优先权日2011年12月14日
发明者严锡君, 张振, 徐立中, 王鑫 申请人:河海大学