3.1. Analysis of Long Time-Series Laser Pointing Change
To verify the accuracy and reliability of this method, the LFIs of the ascending orbit (at night) and complex ground objects were used as experimental data. As shown in
Figure 5a, for night conditions, while avoiding a complex background, the centroid extraction accuracy of this method was within 0.05 pixels and the centroid extraction accuracy of the GCM was approximately 0.08 pixels. As shown in
Figure 5b, under the influence of complex ground objects, the extraction accuracy of the centroid was within 0.08 pixels, which was approximately 2.5 pixels better than the result of the GCM. The error for the GCM mainly derived from exposure dispersion adjacent to the laser spot contour and the influence of background objects; this method removed these interference factors through multiple constraints. Under more complex background object conditions, the gap between the gray barycenter method and algorithm used in this study would be even greater.
To test the accuracy of the TEFM algorithm, 1600 LFIs with similar times and enriched ground object types were randomly selected for accuracy verification. Taking the calibration positions of the two beams in orbit as real values, the mean coordinate error, range, and root-mean-square error (RMSE) were calculated (
Table 3). The statistical results show that the traditional algorithm could not extract the centroid of the laser spots under a complex background. The TEFM algorithm maintained a high level of stability and overall extraction accuracy, which was one order of magnitude better than the traditional algorithm. Thus, the proposed algorithm conformed to the requirements of subpixel centroid extraction under complex background conditions.
From 15 to 30 March 2020, a total of 61 tracks of LFI data were extracted and 86,000 images were obtained. Among them, 27% of the images lost the spot contour information owing to cloud cover, overexposure, or other reasons, which made it impossible to extract the centroid.
Figure 6 shows the experimental results. The abscissa represents the X-direction and the ordinate represents the Y-direction. The left side shows the statistical results for the centroid coordinates of beam #1, with mean coordinates of (120.83, 263.71). The X-direction oscillates at approximately 0.4 pixels, while the Y-direction oscillates at approximately 0.5 pixels. The plane position changes within 1 pixel. This is the statistical result for the centroid coordinates of beam #2 on the right side, with mean coordinates of (216.32, 162.51). The X-direction oscillates at approximately 0.5 pixels while the Y-direction oscillates at approximately 0.5 pixels. The plane position changes within 1 pixel. Therefore, we obtained the following preliminary conclusions: (1) compared with the traditional algorithm, the improved algorithm effectively extracts the centroid of the laser spot under a complex ground object background; and (2) the change in the centroid coordinates is relatively stable every month, in which the amplitude oscillation in the X- and Y-directions is <1 pixel.
The difficulty in establishing a long-term laser pointing monitoring system or extracting the centroid of a long time-series lies in eliminating gross errors and retaining the changing trend(s) of the centroid coordinates. Although many conditions can be used to constrain the spot contour and accurately locate the spot centroid position, the influence of background objects on the spot cannot be completely removed, which leads to gross errors, moves the coordinates of the spot centroid, and destroys the existing coordinate trend. There are two types of gross errors: (1) a high-reflectivity ground object completely covers the spot contour or (2) the amplitude of the laser spot itself is low and dim, such that the features of the light spot are hidden by the features of the ground objects. After a large number of statistical calculations, we found that the spot characteristic parameters, such as the eccentricity and half-axis length, could effectively identify gross errors and improve the overall recognition accuracy. To analyze the change and stability of the centroid of the footprint facula, analyses were carried out from both macroscopic and microscopic perspectives.
As listed in
Table 4, to perform long-term analyses on the stability of the spot centroid in the LFI since its orbital launch, the monthly mean coordinates of the facula centroid from March 2020 to April 2021 were counted. As shown in
Figure 7a,b, for beam #1, the centroid coordinates showed a decreasing trend and changed by approximately 0.4 pixels in the X-direction and 1 pixel in the Y-direction. The plane position changed by approximately 1.1 pixels; the corresponding directional angle changed by approximately 0.341”. As shown in
Figure 7c,d, for beam #2, the centroid coordinates first decreased and then increased by approximately 0.4 pixels in the X-direction. There was an overall decreasing trend in the Y-direction by approximately 1.5 pixels, a change in the plane position by approximately 1.4 pixels, and a change in the corresponding pointing angle by approximately 0.434”. For beams #1 and #2, the monthly change in the centroid of the spot did not always move in a specific direction; the monthly change had a relatively small range.
The laser footprint camera and the laser share the same receiving field of view. When the laser pointing jitter is small, the relative position of the laser footprint image and the laser spot is basically unchanged, and transmitting and receiving are coaxial, which can be regarded as the same reference coordinate system. In order to analyze the change in the plane accuracy of the laser data of the GF-7 satellite in orbit, this section randomly checks the plane accuracy of the long-time laser footprint image data. Google images, airborne aerial photographs (from Germany, China and other regions), and Lidar point clouds were used as reference data for evaluating the positioning accuracy of laser footprints [
31,
32]. As shown in
Table 5, the positioning accuracy of lasers 1 and 2 is within two pixels, corresponding to a ground distance of ~6 m, which shows that the positioning accuracy of the laser footprint of the GF-7 satellite has remained relatively stable since its launch.
In the process of laser data processing and analysis of the GF-7 satellite, we found that the ranging and plane accuracy of the obtained laser spot is degraded to some extent when the satellite is rolling at a large angle; this may be related to: (1) the accuracy of the star sensor being reduced; (2) the angle of view of laser receiving becoming smaller; (3) pointing angle jitter at the time of emergent light; or (4) a change in the internal optical axis frame of the laser system due to temperature change, material thermal deformation, and other factors. We refer to the data obtained when the satellite rolling angle is greater than or equal to 3 as rolling data, as shown in
Table 6, which is the random sampling result for evaluating the plane accuracy of rolling data. When the satellite swings sideways, the plane accuracy of beams 1 and 2 is ~two pixels or more, and the corresponding ground distance is 6–20 m. Compared with the case of no rolling, the plane accuracy of the footprint is greatly degraded, and the elevation value is uncertain. At present, the problem of laser ranging accuracy degradation caused by satellite pendulum measurement has not been solved, and the satellite rolling angle (3) can only be used as a necessary condition for quality control to control the overall accuracy. In this process, we found that satellite rolling has a certain influence on the distribution of laser spot energy (see
Section 3.2).
3.2. Analysis of Long Time-Series Laser Energy Changes
To eliminate the influence of complex background objects on the radiance of the laser spot, we used nighttime data from April to July 2021 for the analysis. There were four tracks (one track sampled every month) and 6216 laser spots.
Figure 8 shows the maximum amplitude of the center point, energy inclusion diagram, and OTF-LESE from top to bottom; beams #1 and #2 are shown from left to right. This analysis provided the following findings:
- (1)
Brightness of center disk: As shown in
Figure 8a,b, for beam #1, the maximum amplitude fluctuated between 5100 and 5600, with an amplitude jitter of 400 (dimensionless amplitude value) during the track crossing stage. The center energy gradually diffused after long-term operation. For beam #2, the maximum amplitude fluctuated between 2000 and 2600, with a jitter of 500 during the track crossing stage. The center energy was stable and tended to increase gradually during long-term operation.
- (2)
Encircled energy diagram: As shown in
Figure 8c,d, the maximum slope represented the spot boundary; the spot radius of beams #1 and #2 was between 8 and 10 pixels. With the centroid coordinate as the center and a radius within 20 pixels, the total energy of the scattered spot of beam #1 was approximately 1,200,000 (dimensionless amplitude value), while that of beam #2 was approximately 600,000 (approximately half that of beam #1). The amplitude value was not equal to the laser emission energy; there was a certain mapping relationship.
- (3)
OTF-LESE: the maximum amplitude of the center point can only evaluate the change in the center value of the laser spot and the amount of energy lost, whereas the energy inclusion diagram can only show the energy dispersion degree adjacent to the spot but cannot fully evaluate the energy change at the time of light emission. As shown in
Figure 8e,f, the range of the OTF-LESE was 0–1 under normal conditions; if it exceeded 1, the transmission waveform gain was too small. Compared with the first two indices, the OTF-LESE indicated that the periodic changes caused by the pointing jitter were considered during energy changes at the time of laser exit; there were notable periodic changes at the time of the crossing orbit for beams #1 and #2. For beam #1, the OTF-LESE changed within 0.7–1.3, with an average value of 0.91. The maximum amplitude changed to 0.3 when crossing the orbit, and the energy decayed to 0.7 with continuous operation. For beam #2, the OTF-LESE varied within 0.7–1.05, with an average value of 0.85. The maximum amplitude changed to 0.3 when crossing the orbit, and the energy decayed to 0.72 with continuous operation. This value was still within the normal working range for beams #1 and #2.
Figure 8.
Analysis of changes in the laser energy at the emission time. (a,b) represent respectively center disk brightness of beam 1&2. (c,d) represent respectively Encircled energy diagram of beam 1&2. (e,f) represent respectively OTF-LESE of beam 1&2.
Figure 8.
Analysis of changes in the laser energy at the emission time. (a,b) represent respectively center disk brightness of beam 1&2. (c,d) represent respectively Encircled energy diagram of beam 1&2. (e,f) represent respectively OTF-LESE of beam 1&2.
OTF-LESE is an index used to evaluate the output state of a satellite laser; using other values as the true values to verify its accuracy remains difficult. Here, we focused on evaluating the sensitivity and comprehensiveness to various changes in the laser state. As shown in
Figure 9 and
Table 7, from left to right, the OTF-LESE was 0.7, 0.8, 0.9, and 1.0. The upper sequence diagram shows the transmission waveform while the lower sequence diagram shows the LCPA data. When the satellite laser was launched, the peak and amplitude of the emission waveform, coordinates for the centroid of the spot in the LCPA, and energy distribution, were in a relatively stable state within a small range. Taking the laser state of the rightmost OTF-LESE as a reference, the plane coordinate for the center of mass of the leftmost laser LCPA changed by approximately 0.8 pixels, the spot radius decreased by approximately 1.2 pixels (with 80% of the central energy amplitude as the constraint based on the energy inclusion diagram), and the central energy amplitude decreased by approximately 80. The peak value of the leftmost emission waveform decreased by approximately 10, the kurtosis coefficient changed by approximately 0.66, the skewness coefficient changed by approximately 0.07, and the waveform width decreased by approximately 0.23. Based on these parameters, the OTF-LESE can represent the typical characteristics of the output state of the satellite laser.
On 3 December 2020, the GF-7 satellite passed over Lake Tanganyika at night, acquiring 18 laser spots, including 10 points from beam #1 and 8 points from beam #2 (
Figure 10). Combined with meteorological data, it is confirmed that there were no waves on the lake when the data were obtained, which would not cause additional observation error. Tanganyika is a freshwater lake in central Africa, with a coastline of 1900 km. Considering the influence from the curvature of Earth, the elevation of the lake presents a linear trend along the track. Therefore, the height profile data for beams #1 and #2 along the rail direction were fitted separately. The height error of the corresponding position was calculated as the true value. The maximum error of beam #1 was 7 cm and that of the OTF-LESE was 0.68. The maximum error of beam #2 was 14 cm and that of the OTF-LESE was 0.63. Regardless of the stability of the working mode or altimetry accuracy, beam #1 was superior to #2.
To further verify whether the OTF-LESE could effectively estimate the emission state of the satellite laser and further estimate the altimetry error caused by the change state, we analyzed the correlation between the altimetry error and OTF-LESE of the data on the lake, collected using the two beams. As shown in
Figure 11, the OTF-LESE had a positive correlation with the altimetry error caused by the jitter of the laser state. The Pearson correlation coefficient between the two variables was 0.78, which shows that the OTF-LESE effectively evaluated the emission state of the satellite laser and provided data quality control for the final altimetry product.