CN105526951A - 一种星敏感器原始观测数据预处理方法及系统 - Google Patents
一种星敏感器原始观测数据预处理方法及系统 Download PDFInfo
- Publication number
- CN105526951A CN105526951A CN201610066655.8A CN201610066655A CN105526951A CN 105526951 A CN105526951 A CN 105526951A CN 201610066655 A CN201610066655 A CN 201610066655A CN 105526951 A CN105526951 A CN 105526951A
- Authority
- CN
- China
- Prior art keywords
- star sensor
- star
- follows
- optical axis
- delta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
- G01C21/165—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
Abstract
本发明提供了一种星敏感器原始观测数据预处理方法及系统,包括基于陀螺时标信息对多星敏感器原始观测数据的时间基准进行统一;计算多星敏感器间光轴夹角时间变化序列;计算星敏感器理论光轴夹角大小,进一步构建星敏感器原始观测数据质量控制判定模型,探测星敏感器中可能存在的粗差,最后将每个星敏感器的时间序列观测值分成正常输出观测值组与非正常输出观测值组;将正常输出观测值作为拟合基准点,基于滑动窗口多项式拟合对非正常输出的星敏感器观测值进行修复。本发明可以实现星敏感器不正常输出观测值的探测与修复,有效削弱星敏感器粗差观测值影响,进一步实现高精度姿态确定,为高分辨率光学影像高精度几何处理提供保证。
Description
技术领域
本发明属于遥感卫星地面预处理领域,特别是涉及到一种星敏感器原始观测数据预处理方法及系统。
背景技术
目前高分辨率光学卫星广泛采用由星敏感器和陀螺构成的卫星姿态确定系统实现精密姿态确定,并用于光学影像的高精度几何处理。星敏感器作为一种常用的高精度姿态测量部件,其主要工作原理是通过恒星相机对恒星进行摄影,利用获取的恒星影像确定卫星的姿态,处理流程包括星点提取、导航星表建立与优化以及星图识别等关键步骤。由于卫星在轨运行过程中所处的空间环境比较复杂以及成像器件本身特性,导致以下干扰因素的产生:(1)星敏感器获取的星图存在各种噪声,例如光子散粒噪声、读出噪声以及暗电流噪声;(2)星图数据纹理信息匮乏、信噪比比较低,存在拖尾现象;(3)星图中导航星分布不均匀;(4)星图匹配识别错误。基于以上因素会导致星敏感器输出的观测值中存在一些不正常观测值,因此需要对星敏感器输出的观测值进行可用性判识,以免将非正常工作中的观测数据引入组合定姿系统中,影响最终精密姿态的确定与光学影像的几何处理。
发明内容
本发明针对星敏感器输出观测值中存在粗差难以探测与修复的问题,提供了一种基于光轴夹角变化检测技术的星敏感器观测数据预处理技术方案。
本发明提供的技术方案为一种星敏感器原始观测数据预处理方法,包括以下步骤:
步骤1,基于陀螺时标信息对多星敏感器原始观测数据的时间基准进行统一;
步骤2,基于步骤1的结果,计算多星敏感器间光轴夹角时间变化序列,多星敏感器间光轴夹角计算实现方式如下,
设在相同时刻t的星敏感器A的四元数观测值为星敏感器B的四元数观测值为得到时刻t的本体到惯性系的星敏感器A、B分别相应旋转矩阵表达式如下,
进一步得到星敏感器A、星敏感器B的光轴矢量ZA,ZB在惯性系中表示如下,
ZA=[2(q1Aq3A+q2Aq0A)2(q2Aq3A-q1Aq0A)-q1A 2-q2A 2+q3A 2+q0A 2]T
ZB=[2(q1Bq3B+q2Bq0B)2(q2Bq3B-q1Bq0B)-q1B 2-q2B 2+q3B 2+q0B 2]T
计算t时刻两光轴的夹角αt,
αt=arccos(ZA·ZB)
步骤3,基于星敏感器地面标定的安装矩阵参数计算星敏感器理论光轴夹角大小,进一步构建星敏感器原始观测数据质量控制判定模型,探测星敏感器中存在的粗差,最后将每个星敏感器的时间序列观测值分成正常输出观测值组与非正常输出观测值组;
所述星敏感器理论光轴夹角计算实现方式如下,
设星敏感器A与星敏感器B的安装矩阵分别表示如下,
则星敏感器A与星敏感器B的光轴在卫星本体坐标系中矢量表示如下,
得到星敏感器A与B的光轴在本体系中的夹角αAB如下,
所述星敏感器原始观测数据的质量控制判定模型如下,
其中,m表示阈值系数,αi表示观测时刻ti的星敏感器A与B的光轴夹角大小,δα为星敏感器A与星敏感器B光轴夹角中误差;
步骤4,根据步骤3的星敏感器所得观测值分组结果,将正常输出观测值作为拟合基准点,基于滑动窗口多项式拟合对非正常输出的星敏感器观测值进行修复。
而且,步骤1中,对多星敏感器原始观测数据的时间基准进行统一采用球面线性插值模型。
而且,步骤3中,误差δα计算如下,
其中,N表示星敏感器观测历元个数。
而且,步骤4中,基于滑动窗口多项式拟合对非正常输出的星敏感器观测值进行修复,实现方式如下,
设某星敏感器A输出的正常输出观测值组包括n个时间序列观测值(q1,q2,q3…qn-1,qn),tk为非正常输出观测值组中的起始时刻,
记n个历元的姿态四元数为(q0i,q1i,q2i,q3i)i=1,2,…,n,相应m-1次最佳正交多项式拟合如下,
其中,t表示时刻,aj,j=0,1,…,m-1表示多项式系数,设上式为各正交多项式δj(t)的线性组合,
其中,cj,j=0,1,…,m-1表示正交多项式系数;
根据最小二乘原理得tk时刻姿态四元数拟合值如下,
其中,表示tk时刻四元数矢量部分拟合值,表示tk时刻四元数标量部分拟合值,表示四元数矢量部分正交多项式拟合系数, 表示四元数矢量部分正交多项式。
本发明相应提供一种星敏感器原始观测数据预处理系统,包括以下模块:
第一模块,用于基于陀螺时标信息对多星敏感器原始观测数据的时间基准进行统一;
第二模块,用于计算多星敏感器间光轴夹角时间变化序列,多星敏感器间光轴夹角计算实现方式如下,
设在相同时刻t的星敏感器A的四元数观测值为星敏感器B的四元数观测值为得到时刻t的本体到惯性系的星敏感器A、B分别相应旋转矩阵表达式如下,
进一步得到星敏感器A、星敏感器B的光轴矢量ZA,ZB在惯性系中表示如下,
ZA=[2(q1Aq3A+q2Aq0A)2(q2Aq3A-q1Aq0A)-q1A 2-q2A 2+q3A 2+q0A 2]T
ZB=[2(q1Bq3B+q2Bq0B)2(q2Bq3B-q1Bq0B)-q1B 2-q2B 2+q3B 2+q0B 2]T
计算t时刻两光轴的夹角αt,
αt=arccos(ZA·ZB)
第三模块,用于基于星敏感器地面标定的安装矩阵参数计算星敏感器理论光轴夹角大小,进一步构建星敏感器原始观测数据质量控制判定模型,探测星敏感器中存在的粗差,最后将每个星敏感器的时间序列观测值分成正常输出观测值组与非正常输出观测值组;
所述星敏感器理论光轴夹角计算实现方式如下,
设星敏感器A与星敏感器B的安装矩阵分别表示如下,
则星敏感器A与星敏感器B的光轴在卫星本体坐标系中矢量表示如下,
得到星敏感器A与B的光轴在本体系中的夹角αAB如下,
所述星敏感器原始观测数据的质量控制判定模型如下,
其中,m表示阈值系数,αi表示观测时刻ti的星敏感器A与B的光轴夹角大小,δα为星敏感器A与星敏感器B光轴夹角中误差;
第四模块,用于根据第三模块的星敏感器所得观测值分组结果,将正常输出观测值作为拟合基准点,基于滑动窗口多项式拟合对非正常输出的星敏感器观测值进行修复。
而且,第一模块中,对多星敏感器原始观测数据的时间基准进行统一采用球面线性插值模型。
而且,第三模块中,误差δα计算如下,
其中,N表示星敏感器观测历元个数。
而且,第四模块中,基于滑动窗口多项式拟合对非正常输出的星敏感器观测值进行修复,实现方式如下,
设某星敏感器A输出的正常输出观测值组包括n个时间序列观测值(q1,q2,q3…qn-1,qn),tk为非正常输出观测值组中的起始时刻,
记n个历元的姿态四元数为(q0i,q1i,q2i,q3i)i=1,2,…,n,相应m-1次最佳正交多项式拟合如下,
其中,t表示时刻,aj,j=0,1,…,m-1表示多项式系数,设上式为各正交多项式δj(t)的线性组合,
其中,cj,j=0,1,…,m-1表示正交多项式系数;
根据最小二乘原理得tk时刻姿态四元数拟合值如下,
其中,表示tk时刻四元数矢量部分拟合值,表示tk时刻四元数标量部分拟合值,表示四元数矢量部分正交多项式拟合系数, 表示四元数矢量部分正交多项式。
本发明提供了一种基于光轴夹角变化检测技术的星敏感器原始观测数据预处理技术方案,实现了星敏感器非正常输出观测值的有效辨识与修复,通过本发明提供的技术方案可以有效削弱星敏感器粗差观测值影响,进一步实现高精度姿态确定,为高分辨率光学影像高精度几何处理提供保证。
附图说明
图1为本发明实施例的流程图。
具体实施方式
以下结合附图和实施例详细说明本发明技术方案。
图1表示的是基于光轴夹角变化检测技术的星敏感器原始观测数据预处理方法实施流程图,以下针对实施例流程中的各步骤,对本发明方法做进一步详细描述。
步骤1,基于陀螺输出观测信息,对多星敏感器观测值的时间基准进行统一,保证多星敏感器输出观测值时刻一致性。
通常包括两个或以上星敏感器,处理方式类似。陀螺输出观测信息即陀螺输出观测数据,多星敏感器观测值即多星敏感器输出观测值。基于这两者实现对多星敏感器观测值的时间基准进行统一,得到星敏感器时间同步数据。陀螺作为星体的短期姿态参考,能够连续提供星体的三轴姿态角速度信息,但由于多台星敏感器与陀螺并不能实现完全时间同步,为了后续在步骤2得到多星敏感器间光轴夹角变化时间序列,本发明基于陀螺的时标信息采用球面线性插值模型实现多星敏感器间原始观测数据时间同步。具体实现方式如下:
球面线性插值不仅会保持其单位长度不变,并且若将q1和q2看为四维空间单位球上的两个点,其将以恒定的速率扫过q1和q2之间的夹角θ。假设t时刻所需插值计算的四元数q(t)位于连接q1和q2的弧上,与q1构成的夹角为tθ(t∈[0,1])与q2构成的夹角为(1-t)θ(t∈[0,1]),于是q(t)可表示为
q(t)=C1(t)q1+C2(t)q2(1)
根据相似三角形关系,可以推导出系数C1(t)和C2(t)的表达式:
即:
步骤2,计算多星敏感器间光轴夹角时间序列数据:基于步骤1的结果,进一步计算多星敏感器间光轴夹角时间变化序列,得到星敏感器间光轴夹角数据。当星敏观测误差满足正态分布,那么时间序列光轴夹角变化也满足正态分布,夹角期望为无偏估计。多星敏感器间光轴夹角数据计算方法的具体实现方式如下:
假设在相同时刻t的星敏感器A的四元数观测值为星敏感器B的四元数观测值为得到时刻t的本体到惯性系的星敏感器A、B分别相应旋转矩阵表达式如下:
进一步得到星敏感器A、星敏感器B的光轴矢量ZA,ZB在惯性系中表示如下:
ZA=[2(q1Aq3A+q2Aq0A)2(q2Aq3A-q1Aq0A)-q1A 2-q2A 2+q3A 2+q0A 2]T(5)
ZB=[2(q1Bq3B+q2Bq0B)2(q2Bq3B-q1Bq0B)-q1B 2-q2B 2+q3B 2+q0B 2]T
可计算t时刻两光轴的夹角αt:
αt=arccos(ZA·ZB)(6)
故基于上述算法,对于时间变化中每个时刻分别计算光轴间的夹角,可以得到多星敏感器间光轴夹角变化时间序列,为步骤3实施的前提条件。
步骤3,星敏感器观测值粗差探测:基于星敏感器地面标定的安装矩阵参数解算星敏感器光轴夹角大小,进一步构建星敏感器原始观测数据质量控制判定模型,探测星敏感器中可能存在的粗差,最终将每个星敏感器的时间序列观测值分成正常输出观测值组与非正常输出观测值组,得到分组后的星敏感器观测值序列。该步骤包括以下子步骤:
步骤3.1,基于地面标定的星敏感器的安装矩阵参数计算星敏感器之间的光轴夹角大小,具体实现方式如下:
设星敏感器A与星敏感器B的安装矩阵分别表示如下:
则星敏感器A与星敏感器B的光轴在卫星本体坐标系中矢量可以表示成:
同样可以得到星敏感器A与B的光轴在本体系中的夹角αAB:
步骤3.2,计算星敏感器A与星敏感器B光轴夹角中误差δα,具体实现方式如下:
其中,αi表示观测时刻ti的星敏感器A与B的光轴夹角大小,通过步骤2得到;N表示星敏感器观测历元个数。
步骤3.3,构建星敏感器观测数据质量控制判定模型,探测星敏感器中可能存在的粗差,最终将每个星敏感器的时间序列观测值分成正常输出观测值组与非正常输出观测值组。具体实现方式如下:
星敏感器原始观测数据质量控制判定模型如下:
其中,m表示阈值系数,且1≤m≤3。故依据上述构建的判定模型对星敏感器输出的观测值时间序列进行判定,将正常与非正常输出的观测值各自构成一组时间序列,即分组后的星敏感器观测值序列。
步骤4,星敏感器观测值数据修复,得到多星敏感器可靠观测值数据:根据步骤3的星敏感器A与B时间序列观测值分组结果,将正常输出观测值作为拟合基准点,基于滑动窗口多项式拟合对非正常输出的星敏感器原始观测值进行修复。具体过程如下:
以星敏感器A为例,设其输出的正常观测值序列组包括n个时间序列观测值(q1,q2,q3…qn-1,qn),非正常输出观测值的时刻序列组根据实际分组情况得到,例如(tk,tk+2,tk+6…,tk+l),tk为非正常输出观测值时刻序列组中的起始时刻,采用的拟合多项式模型如下:
记n个历元的姿态四元数为(q0i,q1i,q2i,q3i)i=1,2,…,n,其中qri,r=1,2,3参数的m-1次最佳正交多项式拟合为:
其中,t表示时刻,aj,j=0,1,…,m-1表示多项式系数。设上式为各正交多项式δj(t)(j=0,1,…,m-1)的线性组合:
其中,cj,j=0,1,…,m-1表示正交多项式系数。
其中δj(t)可用递推公式构造:
δ0(t)=1;δ1(t)=(t-α1);
δj(t)=(t-αj)δj-1(t)-βjδj-2(t)(14)
j=2,…,m-1
根据最小二乘原理可得tk时刻姿态四元数拟合值:
其中,dj表示正交多项式平方和,表示tk时刻四元数矢量部分拟合值,表示tk时刻四元数标量部分拟合值,表示四元数矢量部分正交多项式拟合系数,表示四元数矢量部分正交多项式。故根据上述模型对非正常输出观测值的时刻序列进行修复。
具体实施时,本发明所提供方法可基于软件技术实现自动运行流程,也可采用模块化方式实现相应系统。
本发明相应提供一种星敏感器原始观测数据预处理系统,包括以下模块:
第一模块,用于基于陀螺时标信息对多星敏感器原始观测数据的时间基准进行统一;
第二模块,用于计算多星敏感器间光轴夹角时间变化序列,多星敏感器间光轴夹角计算实现方式如下,
设在相同时刻t的星敏感器A的四元数观测值为星敏感器B的四元数观测值为得到时刻t的本体到惯性系的星敏感器A、B分别相应旋转矩阵表达式如下,
进一步得到星敏感器A、星敏感器B的光轴矢量ZA,ZB在惯性系中表示如下,
ZA=[2(q1Aq3A+q2Aq0A)2(q2Aq3A-q1Aq0A)-q1A 2-q2A 2+q3A 2+q0A 2]T
ZB=[2(q1Bq3B+q2Bq0B)2(q2Bq3B-q1Bq0B)-q1B 2-q2B 2+q3B 2+q0B 2]T
计算t时刻两光轴的夹角αt,
αt=arccos(ZA·ZB)
第三模块,用于基于星敏感器地面标定的安装矩阵参数计算星敏感器理论光轴夹角大小,进一步构建星敏感器原始观测数据质量控制判定模型,探测星敏感器中存在的粗差,最后将每个星敏感器的时间序列观测值分成正常输出观测值组与非正常输出观测值组;
所述星敏感器理论光轴夹角计算实现方式如下,
设星敏感器A与星敏感器B的安装矩阵分别表示如下,
则星敏感器A与星敏感器B的光轴在卫星本体坐标系中矢量表示如下,
得到星敏感器A与B的光轴在本体系中的夹角αAB如下,
所述星敏感器原始观测数据的质量控制判定模型如下,
其中,m表示阈值系数,αi表示观测时刻ti的星敏感器A与B的光轴夹角大小,δα为星敏感器A与星敏感器B光轴夹角中误差;
第四模块,用于根据第三模块的星敏感器所得观测值分组结果,将正常输出观测值作为拟合基准点,基于滑动窗口多项式拟合对非正常输出的星敏感器观测值进行修复。
各模块具体实现可参见相应步骤,本发明不予赘述。
本文中所描述的具体实例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (8)
1.一种星敏感器原始观测数据预处理方法,其特征在于,包括以下步骤:
步骤1,基于陀螺时标信息对多星敏感器原始观测数据的时间基准进行统一;
步骤2,基于步骤1的结果,计算多星敏感器间光轴夹角时间变化序列,多星敏感器间光轴夹角计算实现方式如下,
设在相同时刻t的星敏感器A的四元数观测值为 星敏感器B的四元数观测值为 得到时刻t的本体到惯性系的星敏感器A、B分别相应旋转矩阵表达式如下,
进一步得到星敏感器A、星敏感器B的光轴矢量ZA,ZB在惯性系中表示如下,
ZA=[2(q1Aq3A+q2Aq0A)2(q2Aq3A-q1Aq0A)-q1A 2-q2A 2+q3A 2+q0A 2]T
ZB=[2(q1Bq3B+q2Bq0B)2(q2Bq3B-q1Bq0B)-q1B 2-q2B 2+q3B 2+q0B 2]T
计算t时刻两光轴的夹角αt,
αt=arccos(ZA·ZB)
步骤3,基于星敏感器地面标定的安装矩阵参数计算星敏感器理论光轴夹角大小,进一步构建星敏感器原始观测数据质量控制判定模型,探测星敏感器中存在的粗差,最后将每个星敏感器的时间序列观测值分成正常输出观测值组与非正常输出观测值组;
所述星敏感器理论光轴夹角计算实现方式如下,
设星敏感器A与星敏感器B的安装矩阵分别表示如下,
则星敏感器A与星敏感器B的光轴在卫星本体坐标系中矢量表示如下,
得到星敏感器A与B的光轴在本体系中的夹角αAB如下,
所述星敏感器原始观测数据的质量控制判定模型如下,
其中,m表示阈值系数,αi表示观测时刻ti的星敏感器A与B的光轴夹角大小,δα为星敏感器A与星敏感器B光轴夹角中误差;
步骤4,根据步骤3的星敏感器所得观测值分组结果,将正常输出观测值作为拟合基准点,基于滑动窗口多项式拟合对非正常输出的星敏感器观测值进行修复。
2.根据权利要求1所述星敏感器原始观测数据预处理方法,其特征在于:步骤1中,对多星敏感器原始观测数据的时间基准进行统一采用球面线性插值模型。
3.根据权利要求1所述星敏感器原始观测数据预处理方法,其特征在于:步骤3中,误差δα计算如下,
其中,N表示星敏感器观测历元个数。
4.根据权利要求1或2或3所述星敏感器原始观测数据预处理方法,其特征在于:步骤4中,基于滑动窗口多项式拟合对非正常输出的星敏感器观测值进行修复,实现方式如下,
设某星敏感器A输出的正常输出观测值组包括n个时间序列观测值(q1,q2,q3…qn-1,qn),tk为非正常输出观测值组中的起始时刻,
记n个历元的姿态四元数为(q0i,q1i,q2i,q3i)i=1,2,…,n,相应m-1次最佳正交多项式拟合如下,
其中,t表示时刻,aj,j=0,1,…,m-1表示多项式系数,设上式为各正交多项式δj(t)的线性组合,
其中,cj,j=0,1,…,m-1表示正交多项式系数;
根据最小二乘原理得tk时刻姿态四元数拟合值如下,
其中,表示tk时刻四元数矢量部分拟合值,表示tk时刻四元数标量部分拟合值,表示四元数矢量部分正交多项式拟合系数, 表示四元数矢量部分正交多项式。
5.一种星敏感器原始观测数据预处理系统,其特征在于,包括以下模块:
第一模块,用于基于陀螺时标信息对多星敏感器原始观测数据的时间基准进行统一;
第二模块,用于计算多星敏感器间光轴夹角时间变化序列,多星敏感器间光轴夹角计算实现方式如下,
设在相同时刻t的星敏感器A的四元数观测值为 星敏感器B的四元数观测值为 得到时刻t的本体到惯性系的星敏感器A、B分别相应旋转矩阵表达式如下,
进一步得到星敏感器A、星敏感器B的光轴矢量ZA,ZB在惯性系中表示如下,
ZA=[2(q1Aq3A+q2Aq0A)2(q2Aq3A-q1Aq0A)-q1A 2-q2A 2+q3A 2+q0A 2]T
ZB=[2(q1Bq3B+q2Bq0B)2(q2Bq3B-q1Bq0B)-q1B 2-q2B 2+q3B 2+q0B 2]T
计算t时刻两光轴的夹角αt,
αt=arccos(ZA·ZB)
第三模块,用于基于星敏感器地面标定的安装矩阵参数计算星敏感器理论光轴夹角大小,进一步构建星敏感器原始观测数据质量控制判定模型,探测星敏感器中存在的粗差,最后将每个星敏感器的时间序列观测值分成正常输出观测值组与非正常输出观测值组;
所述星敏感器理论光轴夹角计算实现方式如下,
设星敏感器A与星敏感器B的安装矩阵分别表示如下,
则星敏感器A与星敏感器B的光轴在卫星本体坐标系中矢量表示如下,
得到星敏感器A与B的光轴在本体系中的夹角αAB如下,
所述星敏感器原始观测数据的质量控制判定模型如下,
其中,m表示阈值系数,αi表示观测时刻ti的星敏感器A与B的光轴夹角大小,δα为星敏感器A与星敏感器B光轴夹角中误差;
第四模块,用于根据第三模块的星敏感器所得观测值分组结果,将正常输出观测值作为拟合基准点,基于滑动窗口多项式拟合对非正常输出的星敏感器观测值进行修复。
6.根据权利要求5所述星敏感器原始观测数据预处理系统,其特征在于:第一模块中,对多星敏感器原始观测数据的时间基准进行统一采用球面线性插值模型。
7.根据权利要求5所述星敏感器原始观测数据预处理系统,其特征在于:第三模块中,误差δα计算如下,
其中,N表示星敏感器观测历元个数。
8.根据权利要求5或6或7所述星敏感器原始观测数据预处理系统,其特征在于:第四模块中,基于滑动窗口多项式拟合对非正常输出的星敏感器观测值进行修复,实现方式如下,
设某星敏感器A输出的正常输出观测值组包括n个时间序列观测值(q1,q2,q3…qn-1,qn),tk为非正常输出观测值组中的起始时刻,
记n个历元的姿态四元数为(q0i,q1i,q2i,q3i)i=1,2,…,n,相应m-1次最佳正交多项式拟合如下,
其中,t表示时刻,aj,j=0,1,…,m-1表示多项式系数,设上式为各正交多项式δj(t)的线性组合,
其中,cj,j=0,1,…,m-1表示正交多项式系数;
根据最小二乘原理得tk时刻姿态四元数拟合值如下,
其中,表示tk时刻四元数矢量部分拟合值,表示tk时刻四元数标量部分拟合值,表示四元数矢量部分正交多项式拟合系数, 表示四元数矢量部分正交多项式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610066655.8A CN105526951B (zh) | 2016-01-30 | 2016-01-30 | 一种星敏感器原始观测数据预处理方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610066655.8A CN105526951B (zh) | 2016-01-30 | 2016-01-30 | 一种星敏感器原始观测数据预处理方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105526951A true CN105526951A (zh) | 2016-04-27 |
CN105526951B CN105526951B (zh) | 2018-08-21 |
Family
ID=55769321
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610066655.8A Active CN105526951B (zh) | 2016-01-30 | 2016-01-30 | 一种星敏感器原始观测数据预处理方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105526951B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106052716A (zh) * | 2016-05-25 | 2016-10-26 | 南京航空航天大学 | 惯性系下基于星光信息辅助的陀螺误差在线标定方法 |
CN106289238A (zh) * | 2016-08-02 | 2017-01-04 | 哈尔滨工业大学 | 一种星敏感器数据优化处理方法 |
CN106382928A (zh) * | 2016-08-26 | 2017-02-08 | 北京控制工程研究所 | 一种基于卷帘门曝光星敏感器的动态补偿方法 |
CN106598253A (zh) * | 2016-12-23 | 2017-04-26 | 北京搜狐新媒体信息技术有限公司 | 一种数据预测方法及装置 |
CN109029504A (zh) * | 2018-08-08 | 2018-12-18 | 中国空间技术研究院 | 一种基于光信号的整星级恒星敏感器测试系统及方法 |
CN111174812A (zh) * | 2020-01-19 | 2020-05-19 | 北京空间飞行器总体设计部 | 一种卫星星敏感器常值输出异常自主诊断方法 |
CN111561936A (zh) * | 2020-05-19 | 2020-08-21 | 中国科学院微小卫星创新研究院 | 旋转大幅宽光学卫星精准处理方法及系统 |
CN111637879A (zh) * | 2020-04-22 | 2020-09-08 | 北京控制工程研究所 | 基于多维差异化误差特性分布的双星敏加权姿态确定方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20010045128A1 (en) * | 2000-05-24 | 2001-11-29 | Mccall Hiram | Method of measuring motion |
US20020077768A1 (en) * | 1999-08-18 | 2002-06-20 | Mccall Hiram | Processing method for motion measurement |
CN101893440A (zh) * | 2010-05-19 | 2010-11-24 | 哈尔滨工业大学 | 基于星敏感器的天文自主导航方法 |
US8019209B2 (en) * | 2005-09-30 | 2011-09-13 | Hoya Corporation | Optical axis correction apparatus of an imaging device, and optical axis correction method for an imaging device |
CN102735260A (zh) * | 2012-06-18 | 2012-10-17 | 航天东方红卫星有限公司 | 一种星敏感器在轨测量误差的确定方法 |
CN102175260B (zh) * | 2010-12-31 | 2012-11-14 | 北京控制工程研究所 | 一种自主导航系统误差校正方法 |
-
2016
- 2016-01-30 CN CN201610066655.8A patent/CN105526951B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020077768A1 (en) * | 1999-08-18 | 2002-06-20 | Mccall Hiram | Processing method for motion measurement |
US20010045128A1 (en) * | 2000-05-24 | 2001-11-29 | Mccall Hiram | Method of measuring motion |
US8019209B2 (en) * | 2005-09-30 | 2011-09-13 | Hoya Corporation | Optical axis correction apparatus of an imaging device, and optical axis correction method for an imaging device |
CN101893440A (zh) * | 2010-05-19 | 2010-11-24 | 哈尔滨工业大学 | 基于星敏感器的天文自主导航方法 |
CN102175260B (zh) * | 2010-12-31 | 2012-11-14 | 北京控制工程研究所 | 一种自主导航系统误差校正方法 |
CN102735260A (zh) * | 2012-06-18 | 2012-10-17 | 航天东方红卫星有限公司 | 一种星敏感器在轨测量误差的确定方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106052716A (zh) * | 2016-05-25 | 2016-10-26 | 南京航空航天大学 | 惯性系下基于星光信息辅助的陀螺误差在线标定方法 |
CN106052716B (zh) * | 2016-05-25 | 2019-04-05 | 南京航空航天大学 | 惯性系下基于星光信息辅助的陀螺误差在线标定方法 |
CN106289238A (zh) * | 2016-08-02 | 2017-01-04 | 哈尔滨工业大学 | 一种星敏感器数据优化处理方法 |
CN106382928A (zh) * | 2016-08-26 | 2017-02-08 | 北京控制工程研究所 | 一种基于卷帘门曝光星敏感器的动态补偿方法 |
CN106382928B (zh) * | 2016-08-26 | 2019-05-24 | 北京控制工程研究所 | 一种基于卷帘曝光星敏感器的动态补偿方法 |
CN106598253A (zh) * | 2016-12-23 | 2017-04-26 | 北京搜狐新媒体信息技术有限公司 | 一种数据预测方法及装置 |
CN106598253B (zh) * | 2016-12-23 | 2019-12-10 | 北京搜狐新媒体信息技术有限公司 | 一种数据预测方法及装置 |
CN109029504A (zh) * | 2018-08-08 | 2018-12-18 | 中国空间技术研究院 | 一种基于光信号的整星级恒星敏感器测试系统及方法 |
CN111174812A (zh) * | 2020-01-19 | 2020-05-19 | 北京空间飞行器总体设计部 | 一种卫星星敏感器常值输出异常自主诊断方法 |
CN111637879A (zh) * | 2020-04-22 | 2020-09-08 | 北京控制工程研究所 | 基于多维差异化误差特性分布的双星敏加权姿态确定方法 |
CN111637879B (zh) * | 2020-04-22 | 2021-10-01 | 北京控制工程研究所 | 基于多维差异化误差特性分布的双星敏加权姿态确定方法 |
CN111561936A (zh) * | 2020-05-19 | 2020-08-21 | 中国科学院微小卫星创新研究院 | 旋转大幅宽光学卫星精准处理方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105526951B (zh) | 2018-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105526951A (zh) | 一种星敏感器原始观测数据预处理方法及系统 | |
Kolomenkin et al. | Geometric voting algorithm for star trackers | |
CN103323026B (zh) | 星敏感器和有效载荷的姿态基准偏差估计与修正方法 | |
US6876926B2 (en) | Method and system for processing pulse signals within an inertial navigation system | |
CN105698764B (zh) | 一种光学遥感卫星影像时变系统误差建模补偿方法及系统 | |
CN108592945A (zh) | 一种惯性/天文组合系统误差的在线标定方法 | |
CN105548976A (zh) | 船载雷达海上精度鉴定方法 | |
CN104792340A (zh) | 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法 | |
CN103674021A (zh) | 基于捷联惯导与星敏感器的组合导航系统及方法 | |
Rousseau et al. | Star recognition algorithm for APS star tracker: Oriented triangles | |
CN101968361A (zh) | 基于星光观测的空间绝对定向技术 | |
CN111798523A (zh) | 星相机在轨定标定姿及遥感影像几何定位方法、系统 | |
Gou et al. | INS/CNS integrated navigation based on corrected infrared earth measurement | |
CN112880669B (zh) | 一种航天器星光折射和单轴旋转调制惯性组合导航方法 | |
CN104931071A (zh) | 一种基于选权迭代的星敏感器在轨几何定标方法及系统 | |
CN108072387B (zh) | 一种低精度敏感器在轨偏差的地面校正方法及系统 | |
Fei et al. | A novel approach based on MEMS-gyro's data deep coupling for determining the centroid of star spot | |
Wang et al. | Refraction surface-based stellar atmospheric refraction correction and error estimation for terrestrial star tracker | |
Samaan et al. | Non-Dimensional Star Identification for Un-Calibrated Star Cameras | |
CN115630254A (zh) | 高精度星敏感器光学系统内参数微变化量的在轨标定方法 | |
CN102117199A (zh) | 利用并行流水线机制提高星敏感器数据更新率的方法 | |
Xie et al. | A robust autonomous star identification algorithm for ZY3 satellite | |
CN111649738B (zh) | 微重力场下的加速度计初始姿态解算方法 | |
CN114646993A (zh) | 一种基于gnss、视觉以及imu进行精确定位的数据融合算法 | |
Wang et al. | A factor graph approach for attitude estimation of spacecraft using a stellar gyroscope |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |