Next Article in Journal
Comparison of Five Spatio-Temporal Satellite Image Fusion Models over Landscapes with Various Spatial Heterogeneity and Temporal Variation
Previous Article in Journal
Reconstructing One Kilometre Resolution Daily Clear-Sky LST for China’s Landmass Using the BME Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Constructing a High-Accuracy Geometric Model for Moon-Based Earth Observation

1
Key Laboratory of Digital Earth Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100094, China
2
School of Geography, University of Leeds, Leeds LS2 9JT, UK
3
School of Earth Sciences and Engineering, Nanjing University, Nanjing 210023, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(22), 2611; https://doi.org/10.3390/rs11222611
Submission received: 30 September 2019 / Revised: 1 November 2019 / Accepted: 2 November 2019 / Published: 7 November 2019

Abstract

:
The Moon provides us with a long-term, stable, and unique place for Earth observation. Space agencies of various countries, including the United States, China, and Italy, have made the realization of Moon-based Earth observation an objective of their lunar missions. To date, although some conceptual studies have been presented, an accurate geometric model for Moon-based Earth observation has not yet been described. This paper introduces such a geometric model, which connects the attitude of a Moon-based sensor with a corresponding field of view on the Earth’s surface. The aberration and light time correction are involved. Due to the lack of high-quality experimental data, one qualitative comparison with Chang’E lander images and another quantitative comparison with software output are made. The comparison results show good similarity. The overall model accuracy is evaluated to be better than 4″ at current stage and will be better than 1.5″ if the seleno-graphic position is accurately determined. In direct geolocation process, the aberration and light time will cause a total 0.7 km deviation on the ground near the sublunar point. In SAR range history simulation, the light time effect will lead to a linear error, as large as tens of meters, throughout the integration time.

Graphical Abstract

1. Introduction

No matter how far humans want to travel in the universe, the Moon, the only natural satellite of the Earth and our nearest neighbor, is the first place we need to colonize. Since 1958, satellites and spacecrafts have visited the Moon more than 100 times. Countries and institutions including the Soviet Union-Russia, USA, ESA, China, Japan, and India have been involved in these explorations. In this century, more than 10 missions have been successfully launched, and tens of missions are scheduled. With technological progress and increases in our understanding of the Moon, the colonization of the Moon is no longer an envisaged possibility but is a real challenge humans face. Recently, the ESA proposed a “moon village” program, which is a Moon base under international collaboration among spacefaring nations for science, business, mining, and even tourism [1]. NASA also started its new lunar mission, named Artemis, to return astronauts to the lunar surface by 2024 [2]. Similar plans have emerged over past few decades, but until now, humans have not truly had the ability to realize these plans. Such a permanent settlement on the lunar surface will be significant to scientific studies and deep space exploration. For Earth science, it will also provide an opportunity to complete our Earth observation system. China’s Lunar Exploration Program (CLEP) officially began in 2004, and one of its three overall objectives is to develop the Moon as an Earth observation satellite for monitoring global climate change, eco-environment evolution, pollutant diffusion, natural disaster evolution, and other processes [3]. Recent research has shown that Moon-based Earth observation is characterized by longevity, integrity, stability, and uniqueness [4]. In addition, Moon-based Earth observation is a good complement to existing satellite Earth observation systems, which usually run in low Earth orbit or geostationary orbit. Remote sensors on the lunar surface can be used to monitor solid Earth dynamics, the global energy budget, climate and environmental change, and the near-Earth space environment. An Earth’s radiation budget measurement system is planned to land on the lunar surface with Chang’E follow-up mission in ten years [5]. It consists of a radiometry to measure the outgoing solar and thermal radiation of the earth and an array imaging spectrometry to map the atmospheric and surface composition of the earth disk. Concerning the large distance between the Earth and the Moon, a slight change in the angular parameter is expected to bring in a large deviation on the ground. Therefore, it is necessary to present a high-accuracy geometric model establishing the connection between the Moon-based sensor’s pointing and positioning parameters and a target on the Earth’s surface.
Two scientific Earth observation experiments have been carried out on the lunar surface. During the first Moon-based Earth observation experiment in 1972, Apollo 16 astronauts used a far-ultraviolet camera/spectrograph to observe the Earth’s atmosphere and corona [6]. Forty-one years later, a one-year Earth plasma observation experiment were successfully carried out by the Extreme Ultraviolet (EUV) camera on the Chang’E-3 mission lander [7]. These two experiments were mainly related to space science, independent of geographic coordinates. A simplified geometric model, which only consider the earth’s rotation but ignore the revolution of the moon, was presented for synthetic aperture radar (SAR) system performance analysis [8]. At the same time, ELP2000, a semi-analytic theory, was also selected to construct the geometry for the similar purpose [9]. These two models are good at analyzing the system performance without referring to the sensor’s pointing direction and its position on lunar surface, which are the only designable parameters in the geometry model. Ding et al. introduced the JPL ephemeris in geometry and built connections among sensor’s seleno-graphic position, point direction, and the doppler parameters of moon-based SAR [10]. Recently, a simplified geometric model is widely used for spatial and temporal scope analysis [11,12] and the simulation of earth’s outgoing radiation [13,14,15]. The moon-based SAR imaging geometry is also improved [16,17,18]. However, for the mission design of moon-based Earth’s radiation budget measurement and moon-based SAR imaging simulation, the model is still made with too much simplification. Although the satellite tool kit (STK) provides us with a powerful function library for most space applications, the tools for moon-based sensors is not so fully functional as the spaceborne ones.
Compared with former studies, the geometric model presented in this paper is characterized by: (1) a general geometric frame compatible with both passive optic sensors and SAR system, (2) reversable simulation flow, (3) minimum approximation to fulfill the stringent requirements and give reference to the simplification, (4) the analysis and compensation of the aberration effect and light time effect, and (5) first accuracy evaluation of moon-based earth observation geometric model. Moreover, the geometric model can be used to: (1) reveal the direct geometry relationship for mission design, SAR focusing processing, image correction etc.; (2) simulate data to develop data processing algorithms and software systems; and (3) forecast the sensor’s pointing direction for the mission operation, especially for the non-imaging radiometry and SAR.
A brief discussion about the moon-based sensors will be presented in Section 2. Then, in Section 3, we introduce, in detail, the geometric model and the method of simulation. The definition of spatial coordinate systems and temporal coordinate systems are presented in Section 3.1 and Section 3.2, respectively. The data referring to the position of the Sun, the position and attitude of the Earth, and the position, velocity, and libration of the Moon are introduced in Section 3.3. Section 3.4 and Section 3.5 show the transformation matrices and the process of transformation from a sensor-fixed coordinate to a geocentric coordinate system. The corrections of aberration and light time are shown in Section 3.6 and Section 3.7, respectively, for two passive optic sensors. The geometry simulation, as well as an attempt of validation, is described in Section 4. Then, the model accuracy is evaluated in Section 5. Two examples of application are given before the conclusion. One is predicting the pointing direction of a moon-based sensor and derived constraints on site selection. Another is the simulation of range history of moon-based SAR and its comparison with former approximation.

2. Brief Considerations Concerning Moon-Based Earth Observation Sensors

In 2006, NASA suggested 12 scientific objectives for Earth observation among a total of 160 lunar exploration objectives [19]. Although the ambitious plan to return to the Moon was abandoned due to the transfer of political power, scientists have not lost their enthusiasm for relative studies. In the decade after 2006, scientific objectives were elaborated, and sensor configurations were put forth [8,9,20,21,22]. We chose seven geo-scientific objectives from NASA’s list and rearranged them into three categories: (1) use of multispectral passive sensors at the disk scale for global mapping of atmospheric composition and surface composition, determining the Earth’s bidirectional reflectance distribution function (BRDF) and monitoring the Earth’s hot spots; (2) use of radiometers to measure radiance directed toward the Moon for Earth radiation budget research; and (3) use of SAR for topography, tomography, altimetry, vegetation mapping, and ice dynamics observations in polar regions.
At present, both the radiometry and the array imaging spectrometry are technically feasible to observe the Earth from the Moon. Due to their cost, size, weight, and energy consumption, Moon-based sensors are quite limited. In view of the trade-off between costs and scientific significance, the horizontal resolution in the visible spectrum was initially set at 1–10 km. However, if a permanent manned Moon base is established in the future, the transport, construction, and maintenance costs of the observation system will be greatly reduced by cost allocation, especially for long-lasting missions. A high-resolution optical system for land mapping, source investigation, and target monitoring is also a candidate of the observatory in the future. In long-term planning, a Moon-based SAR can be used to image the global dynamics of the solid Earth. The antenna size and power requirements of such a SAR are very high, so it can only be attached to a permanent Moon base over a certain scale. The suggested spatial resolution is 10–100 m. An L-band system is better for some geoscience applications, but an X-band or C-band system is smaller in size [22].
To improve efficiency and reduce the possibility of contamination, it is better to use a field of view (FOV) just fit to the Earth disk. The diameter of the Earth disk including the atmosphere is slightly smaller than 13,000 km. Therefore, for the closest case, the FOV or footprint must have a diameter of approximately 2.1° with 1%, or, equivalently 75.6″ tolerance. However, if a smaller FOV or footprint is used to cover a target area on the Earth’s surface, especially for SAR, the requirements will become much more stringent. For example, a footprint with a width of approximately 300 km and a 3% tolerance of deviation will allow a maximum deviation of only 9 km, corresponding to 5.2″. The accuracy of presented geometric model is expected to be better than this value.

3. Model and Simulation

3.1. The Spatial Coordinate Systems

The four coordinate systems used in this paper, namely, the International Celestial Reference System (ICRS), International Terrestrial Reference System (ITRS), Barycentric Celestial Reference System (BCRS) and Geocentric Celestial Reference System (GCRS), are rigorously defined by the International Astronomical Union (IAU) and the International Union of Geodesy and Geophysics (IUGG). Their definitions are presented in the Nomenclature for Fundamental Astronomy (NFA) glossary. The BCRS and GCRS are assumed to be oriented with the ICRS axes according to the IAU 2006 Resolution B2. To validate most equations in this paper, the time-space distortion will not be taken into consideration.
In addition, the meanings of several other coordinate systems have to be specified. Similar to the geographic coordinate system (GCS), the seleno-graphic coordinate system (SGS) is set up to define Lat-Lon positions on the lunar surface, and its origin is considered the intersection point of the X-axis of the body-fixed coordinate system and the lunar reference sphere. The lunar reference ellipsoid is a sphere with a radius of 1737.4 km. In addition, the three other right-handed three-dimensional space coordinate systems used in the following sections are:
Sensor Coordinate System (SCS): The SCS is a sensor-fixed reference frame. The origin is the optical center for the camera or the phase center for the antenna. The Z-axis is the optical axis or the direction of the main lobe, namely, the line of sight (LOS) in the following sections. The XY plane is perpendicular to Z. The directions of the X-axis and Y-axis refer to the antenna shape or the image coordinates. In this paper, they are implicitly defined by the zenith angle and azimuth angle.
Observatory Local Coordinate System (OLCS): The origin of the OLCS is an artificially defined center of the observatory with an elevation of H. The Z-axis is the same as the local normal vector, while the XY plane is the local tangent plane. The X-axis points to the north. The zenith angle- and the azimuth angle-related SCS and OLCS will be introduced in Section 3.4.
Selenocentric Body-fixed Coordinate System (SBCS): The coordinate origin is the lunar mass center. Principle axis (PA) coordinates and mean-Earth/rotation-axis (MER) coordinates, which are related by a constant three-angle rotation, are two commonly used body-fixed coordinate systems [23]. The PA coordinates are directly linked with the ICRS by physical librations [24], while lunar cartography customarily refers to MER [25].

3.2. Considerations Concerning the Time Coordinates

The Earth orientation parameter (EOP) data and JPL ephemeris used in this paper are presented in different time coordinates. The terrestrial time (TT) coordinating the Earth’s attitude has a mean rate close to the proper time of the rotating Earth and can accurately be represented by International Atomic Time (TAI), an idealized atomic time system with the SI second on the geoid as the unit of measurement [26]:
T T = T A I + 32.184   second
Barycentric dynamical time (TDB) serves as an independent time argument of barycentric ephemerides, such as the JPL ephemeris. The relationship between TDB and TT is given by [27]:
T D B = T T + P
where P is a periodic term with a main component with an approximately 1.7 ms amplitude and a one-year period. P leads to periodic geometry deviation during the Earth’s orbital cycles; however, its amplitude is too small for the Earth observation application. The maximum space deviation of the lunar mass center is less than 2 m, much smaller than other errors analyzed by Section 5.3. Universal Time Coordinated (UTC) system, which is the standard for all civil time systems, is different from TAI because the former has integer leap seconds forecasted by the International Earth Rotation Service (IERS). The time conversion is realized by the Standards of Fundamental Astronomy (SOFA) time scale and calendar tools in this paper.

3.3. JPL Development Ephemeris and Earth’s Attitude

The JPL ephemeris is in good agreement with lunar laser ranging (LLR), which now achieves millimeter-scale accuracy [28]. In this paper, we use DE430 to calculate the positions and velocities of the helio-center, geocenter, and seleno-center in the ICRS, as well as the amplitude and derivative of the libration. The position vector should be translated from the ICRS to the GCRS for the next step.
Usually, the position of the reflector on the Earth is given in the GCS, defined on the Earth’s reference ellipsoid. The coordinates in the ITRS can be calculated by using the geographic coordinates and parameters of the ellipsoid. Since the center of the Sun, Earth, and Moon are given in the GCRS, we need to understand the orientation matrix between the ITRS and GCRS that describes the Earth’s attitude in the universe. According to IERS conventions, the orientation matrix M I T R S _ G C R S is presented by [29]:
M I T R S _ G C R S = M P N · M E R · M P M
where the M P N is the transform matrix of the precession and nutation, M E R is the transform matrix of the Earth’s rotation, and M P M is the transform matrix of the polar motion. The matrix M P N · M E R · M P M transforms the position vector from the ITRS to the GCRS. The IERS conventions shows how to use the models and data to calculate M I T R S _ G C R S . The SOFA tools for Earth attitude provide the program that implements the official IAU algorithms based on EOP data.

3.4. Orientation Matrices

The relation between the sensor’s position and pointing direction is given by a series of transformation and an intersection. The flow chart of simulation is presented in Section 4. The orientation matrices for transformation are introduced in this section, and the process of transformation will be introduced in the next section. Assuming the coordinate system rotates θ clockwise around the X , Y , and Z axes, the rotation matrices are:
T x ( θ ) = [ 1 0 0 0 cos θ sin θ 0 sin θ cos θ ] , T y ( θ ) = [ cos θ 0 sin θ 0 1 0 sin θ 0 cos θ ] , T z ( θ ) = [ cos θ sin θ 0 sin θ cos θ 0 0 0 1 ]
The vector in the new coordinate system can be calculated by pre-multiplying the unit rotation matrix to the vector in the original coordinates system.
For clarity, we assumed the sensor is on a two-dimensional altazimuth turntable [30,31]. In its initial state, the SCS and OLCS are the same. When controlling the direction of the LOS, we first rotate the azimuth angle α I clockwise around the Z-axis and then rotate the zenith angle θ I clockwise around the Y-axis. In some cases, a third rotation is needed. For example, if we use a rectangular antenna, the long side of the footprint is expected to be parallel with the meridian, and an extra T z should be added. The rotation angle α R is the clockwise angle around the Z-axis, and it is zero for circular aperture. Thus, the orientation matrix M S C S _ O L C S (Figure 1) can be written as:
M S C S _ O L C S = T z ( α I ) T y ( θ I ) T z ( α R )
If the observatory position is ( ξ M , M ) in the SGS, we can get the orientation matrix M O L C S _ M E R (Figure 1).
M O L C S _ M E R = T z ( π + M ) T y ( ξ M π / 2 )
The orientation matrix M M E R _ P A connecting the MER to PA is:
M M E R _ P A = T z ( β 1 ) T y ( β 2 ) T x ( β 3 )
The angles β 1 , β 2 , and β 3 refer to the constant parts of the three libration parameters. Their values compatible with DE430 are 67.573″, 78.580″, and 0.285″ [24]. If the position of the observatory is directly given by PA coordinates ( x P A , y P A , z P A ) T , such as the coordinates of the LLR retroreflectors, the orientation matrix from the OLCS to the PA can be directly expressed by substituting the PA latitude and longitude for the MER latitude and longitude in Equation (6).
The three Euler angles of libration are Ω ^ L , i ^ L , and u ^ L . On the celestial sphere, let N 1 stand for the ascending node of the lunar mantle equator to the ICRS equator (Figure 2). Thus, it can be inferred that the three Euler angles are the angle from the X-axis of the inertial frame along the ICRS equator to the intersection of the lunar mantle equator, the inclination of the lunar mantle equator from the ICRS equator, and the longitude from N 1 to the lunar prime meridian along the lunar mantle equator [25,32]. Thus, the orientation matrix M P A _ G C R S can be presented by:
M P A _ G C R S = T z ( Ω ^ L ) T x ( i ^ L ) T z ( u ^ L )
The orientation matrix M G C R S _ I T R S is already given in Section 3.4 and is calculated by SOFA. Note that the orientation matrix is orthogonal and has the following property:
M 1 = M T

3.5. Transformation from the SCS to the GCRS

In the SCS, the LOS can be presented by:
L S C S = ( 0 , 0 , 1 ) T
In the OLCS, the LOS can be presented by:
L L O C S = M S C S _ O L C S · L S C S
The position of the sensor center (PSC) in the OLCS is S O L C S . This vector relies upon the system design. For simplicity, we assume the two centers coincide in the simulation.
In the MER, the LOS can be presented by:
L M E R = M O L C S _ M E R · L O L C S
And the PSC can be presented by:
S M E R = M O L C S _ M E R · S O L C S + ( R M + H ) ( cos M cos ξ M , sin M cos ξ M , sin ξ M ) T
where the R M is the radius of the Moon, and H is the elevation of the observatory center.
In the PA, the LOS and PSC can be presented by:
L P A = M M E R _ P A · L M E R
S P A = M M E R _ P A · S M E R
In the GCRS, they are presented by:
L G C R S = M P A _ G C R S · L P A
S G C R S = M P A _ G C R S · S P A + M G C R S
where M G C R S is the position of the Moon in the GCRS.

3.6. The Light Aberration

Light takes approximately 1.3 s to travel from the Earth to the Moon, which is much longer than the spaceborne observation. Thus, an obvious geolocation problem for Moon-based Earth observation is that the geometry changes from the time when the light leaves the Earth to the time when the light arrives at the sensor. This phenomenon, known as light aberration in astronomy, results from the relative motion between the reflector and the observer. Recently, high-resolution remote sensing processing has taken aberration compensation into consideration [33,34]. This effect is calibrated by introducing a corrective rotation. Intrinsically and more precisely, a two-step method is developed for moon-based earth observation, and it also can be used in spaceborne cases [35]. For Earth observation, the aberration can be divided into two parts. The part related to the motion of the sensor will be discussed in this section, and the part caused by the Earth’s rotation, called the light time effect, will be discussed in Section 3.7. The aberration angle can be approximately presented by:
ε = v c sin θ S
where θ S is the viewing angle related to the speed vector, c is the speed of light, and v is the speed of the observer. If the observing angle is perpendicular to the relative speed, the angular aberration can easily be derived and is v / c . For example, if the orbital speed of the Moon is 1.02 km/s, the viewing angle is 90°, and the distance between the reflector and the seleno-center is 380,000 km; thus, the light aberration is approximately 3.4 × 10−6 rad (0.7″), corresponding to a 1.3 km deviation on the Earth’s surface.
In the GCRS, the velocity of the sensor is V G C R S , and the real received light is c · L G C R S . The calculated received light is c · L G C R S and is in the moving system that has the same axis directions as the GCRS but moves with the SBCS. Let the received light in the GCRS be u = c · L G C R S , the received light in the moving system be u = c · L G C R S , and the relative speed of the GCRS in the moving system be v . Therefore, the speed addition can be described by the Lorentz transformation:
u = u + [ ( γ 1 ) v · u v 2 + γ ] v γ ( 1 + v · u c 2 )
where v is the magnitude of v = V G C R S and the factor γ = 1 / 1 v 2 / c 2 . Then, the real received light is:
L G C R S = u / c
It can be easily deduced that the sensor speed vectors in the coordinates fixed on the Moon, i.e., V S C S , V L O C S , V M E R , and V P A , are all 0. Therefore, the GCRS speed vector is:
    V G C R S = d d t ( M P A _ G C R S · S P A ) + V M G C R S = d M P A _ G C R S d t · S S B C S + M P A _ G C R S · V P A + V M G C R S = d M P A _ G C R S d t · S S B C S + V M G C R S
where V M G C R S is the velocity of the seleno-center in the GCRS, and:
d M P A _ G C R S d t = M P A _ G C R S Ω ^ L d Ω ^ L d t + M P A _ G C R S i ^ L d i ^ L d t + M P A _ G C R S u ^ L d u ^ L d t
Formulating Equation (22) is a complex task [36,37]. For some applications, d Ω ^ L d t and d i ^ L d t can be considered approximately 0, and d u ^ L d t can be considered 13.2 degrees per day [10]. In this paper, the result is numerically computed by using the Chebyshev polynomial to interpolate the JPL ephemeris.

3.7. The Intersection Point

In the previous calculation, the time-variant parameters refer to the receiving time t 2 , so t 2 is not explicitly written out in the equations. However, the position of the reflector refers to the ‘departure time’ t 1 , namely, the time when light leaves the reflector, in the geometry. Thus, the time difference should be clearly represented in the following equations. The departure time can hardly be measured by an Earth observation system, but the relationship below stands (Figure 3):
A G C R S ( t 1 ) = S G C R S ( t 2 ) + c ( t 2 t 1 ) · L G C R S
where L = c ( t 2 t 1 ) is the length of the LOS.
Meanwhile, at t 1 in the GCRS, A G C R S ( t 1 ) is on the terrestrial ellipsoid, so we get:
[ M G C R S _ I T R S ( t 1 ) A G C R S ( t 1 ) ] [ 1 a e 2 0 0 0 1 a e 2 0 0 0 1 b e 2 ] [ M G C R S _ I T R S ( t 1 ) A G C R S ( t 1 ) ] = 1
where A G C R S ( t 1 ) represents the coordinates ( X A ( t 1 ) , Y A ( t 1 ) , Z A ( t 1 ) ) T of the intersection point of the LOS and the terrestrial reference ellipsoid in the GCRS at time t 1 , and L = c ( t 2 t 1 ) is the length of the LOS. The major semi-axis a e is 6378.137 km and the minor semi-axis b e is 6356.752 km in terms of WGS84 ellipsoid. By solving the equation, we can determine the position of the reflector A I T R S by substituting the estimated t 1 into M G C R S _ I T R S ( t 1 ) A G C R S ( t 1 ) .
Because M G C R S _ I T R S ( t 1 ) is inexplicit, Equation (24) only has numerical solutions. Here, a simple approximation with a max error less than 10 m is presented. Note that if the LOS and the ellipsoid intersect, the range of t 1 is less than 21.5 ms because:
t 2 s 2 + a e 2 c < t 1 < t 2 s a e c
where s is the length of S G C R S ( t 2 ) . We can estimate t 1 by replacing M G C R S _ I T R S ( t 1 ) with M G C R S _ I T R S ( t 2 s 2 + a e 2 c ) or M G C R S _ I T R S ( t 2 s a e c ) in Equation (24). In this time window, the ITRF is under- or over-rotated with a 10 m max error (the Earth’s mean rotation speed is approximately 7.29 × 10 5 rad / s ) corresponding to an error of 3.3 × 10 8   s in the estimation of t 1 . In this paper, M G C R S _ I T R S ( t 2 s a e c ) is chosen because the case wherein the LOS is pointed to the geocenter is mainly discussed.
In particular, imaging spectrometer models are usually described by the colinear equation, which is an effective planar approximation within limits. To address the whole disk, we follow the same procedure but substitute L O S _ X Y for L S C S in Equation (10):
L O S _ X Y = 1 x 2 + y 2 + f 2 ( x , y , f ) T
where f is the focal length. Then, all image coordinates ( x , y ) are connected to geographic coordinates. This method is computationally inefficient but has a high accuracy.

4. Geometry Simulation

The simulation procedure follows the model specified from Section 3.4, Section 3.5, Section 3.6 and Section 3.7. DE430 from JPL and EOP 14 C04 from IERS are two input data tables. Usually, times are initially presented in UTC and should be converted to TDB and TT to interpolate the two input tables. The calculated results are tested by the attached test data and the results from the IERS Rapid Service/Prediction Center (http://maia.usno.navy.mil/t2crequest/t2crequest.html). Due to the libration and the ellipticity of the orbit of the Moon, the LOS will shift out of the Earth disk if the azimuth angle and zenith angle do not change with the revolution of the Moon. To make Equation (24) always solvable, a simple strategy is to point the LOS to the geocenter:
L G C R S = S G C R S ( t 2 ) S G C R S ( t 2 )
Then, the coordinates of the intersection point and corresponding azimuth angle and zenith angle can be calculated by a slight adjustment of the proposed method (Figure 4). The coordinates of the lunar surface are directly given in the PA. The following analyses use these criteria, unless otherwise specified.
A lack of experimental data makes it very hard to validate simulation results. Optically imaging the Earth is not one of the objectives of the Chang’E-3 mission. Although the Terrain Camera (TCAM) lander did take several photos of the Earth, the spatial resolution is only approximately 70 km. Since these images are the only measured data accessible to us, we compromise with a comparison of the imaged Earth disk and the simulated Earth disk (Table 1). The elevation of the lander is 2632.0 m, and the location is ( 44.1206 , 19.5124 ) in MER.DE421 [38], i.e., ( 44.1412 , 19.5382 ) in PA.DE430. The orientation angles differ slightly from the angles in Equation (7) [32]. In addition, the simulation results from the STK 11 are introduced as mutual authentication data. A more credible validation method will only rely on field measurements in future missions.
A qualitative comparison is shown in Table 1. The time is given in UTC +8. The STK results are the views from lunar center, while the images of Chang’E-3 and the simulation results of this paper are viewing from the lander. It can be inferred that the twilight arcs from three sources are largely the same (note that the attitude of the Earth and the rendering effect are slightly different). The disk views from the STK and this paper are similar, but the disk views from the Chang’E-3 lander are unidentifiable due to the low spatial resolution and high cloud coverage. Figure 5 shows infrared global geostationary composite v1 at UTC 2013/12/24 18:15. The cloud texture in the red crescent is very close to that in the Chang’E image acquired four minutes later.
A quantitative comparison between the intermediate result and the STK outputs is also made. Five groups of random numbers corresponding to the timestamp, lunar latitude, lunar longitude, azimuth bias, and zenith bias are combined to generate 100 check points. The aberration correction and light time compensation are not made in this comparison to keep pace with the software. Because the STK does not give built-in functions to report the results for a Moon-based sensor, we estimate the intersection position manually in the 3D interface point by point by using a conic beam as narrow as 2 × 10−9 degrees. This method is very time-consuming. The comparison results are shown in Figure 6. All geo-distance deviations are smaller than 400 m and have no significant correlation to any parameters. The mean value and standard deviation are 108.6 m and 95.6 m, respectively. Interestingly, the data between 2009/01/01 and 2016/01/01 obviously fit much better than the data from other times (between two blue dashed lines). Even from 2010/01/01 to 2015/07/01, all the differences are smaller than 2 m. We think that this result may be related to the release date of the used software version and the operation of the ICRF since 2009/01/01, according to the STK help document.

5. Model Accuracy Evaluation

5.1. The Scale of the Aberration Effect

According to Equation (18), the typical deviation of a moon-based sensor is about 1.3 km on the Earth’s surface. The light time in Equation (23) is approximately 1.27 s and therefore leads to a 0.6 km deviation in geolocation. In most cases, as a result of the homodromy of the Earth’s rotation and the Moon’s revolution, the deviations cancel each other out, leaving a 0.7 km remainder. More detailed results are given in Figure 7. We first use Equation (27) to calculate the zenith angle and azimuth angle on the basis of the flow chart in Figure 4 and then reverse the procedure to recalculate the LOS and intersection points (1) when the aberration is not corrected, (2) when the light time is not compensated, (3) without either treatment, and (4) with both treatments. The sensor is located at ( 60 , 0 ) . The differences between cases (1–3) and case (4) reveal the scale of these two effects. In Figure 7, the red dotted line represents the aberration effect, the blue dashed line represents the light time effect, and the black solid line represents the comprehensive effect. The geo-distance between corresponding points, used as an index, is in line with the analysis in the initial part of this paragraph. It should be noted that the comprehensive effect is weakly related to the seleno-graphic coordinates of the sensor but strongly connected with the position of the reflector. For example, if the reflector is near the edge of the Earth disk, the comprehensive effect can be larger than 10 km. One reason for this is that the same aberration angle will result in a larger projected geo-distance at the edge of Earth disk than at the disk center; in addition, the light time only causes small geo-distance deviations at high latitudes.
The aberration effect and light time effect of Moon-based Earth observation are two orders of magnitude larger than those of satellite imaging. Although the pointing accuracy requirements of passive optic sensors are not very high, they are considerably higher than the model accuracy analyzed in Section 5.3. Moreover, the compensation is essential if the model is used as a direct geo-location or correction model for array spectrometry image.

5.2. Pointing Deviation

There are some parameter uncertainties. We describe the general relationships between the errors of different parameters and their corresponding geo-distance deviations with examples from 2017 and the maximum and minimum deviations in the past 50 years. The LOS is pointed to the geocenter, and the relationship between the parameter errors and geo-distance deviations is approximately linear around the Earth’s disk center.
Roughly, one second angular error will cause a geo-distance deviation of 1.84 km if the distance is 380,000 km. The spatial error will be directly projected on the Earth’s surface with a scaling factor referring to the projection angle and Earth’s curvature. The pointing errors attributed to some sources over time are shown in Figure 8. The figure only shows results from 2017, when the sensor was at ( 60 , 0 ) ; however, the geo-distance deviations of the other five sites over the past 50 years have also been investigated (Table 2). Generally, an angular error of one second is likely to generate a geo-distance deviation less than 1.94 km, while a 10 m positional error is likely to induce a geo-distance deviation less than 10 m. The ground deviation caused by the errors in zenith angle or lunar longitude are relatively not sensitive to time and seleno-graphic coordinate. They have more stable impact than the errors in other designable parameters. The positional errors can be converted to equivalent angular errors. The conversion refers to the nearest Earth-Moon distance to make sure that the maximum value is given. Therefore, the model accuracy can be quantitatively restricted by the sum of all kinds of angular errors.

5.3. Model Accuracy

The error sources are divided into seven categories: Orbit, libration, seleno-graphic, pointing, propagation, EOP, and geographic (Figure 9). The maximum error of each source will be evaluated.
In DE430, libration errors are coupled with orbit errors of approximately 0.1″. For the DE430 gravity field and librations, the uncertainty of the retroreflector PA coordinates is only 0.12–0.27 m [24]. If the error is 0.27 m, the deviation angle is approximately 0.032″. The coordinate accuracy of other sites is much lower. The coordinates of the Apollo lunar surface experiment packages may have errors as large as 30 m [39]. For cartography, the widely used Unified Lunar Control Network 2005 (ULCN 2005) is based on Clementine images acquired in 1994 and the previous ULCN. The solution root mean square is 0.9 pixels, corresponding to 100 m–300 m on the lunar surface [40]. The accuracy of next generation networks may reach 20 m at the resultant tie points [41]. For topography, the Lunar Reconnaissance Orbiter (LRO) Lunar Orbiter Laser Altimeter (LOLA) Digital Elevation Model (DEM) represents more than 6.5 billion measurements, and the average accuracy of each point after crossover correction is better than 20 m horizontally and approximately 1 m radially [42]. Therefore, the lunar map coordinate accuracy can be considered 20 m (2.4″) for our simulation studies. After settlement on the Moon, the observatory position accuracy is expected to reach 0.2 m.
In astronomy, a telescope or radio telescope is required to point to a star or radio source whose right ascension and declination position is already known. The position of a star or radio source seldom changes, so errors are mainly derived from the pointing accuracy of the telescope. At present, the pointing accuracy of a telescope after correction is typically 1″on the Earth’s surface [31,43]. Without atmospheric distortion and with slight crustal deformation variation, the pointing accuracy of a Moon-based sensor may be higher.
The atmospheric refraction is related to wavelength. For visible light, the ground deviation is usually 5 m (0.003″) when the local incident angle is 45° [34]. This deviation will be larger at a large incident angle. For microwaves, the ionospheric effect is more significant. The ground deviation can reach 100 m in the ultrahigh frequency band [44]. The stratospheric effects are much smaller. Thus, we assume the total refractive deviation is 150 m, corresponding to 0.086″.
The EOP and geographic errors usually refer to the geo-center, and should be converted to refer to the moon-based sensors in our analysis. In Moon-based Earth observation, the target position in the right ascension and declination relative to the sensor is always changing, so the transformations introduce extra errors. The model should be reversed to calculate the right ascension and declination target position. Rewrite Equation (23) as:
M I T R S _ G S R S ( t 1 ) A I T R S ( t 1 ) = S G C R S ( t 2 ) + c ( t 2 t 1 ) · L G C R S
where t 1 and L G C R S are unknown but A I T R S ( t 1 ) is already given. An iteration is used to solve the equation. The initial error is 10 m or 0.006″, and almost no estimated error remains after several iterations. However, the SOFA algorithm to calculate the orientation matrix has a max error of approximately 0.9″ in 21st century [45]. Relative to a Moon-based sensor, this error is less than 0.016″.
The impact of the Earth’s DEM can be easily corrected. In the forward model, the intersection with the DEM should be numerically coupled with Equation (24). In the reversed model, the elevation is already contained in A I T R S ( t 1 ) . However, each DEM contains intrinsic errors in relation with a particular terrain and land cover type [46]. The global assessment result shows that for Shuttle Radar Topography Mission (SRTM) DEM the absolute geolocation error is less than 15 m and the absolute height error is less than 10 m [47]. Thus, the geographic error is expected less than 25 m, or 0.015″ referring to the sensors.
The estimated accuracies are summarized in Table 3. All other errors are assumed to be within 0.1″. The main errors are the pointing control error and the position error on lunar surface. A great boost on seleno-graphic accuracy is expectable after the settlement on the Moon. At that time, the lunar surface positioning error will be as small as the position error of the LLR reflectors, and the overall model accuracy will be better than 1.5″. It is high enough for mission operation and image geolocation. At this current stage, if we use the LRO LOLA DEM to simulate data, the overall model accuracy is about 4″. Even in this case, the accuracy is acceptable compared with the limitations presented at the end of Section 2.

6. Examples of Application

6.1. Predicting the Pointing Direction and the Observational Constraints of Site Selection

In passive optic sensor cases, to point the LOS to the geocenter, the sensor attitude is determined by only the location of the observatory and the height of the sensor. Figure 10 shows the change rule of the azimuth angle and zenith angle when the height is set to 0. The seleno-graphic coordinates of the sensors in subfigures (a), (b), and (c) are (0,0), (0,80), and (80,0), respectively. The black solid line shows the zenith angle, the red dotted line shows the azimuth angle and the blue dashed line shows the angle between the LOS and the negative light direction toward the sensor (LOS-SUN angle). The negative zenith angle results from the definition of the rotation direction. As stated in Section 5.3, the overall prediction accuracy is expected to be less than 1.5″ for mission operation.
The LOS-SUN curve in Figure 10 indicates the risk of the sunlight going directly into the sensor. The closer the curve is to 0, the higher the risk is. The similarity of the curves in different subfigures suggests that the LOS-SUN angle is only slightly affected by the sensor’s position. Thus, its monthly approach to 0 is inevitable. To reduce the probability of the Sun entering the FOV, we suggest the use of a small FOV that is slightly larger than the Earth disk.
The irregular relative motion causes the projection of the Earth on the lunar surface to move around the seleno-graphic origin, and the selenocentric angle is typically less than 10°. If the sensor is put in this area (e.g., a circle around the seleno-graphic origin with a 10° radius), the azimuth angle must vary from 0° to 360° in each month to track the Earth (Figure 10a). Otherwise, the azimuth angle will be biased to a smaller range, and the farther the sensor is from the origin, the slower the angle will vary (Figure 10). This necessitates a low-frequency operation of the turntable.
Unlike the azimuth angle, no matter where the sensor is put, the zenith angle always varies within a range of a dozen degrees. However, if the absolute value of the zenith angle is too large, the signal from the Earth will be contaminated by the reflected or emitted energy from the lunar surface. To avoid this contamination, the absolute value of the zenith angle should be smaller than 90 ° FOV / 2 . In Figure 10b,c, the absolute value of the zenith angle is larger than 85° in some time periods. Therefore, places outside the 80° circle should be carefully examined.
Figure 11 shows a 19-year simulation on how the observatory position impacts the observational parameters, including the percentage of full coverage, the mean derivative of elevation angle and azimuth angle, and the number of times when LOS-SUN angle is smaller than FOV/2. As many former studies shows, the full coverage area, where the FOV are not sheltered all the time, is almost in an 81-degree circle around the seleno-graphic origin on lunar surface. At the edge of the Moon’s nearside, the full coverage rate decreases to 0.4 (Figure 11a). That means the sensors should be put inside the 81-degree circle to access to full coverage data at any time. Figure 11b,c shows the derivative of the two pointing angles. The change of elevation angle is slow but the azimuth changes intensively near the seleno-graphic center. The azimuth adjustment rate will increase 10 times if the observatory is moved from the edge of the Moon’s nearside to the center. That will result an extra burden for SAR pointing system. To minimize the pointing adjustment, the observatory should be put far from the center, e.g., areas close to the 81° circle on the lunar near side hemisphere. Compared with the south pole regions, the north polar regions are more likely to be affected by the direct sunlight. The best places to avoid sunlight (the red circle in Figure 11d) is the north tropical zone on the west part of the lunar nearside, about 5% better than the north polar regions. Therefore, after the comprehensive analysis of four subfigures, (15, −75), near the Crater Krafft, may be one of the best sites for our moon-based earth observatory, due to a relatively low possibility of direct sunlit, slow adjustment of pointing direction, and full coverage of the earth disk all through the time.
The south pole area is considered a hot candidate for the permanent Moon base [48,49]. Because of their high altitude, some regions near the south pole have good Earth visibility [50]. On the other hand, since the outer rim of the South Pole–Aitken basin at the near side can be seen from the Earth, this region may provide a suitable place for the attached observatory [51]. According to Figure 11, the area near (−81,0) are not the best but also a quite good place. Another crucial factor in this area for the visibility analysis is the complex terrain. More detailed analyses require high-precision local DEM, and thus will be reported in following papers.

6.2. Simulating the Range History of Moon-Base SAR

The SAR imaging mechanism is quite different from passive optic sensors, but the relative motion patterns are the same. Therefore, we can build the moon-based SAR geometry from modifying the proposed geometric model.
For the sake of light time effect, the ‘stop-go’ hypothesis is no longer valid [18,52]. The SAR baseband signal of a point target echo after demodulation can be written as [39,53]:
s ( τ , t ) = A 0   ω r ( τ R p ( t ) c ) p a 2 ( ϑ ( t ) ) exp ( j 2 π f 0 R p ( t ) c + j π K r ( τ R p ( t ) c ) 2 )
In this equation, τ denotes the fast time taking the center of each pulse as the origin. t denotes the slow time referring to UTC. A 0 is a complex constant related to the properties of the point target. The range window ω r refers to the pulse envelop. The one-way azimuth window p a ( θ ) sin c   ( θ L a f 0 / c ) , where L a is the antenna length. ϑ ( t ) is the angle between the wave and the normal plain of the antenna in range direction. f 0 is the carrier frequency and K r is the chirp rate of the transmitting signal. To simulate the signal, the propagation range history R p ( t ) must be accurately determined. Figure 12 shows the round-trip propagation geometry of an orbit SAR system. In a SAR system, the transmitting time t can be recorded, while the ‘go and back’ propagation time Δ t 1 and Δ t 2 is unknow. Considering the light time effects, we assume that the pulse center leaves the antenna S at t , arrives at the target A at t + Δ t 1 , and reaches the antenna S again at t + Δ t 1 + Δ t 2 . Then the propagation range is:
R p ( t ) = c   Δ t 1 + c   Δ t 2
where
c   Δ t 1 = S G C R S ( t ) A G C R S ( t + Δ t 1 ) c Δ t 2 = A G C R S ( t + Δ t 1 ) S G C R S ( t + Δ t 1 + Δ t 2 )
An iteration method or an approximation similar to Section 3.7 can be used to numerically solve Equation (31). The numerical result from Equations (29)–(31) gives an accurate signal simulation which can be used to develop fast focusing methods.
In SAR footprint geolocation, we need two targets to determine the three pointing angles. The first constrains the LOS, and the second describes the third pointing angle in Equation (5). For example, if the long side of the antenna is required to approximately parallel to the Earth’s longitude, the constraint condition can be written as:
A G C R S B G C R S = M S C S _ G C R S ( x ( 0 0 1 ) + y ( 1 0 0 ) )
where A and B indicate two targets with same longitude; x , y , and α R in M S C S _ O L C S are the three unknown numbers.
The aberration effect will result in a small offset angle in both range and azimuth. Theoretically, the transmitting wave and receiving wave should be translated to SCS to calculate the transmitting and receiving pattern weights. However, concerning the typical aberration angle is only 3.4 μrad, the difference is relatively small. For example, assuming the wavelength is 0.03 m (X-band) and the antenna length is 60 m, the round-trip weight difference is less than 3% in the half power beam width.
The angular aberration effect almost gives no impact on range history, but the light time will cause slight changes in both transmitting and receiving propagation range. The real range history should be the sum of (I) transmitting S ( t ) A ( t + Δ t 1 ) and (II) receiving A ( t + Δ t 1 ) S ( t + Δ t 1 + Δ t 2 ) (Figure 12). The reference range history is the sum of (III) reference transmitting range S ( t ) A ( t ) and (IV) receiving range S ( t + Δ t 1 ) A ( t + Δ t 1 ) used for approximation [18]. In airborne SAR and low earth orbit SAR imaging, the range history is approximated to (III) + (III). In reference [18], it is approximated to (III) + (IV). The differences between (III) + (III) and (I) + (II) reach 44.41 m, 34.12 m, and 7.90 m within a 10-min synthetic aperture when the target is put from low latitude to high latitude places (Figure 13). The variation of difference between (III) + (IV) and (I) + (II) are 22.21 m, 17.06 m, and 3.95 m. Though the residual error reduced by half, it is still greater than the signal wavelength during the long synthetic aperture time (SAT). The detrended results suggest that the variation of residual errors can be described as a linear term (Figure 13d). The high order error is within 1 mm in 200 s SAT and 3 mm in 600s SAT around the time of closest range, at least one order smaller than the wavelength of L, C, and X band. That means though the approximation (III) + (III) or (III) + (IV) is not accurate enough, obvious defocusing will not arise. The linear errors can only cause a mismatch on focusing position.

7. Conclusions

The correspondence between the intersection point and sensor attitude is given by the proposed geometric model. The simulated disk view is in good accordance with the image from the Chang’E lander and STK output. The analysis of the forecasted sensor’s attitude suggests that most places at the near side of the Moon are geometrically suitable for Earth observation, except for the marginal areas of the lunar disk and the neighborhood around the seleno-graphic origin. An observatory at the edge of the disk will likely experience view blocking by the topography or even by the lunar sphere itself, while an observatory near the lunar disk center will have a relatively high frequency of operation. The periodical approach between the Earth and the Sun on the celestial sphere requires a small FOV to prevent direct radiation. The analyses show that the areas near Cater Krafft have relatively low possibility of direct sunlit.
The light aberration scale is approximately 3.4 × 10−6 rad, corresponding to a deviation of 1.3 km on the Earth’s surface, while the light time scale is approximately 0.6 km in the opposite direction. As a result of neutralization, the combined effect is approximately 0.7 km, which is in the same order of magnitude as the pointing error. The light time also causes an error greater that the wavelength in moon-based SAR range history. In a 600s SAT around the time of closest range, the linear error can be as large as 44 m, but the high order error is at least one order smaller than the wavelength. The overall model accuracy is evaluated to be better than 4″ at current stage and will increase to be better than 1.5″ when the lunar positioning error is restricted to the level of the LLR reflector seleno-location accuracy.

Author Contributions

Conceptualization, Y.D., H.G. and G.L.; Data curation, Z.R.; Formal analysis, C.H.; Funding acquisition, H.G. and Y.D.; Investigation, G.L.; Methodology, Y.D.; Project administration, H.G.; Resources, H.S.; Software, Y.D.; Supervision, H.G.; Validation, M.L.; Visualization, H.S. and M.L.; Writing—original draft, Y.D.; Writing—review & editing, Y.D., C.H. and Z.R.

Funding

We thank the National Natural Science Foundation of China [grant numbers 41590853, 41501403] and the Key Project of Frontier Science Research of Chinese Academy of Sciences [QYZDY-SSW-DQC026] for their support.

Acknowledgments

Thanks Yingzhuo Jia for providing the terrain camera photos from Chang’E-3. Software Routines from the IAU SOFA Collection were used. Copyright © International Astronomical Union Standards of Fundamental Astronomy (http://www.iausofa.org). We also thank the JPL for the free access to the DE430 and the IERS for the free access to EOP data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Woerner, J.; Foing, B. The “Moon Village” Concept and Initiative. In Annual Meeting of the Lunar Exploration Analysis Group; USRA: Columbia, MD, USA, 2016. [Google Scholar]
  2. NASA. What Is Artemis? Available online: https://www.nasa.gov/feature/what-is-artemis/ (accessed on 29 June 2019).
  3. Ouyang, Z. Scientific objectives of Chinese lunar exploration project and development strategy. Adv. Earth Sci. 2004, 19, 351–358. [Google Scholar]
  4. Guo, H.; Liu, G.; Ding, Y. Moon-based Earth observation: Scientific concept and potential applications. Int. J. Digit. Earth 2018, 11, 546–557. [Google Scholar] [CrossRef]
  5. Xu, L.; Zou, Y.L.; Jia, Y.Z. China’s Planning for Deep Space Exploration and Lunar Exploration before 2030. Chin. J. Space Sci. 2018, 38, 11–12. [Google Scholar]
  6. Carruthers, G.R.; Page, T. Apollo 16 far-ultraviolet camera/spectrograph: Earth observations. Science 1972, 177, 788–791. [Google Scholar] [CrossRef]
  7. Feng, J.; Liu, J.; He, F.; Yan, W.; Ren, X.; Tan, X.; He, L.; Chen, B.; Zuo, W.; Wen, W.; et al. Data processing and initial results from the CE-3 Extreme Ultraviolet Camera. Res. Astron. Astrophys. 2014, 14, 1664. [Google Scholar] [CrossRef]
  8. Fornaro, G.; Franceschetti, G.; Lombardini, F.; Mori, A.; Calamia, M. Potentials and limitations of Moon-borne SAR imaging. IEEE Trans. Geosci. Remote Sens. 2010, 48, 3009–3019. [Google Scholar] [CrossRef]
  9. Moccia, A.; Renga, A. Synthetic aperture radar for earth observation from a lunar base: Performance and potential applications. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 1034–1051. [Google Scholar] [CrossRef]
  10. Ding, Y.; Guo, H.; Liu, G. Method to estimate the Doppler parameters of moon-borne SAR using JPL ephemeris. J. Beijing Univ. Aeronaut. Astronaut. 2015, 41, 71–76. [Google Scholar]
  11. Guo, H.; Ren, Y.; Liu, G.; Ye, H. The angular characteristics of Moon-based Earth observations. Int. J. Digit. Earth 2019, 1–16. [Google Scholar] [CrossRef]
  12. Ye, H.; Guo, H.; Liu, G.; Ren, Y. Observation duration analysis for Earth surface features from a Moon-based platform. Adv. Space Res. 2018, 62, 274–287. [Google Scholar] [CrossRef]
  13. Duan, W.; Huang, S.; Nie, C. Entrance Pupil Irradiance Estimating Model for a Moon-Based Earth Radiation Observatory Instrument. Remote Sens. 2019, 11, 583. [Google Scholar] [CrossRef]
  14. Nie, C.; Liao, J.; Shen, G.; Duan, W. Simulation of the land surface temperature from moon-based Earth observations. Adv. Space Res. 2019, 63, 826–839. [Google Scholar] [CrossRef]
  15. Duan, W.; Huang, S.; Nie, C. Conceptual design of a Moon-based Earth radiation observatory. Int. J. Remote Sens. 2018, 39, 5834–5849. [Google Scholar] [CrossRef]
  16. Xu, Z.; Chen, K.S. Effects of the Earth’s Curvature and Lunar Revolution on the Imaging Performance of the Moon-Based Synthetic Aperture Radar. IEEE Trans. Geosci. Remote Sens. 2019. [Google Scholar] [CrossRef]
  17. Xu, Z.; Chen, K.S. Temporal-Spatial Varying Background Ionospheric Effects on the Moon-Based Synthetic Aperture Radar Imaging: A Theoretical Analysis. IEEE Access 2018, 6, 66767–66786. [Google Scholar] [CrossRef]
  18. Xu, Z.; Chen, K.S. On signal modeling of moon-based synthetic aperture radar (SAR) imaging of earth. Remote Sens. 2018, 10, 486. [Google Scholar] [CrossRef]
  19. NASA. Lunar Exploration Objectives, Version 1. Available online: http://www.lpi.usra.edu/lunar/strategies/objectives/ex_ob_2006.pdf (accessed on 31 December 2006).
  20. Hamill, P. A Lunar Earth Observatory. Available online: https://www.lpi.usra.edu/meetings/LEA/whitepapers/HamillWhitePaper2.pdf (accessed on 2 March 2007).
  21. Pallé, E.; Goode, P.R. The Lunar Terrestrial Observatory: Observing the Earth using photometers on the Moon’s surface. Adv. Space Res. 2009, 43, 1083–1089. [Google Scholar] [CrossRef]
  22. Guo, H.; Ding, Y.; Liu, G.; Zhang, D.; Fu, W.; Zhang, L. Conceptual study of lunar-based SAR for global change monitoring. Sci. China Earth Sci. 2014, 57, 1771–1779. [Google Scholar] [CrossRef]
  23. Williams, J.G.; Newhall, X.X.; Dickey, J.O. Lunar moments, tides, orientation, and coordinate frames. Planet. Space Sci. 1996, 44, 1077–1080. [Google Scholar] [CrossRef]
  24. Williams, J.G.; Boggs, D.H.; Folkner, W.M. DE430 Lunar Orbit, Physical Librations, and Surface Coordinates; Interoffice Memorandum; Jet Propulsion Laboratory, California Institute of Technology: Pasadena, CA, USA, 2013. [Google Scholar]
  25. Archinal, B.A.; A’Hearn, M.F.; Bowell, E.; Conrad, A.; Consolmagno, G.J.; Courtin, R.; Fukushima, T.; Hestroffer, D.; Hilton, J.L.; Krasinsky, G.A.; et al. Report of the IAU working group on cartographic coordinates and rotational elements: 2009. Celest. Mech. Dyn. Astron. 2011, 109, 101–135. [Google Scholar] [CrossRef]
  26. Seidelmann, P.K.; Fukushima, T. Why new time scales? Astron. Astrophys. 1992, 265, 833–838. [Google Scholar]
  27. Montenbruck, O.; Gill, E. Satellite Orbits: Models, Methods and Applications; Springer: Berlin/Heidelberg, Germany, 2000. [Google Scholar]
  28. Adelberger, E.G.; Battat, J.B.R.; Birkmeier, K.J.; Colmenares, N.R.; Davis, R.; Hoyle, C.D.; Huang, R.L.; McMillan, R.J.; Murphy, T.W.; Schlerman, E.; et al. An absolute calibration system for millimeter-accuracy APOLLO measurements. Class. Quantum Gravity 2017, 34, 245008. [Google Scholar] [CrossRef] [Green Version]
  29. Petit, G.; Luzum, B. (Eds.) IERS Conventions (2010); Verlag des Bundesamts für Kartographie und Geodäsie: Frankfurt am Main, Germany, 2010. [Google Scholar]
  30. Li, Z.; Chen, B. Design and analysis of lunar-based tracking turntable for earth observation. Acta Photonica Sin. 2010, 39, 2215–2219. [Google Scholar]
  31. Hu, Q. Telescope Technology and Celestial Fixing; Southeast University Press: Nanjing, China, 2014. [Google Scholar]
  32. Folkner, W.M.; Williams, J.G.; Boggs, D.H.; Park, R.S.; Kuchynka, P. The planetary and lunar ephemerides DE430 and DE431. Interplanet. Netw. Prog. Rep. 2014, 196, 1–81. [Google Scholar]
  33. Greslou, D.; De Lussy, F.; Montel, J. Light aberration effect in HR geometric model. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2008, XXXVII-B1, 859–864. [Google Scholar]
  34. Lussy, F.; Greslou, D.; Dechoz, C.; Amberg, V.; Delvit, J.M.; Lebegue, L.; Blanchet, G.; Fourest, S. Pleiades HR in flight geometrical calibration: Location and mapping of the focal plane. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2012, XXXIX-B1, 519–523. [Google Scholar] [CrossRef]
  35. Ding, Y.; Guo, H.; Dong, J.; Liu, G.; Han, C.; Ruan, Z. Aberration effects in orbital imaging. Remote Sens. Lett. 2019, 10, 816–825. [Google Scholar] [CrossRef]
  36. Eckhardt, D.H. Theory of the libration of the moon. Moon Planets 1981, 25, 3–49. [Google Scholar] [CrossRef]
  37. Newhall, X.X.; Williams, J.G. Estimation of the lunar physical librations. Celest. Mech. Dyn. Astron. 1996, 66, 21–30. [Google Scholar] [CrossRef]
  38. Huang, Y.; Chang, S.; Li, P.; Hu, X.; Wang, G.; Liu, Q.; Zheng, W.; Fan, M. Orbit determination of Chang’E-3 and positioning of the lander and the rover. Chin. Sci. Bull. 2014, 59, 3858–3867. [Google Scholar] [CrossRef]
  39. Davies, M.E.; Colvin, T.R. Lunar coordinates in the regions of the Apollo landers. J. Geophys. Res. Planets 2000, 105, 20277–20280. [Google Scholar] [CrossRef]
  40. Archinal, B.A.; Rosiek, M.R.; Kirk, R.L.; Redding, B.L. The Unified Lunar Control Network 2005. Available online: https://pubs.usgs.gov/of/2006/1367/ (accessed on 20 December 2006).
  41. Di, K.; Liu, B.; Peng, M.; Xin, X.; Jia, M.; Zuo, W.; Ping, J.; Wu, B.; Oberst, J. An Initiative for Construction of New-Generation Lunar Global Control Network Using Multi-Mission Data. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2017, XLII-3/W1, 29–34. [Google Scholar]
  42. Mazarico, E.; Rowlands, D.D.; Neumann, G.A.; Smith, D.E.; Torrence, M.H.; Lemoine, F.G.; Zuber, M.T. Orbit determination of the lunar reconnaissance orbiter. J. Geod. 2012, 86, 193–207. [Google Scholar] [CrossRef]
  43. Lewis, H.; Lupton, W.; Sirota, M.J.; Mast, T.S.; Nelson, J.E.; Wallace, P.T. Pointing and tracking performance of the WM Keck Telescope. Proc. SPIE 1994, 2199, 117–126. [Google Scholar]
  44. Liu, J.; Kuga, Y.; Ishimaru, A.; Pi, X.; Freeman, A. Ionospheric effects on SAR imaging: A numerical study. IEEE Trans. Geosci. Remote Sens. 2003, 41, 939–947. [Google Scholar]
  45. Wallace, P.T.; Capitaine, N. Precession-nutation procedures consistent with IAU 2006 resolutions. Astron. Astrophys. 2006, 459, 981–985. [Google Scholar] [CrossRef] [Green Version]
  46. Mukherjee, S.; Joshi, P.K.; Mukherjee, S.; Ghosh, A.; Grag, R.D.; Mukhopadhyay, A. Evaluation of vertical accuracy of open source Digital Elevation Model (DEM). Int. J. Appl. Earth Obs. Geoinf. 2013, 21, 205–217. [Google Scholar] [CrossRef]
  47. Rodriguez, E.; Morris, C.S.; Belz, J.E. A global assessment of the SRTM performance. Photogramm. Eng. Remote Sens. 2006, 72, 249–260. [Google Scholar] [CrossRef]
  48. Burke, J.D. Merits of a lunar polar base location. In Lunar Bases and Space Activities of the 21st Century; Lunar & Planetary Institute: Houston, TX, USA, 1985. [Google Scholar]
  49. Speyerer, E.J.; Robinson, M.S. Persistently illuminated regions at the lunar poles: Ideal sites for future exploration. Icarus 2013, 222, 122–136. [Google Scholar] [CrossRef]
  50. Mazarico, E.; Neumann, G.A.; Smith, D.E.; Zuber, M.T.; Torrence, M.H. Illumination conditions of the lunar polar regions using LOLA topography. Icarus 2011, 211, 1066–1081. [Google Scholar] [CrossRef]
  51. Jia, Y.; Zou, Y. Research on lunar site selection for lunar based earth observation. Spacecr. Eng. 2016, 25, 116–121. [Google Scholar]
  52. Hu, C.; Long, T.; Zeng, T.; Liu, F.; Liu, Z. The accurate focusing and resolution analysis method in geosynchronous SAR. IEEE Trans. Geosci. Remote Sens. 2011, 49, 3548–3563. [Google Scholar] [CrossRef]
  53. Cumming, I.G.; Wong, F.H. Digital Processing of Synthetic Aperture Radar Data; Artech house: Norwood, MA, USA, 2005. [Google Scholar]
Figure 1. Coordinate transforms on the lunar sphere (the difference between the mean-Earth/rotation-axis (MER) and principle axis (PA) are not distinguished; thus, they are both presented by the Selenocentric Body-fixed Coordinate System (SBCS) in the figure).
Figure 1. Coordinate transforms on the lunar sphere (the difference between the mean-Earth/rotation-axis (MER) and principle axis (PA) are not distinguished; thus, they are both presented by the Selenocentric Body-fixed Coordinate System (SBCS) in the figure).
Remotesensing 11 02611 g001
Figure 2. The coordinate transforms on the celestial sphere.
Figure 2. The coordinate transforms on the celestial sphere.
Remotesensing 11 02611 g002
Figure 3. The spatial relation in the Geocentric Celestial Reference System (GCRS).
Figure 3. The spatial relation in the Geocentric Celestial Reference System (GCRS).
Remotesensing 11 02611 g003
Figure 4. Simulation flow chart.
Figure 4. Simulation flow chart.
Remotesensing 11 02611 g004
Figure 5. Infrared global geostationary composite v1 at UTC 2013/12/24 18:15 (from EarthData, https://earthdata.nasa.gov/).
Figure 5. Infrared global geostationary composite v1 at UTC 2013/12/24 18:15 (from EarthData, https://earthdata.nasa.gov/).
Remotesensing 11 02611 g005
Figure 6. Geo-distances between the calculated results without aberration and light time correction and the estimated results of the satellite tool kit (STK).
Figure 6. Geo-distances between the calculated results without aberration and light time correction and the estimated results of the satellite tool kit (STK).
Remotesensing 11 02611 g006
Figure 7. Aberration and light time effects.
Figure 7. Aberration and light time effects.
Remotesensing 11 02611 g007
Figure 8. Pointing errors caused by the (a) attitude error, (b) libration error, (c) seleno-graphic error with a magnitude of one second, and (d) lunar center error and elevation error with a magnitude of 10 m. The sensor is located at (−60,0) and is pointed to the geocenter.
Figure 8. Pointing errors caused by the (a) attitude error, (b) libration error, (c) seleno-graphic error with a magnitude of one second, and (d) lunar center error and elevation error with a magnitude of 10 m. The sensor is located at (−60,0) and is pointed to the geocenter.
Remotesensing 11 02611 g008aRemotesensing 11 02611 g008b
Figure 9. There are seven categories of error sources analyzed. The errors from aberration and light time have been corrected, and thus are not included.
Figure 9. There are seven categories of error sources analyzed. The errors from aberration and light time have been corrected, and thus are not included.
Remotesensing 11 02611 g009
Figure 10. Attitude of the sensor.
Figure 10. Attitude of the sensor.
Remotesensing 11 02611 g010
Figure 11. Four observational parameters put constraints on the observatory site selection: (a) the percentage of full coverage, (b) the mean derivative of elevation angle (arcsecond/hour), (c) the mean derivative of azimuth angle (arcsecond/hour), and (d) the number of times when LOS and the negative light direction toward the sensor (LOS-SUN) angle is smaller than FOV/2. The simulation step length is one minute from 2000/1/1 to 2019/1/1. The x and y coordinates indicate seleno-longitude and seleno-latitude, respectively.
Figure 11. Four observational parameters put constraints on the observatory site selection: (a) the percentage of full coverage, (b) the mean derivative of elevation angle (arcsecond/hour), (c) the mean derivative of azimuth angle (arcsecond/hour), and (d) the number of times when LOS and the negative light direction toward the sensor (LOS-SUN) angle is smaller than FOV/2. The simulation step length is one minute from 2000/1/1 to 2019/1/1. The x and y coordinates indicate seleno-longitude and seleno-latitude, respectively.
Remotesensing 11 02611 g011
Figure 12. The moon-based synthetic aperture radar (SAR) geometric model in the GCRS. The dashed line indicates the round-trip propagation range.
Figure 12. The moon-based synthetic aperture radar (SAR) geometric model in the GCRS. The dashed line indicates the round-trip propagation range.
Remotesensing 11 02611 g012
Figure 13. The range difference of (II), (III), and (IV) referred to (I) when the SAR is at ( 80 , 0 ) on lunar surface and the target is at (a) ( 0 , 0 ) , (b) ( 40 , 0 ) , and (c) ( 80 , 0 ) on Earth’s surface. The synthetic aperture center is calibrated to the time of closest range at 2017/10/9. Subfigure (d) indicates the detrended result of subfigure (a).
Figure 13. The range difference of (II), (III), and (IV) referred to (I) when the SAR is at ( 80 , 0 ) on lunar surface and the target is at (a) ( 0 , 0 ) , (b) ( 40 , 0 ) , and (c) ( 80 , 0 ) on Earth’s surface. The synthetic aperture center is calibrated to the time of closest range at 2017/10/9. Subfigure (d) indicates the detrended result of subfigure (a).
Remotesensing 11 02611 g013
Table 1. Data comparison.
Table 1. Data comparison.
UTC+8Chang’E-3 LanderSTK This Paper
2013/12/21
02:46:12
Remotesensing 11 02611 i001 Remotesensing 11 02611 i002 Remotesensing 11 02611 i003
2013/12/22
09:15:50
Remotesensing 11 02611 i004 Remotesensing 11 02611 i005 Remotesensing 11 02611 i006
2013/12/23
10:12:47
Remotesensing 11 02611 i007 Remotesensing 11 02611 i008 Remotesensing 11 02611 i009
+2013/12/25
02:19:33
Remotesensing 11 02611 i010 Remotesensing 11 02611 i011 Remotesensing 11 02611 i012
Table 2. The max and min geo-distance deviation (meters) derived from one second sensor attitude or seleno-graphic coordinate errors over the past 50 years when the line of sight (LOS) is pointed to the geocenter. The simulation step size is one day. The deviations caused by the librations and lunar center coordinates given by DE430 are consistent with Figure 8 and are therefore not shown in this Table.
Table 2. The max and min geo-distance deviation (meters) derived from one second sensor attitude or seleno-graphic coordinate errors over the past 50 years when the line of sight (LOS) is pointed to the geocenter. The simulation step size is one day. The deviations caused by the librations and lunar center coordinates given by DE430 are consistent with Figure 8 and are therefore not shown in this Table.
Location α I θ I ξ M M H (10 m)
( 0 , 0 ) 4~3331689~19321698~19411687~19400~1.8
( 40 , 0 ) 939~14111691~19351698~19411687~19405.5~7.4
( 80 , 0 ) 1630~19371695~19401698–19411687–19409.6–10
( 0 , 40 ) 955–13881691–19341204–15871687–19405.3–7.5
( 0 , 80 ) 1662–19211696–193975–5991687–19409.5–10
Table 3. Estimated accuracies.
Table 3. Estimated accuracies.
CategorySub-CategoryAngular Accuracy
Lunar Orbit and Libration 0.1″
Seleno-graphicLLR0.032″
LRO LOLA DEM2.4″
Pointing 1″
PropagationRefractive0.087″
EOP 0.016″
GeographicGeolocation-Height0.015″
Others 0.1″

Share and Cite

MDPI and ACS Style

Ding, Y.; Guo, H.; Liu, G.; Han, C.; Shang, H.; Ruan, Z.; Lv, M. Constructing a High-Accuracy Geometric Model for Moon-Based Earth Observation. Remote Sens. 2019, 11, 2611. https://doi.org/10.3390/rs11222611

AMA Style

Ding Y, Guo H, Liu G, Han C, Shang H, Ruan Z, Lv M. Constructing a High-Accuracy Geometric Model for Moon-Based Earth Observation. Remote Sensing. 2019; 11(22):2611. https://doi.org/10.3390/rs11222611

Chicago/Turabian Style

Ding, Yixing, Huadong Guo, Guang Liu, Chunming Han, Haolv Shang, Zhixing Ruan, and Mingyang Lv. 2019. "Constructing a High-Accuracy Geometric Model for Moon-Based Earth Observation" Remote Sensing 11, no. 22: 2611. https://doi.org/10.3390/rs11222611

APA Style

Ding, Y., Guo, H., Liu, G., Han, C., Shang, H., Ruan, Z., & Lv, M. (2019). Constructing a High-Accuracy Geometric Model for Moon-Based Earth Observation. Remote Sensing, 11(22), 2611. https://doi.org/10.3390/rs11222611

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop