专利名称：一种星敏感器恒星像点位置提取方法
技术领域：
本发明涉及一种航天测量领域中星敏感器星图中恒星像点位置提取方法；具体地 说是一种基于卡尔曼滤波的星敏感器恒星像点位置提取方法。
背景技术：
已有恒星像点位置提取方法，如质心法、曲面拟合法、神经网络法等。都是直接将 亚像元技术计算得到的结果作为最终的恒星像点位置，对于亚像元技术定位的随机噪声， 未采取有效的抑制措施。因此，虽然恒星像点位置提取的系统误差可以通过算法改进、标定 等方法消除，但随机噪声至今没有很好的解决。此外，恒星像点位置提取是星敏感器工作中 最耗时的程序之一。随着星敏感器工作频率的不断提高，如新一代APS星敏感器工作频率 可达IOHz 100Hz，恒星像点位置提取的速度也面临新的挑战。2002 年，期刊《The Journal of the Astronautical Sciences》第一期刊载 "Predictive Centroiding for Star Trackers with the Effect of Image Smear，，，提出 一种基于预测的恒星像点位置提取方法。该方法采用前一时刻星敏感器姿态矩阵和角速度 信息，预测当前时刻的恒星像点位置坐标，然后，从以这些预测的恒星坐标为中心的窗口内 提取相应的实际星像坐标。由于大大减少了星图阈值扫描的区域(只需扫描以预测的恒星 坐标为中心的窗口)，该方法可以有效提高恒星像点位置提取的速度。但该方法仍然是直 接将亚像元计算得到的恒星像点质心作为恒星像点位置，对随机噪声依然没有采取抑制措 施。

发明内容
本发明要解决的技术问题是针对现有星敏感器恒星像点位置提取方法存在的 问题，提出一种基于卡尔曼滤波的恒星像点位置提取方法。以实现在提高恒星像点位置提 取的速度的同时；消除随机噪声对恒星像点位置提取精度的影响；可进一步提高测量的精度。本发明的技术方案包括以下步骤(1)按现有技术利用全天球识别模式算法，获 得初始时刻星敏感器姿态信息、星图中恒星像点位置和瞬时速度；(2)根据前一帧星图中 恒星像点位置和瞬时速度，预测当前星图中恒星像点位置；(3)读入当前时刻星图；(4)采 用基于小面型的恒星像点位置提取方法，获得恒星像点位置的测量值；(5)利用卡尔曼滤 波更新恒星像点位置；(6)按照现有技术定位视场中新出现的恒星；(7)按照现有技术计算 当前时刻星敏感器姿态矩阵，并由陀螺数据或动力学方程，获得当前时刻的角速度，再回到 步骤2，进入下一循环。所述根据前一帧星图中恒星像点位置和瞬时速度，预测当前星图中恒星像点位 置，包括以下步骤；(1)确定恒星像点瞬时速度利用公式1，求得恒星像点在星敏感器像面上的瞬时运动速度办/Si):
4
式 1 中 其中，f表示星敏感器镜头焦距；(x, y)表示对应恒星的像点在星敏感器坐标系中的位置；qj, j = 1，2，3，4，为星敏感器姿态四元素；(2)预测恒星像点位置利用前一帧星图中恒星像点位置和瞬时速度，按以下公 式近似获得 式3中，是
时刻预测的恒星像点位置；[xtyt]τ和
分别表示t时刻该恒星像点的位置和瞬时速度；δ t星敏感器连续两帧星图的时间间隔。所述采用小面型的恒星像点位置提取方法，获得当前星图恒星像点位置的测量 值，包括以下步骤(1)星体阈值分割在以步骤2获得的恒星预测位置坐标为中心的10 X 10窗口内，提取星体目标
(4)式4中(x，y)为星光照射图像传感器像元的位置，I (x, y)为该像元的信号强度， T为星图噪声域值；(2)利用公式5，计算得到恒星像点位置坐标的测量值
(5)所述利用卡尔曼滤波更新恒星像点位置包括以下步骤(1)构建卡尔曼滤波的状态方程和测量方程
式6中，=4 Λ'，由步骤2获得;
叉⑴^,f，X{t + 5t) = \xl+Sl yl+Slf，分别表示t禾口 t+ δ t时刻的状态量； Z{t+St) = [x,+s, ，表示t时刻恒星像点位置测量值，由步骤3获得； Np和Nm分别表示过程噪声和测量噪声；
(2)通过卡尔曼滤波，获得当前时刻的恒星像点位置，计算方法如式7 X{t + St) = X(t) + H(t)St 式7中，Qe、Re和P分别表示过程噪声、测量噪声和最优估计噪声的协方差矩阵。K 是卡尔曼增益矩阵。所述按照现有技术定位视场中新出现的恒星，包括以下步骤(1)按照现有技术确定新进入视场的恒星像点区域；(2)新进入视场的恒星像点位置提取提取当前时刻新进入视场的恒星像点位置坐标，方法同步骤4 ；本发明的优点(1)利用前一帧星图中恒星像点位置和瞬时速度，预测当前星图 中恒星像点位置；从以这些预测的恒星坐标为中心的窗口内提取相应的星像坐标的测量 值。该方法可以有效提高恒星提取的速度；(2)构建了卡尔曼滤波器，以恒星像点位置的预 测值和测量值作为滤波器输入，最终的恒星像点位置坐标为输出，更新恒星像点位置，可有 效减小随机噪声。进一步提高测量的精度；本发明的测量精度较传统方法可提高2. 4倍。


图1是本发明整体流程图；图2是某恒星像点在星敏感器探测器成像面上运动轨迹示意图；图3是采用目前已有技术，用恒星像点位置测量值直接定位恒星在χ轴和y轴方 向的恒星像点位置误差；图4是采用星敏感器本发明的恒星像点位置提取方法，进行卡尔曼滤波后，χ轴和 y轴方向的恒星像点位置提取误差。
具体实施例方式采用本发明对星敏感器采集的星图中的恒星进行定位，具体步骤如下第一步，利用全天球识别模式算法，获得初始时刻星敏感器姿态信息、星图中恒星 像点位置和瞬时速度。第二步，预测恒星像点位置，方法是2. 1确定恒星像点瞬时速度2. 2预测恒星像点位置
第三步，读入当前时刻星图。第四步，基于小面型的恒星像点位置提取方法，获得当前星图恒星像点位置的测
量值，具体方法4. 1星体阈值分割4. 2利用公式5，计算得到恒星像点位置坐标的测量值第五步，利用卡尔曼滤波更新恒星像点位置第六步，定位视场中新出现的恒星，方法是6. 1确定新进入视场的恒星像点区域6. 2新进入视场的恒星像点位置提取第七步，计算当前时刻星敏感器姿态矩阵，并由陀螺数据或动力学方程，获得当前 时刻的角速度。一旦获得当前时刻星敏感器姿态矩阵和角速度，跳至步骤2，进入下一循环。在轨仿真及结果分析仿真星敏感器的基本参数设计为焦距76· 078mm视场10°X 10°图像传感器像元阵列1024 X 1024像元尺寸13· 3μπιΧ13· 3μπι初始时刻(t = 0)，卫星轨道参数为位置(6862315.59570592m, 4395. 21942281437m, 187. 670712716647m)速度(0.42045091053984m/s, -978. 995452142188m/s, 7562. 64902270816m/s)假设星敏感器在卫星上的方位角和俯仰角分别为-60°和90°。星敏感器三轴姿 态测量精度为(5〃，5〃，25〃)(1σ)，角速度测量精度为(0.05° /h,0. 05° /h,0. 05° / h) (1σ),噪声引起的恒星质心提取的随机误差为17 Urad(Io)0为了进一步验证算法的有效性，对卫星在轨运行0 1000s时间段内出现在星 敏感器视场内的所有恒星，采用本发明的方法定位；可消除随机噪声对恒星像点位置提取 精度的影响；可进一步提高测量的精度。如传统亚像元技术的恒星质心定位的随机误差 为17μ rad(l σ )，而本发明方法的提取误差为5 μ rad (1 σ )，其精度较传统方法提高了 2. 4倍。
权利要求
一种星敏感器恒星像点位置提取方法，其特征在于，包括以下步骤(1)按现有技术利用全天球识别模式算法，获得初始时刻星敏感器姿态信息、星图中恒星像点位置和瞬时速度；(2)根据前一帧星图中恒星像点位置和瞬时速度，预测当前星图中恒星像点位置；(3)读入当前时刻星图；(4)采用基于小面型的恒星像点位置提取方法，获得恒星像点位置的测量值；(5)利用卡尔曼滤波更新恒星像点位置；(6)按照现有技术定位视场中新出现的恒星；(7)按照现有技术计算当前时刻星敏感器姿态矩阵，并由陀螺数据或动力学方程，获得当前时刻的角速度，再返回到步骤2，进入下一循环。
2.根据权利要求1所述的星敏感器恒星像点位置提取方法，其特征在于，所述根据前 一帧星图中恒星像点位置和瞬时速度，预测当前星图中恒星像点位置，包括以下步骤；(1)确定恒星像点瞬时速度利用公式1，求得恒星像点在星敏感器像面上的瞬时运动速度(&/&，办/a) 式1中 预测的恒星像点位置； 分别表示t时刻该恒星像点的位置和瞬时速度； δt星敏感器连续两帧星图的时间间隔。
3.根据权利要求1所述的星敏感器恒星像点位置提取方法，其特征在于，所述采用/J 面型的恒星像点位置提取方法，获得当前星图恒星像点位置的测量值，包括以下步骤 (1)星体阈值分割在以步骤2获得的恒星预测位置坐标为中心的10 X 10窗口内，提取星体目标 其中，f表示星敏感器镜头焦距； (X，y)表示对应恒星的像点在星敏感器坐标系中的位置；态四元素； 一帧星图中恒星像点位置和瞬时速度，按以下公式近Qj, j = 1，2，3，4，为星敏感器姿 (2)预测恒星像点位置利用前 似获得 式3中，[U, λJr是t+St时亥Ij式4中(X，y)为星光照射图像传感器像元的位置，I (x, y)为该像元的信号强度，T为 星图噪声域值； (2)利用公式5，计算得到恒星像点位置坐标的测量值
4.根据权利要求1所述的星敏感器恒星像点位置提取方法，其特征在于，所述利用卡 尔曼滤波更新恒星像点位置包括以下步骤(1)构建卡尔曼滤波的状态方程和测量方程 式6中 由步骤2获得； 分别表示 t 禾Π t+ δ t 时刻的状态量； 表示t时刻恒星像点位置测量值，由步骤3获得； Np和Nm分别表示过程噪声和测量噪声；(2)通过卡尔曼滤波，获得当前时刻的恒星像点位置，计算方法如式7 式7中，Q。民和P分别表示过程噪声、测量噪声和最优估计噪声的协方差矩阵。K是 卡尔曼增益矩阵。
5.根据权利要求1所述的星敏感器恒星像点位置提取方法，其特征在于，所述按照现 有技术定位视场中新出现的恒星，包括以下步骤(1)按照现有技术确定新进入视场的恒星像点区域；(2)新进入视场的恒星像点位置提取提取当前时刻新进入视场的恒星像点位置坐 标，方法同步骤4。
全文摘要
本发明公开了一种星敏感器恒星像点位置提取方法。是根据前一帧星图中恒星像点位置和瞬时速度，预测当前星图中恒星像点位置；采用基于小面型的恒星像点位置提取方法，获得恒星像点位置的测量值；最后以恒星像点位置预测值和测量值为输入，构建卡尔曼滤波器，获得最终的恒星像点位置。本发明利用前一帧星图中恒星像点位置和瞬时速度，预测当前星图中恒星像点位置可以有效提高恒星提取的速度；又利用卡尔曼滤波更新恒星像点位置，可有效减小随机噪声；进一步提高测量的精度；本发明的测量精度较传统方法可提高2.4倍。
文档编号G01C21/02GK101907463SQ20101021664
公开日2010年12月8日 申请日期2010年7月5日 优先权日2010年7月5日
发明者刘海波, 宿德志, 廖家莉, 李修建, 杨俊才, 杨建坤, 王文亮, 谭吉春, 贾辉 申请人:中国人民解放军国防科学技术大学