基于Unity3D的月面复杂地形场景构建及模拟驾驶系统.pdf

第6卷 第1期2 0 2 3年1月飞控与探测F l i g h tC o n t r o l&D e t e c t i o nV o l.6N o.1J a n u a r y2 0 2 3基于U n i t y 3 D的月面复杂地形场景构建及模拟驾驶系统*陈李炜1,曹 涛2,3,易建军1,张雪松1,赵逸凡1,王庆明1(1.华东理工大学 机械工程系上海2 0 0 2 3 7;2.上海航天控制技术研究所上海2 0 1 1 0 9;3.上海市空间智能控制技术重点实验室上海2 0 1 1 0 9)摘 要:为了验证未来月面复杂地形操控技术方案,基于U n i t y 3 D物理引擎对月表模型、月球物理环境以及月面光照环境进行了模拟,设计并开发了一套用于模拟月球自主进行着陆点选址以及在线航天器着陆轨迹规划的系统。该模拟系统的工作过程包含了:基于面阵雷达与立体相机信息融合的局部月表重建、基于3 D重建结果的月表地形着陆代价的快速评估与自动着陆选址,以及使登陆器能够到达选址目标的最优燃料消耗着陆运动规划,实现了月球着陆器短距离自主选址着陆的模拟系统构建,同时也从仿真的角度初步验证了“地形重建-地形评估-自主选址-软着陆”这一自主着陆过程的可行性。关键词:U n i t y 3 D开发;月球登陆器;地形评估;软着陆 中图分类号:V 4 4 8.2 5文献标志码:A文章编号:2 0 9 6-5 9 7 4(2 0 2 3)0 1-0 0 4 8-0 8A u t o n o m o u sC o m p l e xT e r r a i nS o f tL a n d i n gS i m u l a t i o nS y s t e mo fL u n a rL a n d e rB a s e do nU n i t y 3 DCHE NL i w e i1,C AOT a o2,3,Y I J i a n j u n1,Z HANGX u e s o n g1,Z HAOY i f a n1,WANGQ i n g m i n g1(1.D e p a r t m e n to fM e c h a n i c a lE n g i n e e r i n g,E a s tC h i n aU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y,S h a n g h a i 2 0 0 2 3 7;2.S h a n g h a iA e r o s p a c eC o n t r o lT e c h n o l o g yI n s t i t u t e,S h a n g h a i 2 0 1 1 0 9;3.S h a n g h a iK e yL a b o r a t o r yo fA e r o s p a c e I n t e l l i g e n tC o n t r o lT e c h n o l o g y,S h a n g h a i 2 0 1 1 0 9)A b s t r a c t:I no r d e r t ov e r i f yt h e t e c h n i c a l s c h e m eo f a u t o n o m o u s l a n d i n go nc o m p l e x l u n a r t e r r a i n f o r l u n a r l a n d e r s,t h i sp a p e rs i m u l a t e s t h e l u n a rs u r f a c em o d e l,l u n a rp h y s i c a l e n v i r o n m e n ta n dl u n a r i l l u m i n a t i o ne n v i r o n m e n tb a s e do nU n i t y 3 Dp h y s i c a l e n g i n e.A t t h es a m e t i m e,i td e s i g n sa n dd e v e l o p sas i m u l a t i o ns y s t e mf o rs i m u l a t i n gt h ea u-t o n o m o u s l a n d i n gs i t es e l e c t i o na n dr e a l-t i m e l a n d i n gt r a j e c t o r yp l a n n i n go f l u n a r l a n d e r s.T h ew o r k i n gp r o c e s so ft h es i m u l a t i o ns y s t e mi n c l u d e s:l o c a l l u n a r s u r f a c e r e c o n s t r u c t i o nb a s e do n t h e i n f o r m a t i o n f u s i o no f a r e a a r r a y r a d a ra n ds t e r e oc a m e r a,r e a l-t i m e a s s e s s m e n t o f l u n a r s u r f a c e t e r r a i n l a n d i n gc o s t a n da u t o m a t i c l a n d i n g l o c a t i o nb a s e do n3 Dr e c o n s t r u c t i o nr e s u l t s,a n do p t i m a l f u e l c o n s u m p t i o nl a n d i n gm o t i o np l a n n i n gt oe n a b l et h e l a n d e rt or e a c ht h el o c a t i o nt a r g e t.T h es i m u l a t i o n s y s t e m o fs h o r t-r a n g e a u t o n o m o u sl o c a t i o nl a n d i n g o ft h el u n a rl a n d e ri sc o n s t r u c t e d,a n dt h e f e a s i b i l i t yo f t h ea u t o n o m o u s l a n d i n gp r o c e s so f“t e r r a i nr e c o n s t r u c t i o n-t e r r a i nl a n d i n gr i s ke-v a l u a t i o n-a u t o n o m o u s l a n d i n gp o i n t s e l e c t i o n-s o f t l a n d i n g”i sp r e l i m i n a r i l yv e r i f i e df r o mt h ep e r s p e c t i v eo fs i m u l a-t i o n.K e y w o r d s:U n i t y 3 Dd e v e l o p m e n t;l u n a r l a n d e r;t e r r a i ne v a l u a t i o n;s o f t l a n d i n g*收稿日期:2 0 2 2-0 3-1 1;修回日期:2 0 2 2-0 9-2 8基金项目:上海航天创新基金(S A S T-2 0 1 9-0 8 0)作者简介:陈李炜(1 9 9 7),男,博士生。E-m a i l:1 8 9 1 6 3 0 9 2 9 11 6 3.c o m通信作者简介:易建军(1 9 6 9),男,博士,教授。E-m a i l:j j y i e c u s t.e d u.c n第1期陈李炜,等:基于U n i t y 3 D的月面复杂地形场景构建及模拟驾驶系统0 引 言航天飞行器在行星表面的自主软着陆是探索未知行星不可缺少的关键技术。载人或无人飞行器都不可避免地需要经过动力下降进行软着陆,为了着陆的平稳与安全,一般情况下着陆点都选择地势开阔且平坦的地形。1 9 6 9年的载人登月任务A p o l l o1 1,飞行器在动力下降阶段借助发动机反推实现了减速,并通过人工控制进行避障飞行,最后垂直下降并软着陆于月表。然而,随着未来行星探索任务的扩展以及复杂化,飞行器面临不得不在复杂地形区域着陆的情况,对于行星地表的感知、评估选址,以及末端的制导着陆则非常关键。本文通过U n i t y 3 D物理引擎对这类复杂地形软着陆的情况开发了一套模拟系统,用作自主选址着陆方法的测试与验证的平台。本文所构建的复杂地形自主软着陆模拟平台是一个整合了月表地形感知、月表地形评估、自主选址与着陆规划的系统,从其技术构成而言包含了场景三维重建技术、地形评估技术以及飞行器动力下降制导技术。月表地形感知的第一步是获取尽可能稠密的月表三维信息,本文所实现的月面地形表面重建是场景三维重建技术的应用扩展。该技术通常使用立体视觉传感器、激光雷达传感器或多种传感器的混合,并以多视图几何原理或光学成像获取重建对象的深度数据,随后进行数据预处理转化为点云数据,并通过配准算法实现点云信息的配准融合,最终将离散的点云转化为连续表面,把真实场景转换成符合计算机逻辑表达的数学模型。在以往的空间任务中常使用视觉传感器作为主要感知手段,例如小行星表面数字地形建模,行星表面重建工作中存在大量的基于立体视觉和立体光度测量方法。其中,前者基于三角测量方法,得到视野中的像素偏差,随后基于极线约束像点进行标定1-3。该方法模拟了人类视觉的感知原理,但是算法复杂,对参数敏感导致稳定性和误差都会随着运行时间而恶化。后者基于立体成像和摄影测量的原理,从各种机载的传感器得到行星表面的深度数据,由于包含了反照率信息,在处理低入射角的图像时优势明显。但在实际使用中,该方法存在着形状特征被随机消除的缺点,会造成像素误差。另一方面,为了提升重建精度,需要更多的有关场景参数的额外数据量。基于单目序列图像的三维重建,通过几何投影的原理构建数字化三维模型,对场景参数的依赖度低,可以根据实际情况选择拟合方案以降低累积误差,同时随着算法复杂度降低,鲁棒性和实时性也相应增强了4-6。近1 0年,随着激光雷达传感器的兴起,基于激光同步定位与地图构建(S i m u l t a n e o u sL o c a l i z a t i o na n dM a p p i n g,S L AM)方法开展的场景重建表现出了超越以往基于立体视觉重建方法的优秀性能。例如2 0 1 4年,L O AM算法提出了一种利用2轴激光雷达在6自由度范围内进行实时里程测量和地图绘制的方法,通过运动估计进行点云的运动补偿,实现了低漂移和低计算复杂度,无需高精度测距或惯性测量。在后续的工作中,大量学者基于该算法进行了改进,开发出了T-L O AM、L e G O-L O AM等性能更好的方法7-8;同时对于面阵激光雷达的场景重建也设计了L i v o x-L O AM算法,针对扫描方式与传统机械雷达不同的问题,优化了点云特征的提取,适配了固态雷达,并表现出了不错的性能9。由于激光雷达的能量密度高,抗干扰能力更强,是近距离行星表面重建的新选择;另一方面,视觉传感器虽然对参数敏感,易受干扰,但是局部重建结果更为稠密。因此,本文在模拟系统中的月表地形重建采用了视觉-激光融合的月表重建方法,使用了面阵激光雷达的L i v o x-L OAM方法对月面进行重建模拟,同时也使用了视觉传感器作为协同重建的手段。快速的月面地形评估是实现自主软着陆选址的前提条件,在完成月表三维重建后,即可开展月面复杂地形的评估,计算月表各区域的着陆安全性。近年来,随着无人车自主行驶、无人机自主降落任务对环境感知能力的需求增加,大量学者对复杂地形的可通过性和着陆安全性的评估技术开展了研究,这项关键技术对月表地形评估具有很高的参考价值。B.A y h a n和C.Kw a n1 0使用G a b o r特征和支持向量机方法以及S e g-N e t方法,对地形进行有效评估,为起伏地形上的无人机提供着陆安全性评估。2 0 2 0年,Y a nL.等1 1提出了一种从激光雷达点云中自主选择安全着陆点的新方法,采用了主成分分析(P r i n c i p a lC o m p o n e n tA n a l y s i s,P C A)和改进的面积增长算法检测平坦区域后,基于渐进样本一致性(P r o g r e s s i v eS a m p l eC o n s e n s u s,P R O S A C)算法实现平面拟合,以评估地形复杂度。94解决复杂地形的自动驾驶问题也需要先进行地形评估,有学者通过大量的地形样本数据集训练,得到了能够有效评估地形的神经网络模型。Y.K o b a y a s h i等1 2训练了一个高斯过程回归器(G a u s s i a nP r o c e s sR e g r e s s i o n,G P R),以预测车辆在地形上移动时的振动(作为地形可通过性的度量),并结合R G B摄像机图像检测到的地形纹理特征对地形进行评估。M.Q u a n n等1 3提出了一种基于概率能量成本预测的越野导航策略,同样通过高斯过程回归器,实现了从当前机器人姿态、地形坡度(沿机器人运动方向)和灰度图像到运行功耗的映射。能量消耗的地图最终被用来预测路径穿越期间的累计能量消耗,这一类的方法往往需要获取机器人与地形的接触状态,以评判地形的可通过性。完成地形评估并确定着陆选址后,需要进行相对精确的软着陆,为此飞行器需要进行在线实时的软着陆运动规划。大量研究工作已经开发出针对动力下降制导问题的近似解决方案,其中以凸优化为主的在线轨迹优化算法在实际工程中有十分优秀的表现。在文献 1 4的工作中,研究团队将飞行器在火星动力软着陆场景中的轨迹规划问题描述为一个非凸优化问题,并通过引入松弛变量,将松弛凸约束代替非凸约束进行问题凸化,证明了松弛凸约束对于该问题是无损的,最终得 到 一 个 二 阶 锥 规 划 问 题。在NA S A的A D A P T项目中验证了该算法的有效性。随后,该团队又对先前的燃料最优轨迹规划加以改进,优化了着陆误差1 5。随着人工智能算法的兴起,不少学者借助学习算法实现了实时软着陆规划,文献 1 6将迭代线性二次调节器作为控制器产生初始轨迹,并通过多层神经网络拟合制导策略,利用控制器监督策略学习得到可行的策略。考虑到机载计算机的处理能力、地形评估算法的复杂度以及着陆运动规划运行的实时性要求,本文在模拟器中使用了面阵激光雷达与立体相机进行三维重建,并采用了一种复合的地形评估函数对地形着陆代价进行评估,在得到选址结果后,采用了基于凸优化的燃料最优的着陆规划,实现了飞行器在月面高度1 k m以下的软着陆。1 模拟系统构建1.1 软件架构模拟系统的软件架构主要分为两部分:基于U n i t y 3 D开发的着陆器飞行场景,以及基于机器人操作系统(R o b o tO p e r a t i n gS y s t e m,R O S)开发的传感器数据获取与处理系统。基于U n i t y 3 D物理引擎所开发的等比的仿真场景搭建与相关功能(如图1所示),其具体内容包含了:图1 模拟器系统架构F i g.1 T h e s i m u l a t i o ns y s t e ma r c h i t e c t u r e05飞控与探测第6卷第1期陈李炜,等:基于U n i t y 3 D的月面复杂地形场景构建及模拟驾驶系统 1)仿真月表的构建:在U n i t y 3 D中使用地形制作插件绘制长、宽均为6 k m的月面地形,并模拟了月面常见的地形特征,例如:起伏的地形、多种尺寸的环形山、月坑。另一方面,为地形设置了高分辨率的碰撞模型。2)月球环境模拟:在U n i t y 3 D的场景物理属性中设置了月球重力,并将地面摩擦和接触刚性设置为可调参数。3)飞行器特性模拟:在着陆器3 D模型的基础上,开发了其基本功能的仿真与实时动态参数的仿真,具体表现为设置了飞行器所具有的主发动机,调姿发动机的开关功能、可调节的矢量推力仿真、基于比冲参数的发动机燃料消耗计算、总体质量参数的动态变化仿真、飞行器质心变化仿真和飞行器刚体碰撞仿真。基于R O S开发的孪生缩比环境负责对视觉以及激光雷达传感器的数据仿真与处理,其具体内容包含了:1)孪生缩比场景的构建:将U n i t y 3 D中的整体场景进行1 0倍缩比,从而降低运行时的计算机资源占用率,提升孪生场景运行的实时性。2)传感器仿真:使用R O S中的虚拟立体相机仿真以及虚拟激光雷达对飞行器下方的月表进行扫描,通过订阅R O S中的传感器数据话题实时回传仿真的深度图像与点云数据,作为重建算法的输入。3)相关算法的开发:进行月面地形重建、地形评估以及着陆轨迹规划的算法。基于U n i t y 3 D的仿真场景与基于R O S的孪生缩比场景之间通过UD P通信进行信息交互,U n i-t y 3 D将飞行器所在的位置、速度发送给R O S,而R O S将 自 主 选 址 和 着 陆 规 划 的 结 果 发 送 给U n i t y 3 D仿真场景,驱使飞行器通过发动机的控制实现最终的软着陆。1.2 月面自主选址软着陆流程模拟中的着陆流程(如图2所示)分为4个阶段:1)动力下降阶段:飞行器按照既定程序从绕月轨道减速下降,并调整飞行姿态与速度,逐渐接近预定的着陆区域上空。2)月表重建与地形评估选址阶段:在飞行器下降到指定高度后,启用面阵雷达和立体相机的仿真功能,对飞行器下方的月面进行成像,并基于重建算法开展实时地形重建,在完成着陆范围的地形重建后,采用地形评估函数对重建地形的稠密点云进行着陆代价计算,并选择最佳着陆地点,返回着陆点坐标。3)软着陆运动规划阶段:在获取上一个阶段的着陆目标点后,根据飞行器与目标点的相对位置和相对速度,以及当前飞行器性能参数进行着陆规划,获取当前状态下燃料消耗最优的软着陆运动规划结果。4)软着陆控制阶段:控制飞行器遵循上一阶段得到的运动规划结果完成着陆。图2 月面自主选址软着陆流程F i g.2 S o f t l a n d i n gp r o c e s s f o ra u t o n o m o u s s i t e s e l e c t i o no nt h eL u n a r s u r f a c e 第二以及第三阶段是整个着陆流程中最关键的2个过程,也是本文所构建的仿真系统中各关键算法的集中表现,后续2章将分别对月表地形的重建、评估和选址以及针对选址目标的软着陆运动规划进行说明。152 月表的重建与地形评估选址本文的模拟系统进行地形评估所使用的传感器为面阵激光雷达与立体相机,通过对地形的扫描与重建,获取着陆区域局部地形相对稠密的点云。2.1 月面地形的稠密重建本文使用R O S中的面阵激光雷达与立体相机作为主要传感器,对飞行器下方的月表地形进行稠密重建,基础算法使用了由L O AM算法改进而来的L i v o x-L O AM算法进行基于面阵激光雷达的地形重建。同时利用立体相机对点云缺失部分进行补偿,从而获取更为稠密的地形重建结果,算法主要流程为:数据预处理、点云配准、后端优化以及点云集成4个部分。由于立体相机的加入,在原始算法的基础上加入了稠密点云的生成和基于时间戳对齐的配准环节。月面3 D重建流程如图3所示。图3 月面3 D重建流程F i g.3 T h ep r o c e s so fL u n a r s u r f a c e3 Dr e c o n s t r u c t i o n在R O S中,对孪生缩比场景进行重建,得到局部月面地形如图4所示。图4 月面3 D重建效果F i g.4 L u n a r s u r f a c e3 Dr e c o n s t r u c t i o n2.2 月面复杂地形的着陆风险评估将重建获得的点云数据在高度方向进行归一化并以灰度形式进行描述,从而获得高程图(如图5所示)。对地形的评估采用3个参数进行描述:地形梯度、地形起伏程度、局部地形的最大高度差,并采用加权计算的方式得到综合着陆风险模型。图5 月表地形的高程图表示(局部)F i g.5 E l e v a t i o nm a pr e p r e s e n t a t i o no f l u n a rs u r f a c e t o p o g r a p h y(L o c a l)后续的月面复杂地形着陆风险评估均建立在高程地图的栅格化处理基础上。通过在栅格化的高程地图上计算栅格法向量实现对地形梯度的计算:对于每个栅格法向量,取4个样本(当前高程地图像素的左、右、上、下邻接像素),由此获取x与y方向上的梯度。将一个3 D向量投影到x=0,y=0的平面上,与2 D法线的斜率相同,并将其归一化。通过遍历栅格化高程图,计算每个栅格法线的斜率。该方法具有较高的计算速度,在模拟过程中能够实时地对重建地形进行评估。地形起伏程度Lr o u g h n e s s的评估则通过式(1)和式(2)计算,采用区域内高程图高度值的标准差作为衡量方法。Lr o u g h n e s s=1NNi=1(hi-h)2(1)h=1NNi=1hi(2)其中,N表示网格邻域内的网格个数;hi表示第i个栅格内的高程图高度值;h表示区域内的高程图高度均值。地形局部高度差表明了地形区域内的最大高度和最小高度之间的差,这直接关联到飞行器在该区域的地表着陆后的最大可能倾斜程度,因此也将其25飞控与探测第6卷第1期陈李炜,等:基于U n i t y 3 D的月面复杂地形场景构建及模拟驾驶系统作为综合风险模型中的一个参数。综合着陆风险模型对以上多种参数的加权融合,如下所示。综合着陆风险模型r1=1,rs l o p eTs l o p e或rr g sTr g s或rm h dTm h dws l o p ers l o p e+wr g srr g s+wm h drm h d,其他(3)ws l o p e+wr g s+wm h d=1(4)rs l o p e=Ds l o p eTs l o p e(5)rr g s=Lr o u g h n e s sTr g s(6)rm h d=m a x HTm h d(7)其中,r1表示着陆点的风险值,取值范围为0,1;rs l o p e为坡度风险值;rr g s为地面起伏风险值;rm h d为局部地形的最大高度差风险值;ws l o p e表示地形梯度风险的权重;wr g s表示地面起伏程度风险的权重;wm h d表示局部地形的最大高度差风险的权重;Ts l o p e表示地形梯度风险阈值;Tr g s表示地面起伏风险阈值;Tm h d表示地形局部高度差风险阈值;Ds l o p e为地形梯度的计算,由规定范围的栅格区域内的平均梯度计算所得;m a x H为局部地形的最大高度差,由规定范围的栅格区域内的最大高度差计算所得。对于评估得到的风险地图进行阈值分割,得到可以着陆的区域,最后对可以着陆的区域进行连通域的面积评估,确保其有足够的着陆面积提供给飞行器进行软着陆。使用上述的地形评估模型对重建后的地表进行整体的计算,得到了飞行器下方区域可视化的评估结果。本文对两种地形的各个区域进行了评估,获得了如下的效果(图6(a),其中绿色区域为满足着陆安全性的区域,红色区域为存在超过阈值的着陆风险区域。在最终的选址过程中,考虑到飞行器对着陆区域的大小有一定的要求,通过遍历栅格化的着陆代价地图,设置对可安全着陆区域中以各个栅格为中心所绘制的圆形区域中的最大面积区域作为最终着陆选址,如图6(b)所示。3 软着陆运动规划本文所设计的模拟器在进行末端制导软着陆时,采用了燃料消耗最少的最优着陆轨迹规划方法。模拟器选定着陆点后,求解一个使着陆器在初始条件 P,V下最终到达目标着陆点的飞行器(a)着陆风险可视化(b)着陆区域选址可视化图6 局部月表的着陆风险评估与选址方法可视化F i g.6 V i s u a l i z a t i o no f l a n d i n gr i s ka s s e s s m e n t l a n d i n gs i t e s e l e c t i o nm e t h o do n l o c a lL u n a r s u r f a c e推力变化曲线,使得燃料消耗最少,其中P为飞行器的初始位置向量,V为飞行器的初始速度向量。该问题的描述为一个优化问题m i ntf,Tc()tf0 Tc(t)dt(8)优化约束为r(t)=g+Tc(t)m(t)(9)m(t)=-Tc(t)(1 0)01Tc(t)2(1 1)r(t)X,t0,tf(1 2)m(0)=mw e t,m(tf)=md r y(1 3)r(0)=r0,r(0)=r0(1 4)eT1r(tf)=0,e2e3r(tf)=q,r(tf)=0(1 5)其中,Tc(t)表示随时间变化的飞行器推力向量;35tf为着陆过程所消耗的总时间;g表示月球的重力加速度向量;为比冲;1、2为推力大小的上下界;X为轨迹规划的倒锥形约束边界;r(t)为t时刻飞行器的位置向量;r0表示飞行器的初始位置向量;m(t)表示飞行器随时间变化的质量;mw e t表示开始制导时的燃料与飞行器本身质量总和;md r y为完成着陆后的燃料与飞行器本身质量总和;e1、e2、e3为垂直于下滑道约束向量分量的3个单位向量,描述了轨迹的始末状态分别处于锥形约束空间的位置,其中q是着陆目标点在着陆平面上的二维向量表示。在文献 1 4-1 6中,通过引入松弛变量的方式,将以上非凸优化问题转化为凸优化问题,进而使用c v x p y工具求解凸优化问题。如图7所示。图7 文献 1 4-1 6中所描述的制导轨迹在空间中的圆锥约束F i g.7 T h ec o n i cc o n s t r a i n t o fg u i d a n c e t r a j e c t o r y i ns p a c e i nr e f e r e n c e s1 4-1 6在模拟器中,预留了燃料量、最大推力、矢量推力限制角度、比冲、飞行器质量作为用户可调参数,用于模拟不同种类、不同状态的飞行器。实验中,飞行器的着陆轨迹和推力变化如图8所示。(a)着陆轨迹(b)推力变化曲线图8 软着陆规划结果F i g.8 P l a n n i n gr e s u l t o f s o f t l a n d i n g可见,在优化燃料的基础上,规划器也能够给出相对平稳光滑的着陆曲线。4 结 论本文所描述的基于U n i t y 3 D的月球探测器复杂地形软着陆模拟系统,完成了月面地形重建、着陆风险代价评估,以及指定点着陆的轨迹规划任务。模拟系统借助R O S构建了外部孪生缩比场景,模拟了实时面阵激光雷达以及立体相机所获取的数据,基于L i v o x-L OAM算法以及深度图像的点云补全实现了月面地形的稠密重建。在此基础上,本文针对重建所得的月面稠密地形,设计了一种高效的地形评估方法,实现了月表地形的着陆风险评估,开展低风险着陆区域选址,最终基于凸优化方法求解燃料最优的软着陆运动规划问题,使飞行器以较高的精度和最少的燃料消耗降落到指定的低风险着陆区域范围内。本文所设计的模拟系统为激光雷达和立体相机在月面自主软着陆中的应用提供了较为完整的功能模拟和数据处理算法的实现,最终将自主选址结果和着陆运动规划相结合,形成一个凸优化问题进行求解,实现了自主地形探测、自主选址、自主开展着陆运动规划的月球探测器复杂地形软着陆功能。该模拟系统具备一定的动力学仿真能力,并开放了与飞行器自身物理属性、飞行性能相关的可调节参数,使得模拟器具有针对不同的飞行器、不同地形的着陆流程进行模拟的能力。但是由于本模拟系统采用了U n i t y 3 D与R O S的联合仿真,因此这2个孪生环境之间信息传输必然存在着一定的延迟和数据丢失,对仿真的实时性和精度有45飞控与探测第6卷第1期陈李炜,等:基于U n i t y 3 D的月面复杂地形场景构建及模拟驾驶系统一定的影响。因此,在未来的开发中需要提升系统的集成度,在U n i t y 3 D仿真环境中实现多种传感器成像仿真以及传感器数据处理。此外,本系统在仿真细节和功能拓展上还具有优化空间,在未来的研究中需要对系统的动力学和太空环境视觉传感器的成像效果进行仿真,以及对人在回路的介入功能方面开展进一步的优化。参考文献(R e f e r e n c e s)1 I N T I L L ESS,B O B I C K A F.D i s p a r i t y-s p a c ei m a g e sa n dl a r g eo c c l u s i o ns t e r e oC/P r o c e e d i n g so fE u r o p e a nC o n-f e r e n c eo nC o m p u t e rV i s i o n.B e r l i n:S p r i n g e r,1 9 9 4:1 7 9-1 8 6.2S UNJ,S HUM H Y,Z HE N GNN.S t e r e om a t c h i n gu s i n gb e l i e fp r o p a g a t i o nC/P r o c e e d i n g so fE u r o p e a nC o n f e r-e n c eo nC o m p u t e rV i s i o n.B e r l i n:S p r i n g e r,2 0 0 2:5 1 0-5 2 4.3MA R U Y AM,N EMO T OK,T A K A S H I MA Y.T e x t u r eb a s e d3 Ds h a p er e c o n s t r u c t i o nf r o m m u l t i p l es t e r e oi m a g e sC/P r o c e e d i n g so f 1 9 9 21 1t hI A P RI n t e r n a t i o n a lC o n f e r e n c eo nP a t-t e r nR e c o g n i t i o n.T h eH a g u e,N e t h e r l a n d s:I E E E,1 9 9 2:1 3 7-1 4 0.4P A L M E REE,H E A DJN,G A S K E L LRW,e t a l.M e r c a t o r-I n d e p e n d e n tr o v e rl o c a l i z a t i o nu s i n gs t e r e o p h o t o c l i n o m e t r ya n dp a n o r a m i c i m a g e sJ.E a r t ha n dS p a c eS c i e n c e,2 0 1 6,3(1 2):4 8 8-5 0 9.5L O R E N ZDA,O L D SR,MAY A,e ta l.L e s s o n sl e a r n e df r o m O S I R I S-R E xa u t o n o m o u sn a v i g a t i o nu s i n gn a t u r a l f e a-t u r et r a c k i n gC/P r o c e e d i n g so f2 0 1 7I E E E A e r o s p a c eC o n f e r e n c e.B i gS k y:I E E E,2 0 1 7:1-1 2.6S A X E NA A,M I N S,N G A Y.3-Dr e c o n s t r u c t i o nf r o ms p a r s ev i e w su s i n gm o n o c u l a rv i s i o nC/P r o c e e d i n g so fI E E EI n t e r n a t i o n a lC o n f e r e n c eo nC o m p u t e rV i s i o n.R i od eJ a n e i r o,B r a z i l:I E E E,2 0 0 7:1-8.7S HAN T,E N G L O T B.L e GO-L OAM:l i g h t w e i g h ta n dg r o u n d-o p t i m i z e dl i d a ro d o m e t r ya n d m a p p i n go nv a r i a b l et e r r a i nC/P r o c e e d i n g so f2 0 1 8I E E E/R S JI n t e r n a t i o n a lC o n f e r e n c eo nI n t e l l i g e n tR o b o t sa n dS y s t e m s(I R O S).M a c a o:I E E E,2 0 1 9:4 7 5 8-4 7 6 5.8Z HOUP,GUOX,P E IX,e t a l.T-L OAM:t r u n c a t e d l e a s ts q u a r e sL i D A R-o n l yo d o m e t r ya n dm a p p i n g i nr e a l t i m eJ.I E E ET r a n s a c t i o n so nG e o s c i e n c e a n dR e m o t eS e n s i n g,2 0 2 1(9 9):1-1 3.9L I NJ,Z HAN GF.L o a ml i v o x:a f a s t,r o b u s t,h i g h-p r e c i-s i o nL i D A Ro d o m e t r ya n dm a p p i n gp a c k a g ef o rL i D A R so fs m a l lF o VC/P r o c e e d i n g so f2 0 2 0I E E EI n t e r n a t i o n a lC o n f e r e n c eo nR o b o t i c sa n d A u t o m a t i o n(I C R A).P a r i s:I E E E,2 0 2 0:3 1 2 6-3 1 3 1.1 0AYHANB,KWANC.Ac o m p a r a t i v e s t u d yo f t w oa p p r o a-c h e s f o rUAVe m e r g e n c y l a n d i n gs i t e s u r f a c e t y p e e s t i m a t i o nC/P r o c e e d