专利名称：一种动态的月面成像敏感器的目标模拟系统的制作方法
技术领域：
本发明属于航天成像式敏感器测试技术领域，本发明也适用于深空探测、行星探测等空间环境下成像式敏感器的成像性能模拟测试。
背景技术：
随着航天技术的飞速发展，以相机为代表的光学成像式敏感器因其具有信息量大，运用灵活、效率高等特点，在深空探测领域将得到越来越广泛的应用。为确保成像式敏感器在轨可靠工作，需在地面对其各方面功能性能进行充分测试，尤其是其在轨工作环境下的光学成像性能验证。目前，国内外对于空间环境下成像式敏感器光学头部成像性能验证方式主要有两种 I、一是构建模拟场景，对真实景物进行拍摄。对于月面环境下成像式敏感器的模拟成像测试，采用本方式，需要采用灯阵模拟月面光照条件，并构建月面景物模拟试验场，时间长，造价高；同时，灯阵均匀性设计、热设计，月面景物纹理模拟等均为技术难点；此夕卜，若拍摄场景变化需重新测量定位，而且用于成像的月面景物的实际参数需要依靠外测获得，工作量大且精度有限，不利于对成像效果进行定量分析。2、一是采用目标模拟装置，生成目标图像并投影。目前已有多种敏感器的目标模拟装置1)如星模拟器(一种静态多光路星模拟器CN200810057344. 0，包括法兰和至少四个星光管，每个星光管均能够产生一颗模拟星，每个星光管的光谱和亮度均可单独调整；星光管通过角度调整机构安装于法兰上，且任意两个星光管的光轴不平行，任意两个星光管之间的相对角度可在不小于10度的范围内进行调整，可以产生大量不同构型的模拟星座，且每个星点的亮度和光谱特性都可以不相同的，任意两个模拟星座之间的几何相似程度比较低，可以很好地模拟真实星空的物理特性)(一种高精度星模拟器CN200610104766. X，包括光源和依次设置在光路上的红外反射镜、毛玻璃、聚光镜A、胶合棱镜、衰减片组、滤光片、聚光镜B、星点板和平行光管；为了使光源稳定，本发明还可包括光源亮度控制电路、设置在聚光镜A之后的胶合棱镜以及设置胶合棱镜透射光路上的光电二极管)，2)如太阳模拟器(太阳模拟器及其运转方法CN201010214004. 1，包括多个氙灯(41 4η);对所述各氙灯(41 4η)配备的多个光量传感器(SI Sn);和对所述各氙灯(41 4η)配备的、用于控制流经该氙灯(41 4η)的电流或对该氙灯施加的电压多个控制电路(7)，使基于所述各光量传感器(SI Sn)的检测信号反馈到所述各控制电路(7)，控制该控制电路(7)，来控制所述各氙灯(41 4η)的光量。由此提供太阳模拟器，在具备多个氙灯作为电源的太阳模拟器中，能够通过各氙灯稳定地获得所期望的光量，且能够使有效照射面中的照度均匀化)。这些设计均用于模拟行星及太阳等单个或多个点光源的入射及入射角度变化，敏感器本身对成像结果有系统的处理算法(去背景阈值、提取中心点等)，因此设计中不考虑点光源成像后的分辨率(像素个数)，及像源灰度的精确控制输出。而相机等成像式敏感器不进行类似星敏感器、太阳敏感器等的处理算法，其目标模拟装置设计中必须考虑像源的分辨率，及灰度的精确控制输出，以达到模拟景物及纹理的要求。而且，上述这些系统设计远不能够满足鱼眼相机等超广角成像式敏感器的成像视角要求
发明内容
本发明的技术解决问题是克服现有技术的不足，提供一种动态的月面成像敏感器的目标模拟系统。本发明的技术解决方案是一种动态的月面成像敏感器的目标模拟系统，包括目标影像生成单元、影像输出单元和光学传递单元；目标影像生成单元，根据DEM数据和遥感影像数据仿真生成三维月面形貌，并根据月面成像敏感器成像时刻的月面光照条件以及该时刻月面成像敏感器的位置姿态参数动态生成符合月面成像敏感器成像要求的灰度图像数据；影像输出单元将目标影像生成单元生成的灰度图像数据在LCD器件上进行显示；光学传递单元将IXD器件上显示的影像进行光学变换，使所述影像正好成像在月面成像敏感器的成像器件上。所述的仿真生成三维月面形貌过程如下(I)将DEM数据和遥感影像数据进行三维场景的叠加，并将叠加后的三维场景进行剔除处理，剔除场景中的突兀点，得到粗分辨率的月面地貌；(2)对步骤(I)中得到的粗分辨率的月面形貌进行地形分形计算，生成高分辨率的月面地貌；(3)根据模拟测试需要，建立岩石、陨石坑的三维模型，并根据月面的土壤特征，构建月面土壤模型；(4)将岩石、陨石坑的三维模型以及构建的月面土壤模型根据月面形态特征分布在步骤(2)生成的高分辨率的月面地貌上，生成三维月面形貌。所述的生成符合月面成像敏感器成像的灰度图像数据过程如下(1-1)根据月面成像敏感器成像时刻的位置姿态，对生成的三维月面形貌进行中心投影成像；(1-2)通过针孔成像和畸变处理得到符合月面成像敏感器成像要求的灰度图像数据。本发明与现有技术相比有益效果为(I)本发明以DEM数据和遥感影像数据为基础，仿真生成三维月面地貌，并能够根据需要进行景物(岩石、陨石坑)设置，生成特定的三维月面形貌灰度数据，满足月面成像敏感器对多种月面形貌的成像需求，进而可验证图像匹配、障碍识别及路径规划算法的有效性，避免了构建月面模拟场景的花费与耗时。(2)本发明以DEM数据和遥感影像数据为基础构建基本三维月面形貌，并能够根据面成像敏感器成像时刻的位置姿态等参数进行解算，可动态生成符合月面成像敏感器成像要求的灰度图像数据，并进行显示输出及光学变换，生成符合月面成像敏感器光路的虚拟图像，使影像以期望的光路形式成像在月面成像敏感器的成像器件上，在没有实物的情况下实现月面成像敏感器的物理闭环测试。
(3)本发明以高分辨率、多灰度的显示器件作为像源，以满足成像式敏感器成像分辨率及灰度要求。
图I为本发明组成框图；图2为中点分形示意图；图3为旋转曲线形状示意图；图4为简易的柯拉照明结构示意图； 图5为本发明光学传递单元光路结构示意图。
具体实施方式
下面结合附图对本发明做详细的说明，如图I所示，本发明一种动态的月面成像敏感器的目标模拟系统，包括目标影像生成单元、影像输出单元和光学传递单元；月面成像敏感器可以采用APS相机实现。目标影像生成单元，根据DEM(数字高程模型Digital Elevation Model)数据和遥感影像数据仿真生成三维月面形貌，并根据月面成像敏感器成像时刻的月面光照条件以及该时刻月面成像敏感器的位置姿态参数动态生成符合月面成像敏感器成像的灰度图像数据；影像输出单元将目标影像生成单元生成的灰度图像数据在IXD器件上进行显示；光学传递单元将IXD器件上显示的影像进行光学变换，使所述影像正好成像在月面成像敏感器的成像器件上。下面分别对上述三部分进行详细说明。(一)目标影像生成单元目标影像生成单元主要实现两方面的功能，I三维月面形貌的仿真，2成像敏感器成像仿真。I、三维月面形貌的仿真(I)以DEM数据和遥感影像数据为基础，叠加遥感影像纹理，得到粗分辨率的月面地貌；由于DEM数据三维高程数据，而遥感影像数据为纹理数据，因此需要对二者进行叠加处理，即根据生成的随机位置，本例中主要是利用OpenGL实现三维场景的叠加。叠加技术是利用OpenGL接口函数实现，属于公知技术。叠加后的三维场景要进行剔除处理，应用OPENGL中的ALPHA(透明度)来过滤掉纹理的背景、视见体裁剪、地球遮挡剔除、屏幕空间剔除等剔除，调用OpenGL相关函数，设置其参数可实现相关操作；剔除场景中的突兀点，得到粗分辨率的月面地貌；(2)对步骤(I)中得到的粗分辨率的月面形貌进行地形分形计算，生成高分辨率的月面地貌。地形分形计算是将分形几何与分数维布朗运动数学模型相结合，常采用的方法有泊松阶跃法(poisson filtering)、傅立叶滤波法(Fourier filtering)、中点位移法(midpoint displacement)、逐次随机增加法(sue cessive random additions)和带限噪声累积法(summing band limited noises)等 5 类。在上述几种三维分形地形建模方法中，以用于快速地景生成的中点位移法应用最为广泛，它的特点是简洁而快速。本例中拟采用基于中点位移法的Diamond-Square算法来生成高分辨率的月面地貌。Diamond-Square是一种较为常用的中点随机位移算法。这种算法最初是由Fournier, Fussell和Carprnter提出的。如图2所示,取四个点的正方形,并对四个角指定一个高度.然后在每边的中点分割地平面得5个新网格点位置。地面边上的点、f、g、h上的高度可以按最近两个顶点的平均高度加上一随机偏移来计算。如中点e处的高度可用顶点a和b来计算该算法是从由种子点组成的正方形开始，通过采用若干次随机中点位移方法，不断细分，最后获得逼真的三维地形仿真图.具体描述如下初始化二维数组ηΧη，(η为2的整数次幂加I)并将4个角设为相同高度。假定 用一个5X5的数组。图2中，左侧图a的四个角种上了初始高度值，表示为黑点。这是递归细分过程的起点，该过程分两步diamond步取四个点的正方形，在正方形中点生成一个随机值，中点为两对角线交点。中点值是平均四个角值再加上一个随机量计算得到的。这样就得到了一个棱锥。当网格上分布着多个正方形时有点象钻石。square步取每个四点形成的棱锥,在棱锥的中心生成一个随机值。平均角值再加上与diamond步相同的随机量,计算出每条边中点值。这又给你一个正方形。这样,如果已经生成了一个种子正方形并经过单独一次细分过程将得到四个方形。第二次经过该过程得到16个方形，第三次得到64个方形。正方形的个数是2(2+i)，其中i是迭代的次数。Diamond-Square可以用迭代或递归的方法来实现，迭代法的伪代码为
当Square边长大于O时 { 遍历数组，对每个正方形表达执行Diamond步； 遍历数组，对每个棱锥表达执行Square步；
减少随机数范围 }(3)根据模拟测试需要，建立岩石、陨石坑的三维模型，并根据月面的土壤特征，构建月面土壤模型；(3. I)陨石坑的建模由于简单陨石坑形状规则，于是可以采用旋转面的方法生成陨石坑，假设陨石坑中心坐标为(Xl，yi)，点(x，y)的高程偏差我们采用如下计算公式Mi = k-[{d21 rAI r2)(I)其中d= (X-Xl)2+(y-yi)2，也就是离陨石坑中心的距离的平方，r为陨石坑的半径，k是控制参数，用来控制坑的深度。通过上面的高程偏差公式，绘制出曲线如图3 分析上式，当d = 0时，Ah = k· r/4,也就是说在陨石坑的中心位置的深度H =k · r/4 ;在d = r2/2和d = r2时,坑的深度为O,中间形成陨石坑的边缘突起,通过上式绘制出的形状比较符合简单陨石坑的形状和Pike提出的经验公式，而且陨石坑的深度可以根据实际情况进行修正，方便实际应用。在构建好单个简单陨石坑模型之后，使用随机函数将月球陨石坑随机分布在月面地形上。为了尽量符合月球陨石坑的真实分布情况，加入坑密度、深度、大小等参数对陨石坑的数量和形状进行控制，可以根据需要生成疏密不同的陨石坑分布。模拟月球陨石坑算法流程如下首先将要加入陨石坑的原始地形归一化，使之横纵坐标在O到I之间； 然后根据地形的大小生成特定数量的陨石坑并使其随机分布；最后根据陨石坑旋转方程(公式I)遍历整个地形，计算陨石坑内的高程偏差，将计算的偏差和原来的地形融合，这样就生成了带陨石坑的新地形。(3. 2)岩石(石块)的建模根据月表石块数的统计信息，随机生成相应的大小不等的石块。月表每100m2面积范围内的石块数分布如下高度25彡h > 6厘米的石块数100块50彡h > 25厘米的石块数为3-4块；h > 50厘米的石块数为O. 6块。月球石块的最短尺寸与最长尺寸的比值在I : 5至I : I范围内，石块表面可能是圆形、矩形、凹坑形等，建模时采用3DMAX实现可随机确定各种表面类型的石块所占总体石块数的比例。(3.3)月面土壤的构建月面土壤模型的构建主要是根据月面的土壤特征，在3DMAX中采用“小沙子”纹理重复贴图实现月面土壤的构建。(4)通过OpenGL的正射投影方式，将岩石、陨石坑的三维模型以及构建的月面土壤模型根据月面形态特征分布在步骤(2)生成的高分辨率的月面地貌上并进行高程融合，最终生成三维月面形貌的仿真效果。利用函数gl0rtho(l，r, b，t，n, f)完成正射投影变换。函数参数为平行投影的左右、上下和远近裁剪面。利用函数glCopyTexlmage2D(GL_TEXTURE_2D，0，GL_DEPTH_COMPONENT, O, O, demSizeX, demSizeY，0拷贝深度缓冲中的深度信息到纹理，该纹理经过输出并坐标规划后即为叠加陨石坑和岩石的高程图。基于OpenGL三维绘制的月面形(岩石、撞击坑)貌(DEM)高程数据融合算法如下glMatrixMode(GL_PROJECTION);
glLoadldentity();
//设置正射投影方式glOrtho(l,r,b，t，n,f);
glGetFloatv(GL_PROJECTION_MATRIX, IightProjection);
//将视点切换到某个高度位置 glMatrixMode(GL_MODELVIEW); glLoad 丨 dentity();
gluLookAt(0.0f, Height, O.Of, O.Of, O.Of, O.Of, O.Of, 1.0f, O.Of);glGetFloatv(GL_MODELVIEW_MATRIX, IightModeIview);
//设置生成的融合高程数据的大小(demSizeX, demSizeY)glViewport(0, 0, demSizeX, demSizeY);
Il清除深度缓冲
glClear(GL_DEPTH_BUFFER_BIT);
//关闭光照、材质等开关glShadeModel(GL_FLAT); glDisable(GL_LIGHTING); glDisable(GL_COLOR_MATERIAL); glDisable(GL_NORMALIZE); glColorMask(0, 0, 0，0);
//设置多边形偏移，提高绘制精度 glEnable(GL_POLYGON_OFFSET_FILL);
//激活并绑定纹理
glActiveTextureARB(GL_TEXTURE1_ARB);
glBindTexture(GL_TEXTURE_2D,shadowTexturelD);
Il绘制叠加了撞击坑的地形数据 DrawTerrain();
//绘制岩石模型 DrawStoneModels();
//拷贝深度缓冲中的深度信息到纹理，该纹理经过输出并坐标规划后即为叠加撞击坑和岩石的高程图
glCopyTexlmage2D(GL_TEXTURE_2D, O, GL_DEPTH_COMPONENT, O, O, demSizeX, demSizeY, 0);
//恢复显示开关
glShadeModel(GL_SMOOTH);
glEnable(GL_LIGHTING);
glEnable(GL_COLOR_MATERIAL);
glEnable(GL_NORMALIZE);
glColorMask(1, 1，1, 1);
glDisable(GL_POLYGON_OFFSET_FILL);
glActiveTextureARB(GL_TEXTUREO_ARB);2、月面成像敏感器成像仿真(1-1)根据月面成像敏感器成像时刻的位置姿态以及视场角，对仿真生成的三维月面形貌进行中心投影成像；在OpenGL中，模型空间坐标系一般可以定义计算机屏幕中心点为坐标系的原点，人正对屏幕向右为X轴正方向，向上为Y轴正方向，垂直屏幕指向外面的方向为Z轴正方向。在模型空间的三维物体可以通过OpenGL的模型函数进行一系列的平移、旋转和缩放，以确定其大小、位置和形状。在OpenGL中，三维物体的最终定位是通过基本实物模型与模型变换矩阵相乘产生的，即[X y 1]τ = V · P · M · [X Y Z 1]τ。X、Y、Z分别代表三维月面形貌在模型空间坐标系下的点；x、y代表三维月面形貌经中心投影后的屏幕坐标；M、P、V分别代表模型矩阵、投影矩阵和放射变换矩阵；I)定义近平面的左下角和右上角坐标分别为(Xl，yb)和(\，yt)，定义近远平面为Zn 和 Zf。
权利要求
1.一种动态的月面成像敏感器的目标模拟系统，其特征在于包括目标影像生成单元、影像输出单元和光学传递单元； 目标影像生成单元，根据DEM数据和遥感影像数据仿真生成三维月面形貌，并根据月面成像敏感器成像时刻的月面光照条件以及该时刻月面成像敏感器的位置姿态参数动态生成符合月面成像敏感器成像要求的灰度图像数据； 影像输出单元将目标影像生成单元生成的灰度图像数据在LCD器件上进行显示； 光学传递单元将LCD器件上显示的影像进行光学变换，使所述影像正好成像在月面成像敏感器的成像器件上。
2.根据权利要求I所述的一种动态的月面成像敏感器的目标模拟系统，其特征在于所述的仿真生成三维月面形貌过程如下 (1)将DEM数据和遥感影像数据进行三维场景的叠加，并将叠加后的三维场景进行剔除处理，剔除场景中的突兀点，得到粗分辨率的月面地貌； (2)对步骤(I)中得到的粗分辨率的月面形貌进行地形分形计算，生成高分辨率的月面地貌； (3)根据模拟测试需要，建立岩石、陨石坑的三维模型，并根据月面的土壤特征，构建月面土壤模型； (4)将岩石、陨石坑的三维模型以及构建的月面土壤模型根据月面形态特征分布在步骤(2)生成的高分辨率的月面地貌上，生成三维月面形貌。
3.根据权利要求I所述的一种动态的月面成像敏感器的目标模拟系统，其特征在于所述的生成符合月面成像敏感器成像的灰度图像数据过程如下 (1-1)根据月面成像敏感器成像时刻的位置姿态，对生成的三维月面形貌进行中心投影成像； (1-2)通过针孔成像和畸变处理得到符合月面成像敏感器成像要求的灰度图像数据。
全文摘要
一种动态的月面成像敏感器的目标模拟系统，包括目标影像生成单元、影像输出单元和光学传递单元；目标影像生成单元，根据DEM数据和遥感影像数据仿真生成三维月面形貌，并根据月面成像敏感器成像时刻的月面光照条件以及该时刻月面成像敏感器的位置姿态参数动态生成符合月面成像敏感器成像要求的灰度图像数据；影像输出单元将目标影像生成单元生成的灰度图像数据在LCD器件上进行显示；光学传递单元将LCD器件上显示的影像进行光学变换，使所述影像正好成像在月面成像敏感器的成像器件上。
文档编号G01M11/02GK102928201SQ20121041394
公开日2013年2月13日 申请日期2012年10月24日 优先权日2012年10月24日
发明者张志 , 杨洁, 张晋, 李志平, 王磊, 余志鸿 申请人:北京控制工程研究所