一种光学遥感卫星影像时变系统误差建模补偿方法及系统的制作方法
【技术领域】
[0001] 本发明属于遥感卫星地面预处理领域，特别是涉及到一种光学遥感卫星影像时变 系统误差建模补偿方法及系统。
【背景技术】
[0002] 高分辨率光学卫星影像无控制几何处理受到多种误差源影响，其中由于星敏感器 低频误差、空间热环境导致的安装结构变化以及不同星敏感器工作模式切换所导致的姿态 基准的变化是当前影响光学遥感影像高精度几何处理精度的重要因素。星敏感器是一种常 用的精密姿态测量部件，但是由于受到卫星上的冷热交变的空间热环境以及视场变化等因 素影响，其会产生依轨道周期变化的周期性低频误差。安装结构变化以及星敏感器工作模 式切换导致姿态基准发生变化，进而无法保证不同姿态敏感器间姿态基准高精度耦合。因 此，需要解决如何消除或削弱星敏感器低频误差、姿态基准变化对无控制几何定位的影响， 以满足光学遥感影像地面处理要求。
【发明内容】
[0003] 本发明针对星敏感器低频误差参数难以辨识与补偿以及姿态基准时变性变化的 问题，提供了一种光学遥感卫星影像时变系统误差建模补偿技术方案。
[0004] 本发明提供的技术方案为一种光学遥感卫星影像时变系统误差建模补偿方法，包 括以下步骤：
[0005] 步骤1，基于光学遥感卫星对地相机多个成像时间段的多星敏感器观测数据，解算 星敏感器间相对安装参数变化序列，进一步基于加权平均得到安装参数最优估计值；
[0006] 步骤2,根据多个星敏感器成像时间段的观测数据以及步骤1得到的安装参数，实 现多星敏感器最优信息融合，获取高精度姿态数据；
[0007] 步骤3，根据跟踪侦照的光学定标场全色影像数据，定标场D0M/DEM参考数据，采用 严密几何成像模型，实现对地相机精密姿态反演；
[0008] 步骤4,采用傅立叶级数构建时变系统误差补偿模型，进一步根据步骤2和步骤3所 得结果，基于最小二乘原理实现时变系统误差补偿模型参数最优估计。
[0009] 而且，步骤1中，解算星敏感器间相对安装参数变化序列时，设某时刻t的星敏感器 A的四元数观测值为# ，星敏感器B的四元数观测值为 彳A- 知]T，得到时刻t的本体到惯性系的旋转矩阵圮、每，计算A星敏感 器到B星敏感器的旋转矩阵巧=(圮广圯以及欧拉角转换参数a#: 。
[0010] 而且，步骤2中，获取高精度姿态数据实现方式如下，
[0011 ]设有η个星敏感器，多个星敏感器光轴在惯性系下的矢量坐标为Vms，Vms，…… VnCIS，在本体坐标系矢量坐标为VlBody，V2Body，…，VnBody，基于星敏感器光轴矢量的观测方程 如下，Kt? + V3xl = ,? _ Κ,ν,,,Λ / =】.2,…."
[0012] 其中，忠表示本体坐标系到惯性坐标系的旋转矩阵，乃幻为星敏感器测量噪声；
[0013] 当光轴矢量个数大于等于2时，基于最小二乘原理实现姿态参数的最优估计。
[0014] 而且，步骤3中，所述严密几何成像模型为基于探元指向角模型的严密几何成像模 型。
[0015] 而且，步骤4中，采用傅立叶级数构建时变系统误差补偿模型的实现方式如下，
[0016] 设时变系统误差欧拉角表不为~ Δ# .，则Δ k的三个分量米用傅 立叶函数形式进行建模，得到时变系统误差补偿模型如下，
[0018]
，T表示卫星轨道周期，k表示时间步长个数，τ表示时间步长， Μ为正常数，31^々1^，<^，\_/_」々0」表不未知的时变系统误差系数。
[0019] 一种光学遥感卫星影像时变系统误差建模补偿系统，包括以下模块：
[0020] 第一模块，用于基于光学遥感卫星对地相机多个成像时间段的多星敏感器观测数 据，解算星敏感器间相对安装参数变化序列，进一步基于加权平均得到安装参数最优估计 值；
[0021] 第二模块，用于根据多个星敏感器成像时间段的观测数据以及第一模块得到的安 装参数，实现多星敏感器最优信息融合，获取高精度姿态数据；
[0022] 第三模块，用于根据跟踪侦照的光学定标场全色影像数据，定标场D0M/DEM参考数 据，采用严密几何成像模型，实现对地相机精密姿态反演；
[0023] 第四模块，用于采用傅立叶级数构建时变系统误差补偿模型，进一步根据第二模 块和第三模块所得结果，基于最小二乘原理实现时变系统误差补偿模型参数最优估计。
[0024] 而且，第一模块中，解算星敏感器间相对安装参数变化序列时，设某时刻t的星敏 感器A的四元数观测值为# = ，星敏感器B的四元数观测值为 得到时刻t的本体到惯性系的旋转矩阵圮、忠，计算A星敏感 器到B星敏感器的旋转矩阵=?)7 <以及欧拉角转换参数6? 心/。
[0025]而且，第二模块中，获取高精度姿态数据实现方式如下，
[0026]设有η个星敏感器，多个星敏感器光轴在惯性系下的矢量坐标为Vicis，V2cis，...... VnCIS，在本体坐标系矢量坐标为VlBody，V2Body，…，VnBody，基于星敏感器光轴矢量的观测方程 如下，
[0027] VicIS + ν3χ1 = Rg - Vmdy i = 1, % -,,n.
[0028] 其中，患表示本体坐标系到惯性坐标系的旋转矩阵，V3XA星敏感器测量噪声；
[0029]当光轴矢量个数大于等于2时，基于最小二乘原理实现姿态参数的最优估计。
[0030] 而且，第三模块中，所述严密几何成像模型为基于探元指向角模型的严密几何成 像模型。
[0031] 而且，第四模块中，采用傅立叶级数构建时变系统误差补偿模型的实现方式如下，
[0032] 设时变系统误差欧拉角表示为的三个分量采用傅 立叶函数形式进行建模，得到时变系统误差补偿模型如下，
[0034]
.，Τ表示卫星轨道周期，k表示时间步长个数，τ表示时间步长， Μ为正常数，at j，bt j，《 _，ae j，b θ j表不未知的时变系统误差系数。
[0035]本发明提供了一种光学遥感卫星影像时变系统误差建模补偿技术方案，实现了星 敏感器低频误差、姿态基准不统一误差的在轨补偿，通过本发明提供的技术方案可以有效 削弱无控制定位时变系统误差影响，为光学遥感影像几何处理奠定基础。
【附图说明】
[0036]图1为本发明实施例的流程图。
【具体实施方式】
[0037] 以下结合附图和实施例详细说明本发明技术方案。
[0038] 参见图1所示光学遥感卫星影像无控制几何定位的时变系统误差在轨标定与补偿 流程图，以下针对实施例流程中的各步骤，对本发明方法做进一步详细描述。
[0039]步骤1，基于光学遥感卫星对地相机多个成像时间段的多星敏感器观测数据，解算 星敏感器间相对安装参数变化序列，进一步基于加权平均得到安装参数最优估计值，实现 星敏感器安装参数相对标定。
[0040]假设某时刻t的星敏感器A的四元数观测值为# 星敏感 器B的四元数观测值为gli; ,得到时亥Ijt的本体到惯性系的旋转矩阵 圮、埒表达式如下：
[0041]
[0043] 进一步得到A星敏感器到B星敏感器的旋转矩阵^^ = (<)7 <以及欧拉角转换参数 ，按此对多个时刻分别进行处理，可以得到解算星敏感器间相对安装 参数变化序列。基于上述方法对多个成像段观测数据进行处理得到每个成像时间段相对安 装参数序列。
[0044] 对每个成像时间段相对安装参数序列进行误差统计分析得到每个成像时间段安 装参数中误差，最后基于每个成像时间段欧拉角转换参数进行整体加权平均得到最优估计 值。
[0045] 步骤2,根据多个成像时间段的星敏感器原始观测数据以及步骤1得到的安装参数 最优估计值，实现多星敏感器最优信息融合，获取高精度姿态数据。
[0046] 光学遥感卫星在轨成像过程中，由于空间环境复杂性以及星敏感器自身工作特 点，经常出现工作模式切换，导致单一星敏感器定姿、双星敏感器定姿以及多星敏感器定姿 等情况产生，故可以根据步骤1标定的星敏感器间的安装关系（即安装参数最优估计值）以 及星敏感器原始观测数据（即步骤1所述四元数观测值），实现最优光轴矢量观测信息融合， 得到姿态基准高度统一的姿态数据。
[0047] 假设有η个星敏感器，多个星敏感器光轴在惯性系下的矢量坐标为V1CIS，V2CIS，…… VnCIS，在本体坐标系矢量坐标为VlBody，V2Body，…，VnBody，基于星敏感器光轴矢量的观测方程 如下：
[0049] 起表示本体坐标系到惯性坐标系的旋转矩阵，v3X1为星敏感器测量噪声。当光轴 矢量个数大于等于2时，基于最小二乘原理实现姿态参数的最优估计。
[0050] 步骤3,根据长年跟踪侦照的光学定标场全色影像数据，定标场D0M/DEM参考数据， 进行控制点自动量测，一般包括特征点提取、影像模拟、金字塔影像匹配、整像素匹配、子像 素匹配以及粗差点剔除，具体实现可采用现有技术;进一步采用严密几何成像模型，实现对 地相机精密姿态反演解算。
[0051] 即对地相机精密姿态反演数学模型为基于探元指向角模型的严密几何成像模型， 具体实现形式如下：
[0053] 其中，（Xg，Yg，Zg)表示地物点的物方坐标；（Ms)，ilv( S))表示探元号s的指向角大 小；（xgps，Ygps，zgps)与(Bx，By，Bz)分别表示对地相机摄影中心的物方坐标与GPS偏心误差;λ 表示比例系数；、Cy分别表示由WGS84坐标系到J2000坐标系旋转矩阵、由 J2000坐标系到卫星本体坐标系旋转矩阵以及由卫星本体坐标系到相机测量坐标系的旋转 矩阵。
[0054] 由上式可以得到，对于线阵推扫相机，当每个扫描行匹配得到的非共线控制点观 测矢量个数大于等