CN103969024B - 一种双脉冲式泥沙运移信息的图像测量方法 - Google Patents
一种双脉冲式泥沙运移信息的图像测量方法 Download PDFInfo
- Publication number
- CN103969024B CN103969024B CN201410182506.9A CN201410182506A CN103969024B CN 103969024 B CN103969024 B CN 103969024B CN 201410182506 A CN201410182506 A CN 201410182506A CN 103969024 B CN103969024 B CN 103969024B
- Authority
- CN
- China
- Prior art keywords
- image
- imgr
- imgl
- point
- pulse
- 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
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000013049 sediment Substances 0.000 title claims abstract description 10
- 239000013598 vector Substances 0.000 claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 239000000203 mixture Substances 0.000 claims abstract description 6
- 230000000007 visual effect Effects 0.000 claims abstract description 6
- 235000018734 Sambucus australis Nutrition 0.000 claims abstract description 5
- 244000180577 Sambucus australis Species 0.000 claims abstract description 5
- OGFXBIXJCWAUCH-UHFFFAOYSA-N meso-secoisolariciresinol Natural products C1=2C=C(O)C(OC)=CC=2CC(CO)C(CO)C1C1=CC=C(O)C(OC)=C1 OGFXBIXJCWAUCH-UHFFFAOYSA-N 0.000 claims abstract description 5
- 238000009434 installation Methods 0.000 claims abstract description 3
- 238000004458 analytical method Methods 0.000 claims description 35
- 230000008878 coupling Effects 0.000 claims description 15
- 238000010168 coupling process Methods 0.000 claims description 15
- 238000005859 coupling reaction Methods 0.000 claims description 15
- 229910052724 xenon Inorganic materials 0.000 claims description 12
- FHNFHKCVQCLJFQ-UHFFFAOYSA-N xenon atom Chemical compound [Xe] FHNFHKCVQCLJFQ-UHFFFAOYSA-N 0.000 claims description 12
- 238000005286 illumination Methods 0.000 claims description 7
- HUTDUHSNJYTCAR-UHFFFAOYSA-N ancymidol Chemical compound C1=CC(OC)=CC=C1C(O)(C=1C=NC=NC=1)C1CC1 HUTDUHSNJYTCAR-UHFFFAOYSA-N 0.000 claims description 6
- 238000013508 migration Methods 0.000 claims description 4
- 230000005012 migration Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000010191 image analysis Methods 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 7
- 238000002474 experimental method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
一种双脉冲式泥沙运移信息的图像测量方法，属于视觉测量技术领域。两台左右放置的摄像机组成立体视觉系统；其中一台安装点阵光栅片，另一台安装负透镜；在同步控制器的时序控制下，利用点阵光栅片和负透镜脉冲激光器交替打出点阵光与体积椎光，同时利用摄像机的跨帧采集技术在较短的时间内获取不同角度的点阵图像与泥沙图像；通过第一次曝光获取的点阵图像的配对信息来预测第二次曝光泥沙图像的立体匹配位置图像,依照预设定的网格进行归一化协方差相关匹配，最后结合摄像机定标参数重构出泥沙的空间运动矢量；本发明无需停机待时，一次分析便可以得到网格节点处的空间矢量场。
Description
技术领域
本发明属于视觉测量技术领域，涉及到泥沙运移信息的测量问题，特别涉及低纹理目标的三维重建与运动测量方法。
背景技术
在海洋工程中，研究泥沙在风、浪、流作用下的运动以及冲刷时水下设施(如海底管线、水下桩基等)对泥沙的作用时往往需要准确测量泥沙的运移信息。目前常使用的实验方法是：首先在水槽内进行造流、造波，经过一段时间水流与被测目标的相互作用后，停止造流、造波，同时记录时间；然后等水面静止后(为了保证目标不再产生运动)，使用激光扫描法分段对被测目标表面进行重构并记录数据；接下来继续进行造流、造波实验，按上面的操作得到当前时刻的被测目标的表面信息。一直反复如此操作，直到得到所需时间的实验结果。最后，对数据进行三维插值并根据不同时刻的三维信息推测目标表面的空间运动矢量信息。从前面所述的实验情况可以看出，进行一次运动测量需要先后进行“暂停等时”、“激光扫描”、“三维插值”以及“运动矢量合成”等操作。操作步骤较多、运算量较大，因此实验时间较长，效率较低。
综上所述，目前还没有一种高效的三维重构技术可以很好地解决实验水槽中泥沙运移信息的测量问题。
发明内容
针对海洋工程中泥沙运动测量实验过程中因实验等时、计算繁杂、操作步骤较多等原因所引起的实验效率较低的问题，本发明将视觉测量原理与激光技术相结合，提出一种双脉冲式泥沙运移信息的图像测量方法。本发明的技术方案如下：
一种双脉冲式泥沙运移信息的图像测量方法由摄像机、脉冲激光器、和内插图像采集卡的PC机组成；两台左右放置的具有外控双曝光功能的工业摄像机组成立体视觉系统固定在被测地形可以被观测到的地点；两台脉冲激光器配有各自的脉冲激光导光臂用于照明；其中一台安装点阵光栅片，另一台安装负透镜；同步控制器的设置端口与PC机相连，用于设置同步控制器的控制参数，而同步控制器的输出端子分别与摄像机A、摄像机B、激光器A的氙灯与调Q外控端口、激光器B的氙灯与调Q的外控端口相连；在同步控制器的时序控制下，利用点阵光栅片和负透镜，两个脉冲激光器可以交替打出点阵光与体积椎光，对泥沙地形进行照明；同时利用摄像机的跨帧采集技术在较短的时间内获取不同角度的点阵图像与泥沙图像；通过第一次曝光获取的点阵图像的配对信息来预测第二次曝光泥沙图像的立体匹配位置；PC机对不同时刻的图像依照预设定的网格进行归一化协方差相关匹配，从而获取左右像平面上的同名运动矢量，最后结合摄像机定标参数重构出泥沙的空间运动矢量，得到泥沙地形的运移信息；其具体方法如下：
步骤一：设置同步控制器的时序：首先保证摄像机A的外触发时序与摄像机B外触发时序严格一致；其次控制脉冲激光器A的氙灯外控信号与调Q外控触发信号使得脉冲激光器A的脉冲在摄像机A第一次曝光时间间隔内发出，同理控制脉冲激光器B的氙灯外控信号与调Q外控触发信号使得脉冲激光器B的脉冲在摄像机B第二次曝光时间间隔内发出；
步骤二：左右放置两台摄像机，调整角度使得目标区域可以被观测到；调整脉冲激光器A的导光臂和脉冲激光器B的导光臂，使得两台激光器的照明区域均覆盖被测地形；
步骤三：设置两台工业摄像机为同步外控模式，使用张氏定标法进行两台摄像机的定标，完成每台摄像机的内参数与外参数的计算；
步骤四D：设置两台摄像机为双曝光模式；
步骤五：令两台脉冲激光器发出激光脉冲之间的时间间隔为Δt，获取t1时刻的图像，令第一次曝光的左图像为ImgL(t1)右图像为ImgR(t1)，而第二次曝光的左图像为ImgL(t1+Δt)右图像为ImgR(t1+Δt)；
步骤六：分别对ImgL(t1)和ImgR(t1)图像进行二值化处理；然后利用重心法或曲面拟合法对图像中的光斑进行中心提取；之后，在ImgL(t1)图像上建立一个分析窗口，使得该分析窗内包含的光斑数目不小于5个；利用极线约束，寻找分析窗口内每个光斑点在ImgR(t1)图像上的候选匹配位置，即得到光斑点的匹配矢量集；逐一搜索每个匹配矢量集合，选择这些匹配矢量集合中最为相似的匹配矢量作为光斑点的匹配矢量，其所对应的匹配位置为图像ImgL(t1)分析窗口内各个斑点在图像ImgR(t1)中的对应位置；利用同样的方法将图像ImgL(t1)中的每个斑点均在图像ImgR(t1)中找到对应位置，并形成位置对应关系的查找表PosLut；
步骤七：在ImgL(t1+Δt)中划分网格，网格的节点为需要进行分析的位置；然后在ImgL(t1+Δt)中逐一选择分析点P，并以分析点为中心以2×wr+1为边长选择一方形分析区域，与此同时在PosLut查找表中查询与P点最为近邻的匹配位置P′，随后在ImgR(t1+Δt)中的P′位置以4×wr+1为边长选择一分析区域；接下来对这两个区域的数据进行归一化协方差相关匹配，并找出最大值的位置P″；按此过程遍历分析ImgL(t1+Δt)上所有的网格节点，计算得到ImgR(t1+Δt)上的对应点，并将位置对应关系记录在查找表StereoLut_1中；
步骤八：获取t2时刻的图像数据，分别记作ImgL(t2)与ImgR(t2)、ImgL(t2+Δt)与ImgR(t2+Δt)；利用步骤F与步骤G所述的方法，得到ImgL(t2+Δt)与ImgR(t2+Δt)之间对应关系的查找表StereoLut_2；
步骤九：仍使用归一化协方差相关匹配方法在分析网格节点处，分别对ImgL(t1+Δt)与ImgL(t2+Δt)、ImgR(t1+Δt)与ImgR(t2+Δt)进行相关分析，此时t1+Δt时刻的两幅图像中使用的分析窗口边长为2×wr+1而t2+Δt时刻的两幅图像中使用的分析窗口边长为2×Wr+1，其中Wr需同时满足Wr≥wr和Wr≥4×Smax，Smax为t1～t2间隔内被测目标物体图像上的最大移动距离；分析结束得到ImgL(t1+Δt)与ImgL(t2+Δt)之间的对应关系，记录在查找表MatchLut_L中，而ImgR(t1+Δt)与ImgR(t2+Δt)之间的对应关系，记录在查找表MatchLut_R中；
步骤十：遍历ImgL(t1+Δt)上所有分析节点，搜索其在ImgL(t2+Δt)以及ImgR(t2+Δt)上的对应位置：令当前一个分析位置记作P0，然后在MatchLut_L中查找P0点在ImgL(t2+Δt)的对应位置P1，并搜索与P1最近邻的网格节点的位置P1′；在StereoLut_2中查找P1′在ImgR(t2+Δt)上的对应位置P1″，之后将P1′与P1″之间的位置变化关系用于P1点，计算得到P1在ImgR(t2+Δt)上的对应位置P1Z；
步骤十一：在StereoLut_1中查找P0点在ImgR(t1+Δt)上的对应位置P2，同时搜索与P2最近邻的网格节点的位置P2′；在MatchLut_R中查找P2′在ImgR(t2+Δt)上的对应位置P2″；同样将P2′与P2″之间的位置变化关系用于P2点，计算得到P2在ImgR(t2+Δt)上的对应位置P2Z；
步骤十二：当||P1Z-P2Z||＜εT时，认为当前分析点P0处可以得到正确的配对结果，需要进一步计算P0点在ImgR(t2+Δt)上的对应位置P3：以P1Z点与P2Z点之间连线的中点为中心以r为半径，要求半径r不小于P1Z点与P2Z点之间连线长度的一半，划定一个分析区域，记作f3；用r/2为半径分别在P1点与P2点处建立分析区域，记为f1与f2；使用f1与f2分别在f3区域遍历性进行归一化协方差相关计算，得到的相关系数矩阵，记作M1与M2；计算两个相关系数矩阵的乘积，记作M＝M1·M2；在M中寻找最大值的位置，该位置即为所求的P3位置；
步骤十三：利用摄像机定标参数对像面上的点进行三维重构：通过[P0,P2]对应点可以得到t1+Δt的空间位置PW1；通过[P1,P3]对应点可以得到t2+Δt的空间位置PW2，因此可以得到一个空间矢量使用相同的分析方法完成对ImgL(t1+Δt)上所有节点的分析后，便可以得到从t1+Δt时刻到t2+Δt时刻的被测泥沙的运动演变信息。
本发明的效果和益处是使用两路脉冲激光器结合摄像机的双曝光功能实现现场图像的跨帧采集，可以在较短的时间间隔内得到带有光斑点阵的泥沙图像与被体积光照明的泥沙图像。使用点阵图像的特征匹配来预测立体匹配，提高了图像立体匹配的运算效率与正确率。分析时采用网格式图像处理，各个流程之间没有数据的依存关系，便于并行作业。此外，采用网格节点的处理结果作为预则值，大幅度提高了分析速度。另外，图像采集过程与实验过程互不干扰，无需停机待时。且一次分析结束便可以得到网格节点处的空间矢量场，无需插值耗时。
附图说明
附图1是本发明的结构示意图。
附图2是系统的时序控制图。
附图3是图像匹配结构图。
图中：1PC机；2同步控制器；3激光器A；4激光器B；5摄像机A；6摄像机B；7脉冲激光器A的导光臂；8脉冲激光器B的导光臂；9点阵光栅片；10负透镜；11被测地形。
具体实施方式
以下结合附图对本发明的实施方式作进一步描述。
一种双脉冲式泥沙运移信息的图像测量方法，其特征是由摄像机、脉冲激光器、和内插图像采集卡的PC机组成；两台左右放置的具有外控双曝光功能的工业摄像机组成立体视觉系统固定在被测地形11可以被观测到的地点；两台脉冲激光器配有各自的脉冲激光导光臂用于照明；其中一台安装点阵光栅片9，另一台安装负透镜10；同步控制器2的设置端口与PC机1相连，用于设置同步控制器2的控制参数；而同步控制器2的输出端子分别与摄像机A5、摄像机B6、激光器A3的氙灯与调Q外控端口、激光器B4的氙灯与调Q的外控端口相连；在同步控制器2的时序控制下，利用点阵光栅片9和负透镜10，两个脉冲激光器可以交替打出点阵光与体积椎光，对泥沙地形进行照明；同时利用摄像机的跨帧采集技术在较短的时间内获取不同角度的点阵图像与泥沙图像；通过第一次曝光获取的点阵图像的配对信息来预测第二次曝光泥沙图像的立体匹配位置；PC机1对不同时刻的图像依照预设定的网格进行归一化协方差相关匹配，从而获取左右像平面上的同名运动矢量，最后结合摄像机定标参数重构出泥沙的空间运动矢量，得到泥沙地形的运移信息；其具体方法如下：
步骤一：设置同步控制器2的时序：首先保证摄像机A5的外触发时序与摄像机B6的外触发时序严格一致；其次控制脉冲激光器A3的氙灯的外触发时序与调Q外触发时序使得脉冲激光器A3的脉冲在摄像机A5第一次曝光间隔内发出，同理控制脉冲激光器B4的氙灯外触发时序与调Q外触发时序使得脉冲激光器B4的脉冲在摄像机B6第二次曝光间隔内发出；
步骤二：左右放置两台摄像机，调整角度使得目标区域可以被观测到；调整脉冲激光器A的导光臂7和脉冲激光器B的导光臂8，使得两台激光器的照明区域均覆盖被测地形11；
步骤三：设置两台工业摄像机为同步外控模式，使用张氏定标法进行两台摄像机的定标，完成每台摄像机的内参数与外参数的计算；
步骤四：设置两台摄像机为双曝光模式；
步骤五：令两台脉冲激光器发出激光脉冲之间的时间间隔为Δt，获取t1时刻的图像，令第一次曝光的左图像为ImgL(t1)右图像为ImgR(t1)，而第二次曝光的左图像为ImgL(t1+Δt)右图像为ImgR(t1+Δt)；
步骤六：分别对ImgL(t1)和ImgR(t1)图像进行二值化处理；然后利用重心法或曲面拟合法对图像中的光斑进行中心提取；之后，在ImgL(t1)图像上建立一个分析窗口，使得该分析窗内包含的光斑数目不小于5个；利用极线约束，寻找分析窗口内每个光斑点在ImgR(t1)图像上的候选匹配位置，即得到光斑点的匹配矢量集；逐一搜索每个匹配矢量集合，选择这些匹配矢量集合中最为相似的匹配矢量作为光斑点的匹配矢量，其所对应的匹配位置为图像ImgL(t1)分析窗口内各个斑点在图像ImgR(t1)中的对应位置；利用同样的方法将图像ImgL(t1)中的每个斑点均在图像ImgR(t1)中找到对应位置，并形成位置对应关系的查找表PosLut；
步骤七：在ImgL(t1+Δt)中划分网格，网格的节点为需要进行分析的位置；然后在ImgL(t1+Δt)中逐一选择分析点P，并以分析点为中心以2×wr+1为边长选择一方形分析区域，与此同时在PosLut查找表中查询与P点最为近邻的匹配位置P′，随后在ImgR(t1+Δt)中的P′位置以4×wr+1为边长选择一分析区域；接下来对这两个区域的数据进行归一化协方差相关匹配，并找出最大值的位置P″；按此过程遍历分析ImgL(t1+Δt)上所有的网格节点，计算得到ImgR(t1+Δt)上的对应点，并将位置对应关系记录在查找表StereoLut_1中；
步骤八：获取t2时刻的图像数据，分别记作ImgL(t2)与ImgR(t2)、ImgL(t2+Δt)与ImgR(t2+Δt)；利用步骤F与步骤G所述的方法，得到ImgL(t2+Δt)与ImgR(t2+Δt)之间对应关系的查找表StereoLut_2；
步骤九：仍使用归一化协方差相关匹配方法在分析网格节点处，分别对ImgL(t1+Δt)与ImgL(t2+Δt)、ImgR(t1+Δt)与ImgR(t2+Δt)进行相关分析，此时t1+Δt时刻的两幅图像中使用的分析窗口边长为2×wr+1而t2+Δt时刻的两幅图像中使用的分析窗口边长为2×Wr+1，其中Wr需同时满足Wr≥wr和Wr≥4×Smax，Smax为t1～t2间隔内被测目标物体图像上的最大移动距离；分析结束得到ImgL(t1+Δt)与ImgL(t2+Δt)之间的对应关系，记录在查找表MatchLut_L中，而ImgR(t1+Δt)与ImgR(t2+Δt)之间的对应关系，记录在查找表MatchLut_R中；
步骤十：遍历ImgL(t1+Δt)上所有分析节点，搜索其在ImgL(t2+Δt)以及ImgR(t2+Δt)上的对应位置：令当前一个分析位置记作P0，然后在MatchLut_L中查找P0点在ImgL(t2+Δt)的对应位置P1，并搜索与P1最近邻的网格节点的位置P1′；在StereoLut_2中查找P1′在ImgR(t2+Δt)上的对应位置P1″，之后将P1′与P1″之间的位置变化关系用于P1点，计算得到P1在ImgR(t2+Δt)上的对应位置P1Z；
步骤十一：在StereoLut_1中查找P0点在ImgR(t1+Δt)上的对应位置P2，同时搜索与P2最近邻的网格节点的位置P2′；在MatchLut_R中查找P2′在ImgR(t2+Δt)上的对应位置P2″；同样将P2′与P2″之间的位置变化关系用于P2点，计算得到P2在ImgR(t2+Δt)上的对应位置P2Z；
步骤十二：当||P1Z-P2Z||＜εT时，认为当前分析点P0处可以得到正确的配对结果，需要进一步计算P0点在ImgR(t2+Δt)上的对应位置P3：以P1Z点与P2Z点之间连线的中点为中心以r为半径，要求半径r不小于P1Z点与P2Z点之间连线长度的一半，划定一个分析区域，记作f3；用r/2为半径分别在P1点与P2点处建立分析区域，记为f1与f2；使用f1与f2分别在f3区域遍历性进行归一化协方差相关计算，得到的相关系数矩阵，记作M1与M2；计算两个相关系数矩阵的乘积，记作M＝M1·M2；在M中寻找最大值的位置，该位置即为所求的P3位置；
步骤十三：利用摄像机定标参数对像面上的点进行三维重构：通过[P0,P2]对应点可以得到t1+Δt的空间位置PW1；通过[P1,P3]对应点可以得到t2+Δt的空间位置PW2，因此可以得到一个空间矢量使用相同的分析方法完成对ImgL(t1+Δt)上所有节点的分析后，便可以得到从t1+Δt时刻到t2+Δt时刻的被测泥沙的运动演变信息。
Claims (1)
1.一种双脉冲式泥沙运移信息的图像测量方法，其特征是由摄像机、脉冲激光器和内插图像采集卡的PC机组成；两台左右放置的具有外控双曝光功能的工业摄像机组成立体视觉系统固定在被测地形(11)可以被观测到的地点；两台脉冲激光器配有各自的脉冲激光导光臂用于照明；其中一台安装点阵光栅片(9)，另一台安装负透镜(10)；同步控制器(2)的设置端口与PC机(1)相连，用于设置同步控制器(2)的控制参数；而同步控制器(2)的输出端子分别与摄像机A(5)、摄像机B(6)、脉冲激光器A(3)的氙灯与调Q外控端口、脉冲激光器B(4)的氙灯与调Q的外控端口相连；在同步控制器(2)的时序控制下，利用点阵光栅片(9)和负透镜(10)，两个脉冲激光器可以交替打出点阵光与体积椎光，对泥沙地形进行照明；同时利用摄像机的跨帧采集技术在较短的时间内获取不同角度的点阵图像与泥沙图像；通过第一次曝光获取的点阵图像的配对信息来预测第二次曝光泥沙图像的立体匹配位置；PC机(1)对不同时刻的图像依照预设定的网格进行归一化协方差相关匹配，从而获取左右像平面上的同名运动矢量，最后结合摄像机定标参数重构出泥沙的空间运动矢量，得到泥沙地形的运移信息；其具体方法如下：
步骤一：设置同步控制器(2)的时序：首先保证摄像机A(5)的外触发时序与摄像机B(6)的外触发时序严格一致；其次控制脉冲激光器A(3)的氙灯的外触发时序与调Q外触发时序使得脉冲激光器A(3)的脉冲在摄像机A(5)第一次曝光间隔内发出，同理控制脉冲激光器B(4)的氙灯外触发时序与调Q外触发时序使得脉冲激光器B(4)的脉冲在摄像机B(6)第二次曝光间隔内发出；
步骤二：左右放置两台摄像机，调整角度使得目标区域可以被观测到；调整脉冲激光器A的导光臂(7)和脉冲激光器B的导光臂(8)，使得两台脉冲激光器的照明区域均覆盖被测地形(11)；
步骤三：设置两台工业摄像机为同步外控模式，使用张氏定标法进行两台摄像机的定标，完成每台摄像机的内参数与外参数的计算；
步骤四：设置两台摄像机为双曝光模式；
步骤五：令两台脉冲激光器发出激光脉冲之间的时间间隔为Δt，获取t1时刻的图像，令第一次曝光的左图像为ImgL(t1)右图像为ImgR(t1)，而第二次曝光的左图像为ImgL(t1+Δt)右图像为ImgR(t1+Δt)；
步骤六：分别对ImgL(t1)和ImgR(t1)图像进行二值化处理；然后利用重心法或曲面拟合法对图像中的光斑进行中心提取；之后，在ImgL(t1)图像上建立一个分析窗口，使得该分析窗内包含的光斑数目不小于5个；利用极线约束，寻找分析窗口内每个光斑点在ImgR(t1)图像上的候选匹配位置，即得到光斑点的匹配矢量集；逐一搜索每个匹配矢量集合，选择这些匹配矢量集合中最为相似的匹配矢量作为光斑点的匹配矢量，其所对应的匹配位置为图像ImgL(t1)分析窗口内各个斑点在图像ImgR(t1)中的对应位置；利用同样的方法将图像ImgL(t1)中的每个斑点均在图像ImgR(t1)中找到对应位置，并形成位置对应关系的查找表PosLut；
步骤七：在ImgL(t1+Δt)中划分网格，网格的节点为需要进行分析的位置；然后在ImgL(t1+Δt)中逐一选择分析点P，并以分析点为中心以2×wr+1为边长选择一方形分析区域，与此同时在PosLut查找表中查询与P点最为近邻的匹配位置P′，随后在ImgR(t1+Δt)中的P′位置以4×wr+1为边长选择一分析区域；接下来对这两个区域的数据进行归一化协方差相关匹配，并找出最大值的位置P″；按此过程遍历分析ImgL(t1+Δt)上所有的网格节点，计算得到ImgR(t1+Δt)上的对应点，并将位置对应关系记录在查找表StereoLut_1中；
步骤八：获取t2时刻的图像数据，分别记作ImgL(t2)与ImgR(t2)、ImgL(t2+Δt)与ImgR(t2+Δt)；利用步骤六与步骤七所述的方法，得到ImgL(t2+Δt)与ImgR(t2+Δt)之间对应关系的查找表StereoLut_2；
步骤九：仍使用归一化协方差相关匹配方法在分析网格节点处，分别对ImgL(t1+Δt)与ImgL(t2+Δt)、ImgR(t1+Δt)与ImgR(t2+Δt)进行相关分析，此时t1+Δt时刻的两幅图像中使用的分析窗口边长为2×wr+1而t2+Δt时刻的两幅图像中使用的分析窗口边长为2×Wr+1，其中Wr需同时满足Wr≥wr和Wr≥4×Smax，Smax为t1～t2间隔内被测目标物体图像上的最大移动距离；分析结束得到ImgL(t1+Δt)与ImgL(t2+Δt)之间的对应关系，记录在查找表MatchLut_L中，而ImgR(t1+Δt)与ImgR(t2+Δt)之间的对应关系，记录在查找表MatchLut_R中；
步骤十：遍历ImgL(t1+Δt)上所有分析节点，搜索其在ImgL(t2+Δt)以及ImgR(t2+Δt)上的对应位置：令当前一个分析位置记作P0，然后在MatchLut_L中查找P0点在ImgL(t2+Δt)的对应位置P1，并搜索与P1最近邻的网格节点的位置P′1；在StereoLut_2中查找P′1在ImgR(t2+Δt)上的对应位置P″1，之后将P′1与P″1之间的位置变化关系用于P1点，计算得到P1在ImgR(t2+Δt)上的对应位置P1Z；
步骤十一：在StereoLut_1中查找P0点在ImgR(t1+Δt)上的对应位置P2，同时搜索与P2最近邻的网格节点的位置P′2；在MatchLut_R中查找P′2在ImgR(t2+Δt)上的对应位置P″2；同样将P′2与P″2之间的位置变化关系用于P2点，计算得到P2在ImgR(t2+Δt)上的对应位置P2Z；
步骤十二：当||P1Z-P2Z||＜εT时，认为当前分析点P0处可以得到正确的配对结果，需要进一步计算P0点在ImgR(t2+Δt)上的对应位置P3：以P1Z点与P2Z点之间连线的中点为中心以r为半径，要求半径r不小于P1Z点与P2Z点之间连线长度的一半，划定一个分析区域，记作f3；用r/2为半径分别在P1点与P2点处建立分析区域，记为f1与f2；使用f1与f2分别在f3区域遍历性进行归一化协方差相关计算，得到的相关系数矩阵，记作M1与M2；计算两个相关系数矩阵的乘积，记作M＝M1·M2；在M中寻找最大值的位置，该位置即为所求的P3位置；
步骤十三：利用摄像机定标参数对像面上的点进行三维重构：通过[P0,P2]对应点可以得到t1+Δt的空间位置PW1；通过[P1,P3]对应点可以得到t2+Δt的空间位置PW2，因此可以得到一个空间矢量使用相同的分析方法完成对ImgL(t1+Δt)上所有节点的分析后，便可以得到从t1+Δt时刻到t2+Δt时刻的被测泥沙的运动演变信息。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410182506.9A CN103969024B (zh) | 2014-04-30 | 2014-04-30 | 一种双脉冲式泥沙运移信息的图像测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410182506.9A CN103969024B (zh) | 2014-04-30 | 2014-04-30 | 一种双脉冲式泥沙运移信息的图像测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103969024A CN103969024A (zh) | 2014-08-06 |
| CN103969024B true CN103969024B (zh) | 2016-05-11 |
Family
ID=51238796
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410182506.9A Expired - Fee Related CN103969024B (zh) | 2014-04-30 | 2014-04-30 | 一种双脉冲式泥沙运移信息的图像测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103969024B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105806319B (zh) * | 2016-03-11 | 2018-04-10 | 大连理工大学 | 一种用于玻璃管线三维运动分析的跨轴式图像测量方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1563925A (zh) * | 2004-03-23 | 2005-01-12 | 天津大学 | 高速复杂流场的全场光学测量方法与测量装置 |
| CN1664253A (zh) * | 2005-04-21 | 2005-09-07 | 中国科学院力学研究所 | 在海流作用下海底输沙状况的模拟装置 |
| EP1777553A1 (fr) * | 2005-10-19 | 2007-04-25 | Commissariat à l'Energie Atomique | Procédé de mesure de la dynamique sédimentaire d'un site, utilisant un traceur artificiel activable |
| CN2906578Y (zh) * | 2006-01-17 | 2007-05-30 | 常兆丰 | 风沙流输沙量测定仪 |
| EP2339327A1 (en) * | 2009-12-23 | 2011-06-29 | F. Hoffmann-La Roche AG | Automatic test for imaging units |
| EP2594969A1 (fr) * | 2011-11-17 | 2013-05-22 | IFP Energies Nouvelles | Procédé d'exploration d'un bassin sédimentaire à partir d'une simulation stratigraphique |
-
2014
- 2014-04-30 CN CN201410182506.9A patent/CN103969024B/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1563925A (zh) * | 2004-03-23 | 2005-01-12 | 天津大学 | 高速复杂流场的全场光学测量方法与测量装置 |
| CN1664253A (zh) * | 2005-04-21 | 2005-09-07 | 中国科学院力学研究所 | 在海流作用下海底输沙状况的模拟装置 |
| EP1777553A1 (fr) * | 2005-10-19 | 2007-04-25 | Commissariat à l'Energie Atomique | Procédé de mesure de la dynamique sédimentaire d'un site, utilisant un traceur artificiel activable |
| CN2906578Y (zh) * | 2006-01-17 | 2007-05-30 | 常兆丰 | 风沙流输沙量测定仪 |
| EP2339327A1 (en) * | 2009-12-23 | 2011-06-29 | F. Hoffmann-La Roche AG | Automatic test for imaging units |
| EP2594969A1 (fr) * | 2011-11-17 | 2013-05-22 | IFP Energies Nouvelles | Procédé d'exploration d'un bassin sédimentaire à partir d'une simulation stratigraphique |
Non-Patent Citations (2)
| Title |
|---|
| 刘春嵘等.复杂流动下泥沙起动概率的图像测量.《湖南大学学报(自然科学版)》.2008,(第03期), * |
| 李丹勋等.图像处理技术及其在泥沙运动研究中的应用.《水利水电技术》.1999,(第05期), * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103969024A (zh) | 2014-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11334960B2 (en) | Systems and methods for pipelined processing of sensor data using hardware | |
| AU2017324923B2 (en) | Predicting depth from image data using a statistical model | |
| CN109791696B (zh) | 用于场景的3d重建的方法、设备和定位事件摄像机的方法 | |
| US9681150B2 (en) | Optical flow determination using pyramidal block matching | |
| US20220051425A1 (en) | Scale-aware monocular localization and mapping | |
| CN111259706B (zh) | 一种车辆的车道线压线判断方法和系统 | |
| ITRM970530A1 (it) | Procedimento per il riconoscimento stereoscopico di oggetti | |
| CN110060299A (zh) | 基于双目视觉技术的输电通道内危险源识别和定位方法 | |
| CN104517095B (zh) | 一种基于深度图像的人头分割方法 | |
| CN103971405A (zh) | 一种激光散斑结构光及深度信息的三维重建方法 | |
| JP2022511990A (ja) | 情報補完方法、車線認識方法、インテリジェント走行方法及び関連製品 | |
| CN111259704B (zh) | 一种虚线车道线端点检测模型的训练方法 | |
| CN102253057B (zh) | 内窥镜系统及使用内窥镜系统的测量方法 | |
| CN105222725B (zh) | 一种基于光谱分析的高清图像动态采集方法 | |
| CN111026115A (zh) | 一种基于深度学习的机器人避障控制方法及装置 | |
| CN103969024B (zh) | 一种双脉冲式泥沙运移信息的图像测量方法 | |
| JP2024052755A (ja) | 三次元変位計測方法及び三次元変位計測装置 | |
| Elzayat et al. | Real-time car detection-based depth estimation using mono camera | |
| WO2022179124A1 (zh) | 一种图像修复方法及装置 | |
| US20210134053A1 (en) | Three-dimensional point data based on stereo reconstruction using structured light | |
| CN110111332A (zh) | 基于深度卷积神经网络的胶原蛋白肠衣缺陷检测模型、检测方法及系统 | |
| CN113221744B (zh) | 一种基于深度学习的单目图像3d物体检测方法 | |
| JP5642049B2 (ja) | 車両用外界認識装置およびそれを用いた車両システム | |
| IT201900001833A1 (it) | Un procedimento per rilevare oggetti in movimento mediante una telecamera in movimento, e relativo sistema di elaborazione, dispositivo e prodotto informatico | |
| JP2013205867A (ja) | 情報処理装置及び方法、並びにプログラム |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160511 Termination date: 20200430 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |