Next Article in Journal
Satellite Remote Sensing Identification of Discolored Standing Trees for Pine Wilt Disease Based on Semi-Supervised Deep Learning
Previous Article in Journal
Evaluation of LiDAR-Derived Features Relevance and Training Data Minimization for 3D Point Cloud Classification
Previous Article in Special Issue
On-Orbit Radiometric Calibration of Hyperspectral Sensors on Board Micro-Nano Satellite Constellation Based on RadCalNet Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Separate Calibration Method of Laser Pointing and Ranging for the GF-7 Satellite Laser That Does Not Require Field Detectors

1
The School of Earth Sciences and Engineering, Hohai University, Nanjing 211100, China
2
Land Satellite Remote Sensing Application Center, Ministry of Natural Resources of China, Beijing 100048, China
3
College of Geodesy and Geomatics, Shandong University of Science and Technology, Qingdao 266590, China
4
School of Smart City, Chongqing Jiaotong University, Chongqing 400074, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(23), 5935; https://doi.org/10.3390/rs14235935
Submission received: 29 September 2022 / Revised: 16 November 2022 / Accepted: 21 November 2022 / Published: 23 November 2022

Abstract

:
Satellite laser altimeters have been widely used in the surveying, mapping, forestry, and polar regions and by other industries due to their excellent elevation measurement accuracy. Satellite laser on-orbit geometry calibration is a necessary means to ensure elevation accuracy. This study proposes an iterative geometry calibration method for satellite laser altimeter pointing and ranging separation that does not require the use of field detectors. The DSM data were first used to complete the laser pointing calibration, and then the laser footprint elevation was measured accurately to complete the laser ranging calibration. The iterative calibration experiment was repeated until the convergence condition (i.e., the laser point difference was less than 1 × 10-5 degrees and the laser ranging difference was less than 0.01 m) was met, with the calibrated laser pointing angle and ranging separation used as the input parameters. In this work, the GaoFen-7 (GF-7) satellite laser was used as the test object and the actual laser pointing and ranging values derived from ground detector calibrations. The results verified that the pointing accuracy of the GF-7 beam 1 was 2 arcsec and that the ranging accuracy was 2 cm after applying the calibration method presented in this paper. The pointing accuracy of the GF-7 beam 2 was 2.2 arcsec, and the ranging accuracy was approximately 1 cm. This analysis demonstrated that the GF-7 laser mission exceeded its pointing angle requirement of 3 arcsec after laser pointing and ranging separation iterative calibrations were applied. Finally, ground control points were used to verify the calibrated elevation accuracy of the GF-7 satellite laser, and its accuracy on flat terrain was 0.18 m. In summary, it was proven that the satellite laser geometry calibration method proposed in the article is effective.

1. Introduction

The GaoFen-7 (GF-7) satellite was launched on 3 November 2019 and is equipped with two full-waveform profile laser beams for synchronous Earth observation [1]. The operating modes of the two GF-7 laser beams are identical. For example, the wavelength is 1064 nm, the measurement frequency is 3 Hz, and the transmitted and received waveforms are recorded at a sampling rate of 2 GHz. The two laser beams are positioned on either side of the across track at the nadir of the satellite, and each beam is at an angle of 0.7 degrees from the nadir; that is, the distance between the two laser footprints is approximately 12.25 km across track. For more details, please refer to reference [2]. The GF-7 satellite has been in orbit for 3 years and has acquired a large amount of data from the surface of the globe. These data products have been released on the Natural Resources Satellite Remote Sensing Cloud Service Platform (http://sasclouds.com/chinese/normal/, accessed on 20 November 2022) in the *.SLA format. The detailed application process for the GF-7 satellite laser data is in Appendix A. These data have been widely studied by domestic and foreign researchers and engineers. Therefore, ensuring the accuracy of GF-7 satellite laser data is one of the most vital tasks during satellite operation. Thus, the periodic on-orbit geometry calibration of the lasers is an important means to ensure data accuracy.
At present, there are two primary high-precision calibration methods for existing profile satellite lasers: (1) The satellite attitude maneuver calibration method [3,4,5,6], which causes the satellite laser to rotate within a cone over the middle of the ocean. Then, according to the principle of minimum residual error of satellite laser ranging, the geometric calibration of the satellite laser is completed. (2) The calibration method based on laser ground detectors (LGDs) [7,8,9,10]. This method actively captures laser footprints by placing multiple infrared detectors on the surface, and satellite laser calibration is accomplished by triggering the detector. Other satellite laser calibration methods have not been used for high-precision satellite laser calibration due to low accuracy, a low success rate, or difficulty in implementation and so on. For example, the calibration method based on the synchronous imaging measurements through an infrared camera [11] requires the aircraft to fly synchronously with the satellite at night. This method is difficult to implement and has a low success rate. The simulation waveform matching calibration method [12] has some uncertainties with high requirements for the local terrain type, and the reliability is low. The traditional terrain matching calibration method [13] cannot calibrate laser range errors, and the calibration accuracy is low. Additionally, there is an iterative pointing angle calibration method based on small-range terrain matching [14]. This method establishes a calibration error model, and the calibration is carried out by solving the partial derivatives of the pointing angle and ranging. However, this method requires a large number of laser spots and thus cannot be used for profile satellite lasers with low pulse frequencies. However, high-precision calibrations can be accomplished by maneuvering the satellite at attitude, though this method requires the high attitude maneuverability of the satellite platform. For satellites with platform-stabilized lasers, such as the GF-7 satellite, rapid cone measurements cannot be performed. Therefore, only the calibration method based on LGDs is widely used on platforms such as the GF-7 and ZY3-02/03 satellites. The disadvantage of this method is that it requires that staff place a large number of LGDs in the calibration field, and the workload is substantial.
Multibeam satellite lasers are the current trend in satellite laser development, such as the GF-7 dual-beam lasers and the five-beam lasers of terrestrial ecosystem carbon monitoring satellites. The cost of the LGD calibration method is large, and the cost doubles with the increase in beams. Therefore, it is critical to develop a low-cost satellite laser geometry calibration method that does not require ground detectors. Due to this urgent demand, this work proposes an iterative geometry calibration method for satellite laser altimeter pointing and ranging separation that does not require field detectors. First, the residual sum of the elevation of all the GF-7 laser footprints and the reference digital surface model (DSM) were obtained by terrain matching. The satellite laser was then calibrated according to the sum of the minimum elevation residuals. Then, a flat and hard surface area was used as the test area for laser range calibration, and the elevation difference between the laser-derived elevations and the actual elevation was defined as the laser ranging error. The actual elevation of the test location was obtained by accurately measuring the flat surface ground with a precise level. Next, taking the calibrated pointing angle and range of the satellite laser as the initial value, the satellite laser pointing and ranging calibration experiment was repeated until the laser point difference was less than 1 × 10−5 degrees and the laser ranging difference was less than 0.01 m, and the calibrated satellite laser pointing and ranging error values were the output. In this work, the GF-7 satellite laser was used as the test object, and the calibrated pointing and ranging error of the GF-7 satellite laser based on LGDs were used as the true values, which were used to verify the laser pointing and ranging accuracy of the GF-7 satellite after the calibrations based on the method presented in this paper were applied. Additionally, the ground control points (GCPs) were used to verify the laser elevation accuracy of the GF-7 satellite after the calibration.

2. Materials and Methods

2.1. Study Area

2.1.1. Calibration Area

(1)
Point Calibration Area
Considering the minimal time phase difference between the laser data and the reference DSM data, there was a little elevation error between the systems caused by the natural growth of surface vegetation. Even so, this study preferentially selected the bare natural surface as the calibration area. A region in the east Gobi desert near the city of Kuqa, Xinjiang, China, was selected as the GF-7 satellite laser pointing calibration area. The area was 83.30–83.86°E, 41.19–42.03°N and consisted of a total area of approximately 2000 square kilometers. The topography was very complex in this area, the landforms undulated considerably, and there was little tall vegetation on the surface. The elevation changes within the calibration area have been negligible over recent decades; thus, it is an ideal experimental area for satellite laser pointing calibration. This experimental area is shown in Figure 1a.
(2)
Ranging Calibration Area
According to the calibration strategy proposed in this paper, a hard and flat surface was selected for the laser ranging calibration area. An airport runway near Beijing city, as shown in Figure 1b, was selected as the GF-7 laser beam 1 ranging calibration area. The ranging calibration laser footprint time code of the GF-7 beam 1 was 203685561.00 (the time code is the accumulated seconds since 1 January 2014), and there was no elevation difference over approximately 200 m, and thus this location can be regarded as falling within the same elevation plane. The ranging calibration area of GF-7 laser beam 2 was selected along a flat highway in the suburbs of Beijing city, as shown in Figure 1c. The ranging calibration laser footprint time code of GF-7 beam 2 was 203685553.00, and the surface elevation was the same within approximately 100 m of this feature. The satellite laser pointing accuracy is normally approximately 30 m [13] after the first terrain matching calibration. Within this error range, the actual elevation of the ground for the selected laser footprints would not change, and the selected ranging calibration area was reasonable. Finally, the DS05 level was used to accurately measure the surface elevation of the ranging calibration area, which was used as the true elevation for the calibration of the laser ranging error.

2.1.2. Elevation Accuracy Validation Area

The accuracy of satellite-derived laser elevation measurements is a topic that has garnered extensive attention. Surface GCPs are measured by real-time kinematic GPS (RTK) to accurately verify satellite-derived elevation data. The areas used for verification should be flat terrain. Sonid Youqi, Inner Mongolia, China, was selected as the area to verify the GF-7 laser calibrations discussed in this work, as shown in Figure 2. The terrain was relatively flat in this location and tall vegetation was minimal, consisting primarily of xerophytes and small grasses.

2.2. Data

2.2.1. GF-7 Satellite Laser Data

In this article, the GF-7 satellite laser data are divided into three categories: laser pointing calibration data, laser ranging calibration data, and elevation accuracy verification data, as shown in Table 1. The GF-7 pointing calibration data came from the 5278th track on 14 October 2020, the 6160th track on 12 December 2020, and the 8851th track on 7 June 2020, of which beam 1 had a total of 83 laser points and beam 2 had a total of 68 laser points. The laser spatial distribution is shown in Figure 1a. The GF-7 ranging calibration data came from the 3419th track on 15 June 2020. The laser point timecode of beam 1 was 203685561.00 and 203685553.00 for beam 2, and their distributions are shown in Figure 1b,c, respectively. The GF-7 elevation verification data were from the 3326th track on 9 June 2020, the 3402th track on 14 June, the 3478th track on 19 June, and the 3555th track on June 24. Among these, beam 1 and beam 2 both contained 24 laser points used for elevation accuracy verification, as shown in Figure 2.

2.2.2. Reference Terrain Data

(1)
AW3D30 DSM
ALOS World 3D-30 m (AW3D30) is a global DSM dataset with a 30 m grid resolution released by the Japan Aerospace Exploration Agency (JAXA) Earth Observation Research Center in 2016 [15]. Users can download the global dataset after registration. The data download link of the AW3D30 dataset is https://www.eorc.jaxa.jp/ALOS/en/aw3d30/data/index.htm (accessed on 20 November 2022). The global elevation accuracy of the AW3D30 is better than 5 m [16]. We adopted the national control point image database to verify the elevation accuracy of the AW3D30 DSM, and the elevation accuracy of the AW3D30 DSM in the laser pointing calibration area was approximately 2 m [17]. The AW3D30 DSM of the pointing calibration area is shown in Figure 3a.
(2)
LiDAR Data
We used a fixed-wing piloted aircraft to carry the Leica TerrainMapper Light Detection and Ranging (LiDAR) to obtain surface terrain data of the pointing calibration area in September 2021. During the surface terrain data measurement, the aircraft maintained a flight altitude of 5200 m, and the laser emission frequency was 436 kHz. The average point density on the surface was guaranteed to be 1.6 points per square meter when the aircraft was at the maximum flight speed of 309 km/h. After completing field measurements, the filtered LiDAR point cloud data were used to generate a DSM image with a resolution of 1 m and an elevation accuracy of 0.17 m, which is shown in Figure 3. In this article we refer to this DSM image as LiDAR data.

2.2.3. GCPs Dataset

The GCP data in the paper refer to the RTK-derived surface elevation control points set around the GF-7 laser footprint. An example of a GCP measurement is shown in Figure 4b. First, the 81 GCPs were measured around the laser footprint after calibration, such that an array of GCPs was formed around the laser footprint. The GCP array is shown in Figure 4a, where the black dots represent the GCPs and the green dot is the laser point 203168336.67 from GF-7 laser beam 1. The distance between the GCPs in the GCP array was approximately 3 m, and the elevation accuracy of each GCP was better than 5 cm. Subsequently, the GCP arrays surrounding all the GF-7 laser points shown in Figure 2 were combined to form a GCP dataset.

2.3. Methodology

This study proposes an iterative geometry calibration method for satellite laser altimeter pointing and ranging separation that does not require the use of field detectors. The workflow of our proposed method can be summarized as follows (Figure 5):
Step 1:
Experimental data preparation: initial pointing parameters, initial ranging error (0 by default), GF-7 laser data, reference DSM data, and the GCPs around the laser footprint.
Step 2:
Laser pointing calibration: First, a pointing angle grid is established with the initial laser pointing angle as the center. Then, the residual sum of the laser and the reference DSM elevation for each grid is calculated. The laser pointing angle with the smallest elevation residual sum between the laser and the DSM is used as the laser pointing angle after calibration.
Step 3:
Laser ranging calibration: taking the laser pointing angle after calibration as the input parameter, the coordinates of the laser footprint on the flat hard surface are calculated, and the elevation difference between the laser and the surface is used as the ranging error after calibration.
Step 4:
Iterative of laser pointing and ranging calibration: It must be assessed whether the laser point difference (LPD) Δ p between the calibrated pointing and the input pointing in step 2 is less than 1 × 10−5 degrees and whether the laser ranging difference (LRD) Δ ρ between the calibrated ranging error and the input ranging error in step 3 is less than 0.01 m. This is called the iterative convergence condition. When the convergence condition is satisfied, the output pointing and ranging error are the final iterative calibration results. Otherwise, redo steps 2–4 with the calibration result of this iteration as an input parameter. Typically, it takes 2–3 iterations to converge.
Step 5:
Accuracy verification after calibration: the laser pointing and ranging error value after calibration by LGDs is first used as the true value to verify the laser pointing and ranging accuracy following the calibration method outlined in this paper. Then, the GCPs are used to verify the elevation accuracy of GF-7 satellite laser after calibration.

2.3.1. Satellite Pointing Error Calibration

The primary mathematical model for obtaining the satellite calibration data of the laser pointing, ranging error, and laser elevation accuracy verification presented in this paper is the satellite laser geometric positioning model. The geometric relationship between the satellite laser and the surface laser footprint was established, as shown in Figure 6.
The geometric relationship between the laser emission point, the satellite centroid, the GNSS antenna, the center of the Earth, and the laser footprint on the surface is shown in Figure 6. The laser footprint position vector r S is the sum of the GNSS antenna center position vector r G ; the vector r G B from the GNSS antenna center to the satellite centroid; the vector r G L from the satellite centroid to the laser emission point; and the laser ranging vector ρ p , where ρ is the laser ranging value and p is the laser pointing unit vector. Converting it into matrix form, the geometric positioning model of the satellite laser is established, as shown in the following formula [18]. The formula can be used to calculate the three-dimensional coordinates of laser footprints based on satellite orbit, attitude data, and laser ranging value.
X Y Z I T R F = X G Y G Z G I T R F + X G B Y G B Z G B I T R F + R I C R F I T R F R S B F I C R F X B L Y B L Z B L + ρ cos ( α x ) cos ( α y ) cos ( α z )
where X Y Z I T R F T are the surface coordinates of the laser footprint. X G Y G Z G I T R F T are the satellite centroid coordinates in the International Terrestrial Reference Frame (ITRF). R I C R F I T R F is the rotation matrix from the International Celestial Reference Frame (ICRF) to ITRF. R S B F I C R F is the rotation matrix from the satellite body frame (SBF) to ICRF. X G B Y G B Z G B I T R F T are the offset values between the GNSS antenna center and the satellite centroid. X B L Y B L Z B L T are the offset values between the satellite centroid and the laser emission point. ρ is the laser range value, which includes the ranging errors caused by tides and the atmosphere, as well as the ranging errors of the laser itself, where the latter ranging errors need to be corrected by geometric calibration. cos ( α x ) cos ( α y ) cos ( α z ) T = p is the laser pointing unit vector; α x is the intersection angle between the laser optical axis and X axis in the SBF; α y is the intersection angle between the laser optical axis and Y axis in the SBF; and α z is the intersection angle between the laser optical axis and Z axis in the SBF. Among them, α z can be directly calculated from the pointing angles α x and α y , and the calculation formula is shown in Formula (2). Due to the satellite vibration during launch and environmental changes after the satellite reaches orbit, there are pointing differences between the actual alignment α x and α y and the laboratory measurement results that also need to be calibrated.
α z = a cos 1 cos ( α x ) cos ( α x ) cos ( α y ) cos ( α y )
The satellite laser pointing error calibration method adopts the terrain matching method. First, the initial laser pointing angle ( α x , α y ) is used as the center to establish a rectangular pointing grid with the pointing angles α x and α y . It is assumed that the pointing grid range is L x and L y , and grid spacing is Δ L x and Δ L y , respectively. Then, each pointing angle in pointing grid is taken out as the input pointing angle, and the surface coordinates of all the laser footprints are calculated according to the satellite laser geometric positioning model. According to the longitude and latitude of each laser footprint, the elevation of the reference DSM is interpolated to calculate the elevation residual value of all laser footprint elevations and the DSM elevation. All laser footprint elevation residual sums are counted (abbreviated as laser elevation residual sums here). Finally, the previous step is repeated until all value in the pointing grids have been processed. The laser pointing angle with the smallest laser elevation residual sum is the optimal pointing angle, that is, the laser pointing angle after calibration. The detailed implementation process of the Algorithm 1. When the range of the pointing grid is large or the grid spacing is small, the experimental processing speed is slow, and the reference [13] can be used to improve the pointing calibration processing speed by using a pyramid method.
Algorithm 1: Satellite Laser Pointing Calibration
1: Input: The Initial laser pointing α x , α y . The pointing grid range L x , L y . The pointing grid interval Δ L x , Δ L y . The original laser data set P T = P t 1 , P t 2 , P t n .
 2: Output: The pointing after calibration α x o p t , α y o p t .
 3: Generate a laser pointing grid with a length of 2 L x and a width of 2 L y . So, the laser pointing set is L P S = α 1 , α 1 , , α 2 L x / Δ L x , α 1 , , α 2 L x / Δ L x , α 2 L y / Δ L y . There is total m = 2 L x Δ L x × 2 L y Δ L y laser pointing
 4: For  i = 1 to m  do
 5:  For  j = 1 to n  do
 6:   Compute the surface coordinates of the j t h laser footprint X j Y j Z j I T R F T using laser pointing α x i ,     α y i .
 7:   Interpolate the elevation h j of the reference DSM using X j Y j I T R F .
 8:   Calculate the elevation difference Δ h j between h j and Z j , Δ h j = a b s Z j h j .
 9:  end for
 10:  Count the mean of all laser elevation differences Δ H i = Δ h 1 + Δ h 2 + + Δ h n / n .
 11: end for
 12: Find the smallest in the elevation difference from dataset Δ H 1 , , Δ H m , and assume that its corresponding laser pointing is α x m i n , α y m i n . Thus, α x o p t , α y o p t = α x m i n , α y m i n .

2.3.2. Satellite Ranging Error Calibration

Satellite lasers are usually measured in the vertical nadir direction [2,19,20]; that is, the intersection angle between the laser optical axis and the Z axis in the SBF can be ignored. However, for the multibeam lasers of the GF-7 satellite, the intersection angle between the laser optical axis and the Z axis in the SBF is relatively large. That is, α z is approximately 0.7 degrees, so cos α z 1.0 . Therefore, according to Formula (1), it is easy to deduce that the laser footprint elevation error will be completely caused by the ranging error. Based on this theory, we propose a method to directly correct the laser ranging error by using the laser footprint elevation error. The key issue of this method is to accurately measure the surface elevation of the laser footprints. We will address these issues in two aspects. The first approach pertains to laser footprint selection, and the other involves accurately measuring the surface elevation of the laser footprints.
First, laser pointing after calibration is used to calculate the surface coordinates of the laser footprints that may fall on flat and hard ground, such as cement courts and ice. Additional laser footprints within 50 m of the initial laser footprint that fall within the same elevation plane were selected. Then, the DS05 precision level is used to precisely measure the leveling points to establish a fourth order leveling network. Then, the laser footprints are placed on the level network path to measure their elevation values. The elevation accuracies of the laser footprints is guaranteed to be better than 3 cm. Assuming that the elevation of the laser footprint after the pointing calibration is h l a s e r , and assuming that the actual surface elevation of the laser footprint measured by the DS05 level is h s u r f , then the satellite laser ranging error Δ ρ can be directly calculated according to the following formula:
Δ ρ = h l a s e r h s u r f

2.3.3. Elevation Accuracy Verification of Satellite Laser

After iterative calibration, the satellite laser pointing and ranging values are brought into Formula (1) to calculate the initial coordinate X Y Z I T R F T of the laser footprint, which is converted to geodetic coordinates using the WGS84 ellipsoid according to Formula (4). According to the coordinates of the laser footprint and laser emission time, the National Centers for Environmental Prediction (NCEP) data are used to correct the laser ranging delay error caused by atmospheric interference [21,22], and the FES2014 ocean tide model is used to correct the laser elevation tide error [23]. The coordinates of the laser footprint are recalculated using the corrected laser range as the final coordinates of the laser footprint.
L = arctan Y X B = arctan Z ( N + H ) ( X 2 + Y 2 ) [ N ( 1 e 2 ) + H ] H = X 2 + Y 2 cos B N
where N is the radius of the curvature of the prime vertical for the WGS 84 ellipsoid; e is the first eccentricity of the WGS84 ellipsoid; and H is the ellipsoid height of the laser footprint.
Taking the continuously operating reference stations (CORS) as the RTK base station, RTK GPS is used to measure the geodetic elevation value of the laser footprint, H . H is considered to be the true value of the laser footprint elevation in order to evaluate the elevation accuracy of the GF-7 satellite laser H after the calibration. Assuming that there are laser footprints, the elevation residual values d H of n laser footprints are calculated by Formula (5). Finally, the mean value Δ H ¯ and root mean square error (RMSE) σ H of the laser footprint elevation residuals are counted as the evaluation index of the laser elevation accuracy, and the calculation formula of the mean value and RMSE are as follows:
d H = H H
Δ H ¯ = i = 1 i = n d H i n
σ H = i = 1 i = n d H i 2 n

3. Results

3.1. Calibration of Laser Pointing and Ranging

In this work, all laser footprints of GF-7 beams 1 and 2 from the pointing calibration area in Figure 1a were used as the experimental objects, and AW3D30 DSM was used as the reference terrain data to conduct terrain matching experiments according to the method shown in Figure 7. For each beam of the GF-7 laser, the initial laser center pointing angle was measured in the laboratory before the satellite launch. A pointing grid range with a laser pointing angle α x and α y range of ±0.5 degrees and grid spacing of 0.003 degrees was created. This pointing grid had a total of 110,889 pointing angles. Assuming that the satellite orbits at an altitude of 500 km, the laser footprint position error was approximately 2.5 m per arcsec. The resolution of the AW3D30 DSM was 30 m, which is equal to approximately 0.003 degrees of the laser pointing angle. For each pointing angle in the pointing grid, the absolute value of the elevation residuals between the 83 laser footprint elevation of beam 1 and the AW3D30 DSM elevation were calculated. Then, the sum of the absolute values of all the laser footprint elevation residuals of beam 1 was calculated, and the result is referred to as the sum of elevation residuals. Finally, after traversing all the pointing angles in the pointing grid, a set of elevation residuals was generated for the grid pointing angles and the sum of the elevation residuals. The calibration experiment was consistent for beams 2 and 1, and the three-dimensional surface plots of the grid pointing and the sum of the elevation residuals of beam 1 and beam 2 are shown in Figure 7a,b. In the figure, Δ H refers to the elevation residuals between the laser and reference terrain at different pointing angles, and α x and α y are two pointing angles of laser boresight in Formula (1).
From the results of Figure 7a,b, there is an obvious minimum value on the surface between the laser pointing angle and the sum of the elevation residuals. It is easy to understand that as the laser pointing accuracy increases, the sum of the laser elevation residuals decreases. Therefore, the laser pointing angle that corresponds to the minimum sum of the elevation residuals in Figure 7a,b is the optimal pointing angle of GF-7 laser beams 1 and 2 after the first iterative calibration.
After the first calibration was completed, the laser pointing angle was used to calculate the elevations of the laser footprint in Figure 1b,c for the ranging calibration. The elevation difference between the laser elevation and surface elevation was measured by RTK, and this elevation difference value was taken as the laser ranging error after the first calibration and is shown in Table 2. Then, the GF-7 laser beam 1 and beam 2 pointing angle after the first iterative calibration were taken as the center pointing angle, and a new pointing grid with the pointing range of ± 0.03 degrees and the grid spacing of 0.003 degrees was created. The first iterative calibration was repeated using the new pointing grid, and the GF-7 laser pointing after the second iterative calibration was calculated. The sum of the elevation residual surfaces after the second iterative calibration is shown in Figure 7c,d. After the second pointing iterative calibration, the new laser pointing angle was used to recalculate the elevations of the two laser footprints in Figure 1b,c. The elevation difference between the laser elevation and the true surface elevation was used as the laser ranging error to complete the second iterative pointing and ranging calibration. Finally, the third iterative calibration takes the laser pointing angle and ranging error after the second iterative calibration as the input parameters and completely repeats the second iterative calibration experiment. The sum of the residual surface elevation after the third iteration calibration is shown in Figure 7e,f, and the laser ranging error after the calibration is shown in Table 2.
Table 2 shows the coordinates of the laser footprints of GF-7 beam 1 (laser time code is 203685561.00) and beam 2 (laser time code is 203685553.00) after each calibration iteration and true surface elevation value. Note: the true surface elevations of the laser footprints calculated after each iterative calibration were the same because the selected surface was very flat, as shown in Figure 1b,c. Although the position of the footprint changed after iterative calibration, the surface elevation remained the same across the plane. As shown in Table 2, the calculated coordinates of the same laser footprint were completely consistent after the second and third iterations, which indicates that the GF-7 laser pointing and ranging calibration had converged after two iterations using AW3D30 DSM as the reference terrain. Figure 7 also illustrates the same result; for example, the surface plot of the sum of laser elevation residuals was completely consistent after the second and third iterations of calibration.
From the results of the first and second iterative calibrations shown in Figure 7 and Table 2, it can be seen that the laser pointing angle could still be changed after the pointing calibration was repeated by correcting the laser ranging error. However, whether the laser pointing accuracy was improved after the change needs further verification. We used the GF-7 laser pointing after calibration based on LGDs as the true value [10] to verify the GF-7 laser pointing accuracy before and after each iterative calibration. Table 3 shows the difference between the initial pointing of the GF-7 laser, the laser pointing and ranging after iterative calibration based on the AW3D30 DSM, and the real laser pointing and ranging. Among them, Δ α x is the difference between the laser pointing angle α x after the calibration and the real pointing angle. Δ α y is the pointing difference of the pointing angle α y . Δ ρ is the difference between the laser ranging error after the calibration and the real ranging error.
The results in Table 3 show that the calibration method proposed in the paper had a fast convergence speed, and the results converged directly after two iterations. The main reason for this is that the input ranging error of beam 1 changed from 0 to −0.76 m, and that of beam 2 changed from 0 to −0.36 m (as shown in Table 2) during the second iteration. When performing the third iterative calibration, the ranging error of beam 1 changed by 0.11 cm, and the ranging error of beam 2 only changed by 0.08 cm. The change in the ranging error in the third iteration was much smaller than the change in the ranging error from the first to the second calibration. According to the error of the ranging error in Table 3, the ranging error was very close to the true value after the second iteration. Therefore, when there was a small change in the laser ranging, the result of the laser pointing calibration would not change again.
After the first iterative calibration, the laser pointing angle quickly approached the true value, in which the pointing angle accuracy of beam 1 was improved by more than 12 times, from the initial −0.038 degrees to 0.00316 degrees. The beam 2 pointing angle accuracy was improved by an order of 13 from the initial −0.044 degrees to 0.0033 degrees. After iterative calibration and convergence using AW3D30 DSM, the pointing accuracy of GF-7 beam 1 reached 0.00116 degrees (4.2 arcsec), the beam 2 pointing accuracy reached 0.0013 degrees (4.7 arcsec), and the ranging accuracies became 0.03 m and −0.06 m, respectively.

3.2. Elevation Accuracy Validation after Calibration

The accuracy of satellite-derived laser elevation measurements is a topic that has gained extensive interest in various industries. In this study, high-precision GCPs were used as elevation data to verify the elevation accuracy of GF-7 lasers after iterative calibration based on the AW3D30 DSM. First, taking the laser pointing and ranging parameters after the first and second calibration iteration as the input parameters, the coordinates of the 24 laser footprints of beams 1 and beam 2 were calculated, as shown in Figure 2. The four nearest GCPs were chosen to calculate the actual laser footprint elevation using bilinear interpolation. Then, the elevation differences between all laser elevations of beam 1 and beam 2 (as shown in Figure 2) and the actual elevation were calculated to verify the GF-7 laser elevation accuracy after calibration, and the results are shown in Table 4 and Table 5. Because of the large initial pointing error of the GF-7 satellite laser, the GCP array measured by RTK was far from the initial laser footprint position; thus, it was difficult to verify the elevation accuracy of the GF-7 laser prior to calibration. The results in Table 3 show that the pointing and ranging error after the third calibration iteration were consistent with the results after the second iteration. Therefore, the elevation accuracy of the GF-7 laser after the third iteration was no longer verified here.
Note that the GCP elevations in Table 4 and Table 5 were different in the first and second iterations. The reason is that the laser pointing angle changed after the first and second iterations during calibration. This means that the positions of the laser footprints also changed. Thus, the GCPs around the laser footprints changed accordingly, resulting in different GCP elevations. To show the changes more clearly in the elevation accuracy of the GF-7 laser beam 1 and beam 2 after the iterative calibration, Table 4 and Table 5 are presented graphically in Figure 8.
The above figure and table results show that the iterative geometry calibration method for laser altimeter pointing and ranging separation performed without field detectors, as proposed in this paper, obviously improved the elevation accuracy of the GF-7 laser. Among these, the elevation accuracy of beam 1 was improved from 0.35 ± 0.27 m to 0.12 + 0.22 m after the first iteration, and the elevation accuracy of beam 2 was improved from 0.37 ± 0.46 m to −0.11 ± 0.25 m. It is worth noting that the GCP elevation in Table 4 and Table 5 used to verify the laser footprint elevation did change after the first and second iterations of calibration. The reason is that the laser pointing of the GF-7 laser changed after the second iterative calibration, resulting in a change in the surface position of the laser footprint. Therefore, the GCPs used to verify the elevation of the laser footprint were also different.
After iterative calibration using the AW3D30 DSM, the maximum elevation error of GF-7 laser beam 1 was approximately 0.34 m, and the maximum elevation error of beam 2 was approximately −0.36 m. In addition, beam 1 still had a systematic error of 0.12 m, and beam 2 also had a systematic error of −0.11 m. The main reason is that beam 1 and beam 2 of the GF-7 lasers still had pointing and ranging errors after calibration. The maximum pointing error of beam 1 was 4.2 arcsec after calibration, and that of beam 2 was 4.7 arcsec. If a satellite laser pointing error of 1 arcsec will lead to a 2.5 m positioning error in the laser footprint [24], there remains a positioning error of 10.5 m in the laser footprint of beam 1 and a positioning error of 11.8 m in beam 2. Because the grassland in Sonid Youqi, Inner Mongolia, is not as flat as the cement road, the changes to the laser-derived elevation accuracy resulting from a laser footprint positioning error of approximately 11 m cannot be estimated. Therefore, it is reasonable that an elevation error of 0.34 m persists after calibrating the GF-7 satellite.

4. Discussion

When the AW3D30 DSM was used as the reference terrain, the GF-7 satellite laser footprint still had a positioning error of 11 m after the iterative calibration of the pointing and ranging separation. We consider that the error was mainly caused by the accuracy of the reference terrain itself. The resolution of the AW3D30 DSM is 30 m, and the elevation accuracy is approximately 2 m [17]. The 30 m resolution of the reference terrain data resulted in the spacing of the laser pointing grids being only divided to 0.003°. The low elevation accuracy of the DSM also limited the laser pointing accuracy after iterative calibration. For all the above reasons, higher-precision reference terrain data will be used to discuss the calibration accuracy of this method.
The DSM image of Kuqa City in Xinjiang generated by the airborne LiDAR data was used as the reference terrain and had a horizontal 1 m resolution. The methods discussed in this paper were used to repeat the pointing and ranging calibration of the GF-7 beams 1 and 2. Here, the first iterative calibration based on LiDAR took the GF-7 laser pointing angle after the second iterative calibration based on AW3D30 (as presented in Table 3) as the initial pointing, and the grid range was set to three times the grid spacing of the second iterative calibration (which was based on the AW3D30 DSM). The grid spacing was 0.0001 degrees (corresponding to the 1 m resolution of the LiDAR data). The surface plots of the sum of the laser elevation residuals during the first iterative calibration based on the LiDAR data are shown in Figure 9a,b. Similar to the ranging calibration method based on the AW3D30 DSM, the two laser footprint elevations in Figure 1b,c were calculated using the pointing angles derived following after the LiDAR calibration, and the elevation difference between the laser elevation and the true surface elevation was used as the laser ranging error. The ranging error of GF-7 laser beam 1 was −0.65 m, and beam 2 was −0.28 m. Subsequently, the second iterative calibration experiment used the first LiDAR calibration result as the input to completely repeat the previous calibration experiment, including the grid pointing range and spacing settings. The surface plots of the sum of the laser elevation residual were calculated as shown in Figure 9c,d.
In Figure 9, there is still an obvious minimum value in the surface plot of the sum of the elevation residuals, but the grid spacing was 30 times smaller than the grid spacing set during the AW3D30 DSM calibration. That is, with the improvement in the accuracy of the reference terrain, the laser pointing optimal solution still existed. According to Figure 9, however, the surface plots of the sum of the elevation residuals after the second calibration were basically the same as the surface plots of the sum of the elevation residuals after the first calibration. This indicates that the GF-7 laser converged after the first calibration when using LiDAR data for reference. Similarly, the GF-7 laser pointing angle and ranging errors derived from the LGD calibration were used as the true values to verify the pointing and ranging accuracies derived from the LiDAR data calibration. The results are shown in Table 6.
The calibration results in Table 6 also show that the calibration based on LiDAR data converged after one iteration. The reason is that the laser pointing and ranging accuracy were already very high after the first iterate when using LiDAR. Compared with the calibration results based on the AWD3D30 DSM, the method discussed in this paper uses LiDAR data as the reference terrain to further improve the laser pointing and ranging accuracy of the GF-7 satellite. Using this method and LiDAR data nearly doubled the accuracy. Among these, the pointing accuracy of GF-7 beam 1 was −0.00056 degrees (−2 arcsec) and the ranging accuracy was −2 cm after calibration, while the pointing accuracy of GF-7 laser beam 2 was 0.0006 degrees (2.2 arcsec) with a ranging error of −1 cm after calibration. If the pointing error was roughly converted into the GF-7 satellite laser footprint positioning error, the positioning error of beam 1 after calibration was 5 m, and that of beam 2 was 5.5 m.
The above results verified that the laser pointing and ranging accuracy of the GF-7 satellite was improved by using the LiDAR data as the reference terrain. However, whether the GF-7 satellite laser pointing and ranging accuracy can improve the laser elevation accuracy needs further verification. Based on the LiDAR data, the results of the second iterative calibration were consistent with the results of the first iterative calibration. We directly used the GF-7 satellite laser pointing and ranging parameters after the first iterative calibration to calculate the surface coordinates of the same GF-7 laser footprints in Table 4 and Table 5. Finally, the calculated laser elevation was verified with GCPs, and the results are shown in Table 7.
The line graphs shown in Figure 10 are used to visually display the GF-7 satellite laser elevation accuracies after the first and second calibrations using the AW3D30 DSM and the first iteration using the LiDAR reference data.
The results in Figure 10 show that when higher-precision reference terrain data are used for calibration, the elevation accuracy of the GF-7 satellite laser can be further improved. Compared with the elevation error of the GF-7 satellite laser by the AW3D30 DSM, the systematic laser elevation error after calibration using LiDAR data was suppressed and basically fluctuated around the 0 value. Additionally, the trend of the laser elevation accuracy was completely different after calibration with ALOS, as shown in Figure 10a,b. This is because there was still approximately a 10 m positioning error in the GF-7 laser footprints after calibration. In addition, the larger elevation errors for the individual laser footprints that remained after calibration using the AW3D30 DSM became relatively small after the LiDAR-based calibration, such as the 19th point of GF-7 laser beam 1. After calibration based on the LiDAR data, the elevation accuracy of GF-7 beam 1 was improved from 0.12 ± 0.22 m to 0.06 ± 0.11 m and beam 2 was improved from −0.11 ± 0.25 m to −0.05 ± 0.13 m. Based on the calibration method proposed in this paper, the GF-7 satellite laser elevation accuracy was approximately 0.18 m after being calibrated with the high-precision LiDAR data.
In summary, when using higher-precision LiDAR data for terrain matching, the calibration method proposed in the article can double the pointing and ranging accuracy of the satellite laser. The elevation accuracy of the GF-7 satellite was 0.18 m after using the LiDAR data for calibration. The elevation accuracy of the GF-7 satellite laser after being calibrated by LGDs was 0.1 m [10]. The calibration accuracy of the method presented here was lower than that based on LGDs. The main reason is that there was a positioning error of approximately 5 m in the laser pointing angle after applying the calibration method discussed in this paper.

5. Conclusions

This work presents an iterative geometry calibration method for satellite laser altimeter pointing and ranging separation that is not dependent on field detectors. The satellite laser pointing calibration is first completed based on the terrain matching method, and then the satellite laser ranging error calibration is achieved by accurately measuring the ground elevation of the laser footprint. The calibrated laser pointing and ranging errors are used to iteratively calibrate the satellite laser pointing and ranging errors. Then, this work also discussed the calibration accuracy of our calibration method when using the AW3D30 DSM and high-precision LiDAR data. Finally, after the laser calibration of the GF-7 satellite using the method outlined in this paper is complete, the pointing and ranging error accuracy of the GF-7 satellite laser was verified by using the calibration results from LGDs, and the GCPs were used to verify the elevation accuracy of the GF-7 satellite laser. According to the experimental results in this article, the following four main conclusions are drawn.
(1)
The iterative calibration method for the separation of laser pointing and ranging proposed in this paper solves the ranging calibration problem that cannot be solved using traditional terrain matching laser calibration methods. The method proposed here can further improve the calibration accuracy of satellite lasers.
(2)
The convergence speed of this calibration method is fast. Typically, the satellite laser pointing angle converges after 2–3 iterations. Taking the AW3D30 DSM as the reference terrain data, the pointing angle accuracy of GF-7 beam 1 was 4.2 arcsec, and the ranging accuracy was 3 cm after the calibration. The pointing accuracy of the beam 2 was 4.7 arcsec, and the ranging accuracy was 6 cm.
(3)
High-precision DSM data can improve the accuracy of the calibration method. Based on the calibration results from the AW3D30 DSM as the reference terrain, when LiDAR data was used as the reference terrain data, the laser pointing and ranging accuracy of the GF-7 satellite was improved. The improved laser pointing accuracy of GF-7 satellite beam 1 was 2 arcsec, and the ranging accuracy was 2 cm. The GF-7 laser beam 2 pointing accuracy was 2.2 arcsec, and the ranging accuracy was 1 cm.
(4)
Using high-precision LiDAR data as the reference terrain data, the elevation accuracy of the GF-7 satellite laser was 0.18 m in flat terrain after the calibration, which was slightly lower than that of the calibration method based on the ground detector. However, it was still slightly lower than that of the calibration method based on LGDs.
In summary, the application and implementation of the calibration method presented in this paper can provide a new high-precision geometric calibration technique for satellite lasers when ground detector calibration cannot be performed. However, at present, we are only conducting calibration experiments on barren, sloped terrain in northern China. To further examine the applicability of this method, subsequent studies will specifically consider more varied terrain types with different slopes.

Author Contributions

Conceptualization, R.L. and J.X.; formal analysis, R.L. and C.X.; methodology, R.L.; software, C.X. and J.Z.; validation, R.L. and J.Z.; data curation, F.M. and X.Y.; writing—original draft preparation, R.L.; writing—review and editing, J.X.; visualization, F.M. and X.Y.; supervision, J.X. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly supported by the Research and Development of Forest Resources Dynamic Monitoring and Forest Volume Estimation Using LiDAR Remote Sensing (2020YFE0200800), the National Natural Science Foundation of China (41971426), and the Multi-beam Laser Terrain Detection Radar and Application Technology (D040105).

Data Availability Statement

The AW3D30 DSM data are available via the JAXA Earth Observation Research Center website (https://www.eorc.jaxa.jp/ALOS/en/aw3d30/data/index.htm, accessed on 20 November 2022); the GF-7 satellite laser data are available on the Natural Resources Satellite Remote Sensing Cloud Service Platform (http://sasclouds.com/chinese/normal/, accessed on 20 November 2022).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The GF-7 satellite laser data are archived at the Land Satellite Remote Sensing Application Center, Ministry of Natural Resources of P.R. China (LASAC). The application for free GF-7 satellite laser data from the institution requires two steps: (1) GF-7 satellite laser data query; (2) GF-7 satellite laser data request.
(1)
GF-7 Satellite Laser Data Query
The GF-7 satellite laser data query only supports the Chinese version of the Natural Resources Satellite Cloud Service Platform website (http://sasclouds.com/chinese/normal/, accessed on 20 November 2022), not the English version of the website. After accessing the Chinese version of the data website, users need to click “Image query” in the menu bar to enter the data query page, as shown in Figure A1. First, the user sets the data range to be retrieved in the Figure A1 area 1, which supports shp file import, latitude and longitude range setting, etc. The user then needs to select “SLA(scene)” or “SLA(strip)” in area 2 and set the cloud size in area 3. Finally, the dark blue button in area 3 is clicked to complete the data query. The queried data is displayed in the directory page of “query result”. In addition, the “Export Filename” button under the “Export” submenu in this directory can be used to export the queried GF-7 satellite laser data file name to a text file.
Figure A1. GF-7 satellite laser data query web page.
Figure A1. GF-7 satellite laser data query web page.
Remotesensing 14 05935 g0a1
(2)
GF-7 satellite laser data request
After obtaining the file name of the GF-7 satellite laser data, the data can be applied for free. However, the GF-7 satellite laser data application can only be obtained through telephone and email contact via LASAC. The telephone number for the data application is +86-01068412488 and the email address is [email protected]. The applied GF-7 satellite laser data are SLA03 level data, and the data format description can also be downloaded from the data inquiry website (http://sasclouds.com/chinese/applicationExtend/fileShare/application, accessed on 20 November 2022).

References

  1. Tang, X.; Xie, J.; Liu, R.; Huang, G.; Zhao, C.; Zhen, Y.; Tang, H.; Dou, X. Overview of the GF-7 Laser Altimeter System Mission. Earth Space Sci. 2019, 6, e2019EA000777. [Google Scholar] [CrossRef]
  2. Xie, J.; Huang, G.; Liu, R.; Zhao, C.; Dai, J.; Jin, T.; Mo, F.; Zhen, Y.; Xi, S.; Tang, H.; et al. Design and Data Processing of China’s First Spaceborne Laser Altimeter System for Earth Observation: GaoFen-7. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 1034–1044. [Google Scholar] [CrossRef]
  3. Luthcke, S.B.; Rowlands, D.D.; Williams, T.A.; Sirota, M. Reduction of ICESat systematic geolocation errors and the impact on ice sheet elevation change detection. Geophys. Res. Lett. 2005, 32, 312–321. [Google Scholar] [CrossRef] [Green Version]
  4. Luthcke, S.B.; Carabajal, C.C.; Rowlands, D.D. Enhanced geolocation of spaceborne laser altimeter surface returns: Parameter calibration from the simultaneous reduction of altimeter range and navigation tracking data. J. Geodyn. 2002, 34, 447–475. [Google Scholar] [CrossRef]
  5. Luthcke, S.B.; Rowlands, D.D.; Mccarthy, J.J.; Pavlis, D.E.; Stoneking, E. Spaceborne Laser-Altimeter-Pointing Bias Calibration from Range Residual Analysis. J. Spacecr. Rocket. 2000, 37, 374–384. [Google Scholar] [CrossRef]
  6. Rowlands, D.D.; Pavlis, D.E.; Lemoine, F.G.; Neumann, G.A.; Luthcke, S.B. The use of laser altimetry in the orbit and attitude determination of Mars Global Surveyor. Geophys. Res. Lett. 1999, 26, 1191–1194. [Google Scholar] [CrossRef] [Green Version]
  7. Magruder, L.A.; Webb, C.E.; Urban, T.J.; Silverberg, E.C.; Schutz, B.E. ICESat Altimetry Data Product Verification at White Sands Space Harbor. IEEE Trans. Geosci. Remote Sens. 2006, 45, 147–155. [Google Scholar] [CrossRef]
  8. Magruder, L.A.; Schutz, B.E.; Silverberg, E.C. Laser pointing angle and time of measurement verification of the ICESat laser altimeter using a ground-based electro-optical detection system. J. Geod. 2003, 77, 148–154. [Google Scholar] [CrossRef]
  9. Xie, J.; Tang, X.; Mo, F.; Tang, H.; Wang, Z.; Wang, X.; Liu, Y.; Tian, S.; Liu, R.; Xia, X. In-orbit geometric calibration and experimental verification of the ZY3-02 laser altimeter. Photogramm. Rec. 2018, 33, 341–362. [Google Scholar] [CrossRef]
  10. Tang, X.; Xie, J.; Mo, F.; Dou, X.; Li, X.; Li, S.; Li, S.; Huang, G.; Fu, X.; Liu, R.; et al. GF-7 dual-beam laser altimeter on-orbit geometric calibration and test verification. Acta Geod. Cartogr. Sin. 2021, 50, 384–395. [Google Scholar]
  11. Magruder, L.A.; Ricklefs, R.L.; Silverberg, E.C.; Horstman, M.F.; Suleman, M.A.; Schutz, B.E. ICESat Geolocation Validation Using Airborne Photography. IEEE Trans. Geosci. Remote Sens. 2010, 48, 2758–2766. [Google Scholar] [CrossRef]
  12. Liu, R.; Xie, J. Calibration of the laser pointing bias of the GaoFen-7 satellite based on simulation waveform matching. Opt. Express 2021, 29, 21844–21858. [Google Scholar] [CrossRef] [PubMed]
  13. Tang, X.; Xie, J.; Gao, X.; Mo, F.; Feng, W.; Liu, R. The In-Orbit Calibration Method Based on Terrain Matching with Pyramid-Search for the Spaceborne Laser Altimeter. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2019, 12, 1053–1062. [Google Scholar] [CrossRef]
  14. Nan, Y.; Feng, Z.; Liu, E.; Li, B. Iterative Pointing Angle Calibration Method for the Spaceborne Photon-Counting Laser Altimeter Based on Small-Range Terrain Matching. Remote Sens. 2019, 11, 2158. [Google Scholar] [CrossRef] [Green Version]
  15. Takaku, J.; Tadono, T.; Tsutsui, K.; Ichikawa, M. Validation of “Aw3D” Global DSM Generated from Alos Prism. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2016, 3, 25–31. [Google Scholar]
  16. Takaku, J.; Tadono, T.; Tsutsui, K. Generation of high resolution global DSM from ALOS PRISM. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. ISPRS Arch. 2014, 40, 243–248. [Google Scholar]
  17. Mo, F.; Xie, J.; Liu, Y. In Evaluation of AW3D30 Elevation Accuracy in China. In Proceedings of the 4th International Conference on Geographical Information Systems Theory, Applications and Management, Madeira, Portugal, 17–19 March 2018. [Google Scholar]
  18. Schutz, B.E. Laser footprint location (geolocation) and surface profiles. In Geoscience Laser Altimeter System (GLAS): Algorithm Theoretical Basis Document Version 3; Center for Space Research of the University of Texas at Austin: Austin, TX, USA, 2002; pp. 1–24. [Google Scholar]
  19. Abshire, J.B.; Sun, X.; Riris, H.; Sirota, M.; Palm, S.; Ketchum, E.A. Geoscience laser altimeter system (GLAS) for the ICESat mission-pre-launch performance. In Proceedings of the Conference on Lasers and Electro-Optics (CLEO ‘03), Baltimore, MD, USA, 6 June 2003. [Google Scholar]
  20. Abdalati, W.; Zwally, H.J.; Bindschadler, R.; Csatho, B.; Farrell, S.L.; Fricker, H.A.; Harding, D.; Kwok, R.; Lefsky, M.; Markus, T.; et al. The ICESat-2 Laser Altimetry Mission. Proc. IEEE 2010, 98, 735–751. [Google Scholar] [CrossRef]
  21. Ma, Y.; Li, S.; Weng, Y.; Zhou, H. Hydrostatic delay correction for satellite laser altimeter. Infrared Laser Eng. 2013, 42, 909–914. [Google Scholar]
  22. The Algorithm Theoretical Basis Document for the Atmospheric Delay Correction to GLAS Laser Altimeter Ranges. Available online: https://ntrs.nasa.gov/search.jsp?R=20130001652 (accessed on 2 June 2012).
  23. Jentzsch, G. Earth tides and ocean tidal loading. Lect. Notes Earth Sci. 2005, 66, 145–171. [Google Scholar]
  24. Zhao, P.; Li, S.; Ma, Y.; Liu, X.; Yang, J.; Yu, D. A new terrain matching method for estimating laser pointing and ranging systematic biases for spaceborne photon-counting laser altimeters. ISPRS J. Photogramm. 2022, 188, 220–236. [Google Scholar] [CrossRef]
Figure 1. The laser calibration area and GF-7 satellite laser calibration data distribution. (a) The GF-7 satellite laser pointing calibration area. (b) The GF-7 satellite laser beam 1 ranging calibration area. (c) The GF-7 satellite laser beam 2 ranging calibration area.
Figure 1. The laser calibration area and GF-7 satellite laser calibration data distribution. (a) The GF-7 satellite laser pointing calibration area. (b) The GF-7 satellite laser beam 1 ranging calibration area. (c) The GF-7 satellite laser beam 2 ranging calibration area.
Remotesensing 14 05935 g001
Figure 2. Elevation accuracy verification area and data points from the GF-7 satellite laser.
Figure 2. Elevation accuracy verification area and data points from the GF-7 satellite laser.
Remotesensing 14 05935 g002
Figure 3. The reference terrain data of laser pointing calibration in the study. (a) AW3D30 DSM. (b) LiDAR Data.
Figure 3. The reference terrain data of laser pointing calibration in the study. (a) AW3D30 DSM. (b) LiDAR Data.
Remotesensing 14 05935 g003
Figure 4. The GCPs around a GF-7 laser footprint. (a) The GCP array surrounding laser point 203168336.67 of GF-7 laser beam 1. (b) A picture of the GCP measurement.
Figure 4. The GCPs around a GF-7 laser footprint. (a) The GCP array surrounding laser point 203168336.67 of GF-7 laser beam 1. (b) A picture of the GCP measurement.
Remotesensing 14 05935 g004
Figure 5. Flow chart of the proposed satellite laser calibration method without a field detector in the study.
Figure 5. Flow chart of the proposed satellite laser calibration method without a field detector in the study.
Remotesensing 14 05935 g005
Figure 6. Diagram showing the geometric relationship between the satellite laser and the ground laser footprint.
Figure 6. Diagram showing the geometric relationship between the satellite laser and the ground laser footprint.
Remotesensing 14 05935 g006
Figure 7. The laser grid pointing angle and the sum of elevation residual plots after calibration based on AW3D30 DSM. (a,b) are the experimental results of beams 1 and 2 after the first iteration, respectively. (c,d) are the experimental results of beams 1 and 2 after the second iteration, respectively. (e,f) are the experimental results of beams 1 and 2 after the third iteration, respectively.
Figure 7. The laser grid pointing angle and the sum of elevation residual plots after calibration based on AW3D30 DSM. (a,b) are the experimental results of beams 1 and 2 after the first iteration, respectively. (c,d) are the experimental results of beams 1 and 2 after the second iteration, respectively. (e,f) are the experimental results of beams 1 and 2 after the third iteration, respectively.
Remotesensing 14 05935 g007
Figure 8. Elevation accuracy comparison of the GF-7 laser after iterative calibration. (a) The beam 1 result of the GF-7 laser. (b) The beam 2 result of the GF-7 laser.
Figure 8. Elevation accuracy comparison of the GF-7 laser after iterative calibration. (a) The beam 1 result of the GF-7 laser. (b) The beam 2 result of the GF-7 laser.
Remotesensing 14 05935 g008
Figure 9. The laser grid pointing angle and the sum of the elevation residual plots after calibration based on LiDAR data. (a,b) are the experimental results of beams 1 and 2 after the first iteration, respectively. (c,d) are the experimental results of beams 1 and 2 after the second iteration, respectively.
Figure 9. The laser grid pointing angle and the sum of the elevation residual plots after calibration based on LiDAR data. (a,b) are the experimental results of beams 1 and 2 after the first iteration, respectively. (c,d) are the experimental results of beams 1 and 2 after the second iteration, respectively.
Remotesensing 14 05935 g009
Figure 10. Comparison chart of the elevation accuracy of the GF-7 laser after calibration based on AW3D30 and LiDAR. (a) The beam 1 result of the GF-7 laser. (b) The beam 2 result of the GF-7 laser.
Figure 10. Comparison chart of the elevation accuracy of the GF-7 laser after calibration based on AW3D30 and LiDAR. (a) The beam 1 result of the GF-7 laser. (b) The beam 2 result of the GF-7 laser.
Remotesensing 14 05935 g010
Table 1. The basic information of all test laser data of the GF-7 satellite in this study.
Table 1. The basic information of all test laser data of the GF-7 satellite in this study.
TypeFile NameNumber of Laser Points
Beam 1Beam 2
Pointing calibrationKRN-GF7-20201014-005278-00000057033739
KRN-GF7-20201212-006160-0000006634910
KSC-GF7-20210607-008851-00000100083719
Ranging calibrationKSC-GF7-20200615-003419-000000349511
Elevation accuracy verificationSYC-GF7-20200609-003326-000000328975
SYC-GF7-20200614-003402-000000339753
SYC-GF7-20200619-003478-000000352946
KSC-GF7-20200624-003555-0000003615810
Table 2. The ranging error of the GF-7 laser based on three iterations of AW3D30 calibration.
Table 2. The ranging error of the GF-7 laser based on three iterations of AW3D30 calibration.
ExperimentBeam TypeSurface Coordinates of Laser FOOTPRINTSSurface Actual
Elevation (m)
Ranging Error (m)
Latitude (°)Longitude (°)Elevation (m)
First iterationBeam 139.24695115.8277412.1912.95−0.76
Beam 239.72343116.1092835.7836.14−0.36
Second iterationBeam 139.24692115.8279312.3012.95−0.65
Beam 239.72347116.1090835.8636.14−0.28
Third iterationBeam 139.24692115.8279312.3012.95−0.65
Beam 239.72347116.1090835.8636.14−0.28
Table 3. Laser pointing and ranging accuracy before and after three iterative calibrations based on AW3D30.
Table 3. Laser pointing and ranging accuracy before and after three iterative calibrations based on AW3D30.
GF-7 Laser Pointing TypeBeam 1Beam 2
Δ α x   ( ° ) Δ α y   ( ° ) Δ ρ   ( m ) Δ α x   ( ° ) Δ α y   ( ° ) Δ ρ   ( m )
Initial pointing0.031−0.038−1.010.107−0.044−1.26
Pointing after first calibration−0.000430.00316−0.080.000560.00330.14
Pointing after second calibration−0.000430.001160.030.000560.0013−0.06
Pointing after third calibration−0.000430.001160.030.000560.0013−0.06
Table 4. The elevation accuracy of GF-7 laser beam 1 after iterative calibration based on AW3D30.
Table 4. The elevation accuracy of GF-7 laser beam 1 after iterative calibration based on AW3D30.
Point
Number
First IterationSecond Iteration
Latitude (°)Longitude (°)Laser
Elevation
Latitude (°)Longitude (°)Laser
Elevation
Latitude (°)Longitude (°)
143.11004111.78070993.81993.350.46993.51993.180.33
243.06836111.76772999.31998.221.09999.01998.420.59
343.04752111.761241046.361045.910.451046.061045.870.19
443.02668111.754761043.751043.370.381043.451043.230.22
543.00584111.748281029.451028.870.581029.151028.880.27
642.99202112.130911074.511074.320.191074.211074.210
742.97117112.124431084.891084.430.461084.591084.450.14
842.95032112.117971092.551092.270.281092.251092.170.08
942.92947112.111491092.131091.920.211091.831091.790.04
1042.88124112.462411039.261039.81−0.551038.961038.760.2
1142.86038112.455951041.011040.580.431040.711040.510.2
1242.83952112.449491041.041040.450.591040.741040.450.29
1342.81866112.443041051.061050.690.371050.761050.680.08
1442.77694112.430141050.331050.100.231050.031049.950.08
1542.75607112.423691051.811051.450.361051.511051.260.25
1642.73521112.417261067.351066.960.391067.051066.870.18
1742.71435112.410821075.331074.930.41075.031074.850.18
1842.87520112.829201076.661076.420.241076.361076.250.11
1942.72912112.784041086.301085.930.371086.001086.71−0.71
2042.68738112.771171090.221089.780.441089.921089.850.07
2142.66651112.764741089.201088.890.311088.911088.94−0.03
2242.56215112.732641119.021118.760.261118.721118.74−0.02
2342.54128112.726231133.861133.690.171133.561133.530.03
2442.52040112.719821136.731136.400.331136.431136.340.09
Mean value (m)0.35--0.12
RMSE (m)0.27--0.22
Table 5. The elevation accuracy of GF-7 laser beam 2 after iterative calibration based on AW3D30.
Table 5. The elevation accuracy of GF-7 laser beam 2 after iterative calibration based on AW3D30.
Point NumberFirst IterationSecond Iteration
Latitude (°)Longitude (°)Laser
Elevation(m)
GCP
Elevation(m)
Elevation
Difference(m)
Laser
Elevation(m)
GCP
Elevation(m)
Elevation
Difference(m)
143.04254111.913061053.101053.50−0.411053.501053.65−0.15
243.02171111.906531056.061056.84−0.781056.471057.08−0.61
343.00087111.899991065.711066.36−0.651066.121066.58−0.46
443.00798112.288771044.741045.48−0.741045.151045.39−0.24
542.98714112.282241071.271071.65−0.381071.681071.77−0.10
642.96629112.275711078.841079.30−0.461079.241079.38−0.14
742.94545112.269181079.561080.14−0.571079.971080.25−0.28
842.92460112.262661081.051081.82−0.771081.461082.01−0.55
942.90376112.256141087.561088.05−0.491087.971088.16−0.19
1042.87648112.612991075.631074.011.621076.041075.710.33
1142.85562112.606481073.921074.32−0.401074.331074.38−0.05
1242.83477112.599971077.861078.32−0.451078.271078.38−0.11
1342.81391112.593471084.121084.59−0.471084.531084.66−0.14
1442.79305112.586971083.871084.27−0.401084.271084.33−0.06
1542.75134112.573971093.291093.78−0.491093.691093.84−0.15
1642.73048112.567481092.461092.87−0.421092.861092.97−0.11
1742.70962112.560991101.001101.40−0.401101.401101.67−0.26
1842.68877112.554501109.161109.65−0.491109.561109.78−0.22
1942.66791112.548011117.981117.920.061118.381118.030.36
2042.70360112.927411093.631094.06−0.431094.031093.720.31
2142.57840112.888551120.681121.04−0.361121.081121.12−0.04
2242.55754112.882081125.421125.80−0.391125.821125.590.23
2342.51580112.869161128.161128.44−0.281128.561128.59−0.03
2442.49494112.862701124.021124.25−0.231124.421124.370.05
Mean value (m)0.37--−0.11
RMSE (m)0.46--0.25
Table 6. Laser pointing angle and ranging accuracy of GF-7 laser calibrations based on AW3D30 and LiDAR.
Table 6. Laser pointing angle and ranging accuracy of GF-7 laser calibrations based on AW3D30 and LiDAR.
GF-7 Laser Pointing TypeBeam 1Beam 2
Δ α x   ( ° ) Δ α y   ( ° ) Δ ρ   ( m ) Δ α x   ( ° ) Δ α y   ( ° ) Δ ρ   ( m )
Initial pointing0.031−0.038−1.010.107−0.044−1.26
Pointing after second iteration based on AW3D30−0.000430.001160.030.000560.0013−0.06
Pointing after first iteration based on LiDAR0.00033−0.00056−0.020.000360.0006−0.01
Pointing after second iteration based on LiDAR0.00033−0.00056−0.020.000360.0006−0.01
Table 7. The elevation accuracy of the GF-7 laser beams 1 and 2 after iterative calibration based on LiDAR.
Table 7. The elevation accuracy of the GF-7 laser beams 1 and 2 after iterative calibration based on LiDAR.
Beam 1Beam 2
Point NumberLaser
Elevation
GCP
Elevation
Elevation DifferencePoint NumberLaser
Elevation
GCP
Elevation
Elevation
Difference
1993.51993.480.0311053.511053.58−0.07
2999.01998.630.3721056.481056.79−0.31
31046.061045.910.1531066.131066.41−0.29
41043.451043.280.1741045.161045.34−0.18
51029.151028.990.1551071.681071.70−0.02
61074.211074.23−0.0261079.251079.31−0.06
71084.591084.61−0.0271079.981080.16−0.18
81092.251092.26−0.0281081.471081.72−0.25
91091.831091.820.0091087.971088.04−0.07
101038.961038.840.12101076.041076.040.00
111040.711040.590.12111074.341074.300.04
121040.741040.680.07121078.281078.270.01
131050.761050.740.02131084.541084.55−0.01
141050.031049.870.16141084.281084.270.01
151051.511051.360.15151093.701093.80−0.10
161067.051066.930.12161092.871092.91−0.04
171075.031074.920.10171101.411101.340.07
181076.361076.300.06181109.571109.72−0.14
191086.001086.02−0.03191118.391118.160.23
201089.921089.95−0.04201094.041094.020.02
211088.901089.06−0.15211121.091121.040.05
221118.721118.80−0.08221125.831125.810.02
231133.561133.58−0.02231128.571128.550.03
241136.431136.46−0.03241124.431124.370.06
Mean value (m)0.06Mean value (m)−0.05
RMSE (m)0.11RMSE (m)0.13
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, R.; Xie, J.; Xu, C.; Zeng, J.; Mo, F.; Yang, X. A Separate Calibration Method of Laser Pointing and Ranging for the GF-7 Satellite Laser That Does Not Require Field Detectors. Remote Sens. 2022, 14, 5935. https://doi.org/10.3390/rs14235935

AMA Style

Liu R, Xie J, Xu C, Zeng J, Mo F, Yang X. A Separate Calibration Method of Laser Pointing and Ranging for the GF-7 Satellite Laser That Does Not Require Field Detectors. Remote Sensing. 2022; 14(23):5935. https://doi.org/10.3390/rs14235935

Chicago/Turabian Style

Liu, Ren, Junfeng Xie, Chaopeng Xu, Junze Zeng, Fan Mo, and Xiaomeng Yang. 2022. "A Separate Calibration Method of Laser Pointing and Ranging for the GF-7 Satellite Laser That Does Not Require Field Detectors" Remote Sensing 14, no. 23: 5935. https://doi.org/10.3390/rs14235935

APA Style

Liu, R., Xie, J., Xu, C., Zeng, J., Mo, F., & Yang, X. (2022). A Separate Calibration Method of Laser Pointing and Ranging for the GF-7 Satellite Laser That Does Not Require Field Detectors. Remote Sensing, 14(23), 5935. https://doi.org/10.3390/rs14235935

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