发明名称 一种河工模型试验三维地形的量测方法
摘要 发明属于水利量测技术领域,特别涉及三维地形的测量方法。在水利工程的河工模型试验中实时监测河道的变形是河工模型试验的关键技术。基于双目立体视觉原理,利用图像摄取技术、激光矩阵光源示踪技术和三维地形重构技术实时量测河工模型试验中床面的三维形态。集成数码摄像机、激光矩形点阵示踪光斑发射装置和计算机为一体,通过设定的系统安装、相机参数校正、床面上激光点阵光斑投影,经济、安全、高质量地实现河工模型床面示踪、床面图像采集、河床激光示踪点提取及匹配,结合计算公式计算出床面上激光示踪光斑的三维空间坐标。基于计算机,以床面上所有激光示踪光斑的三维空间坐标为媒介生成河工模型的地形图。
申请公布号 CN102538763B 申请公布日期 2014.03.12
申请号 CN201210033188.0 申请日期 2012.02.14
申请人 清华大学 发明人 王兴奎;孟震;任海涛;陈启刚;李丹勋
分类号 G01C11/04(2006.01)I 主分类号 G01C11/04(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 廖元秋
主权项 1.一种河工模型试验三维地形的测量方法,其特征在于,采用两台参数相同的数码相机,一台连续激光矩形点阵示踪光斑发射装置,以及分别与两台数码相机相连的一台计算机;所述连续激光矩形点阵示踪光斑发射装置由连续激光器和微调光栅所组成;该微调光栅在河工模型床面上投射出不同间距及大小的矩形点阵示踪光斑,其中矩形点阵示踪光斑中心点为定位点O; 该方法,包括以下步骤: 1)首先进行相机焦距f、两相机两焦点距离L以及相机偏角φ标定: 1-1)相机焦距f标定计算公式为: <img file="FDA00003610009300011.GIF" wi="284" he="148" />式中,a为CCD较短一边的边长,y<sub>2</sub>指相机初始拍摄垂向最大高度,y<sub>1</sub>为相机沿其光轴轴向向后移动ΔX后所拍摄的垂向最大高度,y<sub>1</sub>>y<sub>2</sub>; 1-2)两焦点距离L标定为:<img file="FDA00003610009300012.GIF" wi="530" he="148" />式中,L<sub>1</sub>为OB<sub>1</sub>的长度,O为定位点,B<sub>1</sub>为与点B在左相机成像面上具有相同像素坐标的床面上的一点,点B为在点O垂直向上y高度处的一点;角度α<sub>1</sub>为点O与左相机焦点F<sub>1</sub>的连线F<sub>1</sub>O与直线F<sub>1</sub>B<sub>1</sub>的夹角,角度β<sub>1</sub>为点O与点B<sub>1</sub>的连线B<sub>1</sub>O与直线B<sub>1</sub>F<sub>1</sub>的夹角,其中α<sub>1</sub>和β<sub>1</sub>在标定时计算出α<sub>1</sub>=arctan(O<sub>1</sub>P<sub>1</sub>/f),β<sub>1</sub>=arctan(y/L<sub>1</sub>),O<sub>1</sub>P<sub>1</sub>为床面上线段OB<sub>1</sub>在左相机像平面上的图像长度; 1-3)相机偏角φ标定为: φ=α<sub>1</sub>+β<sub>1</sub>式中,α<sub>1</sub>为直线F<sub>1</sub>O与直线F<sub>1</sub>B<sub>1</sub>之间的夹角;β<sub>1</sub>为直线B<sub>1</sub>O与直线B<sub>1</sub>F<sub>1</sub>之间的夹角; 2)照相机标定完成后,在两相机焦点中间位置垂直向下安装连续激光矩形点阵示踪光斑发射装置,激光示踪光斑中心点为定位点O,定位点O在左右两相机像平面的成像分别为两相机像平面的中心点O<sub>1</sub>和O<sub>2</sub>;左右两相机成像面以中心点O<sub>1</sub>和O<sub>2</sub>为原点分别建立各自的直角坐标系x<sub>1</sub>O<sub>1</sub>y<sub>1</sub>和x<sub>2</sub>O<sub>2</sub>y<sub>2</sub>,床面上也以点O为中心点建立坐标系XOY,将两相机像平面及床面均分成四个象限,且每个象限所包含或所成像的激光示踪斑点一一对应; 3)对床面进行拍照,分别得到左右两相机所拍摄的系列照片;对所述系列照片进行噪声处理、畸变处理、像素二值化及对示踪光斑边缘检测的常规预处理,以及进行示踪光斑中心坐标的计算,再对示踪光斑进行匹配及特征提取,获得床面上同一个示踪光斑点P分别在左右两相机像平面上的像素坐标P<sub>1</sub>(m<sub>1p</sub>,n<sub>1p</sub>)和P<sub>2</sub>(m<sub>2p</sub>,n<sub>2p</sub>),再根据相机像素与CCD物理尺寸的转换关系,将示踪光斑点P在左右相机像平面上的像素坐标分别转化为像平面上的图像坐标P<sub>1</sub>(x<sub>1p</sub>,y<sub>1p</sub>)和P<sub>2</sub>(x<sub>2p</sub>,y<sub>2p</sub>); 4)以左相机焦点F<sub>1</sub>(0,0,0)为原点建立三维空间坐标系O-XYZ:F<sub>1</sub>F<sub>2</sub>为X轴、从点F<sub>1</sub>向下引垂线为Z轴,过点F<sub>1</sub>沿河床纵向方向为Y轴,根据数码相机的标定参数f、L、φ,将床面上每个示踪光斑点P在左右两相机像平面上的图像坐标P<sub>1</sub>(x<sub>1p</sub>,y<sub>1p</sub>)和P<sub>2</sub>(x<sub>2p</sub>,y<sub>2p</sub>)按转换公式换算为空间坐标系O-XYZ中的三维空间坐标P(X<sub>P</sub>,Y<sub>P</sub>,Z<sub>P</sub>); 激光示踪光斑图像坐标与三维空间坐标的转换公式为: <img file="FDA00003610009300021.GIF" wi="622" he="454" />其中<img file="FDA00003610009300022.GIF" wi="604" he="509" />式中,P<sub>1</sub>(x<sub>1p</sub>,y<sub>1p</sub>)和P<sub>2</sub>(x<sub>2p</sub>,y<sub>2p</sub>)分别为模型床面上激光示踪斑点P在左右两个像平面上的图像坐标;f为左右相机的焦距;L为左右相机的焦点间距;φ为相机偏角;θ<sub>1</sub>和θ<sub>2</sub>为虚拟角度,表示模型床面上的点P在平面I上的投影点P<sub>T</sub>与左右两相机焦点F<sub>1</sub>与点F<sub>2</sub>的连线F<sub>1</sub>P<sub>T</sub>和F<sub>2</sub>P<sub>T</sub>分别与线段F<sub>1</sub>F<sub>2</sub>和F<sub>2</sub>F<sub>1</sub>的夹角,平面I为过直线F<sub>1</sub>F<sub>2</sub>且垂直于模型床面的平面;ψ<sub>1</sub>为P<sub>1</sub>相对于左相机像平面中心点O<sub>1</sub>的偏角∠O<sub>1</sub>F<sub>1</sub>K<sub>1</sub>,<img file="FDA00003610009300023.GIF" wi="55" he="52" />表示角∠P<sub>1</sub>F<sub>1</sub>K<sub>1</sub>,θ<sub>1</sub>表示角∠K<sub>1</sub>F<sub>1</sub>F<sub>2</sub>;同样ψ<sub>2</sub>为P<sub>2</sub>相对于右相机像平面中心点O<sub>2</sub>的偏角∠O<sub>2</sub>F<sub>2</sub>K<sub>2</sub>,θ<sub>2</sub>为∠K<sub>2</sub>F<sub>2</sub>F<sub>1</sub>;K<sub>1</sub>和点K<sub>2</sub>分别为P<sub>1</sub>点和点P<sub>2</sub>投影到I平面的点;5)以床面上所有激光示踪光斑的三维空间坐标<img file="FDA000036100093000212.GIF" wi="298" he="77" />生成河工模型地形图。
地址 100084 北京市海淀区清华园1号