CN115100126A - 一种桥梁结构平面位移场智能感知方法 - Google Patents
一种桥梁结构平面位移场智能感知方法 Download PDFInfo
- Publication number
- CN115100126A CN115100126A CN202210663664.0A CN202210663664A CN115100126A CN 115100126 A CN115100126 A CN 115100126A CN 202210663664 A CN202210663664 A CN 202210663664A CN 115100126 A CN115100126 A CN 115100126A
- Authority
- CN
- China
- Prior art keywords
- grid
- bridge structure
- displacement
- plane
- image
- 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.)
- Pending
Links
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000003384 imaging method Methods 0.000 claims abstract description 10
- 230000009466 transformation Effects 0.000 claims abstract description 6
- 238000006243 chemical reaction Methods 0.000 claims abstract description 5
- 238000005259 measurement Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 abstract description 4
- 230000000007 visual effect Effects 0.000 description 9
- 238000000691 measurement method Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
- G06T2207/30132—Masonry; Concrete
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Quality & Reliability (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提出了一种桥梁结构平面位移场智能感知方法,属于监测领域,包括:在待测桥梁关键位置处设置靶标,利用固定在自动旋转云台上摄像机对桥梁进行扫描式成像,通过透视变换将扫描得到的图像投影至同一个平面并拼接,得到二维全景图像;对相机扰动进行消除,采用GrabCut算法提取待测桥梁结构部位的图像,进行网格离散化;利用靶标对相机进行标定以消除畸变,同时计算出像素‑位移转换系数,采用二次曲面拟合改进的模板匹配算法得出关键节点位移;采用三角形形函数和矩形形函数对关键节点位移信息进行传递,从而获取待测桥梁结构平面位移场。本方法实现了桥梁结构平面位移场的准确快速非接触获取。本发明具有便捷、准确的优点。
Description
技术领域
本发明属于桥梁工程检测与健康监测领域,具体涉及一种桥梁结构平面位移场智能感知方法。
背景技术
变形信息对于评估桥梁结构完整性和状态至关重要,它能反映结构出现的异常或缺陷。在结构检测或健康监测系统中,获取结构变形信息的主要途径是位移测量。因此,位移测量与模态识别、结构损伤识别和结构安全评估等密切相关,是不可或缺的一部分。基于计算机视觉的位移测量方法主要是利用摄像头追踪特征区域,分析其运动以实现结构位移测量。该方法可实现快速低成本的非接触式测量,近年来受到研究人员的广泛关注。
但这些当前最新的方法研究的重点在于结构单点或者多点的位移测量,对于大规模的桥梁结构来说,当前的方法只能获取该结构上的某一个点或者某几个点的位移信息,受限于摄像机的视野局限而无法同时对大规模的桥梁结构进行成像和分析。当前视觉测量方法无法对桥梁结构平面位移场进行测量。因此,亟需基于当前最新的计算机视觉技术发展一套基于视觉的桥梁结构平面位移场智能感知的方法。
发明内容
基于以上不足之处,本发明的目的是提供的一种桥梁结构平面位移场智能感知方法,能够更加直观地反映整个结构的状态,可以为后续的桥梁结构的病害识别和安全评估提供更加丰富的数据基础。
本发明所采用的技术方案如下:一种桥梁结构平面位移场智能感知方法,按以下步骤实现的:
步骤一:在待测桥梁结构关键位置处设置具备明显特征的靶标,关键位置的确定以保证后续扫描式成像得到的相邻两张图像至少包含一个靶标为准;
步骤二:利用固定在自动旋转云台上的可见光摄像机在一定距离内对待测桥梁结构进行平面扫描式成像;
步骤三:通过透视变换将扫描得到的图像投影至同一个平面并进行图像拼接,从而得到待测桥梁结构的二维全景图像;
步骤四:在图像中选取非桥梁部位的背景区域作为静止参考点,对相机的扰动进行消除,采用GrabCut算法提取待测桥梁结构部位的图像,并将其进行网格离散化;
步骤五:利用靶标对相机进行标定以消除畸变,同时计算出像素-位移转换系数,采用二次曲面拟合改进的模板匹配算法得出关键节点位移;
步骤六:采用三角形形函数和矩形形函数对关键节点位移信息进行传递,从而获取待测桥梁结构平面位移场。
进一步的,步骤二对桥梁结构平面扫描式成像过程包括:
固定在自动旋转云台上的可见光摄像机设置在待测桥梁结构的跨中位置,同待测桥梁结构平面的垂直距离根据摄像机的焦距和位移测量的精度要求共同确定;自动旋转云台的初始位置需保证可见光摄像机的镜头所在竖直面同待测桥梁结构平面平行;自动旋转云台的旋转角度需由图像拼接所需的相邻图像正向重叠率和旁向重叠率确定,旋转速度的确定不应产生运动模糊。
进一步的,步骤四对桥梁结构平面二维图像进行网格离散化的过程包括:
在经过相机扰动校正,重新裁剪之后,确定网格大小,在同尺寸全景图像上绘制网格的横线和竖线,此时的网格尚未舍弃背景部分,为全图尺寸的网格,是对全图进行规则的矩形网格离散化;对二维全景图像的待测桥梁结构部位提取图像进行二值化阈值处理,将背景置为0,桥梁结构部位置为1,形成二值化的图像,称为掩膜;将所述的掩膜的每个像素与全尺寸网格图的对应像素进行按位与运算,即1&1=1,0&1=0,从而保留桥梁结构部位范围内的网格;将桥梁结构部位范围内直线的端点首尾相连,将竖直与水平的网格线限定在其中,形成带有边界的网格,此时已经初步形成网格划分,不与边界接触的网格为规则的矩形网格;与边界接触的网格规则不一,需对其完成进一步的精细划分,边缘处的网格形状有梯形、五边形和三角形,为统一形状和精细化,将梯形和五边形分割成若干个三角形,至此完成整个位移场网格的离散划分。
进一步的,步骤六采用三角形和矩形形函数对关键节点位移信息进行传递的过程包括:
(1)已知矩形网格四个关键节点处的位移求得为pi,则网格内所有非节点处的位移p(a,b)则为:
p(a,b)=∑Yipi (1)
其中a和b分别为矩形网格的正则化坐标,Yi为节点i的形函数,其满足以下公式:
Yi=0.25v(a+m)(b+n) (2)
其中,v、m和n均为调整系数,当i为奇数时,v为1,否则为-1,当i为2或3时,m为1,否则为-1,当i大于2时,n为1,否则为-1;
(2)已知三角形网格三个关键节点处的位移求得为qi,则网格内所有非节点处的位移q(x,y)则为:
q(x,y)=∑Xiqi (3)
其中x和y分别为三角形网格的坐标,Xi为节点i的形函数,其满足以下公式:
Xi(x,y)=[(x1y2+x2y3+x3y1)-(x2y1+x3y2+x1y3)](hi+jix+kiy) (4)
其中:
h,j,k为调整系数,(x1,y1),(x2,y2),(x3,y3)分别为三角形网格第1,2,3个节点的坐标。
本发明的有益效果及优点为:本发明实现了基于计算机视觉技术的大规模桥梁结构平面位移场的智能感知;本发明便捷、准确,提升了桥梁结构位移测量的智能化水平。本发明还能满足桥梁结构施工监测。竣工验收等位移测量需求,提高了桥梁建造和管养的信息化和自动化水平。
附图说明
图1桥梁结构二维全景图像。
图2桥梁结构平面部位图像网格离散化结果。
图3桥梁结构平面位移场感知结果。
具体实施方式
本发明针对现有视觉位移测量方法只能针对单点或多点位移进行感知的局限,通过采用图像拼接的方法获取桥梁结构平面的大范围图像后,将结构部位图像分割后进行网格离散化,通过具备物理意义的形函数思想将关键节点位移信息传递到非节点处,从而实现了考虑物理限制的桥梁结构平面位移场的准确快速非接触获取。下面通过举例,通过设置在实验室的大型缩尺桥梁模型试验对本发明做进一步的说明,实验对象为三跨缩尺有机玻璃桥的悬臂跨
(2644mm)。
实施例1:
一种桥梁结构平面位移场智能感知,按照以下步骤实现的:
步骤一:在待测桥梁结构关键位置处设置具备明显特征的靶标,关键位置的确定以保证后续扫描式成像得到的相邻两张图像至少包含一个靶标为准;在表面散布着棋盘格状的靶标,并以桥梁的自由端为起点,每隔100mm进行标记;
步骤二:利用固定在自动旋转云台上的可见光摄像机在一定距离内对待测桥梁结构进行扫描式成像;在桥梁半跨的3.5m垂直距离处架设云台,为尽可能多地获取腹板细节,将相机焦距设置为最大值50mm。由于桥梁高度较小,设定云台的旋转角度幅度为2°,通过仅在桥梁跨度方向左右旋转云台的方式,即可实现对20张全桥局部图像的采集。本步骤针对当前视觉位移测量方法视野有限且多相机同时拍摄同步性差的问题,提出将自动旋转云台应用在视觉位移测量中结合后续的图像拼接技术以获取更大的视野;
步骤三:通过透视变换将扫描得到的图像投影至同一个平面并进行图像拼接,从而得到待测桥梁结构的二维全景图像;结合每张图像拍摄时的倾斜角度,通过透视变换,统一各局部图像的投影平面;最后再基于SIFT特征点实现桥梁全景拼接图;在拼接图中,散布着的棋盘格图案未发生变形,由此得到的转换系数为0.073734mm/pixel;27个100mm间距标记点的视觉测量距离的相对误差在1%以内;两处安装传感器的测点的视觉测量位移误差在2%以内,说明图像的拼接效果良好;
步骤四:在图像中选取非桥梁部位的背景区域作为静止参考点,对相机的扰动进行消除,采用GrabCut算法提取待测桥梁结构部位的图像,并将其进行网格离散化;本步骤针对当前视觉位移测量研究在结构平面位移场方面的技术空白,提出对待测结构平面进行二维全景成像后,使用网格离散化手段对位移场进行描述和求解,既加速了计算过程,又可以根据精度需求自主调节网格尺寸;
步骤五:利用靶标对相机进行标定以消除畸变,同时计算出像素-位移转换系数,采用二次曲面拟合改进的模板匹配算法得出关键节点位移;通过张氏标定法,得到相机内参,完成相机标定,标定误差仅为0.178像素,据此完成畸变校正;利用Harris角点检测法,根据角点的实际距离得到不同拍摄距离对应的转换系数;采用零均值归一化互相关系数,利用模板匹配得到整像素位移;再基于离散的相关系数矩阵,利用二次曲面拟合亚像素算法得到亚像素位移,测量误差在1个像素以内;
步骤六:采用三角形形函数和矩形形函数对关键节点位移信息进行传递,从而获取待测桥梁结构平面位移场;
(1)已知矩形网格四个关键节点处的位移求得为pi,则网格内所有非节点处的位移p(a,b)则为:
p(a,b)=∑Yipi (1)
其中a和b分别为矩形网格的正则化坐标,Yi为节点i的形函数,其满足以下公式:
Yi=0.25v(a+m)(b+n) (2)
其中,v、m和n均为调整系数,当i为奇数时,v为1,否则为-1,当i为2或3时,m为1,否则为-1,当i大于2时,n为1,否则为-1;
(2)已知三角形网格三个关键节点处的位移求得为qi,则网格内所有非节点处的位移q(x,y)则为:
q(x,y)=∑Xiqi (3)
其中(x,y)分别为三角形网格的坐标,Xi为节点i的形函数,其满足以下公式:
Xi(x,y)=[(x1y2+x2y3+x3y1)-(x2y1+x3y2+x1y3)](hi+jix+kiy) (4)
其中:
h,j,k为调整系数,(x1,y1),(x2,y2),(x3,y3)分别为三角形网格第1,2,3个节点的坐标;
本步骤针对桥梁结构平面全景图像离散化后通过上述过程只能求解节点处位移而非节点处位移无法快速求解的问题,创造性地借鉴了有限单元法中的形函数思想,将三角形形函数和矩形形函数对节点-非节点信息传递进行建模,从而实现非节点位移符合物理限制的快速求解。
实施例2
为验证本发明所提算法的准确性,利用Abaqus对实验所用的桥梁模型建立有限元模型,提取竖直方向位移场的数据,同本发明计算得到的位移场数据进行对比,使用归一化相关系数(公式6)对两个位移场的相似度进行描述。
最终求得归一化相关系数达到0.9995,验证了本发明所提算法的准确性。
Claims (4)
1.一种桥梁结构平面位移场智能感知方法,其特征在于,方法按以下步骤实现的:
步骤一:在待测桥梁结构关键位置处设置具备明显特征的靶标,关键位置的确定以保证后续扫描式成像得到的相邻两张图像至少包含一个靶标为准;
步骤二:利用固定在自动旋转云台上的可见光摄像机在一定距离内对待测桥梁结构进行平面扫描式成像;
步骤三:通过透视变换将扫描得到的图像投影至同一个平面并进行图像拼接,从而得到待测桥梁结构的二维全景图像;
步骤四:在图像中选取非桥梁部位的背景区域作为静止参考点,对相机的扰动进行消除,采用GrabCut算法提取待测桥梁结构部位的图像,并将其进行网格离散化;
步骤五:利用靶标对相机进行标定以消除畸变,同时计算出像素-位移转换系数,采用二次曲面拟合改进的模板匹配算法得出关键节点位移;
步骤六:采用三角形形函数和矩形形函数对关键节点位移信息进行传递,从而获取待测桥梁结构平面位移场。
2.根据权利要求1所述的一种桥梁结构平面位移场智能感知方法,其特征在于:步骤二对桥梁结构平面扫描式成像过程包括:
固定在自动旋转云台上的可见光摄像机设置在待测桥梁结构的跨中位置,同待测桥梁结构平面的垂直距离根据摄像机的焦距和位移测量的精度要求共同确定;
自动旋转云台的初始位置需保证可见光摄像机的镜头所在竖直面同待测桥梁结构平面平行;
自动旋转云台的旋转角度需由图像拼接所需的相邻图像正向重叠率和旁向重叠率确定,旋转速度的确定不应产生运动模糊。
3.根据权利要求2所述的一种桥梁结构平面位移场智能感知方法,其特征在于:步骤四对桥梁结构平面二维图像进行网格离散化的过程包括:
在经过相机扰动校正,重新裁剪之后,确定网格大小,在同尺寸全景图像上绘制网格的横线和竖线,此时的网格尚未舍弃背景部分,为全图尺寸的网格,是对全图进行规则的矩形网格离散化;对二维全景图像的待测桥梁结构部位提取图像进行二值化阈值处理,将背景置为0,桥梁结构部位置为1,形成二值化的图像,称为掩膜;将所述的掩膜的每个像素与全尺寸网格图的对应像素进行按位与运算,即1&1=1,0&1=0,从而保留桥梁结构部位范围内的网格;将桥梁结构部位范围内直线的端点首尾相连,将竖直与水平的网格线限定在其中,形成带有边界的网格,此时已经初步形成网格划分,不与边界接触的网格为规则的矩形网格;与边界接触的网格规则不一,需对其完成进一步的精细划分,边缘处的网格形状有梯形、五边形和三角形,为统一形状和精细化,将梯形和五边形分割成若干个三角形,至此完成整个位移场网格的离散划分。
4.根据权利要求3所述的一种桥梁结构平面位移场智能感知方法,其特征在于:步骤六采用三角形和矩形形函数对关键节点位移信息进行传递的过程包括:
(1)已知矩形网格四个关键节点处的位移求得为pi,则网格内所有非节点处的位移p(a,b)则为:
p(a,b)=∑Yipi (1)
其中a和b分别为矩形网格的正则化坐标,Yi为节点i的形函数,其满足以下公式:
Yi=0.25v(a+m)(b+n) (2)
其中,v、m和n均为调整系数,当i为奇数时,v为1,否则为-1,当i为2或3时,m为1,否则为-1,当i大于2时,n为1,否则为-1;
(2)已知三角形网格三个关键节点处的位移求得为qi,则网格内所有非节点处的位移q(x,y)则为:
q(x,y)=ΣXiqi (3)
其中(x,y)为三角形网格的坐标,Xi为节点i的形函数,其满足以下公式:
Xi(x,y)=[(x1y2+x2y3+x3y1)-(x2y1+x3y2+x1y3)](hi+jix+kiy) (4)
其中:
h,j,k为调整系数,(x1,y1),(x2,y2),(x3,y3)分别为三角形网格第1,2,3个节点的坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210663664.0A CN115100126A (zh) | 2022-06-10 | 2022-06-10 | 一种桥梁结构平面位移场智能感知方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210663664.0A CN115100126A (zh) | 2022-06-10 | 2022-06-10 | 一种桥梁结构平面位移场智能感知方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115100126A true CN115100126A (zh) | 2022-09-23 |
Family
ID=83290154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210663664.0A Pending CN115100126A (zh) | 2022-06-10 | 2022-06-10 | 一种桥梁结构平面位移场智能感知方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115100126A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117528065A (zh) * | 2023-12-28 | 2024-02-06 | 深圳市城市公共安全技术研究院有限公司 | 相机扰动效果评估及消除方法、装置、设备及存储介质 |
-
2022
- 2022-06-10 CN CN202210663664.0A patent/CN115100126A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117528065A (zh) * | 2023-12-28 | 2024-02-06 | 深圳市城市公共安全技术研究院有限公司 | 相机扰动效果评估及消除方法、装置、设备及存储介质 |
CN117528065B (zh) * | 2023-12-28 | 2024-04-26 | 深圳市城市公共安全技术研究院有限公司 | 相机扰动效果评估及消除方法、装置、设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102376089B (zh) | 一种标靶校正方法及系统 | |
CN109658457B (zh) | 一种激光与相机任意相对位姿关系的标定方法 | |
JP3735344B2 (ja) | キャリブレーション装置、キャリブレーション方法、及びキャリブレーション用プログラム | |
JP5671281B2 (ja) | 位置姿勢計測装置、位置姿勢計測装置の制御方法及びプログラム | |
US20180075618A1 (en) | Measurement system and method for measuring multi-dimensions | |
WO2021004416A1 (zh) | 一种基于视觉信标建立信标地图的方法、装置 | |
CN102788572B (zh) | 一种工程机械吊钩姿态的测量方法、装置及系统 | |
CN105701827A (zh) | 可见光相机与红外相机的参数联合标定方法及装置 | |
CN109341668B (zh) | 基于折射投影模型和光束追踪法的多相机测量方法 | |
CN107392849B (zh) | 基于图像细分的靶标识别与定位方法 | |
Ding et al. | Crack detection and quantification for concrete structures using UAV and transformer | |
CN110223355B (zh) | 一种基于双重极线约束的特征标志点匹配方法 | |
CN107358628B (zh) | 基于靶标的线阵图像处理方法 | |
CN103743352A (zh) | 一种基于多相机匹配的三维变形测量方法 | |
CN110146030A (zh) | 基于棋盘格标志法的边坡表面变形监测系统和方法 | |
JP2011198330A (ja) | 3次元位置合わせにおける照合方法およびそのプログラム | |
CN110763204A (zh) | 一种平面编码靶标及其位姿测量方法 | |
CN109920009B (zh) | 基于二维码标识的控制点检测与管理方法及装置 | |
CN101487702A (zh) | 基于双目视觉的交通事故现场摄影测量方法 | |
CN111996883B (zh) | 一种检测公路路面宽度的方法 | |
Mi et al. | A vision-based displacement measurement system for foundation pit | |
CN106500729A (zh) | 一种无需控制信息的智能手机自检校方法 | |
CN115100126A (zh) | 一种桥梁结构平面位移场智能感知方法 | |
CN112508885B (zh) | 一种弯管的三维中轴线检测方法及系统 | |
CN114897947A (zh) | 一种基于时空间统一的热红外和可见光图像同步配准方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |