一种火星轨道相机影像仿真方法及系统

1.本发明属于行星摄影测量与遥感领域，特别是涉及一种火星轨道相机影像仿真方法及系统。


背景技术：

2.火星轨道相机影像可以直观的帮助人们了解火星表面形貌特征；利用火星相机立体轨道影像数据产生的全球或局部地形信息可用于火星探测器的着陆选址、精确着陆、着陆器巡视导航任务；通过分析火星形貌特点、特征，可以反映火星地质构造的演化过程、寻找水、火星是否存在生命的证据。火星轨道相机的几何设计及验证至关重要，而火星轨道相机影像仿真是相机验证的重要手段，对优化相机设计参数提高火星立体影像几何定位精度具有重要意义。火星表面独特的形貌特征以及与地球表面相比所具有的弱纹理的特点，使得常规的相机仿真方法难以顾及火星表面的形貌特征、特点。
3.现有的火星轨道器相机成像仿真技术，都是基于地球表面的影像和地形或模拟影像和模拟地形进行仿真，这些方法不能够真实、准确的反应火星表面的形貌特征。
4.本发明提出利用火星全球dem和火星全球数字影像模型(mdim)作为仿真数据源，设计相应的相机影像仿真方法，充分顾及了火星实际形貌特征、特点，能够有效的满足火星轨道相机设计验证需求。


技术实现要素：

5.本发明针对现有技术的不足，提供一种火星轨道相机影像仿真方案。
6.本发明提供的技术方案是一种火星轨道相机影像仿真方法，包括如下步骤：
7.步骤1，输入仿真数据源，包括火星全球dem和火星全球数字影像模型mdim；
8.步骤2，仿真火星轨道相机的位置和姿态参数，包括利用三次多项式建立火星轨道器姿轨模型，并顾及相机安装矩阵，计算曝光时刻相机的位置和姿态参数；
9.步骤3，建立火星轨道相机成像模型，包括根据步骤2计算的曝光时刻相机的位置和姿态参数，建立火星轨道相机成像模型，进一步推导出相机视向量模型；
10.步骤4，计算仿真影像平面上像素对应的物方坐标，包括利用步骤3的相机成像模型，将仿真影像平面上像素逐一投影到火星全球dem表面，通过相机视向量模型计算视向量与火星全球dem表面的交点坐标，获得对应物方点的位置；
11.步骤5，根据步骤4计算的物方点位置，在火星全球影像mdim上内插计算像素的灰度值；
12.步骤6，将步骤5所得火星全球影像mdim上对应像素的灰度值逐一赋值给仿真影像平面上像素。
13.而且，步骤2中，利用以下三次多项式建立火星轨道器姿轨模型，
[0014][0015][0016]
其中，xs(t)、ys(t)、zs(t)表示t时刻轨道器在火星体固坐标系中的位置；ω(t)、κ(t)表示t时刻轨道器在火星体固坐标系中的姿态；ai、bi、ci、di、ei、fi分别表示多项式模型的系数，i＝0,1,2,3；σ
ω
(t)、σ
φ
(t)、σ
κ
(t)分别为位置和姿态随机噪声；
[0017]
设相机安装矩阵为在时刻t轨道器到火星体固坐标系的旋转矩阵为则计算如下，
[0018][0019]
其中，为在时刻t火星轨道相机坐标系到火星体固坐标系的旋转矩阵。
[0020]
而且，步骤3中，建立火星轨道相机成像模型的实现方式如下，
[0021]
根据步骤2中计算的曝光时刻相机的位置和姿态参数，建立t时刻火星轨道相机成像模型，如下，
[0022][0023]
其中，表示火星表面物方点在火星体固坐标系中的坐标；表示火星轨道相机投影中心在火星体固坐标系中的坐标；表示时刻t火星轨道相机坐标系到火星体固坐标系旋转矩阵；λ为比例系数；x,y为仿真影像平面上像元的平面坐标；f为相机主距；
[0024]
设表示t时刻火星轨道相机投影中心位置矢量，表示t时刻火星轨道相机视向量；设有物方点m(x,y,z)，视向量则式(4)进一步表示为式(5)所示相机视向量模型，
[0025][0026]
其中，表示火星表面上m点的向量，μ表示比例系数，分别表示视向量的三个坐标分量。
[0027]
而且，步骤4中，计算仿真影像平面上像素对应的物方坐标的实现方式如下，
[0028]
设火星的椭球方程为：
[0029][0030][0031]
其中，a、b分别表示初始的火星椭球体的赤道半径和极半径；h表示椭球高；a、b分别表示迭代过程中火星椭球体的新的赤道半径和新的极半径；
[0032]
将式(5)代入式(6)得到下列方程式：
[0033][0034]
执行如下步骤，
[0035]
(1)计算视向量与火星椭球面的初始交点坐标
[0036]
令h＝0，解算式(8)，取比例系数μ＝min{μ1,μ2}并代入式(5)，计算交点坐标mk(xk,yk,zk)，下标k表示迭代次数，第一次执行步骤(1)时k＝0；μ1、μ2分别表示式(8)的两个根；
[0037]
(2)将第k次计算的用直角坐标表示的交点坐标mk(xk,yk,zk)转换为经、纬度和椭球高表示的交点坐标根据计算的经、纬度坐标对火星全球dem进行双线性内插，得到新的椭球高hk；
[0038]
(3)重新计算视向量与火星椭球面的交点坐标，
[0039]
令h＝hk，解算式(8)，取μ
min
＝min{μ1,μ2}并代入式(5)，计算新的交点坐标m
k+1
(x
k+1
,y
k+1
,z
k+1
)，令k＝k+1；
[0040]
(4)根据影像分辨率相应的预设阈值ρ，重复(2)、(3)，直至||mkm
k+1
||＜ρ，此时交点收敛；视向量与全球dem交点的直角坐标m(x,y,z)为仿真影像平面上像素对应的物方坐标，转换为经、纬度和椭球高表示的大地坐标作为步骤4最终的结果。
[0041]
而且，步骤5中，根据步骤4计算的仿真影像平面上像素对应的物方点的经、纬度坐标对火星全球影像mdim进行双线性内插，计算对应像素的灰度值。
[0042]
另一方面本发明提供一种火星轨道相机影像仿真系统，用于实现如上所述的一种火星轨道相机影像仿真方法。
[0043]
而且，包括以下模块，
[0044]
第一模块，用于输入仿真数据源，包括火星全球dem和火星全球数字影像模型mdim；
[0045]
第二模块，用于仿真火星轨道相机的位置和姿态参数，包括利用三次多项式建立
火星轨道器姿轨模型，并顾及相机安装矩阵，计算曝光时刻相机的位置和姿态参数；
[0046]
第三模块，用于建立火星轨道相机成像模型，包括根据第二模块计算的曝光时刻相机的位置和姿态参数，建立火星轨道相机成像模型，进一步推导出相机视向量模型；
[0047]
第四模块，用于计算仿真影像平面上像素对应的物方坐标，包括利用第三模块的相机成像模型，将仿真影像平面上像素逐一投影到火星全球dem表面，通过相机视向量模型计算视向量与火星全球dem表面的交点坐标，获得对应物方点的位置；
[0048]
第五模块，用于根据第四模块计算的物方点位置，在火星全球影像mdim上内插计算像素的灰度值；
[0049]
第六模块，用于将第五模块所得火星全球影像mdim上对应像素的灰度值逐一赋值给仿真影像平面上像素。
[0050]
或者，包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如上所述的一种火星轨道相机影像仿真方法。
[0051]
或者，包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如上所述的一种火星轨道相机影像仿真方法。
[0052]
通过以上技术方案，本发明利用火星全球dem和火星全球数字影像模型(mdim)作为仿真数据源；根据火星轨道相机的位置和姿态参数，建立轨道相机模型；利用像方到火星全球dem的投影过程，计算对应物方点的坐标；火星全球影像上内插计算对应像素的灰度值；灰度赋值完成相机影像仿真。
[0053]
与现有技术相比，本发明的有益效果是：
[0054]
火星表面形貌特征丰富，火星表面有陨击坑和盆地、大的盾形火山、峡谷系、千涸河床、崩塌地貌、沙丘、极区沉积层等多种类型形貌，且与地球表面相比具有弱纹理的特点。本发明的优势包括如下四个方面：1.利用火星表面的真实影像及形貌特征生成不同分辨率的仿真影像，能够准确的反应火星表面的形貌特征；2.适用于火星轨道的框幅式相机和线阵推扫式相机等所有类型的成像相机影像的仿真；3.利用火星全球的dem和全球的数字影像，可以根据不同的任务需求，设计火星轨道相机的位置和姿态参数，针对不同火星形貌特征的区域进行相机影像仿真；4.针对不同的相机设计参数、影像的航向旁向重叠度进行影像仿真。本发明所得仿真影像可用于火星轨道相机的几何设计及验证、火星立体影像几何定位分析等。
附图说明
[0055]
图1为本发明实施例的“火星全球dem(mars mgs mola-mex hrsc blended dem global 200m v2)”示意图。
[0056]
图2为本发明实施例的“火星全球数字影像模型mdim”示意图。
[0057]
图3为本发明实施例的“火星坐标系(mars body-fixed frame)”示意图。
[0058]
图4为本发明实施例的“火星轨道相机视向量与火星椭球体”示意图。
[0059]
图5为本发明实施例的“火星全球dem双线性内插”示意图。
具体实施方式
[0060]
下面结合附图和实施例对本发明的技术方案作进一步说明。
[0061]
本发明提供一种火星轨道相机影像仿真方案，利用火星全球dem和火星全球数字影像模型(mdim)作为仿真数据源；根据火星轨道相机的位置和姿态参数，建立轨道相机模型；将仿真影像平面上像元投影至火星全球dem表面，计算对应物方点的坐标；在火星全球影像上内插计算对应像素的灰度值；灰度赋值完成相机影像仿真。
[0062]
实施例提供一种火星轨道相机影像仿真方法，包括如下步骤：
[0063]
步骤1，输入仿真数据源，包括火星全球dem和火星全球数字影像模型mdim。
[0064]
和现有的火星轨道器相机成像仿真技术，基于地球表面的影像和地形或模拟影像和模拟地形进行仿真不同。本发明提出利用火星全球dem和火星全球数字影像模型mdim作为仿真数据源。
[0065]
实施例中，采用火星全球dem(如图1)，采用200m分辨率的dem(mars mgs mola-mex hrsc blended dem global 200m v2)；火星全球数字影像模型mdim(如图2)，采用232m分辨率的数字影像模型(viking mdim2.1 grayscale global mosaic 232m)。火星全球dem和火星全球数字影像模型mdim均采用simple cylindrical地图投影模式。
[0066]
步骤2，仿真设计火星轨道相机的位置和姿态参数。
[0067]
实施例中，利用三次多项式建立火星轨道器姿轨模型，并顾及相机安装矩阵，计算曝光时刻相机的位置和姿态参数，实现如下：
[0068]
如图3为火星体固坐标系(mars body-fixed frame)。火星体固坐标系o-xyz由国际天文联合会和国际大地测量协会iau/iag定义。其中坐标原点o位于火星质心；x轴位于火星平均赤道平面与火星主子午面的交线，z轴为火星自转轴且指向火星北极；y轴按右手法则确定。
[0069]
利用三次多项式(1)、(2)建立火星轨道器姿轨模型，
[0070][0071][0072]
其中：xs(t)、ys(t)、zs(t)表示t时刻轨道器在火星体固坐标系中的位置；ω(t)、κ(t)表示t时刻轨道器在火星体固坐标系中的姿态；ai、bi、ci、di、ei、fi(i＝0,1,2,3)分别表示多项式模型的系数，考虑姿轨模型中的随机噪声的影响，给定任意曝光时刻t按照式(1)(2)可计算轨道器的位置和姿态；
[0073]
σ
ω
(t)、σ
φ
(t)、σ
κ
(t)分别为位置和姿态随机噪声。
[0074]
设相机安装矩阵为在时刻t轨道器到火星体固坐标系的旋转矩阵为可由式(2)计算得到，则在时刻t火星轨道相机坐标系到火星体固坐标系的旋转矩阵计算如下：
[0075][0076]
步骤3，建立火星轨道相机成像模型。
[0077]
实施例中，包括根据步骤2计算的曝光时刻相机的位置和姿态参数，建立火星轨道相机成像模型，进一步推导出相机视向量模型。
[0078]
根据步骤2中式(1)和式(3)计算的曝光时刻相机的位置和姿态参数，建立t时刻火星轨道相机成像模型，如式(4)：
[0079][0080]
其中，表示火星表面物方点在火星体固坐标系中的坐标；表示火星轨道相机投影中心在火星体固坐标系中的坐标；表示时刻t火星轨道相机坐标系到火星体固坐标系旋转矩阵；λ为比例系数；
[0081]
x,y为仿真影像平面(即本发明要生成的火星轨道相机影像)上像元的平面坐标；f为相机主距。
[0082]
如图4所示，o-xyz表示火星体固坐标系；表示t时刻火星轨道相机投影中心位置矢量；表示t时刻火星轨道相机视向量。设有物方点m(x,y,z)，视向量则式(4)可进一步表示为式(5)所示相机视向量模型：
[0083][0084]
其中，表示火星表面上m点的向量；μ表示比例系数；分别表示视向量的三个坐标分量。
[0085]
步骤4，计算仿真影像平面上像素对应的物方坐标
[0086]
利用步骤3的相机成像模型，将仿真影像平面上像素逐一投影到火星全球dem表面，计算对应物方点坐标，实际上就是根据式(5)计算视向量与火星全球dem表面的交点坐标。即通过计算视向量与火星全球dem表面的交点坐标，获得对应物方点的位置。
[0087]
如图4所示，具体计算过程如下：
[0088]
设火星的椭球方程为：
[0089]
[0090]
其中
[0091]
其中，a、b分别表示初始的火星椭球体的赤道半径和极半径；h表示椭球高；a、b分别表示迭代过程中火星椭球体的新的赤道半径和新的极半径。
[0092]
将式(5)代入式(6)得到下列方程式(8)：
[0093][0094]
执行如下步骤，
[0095]
(5)计算视向量与火星椭球面的初始交点坐标
[0096]
令h＝0，解算式(8)，取μ＝min{μ1,μ2}并代入式(5)，计算交点坐标mk(xk,yk,zk)，下标k表示迭代次数，第一次执行步骤(1)时k＝0；μ1、μ2分别表示式(8)的两个根；比例系数μ取两个根中的最小值。
[0097]
(6)将第k次计算的用直角坐标表示的交点坐标mk(xk,yk,zk)转换为经、纬度和椭球高表示的交点坐标如图5所示，h
00
、h
01
、h
10
、h
11
分别表示内插窗口四角高程，h表示待内插高程；δλ、表示待内插点的经、纬度差值。根据计算的经、纬度坐标对火星全球dem进行双线性内插，得到新的椭球高hk。
[0098]
(7)重新计算视向量与火星椭球面的交点坐标
[0099]
令h＝hk，解算式(8)，取μ
min
＝min{μ1,μ2}并代入式(5)，计算新的交点坐标m
k+1
(x
k+1
,y
k+1
,z
k+1
)，令k＝k+1。
[0100]
(8)根据影像分辨率设定相应阈值ρ(小于一个像素)，重复(2)、(3)，直至||mkm
k+1
||＜ρ，此时交点收敛。视向量与全球dem交点的直角坐标m(x,y,z)即仿真影像平面上像素对应的物方坐标，将其转换为经、纬度和椭球高表示的大地坐标作为步骤4最终的结果。
[0101]
步骤5，在火星全球影像mdim上内插计算像素的灰度值。
[0102]
根据步骤4计算的仿真影像平面上像素对应的物方点的经、纬度坐标采用与火星全球dem类似的内插方法，对火星全球影像进行双线性内插，计算对应像素的灰度值。
[0103]
步骤6，仿真影像平面上像素灰度赋值。
[0104]
将步骤5计算的火星全球影像上对应像素的灰度值逐一赋值给对应的仿真影像平面上像素。
[0105]
具体实施时，本发明技术方案提出的方法可由本领域技术人员采用计算机软件技术实现自动运行流程，实现方法的系统装置例如存储本发明技术方案相应计算机程序的计算机可读存储介质以及包括运行相应计算机程序的计算机设备，也应当在本发明的保护范围内。
[0106]
在一些可能的实施例中，提供一种火星轨道相机影像仿真系统，包括以下模块，
[0107]
第一模块，用于输入仿真数据源，包括火星全球dem和火星全球数字影像模型mdim；
[0108]
第二模块，用于仿真火星轨道相机的位置和姿态参数，包括利用三次多项式建立火星轨道器姿轨模型，并顾及相机安装矩阵，计算曝光时刻相机的位置和姿态参数；
[0109]
第三模块，用于建立火星轨道相机成像模型，包括根据第二模块计算的曝光时刻相机的位置和姿态参数，建立火星轨道相机成像模型，进一步推导出相机视向量模型；
[0110]
第四模块，用于计算仿真影像平面上像素对应的物方坐标，包括利用第三模块的相机成像模型，将仿真影像平面上像素逐一投影到火星全球dem表面，通过相机视向量模型计算视向量与火星全球dem表面的交点坐标，获得对应物方点的位置；
[0111]
第五模块，用于根据第四模块计算的物方点位置，在火星全球影像mdim上内插计算像素的灰度值；
[0112]
第六模块，用于将第五模块所得火星全球影像mdim上对应像素的灰度值逐一赋值给仿真影像平面上像素。
[0113]
在一些可能的实施例中，提供一种火星轨道相机影像仿真系统，包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如上所述的一种火星轨道相机影像仿真方法。
[0114]
在一些可能的实施例中，提供一种火星轨道相机影像仿真系统，包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如上所述的一种火星轨道相机影像仿真方法。
[0115]
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

技术特征：
1.一种火星轨道相机影像仿真方法，其特征在于，包括如下步骤：步骤1，输入仿真数据源，包括火星全球dem和火星全球数字影像模型mdim；步骤2，仿真火星轨道相机的位置和姿态参数，包括利用三次多项式建立火星轨道器姿轨模型，并顾及相机安装矩阵，计算曝光时刻相机的位置和姿态参数；步骤3，建立火星轨道相机成像模型，包括根据步骤2计算的曝光时刻相机的位置和姿态参数，建立火星轨道相机成像模型，进一步推导出相机视向量模型；步骤4，计算仿真影像平面上像素对应的物方坐标，包括利用步骤3的相机成像模型，将仿真影像平面上像素逐一投影到火星全球dem表面，通过相机视向量模型计算视向量与火星全球dem表面的交点坐标，获得对应物方点的位置；步骤5，根据步骤4计算的物方点位置，在火星全球影像mdim上内插计算像素的灰度值；步骤6，将步骤5所得火星全球影像mdim上对应像素的灰度值逐一赋值给仿真影像平面上像素。2.根据权利要求1所述火星轨道相机影像仿真方法，其特征在于：步骤2中，利用以下三次多项式建立火星轨道器姿轨模型，次多项式建立火星轨道器姿轨模型，其中，x
s
(t)、y
s
(t)、z
s
(t)表示t时刻轨道器在火星体固坐标系中的位置；ω(t)、κ(t)表示t时刻轨道器在火星体固坐标系中的姿态；a
i
、b
i
、c
i
、d
i
、e
i
、f
i
分别表示多项式模型的系数，i＝0,1,2,3；σ
ω
(t)、σ
φ
(t)、σ
κ
(t)分别为位置和姿态随机噪声；设相机安装矩阵为在时刻t轨道器到火星体固坐标系的旋转矩阵为则计算如下，其中，为在时刻t火星轨道相机坐标系到火星体固坐标系的旋转矩阵。3.根据权利要求1所述火星轨道相机影像仿真方法，其特征在于：步骤3中，建立火星轨道相机成像模型的实现方式如下，根据步骤2中计算的曝光时刻相机的位置和姿态参数，建立t时刻火星轨道相机成像模型，如下，
其中，表示火星表面物方点在火星体固坐标系中的坐标；表示火星轨道相机投影中心在火星体固坐标系中的坐标；表示时刻t火星轨道相机坐标系到火星体固坐标系旋转矩阵；λ为比例系数；x,y为仿真影像平面上像元的平面坐标；f为相机主距；设表示t时刻火星轨道相机投影中心位置矢量，表示t时刻火星轨道相机视向量；设有物方点m(x,y,z)，视向量则式(4)进一步表示为式(5)所示相机视向量模型，其中，表示火星表面上m点的向量，μ表示比例系数，分别表示视向量的三个坐标分量。4.根据权利要求1所述火星轨道相机影像仿真方法，其特征在于：步骤4中，计算仿真影像平面上像素对应的物方坐标的实现方式如下，设火星的椭球方程为：设火星的椭球方程为：其中，a、b分别表示初始的火星椭球体的赤道半径和极半径；h表示椭球高；a、b分别表示迭代过程中火星椭球体的新的赤道半径和新的极半径；将式(5)代入式(6)得到下列方程式：执行如下步骤，(1)计算视向量与火星椭球面的初始交点坐标令h＝0，解算式(8)，取比例系数μ＝min{μ1,μ2}并代入式(5)，计算交点坐标m
k
(x
k
,y
k
,z
k
)，下标k表示迭代次数，第一次执行步骤(1)时k＝0；μ1、μ2分别表示式(8)的两个根；(2)将第k次计算的用直角坐标表示的交点坐标m
k
(x
k
,y
k
,z
k
)转换为经、纬度和椭球高表示的交点坐标根据计算的经、纬度坐标对火星全球dem进行双线性内插，得到新的椭球高h
k
；(3)重新计算视向量与火星椭球面的交点坐标，令h＝h
k
，解算式(8)，取μ
min
＝min{μ1,μ2}并代入式(5)，计算新的交点坐标m
k+1
(x
k+1
,y
k+1
,z
k+1
)，令k＝k+1；
(4)根据影像分辨率相应的预设阈值ρ，重复(2)、(3)，直至‖m
k
m
k+1
‖＜ρ，此时交点收敛；视向量与全球dem交点的直角坐标m(x,y,z)为仿真影像平面上像素对应的物方坐标，转换为经、纬度和椭球高表示的大地坐标作为步骤4最终的结果。5.根据权利要求1所述火星轨道相机影像仿真方法，其特征在于：步骤5中，根据步骤4计算的仿真影像平面上像素对应的物方点的经、纬度坐标对火星全球影像mdim进行双线性内插，计算对应像素的灰度值。6.一种火星轨道相机影像仿真系统，其特征在于：用于实现如权利要求1-5任一项所述的一种火星轨道相机影像仿真方法。7.根据权利要求6所述火星轨道相机影像仿真系统，其特征在于：包括以下模块，第一模块，用于输入仿真数据源，包括火星全球dem和火星全球数字影像模型mdim；第二模块，用于仿真火星轨道相机的位置和姿态参数，包括利用三次多项式建立火星轨道器姿轨模型，并顾及相机安装矩阵，计算曝光时刻相机的位置和姿态参数；第三模块，用于建立火星轨道相机成像模型，包括根据第二模块计算的曝光时刻相机的位置和姿态参数，建立火星轨道相机成像模型，进一步推导出相机视向量模型；第四模块，用于计算仿真影像平面上像素对应的物方坐标，包括利用第三模块的相机成像模型，将仿真影像平面上像素逐一投影到火星全球dem表面，通过相机视向量模型计算视向量与火星全球dem表面的交点坐标，获得对应物方点的位置；第五模块，用于根据第四模块计算的物方点位置，在火星全球影像mdim上内插计算像素的灰度值；第六模块，用于将第五模块所得火星全球影像mdim上对应像素的灰度值逐一赋值给仿真影像平面上像素。8.根据权利要求6所述火星轨道相机影像仿真系统，其特征在于：包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如权利要求1-5任一项所述的一种火星轨道相机影像仿真方法。9.根据权利要求6所述火星轨道相机影像仿真系统，其特征在于：包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如权利要求1-5任一项所述的一种火星轨道相机影像仿真方法。

技术总结
本发明提供一种火星轨道相机影像仿真方法及系统，输入仿真数据源，包括火星全球DEM和火星全球数字影像模型MDIM；仿真火星轨道相机的位置和姿态参数，包括利用三次多项式建立火星轨道器姿轨模型，并顾及相机安装矩阵，计算曝光时刻相机的位置和姿态参数；建立火星轨道相机成像模型和相机视向量模型；将仿真影像平面上像素逐一投影到火星全球DEM表面；在火星全球影像MDIM上内插计算像素的灰度值；将火星全球影像MDIM上对应像素的灰度值逐一赋值给仿真影像平面上像素。本发明能够真实、准确的反应火星表面的形貌特征，适用于火星轨道的框幅式相机和线阵推扫式相机等所有类型的成像相机影像的仿真，可以针对不同火星形貌特征的区域进行相机影像仿真。区域进行相机影像仿真。区域进行相机影像仿真。


技术研发人员：赵双明 喻国荣
受保护的技术使用者：武汉大学
技术研发日：2021.02.26
技术公布日：2022/3/8