Theoretical Feasibility Analysis of Fast Back-Projection Algorithm for Moon-Based SAR in Time Domain

: Nowadays, the Earth observation based on the Moon has attracted attention from many researchers and relevant departments. There also exists a considerable amount of interest in monitor-ing large-scale and long-term geoscience phenomena using the Moon-based SAR (MBS). However, the Earth’s observation from MBS has long transmission time, and the relative motion of MBS with its Earth ground target (EGT) is much different to the space-borne SAR, the above reasons indicate that the traditional stop-and-go model is no longer suitable for MBS in frequency domain imaging. Here a dual-path separate calculation method for single pulse is presented in this paper for a better match of a real scenario, and then the slant range is ﬁtted to a high-order polynomial series. The MBS’s location, the synthetic aperture time and other factors have effects on length of the dual-path and ﬁt bias. Without thoroughly investigated phase de-correlation processing in frequency domain, and to avoid computational costs in traditional back-projection (BP) algorithm, the paper ﬁrst proposes a fast back-projection (FBP) algorithm in time domain for MBS, a platform that has long transmission time and long synthetic aperture time. In the FBP algorithm, the original method, that projected echo on all pixels in the imaging area, is changed to projected echo on a centerline instead. A suitable interpolation for points on the centerline is adopted to reduce the projected error; the synthetic aperture length and imaging area are also divided into subsections to reduce computation cost. The formula indicates that the range error could be control once the product of sub-imaging area’s length and sub-aperture’s length stay constant. Through the theoretical analysis, the detailed range difference mainly at apogee, perigee, ascending, and descending nodes indicate the necessity to separately calculate the dual-path for MBS’s single pulse transmission in Earth-Moon motion, with real ephemeris been adopted; then, the high-order polynomial ﬁtting will better describe the motion trajectory. Lastly, the FBP algorithm proposed is simulated in a speciﬁc scenario under acceptable resolution, and the result shows its feasibility for image compression.


Introduction
In recent years, using the Moon as a platform for remote sensing monitoring, (e.g., the global energy balance, the ocean vortex, the earth's surface tide, etc.) has attract increasing attention across different research fields [1,2]. As the only natural satellite of Earth, the Moon's unique rotation and revolution laws make it as a stable platform for the Earth observation with long-term and large-scale characteristics [3,4], which has an advantage over the space-born SAR platform [5][6][7]. Accordingly, the MBS would be a reliable supplement to the existing remote sensing system [8][9][10].
However, the MBS has several unique features different from space-borne SAR platform. First, the perigee, apogee, and the average distance in a sidereal month between The feasibility of ground-based imaging radar technology on Moon might indicate the mission is reasonable. The size, power, etc. of MBS can use the same scale as that used on Earth. With a primary focus on the renewed algorithm for MBS, the influence of attenuation in signal transmission, the light aberration, and also MBS's and EGT's digital elevation mode (DEM) on imaging are not considered. The remainder of this paper will prove the effectiveness of MBS from two aspects, i.e., the spatial relationship and the establishment of the algorithm model; the framework is organized as below. In Section 2, a description of the FBP algorithm, mainly regarding the division of aperture and imaging areas for MBS, is presented based on real motion between MBS and EGT. The real motion, including the range history, is calculated separately for uplink and downlink, based on ephemeris. In Section 3, a theoretical analysis of the image is provided, and the target responses are simulated for confirming the superiority of FBP. Finally, Section 4 offers a brief conclusion for the simulation.

Algorithm for MBS
Unlike the space-born SAR, the imaging of MBS is much different because of the complex Earth-Moon motion relationship, as well as atmospheric and ionosphere disturbances. The FBP algorithm in this paper mainly focus on the spatial relationship and the establishment of the algorithm model, as seen in the detailed flow chart displayed in Figure 2. Box 1 shows the processing of echo generated from the range history between MBS and EGT; since the Moon revolves around Earth with Earth's rotation, the ex-poser time for MBS' achievable EGT changes significantly. Here, the range history is calculated based on a real scenario, e.g., the JPL ephemeris and EOP at different time epochs, and other MBS parameters, such as carrier frequency, bandwidth, and look angle (here, an angle between a vector, i.e., from MBS to EGT, and another vector, i.e., from MBS to Earth centroid, is calculated as the look angle) mainly contribute to the beam shape on Earth, which is related to the MBS's resolution. The detailed introduction of how to calculate the range history is displayed in Section 2.2. Box 2 shows the progression of the final image generated from the echo in Box 1. The sub-aperture division and sub-imaging area are the keys to the FBP algorithm. In the traditional BP algorithm, the sampling echo is projected on all pixels of the image area; this coherent method incurs huge computational time costs, especially when sampling points are very large. Therefore, for reducing time consumption The feasibility of ground-based imaging radar technology on Moon might indicate the mission is reasonable. The size, power, etc. of MBS can use the same scale as that used on Earth. With a primary focus on the renewed algorithm for MBS, the influence of attenuation in signal transmission, the light aberration, and also MBS's and EGT's digital elevation mode (DEM) on imaging are not considered. The remainder of this paper will prove the effectiveness of MBS from two aspects, i.e., the spatial relationship and the establishment of the algorithm model; the framework is organized as below. In Section 2, a description of the FBP algorithm, mainly regarding the division of aperture and imaging areas for MBS, is presented based on real motion between MBS and EGT. The real motion, including the range history, is calculated separately for uplink and downlink, based on ephemeris. In Section 3, a theoretical analysis of the image is provided, and the target responses are simulated for confirming the superiority of FBP. Finally, Section 4 offers a brief conclusion for the simulation.

Algorithm for MBS
Unlike the space-born SAR, the imaging of MBS is much different because of the complex Earth-Moon motion relationship, as well as atmospheric and ionosphere disturbances. The FBP algorithm in this paper mainly focus on the spatial relationship and the establishment of the algorithm model, as seen in the detailed flow chart displayed in Figure 2.
Box 1 shows the processing of echo generated from the range history between MBS and EGT; since the Moon revolves around Earth with Earth's rotation, the ex-poser time for MBS' achievable EGT changes significantly. Here, the range history is calculated based on a real scenario, e.g., the JPL ephemeris and EOP at different time epochs, and other MBS parameters, such as carrier frequency, bandwidth, and look angle (here, an angle between a vector, i.e., from MBS to EGT, and another vector, i.e., from MBS to Earth centroid, is calculated as the look angle) mainly contribute to the beam shape on Earth, which is related to the MBS's resolution. The detailed introduction of how to calculate the range history is displayed in Section 2.2. Box 2 shows the progression of the final image generated from the echo in Box 1. The sub-aperture division and sub-imaging area are the keys to the FBP algorithm. In the traditional BP algorithm, the sampling echo is projected on all pixels of the image area; this coherent method incurs huge computational time costs, especially when sampling points are very large. Therefore, for reducing time consumption without a loss of generality, the echo projection method used in FBP is changed to a centerline based on the sub-aperture and sub-imaging area division within acceptable error; the range interpolation on the projected centerline contributes to a reduction of range error; the final image will be generated after the superposition processing. The detailed introduction of the core processing will be displayed in Section 2.1.
Appl. Sci. 2022, 12, x FOR PEER REVIEW 4 of 18 without a loss of generality, the echo projection method used in FBP is changed to a centerline based on the sub-aperture and sub-imaging area division within acceptable error; the range interpolation on the projected centerline contributes to a reduction of range error; the final image will be generated after the superposition processing. The detailed introduction of the core processing will be displayed in Section 2.1.

Division Mode in FBP
According to the imaging process of BP algorithm, the dual-delay time from SAR platform to each point in imaging area is calculated; then compensation phase is added to the echo corresponding to each point in time domain regarding the range compressing process, and the final imaging is received by coherently accumulating the projection data consistent with different azimuth moments [16,17]. After the range matching and filtering, the echo of EGT is expressed in Equation (1), as shown below: where λ is the wavelength of SAR platform; c is light speed in vacuum; τ is the fast time at range direction; R(η) is the slant range between MBS and its corresponding EGT in slow time η at azimuth direction; and A is the signal amplitude, which indicates the reflective properties of EGT and is often deemed as a constant of one for simplicity. Unlike the slant range calculated in a low orbit platform, its trajectory can be replaced with second-order hyperbola under a constant speed with little change, the curve motion and the locations' arbitrariness of MBS no longer defines trajectory and speed as platforms in a low orbit; a detail description will be provided in Section 2.2.
After the echo has been received, each point in imaging area is projected to the echo domain, forming the migration trajectory of EGT; then, the complex data of EGT is obtained after the echo been superimposed. The complex data of each pixel (i,k) in the imaging area can be expressed as shown below:

Division Mode in FBP
According to the imaging process of BP algorithm, the dual-delay time from SAR platform to each point in imaging area is calculated; then compensation phase is added to the echo corresponding to each point in time domain regarding the range compressing process, and the final imaging is received by coherently accumulating the projection data consistent with different azimuth moments [16,17]. After the range matching and filtering, the echo of EGT is expressed in Equation (1), as shown below: where λ is the wavelength of SAR platform; c is light speed in vacuum; τ is the fast time at range direction; R(η) is the slant range between MBS and its corresponding EGT in slow time η at azimuth direction; and A is the signal amplitude, which indicates the reflective properties of EGT and is often deemed as a constant of one for simplicity. Unlike the slant range calculated in a low orbit platform, its trajectory can be replaced with second-order hyperbola under a constant speed with little change, the curve motion and the locations' arbitrariness of MBS no longer defines trajectory and speed as platforms in a low orbit; a detail description will be provided in Section 2.2.
After the echo has been received, each point in imaging area is projected to the echo domain, forming the migration trajectory of EGT; then, the complex data of EGT is obtained after the echo been superimposed. The complex data of each pixel (i,k) in the imaging area can be expressed as shown below: where τ(η) = 2·R(η) c indicates the dual-delay time at the azimuth time η, t is the synthetic aperture time of MBS.
Since each point in image area has been projected on echo domain in the BP algorithm, if the sampling numbers in azimuth direction is N a , the pixel numbers along azimuth and range directions in imaging areas are M × N, the computational cost of the projection calculation in imaging area is N a × M × N; the cost of phase compensation is also N a × M × N, and the cost of coherent superposition is (N a − 1) × M × N. Therefore, the total computational cost is the sum of above three computational costs, i.e., (3N a − 1) × M × N. Once the number of sampling points becomes large enough, the computational cost will exhibit undesired exponential growth. However, for MBS located in the near side of Moon, its synthetic aperture time can be hundreds of seconds longer than that on the space-born platform; moreover, the imaging area around the Earth's surface can be larger, e.g., exceeding millions of square kilometers, under a small gaze angle. Therefore, the computational cost of the BP algorithm in the MBS scenario might be intolerable large. Therefore, to bring down the cost, the traditional BP algorithm is modified into the FBP algorithm, shown as follows, based on sub-aperture and sub-imaging area division.
In the proposed FBP algorithm, the synthetic aperture length and corresponding imaging area for MBS will be divided into sub-apertures and sub-imaging areas, as shown in Figure 3a. Here, the lengths in Figure 3a are not shown in the same ratio as the actual scenario, in order to create a clear display. The synthetic aperture length L s is divided into α 1 small non-overlapping sub-apertures; therefore, the actual length of every sub-aperture is d = L s α 1 . The size of the original imaging area, i.e., M × N (N is pixels along the range direction and M is the pixels along the azimuth direction), will be divided into α 2 × α 3 small sub-areas; the size of every sub-area is m × n, where m = M α 2 , m = N α 3 . Here, the dark area in Figure 3a represents the corresponding pixels in the sub-imaging area.
Appl. Sci. 2022, 12, x FOR PEER REVIEW 5 of 17 N × M × N, and the cost of coherent superposition is (N − 1) × M × N. Therefore, the total computational cost is the sum of above three computational costs, i.e., (3N − 1) × M × N. Once the number of sampling points becomes large enough, the computational cost will exhibit undesired exponential growth. However, for MBS located in the near side of Moon, its synthetic aperture time can be hundreds of seconds longer than that on the space-born platform; moreover, the imaging area around the Earth's surface can be larger, e.g., exceeding millions of square kilometers, under a small gaze angle. Therefore, the computational cost of the BP algorithm in the MBS scenario might be intolerable large. Therefore, to bring down the cost, the traditional BP algorithm is modified into the FBP algorithm, shown as follows, based on sub-aperture and sub-imaging area division.
In the proposed FBP algorithm, the synthetic aperture length and corresponding imaging area for MBS will be divided into sub-apertures and sub-imaging areas, as shown in Figure 3a. Here, the lengths in Figure 3a are not shown in the same ratio as the actual scenario, in order to create a clear display. The synthetic aperture length is divided into α 1 small non-overlapping sub-apertures; therefore, the actual length of every subaperture is = α 1 . The size of the original imaging area, i.e., M × N (N is pixels along the range direction and M is the pixels along the azimuth direction), will be divided into α 2 × α 3 small sub-areas; the size of every sub-area is . Here, the dark area in Figure 3a represents the corresponding pixels in the sub-imaging area.
(a) (b) After the sub-aperture and sub-imaging area been divided, the subsequent task is to handle the complex data in each sub-object. As shown in Figure 3b, a group of concentric circles centered on sub-aperture's point can be obtained for lines from every sub-aperture to pixels in the sub-image area. Different circles form their respective sector areas and indicate unequal echo; other back-projection data of points on an identical circle are exactly the same. Thus, within an acceptable range error, it is not necessary to project echo data to all pixels in the imaging area, as the centerline represents all pixels distributed in the sub-image area. For two adjacent sub-apertures, the distance to pixels in the sub-image area can be almost the same. For the α(1 ≤ α ≤ α 1 ) aperture and its adjacent apertures, it can be considered that they exhibit approximately the same concentric circles. Here, the echo data in the sector areas can be considered to be consistent with the data of the line from the sub-aperture's center to the sub-image's area center. Therefore, a centerline can be used as an approximation for the sector areas within an allowable range of error, thus reducing the computational cost required by the FBP algorithm.
For the scenario of MBS, the scale change in slant range undulate insignificantly along with the gaze angle. Figure 3b shows the calculation of range error; the lines in this figure do not have the same scale as the actual scenario, in order to create a clear display. After the sub-aperture and sub-imaging area been divided, the subsequent task is to handle the complex data in each sub-object. As shown in Figure 3b, a group of concentric circles centered on sub-aperture's point can be obtained for lines from every sub-aperture to pixels in the sub-image area. Different circles form their respective sector areas and indicate unequal echo; other back-projection data of points on an identical circle are exactly the same. Thus, within an acceptable range error, it is not necessary to project echo data to all pixels in the imaging area, as the centerline represents all pixels distributed in the sub-image area. For two adjacent sub-apertures, the distance to pixels in the sub-image area can be almost the same. For the α(1 ≤ α ≤ α 1 ) aperture and its adjacent apertures, it can be considered that they exhibit approximately the same concentric circles. Here, the echo data in the sector areas can be considered to be consistent with the data of the line from the sub-aperture's center to the sub-image's area center. Therefore, a centerline can be used as an approximation for the sector areas within an allowable range of error, thus reducing the computational cost required by the FBP algorithm.
For the scenario of MBS, the scale change in slant range undulate insignificantly along with the gaze angle. Figure  is the corresponding projected position; and point E is the intersection between segment AD and equal radius arc L 1 . According to the approximation principle, the points near the centerline use points on the centerline instead, although the approximate processing will thus produce some errors and affect the imaging quality. In this paper, the centerline is generally selected, without loss of generality, from the center of the sub-aperture to the center of the sub-imaging area, but only segments of the centerline are used. On one hand, if the segment is too short, all the pixels in the sub-imaging areas cannot be effectively projected on the line, undoubtedly increasing the range error before and after projection. On the other hand, if the segment is too long, much redundant data will be generated, in opposition to the requirement of less intermediate data. In this paper, the segment equal to the diagonal of the sub-image area is selected to ensure that each pixel in the sub-image area can be projected without generating too much redundant data. After the pixels in the sub-image area have been replaced by points on the centerline, e.g., points in Figure 3b, the range error between the real range AD and the replaced BC on the centerline can be expressed as follows: where θ is angle between the azimuth direction and point E; ∆θ is the angle between the centerline and point E; R is the true slant range (i.e., the length of AC), and ∆R is the slant range error between length of AC and BC (i.e., the error between the true range and the projected range). Because the range error ∆R(θ) is far less than any range from MBS, to achieve EGT, the range error ∆R can be approximately quantified by the equation shown below: According to the cosines law of triangles, the slant range R can be re-expressed as follows: where u is the length from position of the beam to its initial position, in this paper, the initial position is reckoned as the center of sub-aperture, r is the length of centerline, and d is the length of sub-aperture. Therefore, Equation (4) can be reverted as follow.
In actual scenario, u is much smaller than r. Therefore, considering the extreme situation, the maximum value of the range error in Equation (6) can be expressed as shown below: where D is the length of sub-image area in azimuth direction. Equation (7) indicates that the range error will fall within a limit, so long as the product, namely the sub-aperture's length and the sub-image area's length in azimuth direction, remains unchanged. It is generally assumed that the phase error caused by the range error should not be greater than π 4 ; this is to say, the maximum value of the range error should be subject to Equation (8).
4π The maximum range error in Equation (8) will often be less than λ 16 , and the denominator constant can be substituted with a control factor δ. With a combination of Equations (7) and (8), the sub-image area's length can be expressed as shown below: In Equation (9), the product of the sub-image area's length and sub-aperture's length has a negative correlation with the control factor δ. Once δ becomes too large, the subaperture's length will be closer to a pulse, and thus cannot achieve the effect of aperture division and computational cost reduction, so a suitable δ should make d equal to the length of the exponential power of 2, i.e., the sub-aperture's length contains 16, 32, or even 64 pulses.
With suitable parameters before being set, the true position showed in Figure 3b still might not project on an existing sampling point of the centerline, shown as a solid red dot; this situation will cause an error expressed as dr. In order to reduce such an error, a step of different multiples interpolation, shown in Box 2 of Figure 2, will be conducted during the imaging process. Theoretically, dr can be reduced by increasing the interpolation multiple; the higher the interpolation multiple, the smaller the error. However, dr will accumulate with the superposition of sub-apertures; the more sub-apertures, the greater the error accumulated, which may affect the imaging quality.
After the FBP algorithm has been introduced, the next step is to derive the range or phase history of the echo; both are influenced by Earth-Moon motion and the actual location of MBS. For a better description of the processing of echo, the subsequent section introduces a detailed method of range history estimation. In this method, the locations of the MBS during periodic Earth-Moon motion are all taken into consideration.

Range History Calculation
Since the dual-delay time of MBS can exceed 2.3 s, the transmitted and received positions of MBS are not the same; the stop-and-go model applied to a low orbit platform will no longer be applicable for MBS distributed around the near side of Moon. In an actual scenario, the long average distance from MBS to its achievable EGT, the rotation of Earth, the rotation of Moon, and the revolution of Moon, together with DEM of MBS and EGT will have different effects on their slant range. The near side of Moon is distributed with uneven craters, hills, uplifts, and potholes that have unequal DEM, and so does the EGT; the differences in locations of MBS and EGT bring diversity to the ground range and slant range, which have been considered in imaging processing; because of earth tide lock, the revolution and rotation of the Moon's near side periodically faces the Earth. The Earth's unequal rotation in latitude direction makes the exposure time, related to synthetic time and length, of EGT change unevenly. Without loss of generality and to simplify the analysis process, the diverse DEM of MBS and its corresponding EGT are omitted, and the two planets, i.e., the Earth and the Moon, are considered as follows: the Earth is considered to be a standard ellipsoid body, with an equator radius of about 6377 km, and a polar radius of about 6356 km, and the Moon is considered to be a standard sphere, with average radius of about 1737 km, respectively [18]. Then, consider the relative motion between MBS, e.g., its location is about 0 • , 0 • on the Moon's surface, and its corresponding look angle is about 0.1 • , as shown in Figure 4. The distance from MBS to EGT with Earth's rotation (an oscillating smooth curve displayed in the black solid line) shows much different when compared to the results without Earth's rotation (a sine smooth curve displayed in the green solid line). The curve occlusion of the sphere is also ignored to provide a simplified description. Although their overall trends are similar, the actual distance shows that oscillation changes daily, so the Earth's rotation cannot be ignored for actual Earth-Moon SAR imaging, especially when determining MBS and accessible EGT. For a more realistic scenario, the ephemeris are adopted among calculations in this paper.
Earth-Moon SAR imaging, especially when determining MBS and accessible EGT. For a more realistic scenario, the ephemeris are adopted among calculations in this paper. While adopting JPL DE430, or other ephemerides at specific moments for researching, the dual-path method is used, i.e., uplink and downlink components, are calculated separately, with Earth's rotation been considered in every pulse [19,20]. Here, its antenna bore-sight lies in the corresponding zero-Doppler plane after using at least the two-dimensional attitude steering method, and the plane is assumed to remain unchanged during synthetic aperture time. The dual-path method from MBS to achievable EGT at the pulse is shown as Figure 5, where 1 is delay time of the downlink path from MBS to EGT; 2 is delay time of the uplink from EGT to MBS; is pulse repeat time (PRT) of MBS; is number of pulses; , + 1 , + 1 + 2 are three adjacent discrete time variables, indicating the wave propagation of a single pulse; ,   While adopting JPL DE430, or other ephemerides at specific moments for researching, the dual-path method is used, i.e., uplink and downlink components, are calculated separately, with Earth's rotation been considered in every pulse [19,20]. Here, its antenna boresight lies in the corresponding zero-Doppler plane after using at least the two-dimensional attitude steering method, and the plane is assumed to remain unchanged during synthetic aperture time. The dual-path method from MBS to achievable EGT at the n pulse is shown as Figure 5, where τ 1 is delay time of the downlink path from MBS to EGT; τ 2 is delay time of the uplink from EGT to MBS; T is pulse repeat time (PRT) of MBS; n is number of pulses; nT, nT + τ 1 , nT + τ 1 + τ 2 are three adjacent discrete time variables, indicating the wave propagation of a single pulse; r nT t , r Earth-Moon SAR imaging, especially when determining MBS and accessible EGT. For a more realistic scenario, the ephemeris are adopted among calculations in this paper. While adopting JPL DE430, or other ephemerides at specific moments for researching, the dual-path method is used, i.e., uplink and downlink components, are calculated separately, with Earth's rotation been considered in every pulse [19,20]. Here, its antenna bore-sight lies in the corresponding zero-Doppler plane after using at least the two-dimensional attitude steering method, and the plane is assumed to remain unchanged during synthetic aperture time. The dual-path method from MBS to achievable EGT at the pulse is shown as Figure 5, where 1 is delay time of the downlink path from MBS to EGT; 2 is delay time of the uplink from EGT to MBS; is pulse repeat time (PRT) of MBS; is number of pulses; , + 1 , + 1 + 2 are three adjacent discrete time variables, indicating the wave propagation of a single pulse; ,     ; this process can be called an entire delay cycle for a single pulse; here the two velocities v nT t and v nT l are deemed to be unchanged in an entire delay cycle. The coordinate system for the calculation of the dual path involved with the Earth's central rotational coordinate system (ECR), the Earth's central inertial coordinate system (ECI), the Moon's central Moon fixed coordinate system (MCMF), etc. In this study, the position and velocity of MBS and EGT should be represented in the same coordinate system, for example, in the ECR.
From the cosine theorem of trigonometric functions, we can obtain the equations for the signal transmission path of downlink in Equation (10) and uplink in Equation (11), respectively: where the top right subscript "Tr " means matrix transpose, and the symbol "< , >" signifies inner product operator. From Equations (11) and (12), the two delay time variables can be expressed in Equations (12) and (13) as shown below: Then consequently, the total range history at pulse nT that consists of both downlink propagation and uplink propagation, can be expressed as below: If we convert the pulse time nT into azimuth slow time η, then Equation (14) can be transformed to a high-order polynomial series, as shown in Equation (15), to fit the range history and coherent superposition for imaging.
where k is the exponential powers corresponding to coefficientσ, and o η k is the higher order infinitesimal.

Theoretical Analysis
After the FBP and the real range calculation modes have been established, specific MBS on the Moon's surface will be selected for research using the application of FBP. Inspired by feasible ground-based imaging radar technology used on the Moon, and the requirement of computation cost reduction for the traditional BP algorithm, the proposed FBP algorithm aimed to ameliorate the cost. So in the first of this section, different locations of MBS are chosen for research on range history characteristic, and second, images of one specific MBS will be figured to show the performance status at four key moments, i.e., apogee, perigee, ascending node, and descending node, during one specific period.

Range Characteristic between MBS and EGT
Because the revolution speed of Moon's centroid (ranging from about 960 m/s to 1100 m/s) is much faster than its maximum rotation speed (about 4.6 m/s) in the equatorial region, the Moon's revolution plays a major role in slant range compared with its rotation.
Because the maximum rotation speed (about 465 m/s) of Earth in its equatorial region is not identical to Moon's revolution, and the rotation orientation of Earth as well as the Moon's revolution orientation are all from west to east, as observed from the north celestial pole, the delay time from MBS to corresponding EGT will not be equal to delay time from the corresponding EGT to MBS. Take January 2001 as an example: the moon will pass the descending node, perigee, ascending node, and apogee in sequence, and the sub-point of different MBS, as well as the location of the EGT with different look angles, may vary during periodic motion, as shown in Table 1.
In this study, three sample MBS's locations, (30 • W, 3 • S), (0 • , 0 • ), and (30 • E, 30 • N), appear to exhibit a centrosymmetric distribution because the Moon is locked by Earth tide, and the zero-longitude line of the Moon is almost facing the Earth. Here, the EGTs are assumed to be the intersection points of the MBS's antenna bore-sight on the Earth's surface. The antenna has been steered into the zero-Doppler plane, and the look angle is left-looking. As displayed in Table 1, the EGTs' latitude changes are greater than the corresponding longitude changes when the Moon reaches apogee and perigee, indicating that the intersection line of the zero-Doppler plane on Earth's surface is almost parallel to the longitude line when the Moon is at apogee and perigee. However, an angle exists because of the angle between the Moon's rotation axis and revolution orbit when the Moon is at either the ascending node or the descending node.
The influence of MBS's location on the sub-point, and the influence of the look angle on the intersection points suggest that the MBS cannot be deemed to be consistent with the Moon's centroid, although the Moon's radius is far smaller than the two planets' average distance, if a more realistic scenario had been presented.

Apogee
Perigee Ascending Node

Descending Node
The diversity of MBS's location and look angle have different numerical effects on the dual-delay of transmission and the maximum difference; therefore, because of the discrepancy caused by the delay time between uplink and downlink, using MBS and EGT, the results will no longer be constant with the changing of the synthetic aperture time and the look angle, as shown in Tables 2 and 3, suggesting that it is necessary to calculate the two delay time separately for the improvement of accuracy, the discrepancy when synthetic aperture time is 150 s is larger than when the time is 80 s if other parameters are the same, mainly because with an increase in the synthetic aperture time, the fringe angle, i.e., the angle from sampling points at fringe to the zero-Doppler plane, also increases, thereby generating unequal dual-delay time, which is more obvious once MBS is at the ascending node or the descending node. If the load shown in this paper were deployed in a space-born platform, the discrepancy of dual-delay time would be much smaller, mainly because of the shorter distance from the platform and the EGT, as well as the available stop-and-go mode. MBS has many advantages (long-term, large-scale, unique, etc.) over the space-born platform when conducting Earth observations. Compared with the geostationary (GEO) SAR, the average distance of MBS is about 9.5 times longer; the field of the whole disc of the Earth at such a long distance is just a small angle, around 1.8 • , which makes it possible for nearly 50% of the Earth to be observed at a gaze look, when echo reflection characteristics, which can be expressed by incident angle, are ignored. Once incident angle is considered, the field of view (FOV) of MBS is also larger than that of the space-borne platform. Moreover, the synthetic aperture length of MBS is also longer than GEO-SAR's, thus increasing the EGT's exposure time. The Moon is locked by Earth's tides, and almost only one side faces the earth, which leads to a basically fixed line of sight for the Earth's observation. Therefore, it is necessary and beneficial to develop the MBS for Earth observation, although range discrepancy does exist.
Finally, we will describe the necessity of using separate calculations for dual-delay time. Another processing is the fitting of actual trajectories deduced from ephemeris, such as those related to precession, nutation, apsidal motion, etc., in a specific scenario for better imaging. Through Equation (15), the fitting speed of MBS and maximum error in slant range (MESR), i.e., the absolute value of difference between fitting range and the range calculated in Section 3.1, around apogee, perigee, ascending, and descending nodes are displayed in Tables 4 and 5. Generally speaking, the fitting speed is close to the speed of the platform, but the actual fitting speed is significantly different from the speed of the MBS or EGT. Along with an increase in synthetic aperture time and look angle, the MESR will also increase slightly. The reason might be related to MBS's motion status, which can also be shown by the speed components of MBS's sub-point on Earth's East-North-Up (ENU) orientation. The ENU changes continuously when Moon is at four different orbit points shown in Table 6; notably, the North-Up speeds at the ascending node and descending node are faster than the speeds at the apogee and perigee, the MESRs in ascending node and descending node are also larger than those in apogee and perigee. In addition, the limit wavelengths in the ascending node and descending node are longer than those in the apogee and perigee, according to Equation (8). When the Moon is in apogee and perigee, the speed components are more closely related to the orientation of EGT's speeds, i.e., the east orientation, thus decreasing the MESRs under the same parameters.   Table 6. The speed of MBS in ECR, the distance of Moon-Earth centroid, and the speed of corresponding sub-point in ENU orientation. The line parallel to latitude and direct to east is called "E" orientation; the line parallel to longitude and direct to north is called "N" orientation; and the other line that satisfies the right-hand rule with the "E" and "N" is called "U" orientation. The time of apogee, perigee, ascending, and descending nodes are consistent with the time in Tables 2 and 3. For the process of imaging simulation in time domain, the smaller the MESR is, the better the imaging representing the actual scenario. Equation (8) indicates that the higher the carrier frequency f is, the smaller |∆R| max will be. However, Tables 4 and 5 indicates that |∆R| max cannot be too small, because MSER should not exceed |∆R| max in different synthetic aperture time periods and look angles. Therefore, to integrate the experience operation of the space-borne platform and ground-based imaging radar technology, the proposed FBP algorithm will be adopted on MBS in the Earth-Moon motion scenario, as shown in the next section.

Imaging Characteristic between MBS and EGT
Section 3.1 analyzes the range characteristic of MBS at four moments during a specific period. After the real trajectory has been fitted under suitable parameters, the FBP algorithm can be used for final image. Because the load on TerraSAR-X, ground-based imaging radar at Kashi, China, etc., are all X-band continuous waves, high-frequency bands have higher transmission power, and are capable of reaching kilowatts that can compensate for power loss during transmission; a band of high frequency will be less affected by the ionosphere than will the low-frequency P-band [21]. Along With increasing synthetic aperture time, the background ionosphere changes in temporal and spatial domains will also cause more serious azimuth shift and image defocusing for low-frequency bands. Here, the antenna aperture in the range or azimuth direction is up to tens of meters, for a better view in the long Earth-Moon distance; as shown in Table 7, the range width angle with only 0.028 • still shows about a 416 km swath in range when MBS is at apogee. Among periods of the Moon's revolution, slightly considering the influence of the incidence angle and the sub-point on echo, the look angle is selected as a middle number, i.e., 0.3 • . Since the EGT's image of MBS does not require high resolution when considering a large coverage, the range bandwidth can be set to a relatively low value. The basic system parameters of MBS are shown in Table 7, and the processed nine points image of MBS's echo at four moments corresponding to the apogee, perigee, ascending, and descending nodes are expressed in Figure 6. While the final compressed image shows little difference at these four moments, the target echo is still validly compressed using the proposed FBP algorithm. Among the simulation results, the Peak Side Lobe Ratio (PSLR), in row slice and column slice, both exceed −13 dB when MBS is at apogee and perigee, but are shown as somewhat smaller when the MBS is at the ascending node and the descending node, indicating that the images at apogee and perigee are slightly better than those for other two moments.
cending, and descending nodes are expressed in Figure 6. While the final compressed image shows little difference at these four moments, the target echo is still validly compressed using the proposed FBP algorithm. Among the simulation results, the Peak Side Lobe Ratio (PSLR), in row slice and column slice, both exceed −13 dB when MBS is at apogee and perigee, but are shown as somewhat smaller when the MBS is at the ascending node and the descending node, indicating that the images at apogee and perigee are slightly better than those for other two moments. Figure 6. The images generated by using the FBP algorithm when the MBS is at apogee (a), perigee (b), ascending node (c), and descending node (d), and their respective time is consistent with the time in Tables 2 and 3.
Moreover, the FBP's computational cost of projection calculation in the imaging area on the compressed data is N × √M × √N, the computational cost of phase interpolation is also √N × √M × √N, and the cost of coherent superposition is √M × √N. Therefore, the total computational cost is the sum of the above three computational costs, i.e.,(√N + N + 1) • √M × √N, and it is much smaller than the traditional computational cost of BP. Through the analysis of MBS, this FBP algorithm in the time domain based on the division of aperture and image area can be applied to avoid the coupling analysis in the frequency domain and to reduce the large computational cost of the BP algorithm. Figure 6. The images generated by using the FBP algorithm when the MBS is at apogee (a), perigee (b), ascending node (c), and descending node (d), and their respective time is consistent with the time in Tables 2 and 3. Moreover, the FBP's computational cost of projection calculation in the imaging area on the compressed data is N a × √ M × √ N, the computational cost of phase interpolation is also √ N a × √ M × √ N, and the cost of coherent superposition is √ M × √ N. Therefore, the total computational cost is the sum of the above three computational costs, i.e., √ N a + N a + 1 · √ M × √ N, and it is much smaller than the traditional computational cost of BP. Through the analysis of MBS, this FBP algorithm in the time domain based on the division of aperture and image area can be applied to avoid the coupling analysis in the frequency domain and to reduce the large computational cost of the BP algorithm.

Conclusions
The Earth's observation from MBS has a long transmission time; its dual-way path in atmosphere and ionosphere, also appears different when viewed from a space-borne platform; its complex environment undoubtedly increases the complexity in frequency domain. In addition, the stop-and-go mode is no longer suitable for MBS; since positions change during a single pulse, the motion status of MBS and corresponding EGT cannot completely be represented by orbit mode used in satellite, but a method including ephemeris that comprehensively considers the revolution or rotation of Earth-Moon motion, has been used instead. Here, the uplink and downlink components of the pulse transmission path are separately calculated for a better match with reality at four main moments, i.e., the apogee, perigee, ascending and descending nodes, during periodic motion with JPL DE430. Ephemeris has been adopted, and the calculated path method shown in this paper fully represents a real scenario for specific MBS distributed around the near side of Moon. Although the average radius of Moon is much smaller than the average distance between Earth and Moon, this paper no longer considered a point on Moon's surface to be same as the Moon centroid, as adopted by other studies. The results indicate that the longer the synthetic aperture time is, the greater the difference between uplink and downlink will be; therefore, synthetic aperture time is not the longer better for MBS. Meanwhile, the look angle, which determine the position of EGT, also affects the maximum difference, so the dual-path transmission should be calculated separately and then transformed to a high-order polynomial series as a fitting formula. In Earth-Moon motion, the fitting velocity is not close to the real velocity calculated in a specific coordinate system, which is not identical to an air-borne or space-borne platform. For MESR, its value will increase along with the increase in synthetic aperture time and look angle; the smaller its value, the better fitting it will represent.
After system parameters and MESR have been set or calculated, the division mode described in Section 2.1 can be used, and the FBP algorithm can be adopted. In contrast to the traditional BP algorithm, the renewed FBP algorithm will no longer project echo to the imaging area, but instead, a centerline will fully represent the imaging area; meanwhile, the synthetic aperture length and image area are divided into sub-parts for fast calculation. Although the FBP algorithm proposed might have a detrimental effect of reduced resolution, its meter-level resolution still meets the demand of MBS with a small width angle. A large amount of data will be generated along the long synthetic aperture time of MBS at the same sampling rate adopted in other platform; lower PRF is not suitable for the MBS to reduce the data received because of the possible azimuth ambiguity. The range interpolation before superposition indicated in this paper will further decrease the error of true position projected on the assumed sub-image area's diagonal, i.e., a segment of the centerline. The compressed results shown in Section 3.2 indicate that imaging processing is efficient in the FBP algorithm, and the results, when MBS is at apogee and perigee, are slightly better than the results at the ascending node and descending node.
In summary, the paper explains the necessity of calculating a dual path for MBS's single pulse transmission in Earth-Moon motion, with real ephemeris being adopted; meanwhile, a simplified formula for the dual path of a single pulse is derived. Then, a detailed range difference, mainly at apogee, perigee, ascending, and descending nodes is analyzed, also indicating the necessity of high-order polynomial fitting after the dual path of a single pulse been calculated. Finally, the simulation of MBS under acceptable resolution indicates the applicability of the algorithm and the reduction of computational cost. The content described in the paper might help people better understand MBS in an actual scenario, encouraging practical applications for future research.