CN101957203B - 一种星敏感器高精确度的星跟踪方法 - Google Patents
一种星敏感器高精确度的星跟踪方法 Download PDFInfo
- Publication number
- CN101957203B CN101957203B CN2010101931318A CN201010193131A CN101957203B CN 101957203 B CN101957203 B CN 101957203B CN 2010101931318 A CN2010101931318 A CN 2010101931318A CN 201010193131 A CN201010193131 A CN 201010193131A CN 101957203 B CN101957203 B CN 101957203B
- Authority
- CN
- China
- Prior art keywords
- star
- cos
- fixed
- sensor
- sub
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Navigation (AREA)
Abstract
本发明提供一种星敏感器高精确度的星跟踪方法。采用星跟踪方法的跟踪窗口方式,获得星敏感器视场内所有已知恒星星像坐标,根据这些已知恒星的信息,搜索出视场内所有未知恒星,把视场内所有恒星划分为多个不均匀的子区域,其中恒星密度大的子区域面积小,恒星密度小的子区域面积大,这样尽量满足每个子区域能有一定数量的恒星,并把这种子区域的分布、面积大小对应到相应的星图,把星图划分为数量相同、面积比例相同的子星图,根据获得已知恒星的星像坐标,确认是否存在不包含已知恒星的子星图中,而且只识别不包含已知恒星子星图内的未知恒星。从而满足每个区域内都有恒星参与计算姿态,保证了参与计算姿态的恒星在全视场内均匀分布。本发明的优点是减少了计算量,保证了星跟踪模式下的实时性,降低了星敏感器输出的姿态与实际姿态的误差。
Description
(一)技术领域
本发明涉及一种航天领域测量技术,具体涉及一种星敏感器高精确度的星跟踪方法。
(二)背景技术
星敏感器是一种以恒星作为观测基准的高精度飞行器姿态测量敏感器,能够提供角秒级甚至更高精度的姿态信息,被认为是目前精度最高的姿态敏感器。星敏感器有初始姿态捕获和星跟踪模式两种模式。初始姿态捕获就是在全天球范围内对CCD(Charge Coupled Device)或CMOS(Complementary Metal Oxide Semiconductor Transistor)成像器件所得到的观测星图进行匹配、识别。计算初始的精确姿态。随着航天技术的发展,对航天飞行器的自主导航能力提出了越来越高的要求,不仅要求准确的姿态信息,还要有很好的快速性。因此很多学者提出各种全天球星图识别算法或者改进各种全天球星图识别算法来提高星敏感器的快速性。这些算法虽然从某种程度上提高了姿态更新率,但是由于全天球星图识别前,必须从整幅星图中提取星像坐标,而提取星像坐标占了整个全天球星图识别模式的绝大部分时间。随着像平面的面阵增大,每次提取星像坐标过程的时间增加。然而宽视场、大面阵的像平面是新一代星敏感器的发展趋势。所以完全从全天球星图识别算法考虑来提高星敏感器的姿态更新率必然受到限制。
为了满足航天器对星敏感器实时性能的需要,以及随着星敏感器处理能力和存储能力的提高,星敏感器经过一次全天星图识别获得初始姿态后,必然有适当的先验信息,因此星敏感器可以不必连续工作在全天球识别模式,而是充分利用这些先验信息,工作在星跟踪模式。星跟踪模式就是利用这些先验信息可以在星图的适当区域内进行提取星像坐标,从而避免了在很多星图背景处判断是否存在星像的过程,这样进一步节省了星像坐标提取过程中的时间,从而提高了星敏感器数据更新率,提高了星敏感器的实时性。
星跟踪方法主要有跟踪窗口方式和匹配识别方式两种跟踪方式。传统的星跟踪方法大部分时间工作在跟踪窗口方式,不考虑参与计算姿态的恒星在视场内的分布,只有当星敏感器视场内已知恒星较少时,星敏感器才采用匹配识别跟踪方法来获取星敏感器视场内未知恒星信息,如图3。由于飞行器的运动,当星敏感器工作在纯跟踪窗口方式下,随着所有已知恒星从视场中心移到视场边沿,星敏感器输出的姿态与理想姿态的误差增加,而当星敏感器再次工作在匹配识别跟踪方法后,由于对视场内其它未知恒星识别,并参与姿态的计算,星敏感器输出的姿态与理想姿态的误差又减少了。如图1中a所示,随着已知恒星从t1时刻恒星星像位置到tn时刻恒星星像位置的移动,输出姿态从图2的输出姿态曲线的A点移到A1点, 当星敏感器在tn时刻对未知恒星识别后,如图1中b所示,又降低了姿态误差,如A2点。
(三)发明内容
本发明的目的在于提供一种新的星敏感器星跟踪方法,在相同的条件下,使星敏感器在星跟踪模式下不但具有很好的实时性,而且提供更高精度的姿态信息。
1.本发明的目的是这样实现的:方法如下:
(1)采用星跟踪方法的跟踪窗口方式,获得星敏感器视场内所有已知恒星星像坐标;
(2)根据这些已知恒星的信息,从星表中搜索出所有视场内所有恒星,包括已知恒星和未知恒星;
(3)根据星敏感器视场内恒星的分布,把视场内恒星划分为多个不均匀的子区域,其中恒星密度大的子区域面积小,恒星密度小的子区域面积大,这样尽量满足每个子区域能有一定数量的恒星,并把这种子区域的分布、面积大小对应到相应的星图,把星图划分为数量相同、面积比例相同的子星图,根据获得已知恒星的星像坐标,确认是否存在不包含已知恒星的子星图中;
(4)对于不包含已知恒星的子区域,采用从天球坐标系到图像平面坐标系中理想星像质心坐标的公式一和公式二,来计算未知恒星在子星图中的理想星像坐标;
式一
x′=xcosθ-ysinθ
式二
y′=xsinθ+ycosθ
其中:
(A,D)是理想坐标系中心,即光轴指向,在天球坐标系的指向;
θ为星敏感器的滚动角;
这些都可以通过已知恒星信息来计算得到;
(x′,y′)是计算未知恒星在子星图中的理想星像坐标;
(α,δ)是未知恒星在天球坐标系中的赤经和赤纬;
(5)采用星跟踪方法的匹配识别方式验证这些未知恒星的实际星像坐标,即新信息;
(6)保存这些新信息以及原有的信息,并利用这些新信息以及原有的已知信息,计算当前星敏感器姿态。
本发明的优点:
把星敏感器视场恒星的分为已知和未知,把整个视场划分为若干个面积不等的小区域,保证每个小区域的单位面积恒星相等;采用已知恒星的信息,对存在已知恒星的区域内的未知恒星不识别,仅仅识别不存在已知恒星区域内的未知恒星信息,这样减少了计算量,从而保证了星跟踪模式下的实时性;识别不存在已知恒星区域内的未知恒星信息,从而使参与计算姿态的恒星充满整个视场,并且在整个视场内均匀分布,从而降低了星敏感器输出的姿态与理想姿态的误差。
(四)附图说明
图1为传统星跟踪模式下预测、识别未知恒星信息示意图;
图2为传统星跟踪模式下星敏感器输出的姿态与实际姿态的误差曲线(某一轴);
图3为传统星跟踪模式下星敏感器视场内已知恒星数量与跟踪方式的关系;
图4为星敏感器视场内已知恒星和未知恒星的分布示意图;
图5为高精确度星跟踪方法流程图;
图6为传统星跟踪识别模式的三轴欧拉角误差曲线;
图7为高精确度星跟踪方法(即本方案)的三轴欧拉角误差曲线;
图8为本发明应用于某载体的工作流程框图;
图9为一种星敏感器高精确度的星跟踪方法原理图。
(五)具体实施方式
根据星敏感器仿真评价使用的姿态角精度公式:
其中σθ、σψ和σφ分别是星敏感器在俯仰角、偏航角和滚动角方向的姿态测量误差;n是参与计算的恒星数量;σxy是星敏感器像面上所成的像点与焦距的球心焦的测量角度误差;θsep为探测恒星的平均分离角度。
假设星敏感器像元尺寸和图像传感器面阵为正方为a×a,像平面面阵为N×N,那么形,像元尺寸
其中M为像元尺寸与测量误差的比值,f为星敏感器的焦距。
代入上式得:
由上式可知,姿态角的测量精度随着参与计算恒星的数量、星敏感器图像传感器面阵的增大而提高,随着视场的增加而降低,滚动角的精度还与参与计算恒星的分布状况有关。把星图分为若干个区域,由于飞行器的运动,某些区域已知恒星的数量可能增加,同时某些区域已知恒星数量可能减少。当某时刻某些区域已知恒星的数量为0时,立刻对视场内未知恒星进行识别,从使参与计算姿态的已知恒星分布充满整个星敏感器视场,并在星敏感器视场内尽量均匀分布,如图4。
本发明采用星跟踪方法的跟踪窗口方式,获得星敏感器视场内所有已知恒星星像坐标;根据这些已知恒星的信息,从星表中搜索出所有视场内所有恒星,包括已知恒星和未知恒星;根据星敏感器视场内恒星的分布,把视场内恒星划分为多个不均匀的子区域,其中恒星密度大的子区域面积小,恒星密度小的子区域面积大,这样尽量满足每个子区域能有一定数量的恒星,并把这种子区域的分布、面积大小对应到相应的星图,把星图划分为数量相同、面积比例相同的子星图,根据获得已知恒星的星像坐标,确认是否存在不包含已知恒星的子星图中;对于不包含已知恒星的子区域,采用从天球坐标系到图像平面坐标系中理想星像质心坐标的公式一和公式二,来计算未知恒星在子星图中的理想星像坐标;
式一
x′=xcosθ-ysinθ
式二
y′=xsinθ+ycosθ
其中:
(A,D)是理想坐标系中心,即光轴指向,在天球坐标系的指向;
θ为星敏感器的滚动角;
这些都可以通过已知恒星信息来计算得到;
(x′,y′)是计算未知恒星在子星图中的理想星像坐标;
(α,δ)是未知恒星在天球坐标系中的赤经和赤纬;
采用星跟踪方法的匹配识别方式验证这些未知恒星的实际星像坐标,即新信息;保存这些新信息以及原有的信息,并利用这些新信息以及原有的已知信息,计算当前星敏感器姿态。
主要性能指标:
我们选用某型号卫星星敏感器进行试验。试验方法:将星敏感器放在地球表面,初始时随机对准天球某区域与地球相对静止,并随着地球自转旋转,进行长时间运行。星敏感器全天球识别后自主进入星跟踪模式,对于同一幅星图,星敏感器首先采用传统的星跟踪识别方法,即不考虑星敏感器视场内已知恒星的分布,当视场内恒星数量少于一定时才识别未知恒星,输出三轴姿态测量值,然后采用本方法进行星跟踪模式,并输出三轴姿态测量值,把相应测量值与实际值作差,把这些差值实时保存到上位机,如图6和图7相应离线误差曲线。
Claims (1)
1.一种星敏感器高精确度的星跟踪方法,其特征在于,方法如下:
(1)采用星跟踪方法的跟踪窗口方式,获得星敏感器视场内所有已知恒星星像坐标;
(2)根据这些已知恒星的信息,从星表中搜索出所有视场内所有恒星,包括已知恒星和未知恒星;
(3)根据星敏感器视场内恒星的分布,把视场内恒星划分为多个不均匀的子区域,其中恒星密度大的子区域面积小,恒星密度小的子区域面积大,这样尽量满足每个子区域能有一定数量的恒星,并把这种子区域的分布、面积大小对应到相应的星图,把星图划分为数量相同、面积比例相同的子星图,根据获得已知恒星的星像坐标,确认是否存在不包含已知恒星的子星图;
(4)对于不包含已知恒星的子星图,采用从天球坐标系到图像平面坐标系中理想星像质心坐标的公式一和公式二,来计算未知恒星在子星图中的理想星像坐标;
(1)
x′=xcosθ-ysinθ
(2)
y′=xsinθ+ycosθ
其中:
(A,D)是理想坐标系中心,即光轴指向,在天球坐标系的指向;
θ为星敏感器的滚动角;
这些都可以通过已知恒星信息来计算得到;
(x′,y′)是计算未知恒星在子星图中的理想星像坐标;
(α,δ)是未知恒星在天球坐标系中的赤经和赤纬;
(5)采用星跟踪方法的匹配识别方式验证这些未知恒星的实际星像坐标,即新信息;
(6)保存这些新信息以及原有的信息,并利用这些新信息以及原有的已知信息,计算当前星敏感器姿态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101931318A CN101957203B (zh) | 2010-06-07 | 2010-06-07 | 一种星敏感器高精确度的星跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101931318A CN101957203B (zh) | 2010-06-07 | 2010-06-07 | 一种星敏感器高精确度的星跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101957203A CN101957203A (zh) | 2011-01-26 |
CN101957203B true CN101957203B (zh) | 2012-05-23 |
Family
ID=43484672
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101931318A Expired - Fee Related CN101957203B (zh) | 2010-06-07 | 2010-06-07 | 一种星敏感器高精确度的星跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101957203B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102252678B (zh) * | 2011-04-18 | 2013-01-23 | 北京航空航天大学 | 高动态高更新率星敏感器及其实现方法 |
CN102252664B (zh) * | 2011-04-18 | 2013-01-23 | 北京航空航天大学 | 一种恒星陀螺及其实现方法 |
CN102279002B (zh) * | 2011-06-27 | 2013-06-05 | 哈尔滨工业大学 | 星敏感器测量坐标系与载体坐标系转换矩阵的标定方法 |
CN102506829B (zh) * | 2011-11-15 | 2013-08-21 | 中国科学院紫金山天文台 | 快速星图匹配方法 |
CN102749089B (zh) * | 2012-07-13 | 2015-02-11 | 哈尔滨工业大学 | 一种三探头星敏感器姿态的确定方法 |
CN103323027B (zh) * | 2013-05-30 | 2015-07-08 | 北京控制工程研究所 | 基于星点重构的星敏感器动态补偿方法 |
CN103837160B (zh) * | 2014-03-04 | 2016-08-17 | 中国科学院光电技术研究所 | 基于恒星匹配的光电经纬仪参量优化方法 |
CN104202513B (zh) * | 2014-06-30 | 2017-05-31 | 北京控制工程研究所 | 一种基于fpga的星敏感器多模式图像预处理方法 |
CN106382927B (zh) * | 2016-08-19 | 2019-10-18 | 哈尔滨工业大学 | 一种基于卫星识别的星敏感器自主导航方法 |
CN106896827B (zh) * | 2017-04-06 | 2020-03-20 | 中国科学院光电技术研究所 | 一种光电跟踪测量设备指向位置附近恒星实时搜索的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101196398A (zh) * | 2007-05-25 | 2008-06-11 | 北京航空航天大学 | 一种基于Euler-q算法和DD2滤波的航天器姿态确定方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100344937C (zh) * | 2005-06-17 | 2007-10-24 | 北京航空航天大学 | 一种星***快速匹配识别方法 |
CN100348460C (zh) * | 2005-10-12 | 2007-11-14 | 北京航空航天大学 | 一种基于星场的星敏感器校准方法 |
US7451022B1 (en) * | 2006-12-28 | 2008-11-11 | Lockheed Martin Corporation | Calibration of ship attitude reference |
CN100529667C (zh) * | 2007-12-26 | 2009-08-19 | 北京控制工程研究所 | 一种自主恢复轨控故障时的星敏感器定姿方法 |
CN101435704B (zh) * | 2008-12-04 | 2010-06-16 | 哈尔滨工业大学 | 一种星敏感器高动态下的星跟踪方法 |
CN101701822B (zh) * | 2009-11-06 | 2011-08-24 | 中国人民解放军国防科学技术大学 | 基于光学联合变换相关的星敏感器星跟踪方法 |
-
2010
- 2010-06-07 CN CN2010101931318A patent/CN101957203B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101196398A (zh) * | 2007-05-25 | 2008-06-11 | 北京航空航天大学 | 一种基于Euler-q算法和DD2滤波的航天器姿态确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101957203A (zh) | 2011-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101957203B (zh) | 一种星敏感器高精确度的星跟踪方法 | |
CN102496015B (zh) | 二维高斯分布光斑图像中心快速高精度定位方法 | |
Zhang | Star identification | |
Zhang et al. | Vision-based pose estimation for textureless space objects by contour points matching | |
CN105371870B (zh) | 一种基于星图数据的星敏感器在轨精度测量方法 | |
CN101435704B (zh) | 一种星敏感器高动态下的星跟踪方法 | |
US10032273B2 (en) | Machine vision system calibration using inaccurate calibration targets | |
CN106643792B (zh) | 惯性测量单元和地磁传感器整体标定装置及标定方法 | |
CN106529587B (zh) | 基于目标点识别的视觉航向识别方法 | |
CN102155945B (zh) | 一种提高ccd星敏感器动态性能的方法 | |
CN107862719A (zh) | 相机外参的标定方法、装置、计算机设备和存储介质 | |
CN104268935A (zh) | 一种基于特征的机载激光点云与影像数据融合***及方法 | |
CN103697855B (zh) | 一种基于海天线检测的船体水平姿态测量方法 | |
CN103791902B (zh) | 适用于高机动载体的星敏感器自主导航方法 | |
CN106679676B (zh) | 一种单视场多功能光学敏感器及实现方法 | |
CN103136525B (zh) | 一种利用广义Hough变换的异型扩展目标高精度定位方法 | |
CN105444778B (zh) | 一种基于成像几何反演的星敏感器在轨定姿误差获取方法 | |
CN103363987B (zh) | 一种多视场星敏感器的星图识别方法 | |
CN104634248A (zh) | 一种双目视觉下的转轴标定方法 | |
CN102853835A (zh) | 基于尺度不变特征变换的无人飞行器景象匹配定位方法 | |
CN102914306A (zh) | 一种双探头星敏感器及其设计方法 | |
CN104949673A (zh) | 一种基于非视觉感知信息的目标定位方法及装置 | |
CN102944227A (zh) | 一种基于fpga实现实时提取恒星星像坐标的方法 | |
CN113295171B (zh) | 一种基于单目视觉的旋转刚体航天器姿态估计方法 | |
CN110030979A (zh) | 一种基于序列图像的空间非合作目标相对位姿测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120523 |