Next Article in Journal
Vegetation Cover Dynamics in the High Atlas Mountains of Morocco
Previous Article in Journal
The Forest Fire Dynamic Change Influencing Factors and the Impacts on Gross Primary Productivity in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Channel Imbalance Calibration Scheme with Distributed Targets for Circular Quad-Polarization SAR with Reciprocal Crosstalk

1
The Department of Space Microwave Remote Sensing System, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100190, China
2
The School of Electronic, Electrical and Communication Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Submission received: 7 February 2023 / Revised: 26 February 2023 / Accepted: 26 February 2023 / Published: 28 February 2023

Abstract

:
As polarimetric antennas can be isolated through excellent electronic frameworks in circular quad-polarization (CQP) synthetic aperture radar (SAR) systems, cross-polarization (x-pol) and co-polarization (co-pol) channel imbalances are more challenging and essential to calibrate than crosstalk in polarimetric calibration (PolCAL). In uncalibrated CQP SAR images without corner reflectors (CRs), the reciprocity and reflection symmetry assumptions of the distributed targets are commonly used to estimate the x-pol and co-pol channel imbalances, respectively. To suppress the influence of additive noise on determining channel imbalances through distributed targets, high signal-to-noise ratio (SNR) distributed targets should be obtained from the x-pol and co-pol channels of the CQP SAR images: namely, surface-dominated and volume-dominated targets. However, some reflection symmetry assumptions used in the existing calibration literature have poor applicability with volume-dominated targets, resulting in unsatisfactory estimation results for the co-pol channel imbalance phase. In this paper, we assess the priority of the reflection symmetry properties of volume-dominated targets used to calibrate the co-pol channel imbalance phase in CQP SAR data synthesized from linear quad-polarization data of ALOS, GF-3, and RADARSAT-2. In the theoretical part, high-priority reflection symmetry (termed semireflection symmetry) assumptions are confirmed as the most suitable for estimating the co-pol channel imbalance phase, and were selected to develop an algorithm for estimating the co-pol channel imbalance phase. Furthermore, based on the novel method for estimating the co-pol channel imbalance phase, a channel imbalance calibration scheme is proposed for CQP SAR systems with reciprocal crosstalk, including extracting surface-dominated and volume-dominated targets, and estimating and filtering channel imbalances. We demonstrate the effectiveness of our proposed scheme with CRs in simulated CQP SAR images. The experimental results show that the calibration scheme is an effective workflow for estimating channel imbalances in CQP SAR systems with reciprocal crosstalk.

Graphical Abstract

1. Introduction

Quad-polarization synthetic aperture radar (SAR) systems actively transmit and receive different combinations of electromagnetic waves, including (1) linear quad-polarization (LQP) SAR, which alternately transmits and receives horizontal (H) and vertical (V) electromagnetic waves, (2) circular quad-polarization (CQP) SAR, which alternately transmits and receives left-handed (L) and right-handed (R) electromagnetic waves, and (3) hybrid quad-polarization (HQP) SAR, which transmits L and R electromagnetic waves, according to combinations of different phases of H and V electromagnetic waves, and alternately receives H and V waves [1,2]. Due to its high sensitivity to clutter backscatter, quad-polarization SAR provides high-quality data for a wide range of remote applications [3,4,5,6,7]. However, hardware system distortions can lead to unfaithful representations of clutter backscatter. It is very important to calibrate the observation matrix to the undistorted measurement, which is known as polarimetric calibration (PolCAL) [8,9,10]. Polarimetric distortions, which are SAR system parameters, mainly include crosstalk, co-polarization (co-pol) channel imbalance, and cross-polarization (x-pol) channel imbalance [11,12,13,14,15]. For decades, SAR antennas have been highly isolated (better than 32 dB [16,17,18]), with crosstalk that is stable at different incidence angles. Therefore, x-pol and co-pol channel imbalances in PolCAL merit further study.
Although CQP SAR systems have been extensively studied in the imaging radar literature [19], channel imbalance calibration on CQP SAR images has received less attention. Based on the antenna transmit model s-m-a (which minimizes the potential for crosstalk, by ensuring complete physical isolation between channels) and the corresponding receive model, [20] proposed an iterative method to calibrate crosstalk and channel imbalances, using corner reflector (CR) groups (including trihedral, dihedral, and gridded trihedral CRs) in CQP SAR systems with reciprocal crosstalk. In addition, considering that CRs can be used only to determine the distortion of CQP SAR data at specific incidence angles, the reflection symmetry and reciprocity of the distributed targets were also considered when estimating the co-pol and x-pol channel imbalances. As the reciprocity of polarimetric SAR images is generally accepted in polarimetric applications [19,21], and the estimation of the x-pol channel imbalance with reciprocally distributed targets is not perturbed by the reciprocal crosstalk in the CQP SAR calibration, the initial estimate of the x-pol channel imbalance was equal to the true value [20]. However, the statistical reflection symmetry characteristics exhibited by distributed targets in the undistorted CQP SAR images did not fully match the ideal reflection symmetry assumptions, which limited the calibration accuracy of the co-pol channel imbalance. Compared with the CR calibration results, the co-pol channel imbalance phase estimated according to the reflection symmetry properties of the distributed targets differed by 14 degrees ([20], Table 2).In existing LQP systems, the acceptable error of the calibrated co-pol channel imbalance phase is generally within ±10 degrees ([22], Table 4).
The first contribution of this paper is to propose a calibration algorithm for the co-pol channel imbalance phase based on high-priority reflection symmetry (termed semireflection symmetry) assumptions. As additive noise may affect the statistical polarization characteristics of the distributed targets [18,23], high signal-to-noise ratio (SNR) targets in the CQP SAR co-pol channel—namely, volume-dominated targets—should be applied, to calibrate the co-pol channel imbalance. One reflection symmetry condition of LQP SAR images is that the average polarization orientation angle is zero, which can be difficult to achieve for natural targets with large terrain fluctuations and complex scattering characteristics, such as sparse woodland areas [24,25]. Considering that the reflection symmetry of CQP SAR images can be obtained from LQP SAR images through polarimetric basis transformations, the calibration accuracy cannot be guaranteed if the low-priority reflection symmetry properties of volume-dominated targets are used to address the co-pol channel imbalance phase. Therefore, this paper determines the priority of reflection symmetry assumptions in volume-dominated targets with CQP SAR data, and then confirms that semireflection symmetry assumptions are the most suitable for estimating the co-channel imbalance phase. In addition, a calibration method for the co-pol channel imbalance phase is proposed, based on semireflection symmetry properties.
The second contribution of this paper is to provide a channel imbalance calibration scheme for CQP SAR systems with reciprocal crosstalk. Firstly, the equivalent number of looks (ENL) and the power ratio of the co-pol and x-pol channels are utilized, to extract volume-dominated and surface-dominated target regions in uncalibrated CQP SAR images. Secondly, the calibration method, via clutter [20], is introduced, to estimate the x-pol channel imbalance according to the surface-dominated targets, under the condition of reciprocal CQP SAR images and crosstalk. Thirdly, the semireflection symmetry assumptions are implemented, to remove the co-pol channel imbalance. Filter operations are also carried out in the second and third steps, to acquire more robust results than the initial estimates.
The effectiveness of the proposed scheme for calibrating channel imbalances is demonstrated with undistorted CQP SAR images synthesized from LQP SAR images, which were acquired from the Yudaokou airborne flight experiment of the Aerospace Information Research Institute, Chinese Academy of Sciences (AIRCAS). The experiments in this paper show that the proposed scheme can be used to estimate channel imbalances with distributed targets in CQP SAR systems with reciprocal crosstalk. In addition, the proposed method is an effective tool for improving the calibration accuracy of the co-pol channel imbalance phase.
The remainder of this paper is organized as follows. Section 2 introduces the basis of CQP SAR PolCAL, including the distortion model and the priority of the reflection symmetry assumptions. The channel imbalance calibration scheme is presented in Section 3. Section 4 introduces the experimental approach. Some artificially adjusted parameters, and the crosstalk tolerance of the proposed scheme, are discussed in Section 5. Finally, Section 6 provides the conclusions of this paper, and discusses future work.

2. Basis of CQP SAR PolCAL

2.1. Distortion Model

In CQP SAR PolCAL, the undistorted backscatter matrix S relates the observed backscatter matrix O to the transmit distortion matrix T , the receive distortion matrix R , and the random additive thermal noise N , according to the following distortion model [20]:
O l l O l r O r l O r r = r l l r l r r r l r r r S l l S l r S r l S r r t l l t l r t r l t r r + N l l N l r N r l N r r O = R S T + N
where the subscripts l or r indicate left-handed or right-handed electromagnetic waves. S p q or O p q (p or q = l or r, similarly hereinafter) denote the undistorted or distorted complex scattering factor transmitted by antenna q and received by antenna p; t p q denotes the product of the common amplifier factor at antenna q, and the scaling factor when antenna q transmits polarization p; r p q denotes the product of the common amplifier factor at antenna p, and the scaling factor when antenna p receives polarization q; n p q denotes the additive noise in the receive channel, p, when transmitting q. Without loss of generality, we assume that the target backscattering and additive noise are uncorrelated, and that the components in N are independent of one another [20,24]. Considering that the Faraday rotation angle can be well-estimated by external reference [22,24], we do not consider the impact of the Faraday rotation angle on (1). Row-based factorization is applied, to transform O , S , and N to O 4 LR , S 4 LR , and N 4 LR :
O l l O l r O r l O r r = Y k 2 0 0 0 0 k 0 0 0 0 k 0 0 0 0 1 α 0 0 0 0 1 0 0 0 0 α 0 0 0 0 1 1 v w v w z 1 z w w u u v 1 v u z u z 1 S l l S l r S r l S r r + N l l N l r N r l N r r O 4 LR = Y K 4 Q X S 4 LR + N 4 LR
where k, α , and u / v / w / z were the co-pol channel imbalance, x-pol channel imbalance, and crosstalk components, respectively; K 4 , Q , and X are the corresponding co-pol channel imbalance, x-pol channel imbalance, and crosstalk matrices, respectively; and Y represents the absolute calibration factor. The relationship between (1) and (2) can be expressed as
Y = r r r t r r , k = r l l / r r r , α = t l l / t r r r l l / r r r , u = r r l / r r r , v = t r l / t l l , w = r l r / r l l , z = t l r / t r r
In PolCAL,Y is usually considered to not impair the relative relationship of polarimetric channels [20,22,24,26]; therefore, the calibration process of Y is not considered in this article. The coherence information of each polarimetric channel can be adequately described by CQP SAR systems, which provide constraints for distributed target calibration. The corresponding covariance matrix (2) can be written as
M 4 LR = K 4 Q X C 4 LR K 4 Q X + C n
where M 4 LR , C 4 LR , and C n are the covariance matrices corresponding to O 4 LR , S 4 LR , and N 4 LR , respectively, and the superscript † represents the conjugate and transpose matrix operations. It should be stressed that C n is a four-by-four diagonal matrix, in which the first and second elements are equal, and the third and fourth elements are equal [20]. As multiplication of diagonal matrices is commutative, (4) can be transformed as
M 4 LR = Q K 4 X C 4 LR Q K 4 X + C n
In this paper, (2) and (5) are utilized to solve the channel imbalances with reciprocal crosstalk. When reciprocal crosstalk and noise are ignored, the x-pol and co-pol channel imbalances can be estimated by the reciprocity and reflection symmetry of distributed targets [20], which can be written as:
α = M 33 / M 22
α = M 32
k = M 11 / M 44 α 2 4
k = 1 2 M 13 / M 42 + n π n = 0 , ± 1 , ± 2 ,
k = 1 2 M 12 / M 43 α / α * + n π n = 0 , ± 1 , ± 2 ,
k = 1 2 M 14 α + 1 sgn C 14 2 + n π n = 0 , ± 1 , ± 2 ,
where · and · represent the phase and amplitude, respectively, of a complex number; M i j and C i j are the ith row and jth column components of M 4 LR and C 4 LR ; and the superscript ∗ indicates the conjugate operator.

2.2. Priority Verification of Reflection Symmetry in CQP SAR Images

The distributed targets with reflection symmetry are composed of single targets with mean-zero helicity and mean-zero orientation [9,27]. When the repeated equations, due to the reciprocity of the CQP SAR images, are removed, the reflection symmetry assumptions of the distributed targets can be expressed as
Lin 1 : S h h S h v * = 0 Lin 2 : S v v S h v * = 0 Line Basis Cir 1 : S l l S l l * S r r S r r * = 0 Cir 2 : S l l S r r * = 0 Cir 3 : S l r * S l l + S r r = 0 Circular Basis
where · is the spatial average operator of a matrix, which is applied to remove coherent speckle noise; · denotes the imaginary operation; and Lin i / Cir j denote the i/jth reflection symmetry assumption of the LQP/CQP SAR images. The circular basis in (12) can be obtained using a unitary basis-change matrix U l r , via a unitary consimilarity transformation [28], as follows:
S 4 LR = 1 2 1 j j 1 1 2 1 j j 1 T S h h S h v S v h S v v T S 4 LR = U l r U l r T S 4 HV
where ⊗ represents the Kronecker product of a matrix, and the superscript T represents the transpose operator.

2.2.1. Verification of Reflection Symmetry Used in the k Calibration Methods of [20]

In [20], Cir 2 and Cir 3 were used to estimate the phase of k through (9) to (11). Next, we rewrite Cir 2 and Cir 3 as ACir 2 and ACir 3 :
ACir 2 : S l l S r r * = n π , n = 0 , ± 1 , ± 2
ACir 3 : S l l S l r * S r r S l r * = 2 n π , n = 0 , ± 1 , ± 2
to more directly judge whether the reflection symmetry is satisfied in the calibration methods of k.
When calibrating different channel imbalances, high SNR targets must be selected in various channels: thus, the calibration results are less affected by additive noise. By analyzing the process of transforming LQP SAR images into CQP SAR images (13), we find that energy-dominated targets in co-pol channels in CQP SAR images exhibit volume-scattering features; therefore, volume-dominated targets should be used to calibrate k by applying reflection symmetry. Next, nine LQP SAR images from L-band ALOS, C-band GaoFen-3 (GF-3), and C-band RADARSAT-2 (RD-2) are used to analyze the reflection symmetry of the volume-dominated targets. The specific details include the following steps: (1) the volume-dominated targets are obtained by applying the optimal non-negative eigenvalue decomposition (NNED) [29], which is an effective method of reducing the negative power problem and the overestimated contribution of the volume scattering affected by Freeman—Durden decomposition [30]; (2) Cir 1 / 2 Span , ACir 2 , and ACir 3 are calculated by (12), (14) and (15), respectively; (3) histograms of Cir 1 / 2 Span and phase distributions of ACir 2 and ACir 3 are obtained, to confirm the difference between the assumptions exhibited by the actual SAR systems and the ideal theoretical assumptions. The final verification results are shown in Figure 1.
In Figure 1a, we illustrate the Cir 1 / 2 Span histogram of nine images, where Cir 1 values are usually small but nonzero compared to the span value. The ACir 2 and ACir 3 phase distributions of the volume-dominated targets are shown in Figure 1b,c, which are not constant at 0 ° or ±180 ° , but have approximately equal probability at all angles. Therefore, unsatisfactory k-angle calibration results may be obtained if ACir 2 and ACir 3 are used in volume-dominated areas.

2.2.2. Cir 3 Division for Reflection Symmetry Verification

In the common method, Cir 3 , which is represented by (12), can be decomposed into an amplitude and a phase, and the phase of Cir 3 is applied, to determine the k phase ((9) and (10)). However, the amplitude of Cir 3 also presents some information about the reflection symmetry. In Figure 2, Cir 3 / 2 Span (solid lines) is investigated in the volume-dominated targets, and Cir 1 / 2 Span (dotted lines) is also shown, for comparison. Clearly, Cir 3 and Cir 1 have the same order of small amplitude in most volume-dominated targets. Compared with ACir 3 , which has a similar distribution across the entire phase domain, Cir 3 may be a better choice for calculating the k phase. However, when Cir 3 / 2 Span and Cir 1 / 2 Span are close to 0, the frequency of Cir 1 / 2 Span is noticeably greater than that of Cir 3 / 2 Span . Thus, there may be a better method of calibrating the phase of k than directly using Cir 3 .
In addition to dividing a complex number into its amplitude and phase, complex numbers can be decomposed into imaginary and real parts. Compared to Cir 3 , the imaginary and real parts of Cir 3 may be less powerful than those of Cir 3 . Thus, Cir 3 can be divided into an imaginary part, Cir 4 , and a real part, Cir 5 , to achieve a more accurate description of the reflection symmetry in the volume-dominated targets.
The histograms of Cir 4 / 2 Span (dashed line) and Cir 5 / 2 Span (solid line) are shown in Figure 3a. Obviously, in the volume-dominated area, Cir 4 / 2 Span and Cir 5 / 2 Span have the same number level of points close to 0, and more points close to 0 than Cir 3 / 2 Span . Furthermore, the percentages of all volume-dominated points with Cir 4 / 2 Span < 0.1 and Cir 5 / 2 Span < 0.1 in each image are plotted in Figure 3b. Both have high percentages above 98% in all nine images. Additionally, the percentage of Cir 4 / 2 Span < 0.1 is higher than that of Cir 5 / 2 Span < 0.1 . Thus, the simulation experiments with the undistorted CQP SAR images show that Cir 4 and Cir 5 are closer to the ideal reflection symmetry assumptions than ACir 2 , ACir 2 , or Cir 3 . Additionally, while the value of Cir 4 near 0 is slightly higher than that of Cir 5 in the undistorted CQP SAR images, the difference is small.
Based on the above analysis, Cir 4 and Cir 5 are more suitable for estimating the phase of k by analyzing the prioritized reflection symmetry in undistorted volume-dominated areas. Furthermore, the priority between Cir 4 and Cir 5 should be determined in the reflection symmetry and reciprocity simulations of CQP SAR data with distortions (i.e., reciprocal crosstalk and α ), to obtain a higher k phase calibration accuracy. In the following, the influence of α is ignored, because the calibration results of estimating the x-pol channel imbalance with (6) and (7) in the reciprocal targets are the same with true value when the CQP SAR antennas have reciprocal crosstalk [20]. It should be stressed that reciprocity of targets is generally acceptable in polarimetric applications.
For this part, 1000 reflection symmetry points were simulated, to conduct 20,000 Monte Carlo trials. During each trial, we imposed the reciprocal crosstalk unknowns, with random amplitudes from 0 to 0.3 ( 0 u = z , v = w 0.3 10.4576 dB ), and random phases ( 180 ° u = z , v = w 180 ° ) were added to the simulated data first of all. Then, the percentages of Cir 4 / 2 Span < 0.1 and Cir 5 / 2 Span < 0.1 were counted in the presence of reciprocal crosstalk. Finally, we verified the priority of Cir 4 / 2 Span and Cir 5 / 2 Span , according to the maximum crosstalk, ( dB u , dB v , dB w , dB z max ).
The results are shown in Figure 3c,d, where the X axis is the maximum added crosstalk, and the Y axis is the percentages of Cir 4 / 2 Span < 0.1 and Cir 5 / 2 Span < 0.1 . As the crosstalk increases, an overall decrease can be seen in the ratios of Cir 4 / 2 Span < 0.1 and Cir 5 / 2 Span < 0.1 . However, at approximately −60 dB, the drop rate and drop value of Cir 5 / 2 Span are noticeably higher than those of Cir 4 / 2 Span . When the maximum crosstalk increases to approximately −15 dB, the percentage of Cir 5 / 2 Span larger than 0.1 is approximately 0.1 (10%), while the percentage of Cir 4 / 2 Span larger than 0.1 is approximately 0.85 (85%). Therefore, Cir 4 is more important than Cir 5 , when calibrating the k phase in the actual calibrated and simulated uncalibrated data, and Cir 4 is used to develop an algorithm to estimate the k phase.

3. Materials and Methods

Due to the high cost of deploying CRs, it is necessary to use distributed targets to determine the channel imbalances. This section describes the complete scheme of calibrating channel imbalances in CQP SAR images with reciprocal crosstalk. Based on some assumptions of CQP SAR systems and distributed targets, the proposed scheme consists of three main steps. Firstly, ENL and R v b are utilized to extract the surface-dominated and volume-dominated targets as the calibration candidate pixels. Secondly, the reciprocity of the surface-dominated targets is utilized, to determine the x-pol channel imbalance α . Thirdly, the high-priority reflection symmetry assumptions of the surface-dominated targets are utilized, to determine the co-pol channel imbalance k. In the second and third steps, filter operations are performed, to obtain more robust results. The schematic workflow is shown in Figure 4.

3.1. Assumptions of Natural Targets and CQP SAR System

3.1.1. Reciprocity of Natural Targets and Semireflection Symmetry of Volume-Dominated Targets

The first assumption for natural targets is reciprocity, which is a well-recognized assumption that has been applied in many polarization applications [31]. In LQP SAR systems, the reciprocity of natural targets can be expressed as S h v = S v h . Combined with (13), the reciprocity of CQP SAR images can be expressed as S l r = S r l . As described in Section 2.2, Cir 1 and Cir 4 are statistically derived from CQP SAR data (either undistorted CQP SAR data with reciprocity and reflection symmetry or simulated reciprocal CQP SAR with reciprocal crosstalk and reflection symmetry), to determine important reflection symmetry characteristics in volume-dominated targets. Next, we hypothesize that volume-dominated targets satisfy low Cir 1 and Cir 4 power, which is termed the semireflection symmetry of the volume-dominated targets.

3.1.2. Antenna Reciprocity, No Azimuth Drifting, and Range Drifting

Here, reasonable assumptions are made for CQP SAR systems. Firstly, the antennas are passive and can therefore be assumed to behave reciprocally: that is, the crosstalk during transmission and reception is equal [20]. We can derive u = z and w = v by applying the antenna reciprocity assumption to the calibration model (1).
Secondly, the system distortion changes with the flight time, because the system hardware is unstable. In addition, the unknown antenna distortions are reflected in the images at different incident angles, [20,32,33]. Considering that the azimuth time of each CQP SAR image is relatively short, we assume here that the system distortion parameters are constant along the azimuth direction: namely, no azimuth drifting. The distortion parameters are considered to be unstable along the range direction: namely, range drifting. According to these criteria, the uncalibrated CQP SAR images can be divided into bins along the azimuth and range directions, and the distortion parameters are consistent in different azimuth bins along each range direction in the following process.

3.2. Distributed Targets Selection for PolCAL

Before using PolCAL to determine the channel imbalances, it is important to extract distributed targets with appropriate scattering characteristics, so as to determine the accuracy of PolCAL. One criterion for extracting distributed targets is that ground targets with features that are less affected by additive noise in the corresponding calibration channels should be selected. By analyzing the process of transforming from the LQP SAR to the CQP SAR (13), we find that two co-pol channels, S l l and S r r , mainly exhibit volume-scattering features, while the x-pol channels S l r and S r l mainly show surface-scattering features. In natural scattering mechanisms, the volume-dominated target regions are mainly forest areas and complex artificial areas with specific incident angles, while the surface-dominated target regions mainly include bare soil and airport runways. However, identifying different incidence artificial areas is a time-consuming task, which imposes a computational burden on calibration with distributed targets: to address this issue, we select natural targets with strong volume scattering power and strong surface scattering power, after removing complex artificial targets through polarimetric parameters.
It should be pointed out that selecting appropriate distributed targets is challenging, because the scattering properties of the uncalibrated CQP SAR images deviate significantly from the polarization features of the actual undistorted data; therefore, after analyzing many extraction methods in LQP SAR images (such as [24,34,35]), ENL and R v b are introduced, to select the volume-dominated and surface-dominated targets. In general, compared with natural areas, the estimated ENL of complex artificial areas is smaller than the true value [24]; therefore, ENL is selected, to remove the effect of artificial areas. R v b —namely,
R v b = M 44 M 11 M 23 ,
which is related to the ratio of co-pol channel power and x-pol channel power—is used in the remaining regions, to select potential volume-dominated and surface-dominated targets. When we neglect the influences of the second-order items of crosstalk values, and the product of additive noise and crosstalk, R v b is virtually unaffected by distortion factors: namely, M 44 M 11 M 23 C 44 C 11 C 23 . Larger R v b values are used to select volume-dominated targets, while smaller R v b values are used to select surface-dominated targets. Based on the above description, the areas dominated by volume scattering and surface scattering are filtered according to the following criteria:
E N L t h e R v b < t h r l surface dominated R v b > t h r u volume dominated
After setting the ENL threshold, t h e , to discard urban points in the CQP SAR images, the R v b values in the remaining areas are sorted from small to large, and the index of the surface-dominated points ranges from 1 to the total number of points multiplied by t h r l ; the R v b values in the remaining areas are sorted from large to small, and the index of the volume-dominated points ranges from 1 to the total number of points multiplied by t h r u .
As previously mentioned, it is rather difficult to select scatter-dominated areas in uncalibrated CQP SAR images. In Section 4 and Section 5, we demonstrate that the proposed selection strategy (17) is an effective method for selecting volume-dominated and surface-dominated targets to estimate channel imbalances.

3.3. Channel Imbalances Calibration by Surface-Dominated and Volume-Dominated Targets

In the following section, the channel imbalances are calibrated with distributed targets, based on the selected volume-dominated and surface-scattered areas, and the assumptions about natural targets and CQP SAR systems.
Step 1: Considering that reciprocal crosstalk and k will not affect the calibration accuracy of α in reciprocal targets, (6) and (7), in [20], are utilized to calibrate α in the reciprocal surface-dominated targets. In the calibration process, a situation may arise in which there is no selected area along a range direction, resulting in an inability to determine the distortion parameters.
Here, we introduce the filtering operations on the amplitude and phase of α , to extend the amplitude and phase of the finite α values to the entire range of cells. Meanwhile, the filter operations can eliminate the errors caused by the selection of improper scattering areas, to a large extent, to obtain robust results.
The filtering operations can be divided into three steps: (1) we first obtain the α solution corresponding to the extracted surface-dominated areas in the uncalibrated CQP SAR images; (2) considering that the phases and amplitudes of α values obtained at different azimuths along a range direction are approximately the same, according to the no azimuth drift property, we perform first-order polynomial fitting on the azimuths with these values, and the fitting results are averaged to determine the amplitude and phase of α in this range direction; (3) the bin index (horizontal axis)-phase/amplitude of the α (vertical axis) coordinate system is rotated, to determine the main gradient direction, and the histogram statistics of the rotated data are obtained. A total of 70% of the points near the highest value in the histogram are selected for first-order polynomial fitting, to derive the continuous phases and amplitudes of α in the whole CQP SAR image.
For each bin, the corresponding filtered α values are used to calibrate the covariance matrices of all volume-dominated pixels, as follows:
M 4 L R = Q 1 M 4 L R Q 1 = K 4 C 4 L R K 4
Step 2: Semireflection symmetry is used to estimate k by the following process. By expanding (18), Cir 1 and Cir 4 can be expressed as
Cir 1 : p 4 M 11 M 44 = 0
Cir 4 : p 2 p M 12 + p M 24 = 0
where p is the inverse of k. The amplitude of p can be easily obtained from (19):
p = M 44 M 11 4
Compared to the amplitude of p, the phase of p is more difficult to determine accurately. Next, we propose an algorithm based on (20), to estimate the phase of p. We find that p = 0 is a constant solution of (20), which is inconsistent with the actual situation. Therefore, we divide p by (20), to obtain UCir 4 :
UCir 4 : p p M 12 + p p M 24 = 0
Next, Newton’s iterative method is utilized, to solve (22). Assuming that p = a + b j , where a and b are real numbers, the associated differentials of UCir 4 are
a UCir 4 = a b M 24 a 2 + b 2 3 2 + a 2 + b 2 M 12 + a a M 12 + b M 12 a 2 + b 2 a 2 M 24 a 2 + b 2 3 2 + M 24 a 2 + b 2 b UCir 4 = a b M 24 a 2 + b 2 3 2 + a 2 + b 2 M 12 + b a M 12 + b M 12 a 2 + b 2 b 2 M 24 a 2 + b 2 3 2 + M 24 a 2 + b 2
where · denotes the real operation of a complex number. Considering that at least two sets of data are needed to determine the two unknowns in (23), several azimuth bins at the same range are used, to obtain unique p values in the corresponding range bin, according to the azimuth non-drifting property. Assuming that the selected volume-dominated areas are in D D 2 azimuth bins, in the d iteration, the relationship between p and the changing values of a and b—namely, Δ a and Δ b —can be expressed as
p d = p d 1 + Δ a + j Δ b
where
Δ a , Δ b T = UCir 4 T UCir 4 1 UCir 4 T UCir 4
UCir 4 = a UCir 4 1 a UCir 4 2 a U C i r 4 D b UCir 4 1 b UCir 4 2 b U C i r 4 D T
UCir 4 = UCir 4 1 UCir 4 2 UCir 4 D T
The initial p values are set over a wide range, to minimize the effect on the Newton iteration method. The amplitudes are set at −3 dB∼3 dB, and the phases are set at −180 ° ∼180 ° . The final p value is selected according to the minimum residual error, which can be calculated by
error r e = UCir 4 Δ a , Δ b T UCir 4 2
where · 2 denotes the 2-norm of the matrix. Note that p and p e ± j π are the solutions of (22), and further processing is required to eliminate ambiguous values.
Combining the above solutions, based on Cir 1 and UCir 4 , (19) can be used to acquire only the amplitude of p, while (20) can be applied to determine the amplitude and phase of p, but the final phases have ± π ambiguity. Next, we define the filtering operations of k, which are roughly the same as the magnitude and phase filtering of α . The differences can be summarized as follows: (1) in the first amplitude filtering step, the initial value is the amplitude obtained by combining two amplitude solutions of Cir 1 and UCir 4 ; (2) in the phase filtering operation of k, the difference lies in the third step: after deriving the histogram statistics of the rotated coordinates, according to the main gradient, the histogram can be divided into several parts, because of the ± π phase shift. In each part, 80% of the points near the peak are selected to perform first-order polynomial fitting. The fitting results of the ± π shift phases and amplitudes are applied to calibrate k, and the correlation between the orientation angles obtained by the polarization data and the digital elevation model (DEM) decreases [22]. Therefore, the correct phase can be obtained according to the correlation. If the DEM cannot be obtained in the selected area, only the angle fitted by the data in the middle part of the histogram is selected as the final solution.
Finally, the obtained solutions of the phase and amplitude of k are introduced into the calibration matrix through (18), to obtain the final calibrated CQP SAR images.

4. Experiments and Results

In September 2021, AIRCAS carried out L-band airborne LQP SAR flight experiments at Yudaokou Airport, Hebei Province, China. The main SAR system characters are shown in Table 1. The typical features of the ground targets included farmland, forests, weedy fields, and airport runways. The original data size was 2000 × 2000 pixels in the azimuth and range directions, and the resolution was 0.659 × 0.681 m. Several kinds of CRs (see Figure 5) were deployed, to perform different scientific tasks, including trihedral CRs, 0 ° dihedral CRs, and 22.5 ° rotated dihedral CRs. In this section, the L-band CQP SAR data, which are synthesized according to the reciprocal LQP SAR data, are used to estimate the channel imbalances. We add different reciprocal crosstalk and channel imbalances to the reciprocal CQP SAR data, to simulate uncalibrated CQP SAR data. A 7 × 7 averaging window is used, to estimate the covariance matrix. Because the additional noise levels of different SAR systems vary, the following experiments are divided into two parts, to validate the robustness and accuracy of the proposed scheme: (1) channel imbalance calibration without additionally additive noise, and (2) channel imbalance calibration with additionally additive noise.

4.1. Channel Imbalances Calibration without Additionally Additive Noise

In this experiment, we impose co-pol and x-pol channel imbalances, with linear amplitudes ranging from −3 dB to 3 dB, and linear phases ranging from −180 ° to 180 ° , as well as crosstalk, with random amplitudes in the range −50 dB∼−15 dB, and random phases in the range −180 ° ∼180 ° . The crosstalk amplitude is set over a wide range, to verify the robustness of the proposed framework at different levels of crosstalk.
Figure 6a shows the true Pauli RGB image of the CQP SAR data, where the red channel is 1 2 S l l S r r 2 , the green channel is 1 2 S l r + S r l 2 , and the blue channel is 1 2 S l l + S r r 2 . The areas framed by the red line in Figure 6a are the runway, weedy fields, and steel fence of Yudaokou Airport. Moreover, we placed several CRs on the weedy fields, for different experimental purposes, which are shown in Figure 6b. TR.1–TR.5, DR.1–DR.2, and 225DR.1 indicate the trihedral CRs 1–5, the 0 ° dihedral CRs, and the 22.5 ° dihedral CR, respectively. The detailed latitudes, longitudes, and sizes of the CRs are shown in Table 2 and Figure 5a,b According to the flight plan, DR.1–DR.2 and 225DR.1 were prepared for this image. However, the altitude and heading of the route were slightly modified when the airplane took off, which caused the backscatter to deviate slightly from the theoretical values of DR.1–DR.2 and 225DR.1:
S D R l r = U l r U l r T 1 0 0 1 T = 1 0 0 1 T
S 225 DR l r = U l r U l r T 1 1 1 1 T = 1 j 0 0 1 j T
Next, we estimate the channel imbalances in the uncalibrated RGB CQP SAR image shown in Figure 6c, which was derived after adding the above crosstalk and channel imbalances. Upon adding crosstalk and channel imbalance, it becomes apparent that the Pauli RGB plot of the simulated uncalibrated image in Figure 6c exhibits significant color variation, when compared to the true Pauli RGB plot of the undistorted image depicted in Figure 6a. This highlights the necessity of calibrating distorted images, in order to obtain precise feature information for practical applications. Figure 6d shows the ENL images derived from Figure 6c. Heterogeneous scenes, such as steel fences, have small ENL values, while homogeneous areas, such as airport runways and bare soil, have larger ENL values; therefore, complex building areas can be removed after setting the ENL threshold, and the remaining parts are essentially homogeneous areas. R v b is shown in Figure 6e. The R v b values of the 0 ° dihedral CRs and 22.5 ° dihedral CR are larger than those of the remaining ground objects. The numerical relationships among the R v b values shown in Figure 6e can be deduced directly from (29) and (30). To more intuitively reflect the relative relationship among different ground features, in Figure 6f, the color scale is set to 0∼7, and there are slightly bright stripes (in the red box) along the range direction: this is because larger crosstalk may impact R v b . However, the relative relationship between the volume-dominated and surface-dominated R v b values does not change when the area with larger R v b is dominated by volume scatterers and the area with smaller R v b is dominated by surface scatterers.
Therefore, it is effective to select the area where surface scattering and volume scattering dominate according to the R v b histogram. Moreover, the amplitude and phase filtering operations proposed in this paper are performed when the channel imbalance is initially estimated, which reduces the error caused by selecting the wrong area. Here, t h e is set as 0.3, and t h r l and t h r u are set as 0.25 and 0.3, respectively. We note that there is no effective backscatter except noise on the airport runway, so the area selected from the runway cannot be used as effective candidate pixels. When using this selection algorithm to solve the problem, noise-dominated areas, such as calm water and cement floors, should be removed as much as possible.
After setting the number of range and azimuth bins to 100 and 25, respectively, the α and k estimation results are shown in Figure 7. Figure 7a,b represent the preliminary results of the amplitude and phase of α in different azimuth and range bins. By projecting the solution results along the Y axis onto the X-Z plane, we determine that the results of different azimuth bins in the same range bin are approximately the same. According to the α filtering operation shown in Figure 7c,d, the fitting results reduce the error of the initial estimation. After the simulation data are calibrated according to the fitting result of α , k is estimated according to the semireflection symmetry, and the results are shown in Figure 7e–h.
The fitting results shown in Figure 7g,h indicate that by applying the k filtering operations proposed in this paper, the initially estimated errors shown in Figure 7e,f, due to improper selection of the volume-dominated targets and the presence of low crosstalk, can be reduced. As mentioned above, Figure 7f,h show two ambiguous phases, which can be eliminated by using the DEM. If there is no corresponding DEM for the selected area, the middle-fitting line in the obtained result is selected as the final result. In Figure 7e,f, the control group (CG) results, using (8) and (11) to estimate the amplitude and phase of k, are shown in blue lines. Compared to the added channel imbalances, the residual errors of estimating the channel imbalances can be evaluated as
error amp _ d B = mean d B LSL amp k / α / True amp k / α error pha = mean LSL pha k / α True pha k / α
error k / α = mean LSL k / α True k / α / True k / α
where LSL · and True · represent the fitted and true values, respectively.
Compared to the added values (‘True’ in Figure 7c,d) of α , the amplitude and phase errors of α with the proposed scheme are 2.4516 × 10 4 dB/0.4865 ° , respectively, according to (31); according to (32), error α is 0.0085. According to (31), comparing the estimated k of the control group (‘ Cir 1 ’ and ‘ Cir 2 ’ in Figure 7g,h) with the added values (‘True’ in Figure 7g,h), the amplitude and phase errors of k are 0.1783 dB and 12.5263 ° , respectively, and the values determined by the algorithm proposed in this paper (‘K-FO’ in Figure 7g,h) are 0.1458 dB and 2.7558 ° , respectively. In addition, according to (32), the error k of the control group is 0.2215, and that of the algorithm proposed in this paper is 0.0515. It should be stressed that the k errors estimated by the proposed scheme and those of the control group may deviate slightly from the true values, because the originally synthesized CQP SAR data may have residual k. Next, the estimated accuracy is further illustrated by point targets, which is the most accurate error estimation method.
In Figure 8, the co-pol and x-pol signatures of the strong point targets are plotted on a linear polarimetric basis. The original signature diagrams are presented in rows 1 and 5. After adding the channel imbalances and crosstalk, the polarization signatures are distorted, as shown in rows 2 and 6. The results of the proposed scheme and the control group are presented in rows 3 and 7, and in rows 4 and 8, respectively. Compared to rows 1 and 5, the control group results have apparent deviations (marked by the dark red circles in Figure 8), while the results of the proposed scheme do not.
A quantitative evaluation of DR.1–DR.2 and 225DR.1 is presented in Table 3. Columns 1–2, 3–4, and 5–6 show the amplitude ratio (dB) and phase difference (degrees) between the LL and RR in the original synthesized CQP SAR images, the CQP SAR images calibrated by the proposed scheme, and the CQP SAR images calibrated by the control group, respectively. Due to the change in the flight route after deploying the CRs, the values presented in columns 1–2 are slightly different from the theoretical values. According to the evaluation system shown in (31), the mean residual errors of the three CRs (DR.1–DR.2 and 225DR.1) are 0.164 dB/2.590 ° , with the proposed scheme, and 0.292 dB/25.369 ° for the control group. Therefore, the proposed scheme is a more feasible calibration choice than the control group.

4.2. Channel Imbalances Calibration with Additionally Additive Noise

Considering that the additional noise in the SAR images varies, we append random additional noise to the simulated uncalibrated CQP SAR images used in Section 4.1 to validate the accuracy and robustness of the proposed scheme according to (2). The ratios of the additional noise amplitudes in each channel to the span are within −60 dB∼−25 dB, and the noise phases are within −180 ° ∼180 ° . In Figure 9, the co-pol and x-pol signatures of the CRs are plotted, where the meaning of each row is the same as in Figure 8. Although the additional noise affects the CR signatures, the results of the proposed scheme are still similar to the original results, and perform better than the control group. A quantitative evaluation of DR.1–DR.2 and 225DR.1, with additional noise, is presented in Table 4.
Overall, the additional noise impacts the statistical properties of the distributed targets, resulting in larger estimated errors than when additional noise is not included, in most cases. As shown in Table 2, the average amplitude and phase errors of the proposed scheme and control group are 0.614 dB/4.154 ° and 0.927 dB/23.893 ° , respectively. Combined with the results in Section 4.1, these results show that the proposed scheme is a better calibration choice for channel imbalances with distributed targets in CQP SAR images with reciprocal crosstalk.

5. Discussion

The above algorithm includes some factors, such as threshold selection, that may affect accuracy. In addition, the effect of crosstalk on estimating channel imbalances should be further discussed, in order to design SAR systems. The detailed discussions of these problems are presented in this section.

5.1. Number of Range and Azimuth Bins

In this article, the CQP SAR images are divided into different bins along the azimuth and range directions, when estimating the channel imbalances. The number of range bins is strongly related to the drift in the range direction in CQP SAR systems. When the variation in the channel imbalance along the range direction is stable, less range bins are needed. The number of azimuth bins is highly dependent on the Newton iteration calculation process in (22). Because the unknowns in (22) are the real and imaginary parts of p, the number of azimuth bins, D, is set to greater than or equal to 2. When D is set to be larger than 2, the iterative process may be more accurate, but the computational efficiency is reduced. If computational efficiency is the main priority, the readers may prefer to try a smaller D, and compare the corresponding estimated results with the CR results, to determine the value that meets the design specification, before calibrating the other images. If CRs are unavailable, it is recommended to set a larger D, to ensure high estimation accuracy. It is also possible to combine the results estimated under multiple bins, and to filter the results to obtain high estimation accuracy [22].

5.2. Choice of ENL Threshold

When extracting surface-dominated and volume-dominated targets before PolCAL, the ENL is utilized to remove heterogeneous scattering areas, such as cities. The heterogeneous scattering areas are removed for the following reasons: (1) when using R v b to extract volume-dominated areas, targets dominated by secondary scattering are also selected from heterogeneous areas, such as cities, resulting in reduced calibration accuracy when the semireflection symmetry is applied; if the ENL is used to remove as many secondary scattering areas as possible, the accuracy of the selected areas is effectively increased, increasing the accuracy of the final calibration result to some extent; (2) in practical applications, we usually use the multilook operation to estimate the covariance matrix; at the boundary between heterogeneous and homogeneous targets, the multilook operation may make the R v b values of the extracted features differ from those of the required features; as a result, heterogeneous scattering features may be extracted, decreasing the PolCAL accuracy. For example, mixing urban targets and dark targets increases R v b , as the targets may be misinterpreted as volume scattering targets. The heterogeneous property leads to a small ENL for complex scattering areas, such as cities, and a high ENL for naturally homogeneous areas; therefore, urban areas can be eliminated by setting an ENL threshold. Here, we set the ENL threshold to 0.3, to remove complex scattering areas, and select more natural areas dominated by surface scattering and volume scattering.

5.3. Choice of R v b Threshold

After removing complex scattering regions, such as cities, the selection of appropriate scattering properties for channel imbalance calibration remains a great challenge, as many polarization parameters vary with system distortion. In addition to the time-consuming and laborious manual selection, in this paper, R v b is considered an available option to extract surface-dominated and volume-dominated targets in uncalibrated CQP SAR images. In Figure 10, we show the scatter characteristics of our extracted targets in the H / α plane [36], where t h e , t h r u , and t h r l are set to 0.3, 0.3, and 0.25, respectively.
In Figure 10, most surface-dominated points are found at low-entropy, medium-entropy, and low-alpha angles, while volume-dominated points are found at medium-entropy, high-entropy, low-alpha, and medium-alpha angles. However, the selected area, which should be dominated by volume scattering, also contains a smaller number of surface-dominated targets. This is because most of the polarization parameters cannot accurately represent the characteristics of the selected area, due to system distortion. In fact, there are very few algorithms for selecting points with absolutely correct scattering properties in uncalibrated CQP SAR data. It should be stressed that the accuracy in this paper depends not only on the quantity and quality of the selected regions, but also on the channel imbalance filtering operations. After the channel imbalances are initially estimated, the results obtained, based on the inaccurate scattering areas, are filtered out as much as possible, thereby improving the final calibration results.
In Figure 11, we evaluate the errors of estimating the channel imbalances based on (31) by setting different R v b thresholds. In the following section, we set ENL to 0.3, while t h r l and t h r u are uniformly changed from 0.2 to 0.8 in 0.01 steps.
As previously mentioned, in the proposed scheme, unknown reciprocal crosstalk does not significantly impact the estimation of α when the CQP SAR image is reciprocal. In Figure 11a,b, when different t h r l are set to select the surface-dominated targets, the amplitude and phase errors of α are relatively small, within 1.5 × 10 3 dB/0.8 ° . Even if the volume-dominated targets are selected to estimate α , the calibration accuracy of α is ensured by the reciprocity of the volume-dominated targets. In practical PolCAL, when t h r l is set to a larger value, the computational time to estimate α increases, which may not be acceptable in multiple image calibration. Therefore, a lower t h r l value, such as 0.25, is recommended to improve the calibration efficiency when ensuring the calibration accuracy.
The estimated errors of k are presented in Figure 11c,d. In the amplitude errors of Figure 11c, when the wrong semireflection symmetry area is selected, by setting the inappropriate thresholds, the fitting operations will increase the instability of the estimated results, which will lead to a worse estimated result than other situations. However, most of the amplitude errors are stable, below 0.4 dB.
In Figure 11d, as t h r u increases, the estimated phase errors gradually increase up to a certain point, and then start to decrease. The same situation is shown in the reset maximum thumbnail of Figure 11c. As t h r u increases, the energy of the selected targets in the co-pol channel is smaller than that in the volume-dominated area, which causes noise to affect the semireflection symmetry of the selected area, and increases the estimation errors: this is because, as an increasing number of surface-dominated targets are selected that satisfy the semireflection symmetry, the solution errors decrease according to the filter operation. We note that when more semireflective symmetry targets are applied to estimate k, the noise effect on the accuracy calibration of k is attenuated, although noise in the image still affects the semireflection symmetry in the surface-dominated regions. There is a combined area between the surface-dominated and volume-dominated regions, and the reflection symmetry of this area was not analyzed in detail, which seriously affects the accuracy of k, as shown in Figure 11c; therefore, we select only volume-scattering-dominant regions, and 0.3 is recommended as the optimal value of t h r u .

5.4. Crosstalk Effect on Channel Imbalance Estimation

In Section 4, the simulated uncalibrated CQP SAR data are used to verify the capability of the proposed scheme, when the reciprocal crosstalk along the range direction is randomly set within a wide range of −50 dB∼−15 dB, to verify the adaptability of different crosstalk systems: namely, −50 dB u = z , v = w 15 dB and 180 ° u = z , v = w 180 ° . When researchers design CQP SAR systems, the crosstalk (i.e., the power leakage between channels) between the antennas is usually verified by a fixed value, for example, better than −35 dB, better than −30 dB, etc; therefore, in this part, different levels of fixed crosstalk are added to the undistorted CQP SAR images, to estimate the channel imbalance, and to calculate the determination error, with the aim of validating the proposed framework as an effective tool for system design.
In each trial, we impose: (1) system crosstalk with constant amplitudes ranging from −50 dB to −15 dB (−50 dB u = z = v = w 15 dB) and random phases ( 180 ° u = z = v = w 180 ° ); (2) channel imbalances with linear amplitudes ranging from −3 dB to 3 dB and linear phases ranging from −180 ° to 180 ° ; t h e , t h r u and t h r l are set to 0.3, 0.3, and 0.25, respectively. To show the effect of the crosstalk phase on the channel imbalance estimation results, 45 Monte Carlo trials were conducted for each crosstalk amplitude. The results are shown in the boxplots in Figure 12.
As has been emphasized throughout the article, reciprocal crosstalk has little negative effect on the α estimation accuracy. In Figure 11a,b, as the crosstalk increases, the α estimation error remains low, within 1.4 × 10 3 dB and 0.7 ° . However, the crosstalk affects R v b , the selection of the volume-scattering-dominant region, and the establishment of (22). Thus, in Figure 12c,d, different crosstalk unknowns have more pronounced impacts on estimating k than α . As the crosstalk amplitude increases, the accuracy and stability of estimating k decreases significantly. When the crosstalk amplitude is greater than −18 dB, the corresponding estimation accuracy drops directly from less than 0.4 dB to more than 0.5 dB. In addition, the stability effect of the crosstalk phase on the amplitude estimation becomes obvious, and the cabinet length of −17 dB is nearly twice that of −18 dB. If 0.5 dB/5 ° (such as GF-3) is used as the channel imbalance design index, then this algorithm should be used with crosstalk greater than −23 dB.

6. Conclusions

In the existing distributed target calibration for channel imbalances of CQP SAR systems with reciprocal crosstalk, some reflection symmetry assumptions of volume-dominated targets are used to estimate the phase of the co-pol channel imbalance k, i.e., (9)–(11). This paper notes that the reflection symmetry used to estimate the phase of k has poor applicability in volume-dominated target regions, which reduces the estimation accuracy of the k phase. In this paper, we assess the applicability of the reflection symmetry in volume-dominated target regions, and propose an algorithm for k phase calibration based on high-priority reflection symmetry (we name it semireflection symmetry) in volume-dominated targets. Furthermore, as the reciprocity of the surface-dominated targets is sufficient to estimate the x-pol channel imbalance α with [20], a channel imbalance calibration scheme based on crosstalk reciprocity in SAR systems is proposed for polarimetric applications, in which crosstalk can be ignored and CRs are not applied to calibrate CQP SAR images. The framework includes the automatic extraction of surface-dominated and volume-dominated targets as reference pixels, the estimation and filtering of α in the surface-dominated target regions, and the estimation and filtering of k in the volume-dominated target regions. The filtering operations aim to obtain more robust results. We conducted experiments based on CQP SAR images synthesized from reciprocal L-band LQP SAR images. A high calibration accuracy of α was achieved by the experimental part, which was consistent with [20]. For the k calibration, compared to the control group, the algorithm proposed in this paper has an improvement of 0.1 dB/22 ° by CRs without additionally additive noise. The improvement can reach 0.3 dB/19 ° with additionally additive noise. It is recommended that t h e , t h r u and t h r l are set to 0.3, 0.3, and 0.25, respectively, to select the appropriate surface-dominated and volume-dominated targets. In addition, we recommend that researchers use this method when the crosstalk design specification of the system is better than −23 dB. In future work, we will implement the scheme in more measured CQP SAR images, to improve the method.

Author Contributions

Conceptualization, Y.D. and X.L.; methodology, X.Z.; software, X.Z.; validation, X.Z.; formal analysis, H.Z.; investigation, X.Z. and H.Z.; resources, H.Z., Y.D. and X.L.; writing—original draft preparation, X.Z.; writing—review and editing, X.Z., H.Z., Y.D. and X.L.; project administration, Y.D. and X.L.; funding acquisition, Y.D. and X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Beijing Municipal Natural Science Foundation (grant no. 4192065) and the National Natural Science Foundation of China (grant no. 61901445).

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Department of Space Microwave Remote Sensing System, Aerospace Information Research Institute, Chinese Academy of Sciences for providing L-band SAR data. The authors would also like to thank all colleagues who participated in this experiment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, T.; Wang, W.; Yang, Z.; Yin, J.; Yang, J. Ship Detection From PolSAR Imagery Using the Hybrid Polarimetric Covariance Matrix. IEEE Geosci. Remote Sens. Lett. 2021, 18, 1575–1579. [Google Scholar] [CrossRef]
  2. Ainsworth, T.; Kelly, J.; Lee, J.S. Classification comparisons between dual-pol, compact polarimetric and quad-pol SAR imagery. ISPRS J. Photogramm. Remote Sens. 2009, 64, 464–471. [Google Scholar] [CrossRef]
  3. Duysak, H.; Yiğit, E. Investigation of the performance of different wavelet-based fusions of SAR and optical images using Sentinel-1 and Sentinel-2 datasets. Int. J. Eng. Geosci. 2022, 7, 81–90. [Google Scholar] [CrossRef]
  4. Inderkumar, K.; Maurya, H.; Kumar, A.; Bhiravarasu, S.S.; Das, A.; Putrevu, D.; Pandey, D.K.; Panigrahi, R.K. Retrieval of lunar surface dielectric constant using Chandrayaan-2 full-polarimetric SAR data. IEEE Trans. Geosci. Remote. Sens. 2022, 60, 1–12. [Google Scholar] [CrossRef]
  5. Ni, J.; López-Martínez, C.; Hu, Z.; Zhang, F. Multitemporal SAR and Polarimetric SAR Optimization and Classification: Reinterpreting Temporal Coherence. IEEE Trans. Geosci. Remote. Sens. 2022, 60, 1–17. [Google Scholar] [CrossRef]
  6. Park, S.E.; Jung, Y.T.; Kim, H.C. Monitoring permafrost changes in central Yakutia using optical and polarimetric SAR data. Remote. Sens. Environ. 2022, 274, 112989. [Google Scholar] [CrossRef]
  7. Wang, H.; Yang, J.; Lin, M.; Li, W.; Zhu, J.; Ren, L.; Cui, L. Quad-polarimetric SAR sea state retrieval algorithm from Chinese Gaofen-3 wave mode imagettes via deep learning. Remote. Sens. Environ. 2022, 273, 112969. [Google Scholar] [CrossRef]
  8. 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]
  9. Ainsworth, T.L.; Ferro-Famil, L.; Lee, J.S. Orientation angle preserving a posteriori polarimetric SAR calibration. IEEE Trans. Geosci. Remote. Sens. 2006, 44, 994–1003. [Google Scholar] [CrossRef]
  10. Freeman, A.; Van Zyl, J.J.; Klein, J.D.; Zebker, H.A.; Shen, Y. Calibration of Stokes and scattering matrix format polarimetric SAR data. IEEE Trans. Geosci. Remote. Sens. 1992, 30, 531–539. [Google Scholar] [CrossRef]
  11. Freeman, A. SAR calibration: An overview. IEEE Trans. Geosci. Remote. Sens. 1992, 30, 1107–1121. [Google Scholar] [CrossRef]
  12. Quegan, S. A unified algorithm for phase and cross-talk calibration of polarimetric data-theory and observations. IEEE Trans. Geosci. Remote. Sens. 1994, 32, 89–99. [Google Scholar] [CrossRef]
  13. Han, Y.; Liu, X.; Hou, W.; Gao, Y.; Wang, R.; Fan, H.; Liu, D.; Zhao, F. A Distributed Target-Based Calibration Method for Hybrid Quadrature-Polarimetric SAR. In Proceedings of the IGARSS 2022-2022 IEEE International Geoscience and Remote Sensing Symposium, Kuala Lumpur, Malaysia, 17–22 July 2022; pp. 2538–2541. [Google Scholar] [CrossRef]
  14. Kimura, H.; Mizuno, T.; Papathanassiou, K.P.; Hajnsek, I. Improvement of polarimetric SAR calibration based on the Quegan algorithm. In Proceedings of the IGARSS 2004. 2004 IEEE International Geoscience and Remote Sensing Symposium, Anchorage, AK, USA, 20–24 September 2004; Volume 1. [Google Scholar] [CrossRef]
  15. Zhang, H.; Lu, W.; Zhang, B.; Chen, J.; Wang, C. Improvement of polarimetric SAR calibration based on the Ainsworth algorithm for Chinese airborne PolSAR data. IEEE Geosci. Remote. Sens. Lett. 2013, 10, 898–902. [Google Scholar] [CrossRef]
  16. Azcueta, M.; d’Alessandro, M.M.; Zajc, T.; Grunfeld, N.; Thibeault, M. ALOS-2 preliminary calibration assessment. In Proceedings of the 2015 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Milan, Italy, 26–31 July2015; pp. 4117–4120. [Google Scholar] [CrossRef]
  17. Shangguan, S.; Qiu, X.; Fu, K.; Lei, B.; Hong, W. GF-3 polarimetric data quality assessment based on automatic extraction of distributed targets. IEEE J. Sel. Top. Appl. Earth Obs. Remote. Sens. 2020, 13, 4282–4294. [Google Scholar] [CrossRef]
  18. Touzi, R.; Hawkins, R.K.; Cote, S. High-precision assessment and calibration of polarimetric RADARSAT-2 SAR using transponder measurements. IEEE Trans. Geosci. Remote. Sens. 2012, 51, 487–503. [Google Scholar] [CrossRef]
  19. Lee, J.S.; Pottier, E. Polarimetric Radar Imaging: From Basics to Applications; CRC Press: New York, NY, USA, 2017. [Google Scholar]
  20. Pincus, P.; Preiss, M.; Goh, A.S.; Gray, D. Polarimetric calibration of circularly polarized synthetic aperture radar data. IEEE Trans. Geosci. Remote. Sens. 2017, 55, 6824–6839. [Google Scholar] [CrossRef]
  21. Cloude, S.R.; Ossikovski, R.; Garcia-Caurel, E. Bright singularities: Polarimetric calibration of spaceborne PolSAR systems. IEEE Geosci. Remote. Sens. Lett. 2020, 18, 476–479. [Google Scholar] [CrossRef]
  22. Shi, L.; Yang, J.; Li, P. Co-polarization channel imbalance determination by the use of bare soil. ISPRS J. Photogramm. Remote. Sens. 2014, 95, 53–67. [Google Scholar] [CrossRef]
  23. Van Zyl, J.J. Calibration of polarimetric radar images using only image parameters and trihedral corner reflector responses. IEEE Trans. Geosci. Remote. Sens. 1990, 28, 337–348. [Google Scholar] [CrossRef]
  24. Shi, L.; Li, P.; Yang, J.; Zhang, L.; Ding, X.; Zhao, L. Polarimetric SAR calibration and residual error estimation when corner reflectors are unavailable. IEEE Trans. Geosci. Remote. Sens. 2020, 58, 4454–4471. [Google Scholar] [CrossRef]
  25. Dey, S.; Bhattacharya, A.; Frery, A.C.; López-Martínez, C.; Rao, Y.S. A model-free four component scattering power decomposition for polarimetric SAR data. IEEE J. Sel. Top. Appl. Earth Obs. Remote. Sens. 2021, 14, 3887–3902. [Google Scholar] [CrossRef]
  26. Sun, G.; Li, Z.; Huang, L.; Chen, Q.; Zhang, P. Quality analysis and improvement of polarimetric synthetic aperture radar (SAR) images from the GaoFen-3 satellite using the Amazon rainforest as an example. Int. J. Remote. Sens. 2021, 42, 2131–2154. [Google Scholar] [CrossRef]
  27. Huynen, J.R. Phenomenological Theory of Radar Targets. Ph.D. Thesis, Technical University, Delft, The Netherlands, 1970. [Google Scholar]
  28. LÜNEBURG, E. Principles of radar polarimetry. Ieice Trans. Electron. 1995, 78, 1339–1345. [Google Scholar]
  29. Lim, Y.X.; Burgin, M.S.; van Zyl, J.J. An optimal nonnegative eigenvalue decomposition for the Freeman and Durden three-component scattering model. IEEE Trans. Geosci. Remote. Sens. 2017, 55, 2167–2176. [Google Scholar] [CrossRef]
  30. Freeman, A.; Durden, S.L. A three-component scattering model for polarimetric SAR data. IEEE Trans. Geosci. Remote. Sens. 1998, 36, 963–973. [Google Scholar] [CrossRef] [Green Version]
  31. Aubry, A.; Carotenuto, V.; De Maio, A.; Pallotta, L. Assessing reciprocity in polarimetric SAR data. IEEE Geosci. Remote. Sens. Lett. 2019, 17, 87–91. [Google Scholar] [CrossRef]
  32. Fore, A.G.; Chapman, B.D.; Hawkins, B.P.; Hensley, S.; Jones, C.E.; Michel, T.R.; Muellerschoen, R.J. UAVSAR polarimetric calibration. IEEE Trans. Geosci. Remote. Sens. 2015, 53, 3481–3491. [Google Scholar] [CrossRef]
  33. Shi, L.; Li, P.; Yang, J.; Zhang, L.; Ding, X.; Zhao, L. Co-polarization channel imbalance phase estimation by corner-reflector-like targets. ISPRS J. Photogramm. Remote. Sens. 2019, 147, 255–266. [Google Scholar] [CrossRef]
  34. Shi, L.; Li, P.; Yang, J.; Sun, H.; Zhao, L.; Zhang, L. Polarimetric calibration for the distributed Gaofen-3 product by an improved unitary zero helix framework. ISPRS J. Photogramm. Remote. Sens. 2020, 160, 229–243. [Google Scholar] [CrossRef]
  35. Chang, Y.; Li, P.; Yang, J.; Zhao, J.; Zhao, L.; Shi, L. Polarimetric calibration and quality assessment of the GF-3 satellite images. Sensors 2018, 18, 403. [Google Scholar] [CrossRef] [Green Version]
  36. Cloude, S.R.; Pottier, E. An entropy based classification scheme for land applications of polarimetric SAR. IEEE Trans. Geosci. Remote. Sens. 1997, 35, 68–78. [Google Scholar] [CrossRef]
Figure 1. Histograms of (a) Cir 1 / 2 Span , (b) ACir 2 , and (c) ACir 3 of volume-dominated targets in the nine CQP SAR images.
Figure 1. Histograms of (a) Cir 1 / 2 Span , (b) ACir 2 , and (c) ACir 3 of volume-dominated targets in the nine CQP SAR images.
Remotesensing 15 01365 g001
Figure 2. Histograms of Cir 1 / 2 Span (dotted lines) and Cir 3 / 2 Span (solid lines).
Figure 2. Histograms of Cir 1 / 2 Span (dotted lines) and Cir 3 / 2 Span (solid lines).
Remotesensing 15 01365 g002
Figure 3. Comparison of Cir 4 / 2 Span and Cir 5 / 2 Span : (a) histograms of Cir 4 / 2 Span and Cir 5 / 2 Span in the volume-dominated area of nine undistorted CQP SAR images; (b) percentages of all volume-dominated points with Cir 4 / 2 Span < 0.1 and Cir 5 / 2 Span < 0.1 in each undistorted CQP SAR image; (c,d) percentages of all reflection symmetry points with Cir 4 / 2 Span < 0.1 and Cir 5 / 2 Span < 0.1 in the presence of reciprocal crosstalk.
Figure 3. Comparison of Cir 4 / 2 Span and Cir 5 / 2 Span : (a) histograms of Cir 4 / 2 Span and Cir 5 / 2 Span in the volume-dominated area of nine undistorted CQP SAR images; (b) percentages of all volume-dominated points with Cir 4 / 2 Span < 0.1 and Cir 5 / 2 Span < 0.1 in each undistorted CQP SAR image; (c,d) percentages of all reflection symmetry points with Cir 4 / 2 Span < 0.1 and Cir 5 / 2 Span < 0.1 in the presence of reciprocal crosstalk.
Remotesensing 15 01365 g003
Figure 4. Channel imbalance calibration scheme for CQP SAR images with reciprocal crosstalk.
Figure 4. Channel imbalance calibration scheme for CQP SAR images with reciprocal crosstalk.
Remotesensing 15 01365 g004
Figure 5. Ground-deployed CRs at Yudaokou Airport in September 2021: (ac) trihedral CR, dihedral CR, and 22.5 ° rotated dihedral CR, respectively; (d,e) placement process and measuring equipment.
Figure 5. Ground-deployed CRs at Yudaokou Airport in September 2021: (ac) trihedral CR, dihedral CR, and 22.5 ° rotated dihedral CR, respectively; (d,e) placement process and measuring equipment.
Remotesensing 15 01365 g005
Figure 6. L-band CQP SAR images: (a) the true Pauli RGB CQP SAR image, where the red channel is 1 2 S l l S r r 2 , the green channel is 1 2 S l r + S r l 2 , and the blue channel is 1 2 S l l + S r r 2 ; the radar illuminated the scene from top to bottom; the airport location is marked in red; (b) enlarged view of the near range in (a); (c) simulated uncalibrated Pauli RGB image; (d) ENL image; (e,f) R v b images with different scales; the red box in (f) marks bright stripes.
Figure 6. L-band CQP SAR images: (a) the true Pauli RGB CQP SAR image, where the red channel is 1 2 S l l S r r 2 , the green channel is 1 2 S l r + S r l 2 , and the blue channel is 1 2 S l l + S r r 2 ; the radar illuminated the scene from top to bottom; the airport location is marked in red; (b) enlarged view of the near range in (a); (c) simulated uncalibrated Pauli RGB image; (d) ENL image; (e,f) R v b images with different scales; the red box in (f) marks bright stripes.
Remotesensing 15 01365 g006
Figure 7. Estimated PolCAL parameters of the simulated CQP SAR images: (a,b) initial amplitude and phase results of α ; (c,d) alpha filter amplitude and phase; alpha-FO indicates the α filtering results; (e,f) initial amplitude and phase results of k; (g,h) K filter amplitude and phase; K-FO denotes the k filtering results after using Cir 1 & UCir 4 / UCir 4 to estimate the initial amplitude/phase of k.
Figure 7. Estimated PolCAL parameters of the simulated CQP SAR images: (a,b) initial amplitude and phase results of α ; (c,d) alpha filter amplitude and phase; alpha-FO indicates the α filtering results; (e,f) initial amplitude and phase results of k; (g,h) K filter amplitude and phase; K-FO denotes the k filtering results after using Cir 1 & UCir 4 / UCir 4 to estimate the initial amplitude/phase of k.
Remotesensing 15 01365 g007
Figure 8. Polarimetric (co-pol and x-pol) signature responses of the strong point targets in the CQP SAR data set, without additive noise in linear polarimetric basis.
Figure 8. Polarimetric (co-pol and x-pol) signature responses of the strong point targets in the CQP SAR data set, without additive noise in linear polarimetric basis.
Remotesensing 15 01365 g008
Figure 9. Polarimetric (co-pol and x-pol) signature responses of the strong point targets in the CQP SAR data set, with additive noise on linear polarimetric basis.
Figure 9. Polarimetric (co-pol and x-pol) signature responses of the strong point targets in the CQP SAR data set, with additive noise on linear polarimetric basis.
Remotesensing 15 01365 g009
Figure 10. The (a) surface-dominated and (b) volume-dominated targets extracted from the uncalibrated CQP SAR data in Section 4.1. Different point colors represent different areas in the H / α plane.
Figure 10. The (a) surface-dominated and (b) volume-dominated targets extracted from the uncalibrated CQP SAR data in Section 4.1. Different point colors represent different areas in the H / α plane.
Remotesensing 15 01365 g010
Figure 11. Calibration accuracy of the proposed scheme under different R v b thresholds: (a,b) amplitude and phase errors of α ; (c,d) amplitude and phase errors of k.
Figure 11. Calibration accuracy of the proposed scheme under different R v b thresholds: (a,b) amplitude and phase errors of α ; (c,d) amplitude and phase errors of k.
Remotesensing 15 01365 g011
Figure 12. Effect of different crosstalk levels on the channel imbalance estimation: (a,b) amplitude and phase errors of α ; (c,d) amplitude and phase errors of k. The standard boxplot notation is used (lower/upper hinges—first/third quartiles; whiskers extend from the hinges to the largest/lowest values no further than 1.5 × interquartile ranges).
Figure 12. Effect of different crosstalk levels on the channel imbalance estimation: (a,b) amplitude and phase errors of α ; (c,d) amplitude and phase errors of k. The standard boxplot notation is used (lower/upper hinges—first/third quartiles; whiskers extend from the hinges to the largest/lowest values no further than 1.5 × interquartile ranges).
Remotesensing 15 01365 g012
Table 1. LQP SAR system technical specifications.
Table 1. LQP SAR system technical specifications.
Characters Value
Center Frequency (GHz) 1.26
Bandwidth (MHz) 200
Pulse Repeat Frequency (Hz) 1805
Flight Hight (m) 4645.64
Platform Velocity (m/s) 84.93
Incidence Angle (Degree) 55
Table 2. Summery of the CRs. N and E of columns 3 and 4 are north latitude and east longitude, respectively. Detailed schematic diagrams of the CR sizes can be referred to Figure 5a,b.
Table 2. Summery of the CRs. N and E of columns 3 and 4 are north latitude and east longitude, respectively. Detailed schematic diagrams of the CR sizes can be referred to Figure 5a,b.
CR TypesIDLongitude (Degree)Latitude (Degree)Size (m)
Trihedral CRTR.1116.869611E42.139556N1
TR.2116.870528E42.140233N1
TR.3116.872422E42.138997N1
TR.4116.873844E42.138783N1
TR.5116.874867E42.139358N1
0 ° Dihedral CRDR.1116.871078E42.139272N1
DR.2116.873450E42.139661N1
22.5 ° Dihedral CR225DR.1116.872117E42.140114N1
Table 3. Quality evaluation of DR.1–DR.2 and 225DR.1 in the CQP SAR data set without additionally additive noise. The unit of measure in the second row, the fourth row, and the fifth row is dBs, and the unit of measure in the third row is degrees ( ° ). For obvious comparison, the values in brackets in the seventh column are in phase periods of 0 to 360 degrees. The last row shows the errors of the results of the proposed method and the control group compared to the corner reflectors by (31).
Table 3. Quality evaluation of DR.1–DR.2 and 225DR.1 in the CQP SAR data set without additionally additive noise. The unit of measure in the second row, the fourth row, and the fifth row is dBs, and the unit of measure in the third row is degrees ( ° ). For obvious comparison, the values in brackets in the seventh column are in phase periods of 0 to 360 degrees. The last row shows the errors of the results of the proposed method and the control group compared to the corner reflectors by (31).
OriginalProposed SchemeControl Group
DR.10.923172.0290.902174.9460.890−168.993 (191.007)
DR.20.932172.2240.933174.5760.917−160.542 (199.458)
225DR.10.91190.1950.88292.6960.869120.089
Mean error
by (31)
-0.1642.5900.29225.369
Table 4. Quality evaluation of DR.1–DR.2 and 225DR.1 in CQP SAR data set with additionally additive noise. The unit of measure in the second row, the fourth row, and the fifth row is dBs, and the unit of measure in the third row is degrees ( ° ). For obvious comparison, the values in brackets in the seventh column are in phase periods of 0 to 360 degrees. The last row shows the errors of the results of the proposed method and the control group compared to the corner reflectors by (31).
Table 4. Quality evaluation of DR.1–DR.2 and 225DR.1 in CQP SAR data set with additionally additive noise. The unit of measure in the second row, the fourth row, and the fifth row is dBs, and the unit of measure in the third row is degrees ( ° ). For obvious comparison, the values in brackets in the seventh column are in phase periods of 0 to 360 degrees. The last row shows the errors of the results of the proposed method and the control group compared to the corner reflectors by (31).
OriginalProposed SchemeControl Group
DR.10.923172.0290.943178.9530.955−164.271 (195.727)
DR.20.932172.2240.983178.3511.000−158.684 (201.316)
225DR.10.91190.1950.99993.7551.014119.529
Mean error
by (31)
-0.6144.1540.92723.893
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

Zhao, X.; Deng, Y.; Zhang, H.; Liu, X. A Channel Imbalance Calibration Scheme with Distributed Targets for Circular Quad-Polarization SAR with Reciprocal Crosstalk. Remote Sens. 2023, 15, 1365. https://doi.org/10.3390/rs15051365

AMA Style

Zhao X, Deng Y, Zhang H, Liu X. A Channel Imbalance Calibration Scheme with Distributed Targets for Circular Quad-Polarization SAR with Reciprocal Crosstalk. Remote Sensing. 2023; 15(5):1365. https://doi.org/10.3390/rs15051365

Chicago/Turabian Style

Zhao, Xingjie, Yunkai Deng, Heng Zhang, and Xiuqing Liu. 2023. "A Channel Imbalance Calibration Scheme with Distributed Targets for Circular Quad-Polarization SAR with Reciprocal Crosstalk" Remote Sensing 15, no. 5: 1365. https://doi.org/10.3390/rs15051365

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