一种基于蜂巢模型的飞机蒙皮小型损伤定位方法与流程

本发明属于视觉标定测量技术领域，特别是涉及一种基于蜂巢模型的飞机蒙皮小型损伤定位方法。

背景技术：

飞机蒙皮是包围在机身骨架外的围型结构，用粘接剂、铆钉或者螺丝固定于骨架上，飞机蒙皮除了形成和维持飞机的启动外形以外，还需要承受局部气动力。但由于飞机蒙皮长期暴露在空气中，因此是飞机结构中极易产生损伤的结构，在航空史上由于飞机蒙皮损伤而造成的飞机事故屡见不鲜。由于飞机的运行环境恶劣，除正常飞行带来的疲劳损伤外，还有由于环境因素导致的结构腐蚀等，这些缺陷的产生对飞机蒙皮结构的完整性产生了破坏，从而引起结构强度的下降，对于飞行存在一定的安全性隐患。疲劳裂纹对飞机结构强度的影响很大，是导致民用航空器结构发生灾难性破坏的主要原因，所以需定期对飞机蒙皮进行结构检测与修理，以保证飞行安全。但是在对飞机蒙皮进行损伤检测的过程中，对飞机蒙皮复杂曲面的破损点的检测定位成为关键的一步。

无损检测(non-destructivetesting，简称ndt)是一门新兴的综合性应用科学。其通过对物理和化学相关理念的运用，可以实现对各种工业产品、工程材料以及零部件等进行有效的检测和测试，借以评价它们的完整性、连续性、安全可靠性及力学、物理性能等，同时不改变被检测对象的状态和使用性能。由于无损检测具有不破坏试件就可以对其进行百分之百的检测等一系列优点，因而在很多领域上得到了重视并迅速发展。飞机蒙皮的检测要求达到较高的检测灵敏度和较高的缺陷查出概率，也就是要求以较高的检出概率查出较小的缺陷。基于无损检测的不改变检测对象的状态和使用性能优点，利用立体视觉测试技术能够有效地对飞机蒙皮的小型损伤进行检测，该技术不仅适用于对此类大型工件的非接触测量，而且其自动化程度高，针对性强。而在对飞机蒙皮进行损伤检测的过程中，对飞机蒙皮复杂曲面的破损点的检测定位成为关键的一步。

技术实现要素：

为了解决上述问题，本发明的目的在于提供一种基于蜂巢模型的飞机蒙皮小型损伤定位方法。

为了达到上述目的，本发明提供的基于蜂巢模型的飞机蒙皮小型损伤定位方法包括按顺序进行的下列步骤：

步骤1)建立飞机蒙皮小型损伤定位系统；

所述的系统包括测量旋转台、两个摄像机、同步频闪控制装置和计算机；其中：测量旋转台为放置待检测的飞机蒙皮的可控旋转平台；摄像机为包含照明光源和摄像头的图像采集装置，两个摄像机分别置于测量旋转台的两侧上方，且摄像机镜头对准待检测的飞机蒙皮；同步频闪控制装置为图像采集同步控制装置，分别与两个摄像机相连接，用于实现两个摄像机的图像同步采集；计算机分别与两个摄像机及测量旋转台相连接，通过两个摄像机采集待检测的飞机蒙皮的图像信息并进行处理，并控制测量旋转台的旋转速度；

步骤2)利用上述飞机蒙皮小型损伤定位系统中的计算机对两个摄像机进行标定；

步骤3)在待检测的飞机蒙皮表面的破损部位布置蜂巢模型标记点；

步骤4)利用上述飞机蒙皮小型损伤定位系统中的两个摄像机对待检测的飞机蒙皮上的破损部位进行拍摄，获得原始图像；

步骤5)对上述拍摄的原始图像进行亚像素级处理，得到标记点的像素边缘；

步骤6)从上述像素边缘中提取出六边形边界点；

步骤7)根据上述六边形边界点计算出经过标记点中心的直线以及该直线的交点坐标；

步骤8)根据上述标记点的交点坐标计算出标记点的实际中心位置。

在步骤1)中，所述的摄像机采用高速ccd相机。

在步骤2)中，所述的利用上述飞机蒙皮小型损伤定位系统中的计算机对两个摄像机进行标定的方法为非参数模型校准方法，结合摄像机成像原理及垂线法，直接建立空间不同方位与相面位置的对应关系，通过将成像平面细化分割，利用插值解算得到任意被测未知点的偏移量，从而得到高精度标定结果。

在步骤3)中，所述的在待检测的飞机蒙皮表面的破损部位布置蜂巢模型标记点的方法是在待检测的飞机蒙皮表面的破损部位布置蜂巢模型的标记点，所述的标记点采用正六边形反光标记点，由多个正六边形反光标记点构成蜂巢模型。

在步骤5)中，所述的对上述拍摄的原始图像进行亚像素级处理，得到标记点的像素边缘的方法是：对步骤4)采集得到的原始图像进行二值化以及形态学处理，再采用自适应中值滤波的方法滤除噪声，从而得到标记点的像素边缘。

在步骤6)中，所述的从上述像素边缘中提取出六边形边界点的方法是：对步骤5)中经过亚像素级处理得到的像素边缘进行灰度加权计算得到边缘灰度值，从而提取出正六边形反光标记点的中心坐标；并且对于每个像素边缘点，沿着该点到正六边形反光标记点中心的方向计算出径向灰度梯度，然后通过灰度加权质心算法求得六边形边界点。

在步骤7)中，所述的根据上述六边形边界点计算出经过标记点中心的直线以及该直线的交点坐标的方法是：由于边界点符合外接圆投影形成的椭圆曲线的性质，所以利用椭圆拟合与正六边形对角线所在的直线的交点得到方程，计算出经过正六边形反光标记点中心的直线以及该直线的交点坐标。

在步骤8)中，所述的根据上述标记点的交点坐标计算出标记点的实际中心位置的方法是根据调和共轭原理和交比不变性计算出标记点的实际中心位置。

本发明提供的基于蜂巢模型的飞机蒙皮小型损伤定位方法能够实现飞机蒙皮表面标记点间的特征对比与匹配，进而实现对飞机蒙皮破损部位的精确检测。其意义在于利用圆的内接正六边形蜂巢模型的无限与无缝隙外扩性来进行定位设计，以达到飞机蒙皮破损检测的高匹配效率、高检测精度，而且计算量小、计算过程简单。若使用双目匹配方法，则更加方便，效果更加突出。

附图说明

图1为本发明提供的基于蜂巢模型的飞机蒙皮小型损伤定位方法所采用的飞机蒙皮小型损伤定位系统结构示意图。

图2为本发明提供的基于蜂巢模型的飞机蒙皮小型损伤定位方法流程图。

图3为正六边形蜂巢模型示意图。

图4为圆形内接正六边形的投影示意图。

图5为正六边形中某一对角线的交比不变性原理示意图。

具体实施方式

下面结合附图和具体实施例对本发明提供的基于蜂巢模型的飞机蒙皮小型损伤定位方法进行详细说明。

如图1、图2所示，本发明提供的基于蜂巢模型的飞机蒙皮小型损伤定位方法包括按顺序执行的下列步骤：

步骤1)建立如图1所示的飞机蒙皮小型损伤定位系统；

所述的系统包括测量旋转台1、两个摄像机2、同步频闪控制装置3和计算机4；其中：测量旋转台1为放置待检测的飞机蒙皮5的可控旋转平台；摄像机2为包含照明光源和摄像头的图像采集装置，两个摄像机2分别置于测量旋转台1的两侧上方，且摄像机镜头对准待检测的飞机蒙皮5；同步频闪控制装置3为图像采集同步控制装置，分别与两个摄像机2相连接，用于实现两个摄像机2的图像同步采集；计算机4分别与两个摄像机2及测量旋转台1相连接，通过两个摄像机2采集待检测的飞机蒙皮5的图像信息并进行处理，并控制测量旋转台1的旋转速度；

步骤2)利用上述飞机蒙皮小型损伤定位系统中的计算机对两个摄像机进行标定；

将上述系统中的两个摄像机2的主光轴夹角设定为60°；根据双目立体视觉原理，利用计算机4对两个摄像机2进行系统标定，获取其内外参数，得到摄像机投影矩阵；

所述的对两个摄像机2进行系统标定的方法为非参数模型校准方法，结合摄像机成像原理及垂线法，直接建立空间不同方位与相面位置的对应关系，通过将成像平面细化分割，利用插值解算得到任意被测未知点的偏移量，从而得到高精度标定结果。

步骤3)在待检测的飞机蒙皮表面的破损部位布置蜂巢模型标记点；

在对待检测的飞机蒙皮5进行破损部位检测时，由于破损部位分布不均匀以及某些破损部位的隐蔽性，所以首先需要在待检测的飞机蒙皮5表面的破损部位布置蜂巢模型的标记点，如图3所示，所述的标记点采用正六边形反光标记点，由多个正六边形反光标记点构成蜂巢模型，以取代传统的圆形标记点。由于蜂巢模型具有无限与无缝隙外扩性，使得该模型能够满足标记点的对比传递，所以可以达到高匹配效率和高检测精度；其中第一标记点1位于蜂巢模型的中心，通过对第一标记点1进行定位，能够快速地进行第一标记点1的特征提取和识别；第二标记点2为第一标记点1的相邻标记点，通过对第二标记点2进行检测和定位，能够与第一标记点1进行很好的特征匹配与对比；以此类推。

由于采用传统的圆形标记点时投影到摄像机成像平面上的图形为椭圆，在此内接正六边形的六个投影点也是位于椭圆的圆周上，其性质符合椭圆方程。下面我们先对圆形蜂巢模型进行分析。

设定x-y平面为圆平面，在此平面中圆特征满足如下式所示的二维圆方程：

(x-x0)²+(y-y0)²＝r²

其等效矩阵形式为：

p^tbp＝0

式中：为圆的矩阵形式，p＝[xy1]^t。

同理可得椭圆的方程与矩阵形式为：

ax²+by²+cxy+dx+ey+f＝0

p^tcp＝0

式中：p＝[xy1]^t。

当圆在物空间x-y平面上且圆心坐标位于坐标原点处，即为令psc＝[xsys1]^t表示圆周上的一个任意点，得：

psc^tcspsc＝0

根据小孔成像模型以及张正友标定法，得到该任意点在图像坐标系下的坐标pcp：

令结合上两式，得到以下公式：

sci＝(h^-1)^tcsh^-1

式中：pie^tcipie＝0；上式中的ci符合椭圆矩阵的形式这说明空间圆在像平面内的投影是一个椭圆。椭圆方程的系数(a,b,c,d,e,f)是圆半径r的函数：

由椭圆圆心计算公式可以得知，圆心是关于圆半径r的函数。计算公式如下：

式中：m1,m2,m,n1,n2和n是h11,h12,……,h33的函数。

在上式中，当圆半径r置零时，(cx,cy)是实际的圆心：

椭圆圆心和实际圆心的斜率可以表示为：

通过上式可以得出，斜率k不随圆半径r的变化而改变，证明了所有圆的圆心在同一直线上，并且该直线通过实际的圆心。

同样的原理适用于圆的内接正六边形。由于正六边形的六个点都位于圆周上，所以可以取对角的两个点并连线，此线必经过实际的圆心。由于蜂巢模型能够无限外扩，且边缘紧密连接，因此能够快速地进行标记点与标记点间的特征匹配与识别，这也是正六边形有别于其他多边形的特点，所以能够快速地完成对待检测的飞机蒙皮5表面破损部位的定位和检测。同时，由于正六边形反光标记点具有同向高反射效应，因此需要安装照明光源；

步骤4)利用上述飞机蒙皮小型损伤定位系统中的两个摄像机对待检测的飞机蒙皮上的破损部位进行拍摄，获得原始图像；

所述的摄像机2采用高速ccd相机，型号为imperx4m，感光尺寸为2/3，分辨率为2060×2056pixel；设置较小的摄像机曝光时间，以保证得到清晰的图像；

步骤5)对上述拍摄的原始图像进行亚像素级处理，得到标记点的像素边缘；

对步骤4)采集得到的原始图像进行二值化以及形态学处理，再采用自适应中值滤波的方法滤除噪声，从而得到标记点的像素边缘；

步骤6)从上述像素边缘中提取出六边形边界点；

对步骤5)中经过亚像素级处理得到的像素边缘进行灰度加权计算得到边缘灰度值，从而提取出正六边形反光标记点的中心坐标；并且对于每个像素边缘点，沿着该点到正六边形反光标记点中心的方向计算出径向灰度梯度，然后通过灰度加权质心算法求得六边形边界点；

步骤7)根据上述六边形边界点计算出经过标记点中心的直线以及该直线的交点坐标；

由于边界点符合外接圆投影形成的椭圆曲线的性质，所以利用椭圆拟合与正六边形对角线所在的直线的交点得到方程，计算出经过正六边形反光标记点中心的直线以及该直线的交点坐标；具体方法如下：图4和图5分别为圆形内接正六边形的投影示意图和正六边形中某一对角线的交比不变性原理示意图，图中的大小写对应表示同一特征的投影关系，如a是a的投影。如图4所示，点0为发光原点，虚线为光束。在物平面内对于由abcdef所形成的正六边形的投影为像平面内的多边形abcdef。如图5所示，线段eb为正六边形中任意一对角线，并且点e、b位于正六边形的外接圆上，所以说线段eb的中点p即为圆心。同理在像平面内所形成的投影线段eb也符合这一性质，e、b为不规则多边形abcdef中的两点，且abcdef为正六边形的投影，所以e、b即为外接圆所形成的椭圆投影上的两点，符合椭圆拟合方程的性质，点p为线段eb的中点(实际圆心的投影)。点g、h分别为连接线段fd、ac与eb所形成的交点，按照正六边形的性质分析可以得出，点p为线段gh的中点。同理在投影图像中的性质亦是相同。点e、g、h、b的坐标均可通过联立方程求解。

步骤8)根据上述标记点的交点坐标计算出标记点的实际中心位置；

根据调和共轭原理和交比不变性计算出标记点的实际中心位置。具体方法如下：在像平面内，利用调和共轭点的性质以及交比不变性来实现对蜂巢模型的标记点的精确定位。

对于调和共轭点来说，如果p1、p2、pm和p∞是四个共线点且满足方程：

{p1,p2；pm,p∞}＝-1

其中{p1,p2；pm,p∞}＝-1是交比方程，那么称p1和p2为pm和p∞的调和共轭点。如果点pm是p1和p2的中点，那么点p∞是消隐点。由共轭调和点的性质分析可得：

{e,b；p,q}＝{g,h；p,q}＝-1

其中，点e、b是正六边形对应顶点，点p是e、b连线中点，点q是e、b所在直线上的远点，点g、h是e、b连线与其它对角线的交点。由交比不变性可得：

{e,b；p,q}＝{g,h；p,d}＝-1

其中，点e、b、g、h分别是点e、b、g、h在像平面的投影点。p和q是由上式得到的二元二次方程的两个解。根据点p在线段eb上，可以选出点p正确的解，此解就是所求的标记点的实际中心位置。此方法计算量小、计算过程简单。

以上结合附图对本发明的具体实施方式作了说明，但这些说明不能被理解为限制了本发明的范围，本发明的保护范围由随附的权利要求书限定，任何在本发明权利要求基础上的改动都是本发明的保护范围。