A Range Ambiguity Suppression Processing Method for Spaceborne SAR with Up and Down Chirp Modulation

Range ambiguity is one of the factors which affect the SAR image quality. Alternately transmitting up and down chirp modulation pulses is one of the methods used to suppress the range ambiguity. However, the defocusing range ambiguous signal can still hold the stronger backscattering intensity than the mainlobe imaging area in some case, which has a severe impact on visual effects and subsequent applications. In this paper, a novel hybrid range ambiguity suppression method for up and down chirp modulation is proposed. The method can obtain the ambiguity area image and reduce the ambiguity signal power appropriately, by applying pulse compression using a contrary modulation rate and CFAR detecting method. The effectiveness and correctness of the approach is demonstrated by processing the archive images acquired by Chinese Gaofen-3 SAR sensor in full-polarization mode.


Introduction
In Synthetic Aperture Radar (SAR) imaging, range ambiguity is caused by the echoes of the previous and latter transmitted pulses scattered from undesired range zones [1]. The round-trip propagation time of the ambiguous signals differ from that of the desired signal by a multiple of the reciprocal pulse repetition frequency (PRF). It is more severe in spaceborne cases than airborne cases with the increase of PRF, which may dramatically affect the quality of the reconstructed image. Thus, range ambiguity suppression is a highly regarded problem in SAR signal processing.
Compared with a variety of azimuth ambiguity suppression techniques [2][3][4], the study of range ambiguity suppression mainly focused on the radar system. There are mainly two kinds of range ambiguity suppression techniques for SAR systems. One kind of ideas try to modify the way of transmitting and receiving, so as to block the receiving of the ambiguous energy, such as elevation null spacing [5], multiple elevation beams [6], azimuth phase coding(APC) [7,8], and digital beam forming on receiving which can forms an equivalent narrow beam of high antenna gain pointing to the imaging area and suppressing the ambiguity [9]. The others are algorithms which aim to disperse the peak ambiguous energy to reduce the impaction of ambiguity, for example, the finite depth of focus algorithm which automatically defocuses the range ambiguities, thereby reducing their peak levels [10,11], and more commonly modulation of successively transmitted pulses by alternating up and down chirps which spreads the ambiguities in the range direction [12,13]. This technique significantly reduces the ambiguity peaks and performs well for multi-polarization spaceborne SAR where a serious decline in the ambiguity ratio appears in cross polarization mode. With the development of new SAR systems, some new techniques have been put forward in recent years, for example, coding the transmitted pulses using an orthogonal frequency-division multiplexing (OFDM) pulse [14], frequency division segregation scheme of uniformly dividing the entire bandwidth into several smaller bands according to the ground-range resolution [15]. In order to solve the conflict between range and azimuth ambiguities in spaceborne high-resolution and wide-swath (HRWS) SAR, some new techniques based on frequency diverse array (FDA) are proposed [16,17]. Furthermore, the performances of range ambiguity suppression in multichannel ScanSAR with up and down chirp modulation have been studied in [13].
It can be seen that in order to adapt to the new mode SAR, many improvements in radar system design for range ambiguity suppression have been put forward. However, there is a lack of analysis and solutions to the problems that may exist in these system designs. Take the technique of up and down chirp modulation for example, it will show severe destructions to the image quality under certain conditions where targets appearing in the ambiguous area have a stronger backscattering intensity than the imaging targets of interest, such as ships or towns in the ambiguous area for the imaging targets like rivers or mountains. The energy of ambiguous area after SAR imaging process are still too high to reduce the ambiguity ratio which leads to line defocusing ambiguity in the image, having severe impact on visual effects and subsequent applications.
Gaofen-3 satellite is the first full-polarization SAR sensor of China, which implemented alternatively up and down chirp modulation for range ambiguity suppression. The line defocusing ambiguity problem appeared in certain images acquired in full-polarization mode. Therefore, this paper proposed a novel hybrid range ambiguity suppression method for up and down chirp modulation. The method can obtain the ambiguity area image and reduce the ambiguity signal power appropriately, by applying pulse compression using contrary modulation rate and Constant False Alarm Rate (CFAR) detecting method.
In Section 2 the effect of range ambiguity under the condition of transmitting up and down chirp modulation signals is analyzed. In Section 3 the principles of the proposed range ambiguity suppression technique are introduced. Section 4 gives the processing results of Gaofen-3 full-polarization strip mode data and conclusion is summarized in the Section 5.

Range Ambiguity Analysis by Alternatively Up and Down Chirp Modulation
As shown in Figure 1, range ambiguous signals are echoes arrived at the wrong time and scattered from undesired range zones. Considering a target (point A in Figure 1) in the main region with range history R(t) where t presents the azimuth time, the echoes from previous and latter transmitted pulses will arrive as ambiguous signals with the slant range history as [18]: where c is the velocity of light, prf is the pulse repetition frequency, n is the ambiguous number (negative part presents previous pulses when positive part presents latter pulses). The signal echoes consist of mainlobe area parts and ambiguous area parts are shown below: where τ is the range time.
Supposing the system works by sending alternately up and down chirp signals with pulse length T p and modulation rate k r , We can get the echo of ambiguous signals as: where σ is the backscattering intensity, λ is the wave length, W r (·) and W a (·) are the antenna patterns of range and azimuth direction.  Figure 1. This is a figure which shows the principle of range ambiguity.
We can get the signals for main area as: We can see from (3) and (4) that the ambiguous signals with odd n show contrary modulation rate to the main signals, whereas those with even n show the same modulation rate.
The filtering of an up chirp with a down chirp filter will result in defocusing phenomenon along range direction and the pulse compression result as follow [12]: For comparison, the result of the correct filtering of the down chirp is [12]: where krTp 2 is the compression gain.
We can see that the amplitude of mismatched up-chirp ambiguous signals (odd n) is lowered by The signal power of main area is unchangeable once the system antenna pattern and imaging targets are fixed, which leaves ambiguity suppression work to the reduction of ambiguous signal power. Due to the contrary modulation rate to the main signals with odd ambiguous number, the method of alternating up and down chirp modulation is provided with lower ambiguous signal power owing to the defocusing of mismatched range filter. Therefore, it can be widely adopted in the multi-polarization spaceborne SAR where a serious decline in the range-ambiguity-to-signal-ratio appears in cross polarization mode.
We can get the range-ambiguity-to-signal-ratio (RASR) as [18]: We can get the signals for main area as: We can see from (3) and (4) that the ambiguous signals with odd n show contrary modulation rate to the main signals, whereas those with even n show the same modulation rate.
The filtering of an up chirp with a down chirp filter will result in defocusing phenomenon along range direction and the pulse compression result as follow [12]: For comparison, the result of the correct filtering of the down chirp is [12]: where k r T p 2 is the compression gain.
We can see that the amplitude of mismatched up-chirp ambiguous signals (odd n) is lowered by 1/ 2·k r T 2 p . The signal power of main area is unchangeable once the system antenna pattern and imaging targets are fixed, which leaves ambiguity suppression work to the reduction of ambiguous signal power. Due to the contrary modulation rate to the main signals with odd ambiguous number, the method of alternating up and down chirp modulation is provided with lower ambiguous signal power owing to the defocusing of mismatched range filter. Therefore, it can be widely adopted in the multi-polarization spaceborne SAR where a serious decline in the range-ambiguity-to-signal-ratio appears in cross polarization mode. We can get the range-ambiguity-to-signal-ratio (RASR) as [18]: where η represents the incidence angle for a certain ambiguous number, G(·) represents the azimuth antenna patterns for a certain incidence angle, T s represents the data recording time.
As can be seen from (7), the ambiguous signal power is mainly determined by the backscattering intensity of the ambiguous area when the range antenna pattern has already been fixed according to the system design. At certain condition, range antenna pattern and the lower power of the defocusing signal can no longer suppress the ambiguous signal power being lower than the main area signal power, when the backscattering intensity of ambiguous area grows strong enough. This can result in ambiguous range defocusing targets appearing on the product image with a shape of line along range direction, having severe impact on visual effects and subsequent applications.

Ambiguity Suppression Method
In order to solve the problem of poor image quality caused by strong ambiguous signals, we developed an ambiguity suppression method aimed at ambiguous signals with odd ambiguous number. The general flowchart is given in Figure 2. where η represents the incidence angle for a certain ambiguous number, G(·) represents the azimuth antenna patterns for a certain incidence angle, Ts represents the data recording time.
As can be seen from (7), the ambiguous signal power is mainly determined by the backscattering intensity of the ambiguous area when the range antenna pattern has already been fixed according to the system design. At certain condition, range antenna pattern and the lower power of the defocusing signal can no longer suppress the ambiguous signal power being lower than the main area signal power, when the backscattering intensity of ambiguous area grows strong enough. This can result in ambiguous range defocusing targets appearing on the product image with a shape of line along range direction, having severe impact on visual effects and subsequent applications.

Ambiguity Suppression Method
In order to solve the problem of poor image quality caused by strong ambiguous signals, we developed an ambiguity suppression method aimed at ambiguous signals with odd ambiguous number. The general flowchart is given in Figure 2.  From previous analysis results we know that the ambiguous signals with odd n show contrary modulation rate to the main signals, whereas even n show the same rate. Therefore, we can get the image result of ambiguous area with odd n through range compression with contrary modulation rate to main signals and simple azimuth compression. Actually, SAR imaging can be regarded as a process to gather the scattered signal power again. Through ambiguous area imaging, we can get the focused ambiguous signal power rather than that spreads out in echo data. Then the focused From previous analysis results we know that the ambiguous signals with odd n show contrary modulation rate to the main signals, whereas even n show the same rate. Therefore, we can get the image result of ambiguous area with odd n through range compression with contrary modulation rate to main signals and simple azimuth compression. Actually, SAR imaging can be regarded as a process to gather the scattered signal power again. Through ambiguous area imaging, we can get the focused ambiguous signal power rather than that spreads out in echo data. Then the focused ambiguous area can be detected and obtained by means of CFAR. Next, the ambiguity suppression is applied to the detected ambiguous area by simply reducing the amplitude, with no destruction to the phase information. Finally, echo data with suppressed ambiguous signals are obtained by anti-compression in range and azimuth direction. The image after ambiguity suppression are provided with Chirp-Scaling algorithm. The details of each step are as follows.

Ambiguous Area Imaging
The ambiguous area imaging method is based on a classical SAR imaging algorithm called Range-Doppler algorithm (RDA) [19]. In simple terms, RDA includes three steps which are range compression, range cell migration correction (RCMC) and azimuth compression, aiming at main area imaging. To get the image of the ambiguous area, we need to make some changes to the algorithm. The general steps are shown in the right parts of Figure 2. The following contents illustrate the specific process.
• Range compression Firstly, we get the signal in range frequency domain by applying the fast Fourier transform (FFT) along range direction to the echo data as shown in Equations (1)- (4): where f τ is the range frequency, f 0 is the carrier frequency. It is obvious that n = 0 presents the main area signals in range frequency domain and others the ambiguous area signals.
A range matched filter with contrary modulation rate to the main area is required to get the pulse compression result of the ambiguous area. The filter in frequency domain is as follows [19]: We can get the range compression result through the following matching filtering process: where F −1 [·] represents the inverse FFT(IFFT). The specific expression is as follows: where s odd (τ,t,n) represents the compression results of odd n, s even (τ,t,n) the even n (including main area with n = 0). We can see that only the ambiguous area with odd n can obtain the right focused range profiles, contributing to the acquisition of the focused image of this ambiguous area.

• Azimuth compression
There are two steps that need to be followed, including RCMC and azimuth compression, according to the RDA. However, the coupling of range and azimuth should remain unchanged to keep the imaging process reversible, which makes it impossible to add RCMC to the imaging process. Thus, we get the image result of ambiguous area through simple azimuth matched filter compression under the assumption of tiny range migration.
Here we consider s odd (τ,t,n) only since the defocusing range compression results can no longer get the focusing image through azimuth processing. We get the range compression results in Doppler frequency domain by azimuth FFT [19]: where f t is the azimuth frequency, R 0 is the center slant range, V is the equivalent velocity between satellite and target at center time. R ' n (f t ) is the range cell migration in the range direction, which has a complex expression that can be found in [19]. The azimuth matched filter in frequency domain is as follows [19]: where R ref represents the reference slant range. According to the SAR imaging algorithm, R ref and V are calculated in the standard imaging process apparently for main area signals [20]. Considering the long slant range of the spaceborne SAR, V of the ambiguous area varies little to the main area. Thus, we just need an appropriate R ref for the ambiguous area. According to the range ambiguity theory, it can be expressed as: where R main represents the calculated reference slant range from standard imaging process. We can get the ambiguous area image result of certain n through the following matching filtering process: where S rd_even (τ,f t ,n) represents the range compression results of even n in Doppler frequency domain.
• Finding the main ambiguous area For better ambiguity suppression, we need to find the right n and R ref for main ambiguous area which shows strongest intensity in the focused image. Considering the mismatched R ref leads to defocusing targets, we just need to test the quality of images acquired by different R ref (n) and select the best focusing one. Generally, the focusing quality of a SAR image can be measured by image entropy [21]: where I(x,y) represents the gray scale value of the SAR image, M and N represent the total number of sample points in azimuth and range direction. The smaller the image entropy, the better the focused image. What needs to be specified is that larger n can bring about smaller power of the corresponding ambiguous area as a result of farther distance from the antenna center. Thus, in this paper we choose n = 1 and n = −1 to find the location of the main ambiguous area.

• Main ambiguous area location
The right n can not only bring right focused main ambiguous area image, but also the center longitude and latitude of the main ambiguous area.
Supposing R(R x ,R y ,R z ) represents the center position coordinates of the main ambiguous area in the earth coordinate system, α and β represent the longitude (positive the East while negative the West) and latitude (positive the North while negative the South) respectively, we can get the following expression: where R T represents the Earth radius at the location of the imaging target.
According to the geometric relationship of the satellite and Earth, we can get the following equations [22]: where R s and V s represent the position and velocity vector of satellite at the center time, f dc the Doppler center frequency, R equ and R pol the radius of the equator and pole of the Earth. R' cen represents the center slant range of main ambiguous area that can be obtained by: cen R = R cen + n/pr f ·c/2 (20) where R cen represents the center slant range of main imaging area. By solving the Equations (18)-(20), we can get the value of α and β. The optical image can be obtained by the information of longitude and latitude, which can be compared with the acquired SAR image to verify the validity of the imaging process.

Ambiguous Area Detecting Based on CFAR
The imaging result of main ambiguous area includes not only ambiguous targets but also range defocusing results of the main area and other ambiguous areas which show weaker backscattering intensity than focusing ones, bringing about even weaker intensity than the focusing ambiguous targets. The problem of ambiguous area target detection in the background of defocusing targets can be translated into the detection for strong backscattering targets under weak background, which is normally solved by a target detecting algorithm called CFAR [23,24].
Target detection based on CFAR is an algorithm that requires a given detection rate and a strong contrast between target and background, which is extensively used in ship detection of SAR image [25,26]. The excellent results of ship detection demonstrate the ability of CFAR in strong scatter points detecting under weak background of SAR images, which is what we need for ambiguous area detection in this paper.
In this paper we choose the double-parameter CFAR detection algorithm to achieve the detection of main ambiguous area targets [25][26][27][28]. This algorithm sets up three sliding windows as shown in Figure 3 and realize the target detection by comparison of the pixel grayscale in target window and the self-adapting threshold which is acquired by false alarm rate (PFA) and background clutter modeling, which is assumed to be Gaussian distribution with mean value µ and standard deviation σ is usually an integrated area composed of strong scattering points, different from ships which have a uniform shape, making it necessary to do some experiments to determine the reasonable length of the three windows.

Background window
Protection window Target window Figure 3. This is a figure that shows the geometric relation of three sliding windows. They are all squares and the inside is the target window, the middle is the protection window, and the outside is the background window. •

Calculation of threshold including μ and σ for the background window
To set the detection threshold T, we need the distribution function of background clutter and the given PFA. Supposing Φ(·) presents the Gaussian distribution function, we can get the following expression: where μ and σ are calculated according to each background window, making the threshold T changing adaptively. In fact, the background clutter may not be in conformity with the standard Gaussian distribution absolutely, forcing us to determine the reasonable choice of t1 through experiments. •

Target detection
If the gray value of a pixel in the target window is Itar, then the target detection criterion is: As mentioned above, we need some experiments to determine the reasonable parameters of the algorithm. Limited to the characteristics of the ambiguous area image itself, this method cannot evaluate the detection quality by the specific quantity of right detecting targets. We can only find a comparatively good value interval by comparing the results from visual observation.

Ambiguity Suppression
The phase information of the input data should stay complete and not damaged to keep the integrity of main area signals, which makes it necessary to the use of amplitude suppression for ambiguity suppression. We divided the image result acquired in Section 3.1 into two parts according to the detection result: where Cde represents the collection of detected pixel positions and Cre represents the collection of remaining pixel positions. The amplitude of detected pixel is reduced by N times, and we can get the result as: Figure 3. This is a figure that shows the geometric relation of three sliding windows. They are all squares and the inside is the target window, the middle is the protection window, and the outside is the background window.
• Set up for the length of sliding windows The size of the general target window should be the same as the minimum size of targets to be detected. The length of the protection window is generally twice the maximum target size, while the background window being selected as two or four times the protection window size. The step length of the sliding window is in accordance with the size of the target window. The ambiguous area image is usually an integrated area composed of strong scattering points, different from ships which have a uniform shape, making it necessary to do some experiments to determine the reasonable length of the three windows.
• Calculation of threshold including µ and σ for the background window To set the detection threshold T, we need the distribution function of background clutter and the given PFA. Supposing Φ(·) presents the Gaussian distribution function, we can get the following expression: where µ and σ are calculated according to each background window, making the threshold T changing adaptively. In fact, the background clutter may not be in conformity with the standard Gaussian distribution absolutely, forcing us to determine the reasonable choice of t 1 through experiments.

• Target detection
If the gray value of a pixel in the target window is I tar , then the target detection criterion is: As mentioned above, we need some experiments to determine the reasonable parameters of the algorithm. Limited to the characteristics of the ambiguous area image itself, this method cannot evaluate the detection quality by the specific quantity of right detecting targets. We can only find a comparatively good value interval by comparing the results from visual observation.

Ambiguity Suppression
The phase information of the input data should stay complete and not damaged to keep the integrity of main area signals, which makes it necessary to the use of amplitude suppression for ambiguity suppression. We divided the image result acquired in Section 3.1 into two parts according to the detection result: where C de represents the collection of detected pixel positions and C re represents the collection of remaining pixel positions. The amplitude of detected pixel is reduced by N times, and we can get the result as: The extent of suppression is alternative for manually input. Generally, we choose the attenuation multiplier N for the aim that ambiguous signal power lower than main area power after standard imaging by 0~5 dB.

Anti-Compression and Imaging
Since the imaging and ambiguity suppression process are reversible, echo data can be obtained by anti-compression in range and azimuth direction to the data after ambiguity suppression. The image after ambiguity suppression are provided with Chirp-Scaling algorithm.
The proposed ambiguity suppression method has two limitations: • Only for ambiguous areas with odd ambiguous numbers; • Only for strip data with tiny range migration.

Results
The performance of the proposed range ambiguity suppression method has been analyzed by processing full-polarization SAR data from Gaofen-3 sensor strip mode which acquired in the near Argun River of Inner Mongolia (the parameters are listed Table 1). Figure 4 shows the HH polarization and VV polarization SAR images acquired by applying standard Chirp-Scaling algorithm to the echo data with no ambiguity suppression process. There are some obvious bright lines along range direction in the top left of the two images which have serious implications on the visual effect and subsequent applications. Considering the technology of alternatively up and down chirp modulation implemented by Gaofen-3 satellite, the bright lines in the image could be the defocusing range ambiguous signals which hold strong intensity.
The extent of suppression is alternative for manually input. Generally, we choose the attenuation multiplier N for the aim that ambiguous signal power lower than main area power after standard imaging by 0~5 dB.

Anti-Compression and Imaging
Since the imaging and ambiguity suppression process are reversible, echo data can be obtained by anti-compression in range and azimuth direction to the data after ambiguity suppression. The image after ambiguity suppression are provided with Chirp-Scaling algorithm.
The proposed ambiguity suppression method has two limitations: • Only for ambiguous areas with odd ambiguous numbers; • Only for strip data with tiny range migration.

Results
The performance of the proposed range ambiguity suppression method has been analyzed by processing full-polarization SAR data from Gaofen-3 sensor strip mode which acquired in the near Argun River of Inner Mongolia (the parameters are listed Table 1). Figure 4 shows the HH polarization and VV polarization SAR images acquired by applying standard Chirp-Scaling algorithm to the echo data with no ambiguity suppression process. There are some obvious bright lines along range direction in the top left of the two images which have serious implications on the visual effect and subsequent applications. Considering the technology of alternatively up and down chirp modulation implemented by Gaofen-3 satellite, the bright lines in the image could be the defocusing range ambiguous signals which hold strong intensity.      Figure 5 shows the image of ambiguous area targets by processing the HH polarization echo data, where (a) represents the azimuth matched filter R ref with n = −1 when (b) represents the n = +1. To find the right n of the main ambiguous area, entropy values of the two are calculated according to Equation (17). The results are as follows: It's obvious that the focusing performance of the result with n = −1 is better than that with n = +1, with which we can suppose that the main ambiguous area of the bright lines is with n = −1. The center latitude and longitude information can be calculated according to the Equations (18)- (20) and Table 1. Figure 6b shows the optical image of the main ambiguous area with the result above. It is the city of Hulun-Buir near the Argun River, which displays a stronger backscatter intensity than the imaging area (Argun River) indeed. Figure 6a shows the corresponding SAR image which is the detail magnification of the top left part in Figure 5a (marked with a red rectangle). The two show a high degree of consistency, which can verify the validity of the ambiguous area imaging process.
Sensors 2018, 18, x 10 of 14 Figure 5 shows the image of ambiguous area targets by processing the HH polarization echo data, where (a) represents the azimuth matched filter Rref with n = −1 when (b) represents the n = +1. To find the right n of the main ambiguous area, entropy values of the two are calculated according to Equation (17). The results are as follows: It's obvious that the focusing performance of the result with n = −1 is better than that with n = +1, with which we can suppose that the main ambiguous area of the bright lines is with n = −1. The center latitude and longitude information can be calculated according to the Equations (18) Figure 6b shows the optical image of the main ambiguous area with the result above. It is the city of Hulun-Buir near the Argun River, which displays a stronger backscatter intensity than the imaging area (Argun River) indeed. Figure 6a shows the corresponding SAR image which is the detail magnification of the top left part in Figure 5a (marked with a red rectangle). The two show a high degree of consistency, which can verify the validity of the ambiguous area imaging process.   The parameters of CFAR used in target detection process are shown in Table 2. We set up the target window as 2, the protection window as 8 and the background window as 32, and we finally choose t1 = 2.5~3.5 for the calculation of threshold, which is decided after many times of experiments. The parameters of CFAR used in target detection process are shown in Table 2. We set up the target window as 2, the protection window as 8 and the background window as 32, and we finally choose t 1 = 2.5~3.5 for the calculation of threshold, which is decided after many times of experiments. In fact, it can reach almost the same detection rate with few error targets for t 1 in this interval. Figure 7 shows the detecting results of ambiguous area targets where left parts present ambiguous area imaging results which are typical parts of Figure 6a (marked with red rectangles with serial numbers), right parts the detecting results with t 1 = 3. It can be seen from the results that most of the strong scatters in imaging results are detected which presents the composition of the ambiguous area. The results show the accuracy and validity of the detection algorithm. We selected an attenuation multiplier N = 100 according to the power of the ambiguous area and main imaging area. Actually, the performance of ambiguity suppression depends on the detection results more than the attenuation multiplier N. Thus, there is no need to be entangled in the precise choice of N as long as it can suppress the ambiguous power effectively.
The image results after ambiguity suppression and chirp-scaling algorithm imaging are shown in Figure 8 (HH-polarization) and Figure 9 (VV-polarization). Comparing with the results in Figure 4a,b, respectively, the bright lines in the top left part of the image are almost invisible. It indicates that the ambiguity signal power in data have been well suppressed and the effect of ambiguity suppression is remarkable. The results show the accuracy and validity of the proposed ambiguity suppression algorithm.
choice of N as long as it can suppress the ambiguous power effectively.
The image results after ambiguity suppression and chirp-scaling algorithm imaging are shown in Figure 8 (HH-polarization) and Figure 9 (VV-polarization). Comparing with the results in Figures 4a,b, respectively, the bright lines in the top left part of the image are almost invisible. It indicates that the ambiguity signal power in data have been well suppressed and the effect of ambiguity suppression is remarkable. The results show the accuracy and validity of the proposed ambiguity suppression algorithm.

Conclusions
In this paper, we have proposed a new range ambiguity suppression method for up and down chirp modulation. The most important feature is the ambiguous target detection based on pulse compression using a contrary modulation rate and CFAR detection method. The method can obtain the ambiguity area image and reduce the ambiguity signal power appropriately. The effectiveness and correctness of the approach are demonstrated by processing archive images acquired by the Gaofen-3 SAR sensor in full-polarization mode. Therefore, the proposed method can be viewed as an alternative way to improve the image quality for images with range ambiguity. At the same time, the method also has some limitations as mentioned in Section 3.4, and these limitations along with the problem of scientifically evaluating the detection results are the direction of the research in the future.