基于标志点的多相机数字图像相关方法精度分析及土木工程中的应用

第 51卷第 2期 2021年 3月
东 南 大 学 学 报 (自 然 科 学 版 )
JOURNALOFSOUTHEASTUNIVERSITY (NaturalScienceEdition)
DOI:10.3969/j.issn.1001-0505.2021.02.005
Vol.51 No.2 Mar.2021
基于标志点的多相机数字图像相关方法 精度分析及土木工程中的应用
魏 康  员 方  董志强  邵新星  吴 刚  何小元
(东南大学土木工程学院,南京 210096)
摘要:为实现土木工程构件的高精度全场变形测量,分析了基于标志点的多相机测量方法的实 际精度,并将该方法用于 BFRP筋海水海砂珊瑚 UHPC组合梁的抗弯实验.利用编码标志点对 多套双目系统的结果进行坐标统一,并通过引入拼接误差 un对多相机测量结果的坐标统一精度 进行定量评估分析
.实验结果表明:面内位移量较小时,拼接误差 un基本等同于双目系统自身的 匹配误差水平,随着位移量增大,un随之增大;拼接误差的标准偏差 σu始终在 0.01mm内波动, 且位移量最大时拼接误差的均值u-与实际位移的相对误差仍小于 0.03%,多套测量系统的坐标 系统一精度满足工程需求;对于混凝土组合梁的全场位移和应变场的测量,局部位移数据与位移 计数据基本吻合.多相机数字图像相关方法可以大幅提高变形测量分辨率,适合在大型土木工程 构件变形测量中推广运用. 关键词:数字图像相关;非重叠视场;标志点;混凝土组合梁 中图分类号:O34  文献标志码:A  文章编号:1001-0505(2021)02021908
Precisionanalysisandcivilapplicationofmulticameradigitalimage correlationmethodbasedonencodedtargets
WeiKang  YuanFang  DongZhiqiang  ShaoXinxing  WuGang  HeXiaoyuan
(SchoolofCivilEngineering,SoutheastUniversity,Nanjing210096,China)
Abstract:Toachievehighaccuracy,fullfielddefo
rmationmeasurementofcivilengineeringstruc tures,theaccuracyofthetargetsbasedmulticamerameasurementmethodwasanalyzed,andthe methodwasappliedtothebendingtestofseawaterseasandcoralultrahighperformanceconcrete(UH PC)compositebeamsreinforcedwithbasaltfiberreinforcedpolymer(BFRP)bars.Codingtargets wereusedtounifythecoordinatesoftheresultsofmultiplebinocularsystems,andthestitchingerror unwasintroducedtoquantitativelyevaluateandanalyzethecoordinateunificationaccuracyofthe measurementresultsofmultiplecameras.Experimentalresultsshow thatwhentheinplanedisplace mentissmaller,thestitchingerrorunisbasicallyequaltothematchingerrorlevelofthebinocularsys tem,andunincreaseswiththeincreaseofthedisplacement.Thestandarddeviationofthestitchinger rorσualwaysfluctuateswit
hin0.01mm,andtherelativeerrorbetweenthemeanofthestitchingerror u- andtheactualdisplacementisalwayslessthan0.03% whenthedisplacementismaximum,indicating thattheunifiedprecisionofthecoordinatesysteminmultiplemeasurementsystemsmeetstheengineering requirements.Forthemeasurementofthewholefielddisplacementandthestrainfieldoftheconcrete beam,thedisplacementmeasurementdataarebasicallyconsistentwiththedisplacementmeterdata.The multicameradigitalimagecorrelationmethodcanimprovetheresolutionofthedeformationmeasure ment,thuspromotingthedeformationmeasurementoflargecivilengineeringcomponents. Keywords:digitalimagecorrelation(DIC);nonoverlappingfieldofview;encodedtargets;con
cretecompositebeams
收稿日期:20200811. 作者简介:魏康(1995—),男,硕士生;邵新星(联系人),男,博士,讲师,xinxing.shao@seu.edu.cn. 基金项目:国家自然科学基金青年基金资助项目(11902074)、国家自然科学基金重大科研仪器研制资助项目(11827801). 引用本文:魏康,员方,董志强,等.基于标志点的多 相 机 数 字 图 像 相 关 方 法 精 度 分 析 及 土 木 工 程 中 的 应 用 [J].东 南 大 学 学 报 (自 然 科 学
版),2021,51(2):219 226.DOI:10.3969/j.issn.1001-0505.2021.02.005.
220
东南大学学报(自然科学版)             第 51卷
  数字图像相关(digitalimagecorrelation,DIC) 方法由 Peters等[1]在 1982年最早提出,该方法因 其非接触、全场测量、适用范围广等特点在科研和 工程领域得到广泛关注[24].经过近几十年硬件的 快速发展,以 及 匹 配 算 法[57]、标 定 方 法[89]、测 量 精度[1011]等方面的大量研究,目前三维 DIC方法 已被广泛应用于土木工程、航空航天、生物医学等 领域的表面三维形貌和变形测量.但在土木工程领 域,被测对象通常具有结构尺寸大、柱式构件表面 曲率大的特点,此时由常规双相机组成的三维 DIC 测量系统已不
能满足实际测量需求.
为此,学者们在双相机测量系统的基础上提出 多相机三维 DIC测量方法,其中多相机指的是 3个 或更多的相机[12].实现多相机测量的关键是相机外 参的标定.相机外参的标定方式决定了多相机测量 方法应用过程的简便性(如基于重叠视场实现外参 标定的方法,需要相机两两之间连续标定,相机数量 较多时标定步骤繁琐费时)和实际场景应用中的简 便性(如基于经纬仪等大量程测量设备进行外参标 定的方法,需要相机之间利用设备互相观测,无法用 于狭小或有遮挡的测量环境),而且相机外参的标定 精度会直接影响测量结果的准确性.针对多相机外 参标定,学者们做了许多研究工作[1224],其中相机 内参的标定可在实验室内提前完成.
多相机三维 DIC测量方法的外参标定大致可 以分为 2类.第 1类是基于重叠视场[1215]外参标 定方法,该方法借助重叠视场完成多相机系统测量 数据的拼接,其实现需要各相机两两之间存在 1/3 ~1/2左右的重叠区域用以标定和计算.该方法除 了布置方式 受 限 (要 求 各 相 机 视 场 必 须 连 续 )、相 机测量视场不能得到充分利用之外,相机之间的测 量结果互相关联,测量系统稳定性差.第 2类是无 需重叠视场实现外参标定的多相机方法,其主要可 以归为以下几种类型:利用经纬仪或激光跟踪仪等 大量程测量设备 实 [1619] 现全局标定的方法;通过 追踪目标的运动模型[2022]来建立不同相机视场间 对应 关 系 的 方 法;利 用 大 规 模 标 定 物[2324]进 行 标 定的方法;基于标志点重构[2526]的外参统一方法. 其中,基于标志点重构的外参统
一的多相机方法具 有操作简便、相机可灵活布置、各子相机系统之间 基本互相独立、测量结果稳定性高的优点.
因此,本文首先设计基于标志点的多相机三维 DIC方法的精度验证实验,对该方法的可行性和精 度进行了验证.通过寻不同测量系统在重叠区域 下测量点的对应关系,利用对应点拼接误差 un 的
http://journal.seu.edu.cn
均值和标准差来评估多相机数据统一的精度.平移 验证实验表明,经坐标统一后的多套测量结果对应 良好.最后,将该方法用于玄武岩纤维增强树脂基 复合 (basaltfiberreinforcedpolymer,BFRP)筋海 水海砂珊 超高性能混凝土(ultrahighperformance concrete,UHCP)组 合 梁 的 抗 弯 承 载 实 验,得 到 破 坏临界荷载下的全场位移和应变云图,并与位移计 数据进行对比分析,结果显示本文使用的基于标志 点的多相机测量方法操作简便、布置灵活,能精确 地测量混凝土梁的表面位移.
1  基 于 标 志 点 重 构 的 多 相 机 三 维 DIC系统测量原理
  为实现不同坐标系下多相机系统的数据统一, 使用多套独立的三维 DIC测量子系统和一台单反 相
机构成一套多相机三维 DIC测量系统,其中三 维 DIC测量子系统的相机内、外参由标定板获得, 而测量子系统之间的相对外参则借助标志点获得. 由于基于标志点的多相机测量方法无需借助重叠 区域,因此可以根据测量需求和环境限制对相机系 统进行灵活布置,工程适用场景广泛. 1.1 三维 DIC测量原理
三维 DIC测量系统主要由 2个工业相机和采 集设备组成,利用双目立体视觉原理和图像相关匹 配技术实现位移和应变的计算.首先使用 2台相机 从不同方位同时对一目标进行拍摄,将拍摄到的二 维图像进行特征匹配处理,从而寻到各点的对应 关系,再根据标定得到的两相机参数,即可利用立 体视觉原理计算出每个点在定义坐标系下的三维 坐标.双目视觉原理如图 1所示,对于世界坐标系 Owxwywzw 下任意待测点 P,该点分别成像于左、右 相机 A、B像 平 面 坐 标 系 Oi1xi1yi1和 Oi2xi2yi2的 P1、P2处,由于 P、P1、P2 三点共面,可计算出点 P 的三维坐标(X,Y,Z)w,其中 P1、P2 坐标可由标定 和图像匹配计算得到,而 Oc1、Oc2分别为左、右两相
图 1 双目立体视觉原理
第 2期
魏康,等:基于标志点的多相机数字图像相关方法精度分析及土木工程中的应用
221
机坐标系 Oc1xc1yc1和 Oc2xc2yc2的光心.图像相关 匹配技术则包括同一时刻 2幅图的散斑之间和同 一相机图片中散斑变形前后的匹配,即立体匹配和 时序匹配,利用形函数和相关函数寻同一点在不 同图像下的对应关系. 1.2 单相机近景摄影测量
对于标志点三维坐标的重构主要基于近景摄 影测量技术 来 [2728] 实现,即利用不同位置不同角 度拍摄得到的多张贴有编码标志点的被测物体图 片,使用数字图像处理和摄影测量技术处理,进而 精确测定目标在三维空间中位置的方法.由便携式 单反相机拍摄得到的一系列图片用于重构标志点 的三维坐标,借此可以标定各子系统的相对外参, 从而得到坐标转换后相同坐标系下完整的全场三 维数据.本文中使用的圆环形编码标志点编码原理 如图 2所示,编码采用顺时针的 15位二进制数,黑 为 0,白为 1,可以任意位置作为起点读数,选 择对应最小十进制数值作为编码值以保证编码的 唯一性.本文使用长度已知、热膨胀系数较小的高 精度碳纤维刚性标尺进行尺度标定.
摄影测量技术中,利用自标定光束平差法以每 个光 束 作 为 基 本 平 差 单 元 ,将 标 志 点 的 图 像 坐 标
图 2 同心圆环形编码标志点原理图
u~i,j作为测量值,按照共线条件方程列出误差方程 如下:
∑ min [u~i,j-u~(K,D,Rj,Tj,X~i)]2 (1) i,j
式中,下标 i、j分别表示标志点的编号和不同姿态
的图片序号;u(K,D,Rj,Tj,Xi)为标志点的投影图
像坐标,K,D,Rj,Tj,Xi分 别 为 相 机 内 参、镜 头 畸
变、外参以及标志点三维空间坐标.在全区域内统
一进行平差处理,对相机内参 K、镜头畸变 D、不同
姿态下相机坐标系相对于世界坐标系的外参 Rj、
Tj以及空间点坐标 Xi进行优化求解.
1.3 多相机数据统一
由三维 DIC测量子系统可以得到标志点在各
自坐标系 Ocmxcmycmzcm(m=1,2,…,s;s≥2)下的 三维坐标,将世界坐标系定在单反相机坐标系并计
算全部标志点的三维坐标,通过编码值相同的对应
关系可计算不同坐标系的刚性转置信息,实现测量
数据统一,原理如图 3所示.
该过程实质为一个最小二乘求优问题,用公式
可描述为
∑ (Rcmw,Tcmw)=argmin wt(Rcmwpt+Tcmw)-q2t t=1
(2)
式中,m为 DIC子系统序号;t、T分别表示第 m个
DIC子系统识得的标志点序号和总数;wt>0为点 集中各点对的权重;pt和 qt分别为不同坐标系下 的对应点;Rcmw,Tcmw为子系统坐标系转到世界坐 标系的旋转与平移矩阵.
图 3 多相机系统坐标统一原理图
http://journal.seu.edu.cn
222
东南大学学报(自然科学版)             第 51卷
2 精度验证
2.1 评估原理 为实现多 相 机 三 维 DIC系 统 测 量 结 果 的 定
量 评 估 和 误 差 分 析 ,建 立 重 叠 区 域 内 三 维 形 貌 数 据 的 对 应 关 系 ,以 对 应 数 据 差 值 绝 对 值 的 均 值 和 标准差来评估误 差.以 2套 三 维 DIC系 统 为 例, 对应 关 系 的 建 立 具 体
如 下:① 将 点 集 1中 在 (x,y)i1处 的 测 量 值 记 为 D(1)(x,y)i1,利 用 相 机 内、外 参 转 换 到 点 集 2的 像 素 坐 标 系 下 变 为 D(1)(x′,y′)i2,搜索 该 点 在 点 集 2中 距 离 最 近 的 整像素点 D(2)(X,Y)i2,完成初始对应;②为保证 对应精度,将点集 1中实现初始对应的点取 k×k (k可取 3、4或 5)的非线性拟合窗口,得到 测 量 值 D(1) (x′,y′)i2 在 点 集 2 的 亚 像 素 对 应 D(2)(x′,y′)i2.建 立 好 精 确 的 对 应 关 系 后 ,对 点 集 1和点集 2的测量值进行误差分析.
用于评估数据坐标统一精度的均值和标准偏 差公式分别为
u- = n=1
D(1)(x′n,y′n)-D(2)(x′n,y′n) N
(3)
∑ (un -u-)2
槡 σu =
n=1
N-1
(4)
其中
un= D(1)(x′n,y′n)-D(2)(x′n,y′n)
(5)
式中,N为对应的点总数;un为第 n个建立对应关
系的测量值(三维形貌或位移或应变数据)差值的
绝对值(简称拼接误差);u-为 N个对应点的拼接误
差的均值;D(1)(x′n,y′n)表示点集 1在(x′n,y′n)处的
测量值;D(2)(x′n,y′n)表示点集 2由(x′n,y′n)周围测
量值拟合后得到的值,可由下式得到:
D(2)(x′n,y′n)=a0 +a1x′n +a2y′n +a3 (x′n)2 +
  a4x′ny′n+a5(y′n)2
(6)
式中,a0,a1,…,a5 为拟合参数. 2.2 验证实验
为证明本 文 多 相 机 三 维 DIC测 量 方 法 的 可
行性和有 效 性,进 行 散 斑 平 面 平 移 验 证 实 验,实
验方案布置图如图 4所示.使用 4个相机组成 2
套三维 DIC测量系统,为了避免测量结果引入其
他因素的 影 响,采 用 相 同 相 机 型 号 和 镜 头 焦 距,
实验初始保 证 4个 相 机 的 视 场 中 心 均 对 准 待 测
散斑平面 中 心,并 将 2套 测 量 系 统 对 称 布 置.其
http://journal.seu.edu.cn
中,待测 平 面 的 散 斑 区 域 尺 寸 为 300mm ×300 mm,调节镜 头 焦 距 和 光 圈 至 被 测 平 面 表 面 清 晰 成 像 ,清 晰 成 像 时 图 像 中 每 像 素 对 应 的 空 间 尺 寸 约为 0.23mm,而粘贴在待测平面的标志点内圆 心直径为 5mm.
(a)相机系统布置图
(b)实验装置布置图
图 4 平移测量验证实验示意图
平移实验开始前,对静止的待测平面进行一段 时间的噪声采集,用以评估 2套 DIC测量系统的 整体误差水平.由于 2套 DIC测量系统对称布置, 在此以相机系统 A为例进行说明,计算区域取全
部散斑区域,具 体 系 统 误 差 (静 态 噪 声 水 平 )结 果 如图 5和表 1所示.可以看到:该测量系统计算的 面内位移 u和 v均值在 ±0.005像素内波动,离面 位移 w均值在 ±0.02像素内,表面应变 εxx、εxy和 εyy均值范围为 ±25×10-6.因此可以认为待测物 体静止不动时,该测量系统的面内位移误差为0.01 像素(约为 0.0023mm).
平移实验具体步骤:使用电动平移台对散斑平面 给予水平位移并用相机采集图像,每次平移 10mm,
图 5 相机系统 A测量结果的均值
第 2期
魏康,等:基于标志点的多相机数字图像相关方法精度分析及土木工程中的应用
223
统计误差 均值  标准差
u/像素 -0.0006
0.0027
表 1 相机系统 A测量噪声
v/像素 -0.0012
0.0018
w/像素 0.0052 0.0061
εxx/10-6 2.21 22.11
εxy/10-6 0.60 21.59
εyy/10-6 1.93
15.36
一共平移 7次,如此重复 5次.图 6(a)为统计重复 定向位移 的 计 算 结 果 数 据,可 以 发 现 当 位 移 量 为 10mm时,拼接误差的均值u-为 0.0025mm,对应 到像素单位即 0.0011像素.位移量不大时,拼接 误差的均值 u-近似等同双目系统本身的匹配误差 水平 0.01像素,说明 2套相机的面内位移测量结 果经坐标统一后仍对应良好.随着位移量的增大, 拼 接 误 差
的 均 值 u- 呈 线 性 增 长,均 值 u- 从 0.0025mm增加到 0.01mm(对应到像素单位为 0.04像素),而标准偏差 σu从 0.0094mm下降到 最小值 0.0065mm再逐步增大到 0.0084mm.这 是因为随着位移量的增大,系统本身的误差水平增 大,u-也随之增大.
图 6(b)给出了电动平移台移动 70mm 时拼 接误差 un的全场分布情况,图中共有 5966组点 对参与计算.其中,超过 95% 的点对拼接误差 un 在 0.02mm以内,且所有点对的拼接误差最大值 不超过 0.035mm,此时均值 u- 为 0.01mm,标准偏
差 σu为 0.0084mm.这说明随着位移量的增大, 拼接误差 un对应增大,但即使在平移 70mm 时, 拼接误差的均值 u-与实际位移的相对误差仍小于 0.03%,两套相机系统的位移计算结果始终对应良 好.该实验证明了利用标志点实现多相机数字图像 相关测量方法的可行性和精确性.
3 BFRP筋海水海砂珊瑚 UHPC组 合梁的抗弯实验
  为验证基于标志点多相机测量方法的有效性, 本文以 BFRP筋海水海砂珊瑚 UHPC组合梁[29] 为研究对象,测量混凝土梁在加载过程的全场位移 和全场应变情况,并与位移计结果进行对比分析. 3.1 试样准备
混凝土梁试件尺寸均为 2400mm(长)×120 mm(宽)×240mm(高),在净跨 l0 为 2.1m 和剪 切跨度为 0.75m的条件下进行四点弯加载实验, 具体实验加载情况如图 7(a)所示.其中采用的钢 筋均为 HPB235级,混凝土梁底部受拉筋为 2根直 径 13mm 的 BFRP筋,顶部使用 2根直径 10mm 的 BFRP筋和 30mm 厚的 UHPC替换层,而所有 箍筋均为 8mm的 BFRP筋,图 7(b)中给出了混凝 土梁截面的示意图.
在梁底部和两支座处共布置了 9个线性差动 变压 器 式 位 移 计 (linearvariabledifferentialtrans former,LVDT)用于测量梁的竖向位移,其中 7个
(a)不同平移距离时 un的统计误差
(a)实验加载装置示意图
(b)平移 70mm时 un的分布云图
图 6 两子系统测得面内位移结果的对应情况
(b)测量传感器布置图(单位:mm)
图 7 实验布置示意图
http://journal.seu.edu.cn

本文发布于:2024-09-20 11:50:17,感谢您对本站的认可!

本文链接:https://www.17tex.com/tex/3/88067.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:测量   位移   系统   标志点   进行   方法   误差   标定
留言与评论(共有 0 条评论)
   
验证码:
Copyright ©2019-2024 Comsenz Inc.Powered by © 易纺专利技术学习网 豫ICP备2022007602号 豫公网安备41160202000603 站长QQ:729038198 关于我们 投诉建议