CN1038613C - 河工模型试验流场量测中的图像摄取***及图象处理方法 - Google Patents
河工模型试验流场量测中的图像摄取***及图象处理方法 Download PDFInfo
- Publication number
- CN1038613C CN1038613C CN94119429A CN94119429A CN1038613C CN 1038613 C CN1038613 C CN 1038613C CN 94119429 A CN94119429 A CN 94119429A CN 94119429 A CN94119429 A CN 94119429A CN 1038613 C CN1038613 C CN 1038613C
- Authority
- CN
- China
- Prior art keywords
- image
- cloth
- width
- flow field
- point
- 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.)
- Expired - Fee Related
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000005259 measurement Methods 0.000 title claims abstract description 11
- 239000000700 radioactive tracer Substances 0.000 claims abstract description 20
- 238000003672 processing method Methods 0.000 claims abstract description 11
- 238000006243 chemical reaction Methods 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 6
- 239000004744 fabric Substances 0.000 claims description 40
- 239000002245 particle Substances 0.000 claims description 21
- 238000005086 pumping Methods 0.000 claims description 8
- 238000013461 design Methods 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 238000006116 polymerization reaction Methods 0.000 claims description 4
- 238000007639 printing Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 2
- 238000000151 deposition Methods 0.000 claims description 2
- 230000033001 locomotion Effects 0.000 abstract description 6
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000009991 scouring Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 abstract 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 10
- 238000005070 sampling Methods 0.000 description 4
- 230000002596 correlated effect Effects 0.000 description 2
- 239000013049 sediment Substances 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 239000004411 aluminium Substances 0.000 description 1
- 239000013068 control sample Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000009432 framing Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 238000004080 punching Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
- Treatment Of Fiber Materials (AREA)
Abstract
本发明涉及河工模型试验流场量测中的图像摄取***及图像处理方法。该方法是在流场中施放示踪颗粒,采用摄、录像的方法记录颗粒的运动轨迹,经过图像处理得出全试验段的表面流速和流志。摄像***包括摄像机、录像机、控制***和信号传输采集***。分析处理程序包括图像采集、示踪点定位、比尺换算、图象定位定向等。本发明的图象处理方法不干扰流场,不用固定床面,得到的流态图能真实反映试验段的流态,具有较高的精度,对保证模型试验的冲淤相似具有重要意义。
Description
本发明涉及河工模型试验流场量测中的图像摄取***及图象处理方法,属水利量测技术领域。
现有的河工模型试验中,一般采用多点旋桨流速仪量测流速,如代表我国先进水平的长江科学院和南京水科院的三峡坝区泥沙模型均用此法。两家模型用这一方法量测全河段的流速需1~2天时间,而模型容许量测的时间只有1~2小时。长江科学院采用停水固定地形后量测,而固定一次地形至少需要15天时间,耗资约10万元。此外,当泥沙含量较高时,旋桨流速仪的精度很低,甚至难以应用。
另一种测速方法是浮标示踪人工计时或拍照后分析处理,这一方法费时费力且精度不高,难以用于大面积水域的流速量测。
本发明的目的是开发和研制河工模型试验流场量测中的图象摄取装置及图象处理方法,用以实现在短时段内瞬时同步大面积动态量测试验河段的流速和流态。
本发明设计的图象摄取***,包括流场示踪装置、摄像机、控制台、录像机、监视器、计算机图像版和打印机。流场示踪装置为一河工模型和投球架。摄像机置于河工模型上方的云台上,图像信号通过传输电缆进入控制台,由监视器监视并用录像机记录图像信号。该图像送入计算机图像板进行处理后,打印机输出流速流态图。
本发明的图象处理方法包括下列步骤:
(1)用摄像机和录像机记录球粒即示踪颗粒在水中运动的轨迹。
(2)用图象板对录像机的图象作A/D转换,将每幅图像分成512×512个坐标点,每个点的灰度分成256个等级,全黑的灰度为零,全白为255。
(3)采集信号,从上述A/D转换后的图像中,采集四幅,分别存在图像板的四个帧存体中,或用采一存循环方式,即采一幅图后,存入主机内存后再采集下幅图像,一次采集12幅以内的图幅。
(4)在上述采集到的图像中,求出每个示踪颗粒的形心:先取一幅图像中的一条横线,当遇到第一个有效点时,查找其后的连续有效点,求出这些点的平均横坐标,将512条横线逐一处理,得出示踪颗粒横向坐标聚合后的结果,再从第一条横线上第一个示踪点坐标起,求该点与第二条线以下所有点的距离,当所求距离小于某一限定值时,即认为是同一个颗粒,求出同一颗粒各点的平均坐标作为该点的坐标值即示踪颗粒的形心,将已认定过的坐标全部消去。
(5)求出同一示踪颗粒在二幅图中的位移量:采用上一幅图中某一颗粒与下一幅图的所有颗粒相关,求出最小距离的点即认为是同一个颗粒在下幅图中的位置,从而求出位移量。
(6)剔除错误信号:预估计水流的最大流速,将大于该值的信号删去;
(7)比尺标定和坐标变换:在作图象处理时,先求出两个基准点的坐标,由此定出本幅图的位置和比尺,从而得出各示踪颗粒的速度大小和方向,再通过平移和旋转,将量测的结果转换到地形图的坐标体系中;
(8)将上述图形打印、输出,即得流场的流速流态图。
附图说明:
图1是图象摄取装置的结构框图。
图2是摄取装置中所用的投球架的结构示意图。
图3是图2的侧示图。
图4是图象处理方法的流程框图。
图5是用本发明的图像处理方法得到的流场流态图。
图6是用本发明的图像处理方法得到的流场流速图。
下面结合附图,详细介绍本发明的内容。
为了确定所测流场的具***置,需要预先在河工模型上布设基准点,见图1左下图中的园点,在录像的时候,保证每幅图至少要有两上基准点。
对某些水流情况复杂而又需重点研究的河段,在量测时全水面均匀撒上示踪球,用摄像机记录球粒的运动轨迹,见图1左下图中左边的箭头。
为测断面的流速分布而专门制做了投球架,其结构见图2和图3,其释放球粒在水中运动的轨迹见图1左下图中右边的箭头。图2和图3中,1是球槽,2是球孔,3是抽板,4是底部穿孔。在型铝上按要求等间距打孔,底部加活动抽板。使用时将抽板关上,在球槽和每个孔内装入球粒,准备好后拉动抽板,抽板孔与投球架孔对齐时,全部球粒同时落入水中,并随水流运动,同时用录像机记录球粒运动的轨迹。使用投球架时,应将其放在模型上的已知位置以作坐标变换。
示踪球用直径3.0cm的空心塑料球做成。球内加入一定数量的重沙,使球体有一半沉在水中,经过水槽试验检验,其与水流的跟随性很好,与在水面用纸屑示踪的同步观测结果对比,差别在2.0%以内。因试验水体呈乳白色,故示踪球做成黑色。
在流场中施放示踪颗粒,采用摄、录像的方法记录颗粒的运动轨迹,经过图象处理得出全试验段的表面流速和流态,图象摄取***的结构框图见图1。如设置多套摄录像***,则可在要求的时间内完成流场信号的同步记录,再通过图像处理,即能获得需要的结果。
摄录像和图像处理的硬件配套***包括摄像机、录像机、控制***和信号传输***。
图像采集***由图像板和486微机组成,图像板通过A/D转换将一幅图分成512×512个坐标点,每个点的灰度分成256个等级,全黑的灰度为零,全白为255。
分析计算机软件包括图像采集、示踪点聚合、流线相关、错误信号的剔除、比尺换算、图像定位定向、流速分布、流向流态、任意方位的流速剖面、流态图打印等一系列子程序,逻辑流程见图4。这一***能高效、快速地给出试验河段全流场的流速流态分布图。
信号采集指令由图像板经销商提供。将一幅图像转换成512×512个点阵、256级灰度的数字信号。采样程序用C语言编写,调用有关的采样指令。采样分两种方式,一种是连续采集四幅图像,分别存在图像板的四个帧存体中,然后同逐一取出进行处理。这一方法速度快,可在1/25秒内采一幅图。另一种是采一存循环方式,即采一幅图后,存入主机内存后再采集下一幅图像,这种方式能一次采集12幅以内的图幅,但速度较慢,采一幅图需0.6秒。为了调整采样频率,可在程序中加一段瞎时程序来控制。
在一幅图的512×512的坐标点中,只有示踪颗粒占据的点才是有用的信号。对于黑色的示踪颗粒,只需保留灰度小于某一临界值的那些坐标点作为信号的有效点,这一临界值因摄像时的光线明暗会有所不同。在一幅图中,一个示踪颗粒可能占有几个坐标点,这时需要求出这几个点的形心以代表这一颗粒的准确位置,这一过程可用距离相关来完成:
取一幅图像中的一条横线,当遇到第一个有效点时,查找其后的连续有效点,求出这些点的平均横坐标,将512条横线逐一处理,得出示踪颗粒横向坐标聚合后的结果。
从第一条横线上第一个示踪点坐标起,求该点与第二条线以下所有点的距离,当所求距离小于某一限定值时,即认为是同一个颗粒,求出同一颗粒各点的平均坐标作为该点的坐标值,将已认定过的坐标全部消去。将512条横线按顺序处理,最后得出全部示踪点的坐标。这样,一幅图中有多少个示踪颗粒,就只有多少个坐标点。
控制采样频率,保证上一幅图中的颗粒在下一幅图中只移动了一个微小的距离,采用上幅图中某一颗粒与下一幅图的所有颗粒相关,求出最小距离的点即认为是同一个颗粒在下幅图中的位置,从而求出位移量。
当认定下一幅图中的某一点后,将其坐标消去,以避免求下一个颗粒的相关时产生混浠。如在一定的距离内没有点,则用前一幅图中本身点的坐标作为下一幅图中的坐标。用计算机时钟定出采样历时,即可求出颗粒的速度。将各连续图幅逐一处理,最后可得出颗粒运动的连续轨迹,
从而定出流速的大小和方向。
在流线相关的过程中,有可能找错点而产生错误信号,需要剔除这些信号,可用以下判别条件进行处理:预估计水流的最大流速,将大于该值的信号删去。
当水流方向大致相同时,将流向偏离总流向的角度大于某一给定值的信号删去。也可以将全图分成儿个部分,分别给出每一部分的总流向进行处理。最后可用人工操作删除明显不合理的点。
在录像的时候,每幅图至少要有两个基准点,在作图像处理时,先求出两个基准点的坐标,由此定出本幅图的位置和比尺。由于摄像时镜头可位于任一方位,采集的图像即为镜头的方位,故需通过平移和旋转,
将量测的结果转换到地形图的坐标体系中。
得出流场图以后,有时需要知道沿某一方向(如航道中心线)的流速分布,这就需要按要求的方向和距离截出流速剖面图。
打印可分图形和表格两种形式,图形比较直观,表格则给出流速的大小、方向和坐标位置。
本发明设计的图象处理装置和处理方法具有以下技术效果:
1、流速量测精度可达2.0mm/s,比常规方法高一个数量级。
2、能真实反映河床的冲淤变化,不干扰流场,不用固定床面,得出的结果能真实地反映试验段的流态,对保证模型试验的冲淤相似具有重要意义。
3、省时省力省经费,全河段的量测时间只需20分钟。可节省大量的试验时间,具有很大的经济效益。
4、摄录像***包括摄像机、录像机、控制***和信号传输采集***,已实现了遥控自动操作,使用方便灵活。
5、软件***功能齐全,能很快得出全流场表面流速的大小和方向。
Claims (3)
1、一种河工模型试验流场量测中的图象摄取***,其特征在于该摄取***包括流场示踪装置、摄像机、控制台、录像机、监视器、计算机图像版和打印机;所述的流场示踪装置为一河工模型和投球架,摄像机置于河工模型上方的云台上,图像信号通过传输电缆进入控制台,由监视器监视并用录像机记录图像信号,该图像送入计算机图像板进行处理后,打印机输出流速流态图。
2、如权利要求1所述的图象摄取***,其特征在于其中所述的投球架由球架和抽板组成,球架为一凹槽,凹槽中开有等间距的球孔,凹槽底部开有穿孔,抽板插在穿孔中,抽板上开有与凹槽中球孔的孔径、间距相等的圆孔。
3、一种河工模型试验流场量测中的图象处理方法,其特征在于该方法包括下列步骤:
(1)用摄像机和录像机记录球粒即示踪颗粒在水中运动的轨迹;
(2)用图象板对录像机的图象作A/D转换,将每幅图像分成512×512个坐标点,每个点的灰度分成256个等级,全黑的灰度为零,全白为255;
(3)采集信号,从上述A/D转换后的图像中,采集四幅,分别存在图像板的四个帧存体中,或用采一存循环方式,即采一幅图后,存入主机内存后再采集下幅图像,一次采集12幅以内的图幅;
(4)在上述采集到的图像中,求出每个示踪颗粒的形心:先取一幅图像中的一条横线,当遇到第一个有效点时,查找其后的连续有效点,求出这些点的平均横坐标,将512条横线逐一处理,得出示踪颗粒横向坐标聚合后的结果,再从第一条横线上第一个示踪点坐标起,求该点与第二条线以下所有点的距离,当所求距离小于某一限定值时,即认为是同一个颗粒,求出同一颗粒各点的平均坐标作为该点的坐标值即示踪颗粒的形心,将已认定过的坐标全部消去;
(5)求出同一示踪颗粒在二幅图中的位移量:采用上一幅图中某一颗粒与下一幅图的所有颗粒相关,求出最小距离的点即认为是同一个颗粒在下幅图中的位置,从而求出位移量;
(6)剔除错误信号:预估计水流的最大流速,将大于该值的信号删去;
(7)比尺标定和坐标变换:在作图象处理时,先求出两个基准点的坐标,由此定出本幅图的位置和比尺,从而得出各示踪颗粒的速度大小和方向,再通过平移和旋转,将量测的结果转换到地形图的坐标体系中;
(8)将上述图形打印、输出,即得流场的流速流态图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN94119429A CN1038613C (zh) | 1994-12-23 | 1994-12-23 | 河工模型试验流场量测中的图像摄取***及图象处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN94119429A CN1038613C (zh) | 1994-12-23 | 1994-12-23 | 河工模型试验流场量测中的图像摄取***及图象处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1116704A CN1116704A (zh) | 1996-02-14 |
CN1038613C true CN1038613C (zh) | 1998-06-03 |
Family
ID=5039230
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN94119429A Expired - Fee Related CN1038613C (zh) | 1994-12-23 | 1994-12-23 | 河工模型试验流场量测中的图像摄取***及图象处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1038613C (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100343675C (zh) * | 2005-02-01 | 2007-10-17 | 湖南继善高科技有限公司 | 流体内部全流场三维可视流向、流速测量装置 |
CN101408534B (zh) * | 2008-10-28 | 2011-12-07 | 山东科技大学 | 室内定轨移动摄影测量*** |
CN102538763B (zh) * | 2012-02-14 | 2014-03-12 | 清华大学 | 一种河工模型试验三维地形的量测方法 |
CN103149031B (zh) * | 2013-04-02 | 2016-02-10 | 哈尔滨电机厂有限责任公司 | 一种计算机控制的模型水轮机流态观测的同步数字式成像方法 |
CN104374894B (zh) * | 2014-06-13 | 2016-10-26 | 黄河水利委员会黄河水利科学研究院 | 小流域不同地貌单元水蚀过程精细模拟试验方法 |
CN104215795B (zh) * | 2014-09-29 | 2017-07-25 | 南京昊控软件技术有限公司 | 基于昊控表面流场计算算法的大尺度粒子图像测速方法 |
CN104535126B (zh) * | 2015-01-16 | 2018-01-26 | 中国水利水电科学研究院 | 一种河道枯季流量测量***和方法 |
CN106908218B (zh) * | 2017-03-23 | 2018-11-30 | 北京尚水信息技术股份有限公司 | 河道表面流场的处理方法 |
CN107702885B (zh) * | 2017-09-29 | 2019-08-30 | 中国水利水电科学研究院 | 床面颗粒起悬时间的提取方法 |
CN109782019B (zh) * | 2019-03-15 | 2020-04-24 | 中国科学技术大学 | 大气污染物二维运动速度测量方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2569974A (en) * | 1944-04-11 | 1951-10-02 | United Aircraft Corp | Velocity measuring device |
EP0436125A1 (en) * | 1989-12-02 | 1991-07-10 | Ezel Inc. | Flowing velocity measurement system |
EP0505099A2 (en) * | 1991-03-22 | 1992-09-23 | Erin Technologies Inc. | Method and apparatus for determining the flow rate of a viscous fluid stream |
-
1994
- 1994-12-23 CN CN94119429A patent/CN1038613C/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2569974A (en) * | 1944-04-11 | 1951-10-02 | United Aircraft Corp | Velocity measuring device |
EP0436125A1 (en) * | 1989-12-02 | 1991-07-10 | Ezel Inc. | Flowing velocity measurement system |
EP0505099A2 (en) * | 1991-03-22 | 1992-09-23 | Erin Technologies Inc. | Method and apparatus for determining the flow rate of a viscous fluid stream |
Also Published As
Publication number | Publication date |
---|---|
CN1116704A (zh) | 1996-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1038613C (zh) | 河工模型试验流场量测中的图像摄取***及图象处理方法 | |
CN102279973B (zh) | 基于高梯度关键点的海天线检测方法 | |
CN110208859B (zh) | 基于超声成像的油基泥浆井裂缝定量参数智能计算方法 | |
CN107621435A (zh) | 一种骨料在线监测装置 | |
CN100554969C (zh) | 湍动条件下低浓度污泥沉降速率的测量方法和装置 | |
CN109711353B (zh) | 一种基于机器视觉的船舶吃水线区域识别方法 | |
CN106989683A (zh) | 一种盾构机盾尾间隙视觉测量方法 | |
Fox et al. | Detection of changes in ridge‐crest morphology using repeated multibeam sonar surveys | |
CN111476159A (zh) | 一种基于双角回归的检测模型训练、检测方法及装置 | |
CN107462741A (zh) | 一种运动物体速度及加速度测量装置 | |
CN101692070A (zh) | 超声扫描显微镜b扫描图像的构建方法 | |
Li et al. | Large-scale particle tracking velocimetry with multi-channel CCD cameras | |
CN112418028A (zh) | 一种基于深度学习的卫星图像中船舶识别与分割方法 | |
CN107631960A (zh) | 海洋物质输移扩散系数ccd自动监测实验装置及其监测方法 | |
CN107103283B (zh) | 一种sar图像舰船目标几何特征并行提取方法与装置 | |
CN1141553C (zh) | 大范围表面流速场的图像处理***及其同步实时量测方法 | |
Ebersole et al. | DUCK85 photopole experiment | |
CN112834780A (zh) | 红外海藻粒子图像测速装置及其测量方法 | |
CN106991682B (zh) | 自动港口货船的提取方法及装置 | |
CN101852813B (zh) | 一种泥沙沉速的测量装置及其测量方法 | |
CN114155428A (zh) | 基于Yolo-v3算法的水下声呐侧扫图像小目标检测方法 | |
USH1409H (en) | Optical correlation velocity log | |
CN111965628A (zh) | 一种垂直出水航行体瞬时波浪参数估计方法 | |
CN205506347U (zh) | 一种基于高清摄影的薄层水流滚波测量装置 | |
Chen et al. | Research on contour extraction and edge detection of sugarcane image based on matlab |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C06 | Publication | ||
PB01 | Publication | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |