CN114003029B - 一种声光融合导引的自主水下航行器鲁棒对接回收方法 - Google Patents
一种声光融合导引的自主水下航行器鲁棒对接回收方法 Download PDFInfo
- Publication number
- CN114003029B CN114003029B CN202111065374.8A CN202111065374A CN114003029B CN 114003029 B CN114003029 B CN 114003029B CN 202111065374 A CN202111065374 A CN 202111065374A CN 114003029 B CN114003029 B CN 114003029B
- Authority
- CN
- China
- Prior art keywords
- recovery
- auv
- recovery cage
- cage
- docking
- 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
- 238000011084 recovery Methods 0.000 title claims abstract description 170
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000003032 molecular docking Methods 0.000 title claims abstract description 61
- 230000004927 fusion Effects 0.000 title claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 230000008569 process Effects 0.000 claims abstract description 17
- 230000003287 optical effect Effects 0.000 claims abstract description 16
- 238000005259 measurement Methods 0.000 claims abstract description 12
- 238000012937 correction Methods 0.000 claims abstract description 3
- 230000001105 regulatory effect Effects 0.000 claims abstract description 3
- 238000004064 recycling Methods 0.000 claims description 16
- 238000007667 floating Methods 0.000 claims description 8
- 238000013459 approach Methods 0.000 claims description 6
- 230000033001 locomotion Effects 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 claims description 6
- 230000000877 morphologic effect Effects 0.000 claims description 4
- 238000003708 edge detection Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims 1
- 238000005265 energy consumption Methods 0.000 abstract description 5
- 238000012545 processing Methods 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 210000001503 joint Anatomy 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 235000002566 Capsicum Nutrition 0.000 description 3
- 239000006002 Pepper Substances 0.000 description 3
- 241000722363 Piper Species 0.000 description 3
- 235000016761 Piper aduncum Nutrition 0.000 description 3
- 235000017804 Piper guineense Nutrition 0.000 description 3
- 235000008184 Piper nigrum Nutrition 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 150000003839 salts Chemical class 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/02—Control of position or course in two dimensions
- G05D1/0206—Control of position or course in two dimensions specially adapted to water vehicles
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提出一种声光融合导引的自主水下航行器鲁棒对接回收方法，AUV以航路点导引方式依次完成航路点，并到达位于回收笼正前方的最后一个航路点，根据参考航向角以及水下航行器速度，采用模糊控制直线跟踪方法继续行驶，进行水下航行器自适应导引对接回收；并在接近回收笼的过程中，采用声光导引组合，根据超短基线USBL测量结果以及光学手段，结合航行器自身的经纬度，对回收笼的经纬度进行实时更新；当AUV调整到回收笼中心附近时，AUV相机识别到回收笼底部中心目标光源，AUV对回收笼底部中心目标光源进行自适应跟踪与位姿计算结果修正，完成水下航行器的自主对接。本发明能够实现AUV的水下回收，降低AUV能源消耗，同时避免发生不必要的碰撞，提高对接成功率。
Description
技术领域
本发明涉及的是一种自主水下航行器的控制方法，更具体地说是在声光融合导引下的自主水下航行器鲁棒对接回收方法。
背景技术
在过去的几十年里，机器人技术在水下环境中的应用迅速发展，AUV可以在水下进行搜索、探测和侦查等任务。在完成作业任务后都需要能够实现回收处理，人为的吊装回收不仅成本巨大，而且十分的危险，目前AUV远距离通过声学设备导引，近距离通过视觉传感器精准识别水下航行器与对接装置的相对位姿，从而实现AUV自主对接的框架得到了国内外学者的广泛追捧。
对接回收问题是自主水下航行器研究过程中的关键课题，由于被控对象的非线性、耦合性和复杂环境干扰等特性，自主水下航行器对接回收问题的复杂度大大提升。目前为止，国内外先后开展了关于水下航行器自主回收的研究，其中研究的关键就在于如何使水下航行器能够精准的与回收装置实现对接。初期的研究主要是依靠声学设备比如长基线、短基线、超短基线，进行导航定位，但是由于声信号数据更新频率较低，在近距离精准对接的过程中稳定性和精度较差，对接成功率较低。随着计算机视觉的快速发展，视觉传感器具有高帧率、低噪声、延迟低等优点，视觉传感器在水下静距离对接任务中有着独特的优势。
但目前公开的声光融合导航方法还存在控制复杂，AUV能源消耗较大等问题，需要结合工程实际，确定合适的控制逻辑，减少控制复杂程度，降低AUV能源消耗。
发明内容
为解决现有技术存在的问题，本发明提出一种声光融合导引的自主水下航行器鲁棒对接回收方法，既通过声学导航克服了惯性导航存在累计误差的缺点，也解决了单靠声学导航数据更新慢，对接成功率较低的缺点，而且结合工程实际经验，利用测量距离设定声光导航的加权值，控制逻辑简单，在满足系统较高的对接成功率、安全性和对接可靠性的情况下，降低AUV能源消耗。
本发明的技术方案为：
所述一种声光融合导引的自主水下航行器鲁棒对接回收方法，包括以下步骤：
步骤1：按照任务需求设定航路点，并根据回收笼位置及回收方向，将最后一个航路点设置在回收笼正前方；水下航行器以航路点导引方式依次完成航路点，并到达位于回收笼正前方的最后一个航路点；
步骤2：执行回收任务参考深度为回收笼深度，跟踪的航向为回收笼口的朝向，跟踪的终点为回收笼后方；根据参考航向角ψreference以及AUV在大地坐标系下的速度，采用模糊控制直线跟踪方法继续行驶，进行AUV自适应导引对接回收；当AUV与回收笼距离小于设定阈值后未检测到超短基线USBL所需声源，则结束当前对接回收任务，并判断是否达到最大尝试次数，若未达到最大尝试次数则重新回到最后一个航路点执行新一次对接任务，若达到最大尝试次数，则上浮等待回收；
步骤3：在AUV接近回收笼的过程中，采用声光导引组合，根据超短基线USBL 测量结果以及光学手段，结合航行器自身的经纬度loncur和latcur，对回收笼经纬度lond和 latd进行更新；
步骤4：当AUV调整到回收笼中心附近时，AUV相机识别到回收笼底部中心目标光源，AUV对回收笼底部中心目标光源进行自适应跟踪与位姿计算结果修正，完成水下航行器的自主对接；
若AUV未检测到回收笼底部中心目标光源，则结束当前对接回收任务，并判断是否达到最大尝试次数，若未达到最大尝试次数则重新回到最后一个航路点执行新一次对接任务，若达到最大尝试次数，则上浮等待回收。
进一步的，步骤2中，参考航向角ψreference通过以下过程得到：
参考航向角ψreference的计算方式为：
ψreference＝ψcross+ψinit
其中ψcross为直线跟踪角，ψinit为任务初始设定的参考航向角；
ψcross＝arctan(k*L)
k为根据AUV与回收笼距离而确定的直线跟踪参数；L为航行器与回收笼之间的横向距离。
进一步的，步骤2中，在水下航行器距离回收笼100米以外时，选择直线跟踪参数k＝0.03，在水下航行器距离回收站100米及100米以内时，选择直线跟踪参数k＝0.08。
进一步的，步骤2中，L为航行器与回收笼之间的横向距离：
L＝-dE1*cos(ψ)-dN1*sin(ψ)
其中北向距离dN1和东向距离dE1的近似值为
dN1＝(latd-latcur)*R
式中R为地球子午圈的半径。
进一步的，步骤3中，建立声光导引的计算模型：
I(x)＝J(x)f(x)+(1-J(x))g(x)
其中I(x)为声光导引组合下的回收笼经纬度计算值，J(x)为远距离声学系数，取值为1或者0，当检测到回收笼口部位置的光信号个数不足4个时取值为1，反之为0， f(x)为基于USBL的回收笼经纬度测量结果，g(x)为基于光学信号的回收笼经纬度测量结果；
利用上一周期更新后的回收笼经纬度lond和latd计算回收笼与AUV之间的实时距离
在Dis>LL时，AUV沿着直线驶向回收笼的过程中，根据超短基线USBL测量结果对回收笼经纬度lond和latd进行更新：LL为设定距离阈值；
在Dis≤LL时，控制水下航行器速度降低，并增加采用光学手段对回收笼的经纬度进行更新。
进一步的，根据超短基线USBL测量结果对回收笼经纬度lond和latd进行更新过程为：
USBL提供回收笼相对于AUV之间的东向距离dE2和北向距离dN2，得到更新后的回收笼经纬度lond和latd为：
lond＝loncur+dE2*ρ2
latd＝latcur+dN2*ρ1
其中纬度更新率ρ1和经度更新率ρ2根据大地测量学估算。
进一步的，在回收笼口部位置具有至少4个光源，在回收笼底部中心具有一个光源，回收笼口部位置光源颜色与回收笼底部中心位置光源颜色不同；
当Dis≤LL时，若AUV相机识别到回收笼口部位置的光源数小于4个，则航行器仍继续使用USBL获得的信息来进行回收笼经纬度更新，同时打开辅推加大左右横移量寻找剩下的蓝色光源，并且把速度再次降低；
如果当Dis<LT后，AUV相机识别到回收笼口部位置的光源数仍小于4个，则结束当前对接回收任务，并判断是否达到最大尝试次数，若未达到最大尝试次数则重新回到最后一个航路点执行新一次对接任务，若达到最大尝试次数，则上浮等待回收；
如果AUV相机识别到回收笼口部位置的光源数达到4个，则将USBL加权降为零，只利用图像解算结果采用直线跟踪导引算法使AUV朝回收笼中心靠近。
进一步的，LL取为25m，LT取15m。
进一步的，利用图像解算结果采用直线跟踪导引算法使AUV朝回收笼中心靠近过程中，对采集的图像采用最大类间方差法进行自适应阈值分割，采用形态学梯度边缘检测方法识别目标光源的轮廓形状，并采用EPnP算法进行位姿估计，得到AUV与回收笼口中心位置的偏差量，包括横向偏差量、纵向偏差量和轴向偏差量；利用辅推控制，将横向偏差量、纵向偏差量调整至接近0，同时将速度降低至能够保持航行状态的最低速度；当轴向偏差量deltaz小于1500mm时，且横向纵向误差满足需求后，不再调整航行器位姿，并且提高主推转速，维持当前航向前进。
有益效果
本发明所构思的技术方案与现有技术相比，主要具备以下的技术优点如下：
1、本发明使用声光融合导引，采用直线跟踪算法对路径进行规划设计，先采用USBL声学导引对回收对接装置的位置进行更新，不断优化提高测量精度，保证AUV 能安全可靠地靠近回收装置；
2、本发明建立了声光导引的计算模型，结合了声学导引和光学导引的优点，根据情况切换导引方式，大大提高了对接的成功概率，并且在达到条件时才打开光学检测传感器，节约了能源消耗；
3、声光导引过程中，回收装置的位置在不断优化，该方法不仅适用于静态回收装置，也适用于动态回收装置。
4、针对识别到不足四个蓝色光源的情况，利用已有的蓝色光源位置和USBL声学导引信息依靠辅推调整左右偏移量，直至找到四个蓝色光源，解决了笼口方向发生移动以及动态笼口不确定方向的问题。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1回收对接总体示意图；
图2单次对接回收流程图；
图3目标光源安装示意图；
图4图像处理算法流程图；
图5声光融合下的加权系数变化情况；
图6阈值分割结果图；
图7闭运算结果图；
图8 EPnP虚拟控制点模型；
图9航行器速度变化曲线；
图10航行器深度变化曲线；
图11相对角度偏差；
图12相对位置偏差。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图和实施例，对本发明进行进一步详细说明。
步骤1：按照期望的路径设定航路点，航路点的设置应保证避开障碍物，而为了实现水下航行器自适应导引对接回收，设置的离回收笼最近的最后一个航路点在回收笼正前方，本实施例中处于回收笼正前方700米处。
在航路点导引航行器的过程中，航行器通过自身携带的惯性导航系统设备输出解算的自身经度loncur、纬度latcur(本实施例中经纬度均指大地坐标系下的经纬度)和航向角ψ(本实施例中航向角以北偏西为正)，回收笼的经度lond和纬度latd由GPS测得，以此来完成对航路点的跟踪。
步骤2：当水下航行器以航路点导引方式依次完成航路点，并到达离回收笼最近的最后一个航路点后，执行回收任务，参考深度为回收笼深度，跟踪的航向为回收笼口的朝向，跟踪的终点为回收笼(lond,latd)的后方；此时根据参考航向角ψreference以及 AUV在大地坐标系下的速度，采用模糊控制直线跟踪方法继续行驶，进行AUV自适应导引对接回收，其中AUV在大地坐标系下的速度通过惯性导航系统设备输出，参考航向角ψreference通过以下过程得到：
参考航向角ψreference的计算方式为：
ψreference＝ψcross+ψinit (1)
其中ψcross为直线跟踪角，ψinit为任务初始设定的参考航向角。
ψcross＝arctan(k*L) (2)
k为根据AUV与回收笼距离而确定的直线跟踪参数；本实施例中计算直线跟踪角ψcross时，在距离回收笼100米以外时，选择直线跟踪参数k＝0.03，使得航行器直线跟踪运动轨迹较为平滑，超调较少；在距离回收笼100米及100米以内时，选择直线跟踪参数k＝0.08，使得航行器在跟踪水声定位设备更新的回收笼位置时，航行器可以快速且有效地做出相应位姿调整，保证进入回收笼的航向和位置，即
L为航行器与回收笼之间的横向距离：
L＝-dE1*cos(ψ)-dN1*sin(ψ) (4)
而其中北向距离dN1和东向距离dE1的近似值为
dN1＝(latd-latcur)*R (6)
式中R为地球子午圈的半径。北向距离dN1和东向距离dE1两者均为近似值，但都满足导航精度。
本实施例中，当AUV与回收笼距离小于500米后未检测到超短基线USBL所需声源，则结束当前对接回收任务，并判断是否达到最大尝试次数，若未达到最大尝试次数则重新回到最后一个航路点执行新一次对接任务，若达到最大尝试次数，则上浮等待回收。
步骤3：在AUV接近回收笼的过程中，超短基线USBL测量出的信息有：航行器相对于回收笼之间的斜距R，航行器相对于回收笼的方位角θyz(右正左负)，航行器相对于回收笼的俯仰角θxz(上正下负)。由于在水下环境中，回收笼依靠GPS定位的经纬度精度不能满足回收需要，因此本发明采用声光导引组合，根据超短基线USBL测量结果以及光学手段，结合航行器自身的经纬度loncur和latcur，对回收笼经纬度lond和latd进行更新，更新方法如下：
建立声光导引的计算模型：
I(x)＝J(x)f(x)+(1-J(x))g(x) (7)
其中I(x)为声光融合下的回收笼经纬度计算值，J(x)为远距离声学系数，取值为1或者0，当未检测到回收笼口部位置的4个光信号时取值为1，反之为0。f(x)为基于USBL的回收笼经纬度测量结果。g(x)为基于光学信号的回收笼经纬度测量结果。加权系数的变化如图5。
利用上一周期更新后的回收笼经纬度lond和latd计算回收笼与AUV之间的实时距离Dis
在Dis>25m时，AUV沿着直线驶向回收笼的过程中，根据超短基线USBL测量结果对回收笼经纬度lond和latd进行更新：
USBL提供回收笼相对于AUV之间的东向距离dE2和北向距离dN2，得到更新后的回收笼经纬度lond和latd为：
lond＝loncur+dE2*ρ2 (9)
latd＝latcur+dN2*ρ1 (10)
纬度更新率ρ1和经度更新率ρ2能够根据大地测量学估算：
式中地球的第一偏心率e2＝0.0066943799013，地球半径R为6378137.0米。
在Dis≤25m时，控制水下航行器速度降低，并增加采用光学手段对回收笼的经纬度进行更新。
在回收笼口部位置具有四个呈矩形分布的蓝色光源，在回收笼底部中心具有一个绿色光源。而航行器姿态和深度控制也由主推、舵板逐渐切换至主推、舵板、辅推共同作用，最终将切换至主推、辅推共同控制。
当航行器的视觉传感器可以识别到装置口位置的蓝色目标光源时，根据检测到的不同数量的蓝色光源，航行器将作出不同的机动策略。图像处理计算机将相对位姿解算结果传输到航行器的运动控制计算机。检测到四个蓝色光源以后声学导引系统加权降为零，采用辅助推进器进行姿态控制，通过解算结果采用直线跟踪导引算法使航行器朝对接装置中心靠近。
光学系统产生水下图像后，首先进行水下图像滤波处理：
图像中存在着随机的白点或者黑点这是椒盐噪声的常见特征，绘制样本图像的统计直方图发现存在着类似于高斯正态函数的直方图分布噪声。针对水下图像存在的椒盐噪声特性采用中值滤波(设计一个包括像素周围邻近像素的模板；按照灰度值的大小重新排列模板中的像素；将灰度值最大和最小的像素点排列在模板的两边，取模板中间位置的像素值，最后用该像素值替代原本像素值。)处理。针对水下图像存在的椒盐噪声特性采用均值滤波(首先设计一个包含待处理像素的邻近像素模板；计算模板包含的所有像素均值；最后将计算的均值替换待处理的像素。)实现水下图像的复原。
目标光源的安装示意见图3。当航行器的视觉传感器可以识别到回收笼口部蓝色目标光源时，根据检测到的不同的蓝色光源数时，航行器将作出不同的机动策略：
当航行器相机看到小于特征点光源数目(n<4)时，无法解算位姿，航行器仍继续使用USBL获得的信息来进行控制，同时打开辅推加大左右横移量寻找剩下的蓝色光源，并且把速度再次降低，避免过快导致未找到全部光源便驶过终点。
当Dis<15m后，航行器相机仍未能看到4个特征点光源数目时，则结束当前对接回收任务，并判断是否达到最大尝试次数，若未达到最大尝试次数则重新回到最后一个航路点执行新一次对接任务，若达到最大尝试次数，则上浮等待回收。
当航行器相机看到对应的特征点光源数目时(n＝4)时，图像处理计算机将AUV 与四个特征点光源中心的相对位姿解算结果传输到航行器的运动控制计算机。检测到四个蓝色光源以后声学导引系统USBL加权降为零，采用辅助推进器进行姿态控制，通过解算结果采用直线跟踪导引算法使航行器朝对接装置中心靠近。其解算过程中需要自适应阈值分割，此处采用最大类间方差法(OTSU)，按照图像的灰度特征将图像分成背景和目标两部分。由于水下图像背景比较暗，假设图像I(x,y)的大小为M×N，前景(及目标)和背景的分割阈值记为T，图像中像素的灰度值小于阈值T的像素个数记作N1，像素灰度大于阈值T的像素个数记作N2,则有：
N1+N2＝M×N (15)
ω1+ω2＝1 (16)
μ＝μ1×ω1+μ2×ω2 (17)
g＝ω1×(μ-μ1)2+ω2×(μ-μ2)2 (18)
将公式16带入公式17，得到等价公式：
g＝ω1×ω2×(μ-μ2)2 (19)
其中属于前景的像素点数占整幅图像的比例记为ω1，其平均灰度为μ1，背景像素点数占整幅图像的比例为ω2，其平均灰度为μ2，图像的总平均灰度为μ，类间方差记为g，通过上述公式采用遍历的方法得到使类间方差最大的阈值T，记为自适应的阈值分割标准，对采集到的水下图像按照公式18进行分割求得二值化图像g(x,y)：
分割结果见图6。
仅通过灰度的相似性进行分割无法完整的识别目标光源的轮廓形状并且光斑区域会出现大小不一的空洞，这样便无法精确的求解出目标光源的中心点像素坐标，此时采用形态学梯度边缘检测方法。形态学的基本运算有膨胀与腐蚀、开运算与闭运算。采用闭运算处理，闭运算的数学表达：
其中A为水下图像，B为结构元，由于目标光源的真是形状为圆形，因此把结构元B定义为3×3的圆形结构元素形态学的闭运算是对图像先进性膨胀在进行腐蚀，闭运算通过圆形结构元素B与水下图像A进行不同形式的“与”操作将B遍历图像A的所有像素，水下图像中可以使边界向图像中心收缩，去除无意义的水下杂质，缩小前景面积更加逼近目标光源真实形状，并且还可以填充光源内部的小孔以及边缘的小凹陷部分。闭运算后的结果如图7所示，轮廓形状更加逼近光源形状，且识别光斑中的凹陷部分也被很好的填充。
把图像进行处理以后，采用EPnP算法进行位姿估计。
EPnP算法是一种O(n)时间复杂度具有封闭形式解的解算方案，该算法通过引入固定数量的虚拟控制点(4个)表示空间中的所有已知点坐标将复杂的多项式求解问题简化，避免了计算过程中未知量的交互利用线性解法计算虚拟控制点在目标坐标系中的位置对位姿进行最终求解，与最精确的迭代方法相比，精度损失很小，运算速度更快。
EPnP算法采用一组虚拟控制点的加权和来表示参考点，引入
个中间标量，原算法中通过构造四个虚拟控制点的方法，只需要引入4×3一共12个中间变量，进而求得闭式解。EPnP虚拟控制点模型见图8。
n个目标点在世界坐标系中的3D坐标已知，及Pi w，i＝1，2，...，n，可以在世界坐标系中构造四个控制点
上标w表示世界坐标系下的点坐标，上标c表示相机坐标系下的点坐标，将参考点的坐标表示为控制点坐标的加权和：
其中αij是齐次重心坐标，一旦虚拟控制点确定后，且满足4个控制点不共面的前提，αij,j＝1,2,3,4是唯一确定的。在相机坐标系中，存在同样的加权和关系，相机的外参定义为[R t]，那么虚拟控制点
和/>
之间存在如下关系：
当选取的第一个虚拟控制点为目标点的重心时，其余的控制点会都排列在数据的主方向上会使位姿估计算法的稳定性增加。设K是相机的内参矩阵，可以通过标定获得，利用用{ui}i＝1,...,n表示目标点{Pi}i＝1,...,n的2D投影点，则有：
其中，wi是射影深度因子，将每个
控制点用/>
来表示，由相机的针孔模型可知，内参矩阵K可以用焦距系数fu，fν图像坐标原点在像素坐标中的坐标 (uc,vc)。
从上式可以得到两个线性方程：
可以发现wi已经被消掉，联立所有n个目标点，得到线性的方程：
Mx＝0 (29)
其中
是一个由未知数组成的12向量，M是一个2n×12的矩阵
其中，vi是与M的N个非奇异值相关的M的右奇异向量的所有列的集合，可以通过12×12大小的矩阵MTM的零特征向量求得。通过研究发现矩阵M的核空间维数 N＝1,2,3,4，则上式描述的实际上是一个秩亏最小二乘问题，当N取不同的值时，由于相机的外参描述只是坐标变换，虚拟控制点之间的距离已知且没有发生改变，利用
可以得到不同的二次约束，下面通过讨论N＝1，2的情况表述求解过程。
例如当N＝1时，方程(上)简化为x＝βv，由
可得：
当N＝2，x＝β1v1+β2v2，相机坐标系中控制点的距离表述变为：
因为我们有四个控制点，可以转换为一个含有5个方程的线性系统：
Lβ＝ρ (33)
其中L是一个6×3的矩阵，由v1和v2组成，ρ是一个由二次距离
组成的一个六维的向量。通过对不同的二次约束求解出βi＝1,...,N的系数值之后，通过下式求解最小的重投影误差值。
最小的重投影误差值对应的βi＝1,...,N系数值则为EPnP位姿估算最优解。
从而得到AUV与回收笼蓝色光源中心位置的偏差量，包括横向偏差量、纵向偏差量和轴向偏差量。
利用辅推控制，将横向偏差量deltax、纵向偏差量deltay调整至0左右(误差 <＝50mm)同时将速度降低至能够保持航行状态的最低。随着航行器继续前进，当轴向偏差量deltaz小于1500mm时，且横向纵向误差满足需求后，此时不再调整航行器位姿，并且提高主推转速，维持当前航向前进(避免由于相机视野的限制，在距离较近时无法看到全部的蓝色特征光源无法解算位姿而产生错误)。
步骤4：当航行器调整到对接装置中心附近时，相机识别到笼内的绿色目标光源，图像处理计算机对绿色光源进行自适应跟踪与位姿计算结果进行修正，从而完成水下航行器的自主对接。
若未检测到绿色光源，则结束当前对接回收任务，并判断是否达到最大尝试次数，若未达到最大尝试次数则重新回到最后一个航路点执行新一次对接任务，若达到最大尝试次数，则上浮等待回收。
下面通过实验验证本发明设计的声光融合导引的自主水下航行器鲁棒对接回收方法的有效性。
首先根据发明内容中步骤1、2所述，首先航行器解算自身经纬度，读取对接装置的深度为3.2m，跟踪的朝向为笼口朝向。然后按照本发明的步骤3-4进行控制器解算，并获得实验结果。图9为航行器速度变化曲线，图10为航行器深度变化曲线，图11 为相对角度偏差图，图12为相对位置偏差。
由图9航行器的Vx速度曲线可得，在航行器下潜到期望深度前由于存在俯仰角的原因会产生一定的距离倒退，所以任务开始初始航行器的前向速度为负值，当达到期望深度后，主推开始旋转航行器开始加速并且保持0.4m/s的前向航向速度。图10 为截取航行器深度变化曲线中的一部分，在20s左右的时航行器的深度达到期望深度，与航行器开始加速的时间吻合。
图11与图12为航行器与对接装置相对的位姿偏差，即相对位姿估计的结果。由图11可以发现航行器与对接装置的俯仰角度偏差与任务设定的角度偏差相近，初始偏航角度偏差为8度左右，最后的解算角度0.04度满足对接角度要求。由图12可以看出，图中Z(图像深度距离)、X(水平偏差距离)的变化曲线反映了航行器下潜时后退的趋势。由相对位置偏差可以发现，航行器在1.84m左右丢失目标光源信息，此时横向偏差在0.3m，符合对接要求。水池对接实验数据验证了本论文的控制算法的有效性。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (9)
1.一种声光融合导引的自主水下航行器鲁棒对接回收方法，其特征在于：包括以下步骤：
步骤1：按照任务需求设定航路点，并根据回收笼位置及回收方向，将最后一个航路点设置在回收笼正前方；水下航行器以航路点导引方式依次完成航路点，并到达位于回收笼正前方的最后一个航路点；
步骤2：执行回收任务参考深度为回收笼深度，跟踪的航向为回收笼口的朝向，跟踪的终点为回收笼后方；根据参考航向角ψreference以及AUV在大地坐标系下的速度，采用模糊控制直线跟踪方法继续行驶，进行AUV自适应导引对接回收；当AUV与回收笼距离小于设定阈值后未检测到超短基线USBL所需声源，则结束当前对接回收任务，并判断是否达到最大尝试次数，若未达到最大尝试次数则重新回到最后一个航路点执行新一次对接任务，若达到最大尝试次数，则上浮等待回收；
步骤3：在AUV接近回收笼的过程中，采用声光导引组合，根据超短基线USBL测量结果以及光学手段，结合航行器自身的经纬度loncur和latcur，对回收笼经纬度lond和latd进行更新；
步骤4：当AUV调整到回收笼中心附近时，AUV相机识别到回收笼底部中心目标光源，AUV对回收笼底部中心目标光源进行自适应跟踪与位姿计算结果修正，完成水下航行器的自主对接；
若AUV未检测到回收笼底部中心目标光源，则结束当前对接回收任务，并判断是否达到最大尝试次数，若未达到最大尝试次数则重新回到最后一个航路点执行新一次对接任务，若达到最大尝试次数，则上浮等待回收。
2.根据权利要求1所述一种声光融合导引的自主水下航行器鲁棒对接回收方法，其特征在于：步骤2中，参考航向角ψreference通过以下过程得到：
参考航向角ψreference的计算方式为：
ψreference＝ψcross+ψinit
其中ψcross为直线跟踪角，ψinit为任务初始设定的参考航向角；
ψcross＝arctan(k*L)
k为根据AUV与回收笼距离而确定的直线跟踪参数；L为航行器与回收笼之间的横向距离。
3.根据权利要求2所述一种声光融合导引的自主水下航行器鲁棒对接回收方法，其特征在于：步骤2中，在水下航行器距离回收笼100米以外时，选择直线跟踪参数k＝0.03，在水下航行器距离回收站100米及100米以内时，选择直线跟踪参数k＝0.08。
4.根据权利要求2所述一种声光融合导引的自主水下航行器鲁棒对接回收方法，其特征在于：步骤2中，L为航行器与回收笼之间的横向距离：
L＝-dE1*cos(ψ)-dN1*sin(ψ)
其中北向距离dN1和东向距离dE1的近似值为
dN1＝(latd-latcur)*R
式中R为地球子午圈的半径。
5.根据权利要求1所述一种声光融合导引的自主水下航行器鲁棒对接回收方法，其特征在于：步骤3中，建立声光导引的计算模型：
I(x)＝J(x)f(x)+(1-J(x))g(x)
其中I(x)为声光导引组合下的回收笼经纬度计算值，J(x)为远距离声学系数，取值为1或者0，当检测到回收笼口部位置的光信号个数不足4个时取值为1，反之为0，f(x)为基于USBL的回收笼经纬度测量结果，g(x)为基于光学信号的回收笼经纬度测量结果；
利用上一周期更新后的回收笼经纬度lond和latd计算回收笼与AUV之间的实时距离
在Dis>LL时，AUV沿着直线驶向回收笼的过程中，根据超短基线USBL测量结果对回收笼经纬度lond和latd进行更新：LL为设定距离阈值；
在Dis≤LL时，控制水下航行器速度降低，并增加采用光学手段对回收笼的经纬度进行更新。
6.根据权利要求5所述一种声光融合导引的自主水下航行器鲁棒对接回收方法，其特征在于：根据超短基线USBL测量结果对回收笼经纬度lond和latd进行更新过程为：
USBL提供回收笼相对于AUV之间的东向距离dE2和北向距离dN2，得到更新后的回收笼经纬度lond和latd为：
lond＝loncur+dE2*ρ2
latd＝latcur+dN2*ρ1
其中纬度更新率ρ1和经度更新率ρ2根据大地测量学估算。
7.根据权利要求5所述一种声光融合导引的自主水下航行器鲁棒对接回收方法，其特征在于：在回收笼口部位置具有至少4个光源，在回收笼底部中心具有一个光源，回收笼口部位置光源颜色与回收笼底部中心位置光源颜色不同；
当Dis≤LL时，若AUV相机识别到回收笼口部位置的光源数小于4个，则航行器仍继续使用USBL获得的信息来进行回收笼经纬度更新，同时打开辅推加大左右横移量寻找剩下的蓝色光源，并且把速度再次降低；
如果当Dis<LT后，AUV相机识别到回收笼口部位置的光源数仍小于4个，则结束当前对接回收任务，并判断是否达到最大尝试次数，若未达到最大尝试次数则重新回到最后一个航路点执行新一次对接任务，若达到最大尝试次数，则上浮等待回收；
如果AUV相机识别到回收笼口部位置的光源数达到4个，则将USBL加权降为零，只利用图像解算结果采用直线跟踪导引算法使AUV朝回收笼中心靠近。
8.根据权利要求7所述一种声光融合导引的自主水下航行器鲁棒对接回收方法，其特征在于：LL取为25m，LT取15m。
9.根据权利要求7所述一种声光融合导引的自主水下航行器鲁棒对接回收方法，其特征在于：利用图像解算结果采用直线跟踪导引算法使AUV朝回收笼中心靠近过程中，对采集的图像采用最大类间方差法进行自适应阈值分割，采用形态学梯度边缘检测方法识别目标光源的轮廓形状，并采用EPnP算法进行位姿估计，得到AUV与回收笼口中心位置的偏差量，包括横向偏差量、纵向偏差量和轴向偏差量；利用辅推控制，将横向偏差量、纵向偏差量调整至接近0，同时将速度降低至能够保持航行状态的最低速度；当轴向偏差量deltaz小于设定值时，且横向纵向误差满足需求后，不再调整航行器位姿，并且提高主推转速，维持当前航向前进。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111065374.8A CN114003029B (zh) | 2021-09-12 | 2021-09-12 | 一种声光融合导引的自主水下航行器鲁棒对接回收方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111065374.8A CN114003029B (zh) | 2021-09-12 | 2021-09-12 | 一种声光融合导引的自主水下航行器鲁棒对接回收方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114003029A CN114003029A (zh) | 2022-02-01 |
| CN114003029B true CN114003029B (zh) | 2023-06-30 |
Family
ID=79921316
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111065374.8A Active CN114003029B (zh) | 2021-09-12 | 2021-09-12 | 一种声光融合导引的自主水下航行器鲁棒对接回收方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114003029B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106444838A (zh) * | 2016-10-25 | 2017-02-22 | 西安兰海动力科技有限公司 | 一种自主水下航行器的精确航路跟踪控制方法 |
| CN109238291A (zh) * | 2018-10-26 | 2019-01-18 | 河海大学 | 一种水面无人船导引缆回收自主水下航行器的规划方法 |
| JP2021034050A (ja) * | 2019-08-21 | 2021-03-01 | 哈爾浜工程大学 | 強化学習に基づくａｕｖ行動計画及び動作制御方法 |
| CN112987798A (zh) * | 2019-12-17 | 2021-06-18 | 中国科学院沈阳自动化研究所 | 基于声光联合导引重型auv动/静态目标自主对接方法 |
-
2021
- 2021-09-12 CN CN202111065374.8A patent/CN114003029B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106444838A (zh) * | 2016-10-25 | 2017-02-22 | 西安兰海动力科技有限公司 | 一种自主水下航行器的精确航路跟踪控制方法 |
| CN109238291A (zh) * | 2018-10-26 | 2019-01-18 | 河海大学 | 一种水面无人船导引缆回收自主水下航行器的规划方法 |
| JP2021034050A (ja) * | 2019-08-21 | 2021-03-01 | 哈爾浜工程大学 | 強化学習に基づくａｕｖ行動計画及び動作制御方法 |
| CN112987798A (zh) * | 2019-12-17 | 2021-06-18 | 中国科学院沈阳自动化研究所 | 基于声光联合导引重型auv动/静态目标自主对接方法 |
Non-Patent Citations (1)
| Title |
|---|
| 基于模型预测控制的自主水下航行器模糊回坞导引算法;齐贝贝;齐仲里;王林林;;国外电子测量技术(第04期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114003029A (zh) | 2022-02-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sinisterra et al. | Stereovision-based target tracking system for USV operations | |
| CN109239709B (zh) | 一种无人船的局部环境地图自主构建方法 | |
| CN110610134B (zh) | 一种无人艇自主对接方法 | |
| CN113657256A (zh) | 一种无人艇载无人机海空协同视觉跟踪与自主回收方法 | |
| CN109213204A (zh) | 基于数据驱动的auv海底目标搜寻航行系统及方法 | |
| CN110658826A (zh) | 一种基于视觉伺服的欠驱动无人水面艇自主靠泊方法 | |
| CN113778099A (zh) | 基于NDT算法和Hybrid A*算法的无人船路径规划方法 | |
| CN113525631A (zh) | 一种基于光视觉引导的水下终端对接系统及方法 | |
| CN115031718B (zh) | 一种多传感器融合的无人船同步定位与建图方法(slam)及系统 | |
| CN115903800A (zh) | 一种基于领导者协同的多无人船编队的严格安全控制方法 | |
| Lin et al. | Auv docking method in a confined reservoir with good visibility | |
| CN114061565B (zh) | 一种无人船舶slam及其应用方法 | |
| Feng et al. | Automatic tracking method for submarine cables and pipelines of AUV based on side scan sonar | |
| Wu et al. | Autonomous UAV landing system based on visual navigation | |
| Chen et al. | A novel unmanned surface vehicle with 2d-3d fused perception and obstacle avoidance module | |
| CN114119752A (zh) | 基于gnss和视觉的室内外衔接的机器人定位方法 | |
| CN114003029B (zh) | 一种声光融合导引的自主水下航行器鲁棒对接回收方法 | |
| Yin et al. | Study on underwater simultaneous localization and mapping based on different sensors | |
| CN109459046B (zh) | 悬浮式水下自主航行器的定位和导航方法 | |
| Jiang et al. | Underwater loop-closure detection for mechanical scanning imaging sonar by filtering the similarity matrix with probability hypothesis density filter | |
| Calado et al. | Obstacle avoidance using echo sounder sonar | |
| Sang et al. | An autonomous underwater vehicle simulation with fuzzy sensor fusion for pipeline inspection | |
| CN115373383A (zh) | 一种垃圾回收无人艇的自主避障方法、装置及相关设备 | |
| CN112577463B (zh) | 姿态参数修正的航天器单目视觉测距方法 | |
| Qingqing et al. | Detecting water reflection symmetries in point clouds for camera position calibration in unmanned surface vehicles |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |