CN105157592B - 基于双目视觉的柔性后缘可变形机翼的变形形状和速率的测量方法 - Google Patents
基于双目视觉的柔性后缘可变形机翼的变形形状和速率的测量方法 Download PDFInfo
- Publication number
- CN105157592B CN105157592B CN201510530708.2A CN201510530708A CN105157592B CN 105157592 B CN105157592 B CN 105157592B CN 201510530708 A CN201510530708 A CN 201510530708A CN 105157592 B CN105157592 B CN 105157592B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- mtr
- trailing edge
- 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.)
- Active
Links
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于双目视觉的柔性后缘可变形机翼的变形形状和速率的测量方法,属于飞机机翼结构分析和测量技术领域。所述的测量方法安装并标定两台工业相机,获得相机的内参数、外参数和立体参数;然后在可变形机翼的柔性后缘端点和上下蒙皮处粘贴圆形标记点,控制两台相机每间隔一段时间同步拍摄柔性后缘图像;再对柔性后缘图像进行处理,检测圆形标记点圆心的图像平面坐标;计算圆形标记点的三维物方坐标;利用曲线拟合获得柔性后缘上下边缘的形状,计算后缘端点变形速率。本发明解决了柔性后缘结构动态变形中,柔性后缘变形形状和变形速率的测量,避免了接触式测量(如加速度计)对柔性后缘变形形状和速率的影响。
Description
技术领域
本发明属于飞机机翼结构分析和测量技术领域,具体涉及一种柔性后缘可变形机翼的变形形状和速率的测量方法。
背景技术
具有柔性后缘的可变形机翼在飞机的飞行任务中,需要根据不同的飞行环境和任务需求产生主动的柔性变形,改变机翼后缘的弯度,从而改善飞机的气动特性,使飞机在不同的飞行阶段均具有最优的气动性能,其中柔性后缘在主动变形过程中的变形参数,如变形形状和变形速率,对飞机的气动性能产生很大的影响,因此精确地测量具有柔性后缘的可变形机翼的变形形状和速率十分必要。
可变形机翼柔性后缘的变形参数的测量属于大变形动态测量问题,目前,针对可变形机翼变形参数的测量方法主要是利用各种传感器,如加速度计测量变形加速度,但在测量中需与机翼接触带来附加质量的影响;如激光测距传感器,但是只能测量机翼后缘某一点的变形位移,且只能测量小变形的情况;而双目立体视觉测量是一种无接触的适用于大变形的多点测量方法,具有非接触、测量精度高和人工干预少等优点,对于可变形机翼柔性后缘的动态大变形测量十分有效,除了可以测量机翼柔性后缘待测点的三维坐标,提取后缘的变形形状,还可以在测量中补充时间维度的信息,测量后缘变形的速率。
发明内容
本发明的目的在于克服现有技术对可变形机翼变形参数测量的不足,提供一种基于双目视觉的柔性后缘可变形机翼的变形形状和变形速率的测量方法。
为了达到上述目的,本发明采用的技术方案包括以下几个步骤:
步骤一:安装并标定两台工业相机,获得相机的内参数、外参数和立体参数;
步骤二:在可变形机翼的柔性后缘端点和上下蒙皮处粘贴圆形标记点,保持两台工业相机的位置不发生变化,驱动可变形机翼使柔性后缘发生主动变形,控制两台相机每间隔一段时间对柔性后缘进行同步拍摄,得到柔性后缘图像;
步骤三:对同步拍摄的柔性后缘图像进行处理,检测圆形标记点圆心的图像平面坐标。圆形标记点圆心的图像平面坐标的获取包括圆形标记点的边缘检测和hough圆变换,具体方法如下:
(1)将柔性后缘图像进行二值化处理,采用canny算子检测柔性后缘图像中圆形标记点的边缘,并保存边缘点i在像素坐标系的坐标位置(xi,yi);
(2)根据圆的方程:x=a+rcosθ,y=b+rsinθ,其中,a和b分别代表圆心在像素坐标系中的横、纵坐标,r为圆的半径,将圆的边缘点转换到a-b-r的三维参数空间中,每一个边缘点对应三维参数空间中每个半径下的一个圆,实际是一个圆锥,同一个圆形标记点上的所有边缘点对应三维参数空间中的所有三维锥面必然相交于一点,因此检测到这一相交的点即可得到圆形标记点的参数,即圆形标记点圆心的像素坐标(a,b)和圆形标记点的半径r,通过下式可将像素坐标系内的坐标位置转换到图像平面坐标系上的坐标位置(xp,yp):
式中dx、dy分别表示像素坐标系中x和y方向的单位像素的宽度,即一个像素在横轴和纵轴上的物理尺寸,(u0,v0)为图像主点在像素坐标系的坐标位置,图像主点即工业相机中光轴与图像平面的交点;
步骤四:计算圆形标记点的三维物方坐标;
根据图像平面坐标计算圆形标记点三维物方坐标的过程包括图像对的立体校正和对应点的立体匹配,方法如下:
(1)校正立体图像对,将左右相机拍摄的每对柔性后缘图像投影到一个公共平面上,使得两个工业相机的像平面是完全的行对准,采用Bouguet算法进行校正,校正的具体方法如下:
Rl=Rrectrl
Rr=Rrectrr
其中rl和rr为由旋转矩阵R分解得到的合成旋转矩阵,Rrect为变换矩阵。
(2)立体匹配是将左右工业相机同步拍摄的两幅柔性后缘图像(分别称为左图像和右图像)中的对应的圆形标记点进行匹配,具体方法如下:
设左右两台工业相机检测到的圆形标记点的圆心分别为p={p1,p2...,pn}和p'={p1',p2',...,pn'},以每一个圆形标记点的圆心为中心,依据极线约束条件,将其邻域窗口(记为W)的像素值作为该圆形标记点的描述符,对于左图像中的每一个圆形标记点的圆心pi,分别计算其与右图像中各个圆形标记点圆心的邻域像素值差的平方和(SSD),取平方和最小者作为右图像中与左图像中圆形标记点pi的匹配点,平方和的计算公式如下:
步骤五:利用曲线拟合获得柔性后缘上下边缘的形状,得到后缘变形形状;
根据测量的离散点曲线拟合获得柔性后缘上下边缘的形状,方法如下:
曲线拟合采用三次多项式拟合方式,具体的拟合的方法如下:设柔性后缘上边缘测量点的三维坐标为(xi,yi,zi),共有m个点,三次多项式形式为利用曲线拟合的最小二乘法得到参数a0、a1、a2和a3:
步骤六:计算柔性后缘端点变形速率。采用差分原理,具体方法如下:
式中,vx,vy,vz分别为测量点第n秒时在坐标轴三个方向的变形速率,为测量点在第n秒时的三维物方坐标,为测量点在第n+Δt秒时的三维物方坐标,Δt为相邻两张图像拍摄的时间差。
本发明的优点在于:
(1)本发明所述基于双目视觉的柔性后缘可变形机翼的变形形状和速率的测量方法,通过两台工业相机从不同的视角同步拍摄可变形机翼的柔性后缘,柔性后缘粘贴圆形标记点,通过检测圆形标记点圆心的位置,获得圆心的像素坐标,根据投影变换和三角测量原理计算其不同时刻的三维坐标,采用曲线拟合的方法拟合出柔性后缘变形的具体形状,采用差分方法计算柔性后缘端点的变形速率。
(2)本发明解决了柔性后缘结构动态变形中,柔性后缘变形形状和变形速率的测量,避免了接触式测量(如加速度计)对柔性后缘变形形状和速率的影响。
附图说明
图1是双目视觉测量系统与柔性后缘位置关系示意图;
图2是标定板的结构示意图;为圆形阵列标定板,圆形点的个数为14×10,间隔14mm,圆的直径为4mm;
图3是三角测量原理示意图;空间待测点为P,左右相机的光心为Ol和Or,两台工业相机的焦距均为f,两台相机间隔T,Z为空间点的深度信息,视差为d=xl-xr。
图中:1.柔性后缘;2.工业相机A;3.工业相机B;4.铝合金型材。
具体实施方式
下面结合附图对本发明提出的基于双目视觉的柔性后缘可变形机翼的变形形状和速率的测量方法进行详细的说明。
本发明提供的测量方法,对柔性后缘的变形形状和速率进行处理,具体包括以下几个步骤:
步骤一:安装并标定两台工业相机;
如图1所示,将两台工业相机(工业相机A2和工业相机B3)固定在一根铝合金型材4上,保证两台工业相机的相对位置不变。铝合金型材4固定在三脚架上,调整铝合金型材4处于水平位置,调整工业相机A2和工业相机B3的角度使其视野分别覆盖整个柔性后缘,调整两个工业相机焦距f的大小,使得相机的焦点位于柔性后缘1的翼肋所在的平面处,根据环境调整光圈的大小,使得所拍摄的图片亮度合适;手持标定板位于柔性后缘1翼肋所在的平面处,通过单片机控制两台工业相机同步拍摄,然后移动和转动标定板,重复拍摄,每台工业相机所拍摄的图像数量约为15~20张,称为标定板图像,使得标定板与工业相机之间存在足够多的位置关系,以提高标定的准确性,采用的标定板为圆形点阵列标定板,标定板上圆形点的直径为4mm,数量为14×10,间隔(圆心距离)为14mm,如图2所示;然后利用C语言调用OpenCV中的cvstereoCalibrate函数计算两台工业相机的内参数、外参数和立体参数,内参数包括焦距、主点坐标和畸变系数,外参数包括世界坐标系到相机坐标系的旋转矩阵和平移向量,立体参数包括两台工业相机之间的旋转矩阵和平移向量。
所述的铝合金型材4为4040L,与三脚架连接。
步骤二:在可变形机翼的柔性后缘端点和上下蒙皮处粘贴圆形标记点,保持两台工业相机与柔性后缘的相对位置不发生变化,驱动可变形机翼柔性后缘产生主动变形,变形过程中,通过单片机控制两台工业相机同步拍摄柔性后缘图像,保存拍摄的柔性后缘图像。每台工业相机拍摄相邻两张柔性后缘图像的间隔时间Δt=0.02s。
步骤三:对同步拍摄的柔性后缘图像进行处理,检测圆形标记点圆心的图像平面坐标。
双目视觉系统中包含四个坐标系,分别为世界坐标系(也叫物方坐标系)、相机坐标系、图像平面坐标系和像素坐标系,图像平面坐标系与像素坐标系之间的转换关系如下式:
式中(xp,yp)为圆形标记点圆心在图像平面坐标系内的图像平面坐标,(a,b)为圆形标记点圆心在像素坐标系内的像素坐标,dx、dy分别表示像素坐标系中x和y方向的单位像素的宽度,即一个像素在横轴和纵轴上的物理尺寸,(u0,v0)为图像主点的像素坐标,图像主点即工业相机中光轴与图像平面的交点。
圆形标记点圆心的图像平面坐标(xp,yp)的获取包括圆形标记点的边缘检测和hough圆变换,具体方法如下:
(1)计算圆形标记点圆心的图像平面坐标,将柔性后缘图像进行二值化处理,采用canny算子检测圆形标记点的边缘点,并保存边缘点i在像素坐标系中坐标位置(xi,yi);
(2)利用hough圆变换,根据圆的方程:x=a+rcosθ,y=b+rsinθ,其中,a和b分别代表圆心在像素坐标系中的横、纵坐标,r为圆的半径,将圆形标记点的边缘点转换到a-b-r的三维参数空间中,每一个边缘点对应三维参数空间中每个半径r下的一个圆,实际是一个圆锥,同一个圆形标记点上的所有边缘点对应三维参数空间中的所有三维圆锥面必然相交于一点,因此检测到这一相交的点即可得到圆形标记点的参数,即圆形标记点圆心的像素坐标(a,b)和圆形标记点的半径r,通过坐标系的转换关系即式(1),得到圆形标记点圆心的图像平面坐标为(xp,yp)。
步骤四:计算圆形标记点的三维物方坐标;
根据图像平面坐标计算圆形标记点三维物方坐标的过程包括图像对的立体校正和对应点的立体匹配,方法如下:
(1)对左右工业相机拍摄的柔性后缘图像进行立体校正,将左右工业相机得到的每对柔性后缘图像投影到一个公共平面上,使得每对柔性后缘图像的像平面是完全的行对准,采用Bouguet算法进行校正,将右工业相机图像平面旋转到左相机图像平面的旋转矩阵R分解为两部分,分别为合成旋转矩阵rl和rr,取变换矩阵Rrect为:
式中,T为平移向量,Tx和Ty分别为平移向量T在相机坐标系x和y方向的投影向量,通过合成旋转矩阵rl和rr分别与变换矩阵Rrect相乘获得左右工业相机的整体旋转矩阵Rl和Rr,左右工业相机拍摄的柔性后缘图像中圆形标记点圆心的图像平面坐标分别乘以各自的整体旋转矩阵即可完成校正:
Rl=Rrectrl (3)
Rr=Rrectrr
(2)对左右工业相机拍摄的柔性后缘图像进行立体匹配。
立体匹配的目的是将对应的圆形标记点进行匹配,以计算视差;立体匹配的具体方法为:设左右两台工业相机检测到的圆形标记点的圆心分别为p={p1,p2...,pn}和p'={p1',p2',...,pn'},以每一个圆形标记点的圆心为中心,依据极线约束条件,将其邻域窗口(记为W)的像素值作为该圆形标记点的描述符,对于左工业相机拍摄的柔性后缘图像中的每一个圆形标记点的圆心pi,分别计算其与右工业相机拍摄的柔性后缘图像中各个圆形标记点圆心的邻域像素值差的平方和(SSD),取平方和最小者作为右工业相机拍摄的柔性后缘图像中与左工业相机拍摄的柔性后缘图像中圆形标记点pi的匹配点,平方和的计算公式如下:
I1(x,y)和I2(x,y)分别为左工业相机和右工业相机拍摄的柔性后缘图像中圆形标记点邻域窗口W的像素值。
获取圆形标记点圆心的三维世界坐标,将世界坐标系的原点设定在左相机坐标系的原点处,校正后的左右相机的投影矩阵Pl、Pr分别为:
式中,αl和αr是像素畸变比例,一般为0;cx_l,cy_l分别为左工业相机中图像主点与成像中心在图像平面坐标系中两个方向的偏移量,fx_l,fy_l分别为左工业相机在图像平面坐标系中两个方向的焦距,cx_r,cy_r分别代表右工业相机中图像主点与成像中心在图像平面坐标系中两个方向的偏移量,fx_r,fy_r分别代表右工业相机在图像平面坐标系中两个方向的焦距。
进一步计算将平面坐标的二维点投影到三维世界坐标的转换矩阵Q,为:
式中,f表示左工业相机的焦距,c'x表示图像主点在右工业相机拍摄的柔性后缘图像上的x坐标。
根据三角测量原理,如图3所示,在两台相机前向平行的前提下,通过相似三角形计算深度Z,如下式:
若视差取为d=xl-xr,依据三角测量原理,采用齐次坐标的形式,如下式,可得圆形标记点圆心的三维世界坐标
步骤五:利用曲线拟合获得柔性后缘上下边缘的形状,得到后缘变形形状,具有柔性后缘的可变形机翼的变形主要关注翼型的弯曲变形,因此拟合变形形状时,Z方向的位移忽略不计,设柔性后缘上(下)边缘点的三维坐标(Xi,Yi,Zi),共有m个点,采用最小二乘法曲线拟合的方式,拟合采用三次多项式拟合方式,具体的拟合的方法如下:三次多项式形式为利用曲线拟合的最小二乘法得到参数ai:
步骤六:计算柔性后缘端点变形速率,柔性后缘端点不同时刻的变形位移可根据上述方法计算得到,测量过程中,每隔0.02s拍摄一张照片,采用一阶差分方法计算变形速率,具体方法如下:
式中,vx,vy,vz分别为测量点第n秒时在坐标轴三个方向的变形速率,为测量点在第n秒时的三维物方坐标,为测量点在第n+Δt秒时的三维物方坐标,Δt为相邻两张图像拍摄的时间差。
Claims (2)
1.基于双目视觉的柔性后缘可变形机翼的变形形状和速率的测量方法,其特征在于:包括以下几个步骤,
步骤一:安装并标定两台工业相机,获得相机的内参数、外参数和立体参数;
步骤二:在可变形机翼的柔性后缘端点和上下蒙皮处粘贴圆形标记点,保持两台工业相机的相对位置不发生变化,驱动可变形机翼使柔性后缘发生主动变形,控制两台工业相机每间隔一段时间对柔性后缘进行同步拍摄,得到柔性后缘图像;
步骤三:对同步拍摄的柔性后缘图像进行处理,检测圆形标记点圆心的图像平面坐标;
所述的圆形标记点圆心的图像平面坐标的获取包括圆形标记点的边缘检测和hough圆变换,具体方法如下:
(1)将柔性后缘图像进行二值化处理,采用canny算子检测柔性后缘图像中圆形标记点的边缘,并保存边缘点在像素坐标系的坐标位置(xi,yi);
(2)根据圆的方程:x=a+r cosθ,y=b+r sinθ,其中,a和b分别代表圆心在像素坐标系中的横、纵坐标,r为圆的半径,将圆的边缘点转换到a-b-r的三维参数空间中,每一个边缘点对应三维参数空间中每个半径下的一个圆,实际是一个圆锥,同一个圆形标记点上的所有边缘点对应三维参数空间中的所有三维锥面必然相交于一点,因此检测到这一相交的点即得到圆形标记点的参数,即圆形标记点圆心的像素坐标(a,b)和圆形标记点的半径r,通过下式将像素坐标系内的坐标位置转换到图像平面坐标系上的坐标位置(xp,yp):
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mi>p</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<mi>p</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>u</mi>
<mn>0</mn>
</msub>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>d</mi>
<mi>y</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
<mi>d</mi>
<mi>y</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>a</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>b</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中dx、dy分别表示像素坐标系中x和y方向的单位像素的宽度,即一个像素在横轴和纵轴上的物理尺寸,(u0,v0)为图像主点在像素坐标系的坐标位置,图像主点即工业相机中光轴与图像平面的交点;
步骤四:计算圆形标记点的三维物方坐标,包括图像对的立体校正和对应点的立体匹配;具体实现方法如下,
(1)校正立体图像对,将左右相机拍摄的每对柔性后缘图像投影到一个公共平面上,使得两个工业相机的像平面是完全的行对准,采用Bouguet算法进行校正,将右工业相机图像平面旋转到左相机图像平面的旋转矩阵R分解为两部分,分别为合成旋转矩阵rl和rr,取变换矩阵Rrect为:
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>c</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>e</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>e</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mi>T</mi>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>T</mi>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>,</mo>
<msub>
<mi>e</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mi>x</mi>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<msqrt>
<mrow>
<msubsup>
<mi>T</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>T</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mfrac>
<mo>,</mo>
<msub>
<mi>e</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<mo>&times;</mo>
<msub>
<mi>e</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,T为平移向量,Tx和Ty分别为平移向量T在相机坐标系x和y方向的投影向量,通过合成旋转矩阵rl和rr分别与变换矩阵Rrect相乘获得左右工业相机的整体旋转矩阵Rl和Rr,左右工业相机拍摄的柔性后缘图像中圆形标记点圆心的图像平面坐标分别乘以各自的整体旋转矩阵即可完成校正,具体方法如下:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>l</mi>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>c</mi>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>r</mi>
<mi>l</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>c</mi>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>r</mi>
<mi>r</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
(2)立体匹配是将左右工业相机同步拍摄的两幅柔性后缘图像中的对应的圆形标记点进行匹配,具体方法如下:
设左右两台工业相机检测到的圆形标记点的圆心分别为p={p1,p2...,pn}和p'={p1',p2',...,pn'},以每一个圆形标记点的圆心为中心,依据极线约束条件,将其邻域窗口的像素值作为该圆形标记点的描述符,对于左图像中的每一个圆形标记点的圆心pi,分别计算其与右图像中各个圆形标记点圆心的邻域像素值差的平方和SSD,取平方和最小者作为右图像中与左图像中圆形标记点pi的匹配点,平方和的计算公式如下:
<mrow>
<mi>S</mi>
<mi>S</mi>
<mi>D</mi>
<mo>=</mo>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>&Element;</mo>
<mi>W</mi>
</mrow>
</munder>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>I</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>I</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
I1(x,y)和I2(x,y)分别为左工业相机和右工业相机拍摄的柔性后缘图像中圆形标记点邻域窗口W的像素值;
步骤五:利用曲线拟合获得柔性后缘上下边缘的形状,得到后缘变形形状;
步骤六:计算柔性后缘端点变形速率;
采用差分原理,具体方法如下:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>z</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>z</mi>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,vx,vy,vz分别为测量点第n秒时在坐标轴三个方向的变形速率,为测量点在第n秒时的三维物方坐标,为测量点在第n+Δt秒时的三维物方坐标,Δt为相邻两张图像拍摄的时间差。
2.根据权利要求1所述的一种基于双目视觉的柔性后缘可变形机翼的变形形状和速率的测量方法,其特征在于:每台工业相机拍摄相邻两张柔性后缘图像的间隔时间Δt=0.02s。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510530708.2A CN105157592B (zh) | 2015-08-26 | 2015-08-26 | 基于双目视觉的柔性后缘可变形机翼的变形形状和速率的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510530708.2A CN105157592B (zh) | 2015-08-26 | 2015-08-26 | 基于双目视觉的柔性后缘可变形机翼的变形形状和速率的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105157592A CN105157592A (zh) | 2015-12-16 |
CN105157592B true CN105157592B (zh) | 2018-03-06 |
Family
ID=54798536
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510530708.2A Active CN105157592B (zh) | 2015-08-26 | 2015-08-26 | 基于双目视觉的柔性后缘可变形机翼的变形形状和速率的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105157592B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105571491B (zh) * | 2016-02-22 | 2018-01-09 | 麦特汽车服务股份有限公司 | 基于双目视觉的汽车底盘数据测量系统及其方法 |
CN106323587B (zh) * | 2016-08-03 | 2018-07-27 | 中国空气动力研究与发展中心高速空气动力研究所 | 机翼风洞试验模型弹性变形的单目视频高精度测量方法 |
US20180099761A1 (en) * | 2016-10-10 | 2018-04-12 | The Boeing Company | Real Time Aircraft Stress Monitoring System |
CN106289184B (zh) * | 2016-11-01 | 2018-11-13 | 华中师范大学 | 一种无gnss信号和无控制点下无人机协同视觉形变监测方法 |
CN107167092B (zh) * | 2017-05-18 | 2019-12-13 | 上海晶电新能源有限公司 | 一种基于多目图像识别的定日镜面形检测系统及方法 |
CN107576275A (zh) * | 2017-08-11 | 2018-01-12 | 哈尔滨工业大学 | 一种利用摄影测量技术对充气结构进行应变场测量的方法 |
CN109559349B (zh) * | 2017-09-27 | 2021-11-09 | 虹软科技股份有限公司 | 一种用于标定的方法和装置 |
CN108844486B (zh) * | 2018-04-12 | 2021-02-09 | 西安交通大学 | 一种定像距双目仿生三维测量仪器 |
CN108955553B (zh) * | 2018-07-26 | 2023-06-06 | 苏州科技大学 | 基于激光数字散斑干涉定位的多点应变测量装置和方法 |
CN110425996A (zh) * | 2019-07-31 | 2019-11-08 | 工极智能科技(苏州)有限公司 | 基于双目立体视觉的工件尺寸测量方法 |
CN111325802B (zh) * | 2020-02-11 | 2023-04-25 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种直升机风洞试验中圆形标记点识别匹配方法 |
CN111238405B (zh) * | 2020-04-26 | 2020-08-04 | 武汉精立电子技术有限公司 | 一种工件曲率检测系统及方法 |
CN112268548B (zh) * | 2020-12-14 | 2021-03-09 | 成都飞机工业(集团)有限责任公司 | 一种基于双目视觉的飞机局部外形测量方法 |
CN112797912B (zh) * | 2020-12-24 | 2023-04-07 | 中国航天空气动力技术研究院 | 一种基于双目视觉的大柔性无人机机翼翼尖变形测量方法 |
CN112801892B (zh) * | 2021-01-11 | 2022-07-01 | 中国民航大学 | 一种基于FFWA的自适应Canny飞机蒙皮红外图像边缘检测方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102944183A (zh) * | 2012-10-31 | 2013-02-27 | 中国航天空气动力技术研究院 | 一种大展弦比柔性机翼测量方法 |
CN103278138A (zh) * | 2013-05-03 | 2013-09-04 | 中国科学院自动化研究所 | 一种复杂结构薄部件三维位置及姿态的测量方法 |
CN103337094A (zh) * | 2013-06-14 | 2013-10-02 | 西安工业大学 | 一种应用双目摄像机实现运动三维重建的方法 |
CN103575227A (zh) * | 2013-09-23 | 2014-02-12 | 西安新拓三维光测科技有限公司 | 一种基于数字散斑的视觉引伸计实现方法 |
CN103925878A (zh) * | 2014-04-11 | 2014-07-16 | 南京信息工程大学 | 基于Hough圆卷积变换算法的孔组检测方法与系统 |
CN104236478A (zh) * | 2014-09-19 | 2014-12-24 | 山东交通学院 | 一种基于视觉的车辆外廓尺寸自动测量系统及方法 |
CN104853154A (zh) * | 2015-04-29 | 2015-08-19 | 同济大学 | 一种基于运动图像的危岩变形信息提取及报警方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5210203B2 (ja) * | 2009-02-25 | 2013-06-12 | ローランドディー.ジー.株式会社 | 画像の差分による高精度ステレオカメラキャリブレーション |
JP2012058076A (ja) * | 2010-09-09 | 2012-03-22 | 3D Media Co Ltd | 3次元計測装置及び3次元計測方法 |
-
2015
- 2015-08-26 CN CN201510530708.2A patent/CN105157592B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102944183A (zh) * | 2012-10-31 | 2013-02-27 | 中国航天空气动力技术研究院 | 一种大展弦比柔性机翼测量方法 |
CN103278138A (zh) * | 2013-05-03 | 2013-09-04 | 中国科学院自动化研究所 | 一种复杂结构薄部件三维位置及姿态的测量方法 |
CN103337094A (zh) * | 2013-06-14 | 2013-10-02 | 西安工业大学 | 一种应用双目摄像机实现运动三维重建的方法 |
CN103575227A (zh) * | 2013-09-23 | 2014-02-12 | 西安新拓三维光测科技有限公司 | 一种基于数字散斑的视觉引伸计实现方法 |
CN103925878A (zh) * | 2014-04-11 | 2014-07-16 | 南京信息工程大学 | 基于Hough圆卷积变换算法的孔组检测方法与系统 |
CN104236478A (zh) * | 2014-09-19 | 2014-12-24 | 山东交通学院 | 一种基于视觉的车辆外廓尺寸自动测量系统及方法 |
CN104853154A (zh) * | 2015-04-29 | 2015-08-19 | 同济大学 | 一种基于运动图像的危岩变形信息提取及报警方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105157592A (zh) | 2015-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105157592B (zh) | 基于双目视觉的柔性后缘可变形机翼的变形形状和速率的测量方法 | |
CN107449402B (zh) | 一种非合作目标的相对位姿的测量方法 | |
CN105698699B (zh) | 一种基于时间转轴约束的双目视觉测量方法 | |
CN103345737B (zh) | 一种基于误差补偿的uav高分辨率影像几何校正方法 | |
CN105716542B (zh) | 一种基于柔性特征点的三维数据拼接方法 | |
CN104459183B (zh) | 一种基于物联网的单相机车辆测速系统及方法 | |
CN104422425B (zh) | 一种不规则外形物体空间姿态动态测量方法 | |
CN101813465B (zh) | 非接触精密测转角的单目视觉测量方法 | |
CN109685858A (zh) | 一种单目摄像头在线标定方法 | |
CN104034305B (zh) | 一种单目视觉实时定位的方法 | |
CN104240262B (zh) | 一种用于摄影测量的相机外参数标定装置及标定方法 | |
CN106949836B (zh) | 一种立体视觉摄像头同侧目标位置标定装置及方法 | |
CN104240221B (zh) | 一种镜头相背双相机相对方位标定装置与方法 | |
CN103759669A (zh) | 一种大型零件的单目视觉测量方法 | |
CN103278138A (zh) | 一种复杂结构薄部件三维位置及姿态的测量方法 | |
CN106197265B (zh) | 一种空间自由飞行模拟器视觉精密定位方法 | |
CN105469389B (zh) | 一种用于视觉传感器标定的网格球靶标及相应标定方法 | |
CN106338287A (zh) | 基于天花板的室内移动机器人视觉定位方法 | |
CN106971408A (zh) | 一种基于时空转换思想的摄像机标定方法 | |
CN110610518A (zh) | 一种基于特征点位置校正的靶标位姿测量系统及方法 | |
CN104807477B (zh) | 一种基于靶标控制点的卫星ccd阵列影像几何检校方法 | |
CN102589530A (zh) | 基于二维相机和三维相机融合的非合作目标位姿测量方法 | |
CN110470226A (zh) | 一种基于无人机系统的桥梁结构位移测量方法 | |
CN101980292B (zh) | 一种基于正八边形模板的车载摄像机内参数的标定方法 | |
CN104574415A (zh) | 一种基于单摄像机的目标空间定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |