CN107424168B - 连续曝光式的表面流场测量方法 - Google Patents

连续曝光式的表面流场测量方法 Download PDF

Info

Publication number
CN107424168B
CN107424168B CN201710120239.6A CN201710120239A CN107424168B CN 107424168 B CN107424168 B CN 107424168B CN 201710120239 A CN201710120239 A CN 201710120239A CN 107424168 B CN107424168 B CN 107424168B
Authority
CN
China
Prior art keywords
linear information
flow field
information
image
continuous exposure
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
Application number
CN201710120239.6A
Other languages
English (en)
Other versions
CN107424168A (zh
Inventor
罗朝林
林俊
陈若舟
邢方亮
陈俊
王磊
丘瑾炜
王天奕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Pearl River Hydraulic Research Institute of PRWRC
Original Assignee
Pearl River Hydraulic Research Institute of PRWRC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Pearl River Hydraulic Research Institute of PRWRC filed Critical Pearl River Hydraulic Research Institute of PRWRC
Priority to CN201710120239.6A priority Critical patent/CN107424168B/zh
Publication of CN107424168A publication Critical patent/CN107424168A/zh
Application granted granted Critical
Publication of CN107424168B publication Critical patent/CN107424168B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M10/00Hydrodynamic testing; Arrangements in or on ship-testing tanks or water tunnels
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/26Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting optical wave
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Fluid Mechanics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)

Abstract

本发明涉及表面流场测量方法的技术领域,公开了连续曝光式的表面流场测量方法,包括:1)、撒布示踪物;2)、连续曝光,分别形成第一重影图像以及第二重影图像;3)、得到第一线性信息为∑nLines(Start,End,Width)以及第二线性信息为∑n+1Lines(Start,End,Width);4)、以离图像中心点距离升序排列,分别得到∑nSortLines(Start,End,Width)以及∑n+ 1SortLines(Start,End,Width);5)、根据最小领域法,将搜索结果保存在
Figure DDA0001297236180000011
6)、计算速度大小为
Figure DDA0001297236180000012
计算速度方向为
Figure DDA0001297236180000013
计算流场为

Description

连续曝光式的表面流场测量方法
技术领域
本发明涉及表面流场测量方法的技术领域,尤其涉及连续曝光式的表面流场测量方法。
背景技术
表面流场测量技术是基于传统PIV技术衍生系统,把微观的二维/三维流场测量拓展到宏观的二维表面流速场测量,在水力学物理模型领域需求较广。传统技术中,采用自然光或者白炽灯作为光源,以现场水体上漂浮的示踪物作为测量对象,把示踪物的流动速度转换成水体的流动速度,均匀撒布合适规格的示踪物(较好的跟随性、较好的区别度等)则可计算出水体的二维表面流场。
近年来,水力学物理模型表面流场测量技术取得了较快的发展,目前的测量系统一般采用基于灰度值相关的粒子图像测速技术(PIV,Particle Image Velocimetry)和粒子图像跟踪测速技术(PTV,Particle Tracking Velocimetry)。这两种传统的测量方法都没有同时很好地解决高频噪声和低频噪声的干扰。另外,在目前的模型试验场地一般存在光照不均、电磁噪声严重、水面漫反射等不利因素,导致测量数据存在较多误矢量以及测量结果不准确的问题,需要较多的人工干预,从而制约了传统测量技术的发展及推广应用。
针对当前表面流速测量常用的互相关算法存在多峰、多值且易受外界噪声干扰,网格划分不准,粒子分布不均匀,粒子拍摄不完全,提出连续曝光下的新型表面流场测量方法,通过计算连续曝光形成示踪粒子的长度来计算流速大小,运用多帧首尾最近法和移动趋势法确定流速的方向,减少了计算的误矢量,提高了模型试验流场测量的精度。
发明内容
本发明的目的在于提供连续曝光式的表面流场测量方法,旨在解决现有技术中的表面流场测量方法存在较多误矢量以及测量结果不准确的问题。
本发明是这样实现的,连续曝光式的表面流场测量方法,包括以下步骤:
1)、在水力物理模型的水体中撒布示踪物;
2)、在第n时刻和第n+1时刻,利用摄像机对水体的表面进行连续曝光,采集多帧重影图像,分别形成第一重影图像以及第二重影图像;
3)、对第一重影图像以及第二重影图像分别进行图像处理以及Hough变换,提取直(曲)线信息,分别得到第一线性信息为∑nLines(Start,End,Width)以及第二线性信息为∑n+1Lines(Start,End,Width);
4)、分别对所述第一线性信息以及第二线性信息进行排序,以离图像中心点距离升序排列,分别得到第一排列线性信息为∑nSortLines(Start,End,Width),以及第二排列线性信息为∑n+1SortLines(Start,End,Width);
5)、根据最小领域法,将第一排列线性信息在第二排列线性信息中进行搜索,并将搜索结果保存在
Figure BDA0001236797660000021
6)、计算速度大小为
Figure BDA0001236797660000022
计算速度方向为
Figure BDA0001236797660000023
计算流场为
Figure BDA0001236797660000024
进一步地,在所述步骤3)中,对形成的所述第一重影图像及第二重影图像的图像处理包括分别对所述第一重影图像及第二重影图像依序进行滤波处理以及二值化处理。
进一步地,在所述步骤5)的
Figure BDA0001236797660000025
中,当i>1时,则通过示踪物的流向的连续性去剔除,使得i=1。
与现有技术相比,本发明提供的连续曝光式的表面流场测量方法具有以下优点:
1)、通过连续曝光形成第一重影图像以及第二重影图形,构建线性信息库,通过获取重影图像的长度计算流速的大小的标量,解决了单粒子误矢量多的问题;
2)、利用重影图像的长度信息取代粒子中心点信息,减少了粒子光照不均和移动不水平导致图形不完整等外部噪声的影响,提高了计算结果的精度;
3)、连续曝光式的表面流场测量方法摒弃了传统根据网格大小进行相关计算的方法,解决了网格大小不科学引进的计算误差,提高了计算的可靠性,降低了对CCD摄像机的要求,对分辨率和采样频率要求低,性价比较高;
4)、突破了传统依靠搜索单个示踪物的偏移计算速度的方法,转换成通过测量示踪物连续移动形成重影的长度来计算流速,弱化了示踪物形状及边缘和环境光照的要求,减少了窗口设置不佳导致的主观误差,提高了计算的精度,可靠性好。
附图说明
图1是本发明实施例提供的连续曝光式的表面流场测量方法的流程示意图;
图2是本发明实施例提供的连续曝光式的表面流场测量方法的模拟视图;
图3是本发明实施例提供的连续曝光式的表面流场测量方法的模拟视图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
以下结合具体附图对本发明的实现进行详细的描述。
参照图1~3所示,为本发明提供的较佳实施例。
本实施例提供的连续曝光式的表面流场测量方法,用于在水力学物理模型领域中,把微观的二维/三维流场测量拓展到宏观的二维表面流速场测量。
连续曝光式的表面流场测量方法,具体如下步骤:
1)、在水力学物理模型的水体中撒布合适规格的示踪物;
2)、在第n时刻和第n+1时刻,利用摄像机对水体的表面进行连续曝光,采集多帧重影图像,分别形成第一重影图像以及第二重影图像;
3)、对第一重影图像以及第二重影图像分别进行图像处理以及Hough变换,提取直(曲)线信息,分别得到第一线性信息为∑nLines(Start,End,Width)以及第二线性信息为∑n+1Lines(Start,End,Width);
4)、分别对第一线性信息以及第二线性信息进行排序,以离图像中心点距离升序排列,分别得到第一排列线性信息为∑nSortLines(Start,End,Width),以及第二排列线性信息为∑n+1SortLines(Start,End,Width);
5)、根据最小领域法,将第一排列线性信息在第二排列线性信息中进行搜索,并将搜索结果保存在
Figure BDA0001236797660000041
6)、计算速度大小为:
Figure BDA0001236797660000042
计算速度方向为:
Figure BDA0001236797660000043
输出流场
Figure BDA0001236797660000044
上述提供的连续曝光式的表面流场测量方法中,通过计算时刻内连续曝光形成的第一重影图像以及第二重影图像的长度,可计算出该时刻内流速的标量,通过搜索连续多帧图像的首尾起始重影的最小邻域,可计算出该时刻内流速的方向,再通过采用多种粒径的示踪粒子,判断曝光重影的宽度,减少起始点位置判断的误差,可进一步提高计算的准确性。
连续曝光式的表面流场测量方法具有以下优点:
1)、通过连续曝光形成第一重影图像以及第二重影图形,构建线性信息库,通过获取重影图像的长度计算流速的大小的标量,解决了单粒子误矢量多的问题;
2)、利用重影图像的长度信息取代粒子中心点信息,减少了粒子光照不均和移动不水平导致图形不完整等外部噪声的影响,提高了计算结果的精度;
3)、连续曝光式的表面流场测量方法摒弃了传统根据网格大小进行相关计算的方法,解决了网格大小不科学引进的计算误差,提高了计算的可靠性,降低了对CCD摄像机的要求,对分辨率和采样频率要求低,性价比较高;
4)、突破了传统依靠搜索单个示踪物的偏移计算速度的方法,转换成通过测量示踪物连续移动形成重影的长度来计算流速,弱化了示踪物形状及边缘和环境光照的要求,减少了窗口设置不佳导致的主观误差,提高了计算的精度,可靠性好。
本实施例中,在步骤3)中,对形成的第一重影图像及第二重影图像的图像处理包括分别对第一重影图像及第二重影图像依序进行滤波处理以及二值化处理。
具体地,在步骤5)的
Figure BDA0001236797660000051
中,当i>1时,则通过示踪物的流向的连续性去剔除,确保结果的唯一性,也就是确保i=1。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (2)

1.连续曝光式的表面流场测量方法,其特征在于,包括以下步骤:
1)、在水力物理模型的水体中撒布示踪物;
2)、在第n时刻和第n+1时刻,利用摄像机对水体的表面进行连续曝光,采集多帧重影图像,分别形成第一重影图像以及第二重影图像;
3)、对第一重影图像以及第二重影图像分别进行图像处理以及Hough变换,提取直线信息或曲线信息,分别得到第一线性信息为
Figure FDA0002519270590000011
以及第二线性信息为
Figure FDA0002519270590000012
4)、分别对所述第一线性信息以及第二线性信息进行排序,以离图像中心点距离升序排列,分别得到第一排列线性信息为
Figure FDA0002519270590000013
以及第二排列线性信息为
Figure FDA0002519270590000014
5)、根据最小邻域法,将第一排列线性信息在第二排列线性信息中进行搜索,并将搜索结果保存在
Figure FDA0002519270590000015
6)、计算速度大小为
Figure FDA0002519270590000016
计算速度方向为
Figure FDA0002519270590000017
计算流场为
Figure FDA0002519270590000018
在所述步骤5)的
Figure FDA0002519270590000019
中,当i>1时,则通过示踪物的流向的连续性去剔除,使得i=1。
2.如权利要求1所述的连续曝光式的表面流场测量方法,其特征在于,在所述步骤3)中,对形成的所述第一重影图像及第二重影图像的图像处理包括分别对所述第一重影图像及第二重影图像依序进行滤波处理以及二值化处理。
CN201710120239.6A 2017-03-02 2017-03-02 连续曝光式的表面流场测量方法 Active CN107424168B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710120239.6A CN107424168B (zh) 2017-03-02 2017-03-02 连续曝光式的表面流场测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710120239.6A CN107424168B (zh) 2017-03-02 2017-03-02 连续曝光式的表面流场测量方法

Publications (2)

Publication Number Publication Date
CN107424168A CN107424168A (zh) 2017-12-01
CN107424168B true CN107424168B (zh) 2020-09-11

Family

ID=60423089

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710120239.6A Active CN107424168B (zh) 2017-03-02 2017-03-02 连续曝光式的表面流场测量方法

Country Status (1)

Country Link
CN (1) CN107424168B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112734680B (zh) * 2020-12-31 2024-07-05 视涯科技股份有限公司 一种重影测量方法、装置、可读存储介质和计算机设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011058824A (ja) * 2009-09-07 2011-03-24 Kansai Electric Power Co Inc:The 可視化トレーサ
CN102331511A (zh) * 2011-06-17 2012-01-25 河海大学 Piv图像高频采集方法
CN102393473A (zh) * 2011-08-02 2012-03-28 南京理工大学 基于单幅图像的粒子运动速度测试系统
CN103698553A (zh) * 2013-12-26 2014-04-02 天津大学 基于红外摄像的新型表面流场测速系统和测速方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011058824A (ja) * 2009-09-07 2011-03-24 Kansai Electric Power Co Inc:The 可視化トレーサ
CN102331511A (zh) * 2011-06-17 2012-01-25 河海大学 Piv图像高频采集方法
CN102393473A (zh) * 2011-08-02 2012-03-28 南京理工大学 基于单幅图像的粒子运动速度测试系统
CN103698553A (zh) * 2013-12-26 2014-04-02 天津大学 基于红外摄像的新型表面流场测速系统和测速方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
珠江河口整体物理模型流场测量系统研发;罗朝林等;《人民珠江》;20131225(第6期);全文 *

Also Published As

Publication number Publication date
CN107424168A (zh) 2017-12-01

Similar Documents

Publication Publication Date Title
US10682108B1 (en) Methods, systems, and computer readable media for three-dimensional (3D) reconstruction of colonoscopic surfaces for determining missing regions
CN106875415B (zh) 一种动态背景中弱小动目标的连续稳定跟踪方法
CN109059895B (zh) 一种基于手机摄像头和传感器的多模态室内测距及定位方法
CN102789568B (zh) 一种基于深度信息的手势识别方法
CN112258658B (zh) 一种基于深度相机的增强现实可视化方法及应用
CN109791696A (zh) 利用事件摄像机同时进行定位和映射
CN110910447A (zh) 一种基于动静态场景分离的视觉里程计方法
CN110006372B (zh) 一种基于局部优化的三维点云平面度计算方法
CN105578034A (zh) 一种对目标进行跟踪拍摄的控制方法、控制装置及系统
CN110850109B (zh) 一种基于模糊图像测量车速的方法
CN109712247B (zh) 基于混合现实技术的实景训练系统
CN106296736B (zh) 一种仿记忆引导的模式识别方法
CN105912977B (zh) 基于点聚类的车道线检测方法
Li et al. Large-scale particle tracking velocimetry with multi-channel CCD cameras
CN108520550B (zh) 基于噪声分类与mls的点云数据泊松曲面重建方法
CN107424168B (zh) 连续曝光式的表面流场测量方法
US9443312B2 (en) Line parametric object estimation
Yu et al. Correlation analysis of spatio-temporal images for estimating two-dimensional flow velocity field in a rotating flow condition
CN110770793A (zh) 一种三维轨迹数据的异常值处理方法及光学动作捕捉方法
CN114994355A (zh) 一种基于非激光光源的气液两相流流场测量方法和系统
CN110274749A (zh) 基于2维piv的旋流器内部流场测量方法和系统
CN104104911B (zh) 全景图像生成过程中的时间戳消除和重置方法及系统
CN116862922B (zh) 基于图像分割与雷达信息融合的目标定位方法、系统及介质
CN117671033A (zh) 一种基于夜间灯光追踪的相机像主点快速定标方法和系统
CN114972170B (zh) 一种基于鱼眼相机的密集场景下抗遮挡的物体检测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Qiu Jinwei

Inventor after: Chen Ruozhou

Inventor after: Lin Jun

Inventor after: Wang Lei

Inventor after: Xing Fangliang

Inventor after: Chen Jun

Inventor after: Luo Chaolin

Inventor after: Wang Tianyi

Inventor after: Guo Zebin

Inventor before: Luo Chaolin

Inventor before: Lin Jun

Inventor before: Chen Ruozhou

CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Luo Chaolin

Inventor after: Lin Jun

Inventor after: Chen Ruozhou

Inventor after: Xing Fangliang

Inventor after: Chen Jun

Inventor after: Wang Lei

Inventor after: Qiu Jinwei

Inventor after: Wang Tianyi

Inventor before: Qiu Jinwei

Inventor before: Chen Ruozhou

Inventor before: Lin Jun

Inventor before: Wang Lei

Inventor before: Xing Fangliang

Inventor before: Chen Jun

Inventor before: Luo Chaolin

Inventor before: Wang Tianyi

Inventor before: Guo Zebin

GR01 Patent grant
GR01 Patent grant