Next Article in Journal
A Spherical “Earth–Ionosphere” Model for Deep Resource Exploration Using Artificial ELF-EM Field
Next Article in Special Issue
Transformer-Based Global Zenith Tropospheric Delay Forecasting Model
Previous Article in Journal
Fine-Grained Ship Classification by Combining CNN and Swin Transformer
Previous Article in Special Issue
Modeling of Residual GNSS Station Motions through Meteorological Data in a Machine Learning Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved Independent Parameter Decomposition Method for Gaofen-3 Surveying and Mapping Calibration

1
Land Satellite Remote Sensing Application Center, Ministry of Natural Resources, Haidian Dist., Beijing 100048, China
2
Shanghai Institute of Satellite Engineering, Minhang Dist., Shanghai 201109, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(13), 3089; https://doi.org/10.3390/rs14133089
Submission received: 29 March 2022 / Revised: 18 June 2022 / Accepted: 24 June 2022 / Published: 27 June 2022
(This article belongs to the Special Issue Data Science and Machine Learning for Geodetic Earth Observation)

Abstract

:
The Gaofen-3 (GF-3) satellite can provide digital elevation model (DEM) data from its interferogram outputs. However, the accuracy of these data cannot be ensured without applying a surveying and mapping (SAM) calibration process, thus necessitating geometric and interferometric calibration technologies. In this paper, we propose an independent parameter decomposition (IPD) method to conduct SAM calibration on GF-3 data and generate high-accuracy DEMs. We resolved the geometric parameters to improve the location accuracy and resolved the interferometric parameters to improve the height accuracy. First, we established a geometric calibration model, analyzed the Range–Doppler (RD) model and resolved the initial imaging time error as well as the initial slant range error. Then, we established a three-dimensional reconstruction (TDR) model to analyze the height error sources. Finally, the interferometric phase error and baseline vector error were precisely estimated to ensure the vertical accuracy of the interferometric results by establishing the interferometric calibration model. We then used the GF-3 interferometric data derived on the same orbit in a north–south distribution to conduct the calibration experiment. The results show that the plane positioning accuracy was 5.09 m following geometric calibration, that the vertical accuracy of the interferometric results was 4.18 m following interferometric calibration and that the average absolute elevation accuracy of the derived DEM product was better than 3.09 m when using the GF-3 SAR data, thus confirming the correctness and effectiveness of the proposed GF-3 IPD calibration method. These results provide a technical basis for SAM calibration using GF-3 interferograms at the 1:50,000 scale in China.

1. Introduction

The Gaofen-3 (GF-3) satellite is the first civil C-band multipolarization synthetic aperture radar (SAR) satellite in China. This satellite has 12 imaging modes, including strip, spotlight and scanning modes. The image resolution ranges from 1 to 500 m, while the observation width ranges from 10 to 650 km. GF-3 can monitor global ocean and land information during all types of weather at all times and can expand the range of Earth observations through left–right attitude maneuvering. The satellite data can be used in many fields, such asoceanography, disaster reduction, water conservancy and meteorological and mapping research, to improve rapid response capabilities. This satellite has also filled the gap in civil autonomous high-resolution multipolarization SAR remote sensing data in China [1,2]. Topographic mapping is an important application direction in the field of remote sensing. Digital elevation model (DEM) data play important roles in national economic construction and scientific research [3,4]. Precise surveying and mapping (SAM) calibration technology has become one of the most effective means of global mapping; by incorporating SAR interferometry (InSAR) to acquire interferometric phase information of the land surface, high-precision elevation information is retrieved [5,6]. In 2000, the National Aeronautics and Space Administration (NASA), the National Imagery and Mapping Agency (NIMA) and the German Aerospace Center (DLR) jointly conducted the Shuttle Radar Topography Mission (SRTM), successfully completing approximately 80% of the proposed terrestrial topographic mapping mission on Earth; in addition, the absolute elevation accuracy of the DEM obtained by the SRTM was 16 m [7,8]. In 2016, the DLR acquired global TanDEM-X DEM data with an absolute elevation accuracy of 4 m; these data were improved over other global DEM data in their accuracy, thus introducing new opportunities for global topographic mapping [9]. In 2018, the DLR released global TanDEM-X DEM data at a 90 m grid size [10,11]. Therefore, obtaining global DEM data constitutes major engineering projects in the field of InSAR topographic mapping [11]. The abilities of high-precision satellite measurements in this application are important and depend not only on the design of the satellite platform index but also on the calibration of the geometric and interferometric parameters [6,12]. Geometric calibration steps improve the plane positioning accuracy of SAR images by calibrating the geometric parameters, while interferometric calibration technologies improve the vertical accuracy of the interferometric results by calibrating the acquired interferometric parameters [13,14]; together, these calibration steps can be used to determine the accuracy level of a global DEM. At present, many research teams in China and globally have published geometric and interferometric calibration methods that are applicable for mainstream SAR satellites [9,13,15].
After applying the geometric calibration method to the data obtained by the Japanese Advanced Land-Observing Satellite-Phased Array L-band SAR (ALSO-PALSAR), the geometric positioning accuracy of the strip mode reached 9.7 m [16]. A 5.74 m plane positioning accuracy was obtained for the Canada Radarsat-2 satellite after a geometric calibration method using corner reflectors (CRs) was applied [17]. German TerraSAR-X satellite data were subjected to high-precision geometric calibration with CRs and achieved a plane positioning result with an accuracy of 0.3 m [18]. In 2017, the geometric calibration precision of the GF-3 SAR satellite was studied, and the positioning accuracy was found to be better than 3 m [14]. In 2017, a multimode hybrid geometric calibration method was proposed for GF-3 images, and the geometric positioning accuracy was found to be better than 3 m [19]. In 2018, a sparse-control multiplatform geometric calibration method was proposed for TerraSAR-X/TanDEM-X and GF-3 images. The geometric positioning accuracy of the TerraSAR-X/TanDEM-X images found using this method was better than 3 m, while that of the GF-3 images was better than 7.5 m following calibration [15]. In 2020, a wide-area joint geometric calibration method was proposed for GF-3 SAR images, and the positioning error of the calibrated SAR images was better than 9 m [20]. Thus, this geometric calibration method has been successfully applied in engineering and scientific research in the past.
Although interferometric calibration technologies are complex, extensive research work has been conducted on this topic both in China and across the globe, and many interferometric calibration methods have been proposed. In 1999, an airborne interferometric calibration method based on a sensitivity equation was deduced using flat-ground geometric relationships [21]. In 2001, based on this research, an airborne interferometric calibration method was formulated based on the interferometric sensitivity equation under squint angle conditions [22]. In 2003, a three-dimensional reconstruction (TDR) model was introduced to derive the interferometric sensitivity equation, and a CR placement algorithm was proposed to minimize the condition numbers in sensitivity matrices [23]. In 2010, an airborne interferometric calibration model was proposed to accommodate the phase offset, baseline length and baseline angle [24]. Then, an airborne dual-antenna SAR interferometric calibration method based on the three-dimensional information of control points was proposed [25]. In 2011, a distributed satellite SAR interferometric altimetry model was proposed based on the TDR model [26]. An interferometric calibration optimization method for spaceborne distributed SAR systems was proposed to solve the problem of coupling the contribution of interferometric parameters to the elevation error [27]. In 2016, the interferometric calibration method derived based on the InSAR TDR model and sensitivity equation was verified using TanDEM-X data, thus proving the universality of this spaceborne SAR interferometric calibration method. However, in this technology verification, each error was not solved independently [28]. Table 1 further describes the current status of research on geometric and interferometric calibration methods.
To date, geometric calibration technologies have gradually matured, and interferometric calibration technologies have developed from airborne to spaceborne systems. Although extensive research has been conducted on this topic, most analyzed parameter errors exhibit strong coupling effects during the interference process, and no joint solution has been accurately distinguished, leading to inaccurate interferometric parameter corrections [28]. To solve this problem, a precise SAM calibration method is proposed herein based on independent parameter decomposition (IPD). We established a geometric calibration model and resolved the initial imaging time error as well as the initial slant range error to improve the location accuracy. On this basis, the interferometric phase error and baseline vector error were precisely estimated to ensure the vertical accuracy of the interferometric results by establishing the interferometric calibration model. The accuracy and reliability of the proposed method were verified by processing GF-3 SAR data obtained on the same orbit in six scenes, and high-accuracy DEM products were obtained, thus providing a technical reference for the global mapping missions of domestic SAR satellites.

2. Geometric Calibration Model

Due to the measurement errors of satellite orbit and imaging parameters, the plane positioning accuracy of SAR images is always restricted. Geometric calibration technologies are the key to achieving high-precision SAR image positioning information and are also the main method by which the generation of high-precision DEMs can be ensured. The Range–Doppler (R-D) model, the mainstream positioning model used for SAR geometry processing, is a rigorous model that conforms to the SAR imaging mechanism and describes the geometric relationships between satellite sensors and object points in the geocentric coordinate system. The R-D model is a nonlinear system of equations composed of the slant range equation, Doppler equation and Earth ellipsoid equation; these three equations can be expressed as follows [15]:
R 1 = S 1 P
f d o p = 2 λ S 1 P v s v t R 1
X t 2 + Y t 2 R e + H t 2 + Z t 2 R p 2 = 1
where R 1 is the distance between the main satellite sensor and the ground object, S 1 = X s , Y s , Z s T is the position vector of the satellite sensor, P = ( X t , Y t , Z t ) T is the position vector of the ground object, f d o p is the Doppler frequency, λ is the radar wavelength, v s = V x , V y , V z T is the velocity vector of the satellite sensor and v t is the velocity vector of the ground object. Because the ground object point is located on the Earth ellipsoid, the velocity vector of the ground object point in the geocentric coordinate system is 0; that is, v t = 0 . R e and R p are the major semiaxes and minor semiaxes of the Earth reference ellipsoid, respectively, and H t is the elevation of the ground object.
The position and velocity vectors of the satellite sensor can be described as functions of time using cubic polynomials, as shown in Formula (4).
X s = a 0 + a 1 t + a 2 t 2 + a 3 t 3 Y s = b 0 + b 1 t + b 2 t 2 + b 3 t 3 Z s = c 0 + c 1 t + c 2 t 2 + c 3 t 3 V x = a 1 + 2 a 2 t + 3 a 3 t 2 V y = b 1 + 2 b 2 t + 3 b 3 t 2 V z = c 1 + 2 c 2 t + 3 c 3 t 2
According to the R-D positioning model, the variables that are necessary for geometric positioning include the slant range (R), satellite sensor position vector ( S 1 = X s , Y s , Z s T ), satellite sensor velocity vector ( v s = V x , V y , V z T ), Doppler frequency ( f d o p ), radar wavelength ( λ ) and ellipsoid parameters ( R e / R p ). Among them, the radar wavelength is a known value that can be obtained according to the formula λ = c / f , where c is the speed of light and f is the frequency. The ellipsoid parameters can also be obtained after the ellipsoid model is selected.
The slant range can be obtained by measuring the temporal difference between the transmission and reception of the radar pulse; the slant range error causes the position of the field of view of the target to move along the Doppler line, thus resulting in a positioning error. In the geometric positioning model, the slant range can be expressed as a linear function of the image column number as follows:
R 1 = r 0 + m j j
where r 0 is the initial slant range, m j is the range pixel spacing and j is the range pixel coordinate. In Equation (5), only the effect of the initial slant range error is considered when assessing the positioning accuracy in the range direction.
The linear function between the azimuth direction time and image line number can be expressed as follows:
t = t 0 + t 1 i
where t 0 is the initial imaging time, t 1 is the azimuth time spacing and i is the line number of the image along the azimuth. Equation (6) considers only the effect of the initial imaging time error on the positioning accuracy in the azimuth direction.
In summary, we consider only the impacts of the initial slant range and imaging time on the plane positioning accuracy, so the geometric calibration parameters of interest in this paper are ( r 0 , t 0 ) . Then, the geometric calibration model is established according to the R-D model, and the geometric calibration parameters are estimated.
First, a conditional equation based on the slant range equation is established as follows:
F 1 = S 1 P r 0 + m j j
Next, a conditional equation based on the Doppler frequency is established as follows:
F 2 = 2 λ S 1 P v R 1
Then, the geometric calibration parameters are estimated using the nonlinear least squares method by considering several calibration points derived based on the conditional equation. The calibration parameter ( r 0 , t 0 ) is used as an unknown observation value, and the error equation is obtained by linearizing the observation value derivation; these equations can be expressed as follows:
V 1 = F 1 r 0 Δ r 0 + F 1 t 0 Δ t 0 + F 1 0 F 1 V 2 = F 2 r 0 Δ r 0 + F 2 t 0 Δ t 0 + F 2 0 F 2
where the true observation value is F 1 = F 2 = 0 . The observation errors F 1 0 and F 2 0 can be calculated by substituting approximate values of each undetermined geometric calibration parameter ( r 0 , t 0 ) into Equations (7) and (8), respectively. The error equation can be rewritten as a matrix and expressed as follows:
V = B X L
where V = V 1 V 2 , B = F 1 r 0 F 1 t 0 F 2 r 0 F 2 t 0 , X = Δ r 0 Δ t 0 and L = F 1 0 F 2 0 .
The coefficients in matrix B are derived from the partial derivatives of Equations (7) and (8) to obtain the geometric calibration parameter ( r 0 , t 0 ) , which can be expressed as follows:
F 1 r 0 = 1 F 1 t 0 = K R 2 Q 1 R M Q 2 R F 2 r 0 = 0 F 2 t 0 = Q 2 R
where K = 2 / λ , M = X v X s X t + Y v Y s Y t + Z v Z s Z t , and
Q 1 = X v a 1 + 2 a 2 t + 3 a 3 t 2 + X s X t 2 a 2 + 6 a 3 t + Y v b 1 + 2 b 2 t + 3 b 3 t 2 + Y s Y t 2 b 2 + 6 b 3 t + Z v c 1 + 2 c 2 t + 3 c 3 t 2 + Z s Z t 2 c 2 + 6 c 3 t Q 2 = X s X t a 1 + 2 a 2 t + 3 a 3 t 2 + Y s Y t b 1 + 2 b 2 t + 3 b 3 t 2 + Z s Z t c 1 + 2 c 2 t + 3 c 3 t 2
Finally, the normal equation can be established based on the initial value ( X 0 ) of the geometric parameter.
B T B X B T L = 0
The geometric correction parameter can then be corrected based on the normal equation; the result can be written as follows:
Δ X = B T B 1 B T L
The first geometric calibration parameter estimation can be obtained as follows:
X = X 0 + X
The geometric correction parameter solution is obtained in a gradual iterative approach process. First, the initial values of calibration parameters r 0 and t 0 are obtained; then, the error equation is established according to these initial values, and the geometric calibration parameter is corrected using the normal equation. Next, the first approximate value of the calibration parameter is obtained, and this value is taken as the initial value in the next iteration. When the calibration parameter gradually stabilizes after many iterations and the parameter correction is less than the specified threshold ( ε = 10 6 ), the iteration is said to have converged. At this time, the cumulative parameter correction result is considered to be the final error in the geometric calibration parameters.

3. Three-Dimensional Reconstruction Model

The TDR is the basic model applied to spaceborne InSAR. Its basic principle involves obtaining the three-dimensional coordinates of GCPs by combining the slant range equation and Doppler equation in the R-D positioning model through the interferometric equation [29].
Herein, we describe the geometric relationship of the InSAR data in geocentric space rectangular coordinates (O-XYZ), as shown in Figure 1, where S 1 and S 2 are satellites; R 1 and R 2 are the position vectors of the satellite to the GCP; P is the position vector of the GCP; R 1 and R 2 are the slant ranges of S 1 and S 2 , respectively; v is the velocity vector; f d o p is the Doppler frequency; λ is the radar wavelength; φ is the absolute phase; ϕ is the wrapped phase; B is the baseline vector; and Q is the antenna receiving mode. The Q values of 2 and 1 represent the repeated-orbit interferometric mode and the double-antenna interferometric mode, respectively.
We mainly focus on the cross-track interferometry (XTI) mode of the spaceborne InSAR. The slant range equation, the Doppler equation in the R-D model and the interferometric phase equation are combined to establish the TDR model, and the geometric relationships can be expressed as follows:
R 1 = S 1 P
f d o p = 2 λ S 1 P v s v t R 1
φ = ϕ + 2 k π = 2 π Q ( R 1 R 2 ) λ
According to the geometric relationships of the interferometric SAR, the GCP can be expressed in O-XYZ coordinates as follows:
P = S 1 + R 1 = S 1 + R 1 r ^
where R 1 is the look vector and r ^ is the unit look vector. In Equation (18), the look vector decomposition (LVD) method is generally used to obtain the GCP coordinates. In this work, we convert the look vector to the unit look vector during the calculation [29], and the unit look vector needs to be represented in the Madsen moving coordinate (MMC) system [30]. The MMC system is based on applying the phase center of the main antenna as the origin of the coordinate system; the velocity vector direction is set to be equal to the V axis; the velocity and baseline cross product direction are set to be equal to the W axis; and V, N and W are set to meet the requirements of the right-hand coordinate system. To more conveniently describe the system, we refer to the MMC system as VNW in this work. The terms r v , r n and r w are considered components of the unit look vector on the V, N and W axes, respectively, and r ^ can be written as follows:
r ^ = r v r n r w
In the LVD method, the look vector is converted from VNW coordinates to O-XYZ coordinates to obtain the three-dimensional coordinates of the GCPs, as shown in Equation (20) below:
P = S 1 + R 1 = S 1 + R 1 · T v n w 1 · r ^ v n w = S 1 + R 1 · v v v B v v B v v B v B λ f d o p 2 v R 1 2 b p v 1 + B 2 R 1 2 1 λ φ 2 π Q R 1 2 b v r v b p v ± 1 r v 2 r n 2
where T v n w 1 is the transformation matrix, b v is the baseline component in the velocity direction and b p v is the baseline component in the direction perpendicular to the velocity.
According to the TDR expression, the main factors that affect the precision of interferometric SAR data are the interferometric phase error and the baseline vector error. Due to coupling among interferometric parameter error sources, to improve the accuracy of the interferometric parameter calculations during the calibration process, we adopted the IPD method, thus allowing us to precisely estimate the interferometric phase error and baseline vector error, improve the vertical accuracy of interferometric results and ultimately obtain high-precision DEM products.

4. Interferometric Calibration Model

The method utilized in this paper is based on the transmission characteristics of the interferometric phase error and baseline vector error to the elevation error; these characteristics are used herein to establish the interferometric calibration model. In spaceborne InSAR topographic mapping, interferometric phase errors generally consist of random and systematic errors. The random errors are caused by incoherent sources in the InSAR system, and the corresponding phases are randomly distributed; these errors are not involved in the calibration process. The systematic errors consist of phase drift, phase synchronization, temperature and absolute phase deviation errors [12,27]. Because the phase drift, phase synchronization and temperature errors are not the main sources of interferometric phase errors, they are omitted from the calibration process in this work. In a spaceborne InSAR system, when the interferometric phase is filtered and unwrapped, fixed deviations still exist among the flat phase, the unwrapped phase and the absolute interferometric phase; this fixed deviation can be considered the absolute phase deviation [31]. Therefore, the main component of the interferometric phase error is the absolute phase deviation, calculated as the interferometric parameter in this paper.
First, the absolute interferometric phase of each point was obtained from the GCPs. Then, the unwrapped phase and flat phase were subtracted from the absolute phase. Finally, the fixed deviation of the phase data of the SAR image was compensated based on the mean estimation method. In this way, the interferometric phase error of each image was obtained as follows:
φ o f f = φ a b s φ f l a φ u n w = 2 π ( R 1 R 2 ) λ φ f l a φ u n w
where φ a b s is the absolute interferometric phase, φ u n w is the unwrapped phase, φ f l a is the flat phase and φ o f f is the interferometric phase error.
In spaceborne InSAR topographic mapping, the baseline vector error directly affects many factors, such as the elevation ambiguity and the transmission coefficient of the elevation error. The baseline is a three-dimensional vector, and to more conveniently express the baseline parameters, the track–cross–normal (TCN) coordinate system is generally adopted. In this system, the master antenna phase center is set as the origin of the coordinates, and the N axis is derived from the main antenna phase center in the direction pointing toward the geocenter; the cross-product direction between the N axis and the velocity vector is considered the positive direction of the C axis, and T, C and N are set to meet the requirements of the right-hand coordinate system [28]. The baseline vector is then decomposed into the initial baseline and time-related baseline rate in the TCN system. Considering that all calculations are projected to the zero-Doppler plane of the main image during the XTI process conducted in this paper, the baseline along-track direction is always zero ( B t = 0 ) [32]. The baseline vectors are calculated as follows:
B t = 0 B c = b c 0 + t b c v B n = b n 0 + t b n v
where B t , B c and B n are the baseline three-dimensional (3D) vectors of the TCN coordinate system; b c 0 and b c v are the initial baseline and baseline rate of baseline B c in the C direction; and b n 0 and b n v are the initial baseline and baseline rate of baseline B n in the N direction.
In this paper, when the baseline vector parameters are calibrated, we first start from the TDR model expression of the InSAR data; then, the sensitivity equation is established. Finally, we use the adjustment model to correct the interferometric parameters based on the least-squares criterion. To improve the accuracy of the baseline calculation, we adopt the iteration method. The root mean square error (RMSE) of the coordinate difference derived before and after the iteration is defined as the objective function. When the objective function is less than a specified threshold ( ε = 0.05 in this paper), iterative convergence occurs. Finally, high-precision baseline 3D vectors are obtained to ensure the generation of high-precision DEM products. The process by which the baseline vectors are precisely estimated is shown in Figure 2 [13].
First, TDR equations representing the three coordinate axis directions are derived from the spaceborne InSAR model; these equations can be expressed as follows:
F x F y F z = S x S y S z P x P y P z + R 1 a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 r v r n r w
where S x S y S z   T are the coordinates of the master antenna phase center, P x P y P z   T are the coordinates of the GCP, R 1 is the slant range between the phase center of the master antenna and the GCP and a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 is the transformation matrix.
Second, the partial derivative of the baseline vector is calculated based on the TDR equation, and the interferometric sensitivity equation is established; this equation can be expressed as follows:
F x X F y X F z X = S x X S y X S z X + R 1 X a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 r v r n r w + R 1 a 11 X a 12 X a 13 X a 21 X a 22 X a 23 X a 31 X a 32 X a 33 X r v r n r w + a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 r v X r n X r w X
The error equation is then established as follows:
V B = F x 1 b c 0 F x 1 b c v F x 1 b n 0 F x 1 b n v F y 1 b c 0 F y 1 b c v F y 1 b n 0 F x 1 b n v F z 1 b c 0 F z 1 b c v F z 1 b n 0 F x 1 b n v . . . . . . . . . . . . F x n b c 0 F x n b c v F x n b n 0 F x n b n v F y n b c 0 F y n b c v F y n b n 0 F y n b n v F z n b c 0 F z n b c v F z n b n 0 F z n b n v Δ b c 0 Δ b c v Δ b n 0 Δ b n v + ( F x 1 ) ( F y 1 ) ( F z 1 ) ( F x n ) ( F y n ) ( F z n ) F x 1 F y 1 F z 1 F x n F y n F z n
where Δ b c 0 and Δ b c v denote the correction of the initial baseline and the baseline rate in the C direction, respectively, and Δ b n 0 and Δ b n v denote the correction of the initial baseline and the baseline rate in the N direction, respectively. The observation error can be calculated as follows:
( F x 1 ) ( F z n ) T = P x 1 0 P x 1 P z n 0 P z n T
where P x i 0 P y i 0 P z i 0   T are the three-dimensional reconstruction coordinates of the targets, P x i P y i P z i   T are the coordinates of the known GCP, and F x 1 F z n   T = 0 is the true observed value.
In the interferometric calibration process, the observation accuracy usually varies due to different observation conditions; thus, the introduced errors are often different. The coherence coefficient is an important criterion used to evaluate the quality of SAR images. Therefore, in this work, the coherence coefficients of the GCPs are weighted, and the objective function is expressed as follows:
J min = V T W V
where W is the diagonal matrix. The diagonal matrix can be written as follows:
W = w 1 0 0 0 w 2 0 0 0 w n
where w i is the coherence coefficient of each GCP.
Finally, given the initial values of the baseline vector parameters, the normal equation can be established as follows:
A T W A X A T W L = 0
The baseline vector was corrected based on the normal equation, which could be denoted as follows:
Δ X = ( A T W A ) 1 A T W L
The first estimation of the baseline vector was obtained as follows:
X = X 0 + Δ X

5. Experiments and Analyses

5.1. Experimental Data

In this work, the research areas are located in southern Nanyang city in Henan Province and in northern Xiangyang city and western Jingzhou city in Hubei Province. The whole study area covered approximately 55 km in the east–west direction and 411 km in the north–south direction. The topographic relief is very small, and the topographic type is typical flat terrain. The geographical location of the study area is shown in Figure 3. Fifty-six GCPs were acquired from the differential global positioning system (DGPS) in the geodetic coordinate system. The plane positioning error was 15 cm, and the elevation measurement error was 20 cm. The interferometric phase error was less than 1.12°, allowing these data to be used for high-precision measurement calibration tasks. Road intersections, house corners and other feature points were selected as GCPs in this paper [33], and the detailed GCP distribution is shown by the green triangles in Figure 3b,c. Before using the Ice, Cloud and Land Elevation Satellite (ICESat) data, it was necessary to filter out some data with elevation errors of tens or even hundreds of meters due to the presence of clouds, trees or buildings. In the end, a total of 10,438 points were selected, and the detailed point distribution is shown using the red dots in Figure 3b,c.
We chose six FS1 GF-3 images obtained in the descending direction to conduct our interferometric calibration experiment. Table 2 presents the derived statistics on the experimental image parameters. The images are denoted as A–F from top to bottom. SAR images A, C and F were acquired on 29 November 2016, while SAR images B, D and E were acquired on 28 December 2016. In this paper, a total of four pairs of interferometric regions were formed. The number of interferometric pairs, as well as the initial baseline and baseline rates in the TCN coordinate system, are listed in Table 2. The polarization mode was HH, the incident angle was 42.8 degrees, the range resolution was 2248 m and the azimuth resolution was 2865 m. SRTM data at a 90 m grid size were selected as the external reference DEM data.

5.2. Results and Analysis

5.2.1. Plane Positioning Accuracy Analysis

The plane positioning accuracies of the SAR images were evaluated using an indirect positioning algorithm (IPA). First, the image coordinates and geodetic coordinates of the GCPs were obtained. Then, the geodetic coordinates of the GCPs were transformed into image coordinates with the IPA. Finally, the RMSEs of the azimuth and range directions of all GCPs were obtained, and the RMSEs of the image points were calculated to evaluate the plane positioning accuracies of the SAR images. The results are shown in Table 3.
Table 3 shows that the average plane positioning accuracy of the six GF-3 images before the geometric calibration process was 37.706 m. The Δ t and Δ r corrections associated with the geometric parameter errors of each image were obtained after the six images were geometrically calibrated, and the initial time and slant range were both compensated for. The indirect geometric positioning algorithm was used to evaluate the plane positioning accuracies of the calibrated SAR images, and the results are listed in Table 4.
From Table 4, it can be seen that the average plane accuracy of the six GF-3 images after the geometric calibration process was 5.09 m, the average positioning error of the ranges of the image points was 4.22 m (approximately 1.87 pixels) and the average positioning error of the azimuths of the image points was 2.77 m (approximately 0.96 pixels). These results show that the plane positioning accuracies of the SAR image after geometric calibration can ensure the vertical accuracy of the interferometric results in the subsequent calibration process.

5.2.2. Vertical Accuracy Analysis

The elevation information of the GCPs was calculated using the spaceborne SAR interferometric model, and the results were evaluated using the elevations of the known GCPs to obtain the vertical accuracies of the interferometric results, as shown in Table 5.
Table 5 shows that the average elevation accuracy of the four pairs of GF-3 data was 166.159 m before the interferometric calibration process. We performed an interferometric calibration experiment on four pairs of GF-3 data. First, the phase error and baseline vector error were corrected based on IPD. Then, the absolute interferometric phase and baseline vector were compensated for. Finally, the vertical accuracy of the interferometric results was evaluated using the spaceborne SAR interferometric model; the evaluation results are listed in Table 6.
As shown in Table 6, the average vertical accuracy of the four pairs of GF-3 data was 4.18 m following the calibration process, and this accuracy meets the requirements for 1:50000-scale topographic mapping in flat areas. The vertical accuracies of the A–B, B–C, C–D and E–F pairs were 6.31 m, 2.26 m, 4.45 m and 3.70 m, respectively.

5.2.3. Spatiotemporal Analysis

To verify the applicability and stability of the calibration parameters characterizing the same-orbit images obtained at different imaging times, we conducted calibration compensation verification experiments. The interferometric calibration results of the A–B pair were compensated using those of the E–F pair, which had the longest imaging time difference in the same orbit. Finally, the accuracy of the verification data was evaluated; the evaluation results are listed in Table 7.
As Table 7 shows, when the calibration parameters were used to compensate for the verification scenes, the average vertical accuracy of the interferometric results without control points was 4.75 m; this accuracy meets the standard requirements of 1:50,000-scale topographic map mapping above grade 3 in flat areas in China (≤6 m) [34].
To analyze the influence of the spatial distribution of GCPs on the interferometric calibration results, we selected a pair of interferometric data (pair B–C) as an example and designed four schemes characterizing different spatial distributions of GCPs. In schemes (a), (b) and (c), six points were selected, while in scheme (d), 11 points were selected. The GCPs in scheme (a) were distributed along the azimuth, those in scheme (b) were distributed along the range and those in schemes (c) and (d) were randomly distributed, as shown in Figure 4. The number of calibration points included in schemes (a) and (b) were the same, but their spatial distribution modes were differentiated to verify the influence of the spatial distribution mode of the calibration points on the interferometric calibration results. Schemes (c) and (d) had the same calibration point distribution modes, characterized by random distributions, but the numbers of calibration points were differentiated to verify the influence of the number of calibration points on the interferometric calibration results.
Under the guarantee of the plane positioning accuracy of the utilized SAR images, the spatial distribution and number of GCPs directly affect the accuracy of the interferometric calibration results; these effects were explored by establishing four calibration schemes. The vertical accuracies of the interferometric results derived after calibration were evaluated according to the spaceborne SAR interferometric model; the results are listed in Table 8.
As shown in Table 8, when the number of calibration points was consistent, as between schemes (a) and (b), the elevation accuracy of scheme (a) was found to be better than that of scheme (b). Thus, the interferometric calibration accuracy was higher when the GCPs were distributed along the azimuth than along the range. While scheme (c) contained six GCPs and scheme (d) contained 11 GCPs, the GCPs in both schemes were randomly distributed in space; we found that the elevation accuracy of scheme (d) was better than that of scheme (c). Therefore, when GCPs are randomly distributed, the greater the number of points, the higher the interferometric calibration accuracy.

5.2.4. DEM Accuracy Analysis

We proposed a precise SAM calibration method based on GF-3 data. We decomposed the plane and elevation error sources in the spaceborne InSAR system by establishing IPD models; these factors are also the main sources of DEM product errors. In this paper, TanDEM-X DEM and ICESat data were used to evaluate the absolute elevation accuracy of the analyzed DEM after calibration, and the results are listed in Table 9.
Focusing on interferometric pair A–B, following the calibration process, the absolute elevation accuracies of the derived DEM were better than 2.10 m when the TanDEM-X DEM was used and 2.63 m when ICESat data were used. Regarding interferometric pair B–C, the post calibration absolute elevation accuracies of the DEMs were better than 2.73 m when the TanDEM-X DEM was used and 2.78 m when the ICESat DEM was used. For interferometric pair C–D, the post calibration absolute elevation accuracies were better than 2.40 m when TanDEM-X DEM was applied and better than 3.09 m when ICESat was applied. Focusing on interferometric pair E–F, the absolute elevation accuracies of the calibrated DEM were better than 1.44 m when TanDEM-X DEM was used and better than 2.55 m when ICESat data were used.
To quantitatively analyze the elevation variations in the GF-3-derived DEM following the calibration process, we drew a transect in the flat area shown in Figure 5a,b according to the distribution of ICESat points. As shown in Figure 5a, this transect covers pairs A–B, B–C and C–D from the northwest to the southeast corners of the image, with a length of approximately 106 km. As shown by the black line (labeled n ) in the figure, the section map corresponds to Figure 6a. The section line shown in Figure 5b spans from the north to the south regions of the image with a length of approximately 45 km. As shown by the black line (labeled m ) in the figure, the section map corresponds to Figure 6b, and the elevation data of the GF-3 DEM, ICESat DEM and TanDEM-X DEM are represented by green, red and blue lines, respectively, in the figure.
It is obvious that the elevation data contained in the GF-3 DEM and TanDEM-X DEM are basically consistent at interferometric pairs A–B, C–D and E–F. However, the elevations at some ICESat laser points are abnormal, as shown in Figure 6a1–a3,b1,b2. The information in the figure shows that the ICESat laser points in these areas are mainly distributed in building construction areas and vegetation-covered areas, possibly leading to poor elevation results. Therefore, the accuracy of the monitored GF-3 DEM products derived by comparison to the TanDEM-X DEM in this area was higher than that of the ICESat DEM.
Regarding interferometric pair B–C, the elevations of the ICESat DEM and TanDEM-X DEM were basically consistent. However, compared to the elevations of the ICESat DEM and TanDEM-X DEM, some of the elevation data in the GF-3 DEM were obviously abnormal. For further analysis, the coherence distribution was calculated in this area, as shown in Figure 6c; the coherence distribution corresponding to transect m is shown in Figure 6d. According to the statistical results, the coherence of transect m was lower than 0.6, accounting for 26.5%. In other words, 26.5% of the points in this region were located in incoherent areas, mainly due to the abundant ground feature information, including water areas, vegetation areas, and residential and urban areas; this ground feature information is often the main factor leading to incoherence in satellite images. Moreover, large shadows and overlaps are generated in side-looking radar imaging, which may lead to the poor accuracy of the elevation data reflected in the GF-3 DEM obtained in these areas.

6. Conclusions

In this paper, we proposed an IPD method to support the application of SAM approaches using GF-3 SAR data. In the proposed IPD method, geometric and interferometric calibration models are separately utilized to prevent parameter coupling phenomena from arising. We decomposed the sources of geometric and interferometric parameter errors in the GF-3 SAR system by establishing IPD models, and the derived errors could also be considered the main factors affecting the resulting DEM products. Herein, the plane positioning accuracies of the SAR images, the vertical accuracies of the interferometric results and the spatiotemporal distribution characteristics and elevation accuracies of DEM products were analyzed. The correctness, reliability and practicability of the proposed method were verified. The results show that the average plane positioning accuracy of the SAR images was 5.09 m following the geometric calibration process, and the average post-calibration vertical accuracy of the interferometric results was 4.18 m. The corresponding average absolute elevation accuracy of the DEM products was better than 3.09 m. These results confirm that the method proposed herein can improve the elevation accuracies of Gaofen-3 DEM products, which met the 1:50000-scale topographic map SAM precision, or even higher precisions, in the flat areas of China. Moreover, there are still some problems to be further studied in this paper. Due to the lack of high-precision CRs, the method proposed cannot fully verify the on-orbit calibration of satellites. Follow-up research will mainly focus on this aspect.
The follow-up satellites to GF-3, i.e., the 1 m C-SAR 01 and 02, were launched on 23 November 2021 and on 7 April 2022. This ocean-monitoring SAR satellite constellation consists of three satellites. The 1 m C-SAR satellites were designed in consideration of their interferometric capabilities. An increasing number of interferograms can be provided by this new satellite constellation, and these novel interferometric data can also be obtained more easily. We hope that the proposed method will be useful in supporting the usage of these new satellites.

Author Contributions

T.L. and J.F. designed the paper and conducted the experiments described in the manuscript. Y.L. is the corresponding author and edited the paper. J.F. contributed the prototype concept of the paper. R.L. and Y.H. helped edit and improve the manuscript. J.L. helped process the test data and assess the results. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key R&D Programme of China under Grant 2021YFC3000405, the Civil Spaceflight Pre-Research Projects under Grant B0302, and the National Natural Science Foundation of China under Grant 41901303.

Data Availability Statement

The data utilized herein are not publicly available due to privacy restrictions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Q. System Design and Key Technologies of the GF-3 Satellite. Acta Geod. Cartogr. Sin. 2017, 46, 269–277. [Google Scholar]
  2. Sun, J.; Yu, W.; Deng, Y. Design of Working Modes and Payload of SAR for GF-3 Satellite. Spacecr. Eng. 2017, 26, 61–67. [Google Scholar]
  3. Liu, G.X.; Ding, X.L.; Li, Z.L.; Chen, Y.Q.; Li, Z.W. Experimental Investigation on DEM Generation through InSAR. Acta Geod. Cartogr. Sin. 2001, 30, 336–342. [Google Scholar]
  4. Li, Z.; Li, P.; Ding, D.; Wang, H. Research Progress of Global High Resolution Digital Elevation Models. Geomat. Inf. Sci. Wuhan Univ. 2018, 43, 174–189. [Google Scholar]
  5. Tang, X.; Li, T.; Gao, X.; Chen, Q.; Zhang, X. Research on Key Technologies of Precise InSAR Surveying and Mapping Application Using Automatic SAR Imaging. Acta Geod. Cartogr. Sin. 2018, 47, 730–740. [Google Scholar]
  6. Tao, L.; Xinming, T.; Xiaoming, G.; Qianfu, C.; Xiang, Z. Analysis and outlook of the operational topographic surveying and mapping capability of the SAR satellites. Acta Geod. Cartogr. Sin. 2021, 50, 891–904. [Google Scholar]
  7. Bamler, R. The SRTM Mission-A World-Wide 30 m Resolution DEM from SAR Interferometry in 11 Days. Photogramm. Week 1999, 1999, 145–154. [Google Scholar]
  8. Peng, L.; Zhenhong, L.; Chuang, S.; Jingnan, L. Quality Evaluation of 1 Arc Second Version SRTM DEM in China. Bull. Surv. Mapp. 2016, 9, 24–28. [Google Scholar]
  9. Gonzalez, J.H.; Bachmann, M.; Krieger, G.; Fiedler, H. Development of the TanDEM-X Calibration Concept: Analysis of Systematic Errors. IEEE Trans. Geosci. Remote Sens. 2010, 48, 716–726. [Google Scholar] [CrossRef] [Green Version]
  10. Wessel, B.; Huber, M.; Wohlfart, C.; Marschalk, U.; Kosmann, D.; Roth, A. Accuracy assessment of the global TanDEM-X Digital Elevation Model with GPS data. ISPRS J. Photogramm. Remote Sens. 2018, 139, 171–182. [Google Scholar] [CrossRef]
  11. Rizzoli, P.; Martone, M.; Gonzalez, C.; Wecklich, C.; Tridon, D.B.; Bräutigam, B.; Bachmann, M.; Schulze, D.; Fritz, T.; Huber, M.; et al. Generation and performance assessment of the global TanDEM-X digital elevation model. ISPRS J. Photogramm. Remote Sens. 2017, 132, 119–139. [Google Scholar] [CrossRef] [Green Version]
  12. Li, T.; Tang, X.; Gao, X.; Chen, Q. Spaceborne InSAR Topographic Surveying and Mapping Error Sources Analysis. Bull. Surv. Mapp. 2017, 11, 37–41. [Google Scholar]
  13. Fan, J.; Li, T.; Zuo, X.; Chen, Q.; Zhang, X.; Lu, J. Interferometric calibration method for spaceborne SAR based on independent parameter decomposition. Acta Geod. Cartogr. Sin. 2019, 48, 737–746. [Google Scholar]
  14. Chibiao, D.; Jiayin, L.; Bin, L.; Xiaolan, Q. Preliminary exploration of systematic geolocation accuracy of GF-3 SAR satellite system. J. Radar 2017, 6, 11–16. [Google Scholar]
  15. Lv, G.; Tang, X.; Ai, B.; Li, T.; Chen, Q. Hybrid Geometric Calibration Method for Multi-platform Spaceborne SAR Image with Sparse GCPs. Acta Geod. Cartogr. Sin. 2018, 47, 986–995. [Google Scholar] [CrossRef] [Green Version]
  16. Shimada, M.; Isoguchi, O.; Tadono, T. PALSAR Radiometric and Geometric Calibration. IEEE Trans. Geosci. Remote Sens. 2009, 47, 3915–3932. [Google Scholar] [CrossRef]
  17. Zhou, X.; Zeng, Q.; Jiao, J. Field calibration and validation of Radarsat-2. In Proceedings of the 2013 IEEE International Geoscience and Remote Sensing Symposium, Melbourne, Australia, 21–26 July 2013; pp. 4451–4454. [Google Scholar]
  18. Cong, X.; Balss, U.; Eineder, M. Imaging Geodesy—Centimeter-Level Ranging Accuracy With TerraSAR-X: An Update. IEEE Geosci. Remote Sens. Lett. 2012, 9, 948–952. [Google Scholar] [CrossRef] [Green Version]
  19. Zhao, R.; Zhang, G.; Deng, M. Geometric Calibration and Accuracy Verification of the GF-3 Satellite. Sensors 2017, 17, 1977. [Google Scholar] [CrossRef] [Green Version]
  20. Ding, L.; Tao, Q.; Li, T.; Chen, Q.; Chen, Y. A joint geometric calibration technique for GF-3 SAR image in wide area. Acta Geod. Cartogr. Sin. 2020, 49, 598–610. [Google Scholar]
  21. Zink, M.; Geudtner, D. Calibration of the Interferometric X-SAR System on SRTM. In Proceedings of the IEEE 1999 International Geoscience and Remote Sensing Symposium, Hamburg, Germany, 28 June–2 July 1999; Volume 1, pp. 227–229. [Google Scholar]
  22. Mallorqui, J.J.; Bara, M.; Broquetas, A. Calibration requirements for airborne SAR interferometry. Eur. Remote Sens. 2000, 4173, 267–278. [Google Scholar]
  23. Wang, Y. Studies on Calibration Model and Algorithm for Airborne Interferoemtric SAR; Insititute of Electrics, Chinese Academy of Sciences: Beijing, China, 2003. [Google Scholar]
  24. Jin, G.; Zhang, W.; Xiang, M. A New Calibration Algorithm of Interferometric Parameters for Dual-antenna Airborne InSAR. Acta Geod. Cartogr. Sin. 2010, 39, 76–81. [Google Scholar]
  25. Zhang, W.; Xiang, M.; Wu, Y. Using Control Points’ 3D Information to Calibrate the Interferometric Parameters of Dual-antenna Airborne InSAR Systems. Acta Geod. Cartogr. Sin. 2010, 39, 370–377. [Google Scholar]
  26. Hu, J.; Hong, J. Systematic Model and Errors Analysis for Distributed Satellites SAR System Based on Three-dimensional Reconstruction Model. J. Grad. Univ. Chin. Acad. Sci. 2011, 28, 489–497. [Google Scholar]
  27. Zhang, Y. Study on Error Analysis and DEM Precision Improvement Methods of Spaceborne Distributed InSAR; National University of Defense Technology: Changsha, China, 2011. [Google Scholar]
  28. Wu, D. Modeling and Methodology for Airborne and Spaceborne InSAR Calibration; Southwest Jiaotong University: Sichuan, China, 2016. [Google Scholar]
  29. Hua, F. Elevation Inversion through Multi-Baseline InSAR; China University of Mining and Technology: Xuzhou, China, 2014. [Google Scholar]
  30. Madsen, S.N.; Zebker, H.A.; Martin, J. Topographic Mapping Using Radar interferometry: Processing Techniques. IEEE Trans. Geosci. Remote Sens. 1993, 31, 246–256. [Google Scholar] [CrossRef]
  31. Wang, Q. Research on High-Efficiency and High-Precision Processing Techniques of Spaceborne Interferometric Synthetic Aperture Radar; National University of Defense Technology: Changsha, China, 2012. [Google Scholar]
  32. Krieger, G.; Moreira, A.; Fiedler, H. TanDEM-X: A satellite formation for high-resolution SAR interferometry. IEEE Trans. Geosci. Remote Sens. 2007, 45, 3317–3341. [Google Scholar] [CrossRef] [Green Version]
  33. Wang, T.; Zhang, G.; Li, D.; Jiang, W.; Tang, X.; Liu, X. Comparison between Plane and Stereo Block Adjustment for ZY-3 Satellite Images. Acta Geod. Cartogr. Sin. 2014, 43, 389–395. [Google Scholar] [CrossRef]
  34. National Administration of Surveying, Mapping and Geoinformation. Digital Products of Fundamental Geographic Information 1:50001:100001:25000 1:500001:100000 Digital Elevation Models; Surveying and Mapping Press: Beijing, China, 2010. [Google Scholar]
Figure 1. Geometric relationships of the spaceborne InSAR.
Figure 1. Geometric relationships of the spaceborne InSAR.
Remotesensing 14 03089 g001
Figure 2. Methodological flow of the spaceborne SAR interferometric calibration process.
Figure 2. Methodological flow of the spaceborne SAR interferometric calibration process.
Remotesensing 14 03089 g002
Figure 3. The study area considered in this paper: (a) the geographical location of the study area, (b) enlarged views of areas A to D, and (c) enlarged views of areas E to F. In (b,c), the white wireframes indicate the coverage of the SAR images, the green triangles denote GCPs, and the red points are the ICESat track.
Figure 3. The study area considered in this paper: (a) the geographical location of the study area, (b) enlarged views of areas A to D, and (c) enlarged views of areas E to F. In (b,c), the white wireframes indicate the coverage of the SAR images, the green triangles denote GCPs, and the red points are the ICESat track.
Remotesensing 14 03089 g003
Figure 4. Four GCP layout schemes considered herein: (a) denotes the distribution of the 6 points along the azimuth direction, (b) denotes the distribution of the 6 points along the range direction (c) and (d) denote the distribution of 6 and 11 points along the random direction, respectively. Where the blue triangles represent the GCPs.
Figure 4. Four GCP layout schemes considered herein: (a) denotes the distribution of the 6 points along the azimuth direction, (b) denotes the distribution of the 6 points along the range direction (c) and (d) denote the distribution of 6 and 11 points along the random direction, respectively. Where the blue triangles represent the GCPs.
Remotesensing 14 03089 g004
Figure 5. DEM products obtained after the calibration process: (a) the DEM results corresponding to interferometric pairs A–B, B–C and C–D and (b) DEM results corresponding to interferometric pair E–F. In (a,b), the green triangles show the GCPs and the red points show the ICESat laser points.
Figure 5. DEM products obtained after the calibration process: (a) the DEM results corresponding to interferometric pairs A–B, B–C and C–D and (b) DEM results corresponding to interferometric pair E–F. In (a,b), the green triangles show the GCPs and the red points show the ICESat laser points.
Remotesensing 14 03089 g005
Figure 6. Elevation distribution maps of transects m and n: (a) the elevation change along section line m in interferometric pairs A–B, B–C and C–D; (b) the elevation change along section line n in interferometric pair E–F; (c) the coherence distribution maps of interferometric pair B–C; (d) the coherence distribution results of section line m. Panels (a1a3) and (b1,b2) contain enlarged images of regions with abnormal ICEsat points; the labels Pa1–Pb5 represent abnormal ICEsat points.
Figure 6. Elevation distribution maps of transects m and n: (a) the elevation change along section line m in interferometric pairs A–B, B–C and C–D; (b) the elevation change along section line n in interferometric pair E–F; (c) the coherence distribution maps of interferometric pair B–C; (d) the coherence distribution results of section line m. Panels (a1a3) and (b1,b2) contain enlarged images of regions with abnormal ICEsat points; the labels Pa1–Pb5 represent abnormal ICEsat points.
Remotesensing 14 03089 g006
Table 1. Research status of geometric calibration and interferometric calibration methods.
Table 1. Research status of geometric calibration and interferometric calibration methods.
No.YearAuthorsDescription of the Method
12017R. Zhao, G. Zhang, M. Deng [19]A multimode hybrid geometric calibration method was proposed for GF-3 images, and the geometric positioning accuracy was found to be better than 3 m.
22018LÜ Guannan, Tang Xinming, Ai Bo [15]A sparse-control multiplatform geometric calibration method was proposed for Ter-raSAR-X/TanDEM-X and GF-3 images. The geometric positioning accuracy of the Ter-raSAR-X/TanDEM-X images found using this method was better than 3 m, while that of the GF-3 images was better than 7.5 m following calibration.
32020DING Liujian, TAO Qiuxiang, LI Tao [20]A wide-area joint geometric calibration method was proposed for GF-3 SAR images, and the positioning error of the calibrated SAR images was better than 9 m.
41999Zink M, Geudtner D [21]An airborne interferometric calibration method based on a sensitivity equation was deduced using flat-ground geometric relationships. The system parameters that were found to affect the interferometric elevation measurements were summarized, and the interferometric error sensitivity was analyzed.
52001J. J. Mallorqui, M. Bara, and A. Broquetas [22]An airborne interferometric calibration method was deduced based on the interferometric sensitivity equation under squint angle conditions. The simulation data confirmed that the condition numbers of different sensitivity matrices affect the interferometric calibration results differently; however, the identified corner reflector placement problem has not yet been explored.
62003Wang Yanping [23]A three-dimensional reconstruction (TDR) model was introduced to derive the interferometric sensitivity equation, and a CR placement algorithm was proposed to minimize the condition numbers in sensitivity matrices.
72010Jin Guowang, Zhang Wei, Xiang Maosheng [24]An airborne interferometric calibration model was proposed to accommodate the phase offset, baseline length and baseline angle, but the plane information of ground calibration points (GCPs) was not considered in this model.
82010Zhang Wei, Xiang Maosheng, Wu Yirong [25]An airborne dual-antenna SAR interferometric calibration method based on the three-dimensional information of control points was proposed, and this method could verify the plane and elevation accuracies of a DEM following calibration; however, the correction of each interferometric parameter was not solved independently in this method.
92011HU Jiwei, Hong Jun [26]A distributed satellite SAR interferometric altimetry model was proposed based on the TDR model, and the characteristics of the baseline vector and phase error were analyzed using simulated data.
102011Zhang Yongjun [27]An interferometric calibration optimization method for spaceborne distributed SAR systems was proposed to solve the problem of coupling the contribution of interferometric parameters to the elevation error. The simulated results showed that the method was feasible, but the existing research on each error factor has not been sufficiently deep, and studies analyzing measured supporting data have been lacking.
112016Wu Danqin [28]The interferometric calibration method derived based on the InSAR TDR model and sensitivity equation was verified using TanDEM-X data, thus proving the universality of this spaceborne SAR interferometric calibration method. However, in this technology verification, each error was not solved independently.
Table 2. Parameters of experimental images.
Table 2. Parameters of experimental images.
IDInitial Baseline in the TCN Direction (m)Baseline Rate in the TCN Direction (m/s)
A–B0, 1087.691, 419.4820, 0.596, 0.182
B–C0, −1090.332, −419.6500, −0.575, −0.070
C–D0, 1092.000, 420.8470, 0.723, 0.065
E–F0, −1117.443, −423.5860, −0.506, −0.066
Table 3. Plane precisions of the test images before the geometric calibration process.
Table 3. Plane precisions of the test images before the geometric calibration process.
SAR Image NumberNumber of GCPsRMSE
Azimuth (m)Range (m)Plane Precision (m)
A1544.78928.80353.251
B1540.55328.16049.371
C1521.52321.29330.275
D1519.57421.28828.919
E1530.53117.28235.082
F1524.07416.76329.335
Average 30.17422.26437.706
Table 4. Plane precision results of the test images after the geometric calibration process.
Table 4. Plane precision results of the test images after the geometric calibration process.
SAR Image NumberGeometric Parameter CorrectionRMSE
Δ t   ( ms ) Δ r   ( m ) Azimuth (m)Range (m)Plane Precision (m)
A−3.229−19.8433.204.005.12
B−3.227−18.1682.863.044.17
C−3.229−15.5452.975.896.59
D−2.733−21.2472.404.304.92
E−3.401−16.4303.003.244.41
F−4.632−15.3962.204.885.35
Average 2.774.225.09
Table 5. Vertical accuracies of the interferometric results before the interferometric calibration process.
Table 5. Vertical accuracies of the interferometric results before the interferometric calibration process.
Interferometric PairNumber of GCPsVertical Accuracy (m)
A–B15396.682
B–C1116.950
C–D15175.936
E–F1575.069
Average 166.159
Table 6. Vertical accuracies of the interferometric results after the interferometric calibration process.
Table 6. Vertical accuracies of the interferometric results after the interferometric calibration process.
Interferometric PairNumber of GCPsPhase Error (°)Interferometric Parameter CorrectionVertical Accuracy (m)
Δ φ Δ b c 0 Δ b n 0 Δ b c v Δ b n v
A–B15−137.337−0.1940.5580.0113−0.1206.31
B–C11−96.657−0.5340.4120.029−0.0102.26
C–D15−115.9660.816−0.592−0.1260.00414.45
E–F15−97.804−0.1270.5350.0918−0.0343.70
Average 4.18
Table 7. Accuracy evaluation results of the verification scene data.
Table 7. Accuracy evaluation results of the verification scene data.
Interferometric PairImage TypeInterferometric Parameter CorrectionVertical Accuracy (m)
Δ φ   ( ° ) Δ b c 0   ( m ) Δ b n 0   ( m ) Δ b c v   ( m / s ) Δ b n v   ( m / s )
E–FCalibration scene−97.804−0.1270.5350.0918−0.0343.70
E–FVerification scene−137.337−0.1940.558−0.01130.1204.75
Table 8. Evaluation results of the interferometric calibration accuracies derived under 4 calibration schemes.
Table 8. Evaluation results of the interferometric calibration accuracies derived under 4 calibration schemes.
Scheme Δ φ   ( ° ) Δ b c 0   ( m ) Δ b n 0   ( m ) Δ b c v   ( m / s ) Δ b n v   ( m / s ) Elevation
Accuracy (m)
(a)−108.232−0.4230.3330.0760−0.04422.513
(b)−116.711−0.4250.3340.0763−0.04452.640
(c)−96.543−0.2460.2000.775−0.5725.584
(d)−96.657−0.5340.4120.0294−0.01062.260
Table 9. DEM accuracy evaluation results (unit: m).
Table 9. DEM accuracy evaluation results (unit: m).
Interferometric PairTanDEM-X DEMICESat
A–B2.102.63
B–C2.732.78
C–D2.403.09
E–F1.442.55
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, T.; Fan, J.; Liu, Y.; Lu, R.; Hou, Y.; Lu, J. An Improved Independent Parameter Decomposition Method for Gaofen-3 Surveying and Mapping Calibration. Remote Sens. 2022, 14, 3089. https://doi.org/10.3390/rs14133089

AMA Style

Li T, Fan J, Liu Y, Lu R, Hou Y, Lu J. An Improved Independent Parameter Decomposition Method for Gaofen-3 Surveying and Mapping Calibration. Remote Sensing. 2022; 14(13):3089. https://doi.org/10.3390/rs14133089

Chicago/Turabian Style

Li, Tao, Jun Fan, Yanyang Liu, Ruifeng Lu, Yusheng Hou, and Jing Lu. 2022. "An Improved Independent Parameter Decomposition Method for Gaofen-3 Surveying and Mapping Calibration" Remote Sensing 14, no. 13: 3089. https://doi.org/10.3390/rs14133089

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