CN102538763B - 一种河工模型试验三维地形的量测方法 - Google Patents
一种河工模型试验三维地形的量测方法 Download PDFInfo
- Publication number
- CN102538763B CN102538763B CN201210033188.0A CN201210033188A CN102538763B CN 102538763 B CN102538763 B CN 102538763B CN 201210033188 A CN201210033188 A CN 201210033188A CN 102538763 B CN102538763 B CN 102538763B
- Authority
- CN
- China
- Prior art keywords
- camera
- bed surface
- point
- spike
- laser
- 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 30
- 238000012360 testing method Methods 0.000 title claims abstract description 17
- 238000003384 imaging method Methods 0.000 claims description 19
- 238000012876 topography Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 239000000571 coke Substances 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 230000003287 optical effect Effects 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 abstract description 11
- 238000005516 engineering process Methods 0.000 abstract description 10
- 230000004438 eyesight Effects 0.000 abstract description 4
- 239000011159 matrix material Substances 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000004304 visual acuity Effects 0.000 description 2
- 235000012773 waffles Nutrition 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000005137 deposition process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000005329 float glass Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000000608 laser ablation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- 230000009184 walking Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/40—Protecting water resources
- Y02A20/402—River restoration
-
- 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
- Length Measuring Devices By Optical Means (AREA)
- Image Processing (AREA)
Abstract
本发明属于水利量测技术领域,特别涉及三维地形的测量方法。在水利工程的河工模型试验中实时监测河道的变形是河工模型试验的关键技术。基于双目立体视觉原理,利用图像摄取技术、激光矩阵光源示踪技术和三维地形重构技术实时量测河工模型试验中床面的三维形态。集成数码摄像机、激光矩形点阵示踪光斑发射装置和计算机为一体,通过设定的系统安装、相机参数校正、床面上激光点阵光斑投影,经济、安全、高质量地实现河工模型床面示踪、床面图像采集、河床激光示踪点提取及匹配,结合计算公式计算出床面上激光示踪光斑的三维空间坐标。基于计算机,以床面上所有激光示踪光斑的三维空间坐标为媒介生成河工模型的地形图。
Description
技术领域
本发明属于水利量测技术领域,特别是河工模型试验的复杂三维地形的动态量测。
背景技术
在水利工程的河工模型试验中,河床地形将随着河道的冲淤过程而变化,如河岸坍塌、边滩消长等。实时监测河道的变形是河工模型试验的关键技术。
现有的河工模型地形测量技术多为单点测量,常用的方法一般有①钢尺测量法、②测针测量法、③光电式地形仪、④电阻式地形仪、⑤超声波地形仪等,并且前4个测量方法为接触式,操作过程中会对床面造成破坏。另外,这些传统的测量方法费工费时,以常用的④或⑤为例,首先要在模型河道上布置测量横断面,在两岸安装轨道并在轨道上架设测桥、安装测车,采用安装在测车上的超声地形仪、电阻测淤仪等在测桥上自动行走,测量断面上各点的河床高度,继而通过相关软件重构河道的三维地形。北京尚水信息技术股份有限公司的产品“超声地形自动测量分析系统”具有代表性。但是此类单点测量方法工作量过大、数据量过少、测量效率过低,三维重构的精度不高。
近年来,作为非接触式量测技术的激光扫描仪在河工模型地形量测领域内得到了一定的发展,如瑞士的徕卡HDS3000三维激光扫描仪和奥地利的Riegl LPM-2K三维激光扫描仪,可自动进行三维地形的单点连续扫描,并将扫描的散点数据重构生成三维地形。其性能指标:扫描距离2-150米;扫描精度2-6毫米;扫描角度270×360度;扫描速度1000点/秒。该仪器测量精度高,扫描的速度快,不需要轨道、测桥、测车等辅助设施,很适合于固定构造物的一次性精确扫描,如城市勘察、文物建筑测绘等。但在用于河工模型试验中的动态河床测量时,每次移动都需要重新标定仪器的三维坐标和水平方向角;由于一次地形的测量工作量很大,难以实现实时测量;当仪器视线与被测地形的夹角较小时,测量精度较差,甚至会因为地形起伏的遮挡而产生错误;且该仪器的价格昂贵。
近几十年来,双目立体视觉技术在工业、农业、军事、考古、交通等相关领域得到了极大的应用与发展,如机器人快速移动避障碍、植物叶片三维立体重建、月球、火星表面地形探测等。双目立体视觉是通过对同一目标的两幅图像提取、识别、匹配,重建三维环境信息的过程。主要包括:图像捕获、摄像机标定、图像预处理和特征提取、立体匹配以及所测物体表面形状三维重建。该技术在河工模型地形量测方面也得到了初步的尝试,如河海大学毛野等人运用该技术研究河道采砂对长江镇江段河床演变的影响等,但是类似的尝试在河工模型床面上特征点的示踪、匹配、提取及三维空间坐标计算等问题上并没有得到很好的解决。
发明内容
本发明提出了一种河工模型试验的三维地形动态测量方法,用以实时量测河工模型试验中床面的三维形态。本方法基于双目立体视觉原理,利用图像摄取技术、激光矩阵光源示踪技术和三维地形重构技术来实时量测河工模型试验中床面的三维形态。具有经济实用、安全可靠、质量高的优点。
本发明提出的一种河工模型试验三维地形的测量方法,其特征在于,采用两台参数相同的数码相机,一台连续激光矩形点阵示踪光斑发射装置,以及分别与两台数码相机相连的一台计算机;所述连续激光矩形点阵示踪光斑发射装置由连续激光器和微调光栅所组成;该微调光栅在河工模型床面上投射出不同间距及大小的矩形点阵示踪光斑,其中矩形点阵示踪光斑中心点为定位点O;
该方法,包括以下步骤:
1)首先进行相机焦距f、两相机两焦点距离L以及相机偏角φ标定:
1-1)相机焦距f标定计算公式为:
式中,a为CCD较短一边的边长,y2指相机初始拍摄垂向最大高度,y1为相机沿其光轴轴向向后移动ΔX后所拍摄的垂向最大高度(y1>y2);
1-2)两焦点距离L标定为:
式中,L1为OB1的长度,O为定位点,B1为与点B在左相机成像面上具有相同像素坐标的床面上的一点,点B为在点O垂直向上y高度处的一点;角度α1为点O与左相机焦点F1的连线F1O与直线F1B1的夹角,角度β1为点O与点B1的连线B1O与直线B1F1的夹角,其中α1和β1在标定时计算出α1=arctan(O1P1/f),β1=arctan(y/L1),O1P1为床面上线段OB1在左相机像平面上的图像长度;
1-3)相机偏角φ标定为:
φ=180-η=α1+β1
式中,α1为直线F1O与直线F1B1之间的夹角;β1为直线B1O与直线B1F1之间的夹角;
2)照相机标定完成后,在两相机焦点中间位置垂直向下安装连续激光矩形点阵示踪光斑发射装置,激光示踪光斑中心点为定位点O,定位点O在左右两相机像平面的成像分别为两相机像平面的中心点O1和O2;左右两相机成像面以中心点O1和O2为原点分别建立各自的直角坐标系x1O1y1和x2O2y2,床面上也以点O为中心点建立坐标系XOY,将两相机像平面及床面均分成四个象限,且每个象限所包含或所成像的激光示踪斑点一一对应;
3)对床面进行拍照,分别得到左右两相机所拍摄的系列照片;对所述系列照片进行噪声处理、畸变处理、像素二值化及对示踪光斑边缘检测的常规预处理,以及进行示踪光斑中心坐标的计算,再对示踪斑点进行匹配及特征提取,获得床面上同一个示踪光斑点P分别在左右两相机像平面上的像素坐标P1(m1p,n1p)和P2(m2p,n2p),再根据相机像素与CCD物理尺寸的转换关系,将示踪光斑点P在左右相机像平面上的像素坐标分别转化为像平面上的图像坐标P1(x1p,y1p)和P2(x2p,y2p);
4)以左相机焦点F1(0,0,0)为原点建立三维空间坐标系O-XYZ:F1F2为X轴、从点F1向下引垂线为Z轴,过点F1沿河床纵向方向为Y轴,根据数码相机的标定参数f、L、φ,将床面上每个示踪光斑点P在左右两相机像平面上的图像坐标P1(x1p,y1p)和P2(x2p,y2p)按转换公式换算为空间坐标系O-XYZ中的三维空间坐标P(XP,YP,ZP);
激光示踪光斑图像坐标与三维空间坐标的转换公式为:
式中,P1(x1p,y1p)和p2(x2p,y2p)分别为模型床面上激光示踪斑点P在左右两个像平面上的图像坐标;f为左右相机的焦距;L为左右相机的焦点间距;φ为相机偏角;θ1和θ2为虚拟角度,表示模型床面上的点P在平面I上的投影点PT与左右两相机焦点F1与点F2的连线F1PT和F2PT分别与线段F1F2和F2F1的夹角,平面I为过直线F1F2且垂直于模型床面的平面;
本发明的技术特点及有益效果:
1、针对河工模型的大小,由于采用定焦相机,通过调节左右两个相机像平面中垂线与水平方向的夹角或改变两相机与模型床面垂向距离的大小来调节数码相机的拍摄范围。
2、在左右两相机中间处设置连续激光矩形点阵示踪光斑发射装置,该装置将激光示踪光斑呈矩形状排列均匀投射到河床表面。设定激光示踪光斑中心点为定位点O,通过左右两个相机对床面上的激光示踪点进行同步拍照,利用经典的特征点匹配及提取技术(如Zhang提出的核线约束等立体匹配方法),可以得到同一个激光示踪光斑在左右两张照片中的图像坐标。根据图像坐标与空间坐标的换算公式,计算出相应激光示踪点的三维空间坐标。
3、基于计算机,以床面上所有激光示踪光斑的三维空间坐标为媒介生成河工模型地形图。
4、本方法具有经济实用、安全可靠、质量高的优点。
附图说明
图1为本发明的测量系统布置示意图;
图2为本发明的测量系统左右两相机焦点距离L和偏角φ标定示意图。
具体实施方式
本发明提出的一种河工模型试验三维地形的测量方法,其特征在于,采用两台参数相同的数码相机,一台连续激光矩形点阵示踪光斑发射装置,以及分别与两台数码相机相连的一台计算机;所述连续激光矩形点阵示踪光斑发射装置由连续激光器和微调光栅所组成;该微调光栅在河工模型床面上投射出不同间距及大小的矩形点阵示踪光斑,其中矩形点阵示踪光斑中心点为定位点O;
该方法,包括以下步骤:
1)首先进行相机焦距f、两相机两焦点距离L以及相机偏角φ标定:
1-1)相机焦距f标定计算公式为:
式中,a为CCD较短一边的边长,y2指相机初始拍摄垂向最大高度,y1为相机沿其光轴轴向向后移动ΔX后所拍摄的垂向最大高度(y1>y2);
1-2)两焦点距离L标定为:
式中,L1为OB1的长度,O为定位点,B1为与点B在左相机成像面上具有相同像素坐标的床面上的一点,点B为在点O垂直向上y高度处的一点;角度α1为点O与左相机焦点F1的连线F1O与直线F1B1的夹角,角度β1为点O与点B1的连线B1O与直线B1F1的夹角,其中α1和β1在标定时计算出α1=arctan(O1P1/f),β1=arctan(y/L1),O1P1为床面上线段OB1在左相机像平面上的图像长度;
1-3)相机偏角φ标定为:
φ=180-η=α1+β1
式中,α1为直线F1O与直线F1B1之间的夹角;β1为直线B1O与直线B1F1之间的夹角;
2)照相机标定完成后,在两相机焦点中间位置垂直向下安装连续激光矩形点阵示踪光斑发射装置,激光示踪光斑中心点为定位点O,定位点O在左右两相机像平面的成像分别为两相机像平面的中心点O1和O2;左右两相机成像面以中心点O1和O2为原点分别建立各自的直角坐标系x1O1y1和x2O2y2,床面上也以点O为中心点建立坐标系XOY,将两相机像平面及床面均分成四个象限,且每个象限所包含或所成像的激光示踪斑点一一对应;
3)对床面进行拍照,分别得到左右两相机所拍摄的系列照片;对所述系列照片进行噪声处理、畸变处理、像素二值化及对示踪光斑边缘检测的常规预处理,以及进行示踪光斑中心坐标的计算,再对示踪斑点进行匹配及特征提取,获得床面上同一个示踪光斑点P分别在左右两相机像平面上的像素坐标P1(m1p,n1p)和P2(m2p,n2p),再根据相机像素与CCD物理尺寸的转换关系,将示踪光斑点P在左右相机像平面上的像素坐标分别转化为像平面上的图像坐标P1(x1p,y1p)和P2(x2p,y2p);
4)以左相机焦点F1(0,0,0)为原点建立三维空间坐标系O-XYZ:F1F2为X轴、从点F1向下引垂线为Z轴,过点F1沿河床纵向方向为Y轴,根据数码相机的标定参数f、L、φ,将床面上每个示踪光斑点P在左右两相机像平面上的图像坐标P1(x1p,y1p)和P2(x2p,y2p)按转换公式换算为空间坐标系O-XYZ中的三维空间坐标P(XP,YP,ZP);
激光示踪光斑图像坐标与三维空间坐标的转换公式为:
式中,P1(x1p,y1p)和P2(x2p,y2p)分别为模型床面上激光示踪斑点P在左右两个像平面上的图像坐标;f为左右相机的焦距;L为左右相机的焦点间距;φ为相机偏角;θ1和θ2为虚拟角度,表示模型床面上的点P在平面I上的投影点PT与左右两相机焦点F1与点F2的连线F1PT和F2PT分别与线段F1F2和F2F1的夹角,平面I为过直线F1F2且垂直于模型床面的平面;
5),以床面上所有激光示踪光斑的三维空间坐标生成河工模型地形图(可通过计算机内设的绘图程序)。
以下将结合附图及实施例对本发明的一种河工模型试验三维地形的量测方法进行详细说明。
本发明提出的河工模型试验三维地形的测量方法,采用两台参数相同的数码相机1、2,一台连续激光矩形点阵示踪光斑发射装置3和计算机,如图1所示;所述数码相机为高分辨率定焦的两台同参数佳能照相机、所述连续激光矩形点阵示踪光斑发射装置由2W的连续激光器和微调光栅所组成;所述计算机为计算河工模型床面上激光示踪点的三维空间坐标以及实现河工模型地形可视化的平台。
上述装置的各部件的具体实施例说明如下:
两台高分辨率定焦相机参数为:
相机机身采用佳能(canon)EOS 500D:传感器类型:CMOS;传感器尺寸:22.3×14.9mm;有效像素:1510万(最高分辨率4752×3168);影像处理器:DIGIC 4;快门速度:1/60-1/4000秒,采用全自动模式;镜头采用135mm全画幅广角镜头;最大光圈F2.8;最小光圈F22;最近对焦距离0.25m;焦距范围20mm;放大倍率0.14倍;视角范围94度。
连续激光器和微调光栅参数如下:
激光器型号为LW-GL-532nm,功率为2W,工作方式为CW,TTL,Analogue,功率稳定性<5%,发散角<1,束腰直径<3mm。
微调光栅为常用的正交衍射光栅,由通过母板用激光刻蚀到超白浮法玻璃上投射出矩形点阵光斑的两个光学镜片构成,每个镜片厚度1.5mm,尺寸10×10mm,投射点阵100×100(投射光斑数量)。当连续激光经过两光学镜片时,微调两镜片之间的距离可以在河工模型床面上投射出不同间距及大小的矩形点阵示踪光斑,其中矩形点阵示踪光斑中心点为定位点O,其亮度最高,容易识别。
计算机配置参数如下:
CPU性能不低于Intel奔腾4处理器、内存不低于512M、主板配备标准串口卡、操作系统版本不低于Windows XP。
本发明方法,包括以下步骤:
1)首先进行相机焦距f、两相机两焦点距离L以及相机偏角φ的标定:如图1所示,在数码相机安装时,根据河工模型的大小,调整左右两数码相机之间的距离、相机与床面的垂直高度H及其偏角φ(左右两数码相机像平面N1、N2的中垂线与水平线F1F2的夹角),最终实现两台数码相机的拍照范围基本重叠,可采用多套系统覆盖整个河工模型;
其中:
1-1)相机焦距f标定方法:
将数码相机固定在大型车床移动平台上,在数码相机正前方垂直安装一个精密网格板;移动车床的平台,拍摄网格板的图像,记录移动的距离Δx;根据几何条件解算出焦距f,焦距f的计算公式为:
式中,a为CCD较短一边的边长,y2指相机初始拍摄垂向最大高度,y1为相机沿其光轴轴向向后移动ΔX后所拍摄的垂向最大高度(y1>y2);
1-2)两焦点距离L标定:
具体方法如图2所示,对于选定的数码相机和定焦镜头,确立焦距f和CCD的物理尺寸与像素的关系(总像素除以CCD尺寸得到每一像素的实际的物理长度);即在现场安装时,两个相机对称相向安装,两摄像机的像平面中垂线F1O、F2O共面且该面垂直于地面,其中定位点O在两相机像平面的成像均为两相机像平面的中心点;
过点O作一条直线l,该直线在两相机像平面的成像分别为两相机的像平面N1和N2的中心线l1和l2,设两相机焦点F1和F2在直线l上的垂直投影为点T1和T2。在数码相机景深允许的范围内,从点O垂直向上y高度处设置一点B,记录点B在左相机像平面上的像素坐标(X,Y),根据该像素坐标确定出在直线l上的点B1(点B1和点B在左相机像平面的像素坐标相同,均为点P1的像素坐标)。测量O、B1的距离为L1,计算角度β1(∠F1B1O):
点O和点B1左相机像平面的成像分别为O1和P1,根据CCD尺寸及像素坐标的换算公式计算出O1P1的长度,进而计算出α1(∠P1F1O1)的大小(arctan(O1O1/f))。已知α1、β1和L1,则:
η=180-α1-β1 (3)
其中,η为∠F1OB1,L2为OF1的长度。由此可以计算出OT1的长度L3:
L3=L2cos(180-η) (5)
由于相机对称安装,可得两相机焦点之间的距离L为:
1-3)相机偏角φ的标定:φ=α1+β1
由于直线F1F2与直线B1B2平行,可得相机偏角:
φ=180-η=α1+β1 (7)
2)照相机标定完成后,在两相机焦点中间位置垂直向下安装连续激光矩形点阵示踪光斑发射装置,根据拍照范围微调光栅,在床面上形成大小适宜,分布均匀、清晰识别的矩形点阵激光示踪光斑,激光示踪光斑中心点为定位点O,定位点O在左右两相机像平面的成像分别为两相机像平面的中心点O1和O2;左右两相机成像面以中心点O1和O2为原点分别建立各自的直角坐标系x1O1y1和x2O2y2,床面上也以点O为中心点建立坐标系XOY,将两相机像平面及床面均分成四个象限,且每个象限所包含或所成像的激光示踪斑点一一对应;
3)上述所有设备稳定布置后,对床面进行拍照,分别得到左右两相机所拍摄的系列照片;对所述系列照片进行噪声处理、畸变处理、像素二值化及对示踪光斑边缘检测的常规预处理,以及进行示踪光斑中心坐标的计算(可采用如质心法、最小二乘法、椭圆中心修正法),再对激光示踪斑点进行匹配及特征提取(采用经典的特征点匹配及提取技术),获得床面上同一个示踪光斑(如点P)在左右两相机像平面上的像素坐标P1(m1p,n1p)和P2(m2p,n2p),再根据相机像素与CCD物理尺寸的转换关系,将示踪光斑在左右相机像平面上的像素坐标转化为对应的照片的图像坐标P1(x1p,y1p)和P2(x2p,y2p);
4)以左相机焦点F1(0,0,0)为原点建立三维空间坐标系O-XYZ:F1F2为X轴、从点F1向下引垂线为Z轴,过点F1沿河床纵向方向为Y轴,如附图1所示;根据数码相机的标定参数f、L、φ,将床面上每个示踪光斑点P在左右两相机像平面上的图像坐标P1(x1p,y1p)和P2(x2p,y2p)按转换公式(8)换算为空间坐标系O-XYZ中的三维空间坐标P(XP,YP,ZP);
激光示踪光斑图像坐标与三维空间坐标的转换公式为:
式中,P1(x1p,y1p)和P2(x2p,y2p)分别为模型床面上激光示踪斑点P在左右两个像平面上的图像坐标;f为左右相机的焦距;L为左右相机的焦点间距;φ为相机偏角;θ1和θ2为虚拟角度,表示模型床面上的点P在平面I上的投影点PT与左右两相机焦点F1与点F2的连线F1PT和F2PT分别与线段F1F2和F2F1的夹角,平面I为过直线F1F2且垂直于模型床面的平面;
以第四象限为例推求激光示踪斑点图像坐标与空间坐标的转换公式(8)如下:
根据两个相数码机和连连续激光矩形点阵示踪光斑发射装置的安装,左右两相机焦点F1(0,0,0)和F2(L,0,0)与激光示踪斑定位点O三点所在的平面与模型床面垂直,记F1、F2、O三点确定的平面为平面I。选择河工模型床面上第四象限某一个激光示踪光斑点P,点P在左右两相机像平面中的成像分别为点P1和点P2,点P1和点P2投影到I平面分别为点K1(在线l1上)和点K2(在线l2上),床面上的点P在平面I上的投影为点PT,P1相对于左相机像平面中心点O1的偏角∠O1F1K1记为ψ1,同样记∠P1F1K1为记∠K1F1F2为θ1。同样地,仿照上述分析,可以定义ψ2、θ2。由安装时标定的φ角和计算出的角ψ1和角Ψ2,可以求出θ1=φ-ψ1(∠K1F1F2)以及θ2=φ+ψ2(∠K2F2F1)。根据空间解析几何分析(利用三角形正弦定理、三角函数公式)可得:
由此计算出河工模型床面上激光示踪光斑点P在以F1(0,0,0)为原点空间坐标系的三维空间坐标;
在实际应用中,如果测量的河工模型范围较大,可以布设多套本发明的系统,各系统之间有一定的重复范围,通过各系统之间坐标的旋转与平移来对各系统的数据进行集成、最后重构河工模型完整的三维地形图。
Claims (1)
1.一种河工模型试验三维地形的测量方法,其特征在于,采用两台参数相同的数码相机,一台连续激光矩形点阵示踪光斑发射装置,以及分别与两台数码相机相连的一台计算机;所述连续激光矩形点阵示踪光斑发射装置由连续激光器和微调光栅所组成;该微调光栅在河工模型床面上投射出不同间距及大小的矩形点阵示踪光斑,其中矩形点阵示踪光斑中心点为定位点O;
该方法,包括以下步骤:
1)首先进行相机焦距f、两相机两焦点距离L以及相机偏角φ标定:
1-1)相机焦距f标定计算公式为:
式中,a为CCD较短一边的边长,y2指相机初始拍摄垂向最大高度,y1为相机沿其光轴轴向向后移动ΔX后所拍摄的垂向最大高度,y1>y2;
式中,L1为OB1的长度,O为定位点,B1为与点B在左相机成像面上具有相同像素坐标的床面上的一点,点B为在点O垂直向上y高度处的一点;角度α1为点O与左相机焦点F1的连线F1O与直线F1B1的夹角,角度β1为点O与点B1的连线B1O与直线B1F1的夹角,其中α1和β1在标定时计算出α1=arctan(O1P1/f),β1=arctan(y/L1),O1P1为床面上线段OB1在左相机像平面上的图像长度;
1-3)相机偏角φ标定为:
φ=α1+β1
式中,α1为直线F1O与直线F1B1之间的夹角;β1为直线B1O与直线B1F1之间的夹角;
2)照相机标定完成后,在两相机焦点中间位置垂直向下安装连续激光矩形点阵示踪光斑发射装置,激光示踪光斑中心点为定位点O,定位点O在左右两相机像平面的成像分别为两相机像平面的中心点O1和O2;左右两相机成像面以中心点O1和O2为原点分别建立各自的直角坐标系x1O1y1和x2O2y2,床面上也以点O为中心点建立坐标系XOY,将两相机像平面及床面均分成四个象限,且每个象限所包含或所成像的激光示踪斑点一一对应;
3)对床面进行拍照,分别得到左右两相机所拍摄的系列照片;对所述系列照片进行噪声处理、畸变处理、像素二值化及对示踪光斑边缘检测的常规预处理,以及进行示踪光斑中心坐标的计算,再对示踪光斑进行匹配及特征提取,获得床面上同一个示踪光斑点P分别在左右两相机像平面上的像素坐标P1(m1p,n1p)和P2(m2p,n2p),再根据相机像素与CCD物理尺寸的转换关系,将示踪光斑点P在左右相机像平面上的像素坐标分别转化为像平面上的图像坐标P1(x1p,y1p)和P2(x2p,y2p);
4)以左相机焦点F1(0,0,0)为原点建立三维空间坐标系O-XYZ:F1F2为X轴、从点F1向下引垂线为Z轴,过点F1沿河床纵向方向为Y轴,根据数码相机的标定参数f、L、φ,将床面上每个示踪光斑点P在左右两相机像平面上的图像坐标P1(x1p,y1p)和P2(x2p,y2p)按转换公式换算为空间坐标系O-XYZ中的三维空间坐标P(XP,YP,ZP);
激光示踪光斑图像坐标与三维空间坐标的转换公式为:
式中,P1(x1p,y1p)和P2(x2p,y2p)分别为模型床面上激光示踪斑点P在左右两个像平面上的图像坐标;f为左右相机的焦距;L为左右相机的焦点间距;φ为相机偏角;θ1和θ2为虚拟角度,表示模型床面上的点P在平面I上的投影点PT与左右两相机焦点F1与点F2的连线F1PT和F2PT分别与线段F1F2和F2F1的夹角,平面I为过直线F1F2且垂直于模型床面的平面;ψ1为P1相对于左相机像平面中心点O1的偏角∠O1F1K1,表示角∠P1F1K1,θ1表示角∠K1F1F2;同样ψ2为P2相对于右相机像平面中心点O2的偏角∠O2F2K2,θ2为∠K2F2F1;K1和点K2分别为P1点和点P2投影到I平面的点;
5)以床面上所有激光示踪光斑的三维空间坐标生成河工模型地形图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210033188.0A CN102538763B (zh) | 2012-02-14 | 2012-02-14 | 一种河工模型试验三维地形的量测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210033188.0A CN102538763B (zh) | 2012-02-14 | 2012-02-14 | 一种河工模型试验三维地形的量测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102538763A CN102538763A (zh) | 2012-07-04 |
CN102538763B true CN102538763B (zh) | 2014-03-12 |
Family
ID=46346256
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210033188.0A Expired - Fee Related CN102538763B (zh) | 2012-02-14 | 2012-02-14 | 一种河工模型试验三维地形的量测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102538763B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103148842B (zh) * | 2013-02-04 | 2014-11-05 | 国家海洋局第二海洋研究所 | 一种基于遥感图像特征的浅海沙波区多波束测深地形重构方法 |
CN104567796A (zh) * | 2013-10-29 | 2015-04-29 | 江西盛泰光学有限公司 | 一种3d拍摄测距方法 |
CN103617649B (zh) * | 2013-11-05 | 2016-05-11 | 北京江宜科技有限公司 | 一种基于相机自标定技术的河工模型地形测量方法 |
CN104864849B (zh) | 2014-02-24 | 2017-12-26 | 电信科学技术研究院 | 视觉导航方法和装置以及机器人 |
CN103983982A (zh) * | 2014-05-27 | 2014-08-13 | 哈尔滨工业大学 | 车用红外/可见光双摄像头激光雷达 |
CN104267444A (zh) * | 2014-10-23 | 2015-01-07 | 淮海工学院 | 一种孔洞实时三维成像装置 |
CN105423940B (zh) | 2015-12-25 | 2017-12-26 | 同济大学 | 一种地铁隧道结构断面变形快速检测装置 |
CN112710284A (zh) * | 2017-06-16 | 2021-04-27 | 深圳市大疆创新科技有限公司 | 地面形态检测方法、系统和无人机 |
CN110030955A (zh) * | 2018-01-11 | 2019-07-19 | 天津大学 | 一种基于由阴影恢复形状算法的海底粗糙度测量方法 |
DE102018209607A1 (de) * | 2018-06-14 | 2019-12-19 | Volkswagen Aktiengesellschaft | Verfahren und Vorrichtung zum Bestimmen einer Position eines Kraftfahrzeugs |
CN109725640A (zh) * | 2018-12-18 | 2019-05-07 | 天津理工大学 | 一种基于立体视觉与激光光栅的自动驾驶中三维道路重建方法 |
CN111552289B (zh) * | 2020-04-28 | 2021-07-06 | 苏州高之仙自动化科技有限公司 | 检测方法及虚拟雷达装置、电子设备、存储介质 |
CN112066950B (zh) * | 2020-07-24 | 2022-10-14 | 北京空间机电研究所 | 一种多光轴平行的测绘相机单中心投影转换方法 |
CN112461496B (zh) * | 2020-09-09 | 2023-08-18 | 福建省水利水电勘测设计研究院有限公司 | 一种水工及河工模型冲淤试验成果的可视化处理方法 |
CN113280790B (zh) * | 2021-05-07 | 2023-08-15 | 王亚超 | 一种基于激光定位的建筑测绘装置 |
CN113406851B (zh) * | 2021-05-31 | 2023-05-05 | 南方电网调峰调频发电有限公司 | 水工模型安装装置及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1038613C (zh) * | 1994-12-23 | 1998-06-03 | 清华大学 | 河工模型试验流场量测中的图像摄取系统及图象处理方法 |
-
2012
- 2012-02-14 CN CN201210033188.0A patent/CN102538763B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102538763A (zh) | 2012-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102538763B (zh) | 一种河工模型试验三维地形的量测方法 | |
CN102003938B (zh) | 大型高温锻件热态在位检测方法 | |
CN101629806B (zh) | 结合激光发射器的非线性ccd三维定位装置及定位方法 | |
CN100489446C (zh) | 基于相位法的三维轮廓测量方法 | |
CN105931234A (zh) | 一种地面三维激光扫描点云与影像融合及配准的方法 | |
CN104268935A (zh) | 一种基于特征的机载激光点云与影像数据融合系统及方法 | |
CN102519434B (zh) | 一种用于立体视觉三维恢复数据精度测量的试验验证方法 | |
CN102944188B (zh) | 一种点扫描三维形貌测量系统标定方法 | |
JP2009145314A (ja) | 異種センサ統合モデリングによるデジタル写真測量方法及びその装置 | |
CN103971404A (zh) | 一种高性价比的3d实景复制装置 | |
CN102506711B (zh) | 线激光视觉三维旋转扫描方法 | |
CN103617649A (zh) | 一种基于相机自标定技术的河工模型地形测量方法 | |
CN105115560A (zh) | 一种船舱舱容的非接触测量方法 | |
CN102692213A (zh) | 基于主动式全景视觉传感器的交通事故现场测绘仪 | |
CN104034305A (zh) | 一种单目视觉实时定位的方法 | |
CN105466400A (zh) | 一种利用rpc探测多源卫星影像同名像点成像交会角方法 | |
CN105203030A (zh) | 工程场地微小位移量监测方法 | |
Gao et al. | Monocular vision measurement technology applied in dynamic compaction ramming settlement monitoring | |
Grant et al. | The application of an in-line, stereoscopic, PIV system to 3-component velocity measurements | |
CN101770656B (zh) | 一种基于立体正射影像对的大场景立体模型生成及量测方法 | |
Heng Siong et al. | The Application of Smartphone Based Structure From Motion (Sfm) Photogrammetry in Ground Volume Measurement | |
CN116188372B (zh) | 一种基于立体测绘的阻水构筑物识别方法 | |
CN101819038B (zh) | 基于标线参数的交通事故现场二维测量及俯视校正方法 | |
Lingua et al. | Iterative Refraction-Correction Method on Mvs-Sfm for Shallow Stream Bathymetry | |
Stylianidis et al. | A digital close-range photogrammetric technique for monitoring slope displacements |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140312 Termination date: 20180214 |