Next Article in Journal
Research on Target Deviation Measurement of Projectile Based on Shadow Imaging Method in Laser Screen Velocity Measuring System
Previous Article in Journal
Improving Positioning Accuracy via Map Matching Algorithm for Visual–Inertial Odometer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Feasibility of Replacing the Range Doppler Equation of Spaceborne Synthetic Aperture Radar Considering Atmospheric Propagation Delay with a Rational Polynomial Coefficient Model

1
School of surveying and Mapping Science and Technology, Xi’an University of Science and Technology, Xi’an 710054, China
2
State Key Laboratory of Information Engineering in Surveying Mapping and Remote Sensing, Wuhan University, Wuhan 430079, China
3
School of Geomatics, Liaoning Technical University, Fuxin 123000, China
4
China Satellite Navigation Office, Beijing 100034, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(2), 553; https://doi.org/10.3390/s20020553
Submission received: 20 November 2019 / Revised: 8 January 2020 / Accepted: 17 January 2020 / Published: 19 January 2020
(This article belongs to the Section Remote Sensors)

Abstract

:
Usually, the rational polynomial coefficient (RPC) model of spaceborne synthetic aperture radar (SAR) is fitted by the original range Doppler (RD) model. However, the radar signal is affected by two-way atmospheric delay, which causes measurement error in the slant range term of the RD model. In this paper, two atmospheric delay correction methods are proposed for use in terrain-independent RPC fitting: single-scene SAR imaging with a unique atmospheric delay correction parameter (plan 1) and single-scene SAR imaging with spatially varying atmospheric delay correction parameters (plan 2). The feasibility of the two methods was verified by conducting fitting experiments and geometric positioning accuracy verification of the RPC model. The experiments for the GF-3 satellite were performed by using global meteorological data, a global digital elevation model, and ground control data from several regions in China. The experimental results show that it is feasible to use plan 1 or plan 2 to correct the atmospheric delay error, no matter whether in plain, mountainous, or plateau areas. Moreover, the geometric positioning accuracy of the RPC model after correcting the atmospheric delay was improved to better than 3 m. This is of great significance for the efficient and high-precision geometric processing of spaceborne SAR images.

1. Introduction

Since 1999, for the IKONOS satellite, considering technical confidentiality and other factors, American space imaging companies have started to provide rational polynomial coefficient (RPC) models instead of rigorous geometry models to end users as basic imaging products [1]. The RPC model is a universal geometric model of remote sensing satellite sensors, which is usually fitted by a range Doppler (RD) model in spaceborne synthetic aperture radar (SAR) imaging systems. The difference between the RPC and RD models is that the former is a purely mathematical model, is independent of the sensor, and involves simple calculations. Furthermore, with the emergence of various imaging sensors, it is difficult for end users to add new sensor models to existing software systems to process new sensor data, while the RPC model solves this difficulty well. Moreover, from the perspective of the multi-source sensor data application, the RPC model also provides a unified geometric model for the joint adjustment of multi-source, high-resolution remote sensing images. In view of these characteristics and factors, the RPC model has been widely used in photogrammetric processing of remote sensing satellite images.
Many scholars have conducted extensive research on the calculations involved in and application of RPC models. Firstly, the direct and iterative least-squares results of the RPC parameters were derived, and terrain-dependent and terrain-independent solution methods were proposed in a photogrammetric system [1,2]. Subsequently, the terrain-independent RPC solution algorithm was studied in detail, which is based on a global digital elevation model (DEM) to interpolate the maximum and minimum elevation of the study area and does not require initial values. It was noted that the third-order RPC model with unequal denominators achieved the highest replacement accuracy [3].
The replacement accuracy of the third-order RPC model was verified for different remote sensing satellite images, whether it is a high-resolution SAR satellite, such as TerraSAR-X, COSMO-SkyMed, and Radarsat-2; or medium-resolution and low-resolution SAR satellites, such as ERS, ALOS, JERS, and ASAR, it achieved accuracies of better than 0.015 pixels [3,4,5]. In addition, many scholars have used control data and compensation models (such as translation, affine transformation, and polynomial models) to improve the geometric positioning accuracies of RPC models. Consequently, the geometric positioning accuracy of the GF-3 satellite is better than 2 pixels, that of the YaoGan-18 satellite is better than 25 m, and Radarsat can achieve a positioning accuracy of 36.735 m [6,7,8,9].
At present, although the RPC model performs well in terms of substitution and positioning accuracy, the slant range measurement error caused by the two-way atmospheric propagation delay of the SAR signal is still not taken into account in the RPC model, and research has shown that the slant range measurement error will bring the positioning error in the range direction [10]. Meanwhile, some researchers pointed out that the range accuracy of SAR using modern satellites has been verified in centimeter-range after correcting geodynamic, systematic errors, and atmospheric effects [11,12]. Moreover, with the continuous launch of high-resolution SAR satellites, the accuracy and efficiency requirements of SAR image geometry processing have become more stringent. Therefore, in this study, we developed an RPC fitting method that corrects the atmospheric propagation delay, analyzed the feasibility of replacing the RD equation considering atmospheric delay correction with the RPC model, and verified the positioning accuracy of the RPC model of spaceborne synthetic aperture radar (SAR) considering atmospheric propagation delay.
GF-3 is the first C-band multi-polarized SAR satellite with 1 m resolution in China. It has 12 imaging modes, including SpotLight, StripMap, and ScanSAR, along with a wave imaging mode, and can provide data support services for research on water disaster monitoring and assessment, climate change, oceans, agriculture, forestry, and earthquakes. In this study, a series of experiments was conducted using GF-3 SAR satellite data and manually measured real ground control point (GCP) data from several regions in China.
Based on the RD model, this paper describes a method of fitting the RPC model of a spaceborne SAR satellite considering the atmospheric propagation delay and presents an evaluation of the feasibility of this replacement. Section 2 describes three basic models: the atmospheric propagation delay correction model, RD model, and RPC model. The RPC model fitting and accuracy verification methods that corrected the atmospheric propagation delay are presented in Section 3. The atmospheric propagation delay analysis results and RPC model fitting and positioning accuracy verification are discussed in Section 4. Finally, Section 5 presents the conclusions.

2. Models

2.1. Atmospheric Propagation Path Delay Equation

The atmospheric propagation path delay of a radar signal Δ L d e l a y depends on a complex atmospheric spatial distribution pattern, generally at a specific height h = h t , and is represented as the integral of the sum of the dry, wet, and ionosphere components of the refractivity N, from h t to the top of the atmosphere h [13,14]:
N = ( k 1 P T ) d r y + ( k 2 e T + k 3 e T 2 + k 4 W c l o u d ) w e t + ( k 5 n e f 2 ) i o n o s p h e r e = N d r y + N w e t + N i o n o s p h e r e
Δ L d e l a y = 1 0 - 6 c o s θ h t h ( N d r y + N w e t + N i o n o s p h e r e ) d h
where P is the total pressure (mbar), T is the temperature (K), e is the partial pressure of water vapor (mbar), and W c l o u d is the cloud water content (g/m3). These meteorological data can be downloaded from the U.S. National Center for Environmental Prediction (NCEP) website (https://rda.ucar.edu/datasets/ds083.2/). In addition, n e is the electronic density of the ionosphere, which can be obtained by bilinear interpolation in time and space based on the global ionospheric map provided by the European Center for Orbit Determination (CODE) (ftp://cddis.gsfc.nasa.gov/gnss/products/ionex/). Here, f and θ are the frequency and incidence angle of the radar, respectively, which can be obtained from the auxiliary file of the SAR satellite image. The constant coefficients have the following values: k1 = 77.6 K/mbar, k2 = 26.0 K/mbar, k3 = 3.753105 K2/mbar, k4 = 1.45 m3/g, and k5 = 24.033107 m3/s2 [15,16,17].

2.2. Spaceborne SAR RD Equation

The RD equation establishes the correspondence between the image and object point coordinates from the perspective of the spaceborne SAR imaging geometry. It consists of an earth model equation, SAR doppler equation, and SAR range equation [18]:
{ X t 2 + Y t 2 A 2 + Z t 2 B 2 = 1 f D = - 2 λ R ( R s - R t ) ( V s - V t ) R = ( R s - R t ) ( R s - R t ) ,
where R s and R t = [ X t Y t Z t ] are the position vectors of the SAR satellite and target, respectively; V s and V t are the satellite and target velocity vectors, respectively; R is the slant range; λ is the radar wavelength; A = 6378.137 km is the mean equatorial radius; and B = ( 1 - 1 / f ) A is the polar radius with a flattening factor f = 298.255.
Considering the slant range measurement error caused by the atmospheric propagation delay, the slant range can be expressed as
R = ( R s - R t ) ( R s - R t ) + Δ L d e l a y .

2.3. RPC Model

The RPC model uses the ratio polynomial to represent the correspondence between the ground coordinate D ( L a t i t u d e , L o n g i t u d e , H e i g h t ) (the polar stereographic coordinate system is more applicable at 70–75 degrees latitude north or south) and the pixel coordinates d ( l i n e , s a m p l e ) [19]:
Y = N l ( P , L , H ) D l ( P , L , H ) X = N s ( P , L , H ) D s ( P , L , H ) } ,
where N l ( P , L , H ) , N s ( P , L , H ) , D l ( P , L , H ) , and D s ( P , L , H ) are all third-order polynomials. These functions have the form
F ( P , L , H ) = a 1 + a 2 L + a 3 P + a 4 H + a 5 L P + a 6 L H + a 7 P H + a 8 L 2 + a 9 P 2 + a 10 H 2 + a 11 P L H + a 12 L 3 + a 13 L P 2 + a 14 L H 2 + a 15 L 2 P + a 16 P 3 + a 17 P H 2 + a 18 L 2 H + a 19 P 2 H + a 20 H 3
where a i is the RPC model coefficient, i = 1 , 2 , , 20 ; Generally, the first parameter coefficient of the denominator term is set to 1, so there are 78 parameters for RPC model.
Here, ( P , L , H ) and ( X , Y ) are the normalized ground and image coordinates, respectively. These coordinates can be expressed as
P = D l a t - D l a t _ o f f D l a t _ s c a l e L = D l o n - D l o n _ o f f D l o n _ s c a l e H = D h e i - D h e i _ o f f D h e i _ s c a l e }
X = s - s o f f s s c a l e Y = l - l o f f l s c a l e } ,
where D i ( i = l a t _ o f f , l a t _ s c a l e , l o n _ o f f , l o n _ s c a l e , h e i _ o f f , h e i _ s c a l e ) is the normalized parameter of the ground coordinates and s j , l j ( j = o f f , s c a l e ) are the normalized parameters of the image coordinates.
The RPC fitting results obtained based on a variety of spaceborne SAR sensors have shown that RPC models with third-order and unequal denominators have the highest replacement accuracies [7]. The following experiments were based on this model.

3. RPC Model Fitting and Accuracy Verification Methods

The RPC model is essentially a mathematical model. In this study, we used a terrain-independent solution algorithm to calculate the RPC parameters. The core of the algorithm involves establishing a virtual control grid based on the SAR image, DEM data of the coverage area, and forward transformation of the RD model. Then, the RPC parameters are fitted according to the least squares fitting method and the correspondence between the ground coordinates and image pixel coordinates of the virtual control points. However, every certain virtual grid point in space is affected by atmospheric propagation delay, which causes slant range measurement error. Therefore, we used two atmospheric delay correction schemes to correct the slant range values at all virtual grid points, then re-fitted the RPC model and evaluated the fitting and positioning accuracies.
The steps of the algorithm are listed in Figure 1.
The steps of the algorithm can be described as follows.
(1)
Firstly, the input data sets, including the SAR image auxiliary file and DEM of the study area, were prepared;
(2)
Then, a virtual plane grid was built according to the image size, the elevation range of the coverage area was interpolated based on the DEM data and RD model, and layering processing was performed in the elevation direction to obtain a virtual space control grid. The number of layers was required to be greater than three to prevent the design matrix from becoming ill-conditioned [1]. Next, a virtual space check grid was constructed by interpolating between the centers of four adjacent virtual control points and the two adjacent elevation layers;
(3)
The atmospheric propagation delay values were corrected at all virtual grid points using two plans: one in which the atmospheric delay correction value at the point at the center of the scene and the average elevation of the coverage area were employed to correct the atmospheric propagation delay at all virtual points, and one in which the atmospheric delay correction value at each virtual grid point was used to correct its own atmospheric delay;
(4)
The RPC model was fitted separately using the original RD model, RD model modified using atmospheric delay correction plan 1, and RD model modified using atmospheric delay correction plan 2, and then the fitting accuracy was evaluated. A detailed flowchart of the RPC model fitting accuracy evaluation is shown in Figure 2.
(5)
The positioning accuracies of the abovementioned RPC models were evaluated based on the measured GCPs of the SAR image coverage area. To ensure the reliability of the verification accuracy, the root mean square error (RMSE) of multiple control points in one image was calculated using the equation in step 5, ∆X and ∆Y as an evaluation index of the geometric positioning accuracy of the RPC model.

4. Results and Discussion

4.1. Experimental Data

In this study, GF-3 images of four imaging modes, namely Spot-Light (1 m resolution), Ultra-Fine-Strip (3 m resolution), Fine-Strip-I (5 m resolution), and Fine-Strip-II (10 m resolution), were used as experimental data. To verify the positioning accuracy of the RPC model considering the atmospheric propagation delay, GF-3 satellite images and GCP data from Anping, Zhanjiang, Xianning, Zhangye, and Inner Mongolia in China were used. The GCPs were derived from global navigation satellite system (GNSS) receivers, which were selected in locations with obvious features, such as the intersections of pathways. Centimeter-level positioning accuracy could be achieved using real-time kinematic static observation technology. The details of the GF-3 satellite images are provided in Table 1, and the distribution of the experimental data and GCPs is depicted in Figure 3.

4.2. Experimental Results and Analysis

To evaluate the feasibility of replacing the range Doppler equation of spaceborne SAR by considering the atmospheric delay with the RPC model, the following three experiments were conducted. (1) For the terrain-independent calculation algorithm of the RPC model of spaceborne SAR, the effects of the virtual grid pattern, number of elevation layers, and orbit fitting order on the fitting accuracy of the RPC model were analyzed. (2) Based on the appropriate grid pattern, number of elevation layers, and orbital fitting order, the fitting accuracy of the RPC model was assessed using the two proposed atmospheric delay correction plans. (3) Based on the fitting results of the RPC models, with and without considering the atmospheric delay of the radar signal, combined with the GCP data in the image coverage area, the positioning accuracy of the RPC model considering the atmospheric delay was verified.

4.2.1. Impact of RPC Model Fitting Parameter

(1) Grid Size
Assuming there to be five elevation layers, six different grid styles were set (style 1: 4000 pixels × 4000 pixels; style 2: 2000 pixels × 2000 pixels; style 3: 1000 pixels × 1000 pixels; style 4: 500 pixels × 500 pixels; style 5: 200 pixels × 200 pixels; style 6: 100 pixels × 100 pixels), here “N pixels × N pixels” refers to a grid spacing of N pixels, or a grid size. We take N pixels as the sampling interval in the range direction and azimuth direction, and take the sampling points as the virtual control points. Based on the fine-strip product data of the GF-3 satellite, DEM data with 30 m resolution, and the RPC fitting method described in Section 3, the influence of the grid size on the fitting accuracy of the RPC model was analyzed. The experimental results are shown in Figure 4a.
(2) Number of Elevation Layers
Using the appropriate grid pattern identified in the previous step, the number of different elevation layers (two, three, four, five, or six) was set to evaluate the influence of the number of elevation layers on the fitting accuracy of the RPC model. The experimental results are presented in Figure 4b.
(3) Orbit Fitting Order
With the appropriate grid pattern and number of elevation layers, different orders of satellite orbital fitting (third, fourth, or fifth order) were set to analyze the influence of the orbit fitting order on the fitting accuracy of the RPC model. The experimental results are depicted in Figure 4c.
Figure 4 shows that as the control grid becomes smaller and the number of elevation layers increases, the plane RMSE of the checkpoint gradually decreases. Eventually, the accuracy of the RPC model becomes close to that of the rigorous imaging geometry model. As the orbital parameter fitting order increases, the plane RMSE of the checkpoint gradually increases. To balance the fitting accuracy and computational efficiency of the RPC model, we set the grid size to 500 pixels × 500 pixels, the number of elevation layers to five, and the orbit fitting order to third-order for the subsequent study.

4.2.2. RPC Model Fitting Accuracy Evaluation

Based on the atmospheric propagation delay correction model described in this paper combined with the meteorological data stored by the NCEP every 6 h, the global ionospheric map provided by CODE every 1 h, and the average elevation and size of the image coverage area provided by the GF-3 image auxiliary file, the atmospheric propagation delay correction of the radar signal was calculated to analyze the influence of atmospheric propagation delay correction on the RPC fitting accuracy.
In this study, atmospheric propagation delay correction of the radar signal was performed using two methods: one with a GF-3 image of a single scene and unique atmospheric delay correction parameters (plan 1) and one with an image of a single scene and spatially varying atmospheric delay correction parameters (plan 2). The two plans can be described as follows. In plan 1, the atmospheric delay correction at the point at the center of the scene and the average elevation of the coverage area are calculated to correct the atmospheric propagation delay at all virtual grid points. In plan 2, the atmospheric propagation delay correction value is calculated at each virtual grid point to correct its own atmospheric delay. The delay results are presented in Table 2.
In Table 2, the calculation results for plan 1 show that the maximum slant range correction value at the center point among the different images is −3.795 m, and the delay error in the ground range direction was calculated based on the incident angle of the image as approximately −5 m, which indicates that the atmospheric delay error has a significant influence on the positioning accuracy. The calculation results for plan 2 show that the atmospheric delay correction values of the virtual grid points with different spatial distributions on a scene image are different, and the slant range correction values of the virtual grid points in all of the image scenes are between 1.521 m and −3.919 m. Thus, the atmospheric delay correction values differ between imaging angles and imaging regions. It was, therefore, concluded that the atmospheric propagation delay is mainly affected by the radar incidence angle and the topography of the study area.
To analyze the main influencing factors of the atmospheric propagation delay in the SAR image scenes further, we visualized the atmospheric propagation delay correction values of the virtual grid points calculated using plan 2, assuming this to be 0 pixels in the azimuth direction and sampling every 500 pixels in the range direction, with five elevation layers. The trend of the atmospheric delay correction value in the range direction is shown in Figure 5. To analyze the influence of the elevation on the atmospheric propagation delay, we increased the number of elevation layers of the virtual control grid to 300 separately in the plains and mountains area (ZJ-5851) and the plateau area (ZY-7668). The maximum variation ∆φ of the corresponding incident angle at the same sampling point in each elevation layer in the experimental region is shown in Figure 6, where the orbital height of the satellite is known to be approximately 755 km. Here, ∆φ was approximated using Equation (9), and was found to be less than 0.03° in the ZJ-5851 image coverage area and less than 0.15° in the ZY-7668 image coverage area. Thus, the variation of the incident angle in the elevation direction in the study area is negligible, and the same sampling point can be selected in each elevation layer to analyze the atmospheric propagation delay law as a function of elevation, as shown in Figure 7.
Δ φ < H m a x H s × 180 ° π ,
where Hs is the orbital height of the GF-3 satellite; Hmax and Hmin are the maximum and minimum elevations of the image coverage area, respectively.
Figure 5 shows that the atmospheric delay and incident angle are approximately linearly related at the same elevation level. The absolute value of the atmospheric propagation delay correction gradually increases as the range increases.
Some researchers have demonstrated that the elevation–atmospheric delay relationship is usually a linear or exponential function [20,21]. Figure 7a shows that in the plains and mountains area, the elevation–atmospheric delay relationship is approximately linear. Meanwhile, Figure 7b reveals that in the plateau region, the elevation–atmospheric delay relationship is also approximately linear.
According to the RPC model fitting parameter impact analysis results in Section 4.2.1, we established a control grid by setting the grid size to 500 pixels × 500 pixels and the number of elevation layers to five. The checkpoint was determined based on the centers of the four adjacent control points and two adjacent elevation layers, and the fitting accuracy of the RPC model was evaluated after correcting the atmospheric delay grid by grid. The results are presented in Table 3.
In general, as the plane error of the checkpoint does not exceed 5% of the pixels, we believe that the RPC model can be used instead of the RD model for photogrammetric processing of spaceborne SAR satellites. Table 3 shows that in the same scene, both the uncorrected atmospheric delay and atmospheric delay correction with plan 1 yield exactly the same RPC model fitting accuracy, and that the maximum plane RMSE of the checkpoint is 0.00356 pixels. For the atmospheric delay correction with plan 2, the fitting accuracy of the RPC model is very approximate, so both approaches achieve precision better than 1% of a pixel. This indicates that the RPC model still provides high accuracy for the RD model after correcting the atmospheric delay using two plans, in plain, mountainous, and plateau areas.
The reason for obtaining these experimental results can be confirmed by research that has shown that the third-order RPC model can better fit a function whose highest power is no more than 5 [22]. For a single-scene SAR image, plan 1 uses a unique atmospheric delay correction parameter to correct the atmospheric propagation delay at all virtual grid points, which is equivalent to adding the atmospheric propagation delay as a constant function to the slant range (R) term of the RD model. Therefore, the fitting accuracy of the RPC model is unchanged. Meanwhile, plan 2 corrects the atmospheric propagation delay at all virtual grid points with the correction parameters of the corresponding virtual grid points, which is equivalent to adding the atmospheric delay correction model to the R term with the elevation as the dependent variable. Figure 7a shows that the elevation–atmospheric delay function is approximately linear in the plains and mountains area, so it can be correctly fitted with a third-order RPC model. Meanwhile, Figure 7b shows that in the plateau area, the elevation–atmospheric delay function also shows an approximate linear relationship, which indicates that RPC model can not only fit complex terrain, but also fit the terrain-related atmospheric propagation delay.

4.2.3. RPC Model Positioning Accuracy Verification

From the perspective of RPC model application, we carried out RPC model positioning experiments to verify the effectiveness of the RPC model considering atmospheric propagation delay. Based on the experimental data of the GF-3 satellite and GCPs, geometric positioning verification experiments were performed for the four RPC models, and the results are shown in Table 4.
(1)
Firstly, for the RPC model fitted using the original RD model without atmospheric delay correction, positioning accuracy verification was conducted;
(2)
After correcting the systematic error in the slant range direction of the GF-3 SAR satellite, geometric positioning accuracy verification of the RPC model was performed. This systematic error was obtained by geometric calibration of the spaceborne SAR [23];
(3)
After correcting the systematic error in the slant range direction of the GF-3 SAR satellite and correcting the atmospheric propagation delay error using plan 1, geometric positioning accuracy verification of the RPC model was performed;
(4)
After correcting the systematic error in the slant range direction of the GF-3 SAR satellite and correcting the atmospheric propagation delay error using plan 2, geometric positioning accuracy verification of the RPC model was performed.
Table 4 shows that the positioning accuracy of the RPC model is between 21.241 m and 26.004 m without correction of system error and atmospheric delay error. After correcting the system error, the geometric positioning error of the RPC model was between 3.064 m and 5.529 m. After correcting the atmospheric delay error using plans 1 and 2, respectively, the maximum geometric positioning error of the RPC model was found to be 2.948 m and 2.957 m, and the positioning accuracy of a single scene image could be increased by 2.7 m, such as for ZJ-5492. Thus, it is effective to improve the geometric positioning accuracy of the RPC model by correcting the atmospheric propagation delay in the process of fitting the RD model to the RPC model. Furthermore, the positioning accuracy of the four RPC models in the azimuth direction is almost unchanged, which indicates that the slant range error caused by the atmospheric propagation delay mainly affects the positioning accuracy in the range direction.
Figure 8 shows that the positioning error of all GCPs in the SAR image of a single scene is approximately 2.5 m after correcting the atmospheric propagation delay using plans 1 and 2 separately, and this remains relatively stable, which indicates that the geometric positioning results of the RPC model are reliable.
It is worth noting that although both plans are feasible, it is still difficult to produce a RPC model considering atmospheric delay directly and in real time in the SAR system, because sufficiently precise atmospheric models typically have latency. We may be able to choose a general model in which the atmospheric delay is a function of the SAR incidence angle and the elevation of a given point, in order to replace the precise atmospheric model [24]. This will be our next research direction.

5. Conclusions

The atmospheric delays of SAR satellite radar signals can reach several meters, which is non-negligible in high-precision geometric positioning. Thus, two methods of atmospheric delay correction during RPC fitting of an RD model were designed and tested in this study. Based on the GF-3 SAR satellite data, a global digital elevation model, global meteorological data, and ground control data from several regions in China, the feasibility of the two methods was verified by conducting RPC model fitting and positioning experiments. We found that the RPC model exhibited high substitution and positioning accuracies after using plan 1 or plan 2 to correct the atmospheric propagation delay, in both the plains and mountains area where the relative height difference is less than 1000 m, and in the plateau area where the relative height difference is more than 1000 m. Therefore, it is feasible to use plan 1 or plan 2 to correct the atmospheric delay error under any terrain conditions. Moreover, the geometric positioning accuracy of the RPC model after correcting the atmospheric delay was improved to better than 3 m.

Author Contributions

G.Z., Y.H., R.Z., S.H. and P.J. conceived and designed the experiments. G.Z., S.H., and R.Z. performed the experiments. S.H., R.Z. and P.Z. analyzed the data. S.H. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Key Research and Development Program of Ministry of Science and Technology, grant 2016YFB0500801; the National Natural Science Foundation of China, grants 91538106, 41501503, 41601490, and 41501383; and the Open Research Fund of State Key Laboratory of Information Engineering in Surveying, Mapping, and Remote Sensing, grant 15E02. The authors also thank the anonymous reviewers for their constructive comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tao, C.V.; Hu, Y. A comprehensive study of the rational function model for photogrammetric processing. Photogramm. Eng. Remote Sens. 2001, 67, 1347–1357. [Google Scholar]
  2. Toutin, T. Review article: Geometric processing of remote sensing images: Models, algorithms and methods. Int. J. Remote Sens. 2004, 25, 1893–1924. [Google Scholar] [CrossRef]
  3. Zhang, G.; Fei, W.B.; Li, Z.; Zhu, X.; Li, D. Evaluation of the RPC model for spaceborne SAR imagery. Photogramm. Eng. Remote Sens. 2010, 76, 727–733. [Google Scholar] [CrossRef]
  4. Eftekhari, A.; Saadatseresht, M.; Motagh, M. A study on rational function model generation for Terrasar-X imagery. Sensors. 2013, 13, 12030–12043. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Zhang, G.; Zhu, X.Y. A study of the RPC model of Terrasar-X and Cosmo-Skymed SAR imagery. Remote Sens. Spat. Inf. Sci. 2008, 37, 321–324. [Google Scholar]
  6. Wang, T.Y.; Zhang, G.; Yu, L.; Zhao, R.; Deng, M.; Xu, K. Multi-mode gf-3 satellite image geometric accuracy verification using the RPC model. Sensors 2017, 17, 2005. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Jiao, N.G.; Wang, F.; You, H.J.; Qiu, X.L.; Yang, M.D. Geo-positioning accuracy improvement of multi- mode gf-3 satellite SAR imagery based on error sources analysis. Sensors 2018, 18, 2333. [Google Scholar] [CrossRef] [Green Version]
  8. Zhao, R.S.; Jiang, Y.H.; Zhang, G.; Deng, M.J.; Yang, F. Geometric accuracy evaluation of YG-18 satellite imagery based on RFM. Photogramm. Rec. 2017, 32, 33–47. [Google Scholar] [CrossRef] [Green Version]
  9. Huang, Z.W.; He, S. A method for improving positioning accuracy of SAR imagery based on RFM. In Proceedings of the 4th International Conference on Digital Manufacturing & Automation, Qindao, China, 29–30 June 2013; pp. 43–46. [Google Scholar]
  10. Liu, X.F.; Liu, J.Y.; Hong, W. The analysis of the precision in spaceborne SAR image location. J. Remote Sens. 2006, 10, 76–81. [Google Scholar]
  11. Eineder, M.; Minet, C.; Steigenberger, P.; Cong, X.; Fritz, T. Imaging Geodesy—Toward Centimeter-Level Ranging Accuracy with TerraSAR-X. IEEE Trans. Geosci. Remote Sens. 2011, 49, 661–671. [Google Scholar] [CrossRef]
  12. Cong, X.; Balss, U.; Eineder, M.; Fritz, T. 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]
  13. Davis, J.L.; Herring, T.A.; Shapiro, I.I.; Rogers, A.E.E.; Elgered, G. Geodesy by radio interferometry: Effects of atmospheric modeling errors on estimates of baseline length. Radio Sci. 2016, 20, 1593–1607. [Google Scholar] [CrossRef]
  14. Puyssegur, B.; Michel, R.; Avouac, J.P. Tropospheric phase delay in interferometric synthetic aperture radar estimated from meteorological model and multispectral imagery. J. Geophys. Res. Solid Earth. 2007, 112, 12. [Google Scholar] [CrossRef] [Green Version]
  15. Smith, E.K.; Weintraub, S. The constants in the equation for atmospheric refractive index at radio frequencies. Proc. IRE 1953, 41, 1035–1037. [Google Scholar] [CrossRef] [Green Version]
  16. Solheim, F.S.; Vivekanandan, J.; Ware, R.H.; Rocken, C. Propagation delays induced in GPS signals by dry air, water vapor, hydrometeors, and other particulates. J. Geophys. Res. 1999, 104, 9663–9670. [Google Scholar] [CrossRef]
  17. Skone, S.; Cannon, M.E. Ionospheric effects on differential GPS applications during auroral substorm activity. ISPRS J. Photogramm. Remote Sens. 1999, 54, 279–288. [Google Scholar] [CrossRef]
  18. Curlander, J.C. Location of spaceborne SAR imagery. IEEE Trans. Geosci. Remote Sens. 1982, 3, 359–364. [Google Scholar] [CrossRef]
  19. Grodecki, J.; Dial, G. Block adjustment of high-resolution satellite images described by rational polynomials. Photogramm. Eng. Remote Sens. 2003, 69, 59–68. [Google Scholar] [CrossRef]
  20. Cavalie, O.; Doin, M.P.; Lasserre, C.; Briole, P. Ground motion measurement in the Lake Mead area, Nevada, by differential synthetic aperture radar interferometry time series analysis: Probing the lithosphere rheological structure. J. Geophys. Res. Solid Earth 2007, 112, B03403. [Google Scholar] [CrossRef] [Green Version]
  21. Onn, F.; Zebker, H.A. Correction for interferometric synthetic aperture radar atmospheric phase artifacts using time series of zenith wet delay observations from a GPS network. J. Geophys. Res. Solid Earth 2006, 111. [Google Scholar] [CrossRef]
  22. Huang, W.C. Research on Compensation for System Errors of Basic Satellite Products. Ph.D. Thesis, Wuhan University, Wuhan, China, October 2016. [Google Scholar]
  23. Zhao, R.S.; Zhang, G.; Deng, M.J.; Xu, K.; Guo, F.C. Geometric calibration and accuracy verification of the GF-3 satellite. Sensors 2017, 17, 1977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Jehle, M.; Perler, D.; Small, D.; Schubert, A.; Meier, E. Estimation of Atmospheric Path Delays in TerraSAR-X Data using Models vs. Measurements. J. Sens. 2008, 8, 8479–8491. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Step-by-step routine for the rational polynomial coefficient (RPC) model fitting and accuracy verification algorithm. Note: SAR = synthetic aperture radar; DEM = digital elevation model; RD = range Doppler; NCEP = National Center for Environmental Prediction; CODE = European Center for Orbit Determination; GCP = ground control point.
Figure 1. Step-by-step routine for the rational polynomial coefficient (RPC) model fitting and accuracy verification algorithm. Note: SAR = synthetic aperture radar; DEM = digital elevation model; RD = range Doppler; NCEP = National Center for Environmental Prediction; CODE = European Center for Orbit Determination; GCP = ground control point.
Sensors 20 00553 g001
Figure 2. Flowchart of RPC fitting accuracy evaluation.
Figure 2. Flowchart of RPC fitting accuracy evaluation.
Sensors 20 00553 g002
Figure 3. Distribution of the experimental data and GCPs.
Figure 3. Distribution of the experimental data and GCPs.
Sensors 20 00553 g003
Figure 4. Impacts of RPC model fitting parameters: (a) grid size; (b) number of elevation layers; (c) orbit fitting order.
Figure 4. Impacts of RPC model fitting parameters: (a) grid size; (b) number of elevation layers; (c) orbit fitting order.
Sensors 20 00553 g004
Figure 5. Trend of atmospheric delay correction value in the range direction: (a) plains and mountains area (ZJ-5851); (b) plateau area (ZY-7668).
Figure 5. Trend of atmospheric delay correction value in the range direction: (a) plains and mountains area (ZJ-5851); (b) plateau area (ZY-7668).
Sensors 20 00553 g005
Figure 6. Variation of incident angle in the elevation direction in the study area.
Figure 6. Variation of incident angle in the elevation direction in the study area.
Sensors 20 00553 g006
Figure 7. Atmospheric delay correction as a function of elevation: (a) plains and mountains area (ZJ-5851); (b) plateau area (ZY-7668).
Figure 7. Atmospheric delay correction as a function of elevation: (a) plains and mountains area (ZJ-5851); (b) plateau area (ZY-7668).
Sensors 20 00553 g007
Figure 8. Positioning error distributions of all GCPs in one SAR image: (a) atmospheric delay corrected with plan 1 for image AP-6616; (b) atmospheric delay corrected with plan 2 for image AP-6616; (c) atmospheric delay corrected with plan 1 for image NM-3601; (d) atmospheric delay corrected with plan 2 for image NM-3601.
Figure 8. Positioning error distributions of all GCPs in one SAR image: (a) atmospheric delay corrected with plan 1 for image AP-6616; (b) atmospheric delay corrected with plan 2 for image AP-6616; (c) atmospheric delay corrected with plan 1 for image NM-3601; (d) atmospheric delay corrected with plan 2 for image NM-3601.
Sensors 20 00553 g008
Table 1. GF-3 data information of the study area.
Table 1. GF-3 data information of the study area.
Imaging ModeID of ImageImaging RegionImaging DateRrRaNumber of GCPs
SL (1 m resolution)CD-3076Chengdou8 January 20170.56210.3126
DF-6954Dengfeng7 August 20180.56210.3125
LZ-4905Lanzhou6 August 20180.56210.3443
UFS (3 m resolution)BC-1094Beichen18 September 20161.12421.7308
ZJ-8807Zhanjiang12 August 20171.12421.7327
SH-4753Shanghai27 January 20171.12421.7292
FSI (5 m resolution)AP-6616Anping19 February 20172.24842.811012
AP-4130Anping5 May 20171.12422.61836
ZJ-5851Zhanjiang25 January 20172.24842.8109
XN-3579Xianning11 January 20172.24842.81396
NM-3601Neimeng11 January 20172.24842.809813
NM-2661Neimeng28 July 20171.12422.591510
ZY-7668Zhangye18 May 20192.24842.82525
FSII (10 m resolution)CD-7465Chengdou3 October 20172.24844.7566
ZJ-5492Zhanjiang6 October 20172.24844.77108
ZY-0316Zhangye22 January 20192.24844.75767
Table 2. Atmospheric delay correction value calculation results.
Table 2. Atmospheric delay correction value calculation results.
TopographyImage IDIncidence Angle (°)DEM (m)ΔRdelay (m)
Near-RangeFar-RangeMinimumMaximumMeanPlan 1Plan 2
Plains and mountainsSH-475345.69947.142−30360.00−3.611−3.678–−3.544
AP-413026.18229.657−286511.66−2.825−2.825–−2.712
AP-661644.74047.334−277223.54−3.504−3.623–−3.399
BC-109419.75922.199−298113.68−2.813−2.852–−2.768
ZJ-880732.37534.393−3016924.05−3.185−3.258–−3.071
ZJ-585146.51048.971−3535023.03−3.795−3.919–3.587
XN-357944.77747.362−21641111.59−3.393−3.771–3.281
CD-307626.16027.019426975695.23−2.522−2.613–−2.400
PlateauZJ-549231.29138.155−33127122.52−3.531−3.702–2.824
DF-695428.71729.6983121470446.73−2.865−2.943–−2.477
NM-360144.74647.33497619701052.78−3.028−3.516–−2.665
NM-266123.86527.6749522214990.42−2.561−2.617–2.136
ZY-766834.66537.920137132221498.87−2.536−2.636–2.041
ZY-031630.89437.837130047692092.58−2.28−2.569–−1.568
LZ-490541.03541.745147821501606.19−3.012−3.095–−2.780
CD-746531.29437.7183904931810.42−3.111−3.212–1.521
Table 3. Fitting accuracy of the RPC model. Note: Max. = maximum; RMSE = root mean square error.
Table 3. Fitting accuracy of the RPC model. Note: Max. = maximum; RMSE = root mean square error.
TopographyImage IDCompensation MethodControl Point Error (Pixels)Check Point Error (Pixels)
SampleLine2-DSampleLine2-D
Max.RMSEMax.RMSEMax.RMSEMax.RMSEMax.RMSEMax.RMSE
Plains and MountainsSH-4753Uncorrected/Plan 10.000410.000130.000130.000070.000410.000150.000290.000120.000150.000070.00030.00014
Plan 20.000410.000130.000130.000070.000410.000150.000290.000120.000150.000070.00030.00014
AP-4130Uncorrected/Plan 10.01220.003640.00320.00120.012450.003830.009330.003380.002190.001110.00950.00356
Plan 20.012430.003660.00320.00120.012680.003850.009440.00340.002190.001110.009610.00357
AP-6616Uncorrected/Plan 10.001720.000590.000090.000040.001720.000590.001260.000550.000120.000050.001260.00055
Plan 20.001710.000590.000090.000040.001710.00060.001270.000550.000120.000050.001270.00055
BC-1094Uncorrected/Plan 10.000740.000250.000120.000070.000750.000260.000530.000220.000160.000070.000550.00024
Plan 20.000730.000260.000120.000070.000740.000270.000530.000230.000160.000070.000550.00025
ZJ-8807Uncorrected/Plan 10.000770.000240.00120.000470.001320.000520.000590.000220.000770.000430.000840.00048
Plan 20.000780.000240.00120.000470.001320.000520.000590.000220.000770.000430.000840.00048
ZJ-5851Uncorrected/Plan 10.002760.00090.001640.000630.003180.00110.002030.000830.0010.000580.002180.00102
Plan 20.002760.00090.001640.000630.003160.00110.002030.000830.0010.000580.002180.00102
XN-3579Uncorrected/Plan 10.000840.00020.000310.000070.000890.000210.000560.000180.00020.000070.000590.0002
Plan 20.000830.000190.000310.000070.000880.000210.000540.000180.00020.000070.000570.00019
CD-3076Uncorrected/Plan 10.000150.000050.000660.000320.000660.000330.000120.000050.000930.000460.000930.00046
Plan 20.000150.000050.000660.000320.000660.000330.000120.000050.000930.000460.000930.00046
PlateauZJ-5492Uncorrected/Plan 10.008390.002070.000580.000140.008390.002070.006580.001920.000440.000130.006580.00193
Plan 20.008430.002090.000580.000140.008430.002090.00660.001940.000440.000130.00660.00194
DF-6954Uncorrected/Plan 10.001340.000460.000790.000390.001460.000610.001120.000440.001190.000430.001580.00061
Plan 20.001340.000470.000790.000390.001470.000610.001120.000440.001190.000430.001580.00062
NM-3601Uncorrected/Plan 10.000720.000180.000260.000110.000740.000210.000540.000170.000240.00010.000550.0002
Plan 20.000750.000190.000260.000110.000770.000220.000570.000170.000240.00010.000590.0002
NM-2661Uncorrected/Plan 10.008320.001790.001060.000410.008390.001840.005550.001650.000780.000380.005610.00169
Plan 20.008330.001790.001060.000410.008390.001840.005530.001650.000780.000380.005590.00169
ZY-7668Uncorrected/Plan 10.001110.000250.001670.000670.001850.000720.000730.000230.001130.000650.001340.00069
Plan 20.001230.000260.001670.000710.002030.000750.000740.000230.001130.000650.001340.00069
ZY-0316Uncorrected/Plan 10.01060.002160.001410.000330.010610.002190.007560.002010.000940.000310.007560.00203
Plan 20.01060.002150.001420.000330.010610.002170.007520.001990.000940.000310.007520.00202
LZ-4905Uncorrected/Plan 10.00180.000670.001330.000510.002240.000850.001490.000640.001580.000530.002160.00083
Plan 20.001810.000670.001330.000510.002220.000850.00150.000640.001580.000530.002160.00083
CD-7465Uncorrected/Plan 10.010.0020.00670.002410.011130.003710.008260.002610.004690.002230.008970.00344
Plan 20.010.0020.00670.002410.011120.00370.008270.002610.004690.002230.008990.00344
Table 4. Geometric positioning accuracies of the four RPC models.
Table 4. Geometric positioning accuracies of the four RPC models.
Image IDCorrection Method
UncorrectedCorrect System ErrorPlan 1Plan 2
SampleLine2-DSampleLine2-DSampleLine2-DSampleLine2-D
AP-661622.8581.38022.9004.0211.3804.4092.5881.3802.9342.5991.3802.944
AP-413021.6191.00021.6422.7821.0003.0641.8781.0002.1291.8711.0002.123
XN-357922.3971.45722.4463.5611.4574.0042.3821.4572.8012.3801.4572.792
NM-360122.4600.11822.4963.6230.1173.9412.1100.1172.4590.5100.1182.486
NM-266121.2071.18821.2412.7301.1883.0751.6031.1882.0211.6021.1882.021
ZY-766823.1981.37223.2384.3601.3734.5711.4471.3721.9941.3881.3721.952
ZJ-549225.8992.33326.0045.0132.3335.5291.6002.3322.8281.8852.3332.999
ZY-031624.9831.00125.0034.0961.0004.2171.1841.0001.5501.1661.0001.536

Share and Cite

MDPI and ACS Style

Hou, S.; Huang, Y.; Zhang, G.; Zhao, R.; Jia, P. Feasibility of Replacing the Range Doppler Equation of Spaceborne Synthetic Aperture Radar Considering Atmospheric Propagation Delay with a Rational Polynomial Coefficient Model. Sensors 2020, 20, 553. https://doi.org/10.3390/s20020553

AMA Style

Hou S, Huang Y, Zhang G, Zhao R, Jia P. Feasibility of Replacing the Range Doppler Equation of Spaceborne Synthetic Aperture Radar Considering Atmospheric Propagation Delay with a Rational Polynomial Coefficient Model. Sensors. 2020; 20(2):553. https://doi.org/10.3390/s20020553

Chicago/Turabian Style

Hou, Shasha, Yuancheng Huang, Guo Zhang, Ruishan Zhao, and Peng Jia. 2020. "Feasibility of Replacing the Range Doppler Equation of Spaceborne Synthetic Aperture Radar Considering Atmospheric Propagation Delay with a Rational Polynomial Coefficient Model" Sensors 20, no. 2: 553. https://doi.org/10.3390/s20020553

APA Style

Hou, S., Huang, Y., Zhang, G., Zhao, R., & Jia, P. (2020). Feasibility of Replacing the Range Doppler Equation of Spaceborne Synthetic Aperture Radar Considering Atmospheric Propagation Delay with a Rational Polynomial Coefficient Model. Sensors, 20(2), 553. https://doi.org/10.3390/s20020553

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