技术领域
[0001] 本
发明涉及岩石裂缝扩展实验的非
接触式高速视频测量方案，尤其是涉及一种基于高速视频测量的岩石裂缝扩展实验监测方法。
背景技术
[0002] 在工程界，一种材料在投入使用之前，其材料性质和安全系数往往需要通过拉伸、压缩、碰撞等测试性实验进行监测。在岩石的裂缝扩展测量方面，传统接触式测量仪器由于量程有限、测量点位少、增加模型
质量、安装费时费
力等诸多缺点，已渐渐被非接触式测量方法所取代。在光学测量研究工作中，数字相关
算法已成为力学计算的主流解算方法。许多学者已开始运用数字散斑相关技术进行位移场和应变场的解算。但在大部分的实验中，更多的是测量二维平面的位移变化，其中感光平面与被测物面不平行导致的测量误差会严重影像测量结果，而这几乎是无法控制的。此外，三维数字相关技术在近些年也得到了广泛研究，但诸多类似的实验是在普通相机或者较低性能的高速相机下进行拍摄的，无法详尽地记录被测物的突变情况，进而无法提供测量点位完整的空间三维信息变化，影响了实验的测量效果。
发明内容
[0003] 本发明的目的就是为了克服上述
现有技术存在的
缺陷而提供一种基于高速视频测量的岩石裂缝扩展实验监测方法。
[0004] 本发明的目的可以通过以下技术方案来实现：
[0005] 一种基于高速视频测量的岩石裂缝扩展实验监测方法，包括以下步骤：
[0006] 1)通过
双目视觉测量系统中的两台高速相机获取待测岩石在单轴压缩下的破裂过程的序列影像并存储；
[0007] 2)在序列影像解析过程中，通过亚
像素级
跟踪匹配方法计算散斑兴趣区中的序列同名像点坐标，并经过摄影测量解析算法获取待测岩石表面在单轴压缩下的时序三维点
云数据；
[0008] 3)针对序列三维点云的
时空序列分析，获取待测岩石裂缝的三维形变参数，包括位移、速度、
加速度和应变场。
[0009] 所述的步骤1)中，所述的双目视觉测量系统中的高速相机为工业相机，通过同步
控制器保持影像采集同步。
[0010] 所述的两台高速相机
水平放置，并且采用交向摄影方式拍摄，以此到达增加影像重叠
覆盖率的目的。
[0011] 所述的步骤1)中，在待测岩石上
喷涂散斑区域，具体为：
[0012] 待测岩石表面打磨至平整，在观测表面喷涂白色的哑光漆，并进行
风干，在观测表面随机并均匀地喷洒黑色哑光漆或黑色墨水，形成散斑，将整个散斑区域作为散斑兴趣区，将散斑作为目标点。
[0013] 所述的步骤2)具体包括以下步骤：
[0014] 21)高速相机的立体标定：通过基于平面板的标定方法同时获取高速相机的内方位元素和外方位元素；
[0015] 22)散斑影像预处理：在序列影像中选取散斑兴趣区并在兴趣区中确定目标点；
[0016] 23)同名点匹配：将两个高速相机同时采集的影像作为匹
配对象，通过归一化相关系数确定整像素级的粗略点位，并通过最小二乘匹配方法确定亚像素级的精确点位；
[0017] 24)目标跟踪匹配：以高速相机采集的序列影像的前后
帧影像作为匹配对象，获取目标点在序列影像中的时序二维像点坐标；
[0018] 25)三维点云重建：根据匹配后的每对同名点的像点坐标以及标定后的内方位元素和外方位元素，通过基于共线方程的前方交会获取序列影像中同名目标点的时序三维点位坐标。
[0019] 所述的步骤23)中，最小二乘匹配方法以归一化相关系数最大为目标函数，考虑左右影像间的仿射畸变模型，利用窗口内的灰度信息和
位置信息进行平差处理，最终获得精确匹配点位。
[0020] 所述的步骤3)中通过对时序三维点位坐标的差分获取位移数据，将位移数据在时间上进行一次微分和二次微分，分别获得速度数据和加速度数据。
[0021] 与现有技术相比，本发明具有以下优点：
[0022] 本发明通过将数字散斑匹配方法与高速视频测量相结合，可以时刻获取待测岩石在微小时间间隔中的三维形态变化，进而提出了一套完备的非接触式三维形变测量方案，能够为岩石裂缝扩展的定量研究与分析提供多种形变参数。
附图说明
[0023] 图1为本发明的技术路线图。
[0025] 图3为同名匹配示意图。
[0026] 图4为目标点在序列影像中跟踪匹配示意图。
[0027] 图5为高速相机网络布设图
[0028] 图6为归一化相关系数统计直方图。
[0029] 图7为岩石表面的
三维重建，其中，图(7a)为初始时刻的三维点云分布，图(7b)为第8.285s时的三维点云分布。
[0030] 图8为岩石试样在第8.285s时的三维位移场，其中，图(8a)为X方向位移，图(8b)为Y方向位移，图(8c)为Z方向位移。
[0031] 图9为岩石试样在第8.285s时的应变场，其中，图(9a)为Exx应变，图(9b)为Eyy应变，图(9c)为Exy应变。
具体实施方式
[0032] 下面结合附图和具体
实施例对本发明进行详细说明。
[0033] 实施例
[0034] 本
专利针对岩石裂缝扩展实验制订了一套详细的高速视频测量方法。在本实施例中使用了两台高速相机进行立体观测，因此可将该
硬件设备和相关测量理论统称为双目视觉测量系统。该系统的构建需要高速相机、同步控制器、高速
数据采集卡等硬件。其中，高速相机是该系统重要的组成部分，它属于工业相机的一种，具有高
稳定性、高帧频、高传输能力和高抗干扰能力等独特优点，特别适用于
自动光学检测、三维测量、
半导体检测、
机器视觉等领域。此外，同步控制器的作用是使联测的相机保持同步性。在后续的
数据处理中，当解算某一时刻某一目标点的三维空间坐标时，需要获得该时刻拍摄的同名影像。如果不对两台相机进行同步控制，就无法进行后期的三维数据解算，因此同步控制是至关重要的。
[0035] 双目视觉测量系统的网络构建也需要根据实验环境来设计。为了增加影像重叠覆盖率，两台高速相机应使用交向摄影方式。在岩石试样
块的正前方，两台高速相机应尽量水平安放，但在水平面上的相对
姿态不能太过倾斜。这是因为如果倾斜姿态过大，那么左右两组影像的相对透视
变形会相应加大，这会降低后期的亚像素影像匹配
精度。因此在实验实施之前，可根据拍摄影像的
视野范围对两台相机的位置和姿态进行调整。
[0036] 本发明的技术路线，如图1所示。
[0037] 本发明具体包括以下步骤：
[0038] 1、目标匹配与跟踪
[0039] 1.1影像预处理
[0040] 在影像匹配之前，需进行兴趣区提取和目标点位确定，其中目标点选取过程类似于采样过程，而采样的间隔可根据实验的需求而制定。例如，将左相机拍摄的第一张影像(即首帧影像)作为参考影像，先以人工方式选取兴趣区，再通过一定的长度间隔(采样间隔)选取目标点位。采样间隔越小则目标点数量越多，其示意图如图2所示。
[0041] 1.2同名点匹配
[0042] 本实验需要亚像素高精度匹配结果，因此本专利使用最为常规的由粗到精的匹配策略。粗匹配通过计算归一化相关系数来确定整像素匹配粗略点位，精匹配则通过最小二乘匹配方法来确定亚像素精确点位。其中最小二乘匹配方法以归一化相关系数最大为目标并将左右影像变形视为仿射变换，利用窗口内的灰度信息和位置信息进行平差处理，从而可达到1/10甚至1/100像素的匹配精度。为了加快匹配速度，需要确定局部搜索区域。由于本实验中两台相机的相对倾斜姿态较小，致使左右影像中兴趣区的范围相差不大。如图3所示，可先确定左影像某一目标点在左兴趣区中的位置关系，然后再推算出该同名点在右兴趣区中大致的范围，由此来决定影像匹配搜索区域。另外，已经确定了各目标点的点位，因此目标影像块的窗口尺寸和搜索区域的窗口尺寸可以按照实验的需求进行设置。
[0043] 1.3目标跟踪匹配
[0044] 目标跟踪匹配是为了获得各目标点序列影像坐标，其亚像素级匹配方法与同名点匹配相似。不同之处在于匹配对象不再是左右影像，而是各相机存储的序列影像。由于同名点匹配过程已经提供了目标影像块，因此这些影像块也应该在目标跟踪匹配中作为目标影像，而下一帧的搜索区域可由上一帧的目标位置所确定，其跟踪匹配示意图如图4。
[0045] 通过上述流程，每一个目标点都可获得序列影像坐标，由于两台高速相机对所拍摄的影像已经进行了同步采集和存储，因此在首帧影像上获得的同名点在时间序列上依旧保持同名关系。通过三维重建后，岩石表面在每个时间刻度下都可形成三维点云数据。
[0046] 2、双目视觉三维重建
[0047] 2.1立体标定
[0048] 本专利使用了基于平面板的标定方法同步获取相机的内方位元素和外方位元素，其中内方位元素不但要考虑像主点(Cx,Cy)和像距(f)，而且要考虑镜头的畸变参数和像元的实际大小(Sx,Sy)。镜头的畸变参数主要包括径向畸变(K1,K2,K3)与切向畸变(P1,P2)。此外，立体标定不但要确定各相机的内方位元素，同时也要确定相机的外方位元素。外方位元素(α,β,γ,tx,ty,tz)则主要体现在摄影机
坐标系与世界坐标系之间的转换关系。立体标定是本实验最为关键的一步，因为其标定定向的精度影响着最终的实验结果。因此，可在平面标定板上将划取检测点来进一步验证标定精度。
[0049] 2.2三维重建
[0050] 在实验过程中，已标定好的高速相机不可发生移动，否则需要重新标定。在2.1节，立体标定已经确定了各相机的内外方位元素，因此在两台相机采集的序列影像中，每获得一对同名点的像点坐标，便可通过基于共线方程的前方交会解算其三维点位。通过空间和时间上的累积，从而获得大量的点云数据。近景摄影测量的共线条件方程公式如下：
[0051]
[0052] 其中，R＝[a1 b1 c1；a2 b2 c2；a3 b3 c3]。
[0053] 在双目视觉测量系统下，一对同名点可列出4个方程，而需要求解3个未知数，因此可按最小二乘原理进行平差解算。
[0054] 3、形变参数解算
[0055] 在前述的数据处理过程中，已经可以获取任意时刻的三维点云数据。由此可以形成位移、速度、加速度、应变等多种形变参数。
[0056] 3.1位移、速度、加速度
[0057] 位移是指目标点在时间序列上某一时刻的当前位置与该目标点初始时刻的距离差。由此便可知跟踪点初始时刻的位移为0。例如，三维点位移数据的位移公式如下：
[0058]
[0059] 其中， 和 分别表示目标点在X和Y方向在时刻n的位移；X1，Y1和Z1分别表示目标点在X，Y和Z方向于初始时刻的坐标值；Xn，Yn和Zn分别表示目标点在X，Y和Z方向于时刻n的坐标值。
[0060] 由于高速相机的采集帧频是固定不变的，因此可以获得相邻两帧的时间差。若将位移数据分别对时间进行一次微分和二次微分，便可获得相应的速度数据和加速度数据。因此，岩石表面所有目标点都进行以上处理，便可形成位移场、速度场和加速度场。
[0061] 3.2应变值
[0062] 在应变解析中，某一点的应变值可由周围的位移数据计算而成，因而可以将各个目标点为中心，在其周围选取一个位移窗口进行应变值解算。在该位移窗口中，可认为其位移分布是线性的。其位移与坐标的关系可表示为：
[0063]
[0064] 其中，u(i,j)和v(i,j)是位移窗口下的点(i,j)的位移值，ai＝0,1,2和bi＝0,1,2是多项式的待定系数。
[0065] 在较小变形的情况下，应变分量可由下式求解：
[0066]
[0067] 该位移窗口尺寸的大小可根据需求进行确定，如果窗口尺寸较大，也可用高次多项式来表示位移分布。但一般来讲，该位移窗口的尺寸应凭借实验需求选取适中，从公式(3)和公示(4)来看，最少已知三个点的平面坐标便可求解其应变分量。在单轴压缩下，更多关注的是岩石表面内发生的形变。因此，在前期生成的三维点云中，只考虑其空间平面内的应变即可。
[0068] 本实验的观测对象是用来进行裂缝扩展测量的岩石试样块，该试样块是由医用
石膏和水按照一定的比例混合制作而成，经过模具加工成为一个边长为70mm的立方体。为了满足测量的需求，需要喷涂散斑来增加观测点位，其具体过程如下：(1)试样块的观测表面需进行打磨至平整；(2)在观测表面喷涂白色的哑光漆，并进行风干；(3)在观测表面随机并均匀地喷洒黑色哑光漆或黑色墨水，以此来形成散斑影像。此外，如图5所示，两台高速相机形成交向摄影测量方式，且使用大功率
卤素灯对实验对象进行补光，保证影像的拍摄质量。
[0069] 在实验中，两台高速相机形成双目视觉对岩石试样块进行三维测量，其帧频设置为400帧/秒，可精确测量
频率高达40Hz的动态数据，即每10次的影像数据来描述试样1次的形态变化。高速相机所拍影像大小为2304×1720像素，且像元大小为7um。为了精密量测，高速相机配置了50mm的定焦镜头。
[0070] 通过标定，可获得两台高速相机的内方位元素和相对的外方位元素，结果如表1。通过平面板上划定的检测点检核，其标定定向精度可达0.01mm，且反投影误差优于0.2像素。
[0071] 表1高速相机的内方位元素和外方位元素
[0072]
[0073] 按照前述，在匹配之前需要进行影像预处理，本次在兴趣区中选取目标点的采样间隔为5像素，由此便产生了17423个目标观测点，并按131列数和133行数进行规则排列。对所有目标点位确定目标影像块，其中目标窗口的尺寸设置为30像素，同名搜索窗口的尺寸设置为50像素。在同名点匹配之后，可获得归一化相关系数统计图，如图6所示，从图中可以获知，有14084对同名点的相关系数在0.9～1.0这一取值范围，有3315对同名点的相关系数在0.8～0.9这一取值范围。由于相关系数在0.8以上的数值为高度影像相关，因此几乎所有目标点都能匹配到同名点，且匹配精度高。然而部分点相关系数数值偏小，这是由于个别目标影像块的纹理信息较少所造成的匹配精度不高。这种情况可通过数值内插的方式将这些不好的点位解算出来。
[0074] 在所有目标点进行跟踪匹配后，可以获得序列同名像点坐标，再通过三维重建，可以获得任意时刻岩石表面的三维点云数据。如图7，从三维点云数据中也可直观地看出岩石表面发生的三维形变状态。
[0075] 通过形变参数解算可形成三维位移场和应变场，分别如图8和图9。在压缩过程中，岩石的最上方产生了裂缝，进而对岩石产生破坏，发生比较大的形变。此外，从应变场可明显看出，岩石的破裂趋势将向下延伸。
[0076] 本专利使用400帧频/秒的高速相机测量岩石裂缝扩展实验，详细介绍了双目视觉三维重建方法，同时提出了一种散斑影像匹配策略和匹配方法，进而提出了一套完备的三维形变测量方案。
[0077] 1)本方案量测了实验对象在破裂过程中的位移场和应变场动态变化，其中目标点位的空间
定位精度最高可达到0.01mm，验证了整套测量方案在岩石破裂实验中应用的可行性和有效性。
[0078] 2)在散斑影像匹配过程中，目标点的选取间隔和目标影像块窗口可根据实验需求来确定，其增加了数据分析的灵活性。
[0079] 3)高速视频测量技术凭借独有的高帧频特性可测量高速运动物体，采集帧频为400帧/秒的高速相机可精确测量频率高达40Hz的动态数据。
[0080] 4)本专利同样详细阐述了形变参数的求解过程，为进一步研究分析工作提供了可靠的实验数据。
