摘要
桥梁悬臂施工控制过程需进行大量线形测量工作,传统测量方法耗时耗力、难以通测,无人机摄影测量具有高机动性、低成本、不依赖人工等优势,将其应用于桥梁线形测量,可解决自适应施工控制过程中线形通测的难题。考虑桥梁施工测量的特点,提出桥梁线形无人机测量的技术框架和主要控制参数,研究适合施工现场的标记点布置方式与航线规划。通过试验研究无人机飞行高度、投射线夹角、视角数量、照片数量等参数对测量精度的影响,在此基础上,总结出桥梁悬臂施工线形无人机测量方案优化方法,可保证各测点的坐标测量精度达到毫米级。
悬臂拼装法施工速度快、质量易控制,是最常用的大跨径桥梁施工方法。但在施工过程中必须对桥梁线形及受力状态进行测控,使结构实际线形与预期状态之间的误差在容许范围内,以保障拼装精度和成桥质量。大跨度桥梁悬臂施工的线形控制工作量大、精度要求高、难度大,同时施工控制误差的影响因素较多,自适应控制法通过在施工控制中引入参数识别,为主动进行结构施工状态的准确控制提供决策依据,可获得高质量的施工控制结
自适应施工控制过程中,从反馈计算和参数识别的角度出发,希望每个工况内测量的节点数和次数越多越好,大大增加了测量的工作
随着智能全站仪、电子水准仪、全球卫星导航系统(global navigation satellite system, GNSS)、激光扫描、摄影测量等测绘技术的发展,新兴桥梁线形测量方法得到了快速发
摄影测量技术通过量测相机摄取的二维图像,测定物体在三维空间的位置、形状、大小乃至物体的运动。近百年来,摄影测量经历了模拟、解析与数字摄影测量三个阶
无人机(unmanned aerial vehicle, UAV)具有高机动性,可到达人类难以抵达的观测点,快速获取多视角数据,近年来无人机导航、障碍物检测与回避、空中决策等算法和计算机视觉技术发展迅
综上所述,无人机摄影测量具有高机动性、低成本、不依赖人工等优势,将其应用于桥梁线形测量,可解决自适应施工控制过程中线形通测的难题。但对于桥梁施工线形测量这种尺度大、测点多、精度要求高的任务,目前还未有研究探索无人机测量的可行性、方法和精度水平。本文考虑桥梁施工测量的特点,提出桥梁线形无人机测量的技术框架和主要控制参数,研究适合施工现场的标记点布置方式与航线规划,并通过试验定量研究各控制参数对测量精度的影响,在此基础上建立无人机测量方案的优化方法。
摄影测量利用立体视觉原理进行三维重构,结合图像处理、射影几何和统计学方法,解析恢复不同角度图像上特征点的空间交会光线,并通过光束法平差求解特征点的三维坐标。相机成像是摄影测量的第一步,可分为4个步骤:刚体变换、透视投影、畸变校正和数字化图像,涉及到4个坐标系的变化:像素坐标系OUV、图像坐标系OpXpYp、相机坐标系OcXcYcZc、世界坐标系OwXwYwZ

图1 线性相机模型
Fig. 1 Linear camera model
(1) |
式中:(Xw, Yw, Zw)为相机坐标系中空间点的坐标;(u, v)为像素坐标系中像点的坐标;Zc表示空间点到相机坐标系原点距离;K为相机内参矩阵;(u0, v0)为像片中心点;dx和dy为单个像素在图像坐标系u轴和v轴方向的长度;f为相机焦距;R、t为相机坐标系变换到世界坐标系的旋转矩阵和平移向量。
摄影测量常在被测物表面设置专门测量标记点,由外围编码环段和中心圆点组成,编码环段用于计算机快速准确地识别,中心圆点用于坐标定位,如

图2 摄影测量标记点
Fig. 2 Photogrammetric markers
识别并匹配各图像中的标记点后,首先采用五点算法(five-point algorithm)求解基本矩阵,定向每张图像对应的相机位置,然后利用三角测量(triangulation)解算出各标记点的三维坐

图3 根据双目视觉原理计算空间中一点坐标
Fig. 3 Calculation of the coordinates of a point according to the principle of binocular vision
根据相机成像模型,可计算出空间点的三维坐标如下:
(2) |
(3) |
(4) |
式中:(x1,y1)和(x2,y2)为像点1、像点2在图像坐标系O1x1y1和O2x2y2中的坐标;(u,v)为像点1在其像素坐标系中的坐标;dx、dy为单位像素在像素坐标系u轴、v轴方向的长度;Zc为空间点到相机光轴中心的距离;f为相机焦距;r1,r2,…,r9为两个相机坐标系变换矩阵的元素。
为了提高测量精度,通常需要更多角度的图像进行迭代计算以减小误差,如
(5) |
式中:X为满足条件的空间点最优三维坐标;ui为该点在第i幅图像上的像点坐标;为推算出的X投影到图像上的重投影像点坐标。
(6) |
式中:ε为重投影误差的最小值;n为空间点的个数;m为相机位置的个数;Q(bi)为预测的第i个点;xij为第j张图像上的第i个特征点;d(x, y)表示图像上由不同向量x和y表示的点的欧式距离。
根据摄影测量原理,飞行高度H、投射线夹角、视角数量、照片数量等控制参数对精度有较大影响。由式(

图4 不同投射线夹角与点定位精度关系
Fig. 4 Relationship between angle of bundle and point positioning accuracy
摄影测量得到的是局部测量坐标系下的测点坐标,需引入基准点进行刚性坐标变换,将其换算到同一桥梁坐标系中,得到测点的实际标高,才能用于施工控制。坐标系转换可采用间接平差方

图5 引入基准点进行测点坐标系转换
Fig. 5 Introducetion of datum point in transformation of the coordinate system of measuring point
将无人机应用于桥梁悬臂施工控制测量,首先需根据结构特点设计标记点布置方案,然后综合考虑测量精度和效率,确定飞行高度H、投射线夹角、视角数量、照片数量等控制参数,并规划无人机航线。执行无人机测量任务得到图像数据,通过畸变校正、标记点识别和三维坐标求解得到各标记点的坐标,最后关联基准点进行坐标转换,得到各测点在桥梁坐标系下的真实标高,方法框架见

图6 桥梁线形无人机摄影测量方法框架图
Fig. 6 Framework of bridge alignment UAV photogrammetry method
室内摄影测量标记点主要用于设置比例尺或辅助识别被测物轮廓,多密集布置在测区角点、边缘等位置。但测量主梁节段的标高和线形,标记点宜布置在节段表面中轴线和翼缘根部,且布置过密会影响测量便利性,需根据主梁节段的尺度研究合理的标记点布置方式。
针对桥梁施工线形测量的特点,设计了几种标记点布置方案,见

图7 标记点布置方案
Fig. 7 Layout plan of marking points
方案一、二、三均可成功解算出标记点坐标,方案四由于标记点间距过大、数量太少,解算失败。由
预设航线可提高测量效率,对于直线型桥梁,将测区简化为矩形,无人机需到达图中所有预设拍摄位置,效率较高的航线包括巡回航线和S形航线,如

图8 巡回航线
Fig. 8 Itinerant route of UAV

图9 S形航线
Fig. 9 S-shaped route of UAV
从航线长度、转弯次数、镜头旋转次数、飞行时间等方面比较两种航线效果,见
根据无人机摄影测量原理,飞行高度H、投射线夹角、视角数量、照片数量是精度的主要控制参数,为定量分析各控制参数的影响,开展了一系列对照试验。在地面上布置一组标记点,如

图10 标记点排列方式
Fig. 10 Layout plan of marking points in the experiment
控制飞行高度在2.7~7.7 m范围内,分别选择8个视角和4个视角,令无人机在各位置悬停拍摄4张照片,各飞行高度下的测量精度见

图11 无人机飞行高度对精度的影响
Fig. 11 Influence of UAV flight height on accuracy
通过试验发现,无人机摄影测量的精度水平与单位像素真实尺度Zcdx、Zcdy大致相当,当相机分辨率、传感器等参数固定时,精度水平直接与飞行高度相关。对于本试验所用的无人机,要将最大误差控制在1 mm以内,建议的飞行高度为5 m。可以通过提高相机分辨率、传感器画幅尺寸等方式进一步提高适用的飞行高度,拓展应用空间。
两次拍摄相机镜头中心点与物点连线的夹角被称为投射线夹角,如

图12 投射线夹角与云台倾角的关系
Fig. 12 Relationship between angle of bundle and camera inclination
控制投射线夹角从10°增大到120°,测量精度的变化见

图13 摄影测量误差与投射线夹角关系
Fig. 13 Relationship between photogrammetry error and angle of bundle
相机视角数量越多,光束法平差的结果越接近点的真实位置,精度越好。但是减少相机视角数量可提高拍摄效率,且平行或垂直于测量轴线方向的视角更便于无人机规划航线自动拍摄。在同一拍摄高度,保证投射线夹角处于高精度区间,分别采用2~8个相机视角进行拍摄,如

图14 不同相机视角拍摄方案
Fig. 14 Different camera position schemes
各视角选择4张照片进行重建,不同方案的测量精度见

图15 不同视角数量对精度的影响
Fig. 15 Influence of different number of viewing angles on accuracy
通过拍摄多张重复照片可提高测量冗余度,降低无人机悬停振动的影响。控制飞行高度不变,选择8个视角和4个视角,各位置分别选择1~5张重复照片进行重建,测量精度见

图16 照片数量对精度的影响
Fig. 16 Influence of photo quantity on accuracy
各位置只有1张照片时,所选照片质量对结果的影响非常大,标记点识别成功率低,测量误差较大。当照片数量为2张以上时,重建成功率大大提高,最大相对误差也降低到0.10 %以内,RMS误差可控制在0.6 mm以内。当照片数量达到3张后,继续增加照片对精度的影响很小。另外,约3~4 s可完成一张照片的标记点识别、定向和坐标解算,总计算时长和照片数量成正相关,对于4个视角数量方案,取3张照片时,30 s左右可完成坐标解算。
无人机进行线形测量任务时,主要控制参数包括飞行高度、投射线夹角(云台倾角)、相机位置和照片数量。根据第3节的研究,无人机摄影测量的精度水平和单位像素真实尺度大致相当,已知相机参数和测量精度要求时,可根据相机成像原理反推出可行的飞行高度区间。
完成相机标定后,图像分辨率和镜头标定区域占比均为已知,假设镜头焦距为f, mm,图像分辨率为I1×I2(单位:像素),传感器尺寸为W×H(单位:mm),镜头标定区域占比为c,见

图17 无人机飞行高度确定原理示意图
Fig. 17 Key parameters of UAV measurement scheme
假设待测视域为一a×b(单位:m)的矩形,考虑镜头标定区域占比,相机需要达到的视域范围为 (a×b)/c,故飞行高度h需满足
(7) |
根据线性相机模型,飞行高度h和单位像素真实尺度Zcdx、Zcdy的关系见
(8) |
(9) |
(10) |
改变无人机相机云台倾角可调整投射线夹角和投射线倾角(
(11) |

图18 云台倾角与投射线倾角的关系
Fig. 18 Relationship between camera inclination angle and bundle inclination
增加相机位置可提高测量精度,但当增加到4个相机位置之后,提高的幅度较小。综合考虑测量效率、计算时长和无人机航线规划便利性,选择4个不同相机位置采集照片。
为缓解无人机悬停振动对成像质量的影响,保证重建冗余度和测量精度,可在各相机位置采集3张重复照片。
为验证本节所述优化方法的有效性,选择了一段28.5 m的跑道进行了模型试验。测量标记点布置方案如

图19 模型试验标记点布置方案(单位:mm)
Fig. 19 Layout plan of marking points for model experiment (unit: mm)

图20 各测点标高测量误差
Fig. 20 Elevation measurement error of each marking point
本文考虑桥梁施工控制测量的特点和无人机摄影测量原理,提出了桥梁线形无人机测量的技术框架和测量控制参数。主要研究结论包括:
(1)在桥面上布设标记点可提高重构成功率和测量精度,标记点宜布置在节段中轴线和翼缘根部,测点间距可设置为1.5 m。对于直线桥梁,巡回航线耗时短、耗电量低、效率高,可实现自动化测量。
(2)无人机进行桥梁施工线形测量需控制的技术参数包括飞行高度H、投射线夹角θ、视角数量、照片数量等。飞行高度对测量精度有直接影响,对于本无人机,高度不超过5 m时,精度可控制在1 mm以内;投射线夹角与平差效果密切相关,当40°≤θ≤110°时,标记点识别和图像定向成功率都较高;4个视角数量可得到较好的测量精度,也便于无人机航线规划;采集3张重复照片,可保证冗余度和计算速度。
(3)建立了桥梁悬臂施工线形无人机测量控制参数优化方法,无人机摄影测量的精度水平和单位像素真实尺度Zcdx、Zcdy大致相当,已知相机参数和精度要求,可根据成像原理反推出可行的飞行高度区间。确定相机云台最优倾角范围需综合考虑投射线夹角和倾角的要求,控制光束平差和标记点识别精度较高。采用优化后的摄影测量方案进行了模型试验,可将误差控制在2 mm以内。
作者贡献声明
石雪飞:研究理念、研究方法、论文撰写。
许琪:研究方法、数据建模与分析、试验与分析、论文撰写。
马海英:研究理念、数据建模与分析、论文撰写。
参考文献
李乔, 卜一之, 张清华. 基于几何控制的全过程自适应施工控制系统研究[J]. 土木工程学报, 2009, 42(7): 69. [百度学术]
LI Qiao, BU Yizhi, ZHANG Qinghua. Whole-procedure adaptive construction control system based on geometry control method[J]. China Civil Engineering Journal, 2009, 42(7): 69. [百度学术]
石雪飞, 项海帆. 斜拉桥施工控制方法的分类分析[J]. 同济大学学报(自然科学版), 2001, 29(1): 55. [百度学术]
SHI Xuefei, XIANG Haifan. Analysis of construction control systems of cable-stayed bridge[J]. Journal of Tongji University (Natural Science), 2001,29(1): 55. [百度学术]
吴迪军. 桥梁工程测量技术现状及发展方向[J]. 测绘通报, 2016(1): 1. [百度学术]
WU Dijun. Current status and development trend of bridge engineering surveying techniques[J]. Bulletin of Surveying and Mapping, 2016(1):1. [百度学术]
余加勇, 邵旭东, 晏班夫, 等. 基于全球导航卫星系统的桥梁健康监测方法研究进展[J]. 中国公路学报, 2016, 29(4): 30. [百度学术]
YU Jiayong, SHAO Xudong, YAN Banfu, et al. Research and development on global navigation satellite system technology for bridge health monitoring[J]. China Journal of Highway and Transport, 2016, 29(4): 30. [百度学术]
吴家铖, 刘成龙, 畅卫杰, 等. 特大型悬索桥运营期维护测量内容与方法研究[J]. 测绘与空间地理信息, 2020, 43(10): 52. [百度学术]
WU Jiacheng, LIU Chenglong, CHANG Weijie, et al. Research on contents and methods of maintenance survey of extra large suspension bridges in operation period[J]. Geomatics & Spatial Information Technology, 2020, 43(10): 52. [百度学术]
李德仁. 摄影测量与遥感的现状及发展趋势[J]. 武汉测绘科技大学学报, 2000, 25(1): 1. [百度学术]
LI Deren. Towards photogrammetry and remote sensing: status and future development[J]. Journal of Wuhan Technical University of Surveying and Mapping, 2000, 25(1): 1. [百度学术]
OHNISHI Y, NISHIYAMA S, YANO T, et al. A study of the application of digital photogrammetry to slope monitoring systems[J]. International Journal of Rock Mechanics and Mining Sciences, 2006, 43(5): 756. [百度学术]
唐君萍, 张丽艳, 刘胜兰, 等.飞机复杂零件上大量小尺寸导孔的快速视觉检测[J]. 中国机械工程, 2015, 26(18): 2456. [百度学术]
TANG Junping, ZHANG Liyan, LIU Shenglan, et al. Quick vision inspection for large number of small size guide holes on complex aircraft parts[J]. China Mechanical Engineering, 2015, 26(18): 2456. [百度学术]
周颖, 刘彤. 基于计算机视觉的混凝土裂缝识别[J]. 同济大学学报(自然科学版), 2019, 47(9): 1277. [百度学术]
ZHOU Ying, LIU Tong. Computer vision-based crack detection and measurement on concrete structure[J]. Journal of Tongji University (Natural Science), 2019, 47(9): 1277. [百度学术]
AI-KAFF A, MARTÍN D, GARCÍA F, et al. Survey of computer vision algorithms and applications for unmanned aerial vehicles[J]. Expert Systems with Application, 2018, 92, 447. [百度学术]
LIU P, CHEN A Y, HUANG Y N, et al. A review of rotorcraft unmanned aerial vehicle (UAV) developments and applications in civil engineering[J]. Smart Structures and Systems, 2014, 13(6): 1065. [百度学术]
HE F, ZHOU T, XIONG W, et al. Automated aerial triangulation for UAV-based mapping[J]. Remote Sensing, 2018, 10(12), 1952. [百度学术]
DELGADO J P, SORIA P R, ARRUE B C, et al. Bridge mapping for inspection using an UAV assisted by a total station[C]// Iberian Robotics Conference. Seville: Springer, 2017: 309-319. [百度学术]
MORGENTHAL G, HALLERMANN N. Quality assessment of unmanned aerial vehicle (UAV) based visual inspection of structures[J]. Advances in Structural Engineering, 2014, 17(3): 289. [百度学术]
HABIB A, XIONG W, HE F, et al. Improving orthorectification of UAV-based Push-broom scanner imagery using derived orthophotos from frame cameras[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 10(1):1. [百度学术]
ELLENBERG A, BRANCO L, KRICK A, et al. Use of unmanned aerial vehicle for quantitative infrastructure evaluation[J]. Journal of Infrastructure Systems, 2014, 4(1): 40. [百度学术]
刘春, 曾劲涛, 张书航, 等. 面向单体异形建筑的无人机单相机实景三维建模[J]. 同济大学学报(自然科学版), 2018, 46(4): 550. [百度学术]
LIU Chun, ZENG Jintao, ZHANG Shuhang, et al. True 3D modelling towards a special-shaped building unit by unmanned aerial vehicle with a single camera[J]. Journal of Tongji University (Natural Science), 2018, 46(4): 550. [百度学术]
YOON H, SHIN J, SPENCER B F. Structural displacement measurement using an unmanned aerial system[J]. Computer Aided Civil and Infrastructure Engineering, 2018, 33(3): 183. [百度学术]
马颂德, 张正友. 计算机视觉: 计算理论与算法基础[M]. 北京: 科学出版社, 1998. [百度学术]
MA Songde, ZHANG Zhengyou. Computer vision: computational theory and algorithm foundation[M]. Beijing: Science Press, 1998. [百度学术]
NISTER D. An efficient solution to the five-point relative pose problem[J]. IEEE Transactions on Pattern Analysis & Machine Intelligence, 2004, 26(6): 756. [百度学术]
LOURAKIS M. SBA: A software package for generic sparse bundle adjustment[J]. ACM Transactions on Mathematical Software, 2009, 36(1): 2. [百度学术]