CN103578118A - 基于序贯图像矢量平均的时均流场重建方法 - Google Patents
基于序贯图像矢量平均的时均流场重建方法 Download PDFInfo
- Publication number
- CN103578118A CN103578118A CN201310508543.XA CN201310508543A CN103578118A CN 103578118 A CN103578118 A CN 103578118A CN 201310508543 A CN201310508543 A CN 201310508543A CN 103578118 A CN103578118 A CN 103578118A
- Authority
- CN
- China
- Prior art keywords
- vector
- flow field
- overbar
- time
- instantaneous
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种基于序贯图像矢量平均的时均流场重建方法,属于非接触式明渠测流技术领域。所述方法首先载入两帧连续图像作为一组图像对,并采用预先设定的掩膜提取图像对中的有效水面区域;其次采用快速傅立叶变换互相关算法对图像对中的所有分析区域IA进行运动矢量估计;然后采用全局角度直方图方法检测流动主方向并识别瞬时流场中的错误矢量;接下来进行流场时间滤波,舍弃当前错误的瞬时矢量,并用当前正确的瞬时矢量替代此前错误的时均矢量,用当前正确的瞬时矢量平均此前正确的时均矢量;最后进行错误矢量插值和流场平滑滤波得到重建的时均流场。本发明以时均流场的矢量正确率为依据控制求平均的进程,改善了流场重建的时间分辨率。
Description
技术领域
本发明涉及一种二维时均流场重建方法,尤其涉及一种基于序贯图像矢量平均的时均流场重建方法,属于非接触式明渠测流技术领域。
背景技术
天然河流中的水体在河槽中运动时受到断面形状、坡度、糙率、水深、弯道、风、气压、潮汐等因素的影响而产生紊流。紊流内部水质点的瞬时流速在大小和方向上均随时间变化,呈现出脉动现象,但在足够长的时段内其均值保持稳定。因此,对于面向断面流量监测的河流水面成像测速系统而言,时均流场的重建比瞬时流场的获取更为重要。然而在天然河流中,由于水面示踪物密度低、时空分布不均及水面光学噪声的影响,视场中的某些待测区域可能会暂时性地缺乏明显的示踪物,导致无法估计出这些区域的瞬时矢量或出现错误矢量。
对于基于灰度相关匹配法的运动矢量估计,按照求平均操作在处理流程中所处的阶段及作用对象的不同,可以分为图像平均、相关平均和矢量平均三种时均流场重建策略。图像平均策略的优点是仅需计算一次相关测度,运算速度快;缺点是需要一次读入完整的图像序列,占用大量存储空间;并且受到平均帧数的限制,对图像质量有较高的要求。相关平均策略的优点是能够有效提高相关曲面的信噪比,降低出现错误矢量的概率;缺点是对相关曲面的直流偏置噪声敏感,依赖于有效的图像背景抑制方法。矢量平均策略的优点是灵活、高效,可以处理任意方法获得的网格化矢量场;由于是时均流场重建的最高层次,可以获得瞬时流场以及特征级的中间参数,具有较高的容错性;缺点是当示踪物密度较低时,对相关曲面信噪比低下引起的错误矢量敏感。以上策略都在一定程度上提高了位移估计的可靠性,但当错误矢量个数超过一定比例(如30%)时,直接采用线性求平均的方法都将导致时均流场重建精度过低。因此错误矢量的识别和修正对低信噪比条件下的时均流场重建具有重要意义。
发明内容
本发明针对现有方法存在的不足,提供了一种基于序贯图像矢量平均的时均流场重建方法。
该方法通过以下技术手段实现:
方法首先载入两帧连续图像作为一组图像对,并采用预先设定的掩膜提取图像对中的有效水面区域;其次采用快速傅立叶变换互相关算法对图像对中的所有分析区域IA进行运动矢量估计;然后采用全局角度直方图检测流动主方向并识别瞬时流场中的错误矢量;接下来进行流场时间滤波,舍弃当前错误的瞬时矢量,并用当前正确的瞬时矢量替代此前错误的时均矢量,用当前正确的瞬时矢量平均此前正确的时均矢量;最后进行错误矢量插值和流场平滑滤波得到重建的时均流场。
所述的错误矢量识别方法,其实现过程如下:对于i时刻的瞬时流场,首先计算所有矢量Si(x,y)与x方向的夹角θi(x,y),得到瞬时流场中矢量角度的极坐标分布:
θi(x,y)=arctan(Δyi/Δxi) (1)
其中,(x,y)表示矢量的图像坐标,Δx、Δy分别表示矢量在x和y方向的大小;然后以Δθ为角度区间建立矢量方向在0-359°范围内的全局角度直方图;接下来依次搜索各角度区间,将矢量数量最多的区间标记为流动主方向θM;最后以角度阈值θTH为判别窗口识别每个矢量的类型,并建立矢量类型标志如下:
其中,1表示正确矢量,0表示错误矢量。
所述的流场时间滤波过程用下式描述为:
其中,分别为当前和此前时刻的时均矢量;Si(x,y)为当前时刻的瞬时矢量;错误矢量标记为“NULL”,以区别于零矢量;Ni(x,y)为参与当前矢量平均的正确矢量个数;Fi(x,y)分别为此前时均矢量的类型标志和当前瞬时矢量的类型标志,对二者进行或运算得到当前时均矢量的类型标志,如下:
其中,“|”表示或运算符。遍历每个所有IA后统计得到正确和错误矢量的个数N1和N0,则时均流场的矢量正确率为:
所述的错误矢量插值采用错误矢量邻域内正确矢量的均值对错误矢量进行插值:
所述的流场平滑滤波,采用时均流场和滤波器的空域卷积实现:
其中,“*”表示卷积运算,滤波器采用大小为3×3的归一化二维离散高斯窗,表示为:
本发明具有以下有益效果:
方法利用了断面流速方向一致性分布的特点,采用基于非线性统计的全局角度直方图方法检测流动主方向,较好地解决了示踪物密度低、时空分布不均及水面光学噪声影响下区域性错误矢量的识别问题,改善了河流水面流场的估计精度。此外,方法充分利用了运动目标时空分布的冗余信息,以时均流场的矢量正确率为依据控制求平均的进程,有效提高了流场测量的时间分辨率。相比图像平均和相关平均策略,本方法在实时性和硬件开销上具有明显优势,适用于河流水面流场、流量连续在线监测系统。
附图说明
图1是本发明涉及的时均流场重建方法的工作流程图。
具体实施方式
下面结合附图对本发明的技术方案进行详细说明:
本发明的时均流场重建方法的工作流程如图1所示。共分为以下七个步骤:
(1)载入图像对,从相机采集的图像序列中载入两帧连续图像作为一组图像对。
(2)水面区域提取,采用预先设定的掩膜提取图像对中的有效水面区域。
(3)运动矢量估计,采用快速傅立叶变换互相关算法对图像对中的所有分析区域IA进行运动矢量估计,得到当前图像对的瞬时流场;
(4)错误矢量识别,对于i时刻的瞬时流场,首先计算所有矢量Si(x,y)与x方向的夹角θi(x,y),得到瞬时流场中矢量角度的极坐标分布:
θi(x,y)=arctan(Δyi/Δxi) (1)
其中,(x,y)表示矢量的图像坐标,Δx、Δy分别表示矢量在x和y方向的大小;然后以Δθ为角度区间建立矢量方向在0-359°范围内的全局角度直方图;接下来依次搜索各角度区间,将矢量数量最多的区间标记为流动主方向θM;最后以角度阈值θTH为判别窗口识别每个矢量的类型,并建立矢量类型标志如下:
其中,1表示正确矢量,0表示错误矢量。
(5)流场时间滤波,基本思想是:舍弃当前错误的瞬时矢量,并用当前正确的瞬时矢量替代此前错误的时均矢量,用当前正确的瞬时矢量平均此前正确的时均矢量,描述为:
其中,分别为当前和此前时刻的时均矢量;Si(x,y)为当前时刻的瞬时矢量;错误矢量标记为“NULL”,以区别于零矢量;Ni(x,y)为参与当前矢量平均的正确矢量个数;Fi(x,y)分别为此前时均矢量的类型标志和当前瞬时矢量的类型标志,对二者进行或运算得到当前时均矢量的类型标志,如下:
其中,“|”表示或运算符。遍历每个所有IA后统计得到正确和错误矢量的个数N1和N0,则时均流场的矢量正确率为:
(6)错误矢量插值,采用错误矢量邻域内正确矢量的均值对错误矢量进行插值:
(7)流场平滑滤波,采用时均流场和滤波器的空域卷积实现:
其中,“*”表示卷积运算,滤波器采用大小为3×3的归一化二维离散高斯窗,表示为:
至此完成一次河流水面二维时均流场重建过程。
Claims (1)
1.一种基于序贯图像矢量平均的时均流场重建方法,其特征在于:首先载入两帧连续图像作为一组图像对,并采用预先设定的掩膜提取图像对中的有效水面区域;其次采用快速傅立叶变换互相关算法对图像对中的所有分析区域IA进行运动矢量估计;然后采用全局角度直方图方法检测流动主方向并识别瞬时流场中的错误矢量;接下来进行流场时间滤波,舍弃当前错误的瞬时矢量,并用当前正确的瞬时矢量替代此前错误的时均矢量,用当前正确的瞬时矢量平均此前正确的时均矢量;最后进行错误矢量插值和流场平滑滤波得到重建的时均流场;
所述的全局角度直方图方法,对于i时刻的瞬时流场,首先计算所有矢量Si(x,y)与x方向的夹角θi(x,y),得到瞬时流场中矢量角度的极坐标分布:
θi(x,y)=arctan(Δyi/Δxi) (1)
其中,(x,y)表示矢量的图像坐标,Δx、Δy分别表示矢量在x和y方向的大小;然后以Δθ为角度区间建立矢量方向在0-359°范围内的全局角度直方图;接下来依次搜索各角度区间,将矢量数量最多的区间标记为流动主方向θM;最后以角度阈值θTH为判别窗口识别每个矢量的类型,并建立矢量类型标志如下:
其中,1表示正确矢量,0表示错误矢量;
所述的流场时间滤波过程用下式描述为:
其中,分别为当前和此前时刻的时均矢量;Si(x,y)为当前时刻的瞬时矢量;错误矢量标记为“NULL”,以区别于零矢量;Ni(x,y)为参与当前矢量平均的正确矢量个数;Fi(x,y)分别为此前时均矢量的类型标志和当前瞬时矢量的类型标志,对二者进行或运算得到当前时均矢量的类型标志,如下:
其中,“|”表示或运算符。遍历每个所有IA后统计得到正确和错误矢量的个数N1和N0,则时均流场的矢量正确率为:
若小于一个阈值PTH(如90%),则继续处理下一时刻的图像对,直到大于PTH并且求平均的帧数达到预设的阈值,得到一个包含错误矢量的时均流场和对应的矢量类型查找表;反之则认为时均流场重建失效;
所述的错误矢量插值,采用错误矢量邻域内正确矢量的均值对错误矢量进行插值,表示如下:
所述的流场平滑滤波,采用时均流场和滤波器的空域卷积实现,表示为:
其中,“*”表示卷积运算,滤波器采用大小为3×3的归一化二维离散高斯窗,表示为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310508543.XA CN103578118B (zh) | 2013-10-24 | 2013-10-24 | 基于序贯图像矢量平均的时均流场重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310508543.XA CN103578118B (zh) | 2013-10-24 | 2013-10-24 | 基于序贯图像矢量平均的时均流场重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103578118A true CN103578118A (zh) | 2014-02-12 |
CN103578118B CN103578118B (zh) | 2016-02-24 |
Family
ID=50049838
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310508543.XA Active CN103578118B (zh) | 2013-10-24 | 2013-10-24 | 基于序贯图像矢量平均的时均流场重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103578118B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103793924A (zh) * | 2014-01-28 | 2014-05-14 | 河海大学 | 一种基于fht-cc的流场图像自适应运动矢量估计方法 |
CN103996171A (zh) * | 2014-05-05 | 2014-08-20 | 河海大学 | 基于时空图像的流体运动矢量估计方法 |
CN104156920A (zh) * | 2014-08-07 | 2014-11-19 | 江苏大学 | 一种用于离心泵叶轮流场图像处理的方法 |
CN111444614A (zh) * | 2020-03-26 | 2020-07-24 | 西安交通大学 | 一种基于图卷积的流场重构方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6665450B1 (en) * | 2000-09-08 | 2003-12-16 | Avid Technology, Inc. | Interpolation of a sequence of images using motion analysis |
CN102866260A (zh) * | 2012-09-18 | 2013-01-09 | 河海大学 | 非接触式河流表面流场成像量测方法 |
CN103336144A (zh) * | 2013-05-27 | 2013-10-02 | 中国船舶重工集团公司第七○二研究所 | 一种水面微速流场的测试方法 |
-
2013
- 2013-10-24 CN CN201310508543.XA patent/CN103578118B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6665450B1 (en) * | 2000-09-08 | 2003-12-16 | Avid Technology, Inc. | Interpolation of a sequence of images using motion analysis |
CN102866260A (zh) * | 2012-09-18 | 2013-01-09 | 河海大学 | 非接触式河流表面流场成像量测方法 |
CN103336144A (zh) * | 2013-05-27 | 2013-10-02 | 中国船舶重工集团公司第七○二研究所 | 一种水面微速流场的测试方法 |
Non-Patent Citations (3)
Title |
---|
J.LE COZ等: "Performance of image-based velocimetry ( LSPIV ) applied to flash-flood discharge measurements in Mediterranean rivers", 《JOURNAL OF HYDROLOGY》 * |
M.JODEAU等: "Application and evaluation of LS-PIV technique for the monitoring of river surface velocities in high flow conditions", 《FLOW MEASUREMENT AND INSTRUMENTATION》 * |
徐立中等: "非接触式明渠水流监测技术的发展现状", 《水利信息化》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103793924A (zh) * | 2014-01-28 | 2014-05-14 | 河海大学 | 一种基于fht-cc的流场图像自适应运动矢量估计方法 |
CN103996171A (zh) * | 2014-05-05 | 2014-08-20 | 河海大学 | 基于时空图像的流体运动矢量估计方法 |
CN103996171B (zh) * | 2014-05-05 | 2017-02-15 | 河海大学 | 基于时空图像的流体运动矢量估计方法 |
CN104156920A (zh) * | 2014-08-07 | 2014-11-19 | 江苏大学 | 一种用于离心泵叶轮流场图像处理的方法 |
CN104156920B (zh) * | 2014-08-07 | 2017-10-20 | 江苏大学 | 一种用于离心泵叶轮流场图像处理的方法 |
CN111444614A (zh) * | 2020-03-26 | 2020-07-24 | 西安交通大学 | 一种基于图卷积的流场重构方法 |
CN111444614B (zh) * | 2020-03-26 | 2021-12-28 | 西安交通大学 | 一种基于图卷积的流场重构方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103578118B (zh) | 2016-02-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7120582B2 (ja) | オプティカルフロー映像処理を利用する河川流速測定装置および方法 | |
CN102129695B (zh) | 基于遮挡物建模的有遮挡情况下的目标跟踪方法 | |
CN107705318B (zh) | 一种基于边界示踪的湍流边界层流场测速方法 | |
CN103578118A (zh) | 基于序贯图像矢量平均的时均流场重建方法 | |
CN101551901B (zh) | 动态遮挡图像的实时补偿和增强方法 | |
CN103295232B (zh) | 基于直线和区域的sar图像配准方法 | |
CN104732536A (zh) | 一种基于改进形态学的亚像素边缘检测方法 | |
Lin et al. | Applications of image recognition for real-time water level and surface velocity | |
CN104268883A (zh) | 一种基于边缘检测的时频谱曲线提取方法 | |
CN103593849A (zh) | 一种用于影像序列椭圆形人工目标点快速识别与跟踪方法 | |
CN102609945A (zh) | 可见光和热红外图像序列自动配准方法 | |
CN109766838A (zh) | 一种基于卷积神经网络的步态周期检测方法 | |
CN102261916B (zh) | 一种沙质环境中基于视觉的月球车定位方法 | |
CN104680553A (zh) | 基于双边滤波的可变核函数目标跟踪方法 | |
CN104574417A (zh) | 一种图像边缘灰度起伏性度量与自适应检测方法 | |
Zhen et al. | Design and evaluation of an FFT-based space-time image velocimetry (STIV) for time-averaged velocity measurement | |
CN104732530A (zh) | 一种图像边缘的检测方法 | |
CN104392466A (zh) | 一种基于视觉显著原始目标的检测跟踪方法 | |
CN101794395B (zh) | 一种基于Sift算法的图像匹配定位方法 | |
CN104166977B (zh) | 一种图像匹配相似性测度方法及其图像匹配方法 | |
CN103411562B (zh) | 一种基于动态规划和均值漂移的结构光光条中心提取方法 | |
CN117495918A (zh) | 基于光照自适应orb算子的河流水面光流估计方法 | |
CN102881025A (zh) | 一种多运动目标的检测方法 | |
CN104199006A (zh) | 一种基于多假设联合分布滤波器的随机集跟踪方法 | |
Modi et al. | Skew correction for vehicle license plates using principal component of Harris Corner Feature |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |