CN103487055B - 一种星惯组合还原目标空间分布的方法 - Google Patents
一种星惯组合还原目标空间分布的方法 Download PDFInfo
- Publication number
- CN103487055B CN103487055B CN201310392353.6A CN201310392353A CN103487055B CN 103487055 B CN103487055 B CN 103487055B CN 201310392353 A CN201310392353 A CN 201310392353A CN 103487055 B CN103487055 B CN 103487055B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- star sensor
- star
- red used
- impact point
- 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
Classifications
-
- 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/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
一种星惯组合还原目标空间分布的方法属于天文导航技术领域,通过结合航天器的实时姿态和星敏感器的安装矩阵,计算星敏感器坐标系与赤惯坐标系间的姿态转换矩阵;对星敏感器提取的目标点信息构造归一化的观测矢量;根据星敏感器在赤惯坐标系下的姿态转换矩阵,计算目标点在赤惯坐标系下的参考矢量,并计算目标在赤惯坐标系以赤经和赤纬定义的二维坐标系下的坐标;重复上述步骤,将多帧图像的目标点投影到天球二维坐标系,叠加在一起,得到某一时间段内星敏感器视场内目标点在赤惯坐标系下的分布图。本发明解决了事后分析星敏感器在轨期间数据时,判别视场内目标类型工作量大,且结果显示不直观的问题,以简单的方式有效地为结果分析提供基础数据。
Description
技术领域
本发明属于天文导航技术领域,具体涉及一种星惯组合还原目标空间分布的方法。
背景技术
随着科学技术的发展,人类争夺资源的加剧,航天技术的重要性显著增加。可以预见,航天事业将成为支撑人类科学研究、日常生活的支柱产业。在航天技术拥有的众多环节中,航天器姿态确定是其中最基础、最重要的任务之一。
姿态确定由姿态敏感器完成,如经典的星/惯组合导航方式中的星敏感器和惯性导航部件,状态良好的姿态敏感器可以很好地完成姿态确定任务。然而遗憾的是,航天器所处的空间物理环境并非理想的,在利用星敏感器校正惯导***,补偿惯性器件漂移的过程中,其视场中出现的景物除了有效的观测目标——恒星外,还包括如行星、星云、其它空间飞行器、空间碎片、地月形成的杂散光以及航天器自身的喷发物等等,这些都属于非预期的假目标。
星敏感器在轨期间拍摄到的目标属于哪一类,我们无法实时分析,仅能在事后对数据进行分析判别。分析得到的结果,不仅能使我们更了解空间环境,也能够使我们的空间飞行器及星敏感器更好地工作,从而更适应这种空间环境。当然,原始星图对于分析这些目标最为直观,但是由于遥测带宽的限制,无法传输完整的原始星图,仅能够将星敏感器实时提取的目标点信息通过遥测资源下传。
此时,如何利用已有的有限的资源进行分析?分析方法就显得非常重要,通常,可以结合航天器的实时姿态数据,采用模拟星图、速度匹配等方法,但这些方法不直观,而且工作量大。
发明内容
为了解决现有技术中存在的问题,得到直观、有效的分析结果,并减小分析的工作量,本发明提供了一种星惯组合还原目标空间分布的方法,该方法原理简单,使用方便。
本发明解决技术问题所采用的技术方案如下:
一种星惯组合还原目标空间分布的方法,该方法包括如下步骤:
步骤一:结合航天器在赤惯坐标系的姿态和星敏感器在航天器上的安装矩阵,计算星敏感器坐标系与赤惯坐标系间的姿态转换矩阵;
步骤二:对星敏感器提取的目标点信息构造归一化的观测矢量;
步骤三:根据星敏感器在赤惯坐标系下的姿态转换矩阵,计算目标点在赤惯坐标系下的参考矢量,并计算目标在赤惯坐标系以赤经和赤纬定义的二维坐标系下的坐标;
步骤四:重复步骤一至步骤三,将多帧图像的目标点投影到赤惯坐标系,叠加在一幅图像中,得到某一时间段内星敏感器视场内目标点在赤惯坐标系下的分布图。
本发明的有益效果是:本发明将多帧星敏感器提取的目标数据利用姿态变换的方式还原到赤惯坐标系,并转换为可视的图像,结果显示直观,方法简单,易于实现,工作量小,对分析飞行任务中星敏感器视场内的目标具有非常实用的价值和指导性的作用。
附图说明
图1为单帧星敏感器数据变换到赤惯坐标系下的流程图。
图2为星敏感器像平面坐标系和像空间坐标系。
图3为星敏感器视场内的目标还原到赤惯坐标系下的叠加图像。
具体实施方式
下面结合附图和实施例对本发明做进一步详细说明。
如图1所示,为单帧星敏感器数据变换到赤惯坐标系下的流程图,具体包含步骤一~步骤三。
步骤一:
已知航天器在赤惯坐标系下的姿态角和姿态转换矩阵,结合星敏感器在航天器上的安装矩阵,计算星敏感器坐标系与赤惯坐标系间的姿态转换矩阵;
假设星敏感器在航天器上的安装矩阵为Mss,航天器在赤惯坐标系下的姿态转换矩阵为Mst,那么,星敏感器在赤惯坐标系下在姿态转换矩阵为Mxm=Mst×Mss。
步骤二:
对星敏感器提取的目标点信息构造归一化的观测矢量;
星敏感器坐标系如图2所示,图中O-uv平面为星敏感器像平面坐标系,Os-xyz为星敏感器像空间坐标系,OOs为星敏感器焦距f,假设目标点Si在星敏感器像面坐标系中的坐标为(ui,vi),目标点Si在星敏感器像空间坐标系的观测矢量归一化后,即为
步骤三:
根据星敏感器在赤惯坐标系下的姿态转换矩阵,计算目标观测矢量在赤惯坐标系下的参考矢量,并计算目标在赤惯坐标系以(赤经,赤纬)定义的二维坐标系下的坐标;
根据步骤一中得到的星敏感器在赤惯坐标系下的姿态转换矩阵,将星敏感器像空间坐标系中的观测矢量Si转换到赤惯坐标系下,得到参考矢量Ri。
(xiyizi)是恒星在赤道直角坐标系下的三维坐标,(αiδi)是恒星在赤惯坐标系下的二维坐标。则该目标在赤惯二维坐标系下的(赤经,赤纬)坐标为
步骤四:
重复步骤一~步骤三,将多帧目标点投影到赤惯坐标系,叠加在一幅图像中,得到某一时间段星敏感器视场内目标点在赤惯坐标系下的分布图,如图3所示。图中,坐标横轴为赤经,坐标纵轴为赤纬,单位为度,“+”为星表中亮于5.5Mv的恒星,“*”为图像中的目标点还原到赤惯坐标系下的坐标,“◆”表示星敏感器的视轴指向。
图3所示的计算结果已经可以为分析进入星敏感器视场内的目标类型提供基础数据,并作初步判断。若图像中的目标点为恒星,通过该方法还原到赤惯坐标系后,应稳定地位于星表中某一恒星附近;而其它类型的目标,初步结果可以通过叠加后的图像,分析目标点在赤惯坐标系下的运动轨迹来判断,因惯性空间中自由飞行的物体不受或受很小的引力作用,按惯性定律,物体将做匀速直线运动;那些运动轨迹与星敏感器视轴运动轨迹一致的目标,则初步判断为航天器的伴飞物或喷发物。而更进一步的判断则要结合其它数据一起分析。
Claims (1)
1.一种星惯组合还原目标空间分布的方法,其特征在于,该方法包括如下步骤:
步骤一:结合航天器在赤惯坐标系的姿态和星敏感器在航天器上的安装矩阵,计算星敏感器坐标系与赤惯坐标系间的姿态转换矩阵;
步骤二:对星敏感器提取的目标点信息构造归一化的观测矢量;
步骤三:根据星敏感器在赤惯坐标系下的姿态转换矩阵,计算目标点在赤惯坐标系下的参考矢量,并计算目标在赤惯坐标系以赤经和赤纬定义的二维坐标系下的坐标;
步骤四:重复步骤一至步骤三,将多帧图像的目标点投影到赤惯坐标系,叠加在一幅图像中,得到某一时间段内星敏感器视场内目标点在赤惯坐标系下的分布图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310392353.6A CN103487055B (zh) | 2013-09-02 | 2013-09-02 | 一种星惯组合还原目标空间分布的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310392353.6A CN103487055B (zh) | 2013-09-02 | 2013-09-02 | 一种星惯组合还原目标空间分布的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103487055A CN103487055A (zh) | 2014-01-01 |
CN103487055B true CN103487055B (zh) | 2016-03-02 |
Family
ID=49827467
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310392353.6A Expired - Fee Related CN103487055B (zh) | 2013-09-02 | 2013-09-02 | 一种星惯组合还原目标空间分布的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103487055B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108645401B (zh) * | 2018-04-03 | 2020-05-22 | 中国人民解放军国防科技大学 | 基于姿态关联图像叠加的全天时星敏感器星点提取方法 |
CN109724588A (zh) * | 2018-12-25 | 2019-05-07 | 北京控制工程研究所 | 一种白昼星敏感器 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103047985A (zh) * | 2012-12-21 | 2013-04-17 | 中国科学院长春光学精密机械与物理研究所 | 一种空间目标的快速定位方法 |
-
2013
- 2013-09-02 CN CN201310392353.6A patent/CN103487055B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103047985A (zh) * | 2012-12-21 | 2013-04-17 | 中国科学院长春光学精密机械与物理研究所 | 一种空间目标的快速定位方法 |
Non-Patent Citations (2)
Title |
---|
CCD星图模拟器的设计及验证;郭敬明等;《中国光学与应用光学》;20101031;第3卷(第5期);486-493 * |
船用星敏感器星图模拟方法;李辉等;《光电工程》;20130430;第40卷(第4期);72-79 * |
Also Published As
Publication number | Publication date |
---|---|
CN103487055A (zh) | 2014-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Schwadron et al. | Determination of interstellar He parameters using five years of data from the IBEX: beyond closed form approximations | |
Colaninno et al. | Quantitative comparison of methods for predicting the arrival of coronal mass ejections at Earth based on multiview imaging | |
Millward et al. | An operational software tool for the analysis of coronagraph images: Determining CME parameters for input into the WSA‐Enlil heliospheric model | |
CN102879014B (zh) | 深空探测接近过程的光学成像自主导航半物理仿真试验*** | |
CN103344256B (zh) | 一种多视场星敏感器实验室测试方法 | |
McWilliams | Fundamentals of geophysical fluid dynamics | |
CN102354123B (zh) | 一种跨平台可扩展的卫星动态仿真测试*** | |
CN102519455B (zh) | 基于紫外敏感器的自主导航半物理仿真试验*** | |
CN103245364B (zh) | 一种星敏感器动态性能测试方法 | |
Luhmann et al. | Coupled model simulation of a Sun-to-Earth space weather event | |
CN108917772A (zh) | 基于序列图像的非合作目标相对导航运动估计方法 | |
CN106096204A (zh) | 一种基于太阳帆推进技术的航天器日心椭圆悬浮轨道设计方法 | |
CN103487055B (zh) | 一种星惯组合还原目标空间分布的方法 | |
Liu et al. | An improved quaternion Gauss–Newton algorithm for attitude determination using magnetometer and accelerometer | |
Horváth et al. | Evolution of an atmospheric Kármán vortex street from high‐resolution satellite winds: Guadalupe Island case study | |
Hietala et al. | Sequentially released tilted flux ropes in the Earth's magnetotail | |
Zhao et al. | Intelligent recognition of spacecraft components from photorealistic images based on Unreal Engine 4 | |
Pecora et al. | Multipoint turbulence analysis with helioswarm | |
Petrinec et al. | The steepness of the magnetic shear angle “saddle”: A parameter for constraining the location of dayside magnetic reconnection? | |
CN105183959A (zh) | 一种获取卫星在轨进行实传时的数传天线遮挡概率的方法 | |
Wang et al. | The dipole tilt angle dependence of the bow shock for southward IMF: MHD results | |
Amiri Atashgah et al. | An integrated virtual environment for feasibility studies and implementation of aerial MonoSLAM | |
Vernieres et al. | Capturing eddy shedding in the Gulf of Mexico from Lagrangian observations | |
Simonnet et al. | Bifurcation analysis of ocean, atmosphere, and climate models | |
CN104598731A (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 |
Granted publication date: 20160302 Termination date: 20170902 |
|
CF01 | Termination of patent right due to non-payment of annual fee |