Next Article in Journal
A Width Measurement Method of Line Shape Based on Second Harmonic Peak and Modulation Amplitude
Previous Article in Journal
A Near InfraRed Emissive Chemosensor for Zn2+ and Phosphate Derivatives Based on a Di-(2-picolyl)amine-styrylflavylium Push-Pull Fluorophore
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Interpolation Methods in the Validation of Backscattering Coefficient Products

1
Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100101, China
2
College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(1), 469; https://doi.org/10.3390/s23010469
Submission received: 20 November 2022 / Revised: 26 December 2022 / Accepted: 28 December 2022 / Published: 1 January 2023
(This article belongs to the Section Environmental Sensing)

Abstract

:
Validation is the basis of synthetic aperture radar (SAR) image quantification applications. Based on the point target of the field site, the radiation characteristics of the backscattering coefficient image can be used to optimize the SAR imaging, and the product production system can be more closely targeted, to ensure the image product accuracy in the actual quantification application. In this study, the validation of the backscattering coefficient image was examined using calibrators, and the radiometric properties of the image were evaluated by extracting the radar cross-section of each point target. Bilinear interpolation and fast Fourier transform (FFT) interpolation methods were introduced for the local area interpolation of point targets, and the two methods were compared from the perspective of response function imaging and validation accuracy. The results show that the FFT interpolation method is more favorable for validating the backscattering coefficient.

1. Introduction

Satellite load calibration is the key to obtaining high-precision satellite products through which quantitative relationships between load observations and observed physical quantities can be established [1]. Validation uses ground truth values to assess the accuracy and uncertainty of remote sensing products and is an effective way to evaluate the quality, reliability, and applicability of remote sensing products [2,3]. Radiometric calibration is mainly used for single-look complex (SLC) data, which is based on reference targets with a known Radar Cross Section (RCS) to determine the calibration constant and thus obtain the backscattering coefficient value of the ground target. Nevertheless, in some cases only the backscattering coefficient product is available to the user, hence there is the need to evaluate the quality of the backscattering coefficient image. The validation is performed on the backscattering coefficient image using ground reference point targets to quantitatively evaluate the radiometric accuracy of the backscattering coefficient image which is very important for environmental parameter inversion and applications.
The backscattering coefficient refers to the radar scattering cross-section per unit irradiated area and is a measure of the result of the interaction of the incident electromagnetic wave with the ground target [4]. In addition to being related to radar system parameters, it mainly depends on the complex dielectric constant of the ground features, surface roughness, etc. [5]. The backscatter coefficient product is a synthetic aperture radar (SAR) load base radiation characteristic product, which is generally calculated from single-look complex (SLC) data using a calibration constant K [6]. The accuracy of the backscatter coefficient is the basis for carrying out the quantitative application of radar remote sensing, which usually uses calibrators deployed on the ground to extract the point target response energy from the image and compare it with the theoretical value to derive the accuracy of the image.
For point target response energy measurements, because point target peaks account for only a small number of pixels in the processed data, the target response within the peak neighborhood must be interpolated, and the sample window is expanded before analysis [7]. Prasad et al. [8] were the first to propose a fast Fourier transform (FFT) interpolation method and showed that the FFT method can significantly reduce the computational effort and estimate the intermediate values of discrete sequences within a very small error margin. Wang et al. [9] introduced common methods in SAR image interpolation, including linear interpolation, such as nearest-neighbor, bilinear, cubic convolution, FFT, and Sinc interpolation, which are commonly used in signal processing. Chen et al. [10], when using the peak method to measure the response energy, proposed that the FFT complementary zero method has a smaller root mean square error when interpolating Sinc curves compared to the cubic spline interpolation and will more closely fit the Sinc curves.
High-resolution SAR images are significantly affected by speckle noise, and different interpolation methods may significantly affect the extraction accuracy of the point target response energy. Therefore, this study explores the effects of different interpolation methods on the accuracy of the backscattering coefficient and provides suggestions on interpolation methods using experimental data obtained from the Common Application Support Platform of Land Observation Satellite for Country Civil Space Infrastructure (CASPLOS_CCSI) 1-m resolution C-SAR/01 Satellite Radiation Characterization Validation.

2. Methods

2.1. Backscattering Coefficient Product Validation Workflow

The process of validating the SAR image backscattering coefficient products is shown in Figure 1. If this is SLC data, the SLC data are processed to obtain the backscattering coefficient product using the SAR satellite calibration equation, calibration constant, incidence angle, and other related information. First, SAR-ground matching is undertaken. Affine transformation is performed on the latitude and longitude coordinates of the calibrators deployed in the field to obtain the row and column coordinates of the corresponding point target in the image. Subsequently, the local area around the point target is extracted and interpolated [11], and the RCS of the point target is calculated using the integration method [12]. Finally, this was compared with the theoretical RCS of the calibrator to calculate the authenticity check results.
A symmetric imaging window of 2 k × 2 k size was selected, for which the point target position was the center. Because the point target occupies fewer image elements in the SAR image and is stored as discrete data, resulting in the loss of a certain degree of original signal information, the sample window must be extended via interpolation to recover the point target response.
After interpolation, the signal-to-clutter ratio (SCR) value is calculated for each point target using the method described in Equation (1). Point targets with a SCR higher than 20 dB [13] are considered valid and can be used for subsequent validation.
S C R = σ i p sin θ σ i c δ a δ r
where σ i p is the peak value of the backscattering coefficient in the energy integral region of the i-th point target; σ i c is the average backscattering coefficient value of the clutter background around the i-th point target; θ i is the local incidence angle of the i-th point target; δ r and δ a are the range and azimuth pixel spacing respectively.
The measured values of RCS for the window integration method were calculated using Equation (2).
R C S integral = ( i A N A σ i N A N B i B N B σ i ) δ a δ r
where R C S integral is the measured radar scattering cross-sectional area of the point target based on the window integration, the unit is m 2 ; σ i is the calculated backscattering coefficient of the image element of the point target imaging window; A is the actual energy integration region of the point target, including the total number of image elements N A ; B is the background area, including the total number of image elements; and N B , δ r and δ a are the range and azimuth pixel spacing respectively.
The necessary ground reference information, such as the theoretical RCS value of the corner reflector, can be calculated based on the shape and edge length of the corner reflector and waveband information of the satellite-based SAR [14]. The RCS theoretical value of the corner reflector was compared [15] to the actual value extracted from the image to generate the radiation validation index of the SAR image and evaluate the image characteristics.
In this study, the following indicators were extracted to analyze the image radiation characteristics.
(1)
Relative radiation validation accuracy
Ideally, corner reflectors with the same nominal RCS value should have equal RCS measurements; relative radiation validation accuracy is obtained by statistical analysis of the measurements of multiple corner reflectors with the same nominal RCS. This indicator reflects the relative difference between measured values.
Δ R = i = 1 N ( R C S i R C S ¯ ) 2 N 1
where Δ R is the relative radiation validation accuracy, and N is the number of point targets used for validation. R C S i is the RCS measurement value of the i-th point target, and R C S ¯ is the mean RCS measurement value of the point targets.
(2)
Absolute radiation validation accuracy
The absolute radiometric validation explains the absolute difference between the measured and theoretical values. The difference between the measured RCS value of the point target in the backscattering coefficient image and the corresponding theoretical value was calculated, and the maximum value of the absolute value of the result was used as the result of the absolute radiation validation accuracy of the backscattering coefficient image.
Δ A = M a x ( | R C S 1 R C S t 1 | , | R C S 2 R C S t 2 | , , | R C S n R C S t n | )
where Δ A is the absolute radiation validation accuracy; R C S i and R C S t i are the RCS measurement and theoretical value of the target RCS at the i-th point, respectively.

2.2. Interpolation Methods

From the above analysis, it is clear that accurate interpolation is an important step in extracting the point target response energy for backscattering coefficient validation. Linear interpolation is weighted according to spatial proximity, and its error at the interpolation nodes is zero. It is simple and convenient compared to other interpolation methods [16]; therefore, it is widely used in remote sensing image processing. Among them, bilinear interpolation can realistically fit the original image, is computationally faster than cubic convolution [17], and is a common interpolation method in many image processing applications. However, the linear interpolation method only considers spatial proximity and ignores the fluctuating characteristics of the signal. By contrast, the FFT interpolation method increases the sampling frequency of the signal by complementing the high-frequency part of the frequency domain, which can make the signal wave smoother [18]. Therefore, in this study, bilinear and FFT interpolations were compared.

2.2.1. Bilinear Interpolation

Bilinear interpolation, also known as quadratic linear interpolation, is the process of bilinearly interpolating a pixel point by using the four pixel values around it plus the distance weights to sum up the new interpolation result [19]; Figure 2 shows the schematic diagram of bilinear interpolation.
The value taken at the place of ( x + u , y + v ) is calculated using Formula (5):
f ( x + u , y + v ) = ( 1 u ) ( 1 v ) f ( x , y ) + u ( 1 v ) f ( x + 1 , y )                   + v ( 1 u ) f ( x , y + 1 ) + u v f ( x + 1 , y + 1 )
where x and y are the row and column coordinates of the known points in the input image, respectively; u and v are the distance difference between the required point and the known point row and column coordinates respectively.

2.2.2. FFT Interpolation

In digital image processing, an FFT zero-padding operation is often performed on the signal. The FFT interpolation method performs a Fourier transform on the signal in the time domain and then zero padding at a high frequency, that is, the middle part of the signal in the frequency domain, producing the effect of interpolation in the time domain. In Figure 3, (a) identifies the point target before interpolation, (b) shows the result after zeroing the spectrum, and (c) shows the point target after interpolation. The white cross-shaped dashed box in (b) identifies the invalid spectral region that can be zeroed, and the expanded point target can be obtained by inversion of the 2D discrete Fourier transform after zeroing.
The interpolation process based on the FFT is realized using one Fourier transform and one Fourier inverse transform, and its operation process is divided into three main steps. First, a fast Fourier transform is applied to the original sequence X N ( k 1 , k 2 ) of size N × N [17].
X N ( k 1 , k 2 ) = n 1 = 0 N 1 n 2 = 0 N 1 x ( n 1 , n 2 ) exp { j 2 π ( k 1 n 1 N + k 2 n 2 N ) }         n 1 , n 2 , k 1 , k 2 [ 0 , N 1 ]
Then, a new sequence X M ( k 1 , k 2 ) of size M × M is constructed by X N ( k 1 , k 2 ) .
X M ( k 1 , k 2 ) = { L 2 × X N ( k 1 , k 2 ) k 1 [ 0 , N 1 2 ] , k 2 [ 0 , N 1 2 ] 0 k 1 [ N 1 2 + 1 , M N 1 2 1 ] , k 2 [ 0 , N 1 2 ] L 2 × X N ( k 1 M + N , k 2 ) k 1 [ N 1 2 + 1 , M 1 ] , k 2 [ 0 , N 1 2 ] 0 k 1 [ 0 , M 1 ] , k 2 [ N 1 2 + 1 , M N 1 2 1 ] L 2 × X N ( k 1 , k 2 ) k 1 [ 0 , N 1 2 ] , k 2 [ N 1 2 + 1 , M 1 ] 0 k 1 [ N 1 2 + 1 , M N 1 2 1 ] , k 2 [ N 1 2 + 1 , M 1 ] L 2 × X N ( k 1 M + N , k 2 ) k 1 [ N 1 2 + 1 , M 1 ] , k 2 [ N 1 2 + 1 , M 1 ]
After that, the Fourier inversion of the X M ( k 1 , k 2 ) is constructed to get x ^ ( m , m ) .
x ^ ( m 1 , m 2 ) = 1 M 1 M 2 k 1 = 0 M 1 1 k 2 = 0 M 2 1 X M ( k 1 , k 2 ) exp { + j 2 π ( k 1 m 1 M + k 2 m 2 M ) }
In Equations (6)–(8), x ( n 1 , n 2 ) is the original signal sequence, n 1 and n 2 refer to the image coordinates of the original signal sequence. N refers to the row and column size of the original signal. L is the interpolation multiplier. M refers to the row and column size of the sequence after interpolation. X N ( k 1 , k 2 ) is the two-dimensional sequence after fast Fourier transform. X M ( k 1 , k 2 ) is the frequency domain sequence after zero-padding. x ^ ( m 1 , m 2 ) is the interpolated signal sequence. m 1 and m 2 are the coordinates in the interpolated signal sequence.
The key to the FFT zero-padding interpolation method lies in the selection of the sub-window size and setting of the zero-padding multiplier. The response range of the calibrators on the image is different for different edge sizes; specifically, the sub-window needs to have the correct number of sampling points to ensure the effectiveness of the interpolation algorithm, but there should not be too many sampling points, otherwise it will affect computational efficiency [10].

3. Experimental Area and Data

This experiment, supported by CASPLOS_CCSI, relies on the Xilinhot SAR satellite calibration and validation. The surface of the field was mainly covered with kerchief fescue and sheep grass, and the terrain in the study area was flat and open. The scattering characteristics of the background features were uniform and stable without strong targets or electromagnetic interference, which can effectively reduce the interference of background clutter [20]. Seven corner reflectors were deployed for the experiment, and their distributions are shown in Figure 4.
The GaoFen-3 satellite [21] is an important Chinese SAR satellite, and it has played a significant role in applications within the fields of marine monitoring, disaster mitigation, environmental protection, water conservancy, agriculture and meteorology, etc. The C-SAR/01, launched on 23 November 2021, is the first follow-up satellite to GaoFen-3 and is designed to form part of the C-band SAR satellite constellation of China’s sea and land surveillance and monitoring system. The payload of C-SAR/01 is a C-band SAR with a center frequency of 5.4 GHz. It has 12 imaging modes such as spotlight, strip, and TOPSAR, which can acquire SAR images with a resolution of 1 m~500 m, swath of 10 km~650 km, and ranges from single polarization to full polarization to realize the monitoring of ocean and land resources. The imaging modes and parameters of C-SAR/01 are shown in Table 1.
The C-SAR/01 SAR data used in this study are listed in Table 2, and the data are Level-1A SLC images acquired on 12 May 2022. The imaging mode of the image is Ultra-Fine Strip (UFS) with a resolution of 3 m in both range and azimuth direction.
The results were obtained using Equation (9) to calculate the backscattering coefficient products in the experimental area [22], as shown in Figure 5.
σ d B 0 = 10 log 10 ( D N 2 ( Q u a l i f y V a l u e 32,767 ) 2 sin θ ) K d B
where σ d B 0 is the calculated backscattering coefficient value, and the unit is dB. D N is the digital number of the SAR image, D N 2 is obtained by summing the squares of the real and imaginary parts of the SLC data, i.e., D N 2 = I 2 + Q 2 . Q u a l i f y V a l u e and K d B are the image quantization maximum value and calibration constant respectively, which can be obtained from the image metadata file. θ is the local incidence angle, which can be obtained from the incidence angle file of C-SAR/01.
The seven corner reflectors used in this experiment were all developed for CASPLOS_CCSI, as shown in Figure 6. The inner leg length is 1 m. Plate curvature is less than ±2 mm. When the angle is adjusted, the deviation of azimuth angle and elevation angle is less than ±0.1°. Therefore, the corner reflectors used in the experiment can effectively guarantee RCS accuracy for radiometric validation.

4. Results and Discussion

According to Nyquist’s sampling theorem, if the signal is band-limited and the sampling frequency is higher than twice the signal bandwidth [23], the original continuous signal can be completely reconstructed from the sampled values. The SLC used in this experiment consists of a real part and imaginary part, which satisfies Nyquist’s sampling law [24]. Sinc interpolation [25] is a common method for complex signal processing [26,27] and can reconstruct as much continuous information as possible in discrete sampled images [28]. Therefore, for SLC data, the continuous information in the discrete sampled image can be reconstructed by Sinc interpolation.
However, in the backscattering coefficient, the pixel DN value only contains the amplitude information, which no longer satisfies the Nyquist sampling law, and thus the original information cannot be reconstructed by Sinc interpolation method. Therefore, we aim to find a suitable interpolation method, so that the information of the point target in the backscattering coefficient can be restored as much as possible to perform radiometric accuracy validation. Bilinear interpolation is a widely used conventional interpolation method in remote sensing image processing, simple but with high accuracy in most cases [29]. The FFT interpolation method is utilized for interpolating finite-time real sequences with the discrete Fourier transform. Compared to the classical interpolation method [30], the mean square error of FFT interpolation is significantly lower. Besides, the FFT interpolation method requires less computation and has a higher accuracy for long sequence data [8].
In order to evaluate the efficiency of bilinear and FFT interpolation methods, Sinc interpolation on the SLC data was used as a reference. In this experiment, after the backscattering coefficient data to be examined were located at the point target, based on the position information of the calibrator, a local area of 32 × 32 pixels around the corner reflector was selected, and the backscattering coefficient data within the sample window were subjected to 8-fold bilinear interpolation and FFT interpolation, respectively, and the two methods were compared.

4.1. Point Target Response Function Analysis

For the results obtained by the two interpolation methods, the profile lines were extracted along the range and azimuth directions of the peak points and compared with the Sinc interpolation results of the SLC data; the results of CR-1 are shown in Figure 7.
Observing the range and azimuth profile lines in the above figure, the image as a whole shows the same trend, but the results obtained by the bilinear interpolation method present jaggedness, which do not fit the Sinc interpolation results of the SLC image well, which is a defect of the algorithm itself. The discontinuity of the data may introduce errors in the window integration method for RCS calculation, while the FFT interpolation method is smoother in fitting the results, so the effect of the FFT interpolation method is better than the bilinear interpolation from the point target response function profile line.
The two-dimensional interpolation results for the (a) backscattering coefficient products, (b) Sinc interpolation of SLC data, (c) bilinear interpolation, and (d) FFT interpolation are shown in Figure 8.
Visually, the point target behaves as a sine-type function in the processed backscattering coefficient image [7], and the image is symmetric and shows fluctuations. Both the FFT interpolation and bilinear interpolation results of the backscattering coefficient are similar to the Sinc interpolation results of the SLC data, but the FFT interpolation results are smoother and better reflect the characteristics of the fluctuation of the point target signal itself, complementing some of the information loss due to resolution limitations. On the other hand, the bilinear interpolation results are calculated only from the perspective of spatial proximity, and although they fit the original backward scattering coefficient data more closely, the results are more rigid and lose some signal information, which may affect the RCS extraction accuracy.

4.2. RCS Extraction Results

The interpolation multiples affect the sampling resolution and spectral resolution, as well as the radiometric validation accuracy [31]. After comparing the validation accuracy indicators with different interpolation multiples, it shows the validation accuracy is the highest at eight times interpolation. Therefore, in this study, an eight-fold interpolation is uniformly performed on the local area of the point target. The point target SCR and RCS extracted after the Sinc interpolation of the SLC data, and the point target RCS extracted after bilinear and FFT interpolation of the backscattering coefficient are listed in Table 3.
From the graphical results, the FFT interpolation of the backscattering coefficient is closer to the theoretical value of the RCS, and the standard deviation is smaller. Corner reflectors appear as point targets on SAR images. According to Equations (6)–(8), FFT interpolation increases the sampling frequency by complementing zeros in the high frequency portion of the frequency domain, resulting in the signal being interpolated in the time domain [8]. Nevertheless, bilinear interpolation is based on the spatial proximity of the imaging, solving for unknown points based on adjacent known points. As can be seen from Figure 7 and Figure 8, the FFT interpolation method provides a better fit to the fluctuating characteristics of the original signal, and the FFT interpolation method can better restore the signal characteristics of point targets than bilinear interpolation. When the signal recovery in the local area of the point target is greater, the calculated RCS is more accurate. Therefore, FFT is more suitable for the interpolation of point target signals and thus more suitable for backscattering coefficient validation. The RCS values of each corner reflector extracted using Sinc, Bilinear and FFT interpolation methods are shown in Figure 9.
To compare the effects of different interpolation methods on the validation accuracy, the relative validation accuracy and absolute validation accuracy were extracted for Sinc, Bilinear and FFT interpolation methods according to Equations (3) and (4) and listed in Table 4.
By analyzing and comparing the parameters listed in the table, the following conclusions can be drawn: (1) From the relative validation accuracy shown in the table, it can be seen that the results obtained by both FFT interpolation and bilinear interpolation methods are stable, but the FFT interpolation method is slightly better. (2) The absolute validation accuracy is the maximum of all point target errors in each method, and the RCS errors extracted after FFT interpolation are smaller. In summary, the FFT interpolation method has somewhat higher accuracy in the backscattering coefficient validation.

5. Conclusions

Validation is a key step in the quantitative application of the SAR backscattering coefficient products. This paper briefly describes the backscattering coefficient validation process, combines the measured data, extracts the local area of the point target in the C-SAR/01 image for interpolation, calculates the radar scattering cross-sectional area of the point target using the integral method, and generates validation indexes. The interpolation effects of bilinear interpolation and FFT interpolation were compared using the response function and validation accuracy after interpolation.
From the point target response function, using the FFT interpolation method to interpolate the local signals around the calibrators, the interpolation result profile lines were more consistent with the point target characteristics, retaining the original features of the image, and the range and azimuth profile curves were smoother. In contrast, bilinear interpolation results in a jagged image that may cause the loss of some signals. Among the validation indices calculated after using different interpolation methods, the RCS of the calibrators obtained by the FFT interpolation method is closer to the theoretical value and more stable among individual point data, which truly reflects the radiation characteristics of the image. In subsequent studies, in-depth research can be conducted on how to suppress the ringing effect of the interpolated image by adding windows according to the FFT interpolation principle and reducing the influence of background clutter on RCS extraction to improve validation accuracy.

Author Contributions

Y.J. wrote the manuscript and designed the comparative experiments; F.Z. supervised the study and revised the manuscript; Q.H. and L.L. assisted Y.J. in designing the architecture and conducting experiments. X.L. assisted Y.J. in obtaining experimental data. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the China high-resolution earth observation system (21-Y20B01-9003-19/22) and the Common Application Support Platform for Land Observation Satellites of China’s Civil Space Infrastructure (CASPLOS_CCSI).

Acknowledgments

Thanks are given for the data support of all members of the North China Comprehensive Experimental Group of the National Civil Space Infrastructure Terrestrial Observation Satellite Common Application Support Platform Project. In addition, the author would like to thank anonymous reviewers for their constructive comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yuan, X. The Radiometric Calibration of Satellite Borne Synthetic Aperture Radar. Aerosp. Shanghai 1998, 15, 7. [Google Scholar]
  2. Jin, R.; Ran, Y.; Liu, Q.; Zhang, R.; Li, X.; Ma, G.; Ge, Y.; Liu, S.; Xiao, Q. Key Methods and Experiment Verification for the Validation of Quantitative Remote Sensing Products. Adv. Earth Sci. 2017, 32, 630–642. [Google Scholar]
  3. Jiang, X.; Li, Z.; Xi, X.; Li, X. Basic Frame of Remote Sensing Validation System. Arid Land Geogr. 2008, 4, 567–571. [Google Scholar]
  4. Ulaby, F.T. Microwave Remote Sensing; Scientific Publishing House: Beijing, China, 1987. [Google Scholar]
  5. Guo, H. Radar Earth Observation Theory and Applications; Scientific Publishing House: Beijing, China, 2000. [Google Scholar]
  6. Freeman, A.J.I.T. SAR Calibration: An Overview. IEEE Trans. Geosci. Remote Sens. 1992, 30, 1107–1121. [Google Scholar] [CrossRef]
  7. Wong, L.G.C.F.H. Digital Processing of Synthetic Aperture Radar Data Algorithms and Implementation; Publishing House of Electronics Industry: Beijing, China, 2012. [Google Scholar]
  8. Prasad, K.P.; Satyanarayana, P.J.E.L. Fast Interpolation Algorithm Using FFT. Electron. Lett. 1986, 22, 185–187. [Google Scholar] [CrossRef]
  9. Wang, H.; Chen, X. Research on Fast Fourier Transform Algorithm Based on Interpolation Re-sampling Method. Water Power 2012, 38, 3. [Google Scholar]
  10. Chen, J.; Zhang, B.; Wang, C.; Lei, W. Comparative Study on the Applicability of Two Corner Reflector based Radiometric Calibration Methods for High Resolution Airborne SAR Image. Remote Sens. Technol. Appl. 2015, 30, 677–683. [Google Scholar]
  11. Freeman, A.; Shen, Y.; Werner, C.L. Polarimetric SAR Calibration Experiment Using Active Radar Calibrators. IEEE Trans. Geosci. Remote Sens. 1990, 28, 224–240. [Google Scholar] [CrossRef]
  12. Gray, A.L.; Vachon, P.W. Synthetic Aperture Radar Calibration Using Reference Reflectors. IEEE Trans. Geosci. Remote Sens. 1990, 28, 374–383. [Google Scholar] [CrossRef]
  13. Shimada, M.; Oaku, H.; Nakai, M. SAR Calibration Using Frequency-Tunable Active Radar Calibrators. IEEE Trans. Geosci. Remote Sens. 1999, 37, 564–573. [Google Scholar] [CrossRef]
  14. Zhang, T.; Zhang, P.; Zeng, Q. Study on Corner Reflectors in SAR Calibration. Remote Sens. Inf. 2010, 3, 38–42. [Google Scholar]
  15. Zhou, X.; Zeng, Q.; Jiao, J.; Wang, Q.; Xiong, S.; Gao, S. Research on Space-borne SAR Field Calibration Experiment—A Case Study of TerraSAR-X Field Calibration. Remote Sens. Technol. Appl. 2014, 5, 711–718. [Google Scholar]
  16. Wang, J.; Li, H.; Wang, J.; Miao, Z. A Scheme and Its Analysis for Fast Linear Image Interpolation. Acta Electron. Sin. 2009, 37, 1481–1486. [Google Scholar]
  17. Guo, X. A Study of the Filtering Algorithm for Quasi-bilinear Interpolation Based on Curvature. Master’s Thesis, Hebei Normal University, Shijiazhuang, China, 2013. [Google Scholar]
  18. Zhang, Y.; Chen, L.; Hao, H.; Zheng, G.; Liang, J. Improved Algorithm for Interpolation Based on FFT. J. Data Acquis. Process. 2013, 2, 173–177. [Google Scholar]
  19. Bailey, K. A Novel Approach to Real-time Bilinear Interpolation. In Proceedings of the IEEE International Workshop on Electronic Design, Perth, WA, Australia, 28–30 January 2004. [Google Scholar]
  20. Sharma, S.; Dadhich, G.; Rambhia, M.; Mathur, A.K.; Prajapati, R.P.; Patel, P.R.; Shukla, A. Radiometric Calibration Stability Assessment for the RISAT-1 SAR Sensor Using a Deployed Point Target Array at the Desalpar Site, Rann of Kutch, India. Int. J. Remote Sens. 2017, 38, 7242–7259. [Google Scholar] [CrossRef]
  21. Sun, J.; Yu, W.; Deng, Y. The SAR Payload Design and Performance for the GF-3 Mission. Sensors 2017, 17, 2419. [Google Scholar] [CrossRef] [Green Version]
  22. Zheng, M.; Li, X.; Ren, Y. The Method Study on Automatic Sea Ice Detection with Gao Fen-3 Synthetic Aperture Radar Data in Polar Regions. J. Oceanogr. 2018, 9, 113–124. [Google Scholar]
  23. Yu, W.; Wu, S. Comparison of Several Interpolation Algorithms in Range-doppler Method. J. Electron. Inf. Technol. 2001, 3, 308–312. [Google Scholar]
  24. Yang, X.; Lou, C.; Xu, J.; Zhao, Y. Software-Defined Radio: Principles and Practice, 1st ed.; Electronic Industry Publishing House: Beijing, China, 2007; pp. 8–20. [Google Scholar]
  25. Zhang, G.; Jia, X.; Yin, C.; Li, Y. Comparison of Interpolation Algorithms in Synthetic Aperture Radar Imaging Processing. Foreign Electron. Meas. Technol. 2008, 8, 9–11. [Google Scholar]
  26. Oliver, C.; Quegan, S. Understanding Synthetic Aperture Radar Images; Publishing House of Electronics Industry: Beijing, China, 2004. [Google Scholar]
  27. Schanze, T. Sinc Interpolation of Discrete Periodic Signals. EEE Trans. Signal Process 1995, 43, 1502–1503. [Google Scholar] [CrossRef]
  28. Oppenheim, A.V.; Schafer, R.W.; Buck, J.R. Discrete-Time Signal Processing, 3rd ed.; Publishing House of Electronics Industry: Beijing, China, 2015. [Google Scholar]
  29. Han, D. Comparison of Commonly Used Image Interpolation Methods. In Proceedings of the ICCSEE 2013, Hangzhou, China, 22 March 2013. [Google Scholar]
  30. Proakis, J.G.; Manolakis, D.G. Digital Signal Process, 4th ed.; Publishing House of Electronic Industry: Beijing, China, 2007. [Google Scholar]
  31. Li, L.; Zhang, F.; Shao, Y.; Wei, Q.; Huang, Q.; Jiao, Y. Airborne SAR Radiometric Calibration Based on Improved Sliding Window Integral Method. Sensors 2022, 22, 320. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The validation process for the backscattering coefficient.
Figure 1. The validation process for the backscattering coefficient.
Sensors 23 00469 g001
Figure 2. Schematic diagram of the bilinear interpolation method.
Figure 2. Schematic diagram of the bilinear interpolation method.
Sensors 23 00469 g002
Figure 3. Schematic diagram of FFT interpolation. (a) Point target before interpolation. (b) Frequency graph after zero padding. (c) Point target after interpolation.
Figure 3. Schematic diagram of FFT interpolation. (a) Point target before interpolation. (b) Frequency graph after zero padding. (c) Point target after interpolation.
Sensors 23 00469 g003
Figure 4. Study area coverage.
Figure 4. Study area coverage.
Sensors 23 00469 g004
Figure 5. Geocoded backscattering coefficient products.
Figure 5. Geocoded backscattering coefficient products.
Sensors 23 00469 g005
Figure 6. Automatic trihedral corner reflectors.
Figure 6. Automatic trihedral corner reflectors.
Sensors 23 00469 g006
Figure 7. Interpolation results along range and azimuth direction using Sinc, Bilinear and FFT methods for point target (taking CR-1 as an example). (a) Interpolation results along range direction. (b) Interpolation results along azimuth direction.
Figure 7. Interpolation results along range and azimuth direction using Sinc, Bilinear and FFT methods for point target (taking CR-1 as an example). (a) Interpolation results along range direction. (b) Interpolation results along azimuth direction.
Sensors 23 00469 g007
Figure 8. The backscattering coefficient image of CR-1 and the two-dimensional interpolation results using Sinc, Bilinear and FFT method. (a) Raw backscattering coefficient data; (b) Sinc interpolation results of SLC data; (c) Bilinear interpolation results of the backscattering coefficient; (d) FFT interpolation results of the backscattering coefficient.
Figure 8. The backscattering coefficient image of CR-1 and the two-dimensional interpolation results using Sinc, Bilinear and FFT method. (a) Raw backscattering coefficient data; (b) Sinc interpolation results of SLC data; (c) Bilinear interpolation results of the backscattering coefficient; (d) FFT interpolation results of the backscattering coefficient.
Sensors 23 00469 g008aSensors 23 00469 g008b
Figure 9. The extracted RCS values of each corner reflector using Sinc, Bilinear and FFT interpolation methods.
Figure 9. The extracted RCS values of each corner reflector using Sinc, Bilinear and FFT interpolation methods.
Sensors 23 00469 g009
Table 1. Imaging modes and parameters of C-SAR/01.
Table 1. Imaging modes and parameters of C-SAR/01.
Imaging ModeNominal Resolution/mSwath Width/kmPolarization
Spotlight (SL)110HH/VV
StripmapUltra-Fine Strip (UFS)330HH/VV
Fine Strip I (FSI)550HH + HV/VV + VH
Fine Strip II (FSII)10100HH + HV/VV + VH
Standard Strip (SS)25130HH + HV/VV + VH
Quad Polarization Strip I (QPSI)830Quad-polarization
Quad Polarization Strip II (QPSII)2540Quad-polarization
ScanNarrow ScanSAR (NSC)50300HH + HV/VV + VH
Wide ScanSAR (WSC)100500HH + HV/VV + VH
Global (GLO)500650HH + HV/VV + VH
Wave imaging (WAV)820Quad-polarization
Expanded incidence angle (EXT)Low Incidence25130HH + HV/VV + VH
High Incidence2580HH + HV/VV + VH
Table 2. Imaging parameters of C-SAR/01 image.
Table 2. Imaging parameters of C-SAR/01 image.
Date12 May 2022
Ascending or DescendingDescending
Look DirectionRight
Incidence Angle28.43–30.57°
Imaging ModeUFS
Polarization ModeHH
Range Resolution/m3
Azimuth Resolution/m3
Range pixel space/m1.124222
Azimuth pixel space/m1.669818
Table 3. SCR and RCS calculated for different interpolation methods.
Table 3. SCR and RCS calculated for different interpolation methods.
CR
Number
SCR/dBTheoretical Value/dBSinc Interpolation of SLC Data/dBBackscattering Coefficient/dB
Bilinear InterpolationFFT Interpolation
CR-136.90331.33231.12831.04630.951
CR-238.18131.33231.80831.70831.63
CR-338.26231.33231.74831.6531.535
CR-437.39831.33231.70431.60631.492
CR-538.57231.33231.75731.66231.561
CR-637.67831.33231.55931.4631.338
CR-739.90531.33231.87731.79131.672
Table 4. Validation accuracy of different interpolation results.
Table 4. Validation accuracy of different interpolation results.
Relative Validation Accuracy/dBAbsolute Validation Accuracy/dB
Sinc interpolation of SLC data0.2520.544
Bilinear interpolation of backscattering coefficient0.230.459
FFT interpolation of backscattering coefficient0.2280.381
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jiao, Y.; Zhang, F.; Huang, Q.; Liu, X.; Li, L. Analysis of Interpolation Methods in the Validation of Backscattering Coefficient Products. Sensors 2023, 23, 469. https://doi.org/10.3390/s23010469

AMA Style

Jiao Y, Zhang F, Huang Q, Liu X, Li L. Analysis of Interpolation Methods in the Validation of Backscattering Coefficient Products. Sensors. 2023; 23(1):469. https://doi.org/10.3390/s23010469

Chicago/Turabian Style

Jiao, Yanan, Fengli Zhang, Qiqi Huang, Xiaochen Liu, and Lu Li. 2023. "Analysis of Interpolation Methods in the Validation of Backscattering Coefficient Products" Sensors 23, no. 1: 469. https://doi.org/10.3390/s23010469

APA Style

Jiao, Y., Zhang, F., Huang, Q., Liu, X., & Li, L. (2023). Analysis of Interpolation Methods in the Validation of Backscattering Coefficient Products. Sensors, 23(1), 469. https://doi.org/10.3390/s23010469

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