CN103996171B - 基于时空图像的流体运动矢量估计方法 - Google Patents
基于时空图像的流体运动矢量估计方法 Download PDFInfo
- Publication number
- CN103996171B CN103996171B CN201410185108.2A CN201410185108A CN103996171B CN 103996171 B CN103996171 B CN 103996171B CN 201410185108 A CN201410185108 A CN 201410185108A CN 103996171 B CN103996171 B CN 103996171B
- Authority
- CN
- China
- Prior art keywords
- time
- image
- space
- line
- speed
- 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.)
- Active
Links
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种基于时空图像的流体运动矢量估计方法,属于流场测量技术领域。所述方法首先采用单台相机以固定时间间隔采集流体图像序列;然后采用像差修正的直接线性变换法对图像进行光学畸变校正,并沿测流方向设置一组测速线,每条测速线均以像素为横坐标、以时间为纵坐标合成一幅时空图像;接下来对边缘检测后的时空图像进行频谱分析,并采用离散逼近法求解频谱主方向得到与其垂直的纹理主方向;最后根据物方坐标和像方坐标间的映射关系求解测速线的物理尺度,并结合纹理主方向和时间间隔计算一维流速矢量的大小,实现时均流场定标。本发明能够用于河流水面等天然示踪水流的时均流场测量,特别适用于高洪期河流流量的在线监测。
Description
技术领域
本发明涉及一种流体运动矢量估计方法,尤其涉及一种基于时空图像的流体运动矢量估计方法,属于流场测量技术领域。
背景技术
流体中存在着大量不具有固定质量并存在相对运动的流体微团,使得流体的运动形态比刚体和固体更为复杂。在过去的二十多年间,光学、电子及计算机技术的发展促进了粒子图像测速(PIV)这种非接触式瞬时全场流速测量技术的实现和应用。它通过粒子图像序列的分析和计算获得局部流体运动位移及速度的大小、方向、特征和分布情况,极大提高了实验室环境下各种复杂流动的测量能力。其中从粒子图像估计流体的运动矢量是PIV技术的核心和难点,方法的选择不仅依赖于硬件系统而且取决于待测流体的特点。
对于利用树叶、树枝等天然漂浮物以及漩滚、表面波等天然水面模式作为水流示踪物的大尺度粒子图像测速(LSPIV)应用而言,河流水面流场的运动矢量估计又面临着特有的挑战。首先,大多数天然水面模式产生于明渠紊流,它们没有明确的粒子特性,通常存在振荡、扩散及形变等非刚性运动,因而缺乏稳定的几何特征和统计特征,使得经典的粒子跟踪测速(PTV)方法难以直接适用。其次,现场条件下的待测流场面积往往覆盖数百至数干平方米,为满足流速测量的目标精度,图像分辨率需要达到百万至千万像素,远大于目前实验室中十万像素的数量级,使得基于直接互相关(DCC)等空域相关匹配的运动矢量估计方法计算量相当可观,难以满足流场、流量连续监测的实时性要求。再则,河流水面的流速分布在断面方向表现出从中泓区域向两岸递减的趋势,流速梯度加之图像透视畸变和水流示踪物分布不均的影响往往导致水面瞬时流场具有复杂的区域特征。采用固定大小的观测窗口会因为对局部流场的适应性较差而降低运动矢量估计的精度和可靠性。因此研究针对河流水面这类特殊流体的运动矢量估计方法具有重要的理论意义和应用价值。
发明内容
本发明针对现有方法存在的不足,提供了一种基于时空图像的流体运动矢量估计方法。
该方法通过以下技术手段实现:
由于满足质量守恒定律的流体,其运动在短时内通常满足连续性的假设。对于跟随性良好的水流示踪物,其在三维时空域中的位置必然满足某种相关性。这种相关性在一维图像空间和一维序列时间组成的时空图像中表现为较为显著的纹理主方向,如果能够精确估计出时空图像的纹理主方向就能获得流体在指定空间方向的运动矢量。然而从空间域中可以看出,对于天然示踪的流体,水流示踪物的时空分布不均往往造成时空图像受到非规则纹理或噪声的干扰,影响纹理主方向的估计精度。图像的边缘极大程度上反应了图像纹理的方向,并且具有方向性规则纹理的图像在傅立叶频谱中的能量主要聚集在过原点的直线上,直线指向代表的频谱主方向与纹理主方向垂直。因此利用时空图像边缘信息的频谱整体估计纹理主方向可以有效抑制背景噪声,提高估计精度和运算效率。
基于上述原理,方法首先采用单台相机以固定时间间隔采集流体图像序列;然后采用像差修正的直接线性变换法对图像进行光学畸变校正,并沿测流方向设置一组测速线,每条测速线均以像素为横坐标、以时间为纵坐标合成一幅时空图像;接下来对边缘检测后的时空图像进行频谱分析,并采用离散逼近法求解频谱主方向得到与其垂直的纹理主方向;最后根据物方坐标和像方坐标间的映射关系求解测速线的物理尺度,并结合纹理主方向和时间间隔计算一维流速矢量的大小,实现时均流场定标。
所述的流体图像序列为顺序采集的M幅8位灰度图像,M在数值上等于单次时均流场测量的时长T和每帧图像采集的时间间隔Δt之比取整;
所述的一组测速线由L条测速线lij组成,下标i=1,2,…,M表示图像在序列中的编号,下标j=1,2,…,L表示测速线的编号,测速线的宽度为1个像素,测速线的长度为N个像素,N不小于水流示踪物在连续图像间运动位移的2倍;
所述的时空图像的大小由灰度图像的数量M和测速线的长度N决定,为M×N像素,数量和测速线的数量相等,共有L幅;
所述的边缘频谱分析,首先采用Canny算子对时空图像进行边缘检测得到二值化的边缘图像,然后依次进行二维离散傅立叶变换、中心平移、对数变换和灰度拉伸得到直角坐标系下的频谱图像;
所述的频谱主方向定义为经过频谱图像中心的直线中线积分峰值对应的直线角度;
所述的离散逼近法包括如下步骤:
1)以直角坐标为(M/2,N/2)的点为原点、M和N中较小值的一半为半径R=min(M/2,N/2)建立极坐标系;
2)在极坐标系中设置一条过原点且半径为R的搜索线,对于固定的角度γ计算线积分:
3)在0~180°的半圆空间内搜索F(γ)的最大值,首先以1°为步进进行粗定位,再以0.1°为步进在以最大值为中心±0.5°的范围内进行精搜索,并将F(γ)最大值对应的角度γj定义为频谱主方向;最后根据纹理主方向与频谱主方向垂直的关系计算时空图像sTIj的纹理主方向θj,如下:
所述的物方坐标和像方坐标间的映射关系采用像差修正的直接线性变换公式描述如下:
其中,(X,Y,Z)表示物方坐标,(x,y)表示像方坐标,11个未知量li称为直接线性变换系数,光学系统的总像差(δx,δy)通过径向畸变、离心畸变及薄棱镜畸变这三种畸变分量的线性叠加来建模:
其中,xd、yd分别为像点辐射距的水平和垂直分量,k1、k2为径向畸变系数,p1、p2为离心畸变系数,s1、s2为薄棱镜畸变系数。
所述的时均流场定标方法,首先将测速线上首尾像素的物理距离除以像素距离得到单位物理尺度sj:
接下来根据时空图像中位移、时间和流速间的三角关系,对时均流速矢量uj进行定标:
反复执行边缘频谱分析及时均流场定标的步骤直到读入全部L幅时空图像,最后得到时间T内测得的一维时均流速场。
本发明具有以下有益效果:
1、鲁棒性强。利用水流示踪物时空分布的连续性从时空图像中直接估计出一维时均流场,这种时空联合的整体估计方法可以有效克服示踪物非刚性运动引起的不确定性。
2、分辨率高。观测窗口为单像素的测速线,相比空域相关匹配法中的矩形窗口,有效改善了相机在倾斜视角下拍摄引起的图像远场空间分辨率降低,使测量结果更具有代表性。
3、实时性强。方法将流体一维运动矢量估计的问题转化为时空图像频谱主方向的检测问题来解决,由于无需进行穷尽搜索的空域相关运算,计算量仅为传统PIV方法的十分之一。
鉴于以上特点,本方法能够用于河流水面等天然示踪水流的时均流场测量,特别适用于高洪期河流流量的在线监测。
附图说明
图1是本发明涉及的运动矢量估计方法的工作流程图。
具体实施方式
下面结合附图对本发明的技术方案进行详细说明:
本发明的运动矢量估计方法的工作流程如图1所示。共分为以下四个主要步骤:
(1)图像序列采集
首先设置单次时均流场测量的时长T和每帧图像采集的时间间隔Δt。T应当足够大以消除流速脉动效应引起的随机测量误差,同时应当足够小以满足测量的时间分辨率;Δt根据拍摄视场、图像分辨率和待测流进行选取,应当足够大以保证结果的离散化精度,同时应当足够小以保证水流示踪物在图像序列中的分布满足连续性。然后计算图像序列的长度M=T/Δt。最后采用单台相机以时间间隔Δt连续采集M幅8位灰度图像GRIi,i=1,2,…,M表示图像在序列中的编号。
(2)时空图像合成
首先在图像中沿测速方向设置L条测速线lij,j=1,2,…,L表示测速线在图像中的编号;测速线的宽度为1个像素,长度N不小于水流示踪物在连续图像间运动位移的2倍。然后以图像序列的长度M为纵坐标、以测速线的长度N为横坐标,建立L幅大小为M×N像素的时空图像STIj。接下来依次读入一幅采集的灰度图像GRIi,并利用现场布设的控制点采用像差修正的直接线性变换法对图像进行光学畸变校正。最后将无失真图像GRIi中每条测速线lj上的像素载入对应的时空图像STIj,直到读入M幅灰度图像,得到合成的L幅时空图像。
(3)边缘频谱分析
首先依次读入一幅合成的时空图像STIj,采用Canny算子对其进行边缘检测,得到二值化的边缘图像EDIj。然后对EDIj依次进行二维离散傅立叶变换、中心平移、对数变换和灰度拉伸得到直角坐标系下的频谱图像FSIj。接下来将过中心直线中FSIj线积分峰值对应的直线角度定义为频谱主方向γj,采用离散逼近法搜索γj,包括如下步骤:
1)以直角坐标为(M/2,N/2)的点为原点、M和N中较小值的一半为半径R=min(M/2,N/2)建立极坐标系;
2)在极坐标系中设置一条过原点且半径为R的搜索线,对于固定的角度γ计算线积分:
3)在0~180°的半圆空间内搜索F(γ)的最大值,首先以1°为步进进行粗定位,再以0.1°为步进在最大值附近±0.5°的范围内进行精搜索,并将F(γ)最大值对应的角度γj定义为频谱主方向;最后根据纹理主方向与频谱主方向垂直的关系计算时空图像STIj的纹理主方向θj,如下:
(4)时均流场定标
首先根据以下物方坐标(X,Y,Z)和像方坐标(x,y)间的映射关系计算测速线lj上的首尾像素(x1,y1)和(xN,yN)的物方坐标(X1,Y1,Z1)和(XN,YN,ZN):
其中,11个未知量li称为直接线性变换系数,光学系统的总像差(δx,δy)通过径向畸变、离心畸变及薄棱镜畸变这三种畸变分量的线性叠加来建模:
其中,xd、yd分别为像点辐射距的水平和垂直分量,k1、k2为径向畸变系数,p1、p2为离心畸变系数,s1、s2为薄棱镜畸变系数。
然后将测速线上首尾像素的物理距离除以像素距离得到单位物理尺度sj:
接下来根据时空图像中位移、时间和流速间的三角关系,对时均流速矢量uj进行定标:
反复执行步骤(3)和(4)直到读入全部L幅时空图像,最后得到时间T内测得的一维时均流速场。
Claims (5)
1.一种基于时空图像的流体运动矢量估计方法,分为图像序列采集、时空图像合成、边缘频谱分析及时均流场定标四个步骤,其特征在于:方法首先采用单台相机以固定时间间隔采集流体图像序列;然后采用像差修正的直接线性变换法对图像进行光学畸变校正,并沿测流方向设置一组测速线,每条测速线均以像素为横坐标、以时间为纵坐标合成一幅时空图像;接下来对边缘检测后的时空图像进行频谱分析,并采用离散逼近法求解频谱主方向得到与其垂直的纹理主方向;最后根据物方坐标和像方坐标间的映射关系求解测速线的物理尺度,并结合纹理主方向和时间间隔计算一维流速矢量的大小,实现时均流场定标;
所述的流体图像序列为顺序采集的M幅8位灰度图像,M在数值上等于单次时均流场测量的时长T和每帧图像采集的时间间隔Δt之比取整;
所述的一组测速线由L条测速线lij组成,下标i=1,2,…,M表示图像在序列中的编号,下标j=1,2,…,L表示测速线的编号,测速线的宽度为1个像素,测速线的长度为N个像素,N不小于水流示踪物在连续图像间运动位移的2倍;
所述的时空图像的大小由灰度图像的数量M和测速线的长度N决定,为M×N像素,数量和测速线的数量相等,共有L幅;
所述的频谱主方向定义为经过频谱图像中心的直线中线积分峰值对应的直线角度。
2.根据权利要求1所述的基于时空图像的流体运动矢量估计方法,其特征在于:所述的边缘频谱分析首先采用Canny算子对时空图像进行边缘检测得到二值化的边缘图像,然后依次进行二维离散傅立叶变换、中心平移、对数变换和灰度拉伸得到直角坐标系下的频谱图像。
3.根据权利要求1所述的基于时空图像的流体运动矢量估计方法,其特征在于:所述的离散逼近法包括如下步骤:
1)以直角坐标为(M/2,N/2)的点为原点、M和N中较小值的一半为半径R=min(M/2,N/2)建立极坐标系;
2)在极坐标系中设置一条过原点且半径为R的搜索线,对于固定的角度γ计算线积分:
3)在0~180°的半圆空间内搜索F(γ)的最大值,首先以1°为步进进行粗定位,再以0.1°为步进在以最大值为中心±0.5°的范围内进行精搜索,并将F(γ)最大值对应的角度γj定义为频谱主方向;最后根据纹理主方向与频谱主方向垂直的关系计算时空图像STIj的纹理主方向θj,如下:
4.根据权利要求1所述的基于时空图像的流体运动矢量估计方法,其特征在于:所述的物方坐标和像方坐标间的映射关系采用像差修正的直接线性变换公式描述如下:
其中,(X,Y,Z)表示物方坐标,(x,y)表示像方坐标,11个未知量li称为直接线性变换系数,光学系统的总像差(δx,δy)通过径向畸变、离心畸变及薄棱镜畸变这三种畸变分量的线性叠加来建模:
其中,xd、yd分别为像点辐射距的水平和垂直分量,k1、k2为径向畸变系数,p1、p2为离心畸变系数,s1、s2为薄棱镜畸变系数。
5.根据权利要求3所述的基于时空图像的流体运动矢量估计方法,其特征在于:所述的时均流场定标为首先将测速线上首尾像素的物理距离除以像素距离得到单位物理尺度sj:
接下来根据时空图像中位移、时间和流速间的三角关系,对时均流速矢量uj进行定标:
θj如公式(2)所示,反复执行边缘频谱分析及时均流场定标的步骤直到读入全部L幅时空图像,最后得到时间T内测得的一维时均流场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410185108.2A CN103996171B (zh) | 2014-05-05 | 2014-05-05 | 基于时空图像的流体运动矢量估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410185108.2A CN103996171B (zh) | 2014-05-05 | 2014-05-05 | 基于时空图像的流体运动矢量估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103996171A CN103996171A (zh) | 2014-08-20 |
CN103996171B true CN103996171B (zh) | 2017-02-15 |
Family
ID=51310328
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410185108.2A Active CN103996171B (zh) | 2014-05-05 | 2014-05-05 | 基于时空图像的流体运动矢量估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103996171B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104297516B (zh) * | 2014-11-06 | 2017-07-18 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种流体表面二维流速场测量方法 |
CN104777327B (zh) * | 2015-03-17 | 2018-03-20 | 河海大学 | 基于激光辅助标定的时空图像测速系统和方法 |
CN107064546B (zh) * | 2017-05-10 | 2019-09-10 | 东南大学 | 一种基于非圆信号频率估计算法的流体速度估计方法 |
CN108548937A (zh) * | 2018-06-21 | 2018-09-18 | 王开全 | 一种明渠视频流速校正仪 |
CN109977829A (zh) * | 2019-03-18 | 2019-07-05 | 江河瑞通(北京)技术有限公司 | 基于图像识别分析技术的流速检测方法及装置 |
CN110530441B (zh) * | 2019-08-21 | 2020-11-03 | 河海大学 | 一种全天候在线视频测流系统 |
CN111062978B (zh) * | 2019-11-27 | 2022-02-01 | 武汉大学 | 基于频域滤波技术的时空图像测流的纹理识别方法 |
CN112697657B (zh) * | 2021-03-24 | 2021-09-14 | 杭州电子科技大学 | 一种基于气溶胶颗粒物图像的智能测风系统 |
CN113343923A (zh) * | 2021-07-01 | 2021-09-03 | 江苏舆图信息科技有限公司 | 一种基于视频图像的河道排口排水状态实时识别方法 |
CN114397476B (zh) * | 2021-11-15 | 2022-10-14 | 河海大学 | 面向频域时空图像测速的流速有效性识别及修正方法 |
CN115456916B (zh) * | 2022-11-11 | 2023-03-24 | 武汉大水云科技有限公司 | 一种基于图像质量的表面流速计算方法 |
CN116342656B (zh) * | 2023-03-29 | 2024-04-26 | 华北电力大学 | 一种基于自适应边缘检测的时空图像测速方法和装置 |
CN116679080A (zh) * | 2023-05-30 | 2023-09-01 | 广州伏羲智能科技有限公司 | 一种河流表面流速确定方法、装置及电子设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101039183B1 (ko) * | 2011-04-26 | 2011-06-03 | 한국해양대학교 산학협력단 | 보간이 필요없는 어파인변환 기반 하이브리드 피아이브이를 이용한 유동장 계측 시스템 |
KR101305305B1 (ko) * | 2012-06-12 | 2013-09-06 | 동의대학교 산학협력단 | 시공간 영상의 상관 분석을 이용한 표면 유속 측정 시스템 및 방법 |
CN103578118A (zh) * | 2013-10-24 | 2014-02-12 | 河海大学 | 基于序贯图像矢量平均的时均流场重建方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101464471B (zh) * | 2009-01-15 | 2010-09-29 | 清华大学 | 实体模型表面流场实时量测系统及其后处理方法 |
US8339584B2 (en) * | 2010-05-21 | 2012-12-25 | Teledyne Technologies Incorporated | Velocity measuring system |
-
2014
- 2014-05-05 CN CN201410185108.2A patent/CN103996171B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101039183B1 (ko) * | 2011-04-26 | 2011-06-03 | 한국해양대학교 산학협력단 | 보간이 필요없는 어파인변환 기반 하이브리드 피아이브이를 이용한 유동장 계측 시스템 |
KR101305305B1 (ko) * | 2012-06-12 | 2013-09-06 | 동의대학교 산학협력단 | 시공간 영상의 상관 분석을 이용한 표면 유속 측정 시스템 및 방법 |
CN103578118A (zh) * | 2013-10-24 | 2014-02-12 | 河海大学 | 基于序贯图像矢量平均的时均流场重建方法 |
Non-Patent Citations (4)
Title |
---|
A 2D Barcode Extraction Method Based on Texture Direction Analysis;Huaqiao Hu 等;《International Conference on Image and Graphics》;20090923;第759-762页 * |
A novel line detection method in space-time images for microvascular blood flow analysis in sublingual microcirculatory videos;Sheng You 等;《IEEE 10th International Symposium on Biomedical Imaging》;20130411;第828-831页 * |
悬移质运动中水沙两相流的流场测量技术研究综述;陈诚 等;《水利水电科技进展》;20111231;第31卷(第6期);第80-84页 * |
近红外成像的便携式大尺度粒子图像测速表;张振 等;《仪器仪表学报》;20121231;第33卷(第12期);第2840-2850页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103996171A (zh) | 2014-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103996171B (zh) | 基于时空图像的流体运动矢量估计方法 | |
Patalano et al. | Rectification of Image Velocity Results (RIVeR): A simple and user-friendly toolbox for large scale water surface Particle Image Velocimetry (PIV) and Particle Tracking Velocimetry (PTV) | |
Carlier et al. | Experimental study of eddy structures in a turbulent boundary layer using particle image velocimetry | |
CN101297762B (zh) | 在医用诊断超声中的血流特征成像 | |
Ettema et al. | Particle-image velocimetry for whole-field measurement of ice velocities | |
Mikheev et al. | A method for measuring the dynamics of velocity vector fields in a turbulent flow using smoke image-visualization videos | |
Novara et al. | Motion tracking-enhanced MART for tomographic PIV | |
CN108195933B (zh) | 一种检测建筑外墙质量缺陷的检测系统 | |
Sankowski et al. | Estimation of measurement uncertainty in stereo vision system | |
KR101305305B1 (ko) | 시공간 영상의 상관 분석을 이용한 표면 유속 측정 시스템 및 방법 | |
CN108489498A (zh) | 一种基于最大互相关熵无迹粒子滤波的auv协同导航方法 | |
CN104200236A (zh) | 基于dpm的快速目标检测方法 | |
CN108007401A (zh) | 一种基于船载InSAR平台的河湖库沿岸形变检测装置及方法 | |
CN103063143A (zh) | 基于图像识别的测量方法及其系统 | |
Gesemann | From particle tracks to velocity and acceleration fields using B-splines and penalties | |
CN109900331A (zh) | 基于深度相机特征跟踪的河流流量在线检测方法 | |
CN104268880A (zh) | 基于特征和区域匹配相结合的深度信息获取方法 | |
Yeh et al. | Applying adaptive LS-PIV with dynamically adjusting detection region approach on the surface velocity measurement of river flow | |
Herlin et al. | Divergence-free motion estimation | |
Fischer et al. | Acoustic tomography of the atmosphere: Comparison of different reconstruction algorithms | |
CN103793924A (zh) | 一种基于fht-cc的流场图像自适应运动矢量估计方法 | |
Schröder et al. | Near-wall turbulence characterization using 4D-PTV “Shake-The-Box” | |
Douxchamps et al. | Stereoscopic and velocimetric reconstructions of the free surface topography of antidune flows | |
Cai et al. | Estimating small structural motions from multi-view video measurement | |
Lee et al. | The measurement of discharge using a commercial digital video camera in irrigation canals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |