CN104567867B - Method for stably tracking attitude of star sensor - Google Patents

Method for stably tracking attitude of star sensor Download PDF

Info

Publication number
CN104567867B
CN104567867B CN201410838410.3A CN201410838410A CN104567867B CN 104567867 B CN104567867 B CN 104567867B CN 201410838410 A CN201410838410 A CN 201410838410A CN 104567867 B CN104567867 B CN 104567867B
Authority
CN
China
Prior art keywords
star
tracking
node
tln
list
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.)
Active
Application number
CN201410838410.3A
Other languages
Chinese (zh)
Other versions
CN104567867A (en
Inventor
钟红军
杨孟飞
王龙
卢欣
杨君
梁士通
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Control Engineering
Original Assignee
Beijing Institute of Control Engineering
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Control Engineering filed Critical Beijing Institute of Control Engineering
Priority to CN201410838410.3A priority Critical patent/CN104567867B/en
Publication of CN104567867A publication Critical patent/CN104567867A/en
Application granted granted Critical
Publication of CN104567867B publication Critical patent/CN104567867B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/02Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
    • G01C21/025Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means with the use of startrackers

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Astronomy & Astrophysics (AREA)
  • Automation & Control Theory (AREA)
  • General Physics & Mathematics (AREA)
  • Navigation (AREA)

Abstract

The invention discloses a method for stably tracking the attitude of a star sensor. The method comprises steps as follows: firstly, establishing a tracking list structure in star sensor software, and recording tracking and identification conditions of star points; comparing and analyzing a matching identification result of each frame of star map of the star sensor and matching identification results of multiple frames of previous star maps to screen out the star points stably tracked by the star sensor; only adopting star pairs capable of being stably tracked by the star sensor to perform attitude calculation. According to the method, only the star pairs which pass multiple-frame verification and can be stably tracked and identified are adopted for attitude calculation, so that effects caused by output attitudes corresponding to star pairs which cannot be stably identified are eliminated.

Description

A kind of method of star sensor attitude stabilization tracking
Technical field
The invention belongs to spacecraft uses Star Sensor field, it is related to a kind of side of star sensor attitude stabilization tracking Method.
Background technology
Star sensor is a kind of imaging type spatial attitude measurement sensor, is extracted by the asterism to clapped star chart, is recognized And calculating, export star sensor three-axis attitude information, due to its have high accuracy, high reliability the features such as, be widely used in In all kinds of spacecrafts such as satellite, airship.Star sensor adopts the observation star extracted from window star chart by fixed star window trace Point information and navigational star table asterism information carry out star pattern matching identification, and then carry out attitude algorithm.Obtained by star pattern matching identification The star pair of observation star and nautical star composition, is the basis for carrying out attitude algorithm, directly affects the output attitude accuracy of star sensor This important indicator.
Existing star sensor is such as《The celestial body segmented positioning technique study of star map image in star sensor》(author:Wei Xin State, Zhang Guangjun etc.;Source:BJ University of Aeronautics & Astronautics's journal) etc. described in document, directly star pattern matching recognized more obtained by star To being completely used for attitude algorithm, due to nautical star point tolerance, the into/out visual field of asterism, sensitive chip noise, Space Particle spoke Penetrate, the impact of the factor such as veiling glare, there is unstable situation to composition per the star obtained by the identification of frame star pattern matching, cause with its work Attitude algorithm module output for input is unstable, affects the output attitude accuracy of star sensor.
The content of the invention
Present invention solves the technical problem that being:Overcome the deficiencies in the prior art, there is provided a kind of star sensor attitude stabilization with The method of track so that the star for attitude algorithm is stable to composition, so as to improve star sensor output attitude accuracy.
The technical scheme is that:A kind of method of star sensor attitude stabilization tracking, step is as follows:
1) star sensor is obtained after a new frame tracing mode star chart information, carries out asterism extraction and match cognization, is obtained Identification queue IR, described identification queue IR include that the navigational star number GSId and image planes of each navigation asterism extracts coordinate information;
2) the next node TLn in tracking list TL is obtained;Described tracking list TL be before per frame accumulation with Track list, including track star number tsn, determine appearance star number qsn and the navigational star number GSId of each tracking list node, continuous loss Whether number of times LostCycles, tracking number of times TrackedCycles and the node can be used to determine the mark FlgUseToCalc of appearance;
3) the next nautical star IRn in identification queue IR is taken;
4) whether the navigational star number IRn.GSId for comparing the navigational star number TLn.GSId and nautical star IRn of TLn is equal, if phase Deng this nautical star being rejected from identification queue IR, then execution step 6);Otherwise, execution step 5);
5) judge to recognize whether all nautical stars extract and terminate in queue IR, if so, then execution step 7);Otherwise, perform Step 3);
6) the tracking number of times TLn.TrackedCycles for tracking list node TLn adds one, continuously loses number of times TLn.LostCycles sets to 0, if judging TLn.TrackedCycles more than or equal to tracking threshold value minNumForTracked, The appearance mark FlgUseToCalc that determines for putting this node is true, and that otherwise puts this node determines appearance mark FlgUseToCalc for non-real; Judgement terminates rear execution step 8);Wherein minNumForTracked is positive integer;
7) judge the continuous loss number of times TLn.LostCycles for tracking list node TLn whether more than or equal to miss-threshold MaxNumForLosted, if so, empties the tracking list node, i.e., the node is set to NullTLNode;Otherwise, TLn.LostCycles adds one, and that puts this node determines appearance mark FlgUseToCalc for non-real, execution step 8);Wherein MaxNumForLosted is positive integer;
8) judge to track all nodes in list and whether extract to terminate, if so, execution step 9);Otherwise, return to step 2);
9) it is above-mentioned to be filled out with the node being cleared in remaining nautical star information filling tracking list TL in identification queue IR The continuous loss number of times for filling node is set to 0, and tracking number of times is set to 1, determines appearance mark and is set to non-real;
10) calculate and obtain non-NULL nodes A in tracking list TL, tracking star number tsn of tracking list TL is set to into A;Meter Calculation acquisition tracking determines appearance and is masked as genuine nodes B, and appearance star number qsn of determining of tracking list is set to into B;
11) using determining appearance and be masked as the attitude that genuine node nautical star information completes this frame of star sensor in tracking list TL Resolve and tenacious tracking.
Present invention good effect compared with prior art:
1. propose that star chart culminant star point mass multiframe compares evaluation method, realizes the preferred of match cognization result.Prior art Many angles that asterism vector and nautical star point vector are observed by adjusting present frame, reject when the larger identification star pair of frame error; The present invention on the basis of existing technology, increased the multiframe checking of identification star pair, and rejecting can not be by stable between continuous multiple frames Star pair with identification;
2. stably it is used for the star of attitude algorithm to composition, improves star sensor output attitude accuracy.Prior art is more complete Into star pattern matching identification after, directly by star pattern matching recognize obtained by star to being completely used for attitude algorithm, when a certain star is to more When interframe can not be recognized stably, star sensor output attitude error is larger;The present invention is steady only with verify through multiframe, energy The star of Tracking Recognition is determined to carrying out attitude algorithm, eliminating can not stably impact of the identification star to correspondence output attitude.
Description of the drawings
Fig. 1 is the main flow chart of the present invention.
Specific embodiment
The present invention has been applied to domestic certain small-sized star sensor, as shown in figure 1, specific implementation step is as follows:
1) star sensor is obtained after a new frame tracing mode star chart information, carries out asterism extraction and match cognization, is obtained Identification queue IR, described identification queue IR include that the navigational star number GSId and image planes of each navigation asterism extracts coordinate information;
2) the next node TLn in tracking list TL is obtained;Described tracking list TL be before per frame accumulation with Track list, including track star number tsn, determine appearance star number qsn and the navigational star number GSId of each tracking list node, continuous loss Whether number of times LostCycles, tracking number of times TrackedCycles and the node can be used to determine the mark FlgUseToCalc of appearance;
3) the next nautical star IRn in identification queue IR is taken;
4) whether the navigational star number IRn.GSId for comparing the navigational star number TLn.GSId and nautical star IRn of TLn is equal, if phase Deng this nautical star being rejected from identification queue IR, then execution step 6);Otherwise, execution step 5);
5) judge to recognize whether all nautical stars extract and terminate in queue IR, if so, then execution step 7);Otherwise, perform Step 3);
6) the tracking number of times TLn.TrackedCycles for tracking list node TLn adds one, continuously loses number of times TLn.LostCycles sets to 0, if judging TLn.TrackedCycles more than or equal to tracking threshold value minNumForTracked, The appearance mark FlgUseToCalc that determines for putting this node is true, and that otherwise puts this node determines appearance mark FlgUseToCalc for non-real; Judgement terminates rear execution step 8);Wherein minNumForTracked is positive integer;
7) judge the continuous loss number of times TLn.LostCycles for tracking list node TLn whether more than or equal to miss-threshold MaxNumForLosted, if so, empties the tracking list node, i.e., the node is set to NullTLNode;Otherwise, TLn.LostCycles adds one, and that puts this node determines appearance mark FlgUseToCalc for non-real, execution step 8);Wherein MaxNumForLosted is positive integer;
8) judge to track all nodes in list and whether extract to terminate, if so, execution step 9);Otherwise, return to step 2);
9) it is above-mentioned to be filled out with the node being cleared in remaining nautical star information filling tracking list TL in identification queue IR The continuous loss number of times for filling node is set to 0, and tracking number of times is set to 1, determines appearance mark and is set to non-real;
10) calculate and obtain non-NULL nodes A in tracking list TL, tracking star number tsn of tracking list TL is set to into A;Meter Calculation acquisition tracking determines appearance and is masked as genuine nodes B, and appearance star number qsn of determining of tracking list is set to into B;
11) using determining appearance and be masked as the attitude that genuine node nautical star information completes this frame of star sensor in tracking list TL Resolve and tenacious tracking.
Table one
Table two
For example, if tracking threshold value be 5, miss-threshold is 3, previous frame recognition result as shown in Table 1, step 1) in it As shown in Table 2, node 4,12,13 is unidentified due to previous frame, through step 7 for the tracking list of front accumulation) determine appearance be masked as it is non- Very;Node 10 due to tracking number of times less than tracking threshold value, through step 7) determine appearance and be masked as non-real;The tracking number of times of other nodes It is all higher than tracking threshold value, therefore through step 6) determine appearance mark and be very.Step 1) in this frame identification the recognized list such as institute of table three Show.Then through calculating, step 10) tracking list after the renewal that obtains as shown in Table 4, its interior joint 4 is recognized due to present frame And number of times is tracked more than tracking threshold value, therefore through step 6) determine appearance mark and be changed into true;Node 10 recognizes that tracking is secondary due to present frame Number has been equal to tracking threshold value, therefore through step 6) determine appearance mark and be changed into true;Node 11,12 is passed through due to present frame nonrecognition Step 7) determine appearance mark and be changed into non-real, and continuous number of times of losing adds 1;Node 13 is due to present frame nonrecognition and continuously loses number of times Equal to miss-threshold, therefore through step 7) it is cleared, and through step 9) it is filled to be new node.If using tracking list, on One frame step 11) carry out attitude algorithm, present frame step 11 with the interior joint 1,2,3,5,6,7,8,9,11 of table two) use table four Interior joint 1~10 carries out attitude algorithm, wherein the nautical star for determining appearance have three it is different.If not using tracking list, Previous frame and present frame carry out attitude algorithm with table one, recognized list shown in table three respectively, wherein the nautical star for determining appearance is total to Have 5 it is different.Understand, star sensor attitude stabilization tracking proposed by the present invention can be in the unstable situation of recognition result The lower tenacious tracking for realizing attitude.
Table three
Table four
The content not being described in detail in description of the invention belongs to the known technology of those skilled in the art.

Claims (1)

1. a kind of method that star sensor attitude stabilization is tracked, it is characterised in that step is as follows:
1) star sensor is obtained after a new frame tracing mode star chart information, carries out asterism extraction and match cognization, is identified Queue IR, described identification queue IR include that the navigational star number GSId and image planes of each navigation asterism extracts coordinate information;
2) the next node TLn in tracking list TL is obtained;The tracking row that described tracking list TL is accumulated before being per frame Table, including tracking star number tsn, determine appearance star number qsn and each track list node navigational star number GSId, continuously lose number of times Whether LostCycles, tracking number of times TrackedCycles and the node can be used to determine the mark FlgUseToCalc of appearance;
3) the next nautical star IRn in identification queue IR is taken;
4) whether the navigational star number IRn.GSId for comparing the navigational star number TLn.GSId and nautical star IRn of TLn is equal, if equal, This nautical star is rejected from identification queue IR, then execution step 6);Otherwise, execution step 5);
5) judge to recognize whether all nautical stars extract and terminate in queue IR, if so, then execution step 7);Otherwise, execution step 3);
6) the tracking number of times TLn.TrackedCycles for tracking list node TLn adds one, continuously loses number of times TLn.LostCycles sets to 0, if judging TLn.TrackedCycles more than or equal to tracking threshold value minNumForTracked, The appearance mark FlgUseToCalc that determines for putting this node is true, and that otherwise puts this node determines appearance mark FlgUseToCalc for non-real; Judgement terminates rear execution step 8);Wherein minNumForTracked is positive integer;
7) judge the continuous loss number of times TLn.LostCycles for tracking list node TLn whether more than or equal to miss-threshold MaxNumForLosted, if so, empties the tracking list node, i.e., the node is set to NullTLNode;Otherwise, TLn.LostCycles adds one, and that puts this node determines appearance mark FlgUseToCalc for non-real, execution step 8);Wherein MaxNumForLosted is positive integer;
8) judge to track all nodes in list and whether extract to terminate, if so, execution step 9);Otherwise, return to step 2);
9) it is above-mentioned to be filled section with the node being cleared in remaining nautical star information filling tracking list TL in identification queue IR The continuous loss number of times of point is set to 0, and tracking number of times is set to 1, determines appearance mark and is set to non-real;
10) calculate and obtain non-NULL nodes A in tracking list TL, tracking star number tsn of tracking list TL is set to into A;Calculating is obtained Must track and determine appearance and be masked as genuine nodes B, appearance star number qsn of determining of tracking list is set to into B;
11) using determining appearance and be masked as the attitude algorithm that genuine node nautical star information completes this frame of star sensor in tracking list TL And tenacious tracking.
CN201410838410.3A 2014-12-29 2014-12-29 Method for stably tracking attitude of star sensor Active CN104567867B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410838410.3A CN104567867B (en) 2014-12-29 2014-12-29 Method for stably tracking attitude of star sensor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410838410.3A CN104567867B (en) 2014-12-29 2014-12-29 Method for stably tracking attitude of star sensor

Publications (2)

Publication Number Publication Date
CN104567867A CN104567867A (en) 2015-04-29
CN104567867B true CN104567867B (en) 2017-05-10

Family

ID=53084483

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410838410.3A Active CN104567867B (en) 2014-12-29 2014-12-29 Method for stably tracking attitude of star sensor

Country Status (1)

Country Link
CN (1) CN104567867B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109271942A (en) * 2018-09-26 2019-01-25 上海七牛信息技术有限公司 A kind of stream of people's statistical method and system

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106403929B (en) * 2016-08-22 2019-02-05 清华大学 The fast robust method of calibration of star sensor importance in star map recognition and attitude algorithm

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102435763A (en) * 2011-09-16 2012-05-02 中国人民解放军国防科学技术大学 Measuring method for attitude angular velocity of spacecraft based on star sensor
CN103323027A (en) * 2013-05-30 2013-09-25 北京控制工程研究所 Star point reconstruction-based star sensor dynamic-compensation method
CN103337059A (en) * 2013-07-10 2013-10-02 北京控制工程研究所 High-energy particle interference suppression method of star sensor
CN103335648A (en) * 2013-06-27 2013-10-02 北京航天自动控制研究所 Autonomous star map identification method
CN103954280A (en) * 2014-04-08 2014-07-30 北京控制工程研究所 Rapid, high-robustness and autonomous fixed star identification method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102435763A (en) * 2011-09-16 2012-05-02 中国人民解放军国防科学技术大学 Measuring method for attitude angular velocity of spacecraft based on star sensor
CN103323027A (en) * 2013-05-30 2013-09-25 北京控制工程研究所 Star point reconstruction-based star sensor dynamic-compensation method
CN103335648A (en) * 2013-06-27 2013-10-02 北京航天自动控制研究所 Autonomous star map identification method
CN103337059A (en) * 2013-07-10 2013-10-02 北京控制工程研究所 High-energy particle interference suppression method of star sensor
CN103954280A (en) * 2014-04-08 2014-07-30 北京控制工程研究所 Rapid, high-robustness and autonomous fixed star identification method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
星敏感器高动态下自主星跟踪算法;李葆华等;《红外与激光工程》;20120131;第41卷(第1期);190-195 *
高精度星敏感器的导航星识别定位算法研究;李学夔等;《计算机工程与应用》;20101231;第46卷(第21期);220-223 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109271942A (en) * 2018-09-26 2019-01-25 上海七牛信息技术有限公司 A kind of stream of people's statistical method and system

Also Published As

Publication number Publication date
CN104567867A (en) 2015-04-29

Similar Documents

Publication Publication Date Title
CN102607526B (en) Target posture measuring method based on binocular vision under double mediums
CN103697855B (en) A kind of hull horizontal attitude measuring method detected based on sea horizon
CN101839721B (en) Visual navigation method in autonomous rendezvous and docking
CN103617328B (en) A kind of airplane three-dimensional attitude computation method
CN107833249A (en) A kind of carrier-borne aircraft landing mission attitude prediction method of view-based access control model guiding
EP2927870B1 (en) Calibration apparatus, calibration method, and calibration program
CN101364259A (en) Method for extracting road various information of multi-level knowledge driven panchromatic remote sensing image
CN104899590A (en) Visual target tracking method and system for unmanned aerial vehicle
CN105138983B (en) The pedestrian detection method divided based on weighting block model and selective search
CN106340010B (en) A kind of angular-point detection method based on second order profile difference
CN103075998A (en) Monocular space target distance-measuring and angle-measuring method
CN107677274A (en) Unmanned plane independent landing navigation information real-time resolving method based on binocular vision
CN109708643B (en) Evaluation and selection method for asteroid surface optical navigation road sign
CN110095123B (en) Method for evaluating and optimizing observation information of road signs on surface of irregular small celestial body
CN104792299B (en) Asteroid orbit identifying method based on observation angle data
CN104197933A (en) Method for enhancing and extracting high stars and other sliding fixed stars in view field of telescope
CN107798691A (en) A kind of unmanned plane independent landing terrestrial reference real-time detecting and tracking method of view-based access control model
CN106706133A (en) Spot-like target attitude estimation method and system
CN106327536A (en) Neck line measuring method based on sectional point clouds
CN108959182A (en) The small feature loss gravitational field modeling method returned based on Gaussian process
CN104567867B (en) Method for stably tracking attitude of star sensor
CN103854290A (en) Extended target tracking method based on combination of skeleton characteristic points and distribution field descriptors
CN103839274A (en) Extension target tracking method based on geometric proportion relation
CN104407345B (en) Improved dynamic planning method for weak target detection tracking
CN103942786A (en) Self-adaptation block mass target detecting method of unmanned aerial vehicle visible light and infrared images

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