CN112629431A

CN112629431A - 土木结构变形监测方法及相关设备

Info

Publication number: CN112629431A
Application number: CN202011468740.XA
Authority: CN
Inventors: 刘肖琳; 丁晓华; 于起峰; 张跃强; 尹义贺
Original assignee: Shenzhen University
Current assignee: Shenzhen University
Priority date: 2020-12-14
Filing date: 2020-12-14
Publication date: 2021-04-09
Anticipated expiration: 2040-12-14
Also published as: CN112629431B

Abstract

本发明实施例公开一种土木结构变形监测方法及相关设备，土木结构变形监测系统包括标志点，以及设置在同一相机平台的相机单元，标志点包括至少三个基准点、被测土木结构的至少一个待测点；先根据拍摄图像获取至少三个基准点和待测点在第一时间点和第二时间点的像点坐标；再根据至少三个基准点在第一时间点和第二时间点的像点坐标确定相机平台的六自由度变化量；最后，根据六自由度变化量、待测点在第一时间点和第二时间点的像点坐标确定待测点在相机平台坐标系中至少一个方向的位移量。利用上述变形监测方法，不仅可以做到对土木结构的动态变形监测，也可以减少监测成本；而且可以有效提高变形监测测量结果的精度。

Description

土木结构变形监测方法及相关设备

技术领域

本发明涉及监测技术领域，尤其涉及一种土木结构变形监测方法及相关设备。

背景技术

土木结构在长期使用过程中，往往会发生由列车高速行驶、船舶撞击、地震等动荷载作用而引起的高频动态变形；和温度荷载、邻近结构施工的扰动而引起的缓慢准静态变形。因此，对土木结构的动、静态变形进行高精度监测是土木结构健康监测的重要内容之一，也是对土木结构性能进行准确评估和安全预警的必要前提。

现有技术中，结构变形测量是利用专用测量仪器和方法对待测结构的位置或形状变化进行观测，获得待测结构随时间的变形特征，主要有以下几种测量方法。第一种是接触式测量方法，例如位移计、加速度计等。第二种是利用GPS进行变形监测，其优点是可以实现全天候、无人值守，实时或准实时作业。第三种是利用传统的光学测量技术，例如水准仪、光学经纬仪、电子全站仪等仪器，具有非接触测量、精度高等优点。第四种是在普通全站仪上安装自动马达，形成全自动全站仪(测量机器人)，可以实现测量过程的自动化，观测周期也可大幅缩减。

然而，接触式测量方法会对测量对象本身造成干扰，而且测量系统布置繁琐、测量点有限，无法满足对大型结构形貌和变形进行全场、高精度测量以及动态监测的需求。而基于GPS的测量方法存在如下缺点：一是由于每个观测点都需要布设接收机天线，对大量点的测量成本较高；二是测量的前提条件是必须能够接收到GPS信号，所以一般只能应用在露天情况，而难以实现室内或地下作业；三是在铅垂方向上测量精度往往不能满足工程的需要。另外，传统的光学测量技术也有自身的局限性，如难以实现测量过程的自动化，观测周期长，观测受外界环境条件限制多等。最后，常用的全自动全站仪测量手段存在以下缺点：一是测量周期长，需要运动对准，测量一次最快也是按数十分钟或小时计，无法满足高频次的动态测量需求；二是无法实时修正测量平台自身不稳定的影响，无法实现长期有效稳定测量。

发明内容

本发明实施例提供了一种土木结构变形监测方法及相关设备，不仅可以做到对土木结构的动态、静态及准静态等变形监测，也可以减少监测成本；而且可以有效提高变形监测测量结果的精度。

第一方面，本发明实施例提供了一种土木结构变形监测方法，应用于土木结构变形监测系统，所述监测系统包括标志点，以及设置在同一相机平台的用于在第一时间点和第二时间点拍摄所述标志点得到拍摄图像的相机单元，所述标志点包括至少三个基准点、被测土木结构的至少一个待测点；所述方法包括：

根据所述拍摄图像获取所述至少三个基准点和所述待测点在第一时间点和第二时间点的像点坐标；

根据所述至少三个基准点在第一时间点和第二时间点的像点坐标确定所述相机平台的六自由度变化量；

根据所述六自由度变化量、所述待测点在第一时间点和第二时间点的像点坐标确定所述待测点在相机平台坐标系中至少一个方向的位移量。

可选地，所述监测系统还包括设置在所述相机平台上的、用于测量所述相机单元和所述待测点之间的距离的测距模块，所述测距模块与所述相机单元并排设置，且所述测距模块的测距方向和所述相机单元的拍摄方向相同；

所述根据所述六自由度变化量、所述待测点在第一时间点和第二时间点的像点坐标确定所述待测点在相机平台坐标系中至少一个方向的位移量，包括：

所述相机单元的同一相机视场中有且只有一个目标待测点时，根据所述六自由度变化量、所述目标待测点在第一时间点和第二时间点的像点坐标、以及成像几何约束关系，确定所述目标待测点在相机平台坐标系中的第一方向和第二方向上的位移量；

获取对应第一时间点和第二时间点的，所述目标待测点和所述相机单元之间的距离，并根据所述距离和所述六自由度变化量确定所述目标待测点在相机平台坐标系中的第三方向的位移量。

可选地，所述根据所述六自由度变化量、所述待测点在第一时间点和第二时间点的像点坐标确定所述待测点在相机平台坐标系中至少一个方向的位移量，包括：

所述相机单元的同一相机视场中有两个以上的目标待测点时，将所述两个以上的目标待测点的变形运动建模为局部平移运动；

根据所述六自由度变化量、所述目标待测点在第一时间点和第二时间点的像点坐标、以及所述局部平移运动对应的成像几何约束关系，确定所述目标待测点在相机平台坐标系中的第一方向、第二方向、第三方向上的位移量。

所述相机单元的同一相机视场中有三个以上的目标待测点时，将所述三个以上的目标待测点的变形运动建模为局部刚体运动；

根据所述六自由度变化量、所述目标待测点在第一时间点和第二时间点的像点坐标、以及所述局部刚体运动对应的成像几何约束关系，确定所述目标待测点在相机平台坐标系中的第一方向、第二方向、第三方向上的位移量。

可选地，所述根据所述拍摄图像获取所述至少三个基准点和所述待测点在第一时间点和第二时间点的像点坐标，包括：

获取在所述第一时间点和/或所述第二时间点拍摄的包含所述基准点和/或所述待测点的拍摄图像；

对所述拍摄图像进行图像像点定位处理，得到所述基准点和/或所述待测点在所述第一时间点和/或所述第二时间点的像点坐标。

可选地，在对所述拍摄图像进行图像像点定位处理之前，所述方法还包括：

对所述拍摄图像进行滤波处理，以滤除所述拍摄图像中的干扰像素。

第二方面，本发明实施例提供了一种土木结构变形监测系统，包括标志点、处理装置以及设置在同一相机平台的用于获取所述标志点的像点坐标的相机单元，所述标志点包括至少三个基准点、被测土木结构的至少一个待测点；

所述相机单元，用于获取在第一时间点和第二时间点拍摄所述标志点得到拍摄图像；

所述处理装置包括：

获取单元，用于根据所述拍摄图像获取所述至少三个基准点和所述待测点在第一时间点和第二时间点的像点坐标；

第一确定单元，用于根据所述至少三个基准点在第一时间点和第二时间点的像点坐标确定所述相机平台的六自由度变化量；

第二确定单元，用于根据所述六自由度变化量、所述待测点在第一时间点和第二时间点的像点坐标确定所述待测点在相机平台坐标系中至少一个方向的位移量。

可选地，所述相机单元包括至少一个相机。

第三方面，本发明实施例提供了一种土木结构变形监测设备，包括：处理器和存储器；

所述处理器和存储器相连，其中，所述存储器用于存储程序代码，所述处理器用于调用所述程序代码，以执行如第一方面所述的土木结构变形监测方法。

第四方面，本发明实施例提供了一种计算机存储介质，所述计算机存储介质存储有计算机程序，所述计算机程序包括程序指令，所述程序指令当被处理器执行时，执行如第一方面所述的土木结构变形监测方法。

本发明实施例中，土木结构变形监测系统包括标志点，以及设置在同一相机平台的用于在第一时间点和第二时间点拍摄标志点得到拍摄图像的相机单元，其中，标志点包括至少三个基准点、被测土木结构的至少一个待测点；变形监测方法中，先根据拍摄图像获取至少三个基准点和待测点在第一时间点和第二时间点的像点坐标；再根据至少三个基准点在第一时间点和第二时间点的像点坐标确定相机平台的六自由度变化量；最后，根据六自由度变化量、待测点在第一时间点和第二时间点的像点坐标确定待测点在相机平台坐标系中至少一个方向的位移量。利用上述变形监测方法，基于相机单元的拍摄图像进行变形监测，不仅可以做到对土木结构的动态变形监测，也可以减少监测成本；而利用相机平台的六自由度变化量确定待测点在相机平台坐标系中至少一个方向的位移量，可以有效提高变形监测测量结果的精度。

附图说明

为了更清楚地说明本发明实施例技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明实施例提供的一种土木结构变形监测方法的具体流程示意图；

图2是本发明实施例提供的一种土木结构变形监测方法的流程示意图；

图3a、图3b是本发明实施例提供的一种大气抖动对变形监测的影响示意图；

图4是本发明实施例提供的一种土木结构变形监测系统的结构示意图；

图5a、图5b是本发明实施例提供的一种单校正相机和双校正相机的观测误差示意图；

图6a、图6b、图6c、图6d是本发明实施例提供的一种土木结构变形监测系统的结构示意图；

图7是本发明实施例提供的一种单相机摄像测量基本原理示意图；

图8是本发明实施例提供的一种土木结构变形监测系统的处理装置的结构示意图；

图9是本发明实施例提供的一种土木结构变形监测设备的结构示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。

应当理解，本申请的说明书和权利要求书及附图中的术语“第一”、“第二”等是用于区别不同对象，而不是用于描述特定顺序。此外，术语“包括”和“具有”以及它们任何变形，意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可选地还包括没有列出的步骤或单元，或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤或单元。

在本发明中提及“实施例”意味着，结合实施例描述的特定特征、结构或特性可以包含在本发明的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例，也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是，本发明所描述的实施例可以与其它实施例相结合。

一种土木结构变形监测方法，应用于土木结构变形监测系统，本申请中，土木结构包括桥涵、路基、隧道、大坝、高层建筑等。土木结构变形监测系统包括标志点，以及设置在同一相机平台的用于在第一时间点和第二时间点拍摄标志点得到拍摄图像的相机单元，其中，标志点包括至少三个基准点、被测土木结构的至少一个待测点。具体的基准点和待测点的个数和设置位置可以根据实际需要进行设置。

具体地，基准点为位置固定或者运动已知的标志点，而待测点为被测土木结构上需要测量的点位。实际监测时，在被测土木结构的监测作业范围内安装相机平台，以设置相机单元。相机平台也即相机观测站，观测站就是一个用来观测的站点，相机平台就是一个可以放相机的台子，或者墩子，或者一个三脚架等等。而相机单元包括至少一个相机，相机的焦距和数量可根据测量精度要求和实际工况进行灵活调配。调整相机朝向方位和镜头焦距，使得周围的基准点和被测土木结构上的待测点可清晰地出现在相应的相机视场中。另外，参考图1，图1是本发明实施例提供的一种土木结构变形监测方法的具体流程示意图；在对变形监测单元进行初始化之后，相机单元在第一时间点和第二时间点同步跟踪采集基准点和待测点的图像；第一时间点和第二时间点的具体时刻可以根据实际需要进行设定，以充分满足用户的监测要求，例如，第一时间点早于第二时间点，则第一时间点可以为初始时刻或者用户指定的任意时刻，第二时间点也可以为用户指定的晚于第一时间点的任意时刻。

请参见图2，是本发明实施例提供的一种土木结构变形监测方法的流程示意图；变形监测方法包括：

步骤201，根据拍摄图像获取至少三个基准点和待测点在第一时间点和第二时间点的像点坐标；

具体地，利用相机单元在第一时间点和第二时间点同步采集的基准点和待测点的拍摄图像，即在第一时间点同步采集的包括基准点和/或待测点的图像，这些图像是在第一时间点同时采集的，可以为仅包括基准点的图像、仅包括待测点的图像、包括基准点和待测点的图像。同样地，拍摄图像还包括在第二时间点同步采集包括基准点和/或待测点的图像。处理上述拍摄图像以得到监测系统中的基准点和待测点分别在第一时间点和第二时间点时的像点坐标，像点坐标为像点在图像坐标系中的位置。

步骤202，根据至少三个基准点在第一时间点和第二时间点的像点坐标确定相机平台的六自由度变化量；

具体地，本申请中，利用基准点来获得相机平台的六自由度运动的变化量，即六自由度变化量。具体是利用至少三个基准点在第一时间点和第二时间点时的像点坐标来解算确定。

步骤203，根据六自由度变化量、待测点在第一时间点和第二时间点的像点坐标确定待测点在相机平台坐标系中至少一个方向的位移量。

具体地，步骤202得到的六自由度变化量可以用来修正待测点在相机平台坐标系中至少一个方向的位移量，其中，相机平台坐标系中的方向可以包括相机成像的像平面内的两个方向和相机深度方向。而在修正时，需要利用待测点在第一时间点和第二时间点的像点坐标、以及六自由度变化量来确定待测点在相机平台坐标系中至少一个方向上的位移量。

本申请实施例的变形监测方法中，由于相机单元中相机的测量速度是由相机的帧率(采图速率)决定，而且相机的帧率可以根据实际需要进行设定，所以监测方法中的测量速度是不仅可控，而且测量速度快。因此，变形监测方法基于相机单元的拍摄图像进行变形监测，不仅可以做到对土木结构的动态变形监测，也可以减少监测成本；而基于基准点确定了精确的相机平台的六自由度变化量，再利用相机平台的六自由度变化量帮助确定待测点在相机平台坐标系中至少一个方向的位移量，可以有效提高变形监测测量结果的精度。

在一个可能的实施例中，步骤201包括：

步骤2011，获取在第一时间点和/或第二时间点拍摄的包含基准点和/或待测点的拍摄图像；

具体地，以在第一时间点获得的拍摄图像为例进行说明，从相机单元中获得的拍摄图像中，可以是仅包括基准点的图像、仅包括待测点的图像、同时包括基准点和待测点的图像。

步骤2012，对拍摄图像进行图像像点定位处理，得到基准点和/或待测点在第一时间点和/或第二时间点的像点坐标。

具体地，对拍摄图像进行图像像点定位处理，以确定图像中的基准点或待测点在图像中的像点坐标，图像像点定位处理可以采用DIC(Digital Image Correlation，数字图像相关法)技术或最小二乘匹配技术等技术，而DIC技术可以包括自适应模板相关滤波法、自适应阈值重心法、灰度图拟合法等。

在一个可能的实施例中，步骤201中，在对拍摄图像进行图像像点定位处理之前，还包括：

步骤2013，对拍摄图像进行滤波处理，以滤除拍摄图像中的干扰像素。

具体地，由于大气湍流的干扰，当一光束在大气中传过一段距离后，在垂直其传播方向的平面内光束其中心位置将作随机变化。这种发起抖动的现象将对图像中标志点的中心的跟踪提取引入较大误差。因此，本发明实施例还采用数据滤波等数据分析方法抑制大气抖动的影响，对相机采集得到的图像进行滤波，滤除标志点的图像中由于大气抖动造成的干扰像素变化，保留标志点的图像中由于发生真实空间位移产生的像素变化，以保障测量精度。效果对比如图3a和图3b所示，图3a、图3b是本发明实施例提供的一种大气抖动对变形监测的影响示意图，其中，图3a示意了未消除大气抖动影响的位移情况，图3b示意了消除大气抖动影响的位移情况，效果对比结果显示：在中午大气抖动最为显著的环境下，通过消除大气抖动的影响之后，测量精度可达到0.6mm以内，亚毫米级测量精度可以保证。但如果不消除大气抖动的影响，得到的测量结果波动大，相对位移变化量甚至可以达到23mm，导致监测结果不可用。

在一个可能的实施例中，步骤202中，根据至少三个基准点在第一时间点和第二时间点的像点坐标解算相机平台的六自由度变化量；其中，以第一时间点早于第二时间点为例，实时提取基准点在第二时间点时其在图像中的像点坐标，计算此时的基准点在图像坐标系中相对于第一时间点时的像点坐标的水平和竖向像素变化量，再通过成像几何约束关系实时解算出相机平台在第二时间点相对于第一时间点的三维姿态变化(即六自由度位移)，包括了：相机平台在相机平台坐标系中三个方向的平移和三个方向的旋转。

具体地，如果有三个以上的基准点，则根据摄像测量原理可以解算出相机平台在第二时间点相对于第一时间点的六自由度变化量。如果有多个相机用来获得三个以上的基准点时，则可以联立方程进行优化求解。

下面以有两个相机为例，第一时间点t₀时，基准点到相机的成像关系可以表示为：

其中

分别为世界坐标系W到t₀时刻相机坐标系

的投影矩阵，

为相机内参数矩阵，

为刚体变换矩阵,R,T分别为旋转矩阵和平移向量，P为齐次坐标，λ＝(MP)_z为深度因子，如果标定相机相差系数，则p为去畸变后的理想像点。

第二时间点t₁时，基准点到相机的成像关系可以表示为：

引入世界坐标系W到相机平台坐标系B(即测量基准坐标系)变换关系：

P^B＝G_W,BP^W (3)

相机平台坐标系从t₀时刻到t₁时刻的刚体运动可以表示为：

将式(3)带入式(1)有：

将式(3)、(4)带入式(2)有：

如果相机视场内有三个以上的基准点，则可以联立(5)(6)求解出

及

从而可以分解得到相机平台的六自由度变化量

特别说明的是，当有两个以上的相机来获得三个以上的基准点的图像时，只要两个以上的相机获得的图像中包括的基准点的总数大于或等于三个即可，每个相机的视场中并不需要包括三个以上的基准点。而当只有一个相机来获得三个以上的基准点的图像时，在该相机的视场中需要包括有三个以上的基准点，根据方程(1)、方程(2)、方程(5)、方程(6)的第一排即可解算出六自由度变化量。

接着，参考图1，根据相机平台的六自由度变化量的实时解算结果对待测点的位移量进行修正，消除相机平台晃动造成的测量误差，即得到准确的待测点相对于相机平台坐标系的至少一个方向上的位移量。如果在同一个相机视场内只有一个待测点，仅靠相机进行测量只能得到相机成像的像平面内两个方向的位移变化量，此时还需要引入相机深度方向的信息，根据相机深度方向的信息、六自由度变化量、待测点在第一时间点和第二时间点处的像点坐标，才可以获得待测点在相机平台坐标系中三个方向的位移量。其中，可以通过激光测距机等测距模块来提供相机深度方向的信息，即测量待测点沿相机深度方向的位移变化，也即测量待测点到相机单元(具体是指拍摄待测点的相机)的距离。如果在同一个相机视场内有两个以上的待测点，考虑到相机的测量范围相对于被测土木结构的尺寸而言较小，可以将两个以上待测点的运动建模为局部平移运动，则可省去测距模块，直接得到待测局部区域在测量基准坐标系下的三个方向的位移量。如果在同一个相机视场内有三个以上待测点，考虑到相机的测量范围相对于被测土木结构的尺寸而言较小，可以将三个以上待测点的运动建模为局部刚体运动，则可省去测距模块，直接得到待测局部区域在测量基准坐标系下的整体六自由度运动。

在一个可能的实施例中，监测系统还包括设置在相机平台上的、用于测量相机单元和待测点之间的距离的测距模块，测距模块与相机单元并排设置，且测距模块的测距方向和相机单元的拍摄方向相同；具体地，测距模块测量的是拍摄待测点的相机到待测点之间的距离。

步骤203包括：

步骤2031，相机单元的同一相机视场中有且只有一个目标待测点时，根据六自由度变化量、目标待测点在第一时间点和第二时间点的像点坐标、以及成像几何约束关系，确定目标待测点在相机平台坐标系中的第一方向和第二方向上的位移量；

步骤2032，获取对应第一时间点和第二时间点的，目标待测点和相机单元之间的距离，并根据距离和六自由度变化量确定目标待测点在相机平台坐标系中的第三方向的位移量。

具体地，利用六自由度变化量消除相机平台晃动的影响，补偿基于相机单元得到测量结果，提高测量精度。以相机单元的一个相机中的一个待测点为例，有：

其中

对于相机视场内的一个待测点，除了会有相机平台运动带来的待测点运动之外，还会有待测点的实际运动(在测量基准坐标系)，因而有：

其中ΔP＝[Δx,Δy,Δz,0]为待测点在测量基准坐标系下的实际位移量。

对于每个测量点根据式(8)可以提供关于测量基准坐标系下实际位移量的两个约束，因而，联立方程(5)(6)(8)，并可以通过测距模块辅助提供相机深度方向的信息，则可以线性求解出待测点的实际运动ΔP，即在相机平台坐标系中的三个方向的位移量。

在另一个可能的实施例中，步骤203包括：

步骤2033，相机单元的同一相机视场中有两个以上的目标待测点时，将两个以上的目标待测点的变形运动建模为局部平移运动；

步骤2034，根据六自由度变化量、目标待测点在第一时间点和第二时间点的像点坐标、以及局部平移运动对应的成像几何约束关系，确定目标待测点在相机平台坐标系中的第一方向、第二方向、第三方向上的位移量。

具体地，如果在同一个相机视场内有两个以上的待测点，考虑到测量相机的测量范围相对于被测土木结构的尺寸较小，则可以将两个以上待测点的运动建模为局部平移运动，即：

根据式(9)，利用两个或者两个以上待测点，忽略

中三个方向的转角，可以求得

中三个方向位移分量，也即待测局部区域在测量基准坐标系下的三个方向的位移量。

在又一个可能的实施例中，步骤203包括：

步骤2035，相机单元的同一相机视场中有三个以上的目标待测点时，将三个以上的目标待测点的变形运动建模为局部刚体运动；

步骤2036，根据六自由度变化量、目标待测点在第一时间点和第二时间点的像点坐标、以及局部刚体运动对应的成像几何约束关系，确定目标待测点在相机平台坐标系中的第一方向、第二方向、第三方向上的位移量。

具体地，如果在同一个相机视场内有三个以上待测点，考虑到测量相机测量范围相对于被测土木结构的尺寸较小，则可以将三个以上待测点的运动建模为局部刚体运动。根据式(9)利用三个以上待测点，可以求得待测局部区域在测量基准坐标系下的六自由度运动

特别指出的是，对于相机平台坐标系中一个方向、两个方向的修正方法与上述修正三个方向是一样的，都是利用上面的公式。上述步骤2031-步骤2036的步骤编号仅为区分不同步骤，不对步骤的执行顺序造成限定。上述对应步骤2034-步骤2036的两个实施例的方法，虽然可以获得待测点在相机平台坐标系中至少三个方向的位移量，但是，在相机深度(即相机光轴方向)这一方向上的精度较低，为了提高这一方向上的位移测量精度，可以选配激光测距机等测距模块同相机单元固联安装在相机平台上，以提高待测点沿测量相机光轴方向的变形测量精度。同样地，对于相机平台的六自由度变化量的精度，也可以利用相机深度的信息来提升对应方向上的精度。

在一个可能的实施例中，相机单元至少包括一个相机，用于获得包括基准点和/或待测点的图像。相机单元中的相机可以是同时用于拍摄基准点和待测点的相机，也可以是按照拍摄对象的不同进行划分，如划分为拍摄基准点的校正相机，以及拍摄待测点的测量相机，此时的相机单元的相机数目需要至少两个，其中，至少一个为测量相机，至少一个为校正相机。在实际监测过程中，参考图4，图4是本发明实施例提供的一种土木结构变形监测系统的结构示意图；各校正相机404和测量相机403严格固连在同一个相机平台402并同步触发拍摄，相机帧率可根据用户需求选取，可满足对大型结构的静态、准静态、动态变形进行高精度的监测需求。相机平台上固连安装的校正相机或测量相机的视场需要能覆盖到位于不同方位的基准点和/或待测点，例如，测量相机403的视场覆盖到了待测点401，而校正相机404的视场覆盖到了基准点405。

其中，监测系统中，若干校正相机观测不同位置处的若干个稳定基准点或运动已知的基准点(校正相机的数量配置可根据测量需求进行调整)，可实现对观测相机平台的六自由度姿态变化的高精度实时解算，从而消除相机观测站自身不稳定的影响，可满足对土木结构变形进行高精度(亚毫米级)且长期稳定动态监测的需求。另外，校正相机拍摄的基准点，其位置不限，可以根据测量现场环境选定。

特别地，在相机固连安装完成后，需要先对相机进行标定，且标定后，在监测工作过程中相对安装关系不发生变化。在相机标定的时候，需要足够的标志点数量，所以基准点和待测点都属于标志点，但是用于标定的标志点不仅仅包含基准点和待测点。当基准点和待测点数量不够的时候，还需要再加一些临时的标志点，这些标志点仅用来做相机标定，标定完了就可以撤离。

以结合标志点的空间坐标和像点坐标对相机进行标定的方法为例，在标定之前，需要获取用于标定的标志点的空间坐标，可以通过全站仪等手段测得，测得手段还包括GPS，或者，经纬仪配合光电测距仪，或者，三维激光扫描仪等可获取标志点空间坐标的手段。相机参数标定时，以全站仪为例，对标定某一个相机进行说明：利用全站仪依次对该相机视场中的标志点进行测量，获得各标志点的三维空间坐标；与此同时，该相机对标志点进行同步图像采集，提取得到各标志点的像点坐标；根据标志点在全站仪坐标系中的空间坐标和对应像点的图像坐标，利用光束法平差进行相机内外参数的优化求解，完成对相机的标定。

还可以利用其他计算机视觉标定方法对相机进行标定，本申请不做特别限定。

在一个可能的实施例中，本申请中的标志点既可以是被测土木结构上的自然特征，也可以为固定安装的合作标志，合作标志可以是圆形，也可以是对顶角，或十字丝或其他易于识别的形状。标志可以主动发光，也可以依靠反射日光或其他光源进行成像。标志优选为红外发光标志，以满足全天时测量需求。

在一个可能的实施例中，测量相机的个数由待测点的个数和分布来确定。例如：若待测点只有一个，那么只要一个测量相机即可拍到这个点。如果有两个待测点，并且这两个点距离比较近，可以用一个相机拍到，那么只用一个测量相机即可，但如果这两个点距离比较远，用一个相机不能同时拍到，那么就用两个测量相机。依次类推。

而校正相机的个数，是根据待测点监测的精度需求和待测点监测的方向需求确定的。例如，如果只关注待测点的一个方向的位移。那么至少只要一个校正相机即可满足，且精度也可满足。例如，如果测量需求只关注待测目标的竖向沉降，则每套测量系统可仅配置一个校正相机。

如果关注待测点的二维位移(测量相机面内两个方向位移)，则需要至少一个校正相机。若测量精度要求高，需要两个或两个以上的校正相机。

如果关注待测点的三个方向的位移，那么只要一个校正相机也可以，但是精度不够，如果对精度要求高，此时需要用两个或两个以上的校正相机。

高精度解算相机平台的六自由度变化量是测量相机自校准高精度获得待测点三个方向位移变化量的关键。通过数值仿真，对比分析了单校正相机和双校正相机解算观测平台六自由度晃动量的精度。两种模型中除校正相机个数不同外，其余条件均相同。数值仿真流程如下：首先设定相机参数、图像分辨率大小、拍摄距离，并随机给定相机姿态角，生成标志物理想图像。在理想图像点上加入提取误差来模拟实际标志定位误差，通过带有误差的图像标志定位点求解相机位姿，计算姿态角解算误差，分析误差的均值、标准差、均方根值等。

表1均方根误差对比表格

	单校正相机	双校正相机
			相机坐标系x轴方向位移ΔX(mm)	1.444	0.411
相机坐标系y轴方向位移ΔY(mm)	1.524	0.296
			相机光轴z方向位移ΔZ(mm)	61.866	0.396
绕相机坐标系x轴旋转Δα(角秒)	2.516	0.544
			绕相机坐标系y轴旋转Δβ(角秒)	2.407	0.684
绕相机坐标系z轴旋转Δγ(角秒)	46.218	0.536

数值仿真结果如图5a、图5b和表1所示，其中，图5a示意了单校正相机解算观测平台姿态的误差情况，图5b示意了双校正相机解算观测平台姿态的误差。仿真结果表明：单个校正相机能解算出相机平台的六个自由度位移量(但沿相机深度方向的位移和转角相对于其他两个方向的位移和转角精度低一到两个数量级)，而两个呈一定夹角放置的校正相机，可补充单个相机深度方向的解算结果，可高精度解算出相机平台的六自由度变化量。并且单个校正相机对滚转角的解算精度较差，解算误差较其他两个方向的旋转角高了一到两个数量级，其原因在于单个相机滚转量引起的图像变化量不敏感；而对于两个呈一定夹角的校正相机，其中一个相机的滚转角，可由另一相机的俯仰角补偿求解，从而可大大提高对观测平台滚转角的求解精度。因此，在实际工程应用中可根据测量需求调整校正相机的个数，以达到最佳配置。

在一个可能的实施例中，参考图6a、图6b、图6c、图6d，是本发明实施例提供的一种土木结构变形监测系统的结构示意图；其中，第一种可能，参考图6a，监测系统中，相机平台602上，校正相机606配合云台605进行使用，通过云台的运动观测不同位置的基准点601；而测量相机604通过在固定安装时对准控制点603而不需要用云台。第二种情况，参考图6b，测量相机604配合云台605进行使用，通过云台的运动观测不同位置的待测点603；而相反，校正相机606通过在固定安装时对准基准点601而不需要用云台。第三种情况，参考图6c，校正相机606和测量相机604分别设置一个云台进行使用，通过云台的运动，可以观测不同位置的基准点和待测点。第四种情况，相机单元中，仅由一个相机607配合云台605进行使用，通过云台的运动，仅由一个相机观测不同位置的基准点和待测点。

另外，如果某测量相机或校正相机视场内既有若干待测点也有若干不动基准点，则该相机即可作为测量相机也可作为校正相机，从而适当减少相机观测站中的相机个数。

在一个可能的实施例中，通过合理设置测量目标的图像物面分辨率，以及合理选择标志定位提取算法，可以保障亚毫米级的测量精度。图像物面分辨率是决定摄像测量精度的最重要因素。如果图像物面分辨率达到了毫米量级，即一个像素代表物面的几个平方毫米，配合亚像素的定位提取算法，则其测量精度就可能达到亚毫米量级。或者，如果图像物面分辨率达到了亚毫米量级，则配合整像素的提取算法，其测量精度就可能达到亚毫米量级；如果配合亚像素的提取算法，则测量精度可更高，以此类推。因此，设计一个具体测量任务，首先要根据测量精度的需求，对图像物面分辨率和合作标志定位提取算法进行综合考虑，选取合适的合作标志定位提取算法，并设定相应的图像物面分辨率，进而确定要达到此物面分辨率所需要的相机分辨率、镜头焦距和视场大小等。

具体地，参考图7，图7是本发明实施例提供的一种单相机摄像测量基本原理示意图；设拍摄视场大小为W×H，设相机分辨率为M×N，则图像物面分辨率的计算公式为：

R_x＝W/M，R_y＝H/N

视场是与物距相关的量，物距越远视场越大，如图，因此图像物面分辨率与物距有关，又称为景物对像面的放大系数。设感光像元尺寸为dx×dy，焦距为f，物距为D，则根据相似关系得：

本发明实施例中，通过安装若干个校正相机实现相机平台自身的六自由度姿态变化的高精度解算，同时通过安装若干个测量相机实现大尺度范围变形同步高精度测量，有效拓展了摄像测量技术的应用范围和解决实际工程问题的能力，可用于桥涵、路基、隧道、大坝、高层建筑等结构的变形监测，具有重要的理论研究意义和广泛的应用前景。基于摄像测量原理的静动态亚毫米级变形自校准摄像监测方法，可监测各类结构变形。摄像测量方法理论成熟、测量精度高、测量频率高、设备成本低、操作方便、可对任意多点进行监测、无需棱镜也无需年检，系统数字化程度高，可满足军用和民用领域中对大型结构变形进行长时间、多点、动态、非接触、实时、亚毫米级精度测量的需求，极大提高了大型结构变形测量水平。

基于上述土木结构变形监测方法实施例的描述，本发明实施例还公开了一种土木结构变形监测系统，土木结构变形监测系统包括标志点、处理装置以及设置在同一相机平台的用于获取标志点的像点坐标的相机单元，标志点包括至少三个基准点、被测土木结构的至少一个待测点；

相机单元，用于获取在第一时间点和第二时间点拍摄标志点得到拍摄图像；

参考图8，图8是本发明实施例提供的一种土木结构变形监测系统的处理装置的结构示意图；处理装置包括：

获取单元801，用于根据拍摄图像获取至少三个基准点和待测点在第一时间点和第二时间点的像点坐标；

第一确定单元802，用于根据至少三个基准点在第一时间点和第二时间点的像点坐标确定相机平台的六自由度变化量；

第二确定单元803，用于根据六自由度变化量、待测点在第一时间点和第二时间点的像点坐标确定待测点在相机平台坐标系中至少一个方向的位移量。

在一个可能的实施例中，相机单元包括至少一个相机。

在一个可能的实施例中，获取单元801包括：

第一模块，用于获取在第一时间点和/或第二时间点拍摄的包含基准点和/或待测点的拍摄图像；

第二模块，用于对拍摄图像进行图像像点定位处理，得到基准点和/或待测点在第一时间点和/或第二时间点的像点坐标。

在一个可能的实施例中，获取单元801还包括：

第三模块，用于在对拍摄图像进行图像像点定位处理之前，对拍摄图像进行滤波处理，以滤除拍摄图像中的干扰像素。

在一个可能的实施例中，监测系统还包括设置在相机平台上的、用于测量相机单元和待测点之间的距离的测距模块，测距模块与相机单元并排设置，且测距模块的测距方向和相机单元的拍摄方向相同。另外，第二确定单元803包括：

第四模块，用于相机单元的同一相机视场中有且只有一个目标待测点时，根据六自由度变化量、目标待测点在第一时间点和第二时间点的像点坐标、以及成像几何约束关系，确定目标待测点在相机平台坐标系中的第一方向和第二方向上的位移量；获取对应第一时间点和第二时间点的，目标待测点和相机单元之间的距离，并根据距离和六自由度变化量确定目标待测点在相机平台坐标系中的第三方向的位移量。

在另一个可能的实施例中，第二确定单元803包括：

第五模块，用于相机单元的同一相机视场中有两个以上的目标待测点时，将两个以上的目标待测点的变形运动建模为局部平移运动；根据六自由度变化量、目标待测点在第一时间点和第二时间点的像点坐标、以及局部平移运动对应的成像几何约束关系，确定目标待测点在相机平台坐标系中的第一方向、第二方向、第三方向上的位移量。

在又一个可能的实施例中，第二确定单元803包括：

第六模块，用于相机单元的同一相机视场中有三个以上的目标待测点时，将三个以上的目标待测点的变形运动建模为局部刚体运动；根据六自由度变化量、目标待测点在第一时间点和第二时间点的像点坐标、以及局部刚体运动对应的成像几何约束关系，确定目标待测点在相机平台坐标系中的第一方向、第二方向、第三方向上的位移量。

值得指出的是，其中，土木结构变形监测系统的具体功能实现方式可以参见上述土木结构变形监测方法的描述，这里不再进行赘述。土木结构变形监测系统中的各个单元或模块可以分别或全部合并为一个或若干个另外的单元或模块来构成，或者其中的某个(些)单元或模块还可以再拆分为功能上更小的多个单元或模块来构成，这可以实现同样的操作，而不影响本发明的实施例的技术效果的实现。上述单元或模块是基于逻辑功能划分的，在实际应用中，一个单元(或模块)的功能也可以由多个单元(或模块)来实现，或者多个单元(或模块)的功能由一个单元(或模块)实现。

基于上述方法实施例以及装置实施例的描述，本发明实施例还提供一种土木结构变形监测设备。

请参见图9，是本发明实施例提供的一种土木结构变形监测设备的结构示意图。如图9所示，上述的土木结构变形监测系统可以应用于所述土木结构变形监测设备900，所述土木结构变形监测设备900可以包括：处理器901，网络接口904和存储器905，此外，所述土木结构变形监测设备900还可以包括：用户接口903，和至少一个通信总线902。其中，通信总线902用于实现这些组件之间的连接通信。其中，用户接口903可以包括显示屏(Display)、键盘(Keyboard)，可选用户接口903还可以包括标准的有线接口、无线接口。网络接口904可选的可以包括标准的有线接口、无线接口(如WI-FI接口)。存储器905可以是高速RAM存储器，也可以是非不稳定的存储器(non-volatile memory)，例如至少一个磁盘存储器。存储器905可选的还可以是至少一个位于远离前述处理器901的存储装置。如图9所示，作为一种计算机存储介质的存储器905中可以包括操作系统、网络通信模块、用户接口模块以及设备控制应用程序。

在图9所示的土木结构变形监测设备900中，网络接口904可提供网络通讯功能；而用户接口903主要用于为用户提供输入的接口；而处理器901可以用于调用存储器905中存储的设备控制应用程序，以实现上述土木结构变形监测方法的步骤。

应当理解，本发明实施例中所描述的土木结构变形监测设备900可执行前文所述土木结构变形监测方法，也可执行前文所述土木结构变形监测系统的描述，在此不再赘述。另外，对采用相同方法的有益效果描述，也不再进行赘述。

此外，这里需要指出的是：本发明实施例还提供了一种计算机存储介质，且所述计算机存储介质中存储有前文提及的土木结构变形监测系统所执行的计算机程序，且所述计算机程序包括程序指令，当处理器执行所述程序指令时，能够执行前文所述土木结构变形监测方法的描述，因此，这里将不再进行赘述。另外，对采用相同方法的有益效果描述，也不再进行赘述。对于本发明所涉及的计算机存储介质实施例中未披露的技术细节，请参照本发明方法实施例的描述。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存储记忆体(Random Access Memory，RAM)等。

以上所揭露的仅为本发明较佳实施例而已，当然不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。

Claims

1.一种土木结构变形监测方法，其特征在于，应用于土木结构变形监测系统，所述监测系统包括标志点，以及设置在同一相机平台的用于在第一时间点和第二时间点拍摄所述标志点得到拍摄图像的相机单元，所述标志点包括至少三个基准点、被测土木结构的至少一个待测点；所述方法包括：

2.根据权利要求1所述的方法，其特征在于，所述监测系统还包括设置在所述相机平台上的、用于测量所述相机单元和所述待测点之间的距离的测距模块，所述测距模块与所述相机单元并排设置，且所述测距模块的测距方向和所述相机单元的拍摄方向相同；

3.根据权利要求1所述的方法，其特征在于，所述根据所述六自由度变化量、所述待测点在第一时间点和第二时间点的像点坐标确定所述待测点在相机平台坐标系中至少一个方向的位移量，包括：

4.根据权利要求1所述的方法，其特征在于，所述根据所述六自由度变化量、所述待测点在第一时间点和第二时间点的像点坐标确定所述待测点在相机平台坐标系中至少一个方向的位移量，包括：

5.根据权利要求1至4任一项所述的方法，其特征在于，所述根据所述拍摄图像获取所述至少三个基准点和所述待测点在第一时间点和第二时间点的像点坐标，包括：

6.根据权利要求5所述的方法，其特征在于，在对所述拍摄图像进行图像像点定位处理之前，所述方法还包括：

7.一种土木结构变形监测系统，其特征在于，包括标志点、处理装置以及设置在同一相机平台的用于获取所述标志点的像点坐标的相机单元，所述标志点包括至少三个基准点、被测土木结构的至少一个待测点；

所述处理装置包括：

8.根据权利要求7所述的系统，其特征在于，所述相机单元包括至少一个相机。

9.一种土木结构变形监测设备，其特征在于，包括：处理器和存储器；

所述处理器和存储器相连，其中，所述存储器用于存储程序代码，所述处理器用于调用所述程序代码，以执行如权利要求1-6任一项所述的土木结构变形监测方法。

10.一种计算机存储介质，其特征在于，所述计算机存储介质存储有计算机程序，所述计算机程序包括程序指令，所述程序指令当被处理器执行时，执行如权利要求1-6任一项所述的土木结构变形监测方法。