专利名称：一种用于深空探测的彗星光学面信号模拟方法
技术领域：
本发明属于航天导航技术领域，尤其涉及一种用于深空探测的彗星光学面信号模 拟方法。
背景技术：
2005年美国“深度撞击号”(De印Impact)飞船发射的“撞击器”(Impactor)成功 撞击了“坦普尔1号”(Tempell)彗星，从此拉开了人类撞击彗星探索的序幕。这次撞击探 测中光学导航是撞击器能够准确击中目标的重要保障，在光学导航中，无论是导航相机的 测试，还是光学导航中数据处理方案的验证，都离不开星空及目标彗星的光学信号，在探测 器实际拍得真实的星空及目标彗星图像前，唯一获得光学信号的途径就是地面实验室环境 下的模拟仿真。目前，国内外已有多家单位成功研制了恒星模拟器，较好的满足了利用恒星的姿 态进行测量的测试实验工作，但是导航相机单纯获取恒星信息只能用于姿态确定，要确定 探测器与目标彗星的相对位置，还必须获得目标彗星的准确光学信息。在彗星距离探测器较近时，彗星中的慧核会完全暴露出来，此时探测器观测到的 是面信号，目前还没有较为有效的方法对彗星的面信号进行模拟，以获得目标彗星的准确 光学fe息。
发明内容
本发明的目的是提供一种用于深空探测的彗星光学面信号模拟方法，以获得目标 彗星的准确光学信息。为实现上述目的，本发明采用如下技术方案一种用于深空探测的彗星光学面信 号模拟方法，该方法包括以下步骤
(1)利用现有的成熟技术绘制出具有准确时空关系和亮度等级的星空背景；
(2)依照彗星的面信号光学信号模型，计算出在不同方位，不同距离上，探测器所能观 测到的亮度等级，并运用三角面拟合几何体方法对目标彗星的面信号进行重建；
(3)根据模拟器输入的时间、观测平台的位置以及导航相机的指向，生成目标彗星相应 方向的模拟信号；
(4)将投影到高性能显示器上，完成光学信号在空间和辐射亮度上的数/模转换，最终 形成模拟的光学信号。对目标彗星的面信号进行重建的具体步骤如下；
①利用三角面拟合几何体的技术来实现慧核的几何模型重建三角面拟合几何体即用 三角面无限逼近几何体的真实几何形状；
②根据彗星面信号的反射度模型生成相应的反射度纹理；
③将步骤②所得反射度纹理叠加到彗星几何模型上；
④在步骤③的基础上添加太阳光照模型。
E1" _ ^ C" -T(Z1^Z2)-SecQ'
JJ J ''J X -J _ Gd1
τ r A τ 其中，α为探测器相对彗核的天顶角；
.T(ZpZ2) = ^KiPdZ为光线通过厚度；至高度&的彗发区段时该区段彗发的光学
厚度；
Ex = E0e~T(Zl'Z2>sece为在n波长范围内太阳在行星表面产生的照度，0为此 时的太阳天顶角；
^e = 1 Oia KcID'1 Γ r [ Φ (s Z螂)-1Γ以为太阳的辐射照度，其中太阳半径
li2I
R = 6.9599 χIO8 m ,义为波长(辨〉，T0 为黑体的温度(ζ) , C1 = 3.742χIO^r-^m2
为第一辐射常数，C2 = 14388/ · K为第二辐射常数，T0 = 5900Κ, D为太阳与行星 的距离。太阳光照模型中太阳对光学探测器的简单影响模型如下式所示 je>r#不受太阳影响
受太阳影响 °在日心黄道坐标系中，探测器对太阳的观测相角θ力导航相机成像方向与太阳方 向直接的夹角，其公式如下
r P
θ = arccosf-)
ηιρΓ
其中，探测器到太阳的距离为r，敏感器成像方向为矢量P。本发明用于深空探测的彗星光学面信号模拟方法运用OpenGL三维图形库编程实 现了光学信号的模拟，对彗星的面信号运用粒子系统进行重建接着根据模拟器输入的时 间、观测平台的位置以及相机的指向，生成相应方向的模拟信号，并投影到高性能显示器 上，完成光学信号在空间和辐射亮度上的数/模转换，最终形成模拟的光学信号。该方法绘
制简单直观，效果较好，为深空探测尤其是彗星撞击探测中光学导航相机提供了可靠的动 态光学信号。
图1是彗星形态示意图2是探测器敏感器成像方向与太阳方向之间的关系图； 图3是模拟器输入参数生成光学信号的流程图； 图4是彗星面信号生成流程图； 图5是探测器实拍“坦普尔1号”面信号影像图；图6是“坦普尔1号”实拍时本发明模拟器彗星面模拟信号图。
具体实施例方式下面结合具体的实施方式对本发明做进一步介绍 1.彗星光学信号模型
如图1所示是彗星形态示意图，典型的彗星形态包括四个部分即彗核、彗发、彗尾及 彗云，其中彗核与彗发合起来称彗头，但只是对多数彗星而言。有些彗星形态并不一定四部 分都有，有些不仅没有彗尾，连彗发都没有。就是同一颗彗星，在绕日运行过程中，它的形 态也是不断变化的，在接近太阳时，逐渐出现彗发和彗尾；而远离太阳时，彗发和彗尾又逐 渐消失了。在彗星距离探测器较近时，慧核将完全暴露出来，此时探测器观测到的将是面信 号，慧核本身不发光，其入射到光学导航相机的光信号来自对太阳直接辐射能量的反射。 但是彗核表面包裹着一层大气，近于球形，也即慧发，是彗星接近太阳出彗核中蒸发出来 的气体和微小尘粒，由于这些气体分子和微小尘粒的影响，太阳光线在透过慧发时，会产生
相应的衰减，表示衰减程度的系数可以称为消光系数，若入射光亮度为L，在吸收和散射 物质密度为P的介质中经过光路式后，光亮度减弱了 ，则消光系数表示为
产令 (1) pi,,
当光线通过厚度Z1至高度Z2的慧发区段时，该区段慧发的光学厚度TXZ1,Z2)可由消 光系数来定义，即
mZj
y^i^i，Ζ··η,》—J ^ ^ydjLi (2)
当太阳天顶角为5时，光线通过的路径会产生变化，这时其透过率Tt而.’巧可表达为 TiZl^e) =(3)
又因为太阳所发出的总辐射量在空间方向上的分布是均勻的，在為 為波长范围内太 阳在行星表面产生的照度尾则可表示为
E1 =E0e~T(Zl'E2>seca (4)
其中馬如下式(5)所示，为太阳的辐射照度；到Z2为包裹慧核的大气(即慧发)的厚 度；0为慧核上相应的太阳天顶角。E1 = 10" E'q ZTj J^ ；T5 [exp (c / λΤ ) -1]"1 λ (5)
11JliJ
式(5)中，太阳半径R. =6,9599x10' m , 1为波长U— , 为黑体的温度(幻， =周=为第一辐射常数， = 4388/ .|为第二辐射常数，ΓΛ = 5900Κ (通常太阳辐射可以认为是温度为5900K的黑体辐射)￡力太阳与行星的距离。照度p的单位为 WZm1 ο彗核反射后到达探测器的辐射照度为(在垂直传感时) J ^ = ρβτβ~Τ(Ζι'Ζ2) (6)
当探测器相对彗核的天顶角力”时，探测器的接收到的辐射照度则为 Ef = pzET�-T(Zl'Z2>seca (7)
以上就是探测器与彗星较近时，慧核完全暴露的情况下，慧核的面信号光学模型，也就 是彗星彗核面信号的反射度模型。2.太阳干扰模型
在彗星面信号的模拟过程中，太阳的干扰也是不可忽视的。太阳对各类深空光学探测 器的工作都存在重大影响，因为一旦太阳进入探测器视�。岬贾绿讲馄鰿CD曝光饱和， 无法分辨出任何天体有效信号，严重时还可能导致CCD的损坏。因此，信号模拟时必须对太 阳的干扰进行建模。如图2所示，定义导航相机成像方向与太阳方向直接的夹角为观测相角3，假设在 日心黄道坐标系中，探测器到太阳的距离为r，敏感器成像方向为矢量P，那么，
Γ ·ρ
= arccQs(——)(8)
HIpI
通常，太阳对光学探测器的简单影响模型如下式所示， r,不受太阳影响
(
(9)
5 <7；受太阳影响
探测器对观测相角6都有一定限制，只有观测相角小于等于某一阈值指标？；时，太阳
才会对成像产生干扰。 3.彗星光学信号的模拟实现
如图3所示，本发明用于深空探测的彗星光学信号模拟方法以彗星的光学信号为基 �。赩C++6. 0中，运用OpenGL三维图形库编程实现了光学信号的模拟，首先利用现有的成 熟技术绘制出具有准确时空关系和亮度等级的星空背景，具体过程可描述如下在利用星 历表中的数据计算得到恒星的精确位置，以及恒星的星等值后，利用OpenGL三维图形库对 星空背景进行三维重建，恢复恒星在空间中的位置，将恒星星等线性映射到0，1区间，利 用OpenGL绘制三维点来表示恒星，利用函数指定不同亮度(颜色)来区分不同的照在星空 背景绘制完成后，依照彗星的面信号光学信号模型，对彗星的面信号进行重建，接着根据模 拟器输入的时间、观测平台的位置以及相机的指向，计算出需要显示的区域，生成相应方向 的模拟信号(该计算过程相当于将照相机指向要显示的方向拍照)，其具体过程如下根据 模拟器输入的时间，观测平台的位置和导航相机的指向计算出彗星在导航相机中的位置， 其中位置计算公式如(10):
上式中坐标单位为光年，(；动的为彗星的赤经、赤纬，单位为弧度。r为彗星到地球
的距离。利用OpenGL图形开发包，采用三角面拟合几何体技术，完成相应得彗星面信号模 拟。模拟信号生成后投影到高性能显示器上，完成光学信号在空间和辐射亮度上的数/模 转换，最终形成模拟的光学信号。慧核的几何建�？梢岳米匀惶焯灞砻娴娜７椒ǎ腔酆擞衅渥陨淼奶� 殊性，如无法得到其DEM和表面的影像数据，慧核不一定类似于球体，很可能是不规则的几 何体，采用自然天体几何维建模方法工作量将非常巨大，而得到的绘制效率和模拟效果不 一定能满足要求，因此这里采用一种三角面拟合几何体的技术来实行慧核的几何模型重建。为了进一步理解使用三角面拟合几何体的概念，本文在此使用称为顶点数组的数 据结构，用此数据结构去描述被拟合的几何体。顶点数组中每一个元素代表一个顶点，其 结构中包含有该顶点的纹理坐标(a,ν)、法向量值(《m )以及顶点坐标，每一个 三角面由三个顶点构成，在三角面结构中不存放构成三角面的三个顶点的坐标，而是存放 了三个顶点的索引值po，pi, p2。对于任意一个几何体，可用有限个三角面拟合，并配以 相应的纹理信息数据。几何体(Geometry)由三角面数组构成，并且包含该几何体的纹理名 (texturename )以及纹理物体号(texturelD )。如图4所示，几何建模完成后并不能完成面信号的模拟，这里实现面信号模拟的 方法是反射度纹理叠加光照模型，首先根据面信号的反射度模型生成相应的反射度纹理， 然后加上太阳光照模型，最终实现面信号的模拟，其具体的实现过程如下
第一步利用三角面拟合几何体的技术来实现慧核的几何模型重建。这步的作用主要是完成慧核的几何模型的构建，主要方法是用多个三角面来无限 逼近慧核的真实几何形状。逼近的程度由慧核与视点的远近来决定，慧核与视点较远时，逼 近的程度较低；当慧核与视点较近时，逼近的程度较高。第二步根据面信号的反射度模型生成相应的反射度纹理。首先，根据慧核的真实面貌给定一个基础纹理，该纹理主要反映慧核的基本面貌， 然后再利用面信号反射度模型生成各点处的反射度，然后将反射度结果叠加到基础纹理 上，得到新的反射度纹理。第三步将反射度纹理叠加到几何模型上。计算出各点处的纹理坐标，利用OpenGl提供的纹理映射机制，将得到的反射度纹 理叠加到几何模型上。第四步添加太阳光照模型，最终实现面信号的模拟根据公式(8)和(9)，当探测 器对太阳的观测相角0小于等于某一阈值指标4时，太阳会对成像产生干扰。这时，相机的 镜头会曝光过度，所以此时，面成像为一片�：拿妫坏碧讲馄鞫蕴舻墓鄄庀嘟�$大于该
8
4.本发明用于深空探测的彗星面光学信号模拟方法的实验验证
为了验证本发明提出的彗星面光学信号方法的可行性和适用性，本方法以“深空撞击” 探测任务为背景，模拟探测器在接近目标天体时对目标彗星成像信号在以本发明方法为 核心的模拟器系统中输入时间、探测器的位置、光学相机指向等参数后，由工作站生成相应 的光学信号并显示在高分辨率显示器上，其比较结果如下2005年7月4号下午1点52 分左右，深度撞击号成功撞击“坦普尔1号”彗星，图5是撞击舱在撞击发生前6分钟拍摄 的“坦普尔1号”彗星影像图，图6为相应时间的彗星模拟信号。从试验结果可以看出，本发明用于深空探测的彗星光学信号模拟方法所生成的光 学信号与深度撞击中探测器实际拍摄到的信号强度和相对位置关系正确，满足精度要求， 可以为彗星撞击实验中光学导航相机提供准确的光学模拟信号。
权利要求
一种用于深空探测的彗星光学面信号模拟方法，其特征在于，该方法包括以下步骤(1)利用现有的成熟技术绘制出具有准确时空关系和亮度等级的星空背景；(2)依照彗星的面信号光学信号模型，计算出在不同方位，不同距离上，探测器所能观测到的亮度等级，并运用三角面拟合几何体方法对目标彗星的面信号进行重建； (3)根据模拟器输入的时间、观测平台的位置以及导航相机的指向，生成目标彗星相应方向的模拟信号；(4)将投影到高性能显示器上，完成光学信号在空间和辐射亮度上的数/模转换，最终形成模拟的光学信号。
2.权利要求1所述的用于深空探测的彗星光学面信号模拟方法，其特征在于对目标 彗星的面信号进行重建的具体步骤如下；①利用三角面拟合几何体的技术来实现慧核的几何模型重建三角面拟合几何体即用 三角面无限逼近几何体的真实几何形状；②根据彗星面信号的反射度模型生成相应的反射度纹理；③将步骤②所得反射度纹理叠加到彗星几何模型上；④在步骤③的基础上添加太阳干扰模型。
3.权利要求1所述的用于深空探测的彗星光学面信号模拟方法，其特征在于彗星的 面信号光学信号模型如下 其中，《为探测器相对彗核的天顶角； 为光线通过厚度至高度^的彗发区段时该区段彗发的光学厚度； 为在;波长范围内太阳在行星表面产生的照度、为此时的太阳天顶角；民= U为太阳的辐射照度，其中太阳半径 为彼长、,、，T0 为黑体的温度 为第一辐射常数，C2 143 /￡Μ· Γ为第二辐射常数，；=5900Κ, D为太阳与行星 的距离。
4.根据权利要求1所述的用于深空探测的彗星光学面信号模拟方法，其特征在于太 阳光照模型中太阳对光学探测器的简单影响模型如下式所示 不受太阳影响 受太阳影响 °
5.权利要求3所述的用于深空探测的彗星光学面信号模拟方法，其特征在于在曰心 黄道坐标系中，探测器对太阳的观测相角4力导航相机成像方向与太阳方向直接的夹角， 其公式如下 其中，探测器到太阳的距离为Γ，敏感器成像方向为矢量P。
全文摘要
本发明提供了一种用于深空探测的彗星光学面信号模拟方法，运用OpenGL三维图形库编程实现了光学信号的模拟，对彗星的面信号运用粒子系统进行重建接着根据模拟器输入的时间、观测平台的位置以及相机的指向，生成相应方向的模拟信号，并投影到高性能显示器上，完成光学信号在空间和辐射亮度上的数/模转换，最终形成模拟的光学信号。该方法绘制简单直观，效果较好，为深空探测尤其是彗星撞击探测中光学导航相机提供了可靠的动态光学信号。
文档编号G01C21/24GK101915581SQ20101023287
公开日2010年12月15日 申请日期2010年7月21日 优先权日2010年7月21日
发明者何钰, 周杨, 孙伟, 张衡, 徐青, 施群山, 李建胜, 王栋, 蓝朝桢, 邢帅, 马东洋 申请人:中国人民解放军信息工程大学
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
