一种物体表面高精度定位方法及系统与流程

1.本发明涉及测量技术领域，尤其涉及一种物体表面高精度定位方法及系统。


背景技术：

2.对物体表面进行高精度定位是测量领域的重要内容。通过物体表面的高精度定位，可以直接高精度测量该物体的位移或转动，进而可以为精密仪器、精密设备或高精度武器装备等提供高精度的位移或转角数据。此类技术可以广泛地应用于机器人、航空航天装备、光电经纬仪、雷达、数控机床及各种工业自动化设备等领域，对国防、航空航天、精密制造、自动化产业等行业都具有重要意义。
3.目前，物体表面高精度定位的主要方式大多需要依靠光栅来实现。光栅测量技术是一种利用莫尔条纹进行精密测量的技术。其测量原理是：两个光栅重叠时会产生莫尔条纹，莫尔条纹能够对两个光栅的相对微小位移进行放大，通过测量莫尔条纹的相关信息即可实现对光栅的高精度定位，从而实现位移或角度的高精度测量。虽然光栅测量技术具有体积小、精度高、抗干扰能力强等优点，但对材质、工艺的要求非常高，制造加工很困难；即使完成了高精度光栅的加工制造，在进行实际应用时，光栅测量技术依然会受到多种因素制约，实际测量精度受到很大影响，例如码盘的母板误差和刻线误差等编码器光学部分的影响、轴承与结构等机械部分的影响、光发射源与接收单元等电气部分影响以及编码器使用中的信号延迟的影响等。


技术实现要素：

4.为了解决上述技术问题，本发明的目的是提供一种物体表面高精度定位方法及系统，此算法为半监督学习方法，能够通过单台相机拍摄物体表面的自带特征或制作的纹理特征，实现物体二维平移、三维平移、二维转角等参数的高精度测量。
5.本发明所采用的第一技术方案是：一种物体表面高精度定位方法，包括以下步骤：
6.在物体表面制作预设分布密度的标记，得到带标记的物体；
7.制作的标记需要符合如下设置：1.确保在对物体表面成像后，图像中不同标记的背景图有一定区分度，从而可以快速可靠地从图像中检测和识别标记；2.标记在物体表面需要达到一定的分布密度，保证相机在对物体表面的任何待测局部拍摄时，成像视场内至少存在1个完整标记。
8.对带标记的物体进行拍摄，并编号图像中的标记，得到带编号标记的图像；
9.有3种不同命名编号的方式，可以根据实际情况选择其一：1.根据标记自身的图像特征信息进行命名编号；2.根据标记周围的随机纹理图像信息进行命名编号；3.标记自身信息结合其周围纹理的图像信息进行命名编号。
10.根据带编号标记的图像获取物体表面各个标记的相对位置信息并建立标记编号与标记位置信息的对应关系；
11.获取实时待测图像并检测实时待测图像内的标记，得到待测标记；
12.根据待测标记在实时待测图像中的位置以及该标记对应的物体表面位置信息，计算得到当前实时待测图像所对应的物体表面位置。
13.进一步，所述在物体表面制作预设分布密度的标记，得到带标记的物体这一步骤，其具体包括：
14.采用预设样式的标记，通过标记制作方法表面刻制标记，得到带标记的物体；
15.所述预设样式的标记包括线段、点、圆、方形、十字及其组合形状；
16.所述标记制作方法包括激光雕刻、印刷和刻蚀；
17.所述标记的刻制间距小于图像画幅的一半。
18.进一步，所述对带标记的物体进行拍摄，并编号图像中的标记，得到带编号标记的表面图像这一步骤，其具体包括：
19.对带标记的物体的表面拍摄成像，得到表面图像；
20.通过阈值分割对表面图像进行处理，检测得到表面图像内的所有标记的目标中心；
21.以目标中心为参考点，截取预设像素大小的图像作为对应标记的标记特征图像；
22.采用主成分分析法建立标记特征图像与标记编号的对应关系，将标记特征图像降至低维，并取主成分矩阵前10个主成分为表示该标记的一维向量，确定唯一编号；
23.得到带编号标记的表面图像。
24.进一步，所述根据带编号标记的图像获取物体表面各个标记的相对位置信息并建立标记编号与标记位置信息的对应关系这一步骤，其具体包括：
25.通过图像拼接方法将带编号标记的图像进行图像拼接，得到完整表面图像；
26.以完整表面图像为基准图并建立基准图的图像像素坐标与实际物体表面坐标的关系；
27.在基准图像上进行标记的识别和定位，并根据基准图像像素坐标与实际物体表面坐标的关系，确定各个标记与物体表面坐标的对应关系。
28.进一步，所述获取实时待测图像并检测实时待测图像内的标记，得到待测标记这一步骤，其具体包括：
29.根据变换公式确定各标记中心点在物体表面的精确位置；
30.对物体表面进行实时成像，得到实时待测图像；
31.对实时待测图像进行标记检测，基于向量比对确认实时待测图像内标记对应的标记编号。
32.进一步，所述变换公式如下式：
[0033][0034]
上式中，s为中心点在摄像机坐标系中的z坐标，(u,v)为中心图像坐标，(f
x
,fy)为相机的等效焦距，(c
x
,cy)为图像主点，r
ij
(i,j＝1,2,3)为旋转矩阵r元素，ti(i＝1,2,3)为平移向量，(xw,yw,zw)代表物体表面与标记中心对应的三维点坐标。
[0035]
进一步，所述根据待测标记在实时待测图像中的位置以及该标记对应的物体表面位置信息，计算得到当前实时待测图像所对应的物体表面位置这一步骤，其具体包括：
[0036]
根据标记编号获取实时图像中心在基准图像上的对应坐标；
[0037]
根据基准图像像素坐标与物体表面参数的对应关系、实时图像中心在基准图像上的坐标，确定实时图像中心对应的物体表面点；
[0038]
根据实时图像中心对应的物体表面点，得到当前实时待测图像所对应的物体表面位置；
[0039]
所述当前实时待测图像所对应的物体表面位置包括一维坐标、二维坐标、三维坐标、一维转角和二维转角。
[0040]
进一步，所述根据待测标记在实时待测图像中的位置以及该标记对应的物体表面位置信息，计算得到当前实时待测图像所对应的物体表面位置这一步骤，其具体包括：
[0041][0042]
上式中，θ为待测旋转体的旋转角度，y0为旋转体外表面展开的周线原点，yc为旋转体外表面展开的标识中心坐标点，l为旋转体外表面展开的周线长。
[0043]
进一步，所述二维转角的计算公式如下：
[0044][0045]
上式中，图像中心坐标为(u,v)对应在待测球体表面的世界坐标记为(xw,yw,zw)，α为待测球体的经度，β为待测球体的纬度。
[0046]
本发明所采用的第二技术方案是：一种物体表面高精度定位系统，包括：
[0047]
标记制作模块，用于在物体表面制作预设分布密度的标记，得到带标记的物体；
[0048]
编号模块，用于对带标记的物体进行拍摄，并编号图像中的标记，得到带编号标记的图像；
[0049]
标记位置关系模块，用于根据带编号标记的图像获取物体表面各个标记的相对位置信息并建立标记编号与标记位置信息的对应关系；
[0050]
待测图像模块，用于获取实时待测图像并检测实时待测图像内的标记，得到待测标记；
[0051]
物理量结算模块，用于根据待测标记在实时待测图像中的位置以及该标记对应的物体表面位置信息，计算得到当前实时待测图像所对应的物体表面位置。
[0052]
本发明方法及系统的有益效果是：本发明通过单台相机拍摄物体表面的自带特征或制作的纹理特征，实现物体二维平移、三维平移、二维转角等参数的高精度测量。通过对制作完成后的标记进行测量，高精度得到各标记相对于待测物体的实际位置，再根据标记的实际位置(而不是理论位置)来进行物体表面的定位和测量，从而避免了标记的加工误差、安装误差等因素对测量精度的影响；且本发明无需制作和安装信息环，无需高精度加工和制作编码标志，方法更为简单经济。
附图说明
[0053]
图1是本发明一种物体表面高精度定位方法的步骤流程图；
[0054]
图2是本发明一种物体表面高精度定位系统的结构框图；
[0055]
图3是本发明具体实施例相机成像示意图；
[0056]
图4是本发明具体实施例标记制作示意图；
[0057]
图5是本发明具体实施例命名编号示意图；
[0058]
图6是本发明具体实施例物理量计算示意图。
具体实施方式
[0059]
下面结合附图和具体实施例对本发明做进一步的详细说明。对于以下实施例中的步骤编号，其仅为了便于阐述说明而设置，对步骤之间的顺序不做任何限定，实施例中的各步骤的执行顺序均可根据本领域技术人员的理解来进行适应性调整。
[0060]
本发明以单台线阵相机对圆台侧面进行扫描成像、且圆台表面缺少明显随机自然纹理为例，对本发明的实施方式进行进一步的详细说明。
[0061]
参照图1，本发明提供了一种物体表面高精度定位方法，该方法包括以下步骤：
[0062]
s1、在物体表面制作预设分布密度的标记，得到带标记的物体；
[0063]
具体地，如图4所示，通过激光雕刻技术，采用线段样式的标记，在已经制作了随机纹理的物体表面刻制线段，线段的刻制间距小于图像画幅的一半，保证在对待测物体表面成像时图像内至少有1个完整标记。
[0064]
s2、对带标记的物体进行拍摄，并编号图像中的标记，得到带编号标记的图像；
[0065]
具体地，对物体表面所有的标记成像，然后对所有图像通过阈值分割检测标记，使用质心法确定标记图像坐标；并以标记中心为参考点截取一定范围(如200
×
200像素)的图像作为该标记对应的特征图像(由于纹理信息的随机性，每个标记的特征图像都不相同)；最后采用主成分分析法(pca)区分标记的特征图像并编号，编号示意图参照图5。
[0066]
s3、根据带编号标记的图像获取物体表面各个标记的相对位置信息并建立标记编号与标记位置信息的对应关系；
[0067]
具体地，如图3所示，对于表面缺少明显随机自然纹理的圆台，首先在待测圆台表面喷涂金属漆，金属颗粒的大小根据图像像素的物理分辨率来确定，使得物体表面成像后产生散斑形状的随机图案。然后使用对圆台侧面的各个局部进行具有一定重叠度的成像，并通过图像拼接技术得到圆台一周完整的图像，以此作为基准图。建立基准图的图像像素坐标与圆台表面局部区域二维坐标的对应关系。
[0068]
可使用三坐标测量仪或者立体视觉测量的技术，直接测量物体表面各标记的坐标，确定各个标记与物体表面坐标的对应关系。
[0069]
s4、获取实时待测图像并检测实时待测图像内的标记，得到待测标记；
[0070]
s5、根据待测标记在实时待测图像中的位置以及该标记对应的物体表面位置信息，计算得到当前实时待测图像所对应的物体表面位置。
[0071]
具体地，物理量解算参照图6。
[0072]
进一步作为本方法优选实施例，所述s2、对带标记的物体进行拍摄，并编号图像中的标记，得到带编号标记的图像这一步骤，其具体包括：
[0073]
s21、对带标记的物体的表面拍摄成像，得到表面图像；
[0074]
s22、通过阈值分割对表面图像进行处理，检测得到表面图像内的所有标记的目标中心；具体地，通过阈值对图像进行二值化分割：
[0075][0076]
式中，f(x,y)为图像灰度值，x,y为像素坐标，th为阈值。然后对灰度值为255的像素点进行区域增长，并根据目标的面积确认目标，如果当前像素点周围相连的像素值为255的点个数小于阈值ts，则当前点不是目标；相反，如果像素值为255的点个数大于阈值ts，则当前点确认是目标：
[0077][0078]
在对确认的目标进行质心定位得到目标中心，质心为各个像素值为255的点的坐标中心：
[0079][0080]
其中，n为当前目标中像素值为255的点的个数，(xi,yi)为各点的图像坐标。
[0081]
s23、以目标中心为参考点，截取预设像素大小的图像作为对应标记的标记特征图像；
[0082]
具体地，以目标中心为参考点，截取200
×
200像素大小的图像作为该标记的特征图像，可以通过特征图像来区分不同的标记。
[0083]
s24、采用主成分分析法建立标记特征图像与标记编号的对应关系，将标记特征图像降至预设的10维，并取主成分矩阵前10个主成分为表示该标记的一维向量，确定唯一编号；
[0084]
具体地，采用主成分分析法建立标记特征图像与标记编号的对应关系，将200
×
200像素大小的标记特征图像降维成10个分量表示的一维向量，并以该一维向量作为标记的位移编号。
[0085]
s25、得到带编号标记的表面图像。
[0086]
进一步作为本方法优选实施例，所述s24这一步骤，其具体包括：
[0087]
s241、将200
×
200的矩阵的每一个行向量(每个变量)都减去该行向量的均值，从而使得新行向量的均值为0，得到新的数据集矩阵x；
[0088]
s242、求x的协方差矩阵，并求出协方差矩阵的特征值λ和单位特征向量e；
[0089]
s243、按照特征值从大到小的顺序，将单位特征向量排列成矩阵，得到转换矩阵p，并按px计算出主成分矩阵；
[0090]
s244、将特征图像降至特定的10维，直接取主成分矩阵前10个主成分；如果对于表示图像的目标维度无固定要求，可以用特征值计算方差贡献率和方差累计贡献率，取方差累计贡献率最大的前10个主成分，这10个主成分即为表示该标记的一维向量，以此向量作为标记的特征信息，确定唯一编号。
[0091]
进一步作为本方法优选实施例，所述根据带编号标记的图像获取物体表面各个标记的相对位置信息并建立标记编号与标记位置信息的对应关系这一步骤，其具体包括：
[0092]
s31、通过图像拼接方法将带编号标记的图像进行图像拼接，得到完整表面图像；
[0093]
s32、以完整表面图像为基准图并建立基准图的图像像素坐标与实际物体表面坐标的关系；
[0094]
s33、在基准图像上进行标记的识别和定位，并根据基准图像像素坐标与实际物体表面坐标的关系，确定各个标记与物体表面坐标的对应关系。
[0095]
具体地，假设圆台表面局部区域为近似平面，则图像平面点和物体表面局部区域上二维点的关系可表示为：
[0096][0097]
其中，a,b,c,d,e,f为描述两个平面之间变形的6个参数，(u,v)为图像像素坐标，(x,y)为圆台表面近似平面上的坐标。变形参数可以通过物体表面控制点(ui,vi)与图像对应点(xi,yi)联立如下方程组计算获取：
[0098][0099]
以基准图v方向对应于回转体的转轴高度方向，以基准图的u方向对应于回转体的转角，圆台表面点用其对应的转轴高度和转角表示为(θ,h)，图像坐标(u,v)与圆柱侧表面坐标(θ,h)的对应关系为：
[0100][0101]
其中，a,b为比例参数，可以通过相机标定获取，也可以通过物体表面控制点(θi,hi)与图像对应点(ui,vi)联立得到下列方程组计算获取：
[0102][0103]
在基准图像上进行各个标记的识别和定位，并根据基准图像像素坐标与实际物体表面坐标的关系，确定各个标记与物体表面坐标的对应关系。
[0104]
进一步作为本方法优选实施例，所述获取实时待测图像并检测实时待测图像内的标记，得到待测标记这一步骤，其具体包括：
[0105]
s41、根据变换公式确定各标记中心点在物体表面的精确位置；
[0106]
具体地，测量得到各标记中心点对应的圆台表面精确位置，如二维坐标或者三维坐标。其中，标志中心在圆台坐标系下的二维坐标(xc,yc)和三维坐标(xw,yw,zw)可以通过二维或三维坐标测量仪直接测量得到；也可以使用三维扫描仪，扫描得到标志点及圆台表面其它部位的三维点云，从而得到标志点相对于圆台坐标系的三维坐标(xw,yw,zw)。
[0107]
采用图像测量技术，可以根据成像参数确定标志中心点对应的圆台表面位置：设标志中心图像坐标为(u,v)，则其与世界坐标(xw,yw,zw)的变换公式如下：
[0108]
[0109]
如上式所示，右边第一个矩阵为相机内参数矩阵，第二个矩阵为相机外参数矩阵。相机的内参矩阵和外参矩阵可以通过张正友标定法等相机标定方法获取。使用上式中的(xw,yw,zw)代表圆台表面与标志中心对应的三维点坐标。s为中心点在摄像机坐标系中的z坐标，(u,v)为中心图像坐标，(f
x
,fy)为相机的等效焦距，(c
x
,cy)为图像主点，r
ij
(i,j＝1,2,3)为旋转矩阵r元素，ti(i＝1,2,3)为平移向量，(xw,yw,zw)代表物体表面与标记中心对应的三维点坐标。
[0110]
s42、对物体表面进行实时成像，得到实时待测图像；
[0111]
s43、对实时待测图像进行标记检测，基于向量比对确认实时待测图像内标记对应的标记编号。
[0112]
具体地，对圆台表面进行实时成像，对实时图像进行标记检测和定位(检测过程同s3)，确认当前标记对应的标记编号，并计算出实时图像中心在基准图像上的对应坐标。具体步骤如下：
[0113]
标记检测确认：根据pca降维准则对当前标记的特征图像进行降维，得到10分量的一维向量，将该向量与基准图中所有标记的一维向量分别进行比对，以基准图中相似度最大的标记编号作为当前标记编号。两个向量{ai}、{bi}的相似度c表示为：
[0114][0115]
根据标记定位，确定实时图中心在基准图上的位置：设实时图像上检测出的标记在实时图上的坐标为(ur,vr)，根据标记向量的相似度判断为基准图上的第n号标记，第n号标记(中心)在基准图上的坐标为(xn,yn)，则根据该标记得到实时图中心在基准图像上的坐标(uc,vc)为
[0116][0117]
其中，width和height分别为实时图像的像素宽度和高度。
[0118]
进一步作为本方法优选实施例，所述根据待测标记在实时待测图像中的位置以及该标记对应的物体表面位置信息，计算得到当前实时待测图像所对应的物体表面位置这一步骤，其具体包括这一步骤，其具体包括：
[0119]
s51、根据标记编号获取实时图像中心在基准图像上的对应坐标。
[0120]
s52、根据基准图像像素坐标与物体表面参数的对应关系、实时图像中心在基准图像上的坐标，确定实时图像中心对应的物体表面点；
[0121]
s53、根据实时图像中心对应的物体表面点，得到当前实时待测图像所对应的物体表面位置；
[0122]
具体地，得到实时图像中心在基准图像上的坐标(uc,vc)后，根据基准图像像素坐标与物体表面参数的对应关系，就可以确定实时图中心对应的物体表面点。
[0123]
当物体表面为平面时，实时图像中心(uc,vc)对应的物体表面的二维坐标(xc,yc)可以表示为：
[0124]
[0125]
根据物体表面的二维坐标(xc,yc)，可以直接计算实时图中心对应点的位移等参数。
[0126]
对于圆台，以基准图v方向对应于圆台的转轴高度方向，以基准图的u方向对应于圆台的转动角度，实时图中心(uc,vc)对应的物体表面点(θc,hc)表示为：
[0127][0128]
所述当前实时待测图像所对应的物体表面位置包括一维坐标、二维坐标、三维坐标、一维转角和二维转角。
[0129]
根据实时图中心(uc,vc)对应的物体表面点(θc,hc)，可以计算一维角度和二维角度的变化：
[0130]
图像位置与物体一维转角的对应关系：设圆台周长为l，圆台轴线方向为x方向，圆台周线方向为y方向，设图像中心的像素坐标为(u,v)，其对应的物体表面二维坐标为(xc,yc)，则该点对应的旋转角度为：
[0131][0132]
上式中，θ为待测旋转体的旋转角度，y0为旋转体外表面展开的周线原点，yc为旋转体外表面展开的标识中心坐标点，l为旋转体外表面展开的周线长。
[0133]
图像位置与球面二维转角(经纬度)的对应关系：图像中心坐标为(u,v)对应在球面上的世界坐标为(xw,yw,zw)，假设世界坐标系是以球心为原点，则经纬度为：
[0134][0135]
上式中，图像中心坐标为(u,v)对应在待测球体表面的世界坐标记为(xw,yw,zw)，α为待测球体的经度，β为待测球体的纬度。
[0136]
一维角度变化：两次实时成像时图像中心对应的物体一维角度变化计算如下：
[0137]
δθ＝θ
c1-θ
c2
＝au
c1-au
c2
[0138]
二维角度变化：两次实时成像时图像中心对应的球体表面二维角度变化(假设原点在球心，且球半径为r)计算如下：
[0139][0140]
本发明具有如下优势：
[0141]
1、相对于光栅技术和成像式编码器技术依据刻度和标识的理论位置进行测量，本发明的测量依据是待测物体表面标记的真实位置。本发明通过对制作完成后的随机纹理信息和标记进行测量，高精度得到各标记相对于待测物体的实际位置，再根据标记的实际位置(而不是理论位置)来进行物体表面的定位和测量，从而避免了标记的加工误差、安装误差等因素对测量精度的影响。这样测量不仅精度高、可靠性好，而且大幅降低了对加工工艺、器件材质的要求。
[0142]
2、相对于信息环识别的技术，本发明可以直接在物体表明制作随机散斑信息和标记(喷涂、粘贴)，无需制作和安装信息环，无需高精度加工和制作编码标志，方法更为简单经济。
[0143]
3、相对于信息环只能针对圆柱体测量，通过测量圆柱侧面一维位移进行转角测量。本发明可以直接应用于圆台(侧面与上下底面均可)、圆锥、圆盘以及任意不规则回转体的表面位移测量与转角测量。
[0144]
4、相对于信息环识别和双像机的测量技术，其采用的遍历搜索式图像定位方法效率低下，很难用于实时快速测量。本发明提出了直接检测和识别标识以进行图像定位的方法，无需在图像定位的时候对基准图像的所有像素进行遍历搜索，具有执行效率高、实时性好、简单可靠的优点。
[0145]
5、相对于成像式编码器和信息环标识采用的单一标识编码测量方式，本发明采用无编码的简单图案(如点、线、十字丝等)，也可以兼容编码标识，这样不仅使标识的加工制作更为简单方便，彻底摆脱了高精度加工要求，还可以使像机成像的范围更小，从而获得更高测量精度。
[0146]
6、相对于双像机的测量方案，本发明只需要采用单像机，即可实现对物体表面标记点的识别定位，使得整个测量过程更为简单高效。
[0147]
7、相对于现有的只能进行一维转角和一维位移高精度测量的技术，本发明在此基础上还可以用于二维平移、三维平移、二维转角等参数的测量。
[0148]
如图2所示，一种物体表面高精度定位系统，包括：
[0149]
标记制作模块，用于在物体表面制作预设分布密度的标记，得到带标记的物体；
[0150]
编号模块，用于对带标记的物体进行拍摄，并编号图像中的标记，得到带编号标记的图像；
[0151]
标记位置关系模块，用于根据带编号标记的图像获取物体表面各个标记的相对位置信息并建立标记编号与标记位置信息的对应关系；
[0152]
待测图像模块，用于获取实时待测图像并检测实时待测图像内的标记，得到待测标记；
[0153]
物理量解算模块，用于根据待测标记在实时待测图像中的位置以及该标记对应的物体表面位置信息，计算得到当前实时待测图像所对应的物体表面位置。
[0154]
上述方法实施例中的内容均适用于本系统实施例中，本系统实施例所具体实现的功能与上述方法实施例相同，并且达到的有益效果与上述方法实施例所达到的有益效果也相同。
[0155]
以上是对本发明的较佳实施进行了具体说明，但本发明创造并不限于所述实施例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替换，这些等同的变形或替换均包含在本技术权利要求所限定的范围内。