CN101236090A - Moon-tracking segment ultraviolet sensor earth and moon photography time predication method - Google Patents

Moon-tracking segment ultraviolet sensor earth and moon photography time predication method Download PDF

Info

Publication number
CN101236090A
CN101236090A CNA2008100573474A CN200810057347A CN101236090A CN 101236090 A CN101236090 A CN 101236090A CN A2008100573474 A CNA2008100573474 A CN A2008100573474A CN 200810057347 A CN200810057347 A CN 200810057347A CN 101236090 A CN101236090 A CN 101236090A
Authority
CN
China
Prior art keywords
moon
earth
vector
time
heart
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
Application number
CNA2008100573474A
Other languages
Chinese (zh)
Other versions
CN100585334C (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 CN200810057347A priority Critical patent/CN100585334C/en
Publication of CN101236090A publication Critical patent/CN101236090A/en
Application granted granted Critical
Publication of CN100585334C publication Critical patent/CN100585334C/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a time prediction method for earth-moon picture shooting by an ultraviolet sensor at cislunar fight, which comprises steps of (1) selecting an orbit arc meeting the light condition for earth-moon shooting, (2) selecting an orbit arc, in which earth-moon appears, within viewing field area of the ultraviolet sensor, (3) taking the orbit position simultaneously meeting the step (1) and the step (2) as a selectable orbit for shooting, inputting an orbit prediction position corresponding to the selectable orbit, attitude matrix at a certain time before the selectable orbit, sun-moon-earth ephemeredes and a satellite attitude motion model after launching the satellite and before the coming of the selectable orbit, calculating the attitude matrix and a selenocentric/geocentric vector within the prediction time period and making a curve, changing along the time, of an included angle of the selenocentric/geocentric vector and a viewing field central vector, then obtaining continuous visible time range(t1-t2) of earth/moon by the curve and the size of the viewing field of the ultraviolet sensor and calculating the best time according to the obtained continuous visible time range.

Description

The moon-tracking segment ultraviolet sensor earth and moon photography time predication method
Technical field
The present invention relates to a kind of Benyue section and adopt ultraviolet sensors to carry out the method for ground month time predication, belong to optical imagery attitude sensor field.
Background technology
Along with the quick progress of imaging detection device and processor technology, the spacecraft attitude sensor is developed to imaging type by the unit scan formula gradually, and ultraviolet sensors is exactly a kind of big view field imaging formula attitude sensor that is different from traditional horizon instrument.
Ultraviolet sensors possesses the ability that the Benyue section is obtained the ultraviolet earth, ultraviolet moon image, not only important value but also utilization are obtained on science data the ground moon image information and cooperate the inertia attitude for realizing that independent navigation has important value.
At first ultraviolet sensors has special combined type visual field distribution, and the celestial sphere that forms distinguish covers, and therefore needs analysis ultraviolet sensors visual field model; The earth, the moon belong to the reflected sunlight imaging in addition, and shooting condition also needs certain sun, the earth/moon, satellite position demand; Also have, predict that the moon/earth appears at the imaging that the short accurate prediction that needs the shooting point moment of time visual field in could realize the earth/moon certainly not influencing under the satellite motion mode state to take to change next step motion of satellite and attitude.
Ultraviolet sensors itself is a brand-new imaging type attitude sensor, and this patent content is the expanded application research of carrying out on the ultraviolet sensors basis.In addition, the Benyue section is utilized ultraviolet sensors to carry out the moon/earth technique for taking and is verified it also is a new task, and has passed through in the rail flight test.Utilize ultraviolet sensors that the accurate prediction algorithm of the moon/earth shooting time be there is no domestic and international report.
External techwatch [J], 2006 the 9th phases, talk about Smart-1 in " No. 1, first lunar excursion module wisdom of Europe " and utilize the high precision camera moon to be carried out imaging, but do not specify segmental arc selection criterion, the best prediction mode method constantly of taking in the Benyue section.
A elfving, L Stagnearo, A Winton, introduced the situation that the Benyue section moon is taken camera that is used among the SMART-1:Key technologies andautonomy implementations.Acta Astronautic 52 (2003), not enough system of selection at undeclared shooting segmental arc, best shooting time.
U.S. Pat 5319969 has been introduced a kind of three-axis attitude of ultraviolet spectral coverage attitude sensor that utilizes in the title " Method for determining 3-axis spacecraftattitude " and has been determined method, does not wherein relate to the shooting forecasting problem of Benyue section.
Summary of the invention
Technology of the present invention is dealt with problems and is: overcome the deficiencies in the prior art, a kind of moon-tracking segment ultraviolet sensor earth and moon photography time predication method is provided, this method is taken ground month ultraviolet image according to the ground transmission instruction that predicts the outcome in the best time, reaches not change the shooting of satellite motion pattern realization to the earth/moon.
Technical solution of the present invention is: the moon-tracking segment ultraviolet sensor earth and moon photography time predication method comprises the following steps:
(1) carries out the illumination analysis according to a solar vector and a month heart vector in nominal track, the satellite body system, select the track segmental arc of taking illumination condition by the moon contentedly;
(2) obtain to satisfy the earth under the orbital arc fragment position in the step (1)/moon heart vector set according to the Satellite Attitude Movement of design in advance, whether in described vector set have vector be present in ultraviolet sensors field of view in, select to satisfy the track segmental arc in the ultraviolet sensors field of view if judging;
(3) will satisfy the orbital position of step (1) and step (2) simultaneously as taking optional track, before described optional track arrives, import the predicted orbit position of described optional track correspondence, optional track a certain moment attitude matrix, life ground ephemeris and Satellite Attitude Movement model before after the satellites transmits; Calculate attitude matrix, month heart/the earth's core vector in the predicted time section, and make a moon heart/the earth's core vector and the time dependent curve of visual field center vector angle, obtain the visible time range [t of the continuous earth/moon by curve and ultraviolet sensors visual field size 1~t 2], according to the continuously visible time range calculating optimum that obtains constantly.
Illumination analytic process in the described step (1) is: calculate solar vector V SWith moon heart vector V MOr the earth's core vector V EBetween angle , computing formula is:
=acos(V S·V i),
As > T, think under this orbital position that contentedly month taking an illumination condition obtains the better phases of the moon; Wherein, V iRepresent a month heart vector V MOr the earth's core vector V E
TBe the angle threshold value, size is determined according to the shooting demand.
The computing formula of attitude matrix is in the described step (3):
C t=R t·C 0
Wherein, R tBe the posture changing matrix, determine according to the Satellite Attitude Movement model;
C 0A certain moment attitude matrix for input.
Calculating month heart/the earth's core vector formula in the described step (3) is:
M t=R t·M 0
E t=R t·E 0
Wherein, M t-t moon heart vector constantly;
E t-t the earth's core vector constantly;
R tBe the posture changing matrix, determine according to the Satellite Attitude Movement model;
M 0The a certain moment t of-input 0Corresponding moon heart vector.
Visible time range deterministic process is as follows in the described step (3):
At first, calculate t moon heart vector M constantly tWith the earth's core vector E tThe angle  of the center vector Fs of ultraviolet sensors visual field, computing formula is:
M(t)=acos(Fs·M t)
E(t)=acos(Fs·E t)
Then, with the value and the visual field threshold value T of curve correspondence FOVCompare, work as  M<T FOV, the moon is described as seen, otherwise invisible; Work as  E<T FOV, the earth is described as seen, otherwise invisible; Continuous visible time border is visible time range; Wherein, T FOVEqual the visual field size of ultraviolet sensors optical system for the visual field threshold value.
The best time is got described visible time range [t in the described step (3) 1~t 2] intermediate value, or [t 1~t 2] the angle  moment corresponding of time range minimum.
The present invention compared with prior art beneficial effect is:
(1) the technology of the present invention is divided into ground moon bat figure prediction the preceding track selection work of emission and launches back prediction work in time, and according to field angle optimum shooting point is predicted, solved the problem of utilizing the success of ultraviolet sensors Benyue stage to take the earth, the moon at rail.
(2) the present invention carries out the illumination analysis by the angle that calculates between solar vector and month heart vector/the earth's core vector, has guaranteed that the moon/earth image possesses higher image-forming condition and science construction value in the selected shooting segmental arc.
(3) the present invention utilizes the method in a moon heart/the earth's core vector and visual field center vector angle analysis means acquisition continuous imaging interval to have and calculates simple, complete advantage.
(4) the technology of the present invention angle method placed in the middle or minimum service time can obtain optimum shooting moment point, guarantees to photograph the moon/earth and can photograph optimal image.
Description of drawings
Fig. 1 is the technology of the present invention algorithm flow chart;
Fig. 2 is an embodiment of the invention moon heart vector and visual field center vector angle curve;
Fig. 3 for the embodiment of the invention moon heart vector be in the visual field inner curve continuously;
Fig. 4 takes emulating image for our bright embodiment.
Embodiment
It is disclosed a kind of three-axis attitude sensor that utilizes the ultraviolet spectral coverage in the US5837894 title " Wide Field of View Sensor with diffractive OpticalCorrector " that the ultraviolet sensors that relates among the present invention can adopt the patent No. of Honeywell company application.Can also adopt disclosed a kind of three-axis attitude sensor that utilizes the ultraviolet spectral coverage in U.S. Pat 5319969 titles " Method for determining 3-axis spacecraft attitude ".
Below the moon is taken pictures and specifically describe for example.
As shown in Figure 1, be method flow diagram of the present invention, following mask body is introduced the implementation procedure of this method.
(1) according to the nominal track (X, Y, Z), solar vector and a month heart vector carry out the illumination analysis in the satellite body system, select to satisfy the track segmental arc that the earth/moon is taken illumination condition;
The illumination analytic process is: according to solar vector V S, month heart vector V MCalculate angle β:
β=acos(V S·V M)
As β>β T, think that satisfying illumination condition under this orbital position obtains better lunar map picture, β TFor the phase angle threshold value is selected according to demand, 90 ° of corresponding crescent or half moons wherein, and generally require greater than 30 ° otherwise can only see very little crescent moon, the embodiment threshold value is selected 30 °.
(2) obtain to satisfy the earth under the orbital arc fragment position in the step (1)/moon heart vector set according to the Satellite Attitude Movement of design in advance, whether in described vector set have vector be present in ultraviolet sensors field of view in, select to satisfy the track segmental arc in the ultraviolet sensors field of view if judging;
Whether in described vector set have vector be present in ultraviolet sensors field of view in, method is for judging vector and visual field center vector F if judging sAngle  whether less than field angle, select the orbital position of the vector correspondence in the ultraviolet sensors field of view;
(3) orbital position that will satisfy step (1) and step (2) simultaneously before described optional track arrives, is imported the predicted orbit position of described optional track correspondence, optional track a certain moment t before as taking optional track after the satellites transmits 0Attitude matrix, life ground ephemeris and Satellite Attitude Movement model; Calculate attitude matrix, month heart/the earth's core vector in the predicted time section, and make a moon heart/the earth's core vector and the time dependent curve of visual field center vector angle, obtain the visible time range [t of the continuous earth/moon by curve and ultraviolet sensors visual field size 1~t 2], according to the continuously visible time range calculating optimum that obtains constantly.
Moon heart vector V under the body series M:
x_moon=-0.6802; y_moon=0.3213; z_moon=-0.6589;
Solar vector V S:
x_sun=1; y_sun=0; z_sun=0;
A solar vector and a month heart vector angle β:
β=acos(V M·V S)=132.86°
As seen angle is much larger than 30 °, is between crescent or half moon and the full moon to show as gibbous moon.
According to t 0Attitude matrix C 0Calculate predicted time interval [t 0~t n] corresponding attitude matrix C t:
C t=R t·C 0
R wherein tBeing the posture changing matrix, is to be determined by the Satellite Attitude Movement model, sets satellite in the utilization example to center on+the X-axis rotation, and speed is that 0.1 °/s is designated as V, selects predicted time length t L=t n-t 0Be 1 hour, remember t among the embodiment n=3600s, t 0=Os, R tFor:
R t = 1 0 0 0 cos θ sin θ 0 - sin θ cos θ
Rotate then around+Y-axis if select:
R t = cos θ 0 - sin θ 0 1 0 sin θ 0 cos θ
Rotate then around+Z axle if select:
R t = cos θ sin θ 0 - sin θ cos θ 0 0 0 1
In the formula: θ=(t-t 0) * V
According to posture changing matrix R tObtain the moon heart vector V under the moment t MtWith the earth's core vector V Et
M t=R t·M 0
E t=R t·E 0
Calculate a month heart vector M tThe earth's core vector E tAngle  with the centre normal Fs of visual field
M=acos(Fs·M t)
E=acos(Fs·E t)
M<T FOVThe moon as seen, and is otherwise invisible
E<T FOVThe earth as seen, and is otherwise invisible
T FOVVisual field size for visual field threshold value Here it is sensor optical system is set at 10 ° in the example.The Fs value is as follows among the embodiment:
x_Fs=0 y_Fs=0 z_Fs=1
Predicted time is selected 1 hour, and the angle curve at month heart vector and center, visual field is seen Fig. 2 in the predicted time that calculates so, and Fig. 3 is seen less than 10 ° curve in the visual field continuously in the search back, obtains the time:
t 1=2510(s) t 2=2760(s)
The best time t then Optimal:
Have according to time average:
t optimal = t 2 + t 1 2 t optimal=2635(s)
[t 1~t 2] minimum angle  is constantly corresponding for time range:
t optimal={t| t=min()} t optimal=2640(s)
As seen predicted time is very approximate under two standards, and the corresponding analogous diagram of taking of time average is seen Fig. 4.
Inventive concept and method can be applied in the time prediction of other space sensitive device, as long as the method for utilizing thought of the present invention to realize all falls into protection scope of the present invention, the unspecified part of the present invention belongs to general knowledge as well known to those skilled in the art.

Claims (6)

1, moon-tracking segment ultraviolet sensor earth and moon photography time predication method is characterized in that comprising the following steps:
(1) carries out the illumination analysis according to a solar vector and a month heart vector in nominal track, the satellite body system, select the track segmental arc of taking illumination condition by the moon contentedly;
(2) obtain to satisfy the earth under the orbital arc fragment position in the step (1)/moon heart vector set according to the Satellite Attitude Movement of design in advance, judge whether described vector set is present in the ultraviolet sensors field of view, select to satisfy the track segmental arc in the ultraviolet sensors field of view;
(3) will satisfy the orbital position of step (1) and step (2) simultaneously as taking optional track, before described optional track arrives, import the predicted orbit position of described optional track correspondence, optional track a certain moment attitude matrix, life ground ephemeris and Satellite Attitude Movement model before after the satellites transmits; Calculate attitude matrix, month heart/the earth's core vector in the predicted time section, and make a moon heart/the earth's core vector and the time dependent curve of visual field center vector angle, obtain the visible time range [t of the continuous earth/moon by curve and ultraviolet sensors visual field size 1~t 2], according to the continuously visible time range calculating optimum that obtains constantly.
2, moon-tracking segment ultraviolet sensor earth and moon photography time predication method according to claim 1 is characterized in that: the illumination analytic process in the described step (1) is: calculate solar vector V SWith moon heart vector V MOr the earth's core vector V EBetween angle , computing formula is:
=acos(V S·V i)
As > T, think under this orbital position that contentedly month taking an illumination condition obtains the better phases of the moon;
Wherein, V iRepresent a month heart vector V MOr the earth's core vector V E
TBe the angle threshold value, size is determined according to the shooting demand.
3, moon-tracking segment ultraviolet sensor earth and moon photography time predication method according to claim 1 is characterized in that: the computing formula of attitude matrix is in the described step (3):
C t=R t·C 0
Wherein, R tBe the posture changing matrix, determine according to the Satellite Attitude Movement model;
C 0A certain moment attitude matrix for input.
4, moon-tracking segment ultraviolet sensor earth and moon photography time predication method according to claim 1 is characterized in that: calculating month heart/the earth's core vector formula in the described step (3) is:
M t=R t·M 0
E t=R t·E 0
Wherein, M t-t moon heart vector constantly;
E t-t the earth's core vector constantly;
R tBe the posture changing matrix, determine according to the Satellite Attitude Movement model;
M 0The a certain moment t of-input 0Corresponding moon heart vector.
5, moon-tracking segment ultraviolet sensor earth and moon photography time predication method according to claim 1 is characterized in that: visible time range deterministic process is as follows in the described step (3):
At first, calculate t moon heart vector M constantly tWith the earth's core vector E tWith the angle  of the center vector Fs of ultraviolet sensors visual field, computing formula is:
M(t)=acos(Fs·M t)
E(t)=acos(Fs·E t)
Then, with the value and the visual field threshold value T of curve correspondence FOVCompare, work as  M<T FOV, the moon is described as seen, otherwise invisible; Work as  E<T FOV, the earth is described as seen, otherwise invisible; Continuous visible time border is visible time range; Wherein, T FOVEqual the visual field size of ultraviolet sensors optical system for the visual field threshold value.
6, moon-tracking segment ultraviolet sensor earth and moon photography time predication method according to claim 1 is characterized in that: the best time is got described visible time range [t in the described step (3) 1~t 2] intermediate value, or [t 1~t 2] the angle  moment corresponding of time range minimum.
CN200810057347A 2008-01-31 2008-01-31 The moon-tracking segment ultraviolet sensor earth and moon photography time predication method Active CN100585334C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200810057347A CN100585334C (en) 2008-01-31 2008-01-31 The moon-tracking segment ultraviolet sensor earth and moon photography time predication method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200810057347A CN100585334C (en) 2008-01-31 2008-01-31 The moon-tracking segment ultraviolet sensor earth and moon photography time predication method

Publications (2)

Publication Number Publication Date
CN101236090A true CN101236090A (en) 2008-08-06
CN100585334C CN100585334C (en) 2010-01-27

Family

ID=39919870

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200810057347A Active CN100585334C (en) 2008-01-31 2008-01-31 The moon-tracking segment ultraviolet sensor earth and moon photography time predication method

Country Status (1)

Country Link
CN (1) CN100585334C (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102175258A (en) * 2010-12-31 2011-09-07 北京控制工程研究所 Correction method of moon center direction of moon sensor
CN104567881A (en) * 2014-12-26 2015-04-29 北京控制工程研究所 Patrol device position posture determining method based on gravity vectors of sun, earth core and lunar surface
CN109724590A (en) * 2019-01-04 2019-05-07 中国科学院遥感与数字地球研究所 A kind of calculation method of month the based platform optical sensor whole world video imaging frequency and time interval
CN111366149A (en) * 2020-03-25 2020-07-03 中国科学院微小卫星创新研究院 Method for extracting sun vector direction based on lunar phase

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102175258A (en) * 2010-12-31 2011-09-07 北京控制工程研究所 Correction method of moon center direction of moon sensor
CN102175258B (en) * 2010-12-31 2012-08-22 北京控制工程研究所 Correction method of moon center direction of moon sensor
CN104567881A (en) * 2014-12-26 2015-04-29 北京控制工程研究所 Patrol device position posture determining method based on gravity vectors of sun, earth core and lunar surface
CN109724590A (en) * 2019-01-04 2019-05-07 中国科学院遥感与数字地球研究所 A kind of calculation method of month the based platform optical sensor whole world video imaging frequency and time interval
CN111366149A (en) * 2020-03-25 2020-07-03 中国科学院微小卫星创新研究院 Method for extracting sun vector direction based on lunar phase
CN111366149B (en) * 2020-03-25 2022-02-15 中国科学院微小卫星创新研究院 Method for extracting sun vector direction based on lunar phase

Also Published As

Publication number Publication date
CN100585334C (en) 2010-01-27

Similar Documents

Publication Publication Date Title
WO2013070518A1 (en) Solar timer using gps technology
CN100533065C (en) Interstellar cruising self-navigation method based on multi-star road sign
Buie et al. Size and shape from stellar occultation observations of the double Jupiter Trojan Patroclus and Menoetius
Miyake et al. A sub-Saturn mass planet, MOA-2009-BLG-319Lb
CN100585334C (en) The moon-tracking segment ultraviolet sensor earth and moon photography time predication method
CN103837150A (en) Method for performing rapid celestial fix through CCD (charge coupled device) zenith telescope on ground
Santos-Sanz et al. James Webb Space Telescope observations of stellar occultations by solar system bodies and rings
Flohrer et al. Proposed strategies for optical observations in a future European Space Surveillance network
Fujita et al. An Overview of JAXA's Ground-Observation Activities for HAYABUSA Reentry
Helin et al. Search techniques for near-Earth asteroids
Freeman et al. Can the masses of isolated planetary-mass gravitational lenses be measured by terrestrial parallax?
Weaver et al. High-resolution Search for Kuiper Belt Object Binaries from New Horizons
Sydney et al. Raven automated small telescope systems
Liebe Solar compass chip
Plotke et al. Dual use star tracker and space domain awareness sensor in-space test
Naidu et al. Orbital and physical characterization of asteroid dimorphos following the DART impact
Leonard et al. UMaMI: a new frontiers-style mission concept to explore the Uranian system
Nelson et al. Navigation and Orbit Estimation for New Horizons’ Arrokoth Flyby: Overview, Results and Lessons Learned
Sekanina Unprecedented Daylight Display of Kreutz Sungrazers in AD 363?
Betto et al. Advanced stellar compass deep space navigation, ground testing results
Dombrowski Investigation of Mirror Satellite Concept to Augment Natural Lighting Conditions on Orbit
Umehara et al. Scan by Monitoring a Pair of Points Optical Survey Method for Near-Geosynchronous Orbits
RU2576643C2 (en) Ultrahigh altitude sun synchronous orbit satellite system
Buchheim et al. The Planets
Vaughan et al. MESSENGER Spacecraft Pointing Performance during the Mission’s Mercury Orbital Phase

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