协同多载荷数据的高分辨率遥感卫星颤振探测与建模的方法

技术领域

本发明涉及高分辨率遥感卫星平台颤振探测与建模技术领域，尤其涉及 一种协同多载荷数据的高分辨率遥感卫星颤振探测与建模的方法。

背景技术

颤振是在轨卫星普遍存在的物理现象，是对地观测卫星的主要研究内容 之一。随着高分辨率遥感卫星的质量不断增高，地面分辨率等参数得到大幅 优化，特定频率的颤振逐渐成为影像及其产品质量的重要影响因素。为了降 低颤振的影响，卫星颤振探测与建模的研究逐渐成为热点。

国外在颤振探测与建模领域的研究较早，具有一定的技术积累。国内起 步稍晚，但是由于国产遥感卫星在轨数量增速较快，相应的颤振探测与建模 技术发展迅速。分析不同载荷在平台颤振期间采集的数据来反映平台颤振， 是颤振探测的基本策略。卫星颤振探测与建模的重点和难点在于探测，建模 通常采用一般多项式和三角函数，因此这里主要分析颤振探测的研究进展。 根据所依赖的载荷及其数据，可将颤振探测方法归纳为四类：①基于高频传 感器的颤振探测，如角度传感器和角位移传感器等；②基于线阵影像匹配的颤振探测，如全色影像对和多光谱分时成像影像等；③基于影像定位与地面 控制对比分析的颤振探测，如控制点和控制影像产品等；④基于卫星其他载 荷数据的颤振探测，如星图和激光测高数据等。基于上述方法，不同卫星颤 振探测结果情况如表1所示：

表1

振幅a参考公开文献，结合卫星参数统一转换到地面单位(米)，没有 公开或者无法转换的用“-”表示；方法b表示上述总结的四种方法。

根据各载荷的性质，结合表1可以得出四种卫星颤振探测方法的特：

方法①探测的频率范围较宽，特别针对高频颤振，具有很好的识别效果， 精度较高，但是受到卫星平台空间和卫星系统预算等因素限制，较少的卫星 搭载了高频传感器。方法②是目前高分辨率遥感卫星颤振探测与建模的主流 方法，对低频颤振敏感度较高，也可用于高频颤振探测，但是只有特定的相 机安装关系才适用该方法，且受影像匹配精度制约。方法③主要针对低频颤 振，具有较好的可靠性，但是成本较高，需要大量的地面控制数据。方法④ 灵活性高，根据载荷的特性，设计颤振探测方案，但是精度不稳定，且针对 性较强。

虽然现有颤振探测方法较多，但是均存在局限性，单一方法难以完成高精度 的颤振探测与建模。因此，本领域需要一种能够准确可靠探测卫星平台颤振并建 模的方法。

发明内容

为解决上述技术问题，本发明的目的是提供一种协同多载荷数据的高分 辨率遥感卫星颤振探测与建模的方法，该方法基于多载荷数据探测颤振与融 合建模的研究，利用多载荷数据探测颤振，融合得到可靠性更高、精确性更 好和频域覆盖更全面的颤振结果，实现多种方法优势互补。在此基础上，构 建高精度颤振模型，为卫星颤振辅助姿态确定奠定基础。

本发明的目的通过以下的技术方案来实现：

协同多载荷数据的高分辨率遥感卫星颤振探测与建模的方法，包括：

步骤1基于多载荷数据进行卫星平台颤振探测；

步骤2采用自适应贝叶斯算法融合卫星平台颤振结果；

步骤3顾及惩罚因子对融合后的颤振结果进行建模。

与现有技术相比，本发明的一个或多个实施例可以具有如下优点：

基于高分辨率遥感卫星平台搭载的多种载荷，从多个角度探测卫星平台 颤振，规避了单一算法的局限性，充分利用了多载荷数据的宽频域探测特性； 顾及了长时序姿态频率的变化，对低频误差具有很好的敏感性。

附图说明

图1是协同多载荷数据的高分辨率遥感卫星颤振探测与建模的方法的流程 图；

图2是资源三号卫星多光谱相机多片CCD安装关系示意图；

图3a和图3b是资源三号卫星陀螺和星敏与平台的安装关系示意图；

图4是能量中心计算原理示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合实施例及附图 对本发明作进一步详细的描述。

如图1所示，为协同多载荷数据的高分辨率遥感卫星颤振探测与建模的方法 的流程，包括以下步骤：

步骤1基于多载荷数据开展卫星平台颤振探测；

步骤1.1采用高精度影像匹配方法，对存在不同时成像的影像进行匹配，获 得短时间内的片间视差，从而基于光学相机探测卫星平台颤振；

资源三号卫星多光谱相机包含蓝、绿、红和近红外4个波段，依次记为B1、 B2、B3和B4，如图2所示，B1表示蓝光谱段及三片CCD，B2表示绿光谱段及 三片CCD，B3表示红光谱段及三片CCD，B4表示近红外谱段及三片CCD，多 光谱各波段共视场成像，即同一时刻各波段CCD拍摄区域对应地面相隔一定距离 的地面区域，同一区域则被各波段CCD在不同的连续时刻拍摄到。每个波段由3 片CCD组成，每片CCD包含3072个探测元器件，为了便于后期的影像拼接， 波段内相邻CCD之间具有195个像素的重叠，每个像素大小为0.02mm。在相机 焦平面上，4个波段线阵CCD器件在沿轨方向依次平行摆放。

采用高精度亚像素相位相关匹配方法，对资源三号卫星多光谱分片影像进行匹配，得到片间视差，通过波形叠加理论得到卫星平台的颤振，波形分解公式如下；

g(t

式中，g(t

步骤1.2采用波形分析算法，对陀螺、星敏等测量载荷获取的姿态数据进行分 析，从而基于陀螺、星敏等姿态测量载荷探测卫星平台颤振；

资源三号卫星安装了4组陀螺(G1G3G5、G2G4G6、G7G8G9、G10G11G12) 和3个星敏(ST1、ST2、ST3)，陀螺和星敏与平台的安装关系如图3a和3b所 示。

采用傅里叶变换算法，对资源三号卫星陀螺和星敏数据进行处理，得到卫星平 台颤振信息，傅里叶变换公式如下：

式中，F(u)表示频域值，f(t)表示时域陀螺、星敏数据，t表示时间，u表示 频率，j表示复数的虚部；

步骤1.3采用质心提取算法，对恒星相机获取的星图数据进行分析，结合波 形分析算法，从而基于恒星相机探测卫星平台颤振；

资源三号03星的姿态测量系统包含了3台恒星相机，恒星相机成像得到连续 的星图，通过质心提取算法得到恒星在星图上变化，从而基于恒星相机载荷探测卫 星平台颤振，考虑到星图质心的能量分布情况，采用能量中心计算方法；能量中心 计算原理示意图如图4所示，峰值点B(x

根据图4可知，存在两个几何关系，利用公式表示如下：

通过对公式(3)进行推导化简，可得一元二次方程：

2(y

上式可以简化表示为：

ax

式中，a＝2(y

则公式(5)根据一元二次方程解可得：

式中，满足x

同理，当峰值点B(x

2(y

利用公式(5)(6)可求得公式(7)的x值，其中满足x

步骤1.4采用波形分析算法，对角位移传感器等微振动测量载荷获取的数据 进行分析，从而基于角位移传感器等微振动测量载荷探测卫星平台颤振；

资源三号卫星未装载角位移传感器等高频角动量测量传感器，如在实际作业过程中，目标卫星安装了类似载荷，可采用步骤1.2方法进行处理；

步骤1.5采用地形匹配分析算法，对激光测高仪获取的高程数据进行分析， 从而基于激光测高仪探测卫星平台颤振；

资源三号卫星(02星和03星)装载了激光测高仪，地表高程测量精度分别达 到1.0米和0.3米，与已有高精度地形数据匹配分析，可以通过匹配结果分析卫星 平台的颤振情况，地形匹配的核心是激光严格几何测高模型，公式如下：

式中，

步骤2采用自适应贝叶斯算法融合卫星平台颤振结果；

步骤2.1顾及载荷获取数据的精度和可靠度等，设定自适应的先验概率；

综合考虑所依赖载荷的测量精度、测量可靠度设定先验概率，根据颤振探测的 结果对设定的概率进行调整，以资源三号卫星为例，根据经验设定光学相机的先验 概率为0.3，陀螺星敏0.5，恒星相机0.1，激光测高仪0.1；

步骤2.2采用贝叶斯算法融合多载荷颤振数据，以得到可靠性更高的颤振结 果。

在基于多种载荷探测卫星平台颤振的基础上，根据已设定的先验概率，结合贝 叶斯算法融合多载荷颤振数据，贝叶斯的判定准则如下：

式中，P(O

步骤3顾及惩罚因子对融合后的颤振结果进行建模。

步骤3.1采用一般多项式和三角函数作为颤振预估模型，设定惩罚因子，在区 段时间内建立颤振模型；

针对资源三号卫星颤振建模，采用一般多项式和三角函数拟合，通过引入惩罚 函数估计模型参数，设定惩罚因子，获得最优化目标函数参数。待获取参数的卫星 颤振模型为目标函数，如下所示：

式中，Att(t)为构建的连续颤振模型，f

由于颤振测量受到测量载荷的采样频率限制，导致颤振数据是离散数据，定义 为ATT，为了最优化拟合模型参数，惩罚因子定义为：

式中，

步骤3.2考虑长时序模型的参数变化，优化颤振模型。

为了构建高精度颤振模型，考虑低频部分，因此需要分析卫星颤振模型参数变 化情况。依赖长时间序列的模型参数，构建参数变化模型。

轨道周期、日周期和年周期变化通常为周期性变化，颤振模型参数变化采用三 角函数来表示：

F

式中，i为颤振模型包含的n个参数，F为参数变化函数，a为振幅，ω为角频， θ为初相，l为常数项。

虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本 发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内 的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的 形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所 附的权利要求书所界定的范围为准。