CN107421504B - The shooting time calculation method of month base earth observation electro-optical photo - Google Patents
The shooting time calculation method of month base earth observation electro-optical photo Download PDFInfo
- Publication number
- CN107421504B CN107421504B CN201710650254.1A CN201710650254A CN107421504B CN 107421504 B CN107421504 B CN 107421504B CN 201710650254 A CN201710650254 A CN 201710650254A CN 107421504 B CN107421504 B CN 107421504B
- Authority
- CN
- China
- Prior art keywords
- moon
- observed
- target
- earth
- shooting time
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electric Clocks (AREA)
Abstract
The invention discloses the shooting time calculation methods of a kind of month base earth observation electro-optical photo, it is calculated for shooting time of the photoelectric sensor laid on moon during being observed to tellurian target, this method comprises: the location information based on target being observed on position on the moon of ephemeris, moon based sensor and the earth, calculate under same referential the moon based sensor, the time-varying position of target being observed and the sun on the earth；The light conditions of target are obtained based on the position of target being observed and the sun；Position based on position of the earth, target being observed and sensor obtains the observed altitude angle of sensor；If target being observed meets solar irradiation and observed altitude corner condition simultaneously, this moment is can shooting time；The moon base observation shooting time that can be obtained in the period is calculated by each moment of the period to selection.Computational efficiency of the present invention is high, visual result is illustrated, and has practical value in moon base earth observation field.
Description
Technical field
The invention belongs to empty day earth observation message areas, relate to the use of the photoelectric sensor laid on the moon to the earth
On target its shooting time when being observed calculation method.
Background technique
It ascends a height and enjoys a distant view from initial to utilization balloon, dirigible, aircraft, then arrives satellite, space shuttle and space station, see over the ground
The platform and sensor of survey are all in development at full speed.Phase of the moon as the unique natural satellite of the earth, with the sun and the earth
Mutual relation makes it have unique advantage in terms of earth observation.The earth observation potentiality of month base earth observation sensor will be complete
Ball scale multi-turn layer Integrated research provides safeguard, it is possible to really realize the angle from earth system science to global multi-turn layer phase
A series of problem in science of keys of mutual coupling provides completely new answer, including global ocean-atmosphere interaction, Land-air interaction
It is existing with boundary layer Atmospheric processes, sea Lu Guanlian variation and littoral zone process etc., especially solid tide, global energy balance etc.
As.
Shooting time during being observed using the photoelectric sensor laid on the moon to tellurian target
Selection mainly with day the moon three's positional relationship and the sensor position and target being observed laid on the moon on earth
Position it is related.1972, spacefarer utilized Extreme Ultraviolet Telescope on the moon when No. 16 airship manned moon landings of american apollo
The earth is shot, obtains the far ultraviolet photo of the earth and being shown on the net in Nasa, but do not have
Illustrate the calculation method of shooting time selection；Optical precision engineering [J], the 12nd phase of volume 18 in 2010, " moon base observes earth etc.
Observation time of the extreme ultraviolet observation instrument in a lunar cycles is mentioned in gas ions layer extreme ultraviolet radiation characteristic " to need to meet
Condition and carried out preliminary analysis, but do not provide the calculation method of specific shooting time selection.
It is the weight for carrying out the earth observation of moon base that the photoelectric sensor laid on the moon, which is observed tellurian target,
Component part is wanted, the selection to shooting time is the important content of research required for playing moon base earth observation advantage, therefore, such as
What, which finds a kind of accurate and high-precision shooting time selection method, is necessary.
Summary of the invention
The present invention solves the technical problem of be directed to the photoelectric sensor laid on the moon to tellurian target
During being observed, the problem of how selecting the time shot.
The present invention is realized by following steps:
Step 1: the moon based sensor and target being observed position resolving
(1) the initial time t calculated as needed1With end time t2, first by moon based sensor longitude and latitude positionIt transforms to and consolidates three-dimensional coordinate by moon heart moonThen will consolidate three-dimensional coordinate transformation the heart moon moon to arrive
Month heart Celestial Reference System coordinate, further by moon heart Celestial Reference System coordinate transform to the earth's core Celestial Reference System coordinate, further according to
International Astronomical Union (IAU) (IAU) resolution and International Earth Rotation service (IERS) specification will be consolidated three-dimensional coordinate transformation the heart moon moon and arrived
Ground heart is consolidated under referential, and above-mentioned transformation relation can be calculate by the following formula:
Wherein [W] is Ghandler motion matrix, and [R] is earth rotation matrix, and [N] is nutating matrix, and [P] is precession of the equinoxes matrix, and [B] is
Frame bias matrix, these matrixes can be resolved according to International Astronomical Union (IAU) (IAU) and International Earth Rotation service
(IERS) proper calculation obtains；For the trivector of the moon heart and the earth's core, [L] is lunar libration spin matrix, can pass through
Almanac data, which calculates, to be obtained；
(2) the initial time t calculated as needed1With end time t2, by the longitude and latitude position of target being observed on the earthIt transforms to ground heart and consolidates target being observed coordinate under referential
(3) the initial time t calculated as needed1With end time t2, the almanac data based on U.S. jet laboratory with
And the coordinate position of International Astronomical Union (IAU) (IAU) resolution and International Earth Rotation service (IERS) proper calculation sun
Step 2: the solar irradiation condition of target being observed is calculated
Assuming that the earth is sphere, sunlight is incident in parallel, consolidates the coordinate position under referential in ground heart based on the sunCalculate substar positionThen turn
It is changed to latitude and longitude coordinates positionBy rotation transformation by latitude and longitude coordinates position
(φt；λt) it is transformed to pole P, the longitude and latitude position of target being observed at this timeIt can be with the difference of latitude of pole P
It indicates are as follows:
φT=arccos (cos φt cosφcosλt cosλ+cosφt cosφsinλtsinλ+sinφt sinφ)
(2)
Transformed difference of latitude is then set less than 90 degreeIt is 1, is otherwise set as 0；
Step 3: the observed altitude angle to target of sensor is calculated
Coordinate position based on target being observedAnd sensor coordinatesMeter
Calculate the observed altitude angle θ (t) to target at the earth's surface of sensor, the calculation formula are as follows:
Step 4: if target being observed meets solar irradiation simultaneouslyCondition and observed altitude angle θ (t) item
Part, this moment is can shooting time.It calculates and can obtain by each moment of the starting and ending time to selection
Moon base earth observation shooting time in this period.Also it can choose the smaller restriction range in observed altitude angle.
The beneficial effects of the present invention are:
(1) the method for the present invention is based on almanac data, while considering to be seen on moon position of based sensor laying and the earth
With the surveying target position and day correlation between the moon, it is established that the moon based sensor sat admittedly from lunar surface latitude and longitude coordinates to ground heart
The calculation formula of mark system conversion, can accurately predict booting shooting time.
(2) the present invention is based on the observation conditions and sensor observed altitude corner condition that calculate the target being observed obtained, come
Continuous shooting time section is obtained, has the advantages that calculating is simple, complete, can satisfy the base earth observation of the photoelectric sensor moon
Needs.
Detailed description of the invention
The above is merely an overview of the technical solutions of the present invention, in order to better understand the technical means of the present invention, below
In conjunction with attached drawing, the present invention is described in further detail with specific embodiment.
Fig. 1 is the techniqueflow chart of the method for the present invention；
Fig. 2 is three-dimensional position schematic diagram of month based sensor under ECEF coordinate system；
Fig. 3 is the difference of latitude schematic diagram between the position of target being observed and pole P；
Fig. 4 is a month elevation angle schematic diagram for based sensor observation；
Fig. 5 be photoelectric sensor can shooting time schematic diagram.
Specific embodiment
Below by 1 day 12 January of 2018 Coordinated Universal Time(UTC) for 31 Time of Day section in January, in conjunction with attached drawing to the present invention
Be described in detail, so that the advantages and features of the present invention can be more easily understood by those skilled in the art, thus to this hair
Bright protection scope is more clearly defined.
Refering to fig. 1, the invention discloses the shooting time calculation methods of a kind of month base earth observation electro-optical photo, including with
Lower step:
Step 1: the resolving of moon based sensor and target being observed position
Moon based sensor is laid in No. three landing positions of the goddess in the moon, 19.51 degree of moonscape west longitude, 44.12 degree of north latitude
Region, target being observed are located at 60 degree of earth surface east longitude, 30 degree of north latitude.
(1) the initial time Coordinated Universal Time(UTC) t calculated as needed1On January 1st, 1 and end time t2=2018
On January 31, in, first by moon based sensor longitude and latitude positionTransform to moon heart moon solid three-dimensional
CoordinateUnit is rice, then will consolidate three-dimensional coordinate transformation the heart moon moon and arrives
Month heart Celestial Reference System coordinate, further by moon heart Celestial Reference System coordinate transform to the earth's core Celestial Reference System coordinate, further according to
International Astronomical Union (IAU) (IAU) resolution and International Earth Rotation service (IERS) specification will be consolidated three-dimensional coordinate transformation the heart moon moon and arrived
Ground heart is consolidated under referential, such as t1Moment sensor coordinates areUnit is
Kilometer.Fig. 2 is that three-dimensional position diagram of month based sensor under ECEF coordinate system is intended to；
(2) by the longitude and latitude position of target being observed on the earthIt transforms to ground heart and consolidates referential
Lower target being observed coordinate isUnit is rice；
(3) almanac data and International Astronomical Union (IAU) (IAU) resolution and international ground based on U.S. jet laboratory
The coordinate position of revolutions service (IERS) proper calculation sun, such as t1Moment isUnit is kilometer；
Step 2: the solar irradiation condition of target being observed is calculated
Assuming that the earth is sphere, sunlight is incident in parallel, consolidates the coordinate position under referential in ground heart based on the sunCalculate substar positionThen
It is converted into latitude coordinate locationBy rotation transformation by latitude and longitude coordinates position
Set (φt；λt) it is transformed to pole P, the longitude and latitude position of target being observed at this timeDifference of latitude with pole P can
To indicate are as follows:
φT=arccos (cos φt cosφcosλt cosλ+cosφt cosφsinλtsinλ+sinφt sinφ)
(1)
Transformed difference of latitude is then set less than 90 degreeIt is 1, is otherwise set as 0.Such as t1The difference of latitude at moment is
76.7 degree, ifIt is 1.Fig. 3 is the difference of latitude schematic diagram between the position of target being observed and pole P；
Step 3: the observed altitude angle to target of sensor is calculated
Coordinate position based on target being observedAnd sensor coordinatesMeter
Calculate the observed altitude angle θ (t) to target at the earth's surface of sensor, the calculation formula are as follows:
Such as t1Moment elevation angle θ (t1) it is 99.3 degree, Fig. 4 is a month based sensor observed altitude angle schematic diagram；
Step 4: if target being observed meets solar irradiation simultaneouslyCondition and observed altitude angle θ (t) are less than
90 degree, this moment is can shooting time.It calculates and can obtain by each moment of the starting and ending time to selection
Obtain the moon base earth observation shooting time in this period.Such as t1Moment due to observed altitude angle be greater than 90 degree, so
Effective shooting is not can be carried out.Fig. 5 be photoelectric sensor can shooting time schematic diagram.
It is sub that the above description is only an embodiment of the present invention, is not intended to limit the scope of the invention, all to utilize this
Equivalent structure or equivalent flow shift made by description of the invention and accompanying drawing content, it is relevant to be applied directly or indirectly in other
Technical field is included within the scope of the present invention.
Claims (4)
1. the shooting time calculation method of a kind of month base earth observation electro-optical photo, it is characterised in that: for what is laid on the moon
Photoelectric sensor is calculated to tellurian target progress shooting time, and specific step is as follows for this method:
Step 1: calculative initial time t is selected1With end time t2, based on almanac data, moon based sensor in the moon
The longitude and latitude position of upper layingAnd on the earth target being observed longitude and latitude positionIt is logical
Coordinate transform is crossed, heart consolidates the sensor coordinates under referential with obtaining each momentIt is observed mesh
Mark coordinateAnd the coordinate position of the sun
Step 2: assuming that the earth is sphere, the sunlight incident earth in parallel, the coordinate position based on the sun
And the coordinate position of target being observedThe light conditions of target being observed are calculated, if being observed mesh
Mark is under solar irradiation, ifIt is 1, is otherwise set as 0；
Step 3: the coordinate position based on target being observedAnd sensor coordinatesCalculate the observed altitude angle θ (l) to target at the earth's surface of sensor, the calculation formula
Are as follows:
Step 4: if target being observed meets solar irradiation simultaneouslyCondition and observed altitude angle θ (t) condition are for example
Less than 90 degree, this moment is can shooting time；Being calculated by each moment of the starting and ending time to selection can
To obtain the moon base earth observation shooting time in this period.
2. the shooting time calculation method of moon base earth observation electro-optical photo according to claim 1, it is characterised in that: described
The longitude and latitude position that moon based sensor in step 1 is laid on the moonConsolidate under referential to ground heart
Sensor coordinatesConversion process in needed first by moon based sensor longitude and latitude position
It transforms to and consolidates three-dimensional coordinate by moon heart moonThen will consolidate three-dimensional coordinate transformation the heart moon moon to refer to moon heart celestial sphere
It is coordinate, further by moon heart Celestial Reference System coordinate transform to the earth's core Celestial Reference System coordinate, further according to International Astronomical student's federation
Credit union (IAU) resolution and International Earth Rotation service (IERS) specification will be consolidated three-dimensional coordinate transformation the heart moon moon and referred to admittedly to ground heart
Under system, above-mentioned transformation relation can be calculate by the following formula:
Wherein [W] is Ghandler motion matrix, and [R] is earth rotation matrix, and [N] is nutating matrix, and [P] is precession of the equinoxes matrix, and [B] is frame
Bias matrix, these matrixes can be resolved according to International Astronomical Union (IAU) (IAU) and International Earth Rotation service (IERS)
Proper calculation obtains；For the trivector of the moon heart and the earth's core, [L] is lunar libration spin matrix, can pass through ephemeris number
It is obtained according to calculating.
3. the shooting time calculation method of moon base earth observation electro-optical photo according to claim 1, it is characterised in that: described
The process that the light conditions of target being observed are calculated in step 2 is as follows:
It is primarily based on the sun and consolidates the coordinate position under referential in ground heartCalculate its corresponding substar position
It setsIt is then converted into latitude and longitude coordinates positionBy rotation transformation by coordinate position (φt, λt) it is transformed to pole P, at this time
The longitude and latitude position of target being observedIt will be converted with the difference of latitude of pole P are as follows:
φT=arccos (cos φtcosφcosλtcosλ+cosφtcosφsinλtsinλ+sinφtsinφ) (3)
Transformed difference of latitude is less than such as 90 degree of certain angle, then setsIt is 1, is otherwise set as 0.
4. the shooting time calculation method of moon base earth observation electro-optical photo according to claim 1, it is characterised in that: described
In step 4 can shooting time be between the starting and ending time of selection, and meet solar irradiation simultaneouslyItem
The period of part and observed altitude angle θ (t) condition.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710650254.1A CN107421504B (en) | 2017-08-02 | 2017-08-02 | The shooting time calculation method of month base earth observation electro-optical photo |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710650254.1A CN107421504B (en) | 2017-08-02 | 2017-08-02 | The shooting time calculation method of month base earth observation electro-optical photo |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107421504A CN107421504A (en) | 2017-12-01 |
| CN107421504B true CN107421504B (en) | 2019-11-12 |
Family
ID=60436573
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710650254.1A Active CN107421504B (en) | 2017-08-02 | 2017-08-02 | The shooting time calculation method of month base earth observation electro-optical photo |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107421504B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108896038B (en) * | 2018-05-15 | 2022-03-25 | 中国科学院遥感与数字地球研究所 | Imaging method of lunar-based optical sensor |
| CN108931259B (en) * | 2018-05-25 | 2021-02-05 | 中国科学院遥感与数字地球研究所 | Method for simulating lunar-base earth observation image |
| CN109724590B (en) * | 2019-01-04 | 2020-07-24 | 中国科学院遥感与数字地球研究所 | Method for calculating global image imaging frequency and time interval of lunar-based platform optical sensor |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101424530A (en) * | 2008-12-09 | 2009-05-06 | 武汉大学 | Method for generating approximate kernel line of satellite stereo image pairs based on projection reference surface |
| CN104931022A (en) * | 2015-04-21 | 2015-09-23 | 国家测绘地理信息局卫星测绘应用中心 | Satellite image three-dimensional area network adjustment method based on satellite-borne laser height measurement data |
| CN105550639A (en) * | 2015-12-07 | 2016-05-04 | 国家测绘地理信息局卫星测绘应用中心 | Automatic extraction method for Earth observation laser height measurement satellite elevation control points and data processing method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2935579B1 (en) * | 2008-08-28 | 2010-11-05 | Centre Nat Etd Spatiales | METHOD OF ACQUIRING, REDUCING AND TRANSMITTING SATELLITE IMAGES |
-
2017
- 2017-08-02 CN CN201710650254.1A patent/CN107421504B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101424530A (en) * | 2008-12-09 | 2009-05-06 | 武汉大学 | Method for generating approximate kernel line of satellite stereo image pairs based on projection reference surface |
| CN104931022A (en) * | 2015-04-21 | 2015-09-23 | 国家测绘地理信息局卫星测绘应用中心 | Satellite image three-dimensional area network adjustment method based on satellite-borne laser height measurement data |
| CN105550639A (en) * | 2015-12-07 | 2016-05-04 | 国家测绘地理信息局卫星测绘应用中心 | Automatic extraction method for Earth observation laser height measurement satellite elevation control points and data processing method |
Non-Patent Citations (1)
| Title |
|---|
| 面向全球变化探测的月基成像雷达概念研究;郭华东，丁翼星，刘广，张道卫，傅文学，张露;《中国科学:地球科学》;20131030;第43卷(第11期);1760-1769 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107421504A (en) | 2017-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9551980B2 (en) | Solar timer using GPS technology | |
| Zhang | Star identification | |
| CN103047985B (en) | A kind of method for rapidly positioning of extraterrestrial target | |
| CN101893440B (en) | Celestial autonomous navigation method based on star sensors | |
| CN102175241B (en) | Autonomous astronomical navigation method of Mars probe in cruise section | |
| CN107421504B (en) | The shooting time calculation method of month base earth observation electro-optical photo | |
| CN102706363B (en) | Precision measuring method of high-precision star sensor | |
| CN104462776A (en) | Method for absolutely radiometric calibration of low orbit earth observation satellite with moon as reference | |
| CN103063217A (en) | Deep space detector astronomy/radio combination navigation method based on ephemeris correction | |
| CN108317993A (en) | A kind of deviation of plumb line measuring device and method of integrated GNSS and laser tracker | |
| CN103837150A (en) | Method for performing rapid celestial fix through CCD (charge coupled device) zenith telescope on ground | |
| CN102288201A (en) | Precision measurement method for star sensor | |
| CN102288200A (en) | Accuracy measurement system for star sensor | |
| CN105446346A (en) | Moon relative calibrating posture adjustment method by remote sensing satellite | |
| Wang et al. | Image and data processing of digital zenith telescope (DZT-1) of China | |
| CN105486315A (en) | Method for adjusting remote sensing satellite absolute calibration attitude relative to moon | |
| Stooke | The international atlas of Mars exploration: volume 1, 1953 to 2003: the first five decades | |
| CN110146052A (en) | A kind of plane normal astronomical orientation measurement method and system based on total station | |
| Zeng et al. | Pan-location mapping and localization for the in-situ science exploration of Zhurong Mars rover | |
| CN108896038A (en) | Month base optical sensor imaging method | |
| CN102607597B (en) | Three-axis precision expression and measurement method for star sensor | |
| CN104458653B (en) | Method and system for measuring atmospheric refraction value at large zenith distance | |
| Liebe | Solar compass chip | |
| Shaltout et al. | Studying the orientations of luxor ancient Egyptian temples using QuickBird images | |
| CN108931259B (en) | Method for simulating lunar-base earth observation image |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |