In this paper, the effectiveness and feasibility of the proposed method were verified through the following series of experiments. The robustness of the proposed method was verified by applying it to various data sources, especially on P/L/S/C/X multi-band data of airborne data. For the 180-degree phase ambiguity problem, two implementation examples realized by specific ground buildings are given in
Section 3.2 and
Section 3.3. Through these two complete solution examples, the correctness of the method’s principle was verified, and a deeper understanding of the double-bounce reflections of the ground surface was obtained. Since the method can be applied to estimate the phase imbalance of full-polarization data, we used other polarization calibration algorithms to verify the accuracy of the estimated phase imbalances.
3.1. Phase Imbalance Monitoring Using Sentinel-1 Dual-Pol Data
As seen in the flow chart of
Figure 3, the 180-degree phase ambiguity determination operation can be omitted to obtain a fast phase imbalance monitoring for the calibrated SAR images. Because the residual phase imbalance is unlikely to reach or exceed 90°, the peak close to 0° should be the target peak of the two estimated phase-imbalance peaks. The experiment used Sentinel-1 calibrated dual-polarization single look complex (SLC) data of the Stripmap mode. The Stripmap mode acquires data with an 80 km swath and 5 m by 5 m spatial resolution. The characteristics of the Stripmap products are a phase error within 5°, radiometric accuracy with 1 dB (3σ), and maximum NESZ with −22 dB [
23]. The operation of the experiment data is the same as shown in
Figure 3, that is, first, extract the RDBs through polarimetric strong scattering points filtering and coherence filtering, then, calculate each RDB’s phase imbalance. Finally, the histogram statistics of the results are performed. In the specific parameters in the experiment, K was 2, and the threshold was set to 0.8 for local 7-by-7 pixel areas. The results based on Sentinel-1 data are shown in
Figure 4.
This scene locates in the city of Houston and contains the entire urban area. One local area is shown in
Figure 4a, where the yellow crosses mark the extracted RDB targets. There are a total of 186,474 RDBs extracted in this scenario. The histogram statistical results are shown in
Figure 4b. As can be seen, the peak of the distribution histogram occurs at the position of about 0°, which is in line with expectation, and the other peak that locates around 180° can be ignored. Meanwhile,
Figure 4b shows the normal fitting curve of the data within the range of plus or minus 40° centered at 0°, and the corresponding peak position can be obtained as 2.55° according to the normal fitting results. Therefore, the phase imbalance of this scene is estimated to be 2.55°.
As seen from the experimental results, there were enough RDB targets extracted in such an urban located scene, and its statistical distribution presented an obvious characteristic similar to the normal distribution, which can ensure that the final results obtained are accurate enough in the statistical sense. The results obtained in the experiments prove an excellent polarization phase performance of the Sentinel-1 data, which is consistent with its technical index that phase error within 5°. At the same time, the experimental results also prove the effectiveness of the proposed method for phase imbalance monitoring applications. Although the quantity of RDBs is large, the processing algorithm of individual RDB is not complicated, so the method can realize the phase imbalance monitoring mission quickly.
3.2. Phase Imbalance Estimation Using Uncalibrated GF-3 Quad-Pol Data
In this part, we use GAOFEN-3 (GF-3) quad-pol data for the estimation of the unknown polarimetric phase imbalances at the transmitter and receiver. GF-3 is the first full polarization SAR satellite of China launched in August 2016. It is designed so that the channel isolation is better than −35 dB, and the channel imbalance is within 0.5 dB/10°. It provides full-polarized data with swaths of at least 20 km with a resolution of about 8 m (QPSI mode), and a 35 km swath with a resolution of about 25 m (QPSII mode).
GF-3 has been polarimetric-calibrated using polarimetric active radar calibrators’ methods. The calibration of GF-3 quad-pol data was implemented in July 2017 [
2]. However, to validate the method for the estimation of unknown and possibly large phase imbalances, in this experiment, we specifically chose the SAR data that was imaged earlier than July 2017. It means that the used quad-pol SAR data was not calibrated.
In the case of full-polarization data, the comparison object of experimental results is a full-polarization distortion parameters’ estimation method based on common distributed targets [
11], which utilizes the forest area of volume scattering to effectively obtain the channel imbalances with an accuracy of about 0.3 dB/4°. To facilitate the comparison of the two methods, the experimental scene was selected to be located in Jiangmen, Guangdong Province, imaged on 30 December 2016. This scene has both the urban area and forest area. Meanwhile, the quantization between polarization channels should be considered when processing the SLC data of GF3, the details of which can be found in the literature [
12]. The Pauli polarization pseudo-color image of this scene and the phase imbalance estimation results of the two methods are shown in
Figure 5.
As shown in
Figure 5a,b, there were a lot of factories and residential areas in this complex scene. In addition, the area marked by the white dashed box in the upper right corner of
Figure 5a was a densely forested mountainous area, which satisfied the requirements of the quad-pol estimation method for distribution targets. The extraction details of RDBs were the same as in the previous experiment, and the RDB targets obtained are marked as yellow crosses in
Figure 5a. It can be seen that they were densely distributed in areas of buildings, and a small number of sporadic targets in areas such as farmland and mountain forests could be considered discrete values, which had little effect on the final statistical results.
Figure 5c,d shows the results obtained by the distributed-target-based method, for which the transmitting phase imbalance was −91.1°and the receiving phase imbalance was 52.3°.
Figure 5e,f shows the estimation results that come from RDBs. HH and VH were used to obtain the transmitting phase imbalance, from
Figure 5f, the peak position of the normal fitting of the left peak was −94.5°. Similarly, it can be seen in
Figure 5e that the receiving phase imbalance at the right peak was 49.3°.
As seen in
Figure 5e,f, there were two peaks in one distribution histogram. It was the representation of the 180-degree phase ambiguity problem, which was caused by the fact that the phase relations of Equations (5a) and (5b) are not distinguished in the RDBs extraction. The difference of about 180° between the two peaks was consistent with the analysis in
Section 2. For the case of estimating the receiving phase imbalance, we marked these RDB targets of the two peaks separately in the SAR scene, as shown in
Figure 5b. The yellow rings correspond to the right peak, and the blue triangles correspond to the left peak in
Figure 5e. As the orientations of the buildings in the same residential area tend to be consistent, the rotation angles of these effective dihedrals were likely to be consistent, which caused the RDB of the same peak to cluster according to geographic location.
Next, the 180-degree phase ambiguity determination needed to be executed to determine the final phase imbalance estimation result. The experiment was based on the case of receiving phase imbalance estimation, which is shown in
Figure 6. Some RDBs corresponding to the left peak in
Figure 5e appeared on the factory circled in
Figure 6a. It is an L1-level image, where the radar illuminates from the left to the right on the image.
Figure 6b is a Google optical image of the corresponding area, where the circled building is the research target. It can be seen that the
effective dihedral of the factory building should be composed of protruding bars on the roof and the roof plane. Further, the orientations of the walls of the buildings in this area were all the same, facing the direction perpendicular to the road. From the SAR image, we could obtain about 96 m in the direction of
(the direction is shown in
Figure 2) and about 42 m in the azimuth direction. Then,
was estimated to be about 23.6°. The incidence angle of this scene was 26.17°. Assuming that the scene had no obvious slope, the local incidence angle was considered to be 26.17°. Using Equation (11), we could obtain the rotation angle
of these
effective dihedrals, which was about −25.9 degrees. From
Figure 1b, we can see that the phase relationship between HH and HV of these
effective dihedrals was opposite, which was not the assumed case in the method (Equation (5a)). Therefore, the corresponding peak (left peak in
Figure 5e) was not the target peak. The receiving phase imbalance of the scene finally obtained with the proposed method was 49.3°, which was consistent with the result of the distributed-targets-based method.
It can also be proved that the left peak was the target peak in the estimation of transmitting phase imbalance; the similar derivation process was not repeated here. From the experimental results, the phase imbalance estimation differences between the two methods turned out to be 3.4 degrees and 3.0 degrees. Considering both methods’ estimation errors, the experimental results were reasonable and in line with expectations. This experiment fully verified the correctness, accuracy, and completeness of the proposed method.
3.3. The Receive Phase Imbalance Estimation Using Uncalibrated GF-3 Dual-Pol Data
In this part, we apply the proposed phase imbalance estimation method to the uncalibrated GF-3 dual-polarization data. The scene was imaged in the urban area of Beijing, which was the FSI mode, ascending, and with an incidence angle of 30.77°. The operation of the experiment was the same as in the previous two experiments. The result is shown in
Figure 7.
Figure 7a shows the polarization pseudo-color image of a local area in the Beijing scene, where the RDB targets marked in blue triangles correspond to the left peak in
Figure 7b, and the RDBs marked in yellow circles correspond to the right peak.
Figure 7b shows the statistical results of the phase imbalance estimation of the extracted total of 528,819 RDB targets from the entire scene. It can be seen that the peak on the left was dominant, and its phase imbalance was estimated to be −108.1°. The peak on the right was relatively weak, it corresponded to an estimated phase imbalance of about 70°.
The phenomenon that one peak was much stronger than the other peak was seen in this experiment. When a large number of buildings in an urban area have the same orientation (i.e., most buildings in Beijing have a north-south orientation), the orientation of the walls corresponding to the effective dihedrals is the same. Then, the extracted RDBs correspond to the same type mostly, resulting in a strong peak.
However, since the data was not calibrated, we cannot judge whether the dominant peak was the target peak. The 180-degree phase ambiguity determination still needed to be performed.
The experiment is shown in
Figure 8. It was a local area of this scene, in which there were a large number of RDBs marked in blue and a small number of RDBs marked in yellow, which corresponds to the two peaks, respectively, in
Figure 7b. Two typical architectural targets were selected as the research objects. One was the building cluster corresponding to the circled blue RDBs, which was the Beijing West Railway Station and nearby buildings, almost facing due south, as shown in
Figure 8c. The other was the building corresponding to the yellow RDBs, which was the Beijing Electric Power Hospital, as shown in
Figure 8d.
The imaging range of the entire scene is displayed in
Figure 8b, the latitude and longitude of the upper-left corner were (40.150°N, 115.912°E), and the lower-left corner was (39.669°N, 116.031°E). After conversion, the difference between the two points was 53 km in the north-south direction and 10 km in the east-west direction. The angle about the north-south direction was about 10.7°.
For the wall of the buildings facing south, the angle
was equivalent to about 79.3°. Meanwhile, the incidence angle of this scene was 30.77°. Assuming that the scene had no obvious slope, the local incidence angle should be about 30.77°. Using Equation (11), we could obtain the rotation angle
of these
effective dihedrals was about −80.78 degrees. However, note that the only meaningful range for
is between −45°and 45° [
24]. Since urban blocks are usually rectangular, rotations exceeding ±45° begin to present the orthogonal sides as the stronger backscatter target [
19]. Therefore, the value of
should be about –10.7°, for which the western wall should be the main backscatter target. Then, the rotation angle
of these
effective dihedrals should be about 12.4°.
From
Figure 1b, it can be seen that these RDBs should belong to the type of phase relationship of Equation (5a), which was exactly the condition used in the method. Therefore, the left peak in
Figure 7b should be the target peak. The phase imbalance of this scene was eventually estimated to be −108.1°.
While the building shown in
Figure 8d was located along the road, the wall of the main building facing the road was the structure corresponding to the strong scatter targets, whose orientation was not the north-south direction. It can be seen that it was rotated about 45 degrees from the north-south direction, then, the phase characteristic of HH and HV changed. For the corresponding RDBs, the yellow marks in the figure indicate that they belonged to the right peak in
Figure 7b, which was not the target peak.
In general, the estimated rotation angle
of the
effective dihedral may not have been a completely accurate value, because the estimation of
was affected by image parameters, geometric accuracy, etc., and the local incidence angle
was also based on the assumption that the ground was flat. However, the phase difference of
from one peak to another peak in
Figure 1b was 45 degrees, which provided enough redundancy and insurance for the estimation errors of the rotation angle
. Therefore, the judgment of the 180-degree phase ambiguity determination of the two experiments should have been correct and valid.
Since GF-3 dual-polarization data were not calibrated, the result could not be verified temporarily, but a meaningful reference value has been obtained for the first time.
3.4. The Phase Imbalance Estimation Experiments on Multiple Band Airborne SAR Data
The experiments in this subsection apply the method to airborne SAR data of different bands to verify the universality and robustness of the method.
The Aerospace Information Research Institute, Chinese Academy of Sciences (AIRCAS) led a development of an airborne multidimensional space joint-observation (SARMSJosSAR) system and carried out data acquisition experiments [
25]. The airborne multidimensional SAR system includes six bands in total, among which the Ka-band has some problems, and here, a total of five bands of P, L S, C, and X data were used accordingly in this experiment. The system parameters of this airborne multi-band data are shown in the
Table 1. “Br” is bandwidth, “Fsr” is the sampling rate, “Xbin” is the pixel interval in the azimuth direction, and “Rbin” is the pixel interval in the range direction.
The data set used in the experiment was imaged in the Houhai area, east of Wanning, Hainan Province, on 25 December 2020. This data set has not undergone a calibration operation, the distortion parameters of each band are unknown. We chose a scene that contained a village and a forest area, as seen in
Figure 9. The airborne SAR image range was small, encompassing only a few towns and villages. As can be seen from the figure, the color of some ground surface changed significantly with the change of radar wavelength.
The experiment was conducted in the same way as in
Section 3.2. The reference values were obtained by the distributed-targets-based polarimetric distortion estimation method. In the experiment, we chose the mountain forest region on the upper right of this scene as the object. The method proposed in the paper estimated the phase imbalance parameters using a limited number of village buildings in the scene. Meanwhile, the points marked in yellow in
Figure 9 are the extracted RDB targets, for which we appropriately relaxed the K value to 2 or 1.5 to ensure that enough point targets were extracted. With the results of a forest-based quad-pol calibration method as a reference, we skipped the determination for 180-degree phase ambiguity. The results are shown in
Table 2.
From
Table 2, it can be seen that the estimation differences between the two methods in most cases were within the range of about 5 degrees. The S-band transmitting result that reached a level of 10 degrees was regarded as a discrete value, which may have been a statistical error caused by insufficient target points, for which we think it had no significant representative significance. Interestingly, the phase imbalance results estimated by RDBs were generally larger than those estimated by the distributed targets in this experiment. The systematic difference may have come from the unknown crosstalk of the airborne polarimetric system.
As for the reasons why the results were not ideal, firstly, the SAR data quality of the airborne data was not very good, including the poor radiation quality and unsatisfactory polarization isolation. In addition, due to the small image scene of the airborne data, the quantity of the RDBs obtained in this experiment was small (about 20,000, and even less for some bands’ data), which reduced the accuracy of the estimation results.
Nevertheless, the estimation differences were still within an allowable range in such complex multi-band scenes with poor data quality. From the results, the estimated errors did not show a clear correlation with the band, and it was seen that the method worked effectively in these bands. Therefore, it was shown that the method could be used for a rough estimation in complex airborne SAR data conditions. The method had a significant contribution to the quality assurance of full-polarization data and dual-polarization data.