CN107421504A - The shooting time computational methods of month base earth observation electro-optical photo - Google Patents
The shooting time computational methods of month base earth observation electro-optical photo Download PDFInfo
- Publication number
- CN107421504A CN107421504A CN201710650254.1A CN201710650254A CN107421504A CN 107421504 A CN107421504 A CN 107421504A CN 201710650254 A CN201710650254 A CN 201710650254A CN 107421504 A CN107421504 A CN 107421504A
- 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.)
- Granted
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 computational methods of a kind of month base earth observation electro-optical photo, calculated for the shooting time of the photoelectric sensor laid on the moon during being observed to tellurian target, this method includes：Based on the positional information of 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 understands there is practical value in moon base earth observation field.
Description
Technical field
The invention belongs to empty day earth observation message area, relates to the use of the photoelectric sensor laid on the moon to the earth
On target its shooting time when being observed computational methods.
Background technology
Ascend a height and enjoy a distant view to utilization balloon, dirigible, aircraft from initial, then to satellite, space shuttle and space station, see over the ground
The platform and sensor of survey are all in development at full speed.The moon is as the unique natural satellite of the earth, its phase 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
The more ring layer Integrated research of ball yardstick provide safeguard, it is possible to really realize the angle from earth system science to the more ring layer phases in the whole world
A series of crucial problem in science of mutual coupling provide brand-new answer, including global ocean-atmosphere interaction, Land-air interaction
It is existing with boundary layer Atmospheric processes, extra large Lu Guanlian changes and littoral zone process etc., particularly 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 position relationship and the sensing station and target being observed laid on the moon on earth
Position it is relevant.1972, spacefarer utilized Extreme Ultraviolet Telescope on the moon during No. 16 airship manned moon landings of american apollo
The earth is shot, the far ultraviolet photo of the earth is obtained and is shown in the online of Nasa, but do not have
Illustrate the computational methods of shooting time selection；Optical precision engineering [J], the 12nd phase of volume 18 in 2010, " moon base observation 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 be needed to meet
Condition and carried out initial analysis, but do not provide the computational methods of specific shooting time selection.
It is the weight for carrying out moon base earth observation that the photoelectric sensor laid on the moon is observed to tellurian target
Part is wanted, the selection to shooting time is to play the important content studied required for moon base earth observation advantage, therefore, such as
It is necessary that what, which finds a kind of accurate and high-precision shooting time system of selection,.
The content of the invention
The present invention solves the technical problem of the photoelectric sensor for being 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：Month based sensor and the resolving of target being observed position
(1) the initial time t calculated as needed1With end time t2, first by moon based sensor longitude and latitude positionTransform to moon heart moon solid three-dimensional coordinateThen the heart moon moon solid three-dimensional coordinate transformation is arrived
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) is resolved and the heart moon moon is consolidated three-dimensional coordinate transformation by International Earth Rotation service (IERS) specification
To under the solid referential of ground heart, above-mentioned transformation relation can be calculated by 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, [B]
For framework 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 moon heart and the trivector in the earth's core, [L] is lunar libration spin matrix, can passed 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 earthTransform to target being observed coordinate under the solid referential of ground heart
(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：Calculate the solar irradiation condition of target being observed
Assuming that the earth is spheroid, sun parallel light is incident, the coordinate position based on the sun under the solid referential of ground heartCalculate substar positionThen turn
It is changed to latitude and longitude coordinates positionBy rotation transformation by latitude and longitude coordinates position
(φt；λt) limit P is transformed to, the now longitude and latitude position of target being observedCan be with limit P difference of latitude
It is expressed as：
φT=arccos (cos φt cosφcosλt cosλ+cosφt cosφsinλtsinλ+sinφt sinφ)
(2)
Difference of latitude after conversion is less than 90 degree, then setsFor 1,0 is otherwise set to；
Step 3：Calculate the observed altitude angle to target of sensor
Coordinate position based on target being observedAnd sensor coordinatesMeter
The observed altitude angle θ (t) to target at the earth's surface of sensor is calculated, described calculation formula is：
Step 4：If target being observed meets solar irradiation simultaneouslyCondition and observed altitude angle θ (t) bars
Part, this moment is can shooting time.Carrying out calculating by each moment of the starting and ending time to selection can obtain
Moon base earth observation shooting time in this period.The restriction scope that observed altitude angle is smaller can also be selected.
The beneficial effects of the invention are as follows：
(1) the inventive method is based on almanac data, while considers to be seen on moon position of based sensor laying and the earth
With the surveying target location 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 start shooting time.
(2) present invention is come based on the observation condition and sensor observed altitude corner condition for calculating the target being observed obtained
Continuous shooting time section is obtained, has and calculates the advantages of simple, complete, the base earth observation of the photoelectric sensor moon can be met
Needs.
Brief description of the drawings
Above-mentioned is only the general introduction of technical solution of the present invention, in order to better understand the technological means of the present invention, below
With reference to accompanying drawing, the present invention is described in further detail with embodiment.
Fig. 1 is the techniqueflow chart of the inventive method；
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 limit P；
Fig. 4 is a month elevation angle schematic diagram for based sensor observation；
Fig. 5 be photoelectric sensor can shooting time schematic diagram.
Embodiment
Below by 1 day 12 January of 2018 Coordinated Universal Time(UTC) exemplified by January 31 Time of Day section, with reference to accompanying drawing to the present invention
Be described in detail so that advantages and features of the invention can be easier to be readily appreciated by one skilled in the art, so as to this hair
Bright protection domain, which is made, apparent clearly to be defined.
Refering to Fig. 1, the invention discloses the shooting time computational 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 arrives the heart moon moon solid three-dimensional coordinate transformation
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) resolves and International Earth Rotation service (IERS) specification arrives the heart moon moon solid three-dimensional coordinate transformation
Under the solid referential of ground heart, 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 earthTransform to the solid referential of ground heart
Lower target being observed coordinate isUnit is rice；
(3) almanac data and International Astronomical Union (IAU) (IAU) resolution based on U.S. jet laboratory
With the coordinate position of International Earth Rotation service (IERS) proper calculation sun, such as t1Moment isUnit is kilometer；
Step 2：Calculate the solar irradiation condition of target being observed
Assuming that the earth is spheroid, sun parallel light is incident, the coordinate position based on the sun under the solid referential of ground heartCalculate substar positionThen turn
It is changed to latitude and longitude coordinates positionBy rotation transformation by latitude and longitude coordinates position
(φt；λt) limit P is transformed to, the now longitude and latitude position of target being observedCan be with limit P difference of latitude
It is expressed as：
φT=arccos (cos φt cosφcosλt cosλ+cosφt cosφsinλtsinλ+sinφt sinφ)
(1)
Difference of latitude after conversion is less than 90 degree, then setsFor 1,0 is otherwise set to.Such as t1The difference of latitude at moment is
76.7 degree, ifFor 1.Fig. 3 is the difference of latitude schematic diagram between the position of target being observed and limit P；
Step 3：Calculate the observed altitude angle to target of sensor
Coordinate position based on target being observedAnd sensor coordinatesMeter
The observed altitude angle θ (t) to target at the earth's surface of sensor is calculated, described calculation formula is：
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 small
In 90 degree, this moment is can shooting time.Carrying out calculating by each moment of the starting and ending time to selection can be with
Obtain the moon base earth observation shooting time in this period.Such as t1Moment due to observed altitude angle be more than 90 degree, institute
Can not effectively be shot.Fig. 5 be photoelectric sensor can shooting time schematic diagram.
Embodiments of the invention is the foregoing is only, is not intended to limit the scope of the invention, it is every to utilize this
The equivalent structure or equivalent flow conversion that description of the invention and accompanying drawing content are made, or directly or indirectly it is used in other correlations
Technical field, it is included within the scope of the present invention.
Claims (4)
1. the shooting time computational methods 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 calculates to tellurian target progress shooting time, and this method comprises the following steps that：
Step 1：Select calculative initial time t1With 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, with obtaining each moment the sensor coordinates under the solid referential of heartIt is observed mesh
Mark coordinateAnd the coordinate position of the sun
Step 2：Assuming that the earth is spheroid, and the sun parallel light incidence earth, 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, ifFor 1,0 is otherwise set to；
Step 3：Coordinate position based on target being observedAnd sensor coordinatesCalculate the observed altitude angle θ (t) to target at the earth's surface of sensor, described calculation formula
For：
Step 4：If target being observed meets solar irradiation simultaneouslyCondition and observed altitude angle θ (t) conditions are for example
Less than 90 degree, this moment is can shooting time.Carrying out calculating 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 computational methods of moon base earth observation electro-optical photo according to claim 1, it is characterised in that：It is described
The longitude and latitude position that moon based sensor in step 1 is laid on the moonBiography under the solid referential of ground heart
Sensor coordinateConversion process in need moon based sensor longitude and latitude position firstBecome
Change to moon heart moon solid three-dimensional coordinateThen three-dimensional coordinate transformation will be consolidated the heart moon moon to moon heart Celestial Reference System
Coordinate, further by moon heart Celestial Reference System coordinate transform to the earth's core Celestial Reference System coordinate, combine further according to International Astronomical
Meeting (IAU) resolve and International Earth Rotation service (IERS) specification by the heart moon moon solid three-dimensional coordinate transformation to ground heart referential admittedly
Under, above-mentioned transformation relation can be calculated by 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 framework
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 moon heart and the trivector in the earth's core, [L] is lunar libration spin matrix, can pass through ephemeris number
Obtained according to calculating.
3. the shooting time computational methods of moon base earth observation electro-optical photo according to claim 1, it is characterised in that：It is described
The process that the light conditions of target being observed are calculated in step 2 is as follows：
It is primarily based on coordinate position of the sun under the solid referential of ground heartCalculate its corresponding substar position
PutIt is then converted into latitude and longitude coordinates positionBy rotation transformation by coordinate position (φt；λt) limit P is transformed to, now
The longitude and latitude position of target being observedIt will be transformed to limit P difference of latitude：
φT=arccos (cos φtcosφcosλtcosλ+cosφtcosφsinλtsinλ+sinφtsinφ) (3)
Difference of latitude after conversion is less than such as 90 degree of certain angle, then setsFor 1,0 is otherwise set to；
4. the shooting time computational methods of moon base earth observation electro-optical photo according to claim 1, it is characterised in that：It is described
In step 4 can shooting time be between the starting and ending time of selection, and meet solar irradiation simultaneouslyBar
The period of part and observed altitude angle θ (t) conditions.
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 true CN107421504A (en) | 2017-12-01 |
| CN107421504B 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) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108896038A (en) * | 2018-05-15 | 2018-11-27 | 中国科学院遥感与数字地球研究所 | Month base optical sensor imaging method |
| CN108931259A (en) * | 2018-05-25 | 2018-12-04 | 中国科学院遥感与数字地球研究所 | A kind of method of month base earth observation image simulation |
| 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 |
Citations (4)
| 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 |
| US20110187860A1 (en) * | 2008-08-28 | 2011-08-04 | Centre National D'etudes Spatiales (Cnes) | Compression of earth observation satellite images |
| 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 |
-
2017
- 2017-08-02 CN CN201710650254.1A patent/CN107421504B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110187860A1 (en) * | 2008-08-28 | 2011-08-04 | Centre National D'etudes Spatiales (Cnes) | Compression of earth observation satellite images |
| 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 |
|---|
| 郭华东，丁翼星，刘广，张道卫，傅文学，张露: "面向全球变化探测的月基成像雷达概念研究", 《中国科学:地球科学》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108896038A (en) * | 2018-05-15 | 2018-11-27 | 中国科学院遥感与数字地球研究所 | Month base optical sensor imaging method |
| CN108896038B (en) * | 2018-05-15 | 2022-03-25 | 中国科学院遥感与数字地球研究所 | Imaging method of lunar-based optical sensor |
| CN108931259A (en) * | 2018-05-25 | 2018-12-04 | 中国科学院遥感与数字地球研究所 | A kind of method of month base earth observation image simulation |
| CN108931259B (en) * | 2018-05-25 | 2021-02-05 | 中国科学院遥感与数字地球研究所 | Method for simulating lunar-base earth observation image |
| 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107421504B (en) | 2019-11-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Tang et al. | Overview of the GF‐7 laser altimeter system mission | |
| US9551980B2 (en) | Solar timer using GPS technology | |
| CN102706363B (en) | Precision measuring method of high-precision star sensor | |
| CN102288201B (en) | Precision measurement method for star sensor | |
| CN107421504B (en) | The shooting time calculation method of month base earth observation electro-optical photo | |
| CN105486314A (en) | Lagrange navigation constellations for seamless coverage of moon space, and construction method thereof | |
| CN104462776A (en) | Method for absolutely radiometric calibration of low orbit earth observation satellite with moon as reference | |
| CN102288200B (en) | Accuracy measurement system for star sensor | |
| CN103837150A (en) | Method for performing rapid celestial fix through CCD (charge coupled device) zenith telescope on ground | |
| CN107144283A (en) | A kind of high considerable degree optical pulsar hybrid navigation method for deep space probe | |
| CN105446346A (en) | Moon relative calibrating posture adjustment method by remote sensing satellite | |
| CN104236553B (en) | Autonomous all-weather stellar refraction satellite location method | |
| CN105486315A (en) | Method for adjusting remote sensing satellite absolute calibration attitude relative to moon | |
| 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 | |
| Hawkins | Callanish, a Scottish Stonehenge: A group of standing stones was used by Stone Age man to mark the seasons and perhaps to predict eclipse seasons. | |
| CN109099911A (en) | A kind of method and system of air line navigator fix | |
| Wang et al. | Cloud‐tracked winds for the first Mars Global Surveyor mapping year | |
| CN108931259B (en) | Method for simulating lunar-base earth observation image | |
| Shaltout et al. | Studying the orientations of luxor ancient Egyptian temples using QuickBird images | |
| CN106767824A (en) | A kind of method for calculating double detector in objects outside Earth surface relative position | |
| Tao | Space Optical Remote Sensing: Fundamentals and System Design | |
| Zhang et al. | Random Star Recognition Algorithm Based on Image Total Station and Its Application to Astronomical Positioning |
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 |