专利摘要

本发明公开了一种基于虚拟立体视觉的非合作目标协同测量方法，包括以下步骤：获取非合作翻滚目标的三维点云图；提取出卫星的喷嘴中心与三角支架内切圆中心的相对位姿关系；对获取的卫星喷嘴与三角支架的图像进行图像预处理操作；通过坐标转换获得在特定坐标系下的位姿关系；求得所述卫星喷嘴与三角支架内切圆的半径，以得到喷嘴与三角支架中心的位姿信息。以及一种基于虚拟立体视觉的非合作目标协同测量系统。其降低了测量系统的硬件数量，扩展了测量系统的能力，具有较高的计算效率，减少了传感器的数量，降低了空间发射的成本，且拓展了空间测量的能力，广泛应用于空间非合作卫星视觉测量的技术领域。

权利要求

1.一种基于虚拟立体视觉的非合作目标协同测量方法，其特征在于，包括以下步骤：

获取非合作翻滚目标的三维点云图；

提取出所述三维点云图中卫星的喷嘴中心与三角支架内切圆中心的相对位姿关系；

对获取的卫星喷嘴与三角支架的图像进行图像预处理操作；

对所述卫星喷嘴与三角支架的图像分别进行处理，通过空间圆形目标成像的坐标变换方法，分别求得所述喷嘴中心与三角支架内切圆中心在各自相机坐标系下的位姿关系；

具体包括：将所述三角支架内切圆中心的位姿转换到左相机坐标系下，组成封闭解，求得卫星喷嘴与三角支架内切圆的半径，从而求得喷嘴与三角支架中心的位姿信息；

利用最小二乘法求得所述卫星喷嘴与三角支架内切圆的半径，以得到喷嘴与三角支架中心的位姿信息。

2.根据权利要求1所述的基于虚拟立体视觉的非合作目标协同测量方法，其特征在于，通过采用双目相机来获取所述三维点云图。

3.根据权利要求1所述的基于虚拟立体视觉的非合作目标协同测量方法，其特征在于，所述方法包括：通过Hough检测获取到卫星的喷嘴中心及三角支架。

4.根据权利要求1至3任一项所述的基于虚拟立体视觉的非合作目标协同测量方法，其特征在于，所述方法还包括：对获取的卫星喷嘴与三角支架的图像进行图像增强处理，并对所述图像进行滤波和边缘检测。

5.根据权利要求4所述的基于虚拟立体视觉的非合作目标协同测量方法，其特征在于，所述方法包括：通过分别设置于机器人左、右机械臂的单目相机来分别采集卫星喷嘴与三角支架的图像。

6.根据权利要求5所述的基于虚拟立体视觉的非合作目标协同测量方法，其特征在于，所述方法还包括：对左臂手眼相机采集的图像进行Hough椭圆检测，得到表示椭圆的相关参数；对右臂手眼相机进行Hough直线检测，并提取三条直线的交点，从而获取三角支架的边缘特征。

7.根据权利要求6所述的基于虚拟立体视觉的非合作目标协同测量方法，其特征在于，所述方法还包括：获取所述三角支架的内切椭圆的相关参数。

8.一种基于虚拟立体视觉的非合作目标协同测量系统，其特征在于，包括：图像采集模块，用于执行步骤获取非合作翻滚目标的三维点云图；

图像处理模块，用于执行步骤提取出所述三维点云图中卫星的喷嘴中心与三角支架内切圆中心的相对位姿关系；

对获取的卫星喷嘴与三角支架的图像进行图像预处理操作；

对所述卫星喷嘴与三角支架的图像分别进行处理，通过坐标转换获得在特定坐标系下的位姿关系；

位姿求解模块，用于执行步骤利用最小二乘法求得所述卫星喷嘴与三角支架内切圆的半径，以得到喷嘴与三角支架中心的位姿信息。

说明书

技术领域

本发明涉及空间非合作失效卫星视觉测量的技术领域，具体为基于非合作目标协同测量方法及系统。

背景技术

随着空间探索的不断深入，在轨航天器数量越来越多，除了少部分仍然正常工作外，大部分已到寿命末期而被废弃，且部分因为出现故障而无法正常工作，使得轨道垃圾不断累积。利用空间机器人进行卫星维修、轨道垃圾清理成为近年来热门热点的研究方向。然而，作为在轨维护的非合作目标，大都处于翻滚状态，动力学特性极其复杂。

为了实现在轨服务，首先要解决的技术难点就是非合作目标的视觉测量问题。由于失效目标上没有任何的合作标识，且目标又处于翻滚的状态，这对非合作目标的视觉测量造成了很大的影响。以往的视觉测量研究中，如美国的“通用航天器在轨修正系统”以及日本的ETS-VII项目等主要针对的是合作目标的视觉测量。目前的研究，针对非合作目标的视觉测量机械臂末端都需要携带双目相机，在双臂协同抓捕非合作目标前，传统的研究需要两对双目相机进行双臂的独立测量(例如圆形的卫星喷嘴，三角形的帆板支架等)，很难满足空间视觉测量的实时性要求，同时图像数据量很大，对星载系统的存储能力和计算能力提出了更高的要求，不利于空间机器人的视觉测量，且为双臂的协同捕获带来诸多不利因素。

为了减少双臂测量过程中视觉传感器的数量，降低图像存储量，提高视觉测量算法的效率，在保证测量精度的前提下有必要开展提高视觉计算速度的同时又能降低空间测量的成本、功耗与质量等相关视觉测量技术的研究。

发明内容

本发明旨在至少在一定程度上解决相关技术中的技术问题之一。为此，本发明的一个目的是提供一种提高视觉计算速度、节省成本的基于非合作目标协同测量方法及系统。

本发明所采用的技术方案是：

本发明提供一种基于虚拟立体视觉的非合作目标协同测量方法，包括以下步骤：

获取非合作翻滚目标的三维点云图；

提取出所述三维点云图中卫星的喷嘴中心与三角支架内切圆中心的相对位姿关系；

对获取的卫星喷嘴与三角支架的图像进行图像预处理操作；

对所述卫星喷嘴与三角支架的图像分别进行处理，通过坐标转换获得在特定坐标系下的位姿关系；

利用最小二乘法求得所述卫星喷嘴与三角支架内切圆的半径，以得到喷嘴与三角支架中心的位姿信息。

作为该技术方案的改进，通过采用双目相机来获取所述三维点云图。

作为该技术方案的改进，所述方法包括：通过Hough检测获取到卫星的喷嘴中心及三角支架。

作为该技术方案的改进，所述方法还包括：对获取的卫星喷嘴与三角支架的图像进行图像增强处理，并对所述图像进行滤波和边缘检测。

进一步地，所述方法包括：通过分别设置于机器人左、右机械臂的单目相机来分别采集卫星喷嘴与三角支架的图像。

进一步地，所述方法还包括：对左臂手眼相机采集的图像进行Hough椭圆检测，得到表示椭圆的相关参数；对右臂手眼相机进行Hough直线检测，并提取三条直线的交点，从而获取三角支架的边缘特征。

进一步地，所述方法还包括：获取所述三角支架的内切椭圆的相关参数。

进一步地，通过空间圆形目标成像的坐标变换方法，分别求得所述喷嘴中心与三角支架内切圆中心在各自相机坐标系下的位姿关系。

进一步地，所述方法还包括：将所述三角支架内切圆中心的位姿转换到左相机坐标系下，组成封闭解，求得卫星喷嘴与三角支架内切圆的半径，从而求得喷嘴与三角支架中心的位姿信息。

另一方面，本发明还提供一种基于虚拟立体视觉的非合作目标协同测量系统，包括：

图像采集模块，用于执行步骤获取非合作翻滚目标的三维点云图；

图像处理模块，用于执行步骤提取出所述三维点云图中卫星的喷嘴中心与三角支架内切圆中心的相对位姿关系；

对获取的卫星喷嘴与三角支架的图像进行图像预处理操作；

对所述卫星喷嘴与三角支架的图像分别进行处理，通过坐标转换获得在特定坐标系下的位姿关系；

位姿求解模块，用于执行步骤利用最小二乘法求得所述卫星喷嘴与三角支架内切圆的半径，以得到喷嘴与三角支架中心的位姿信息。

本发明的有益效果是：

本发明提供的基于非合作目标协同测量方法、系统，基于多视点观测重构出翻滚目标的三维表面信息，进一步地提取双臂抓捕手柄的相对位姿关系，作为后续手眼视觉测量的先验知识，圆形的卫星喷嘴和三角形的帆板支架常作为空间抓捕的对象，根据双臂末端携带的单目手眼相机分别独立地观测各自的抓捕位置，通过双目的协调测量实现喷嘴和三角支架的位姿测量，实现了双臂空间机器人系统采用每个臂携带一个相机并观测不同对象下(每个臂“单臂-单目”配置)的双臂协同测量方法，降低了测量系统的硬件数量，扩展了测量系统的能力。通过求解多特征的几何约束方程，从而确定位姿信息，具有较高的计算效率，减少了传感器的数量，降低了空间发射的成本，且拓展了空间测量的能力，广泛应用于空间非合作卫星视觉测量的技术领域。

附图说明

图1是本发明一实施例的虚拟立体视觉测量模型图；

图2是本发明一实施例的空间非合作目标的虚拟立体视觉协同测量的结构模型示意图；

图3是本发明一实施例的基于虚拟立体视觉的非合作目标协同测量流程图；

图4是本发明一实施例的多视点三维重建示意图；

图5是本发明一实施例的三角支架像平面投影关系示意图；

图6是本发明一实施例的三角支架内切圆成像示意图；

图7是本发明一实施例的虚拟立体视觉位姿解算结构图。

具体实施方式

需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。

本发明提供一种基于虚拟立体视觉的非合作目标协同测量方法，包括以下步骤：

获取非合作翻滚目标的三维点云图；

提取出所述三维点云图中卫星的喷嘴中心与三角支架内切圆中心的相对位姿关系；

对获取的卫星喷嘴与三角支架的图像进行图像预处理操作；

对所述卫星喷嘴与三角支架的图像分别进行处理，通过坐标转换获得在特定坐标系下的位姿关系；

利用最小二乘法求得所述卫星喷嘴与三角支架内切圆的半径，以得到喷嘴与三角支架中心的位姿信息。

作为该技术方案的改进，通过采用双目相机来获取所述三维点云图。

作为该技术方案的改进，所述方法包括：通过Hough检测获取到卫星的喷嘴中心及三角支架。

作为该技术方案的改进，所述方法还包括：对获取的卫星喷嘴与三角支架的图像进行图像增强处理，并对所述图像进行滤波和边缘检测。

进一步地，所述方法包括：通过分别设置于机器人左、右机械臂的单目相机来分别采集卫星喷嘴与三角支架的图像。

进一步地，所述方法还包括：对左臂手眼相机采集的图像进行Hough椭圆检测，得到表示椭圆的相关参数；对右臂手眼相机进行Hough直线检测，并提取三条直线的交点，从而获取三角支架的边缘特征。

进一步地，所述方法还包括：获取所述三角支架的内切椭圆的相关参数。

进一步地，通过空间圆形目标成像的坐标变换方法，分别求得所述喷嘴中心与三角支架内切圆中心在各自相机坐标系下的位姿关系。

进一步地，所述方法还包括：将所述三角支架内切圆中心的位姿转换到左相机坐标系下，组成封闭解，求得卫星喷嘴与三角支架内切圆的半径，从而求得喷嘴与三角支架中心的位姿信息。

参照图1，是本发明一实施例的虚拟立体视觉测量模型图。通过对空间非合作目标的三维重构实现卫星喷嘴中心与三角支架内切圆中心的相对位姿关系，这个位姿关系作为后续协同测量的先验知识。

参照图2，是本发明一实施例的空间非合作目标的虚拟立体视觉协同测量结构示意图。该服务卫星由基座和搭载基座上两个串联空间机械臂组成。在机械臂末端分别安装有测量目标卫星位姿信息的单目手眼相机，为了便于讨论定义了如下符号：

Ot-XtYtZt：为空间目标的轨道坐标系；

Obt-XbtYbtZbt：为非合作目标的质心坐标系；

Obc-XbcYbcZbc：为空间机器人基座质心坐标系；

Oea-XeaYeaZea：为左臂末端坐标系；

Oeb-XebYebZeb：为右臂末端坐标系；

OL-XLYLZL：为左臂手眼相机坐标系；

OR-XRYRZR：为右臂手眼相机坐标系；

分别表示喷嘴圆心在左臂手眼相机坐标系下的位置矢量与法向量；

分别表示三角支架内切圆圆心在右臂手眼相机坐标系下的位置矢量与法向量。

参照图3，是本发明一实施例的基于虚拟立体视觉的空间非合作目标协同测量方法的流程图。包括步骤如下：

S1：根据安装在基座上的双目相机测量，得到非合作翻滚目标的三维点云图。

S2：根据S1中的点云图通过Hough检测找出圆形卫星喷嘴以及三角形的帆板支架，提取喷嘴中心与三角形支架内切圆中心的相对位姿关系

S3：对左右臂的单目相机获取的图像进行图像增强，再对该图像进行滤波和边缘检测等图像预处理操作。

S4：由S3得到的左臂手眼相机图像进行Hough椭圆检测，得到表示椭圆的5个参数，对右臂手眼相机进行Hough直线检测，并提取三条直线的交点P1,P2,P3，从而获取三角支架的边缘特征。

S5：根据S4中得到的三角形特征提取表示该三角形的内切椭圆的5个参数(xc,yc,ac,bc,θ)。

S6：由S5与S4中得到的两个椭圆特征的5个参数，通过空间圆形目标成像的坐标变换分别计算圆心在各自相机坐标系下的位姿表示，记为 r)。

S7:由S2得到 将S6中得到的三角形内切圆中心位姿转换到左相机坐标系下，组成封闭解，利用最小二乘法求得卫星喷嘴与三角支架内切圆的半径，从而求得喷嘴与三角型支架中心的位姿信息。

参照图4，作为本发明的一个实施例的多视点三维重建示意图。将点云邻域点的平均距离、曲率和法向夹角来表达点云模型的特征。

首先，将待匹配的两帧点云记为PA和PB，计算两组点云的特征参数f(pi)，判断所选特征点是否满足f(pi)＞τf(τf为特征参数的阀值)。

设点 点 如果满足 和 则点 是点 的对应点。

其次，把对应点对之间的欧氏距离作为目标的优化指标，则配准的目标函数可定义为：

采用单位四元数法来求解变换矩阵，根据式(1)可以求得旋转矩阵ARB，于是，平移向量AtB可表示为：

其中， 分别为点集PA，PB的均值。

进一步地，将几何特征引入到误差函数可得：

重复步骤(2)与(3)，直至迭代误差d＜τ或迭代次数超过阈值才退出。

进一步地，将求得的位姿变换矩阵为依据将目标点云转换到参考点云所在的坐标系中，实现对待配卫星点云的配准。

卫星三维重建后的三维无序点云只包含三维坐标信息，无任何拓扑信息，它的边界线是由明显曲率特征的特征点组成。假设三维点云的数据集为{P}，对于任意点Pi的邻域内n个点，如果分布在一侧，则可认为点Pi为边界特征点。反之，认为Pi是内部点。基于此可以利用数据点Pi及其周围点的分布来判断是否为边界点。令任意点Pi的坐标为(xi,yi,zi)T，当找出了所有的边界点，Pi邻域内点形心的坐标可表示为

假设Pi邻域内点拟合的平面方程π为：

a1x+a2y+a3z+a4＝0(5)

以Pi邻域内n个点到π距离的平方和最小建立目标函数：

为了使f(a1,a2,a3,a4)最小，令 可以得到：

a4＝-(a1xc+a2yc+a3zc)(7)

于是，式(6)可进一步化简为

由于 根据(8)可以建立拉格朗日优化方程：

令

对式(8)求偏导可以得到：

(J-κI33)·[a1,a2,a3]T＝0 (10)

根据式(10)可以得到：

于是，最小特征值λmin(J)对应矩阵J的特征向量就是所求切平面的法向量 根据求得的形心Pcenter和法向量(a1,a2,a3)T即可确定切平面的方程。设Pi点在所求切平面上投影点Pi'的坐标为 由于该投影点在过点Pi且方向向量为 的空间直线上，即有：

由式(12)可解出投影点的坐标为：

提取到切平面π之后，把点集{Pi}向切平面π投影，得到二维投影点集{Pi'}。当提取到了二维离散边界曲线后，就可以根据霍夫变换获得二维点集的椭圆特征与三角支架特征，进而获得喷嘴圆心与三角支架内切圆圆心的位置矢量以及法向量。

参照图5，是本发明实施例的等腰三角形像平面投影关系示意图。

根据等腰三角形内切圆在像平面的仿射变换可知，原三角支架的内切圆圆心经过仿射变换在像平面内变成了椭圆的中心，切点PE的坐标可表示为：

假设平面椭圆的方程为：

Au2+Buv+Cv2+Du+Ev+F＝0(15)

式中，A,B,C,D,E,F为椭圆方程二次型的6个参数。

进一步地，椭圆二次型可以表示为：

联立式(16)与(15)，椭圆的几何形式的5个参数可以表示为：

显然，点O'在AE两点组成的直线上，即有：

(vE-v1)uc-(uE-u1)vc＝(vE-v1)u1-(uE-u1)v1 (18)

在像平面内，过椭圆上任意点(up,vp)的切线方程可以表示为：

2Auup+B(upv+uvp)+2Cvvp+D(u+up)+E(v+vp)＝0 (19)

由于P2'与P3'在过PD点的椭圆切线上，于是有：

根据方程组(20)，PD点的像素坐标可表示为:

同理，PE与PF点的像素坐标可分别表示为:

根据式(14)与式(22)，有如下关系：

另一方面，三个切点(k＝D,E,F)都在椭圆上，即有

联立式(24)、(25)以及式(18)即可通过数值法求得椭圆方程的6个参数。

参照图6，是本发明实施例的三角支架的位姿等效测量示意图。三角形卫星支架的内切圆经过相机的仿射变换变成了相平面里面的内切椭圆。

三角支架内切圆在像平面内呈现的结果是椭圆，根据求解得到的椭圆六个参数A,B,C,D,E,F，结合公式(17)可以得到像平面上的椭圆参数(uc,vc,ac,bc,θ)。

把椭圆平面转换到归一化焦平面上有：

式中，fu,fv,u0,v0为相机标定参数。

将相机模型(26)代入(16)可得空间斜圆锥在相机坐标系下的方程：

式中

令Xc＝[xc,yc,zc]T，则斜圆锥方程(27)又可以表示为：

式中，

假设矩阵A的特征值为λ1(A),λ2(A),λ3(A)，对应的规范化特征向量为 把P写成列向量的形式[r1r2r3]，不失一般性，令λ1(A)≥λ2(A)＞0，则有

参照图7，是本发明实施例的虚拟立体视觉位姿解算结构图。 与 分别表示左、右相机的图像坐标系， 表示右相机相对于左相机坐标系的齐次变换矩阵，所测量的空间圆平面分别为 与

根据双目标定结果，右抓捕圆面的位姿在左相机坐标系下的表示为：

根据远距离多视点协同测量的结果，左右抓捕圆面的圆心位置以及法向量之间的相对关系可以用齐次变换矩阵 来表示，于是右抓捕面相对于左抓捕面的位置和姿态可表示为：

比较公式(31)与公式(33)可得：

式中

对于公式(35)的求解可以采用最小二乘法，于是

式中，

于是，喷嘴中心的位姿为：

同理，三角支架内切圆中心的位姿为：

由于每次测量都可以得到两组解，为了获得目标抓捕面的法向量，可以通过比较左右手眼相机在左相机坐标系下两个法向量的夹角来衡量，取它们之间夹角最小的那组解作为最终视觉测量的位姿解。

另一方面，本发明还提供一种基于虚拟立体视觉的非合作目标协同测量系统，包括：

图像采集模块，用于执行步骤获取非合作翻滚目标的三维点云图；

图像处理模块，用于执行步骤提取出所述三维点云图中卫星的喷嘴中心与三角支架内切圆中心的相对位姿关系；

对获取的卫星喷嘴与三角支架的图像进行图像预处理操作；

对所述卫星喷嘴与三角支架的图像分别进行处理，通过坐标转换获得在特定坐标系下的位姿关系；

位姿求解模块，用于执行步骤利用最小二乘法求得所述卫星喷嘴与三角支架内切圆的半径，以得到喷嘴与三角支架中心的位姿信息。

再一方面，本发明还提供一种基于虚拟立体视觉的非合作目标协同测量设备，包括：

至少一个处理器；以及，

与所述至少一个处理器通信连接的存储器；其中，

所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行如所述的方法。

本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机可执行指令，所述计算机可执行指令用于使计算机执行如所述的方法。

本发明实现了双臂空间机器人系统采用每个臂携带一个相机并观测不同对象下(每个臂“单臂-单目”配置)的双臂协同测量方法，降低了测量系统的硬件数量，扩展了测量系统的能力。通过求解多特征的几何约束方程，从而确定位姿信息，减少了传感器数量和计算复杂度，降低了空间测量的成本、功耗与质量，对于太空应用而言具有极其重要的意义。

以上是对本发明的较佳实施进行了具体说明，但本发明创造并不限于所述实施例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做出种种的等同变形或替换，这些等同的变形或替换均包含在本申请权利要求所限定的范围内。

基于非合作目标协同测量方法及系统专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。