[发明专利]一种地球静止轨道光学相机热变形误差修正方法

权利要求书

1.一种地球静止轨道光学相机热变形误差修正方法，其特征在于包括以下步骤:

步骤1：通过严格几何定位模型，获取空间热变形引起的地球静止轨道光学相机在赤经赤纬方向的定位误差在连续周期内的观测值，分析该定位误差观测值随在轨时间的变化规律；

步骤2：根据最小二乘原理和信赖域收敛算法，采用傅里叶级数模型对地球静止轨道光学相机热变形定位误差观测值随时间的变化趋势进行高精度拟合；设Data1＝{(t_1i,error_1i)|i＝1,2…N₁}为第一天的误差数据集，error_1i为t_1i时刻对应的定位误差值，根据三阶傅里叶级数模型：

其中，a₀，a₁，b₁，a₂，b₂，a₃，b₃，ω为三阶傅里叶级数模型的参数，Time为时间，FError为拟合误差值；根据最小二乘原理：

带入第一天的数据Data1＝{(t_1i,error_1i)|i＝1,2…N}即可求得拟合模型的参数；

步骤3：根据地球静止轨道光学相机热变形定位误差在连续周期内的观测数据，通过Levenberg-Marquardt算法，采用第n+1天的定位误差数据解算由第1，2，...，n天误差拟合模型构成的修正模型的权系数，并保存该权系数矩阵，n为天数；

根据步骤2方法可求得第一、二、三天的误差拟合模型分别为FError1、FError2、FError3，采用第四天的定位误差数据：

Data4＝{(t_4i,error_4i)|i＝1,2…N₄} (3)

解算由第一、二、三天误差拟合模型构成的修正模型如下：

F＝q₀+q₁FError1+q₂FError2+q₃FError3 (4)

式中，q₀，q₁，q₂，q₃为修正模型的权系数，F为修正模型，解算并保存该权系数矩阵；

步骤4：将步骤3获得的权系数矩阵应用到第2，3，...，n+1天的误差拟合模型上，得到第n+2天的误差修正模型，通过该模型对第n+2天的定位误差进行修正；

步骤5：由第n+2天的地球静止轨道光学相机热变形误差观测值减去修正模型的预测值，得到修正后的定位误差；

步骤6：重复执行步骤2、步骤3、步骤4、步骤5，并以卫星周期为单位随时间逐步移动，从而修正每一天的热变形定位误差。