Next Article in Journal
Filter Selection for Optimizing the Spectral Sensitivity of Broadband Multispectral Cameras Based on Maximum Linear Independence
Next Article in Special Issue
Speckle Filtering of GF-3 Polarimetric SAR Data with Joint Restriction Principle
Previous Article in Journal
Privacy-Preserving Authentication Using a Double Pseudonym for Internet of Vehicles
Previous Article in Special Issue
Water Detection in Urban Areas from GF-3
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Laboratory of spatial information intelligent processing system, Institute of Electronics, Chinese Academy of Sciences, Suzhou 215000, China
2
Institute of Electronics, Chinese Academy of Sciences, Beijing 100190, China
3
China Center for Resources Satellite Data and Application, Beijing 100094, China
*
Author to whom correspondence should be addressed.
Submission received: 12 March 2018 / Revised: 30 April 2018 / Accepted: 4 May 2018 / Published: 7 May 2018
(This article belongs to the Special Issue First Experiences with Chinese Gaofen-3 SAR Sensor)

Abstract

:
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.

1. 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.

2. 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]:
R n ( t ) = R ( t ) + c / 2 · n / p r f ( n = ± 1 , ± 2 , ± 3 )
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:
s ( τ , t ) = s m a i n ( τ , t ) + n s a m b i g u i t y ( τ , t , n )
where τ is the range time.
Supposing the system works by sending alternately up and down chirp signals with pulse length Tp and modulation rate kr, We can get the echo of ambiguous signals as:
s a m b i g u i t y ( τ , t , n ) = σ n · W a ( t , n ) · W r [ τ 2 R n ( t ) / c ] · exp [ j 4 π / λ · R n ( t ) ] · exp { ( 1 ) | n | j π k r · [ τ 2 R n ( t ) / c ] 2 }
where σ is the backscattering intensity, λ is the wave length, Wr(·) and Wa(·) are the antenna patterns of range and azimuth direction.
We can get the signals for main area as:
s m a i n ( τ , t ) = σ · W a ( t ) · W r [ τ 2 R ( t ) / c ] · exp [ j 4 π / λ · R ( t ) ] · exp { j π k r · [ τ 2 R ( t ) / c ] 2 }
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]:
s d e ( τ ) = 1 / 2 · r e c t [ τ / 2 · T p ] · exp ( j π k r τ 2 / 2 )
For comparison, the result of the correct filtering of the down chirp is [12]:
s d o w n ( τ ) = exp ( j π / 4 ) · k r T p 2 sin c [ π k r T p τ ]
where krTp2 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 p 2 .
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]:
R A S R = n = T s / 2 T s / 2 G 2 [ η ( t , n ) ] · σ ( t , n ) sin [ η ( t , n ) ] / R n 3 ( t ) d t T s / 2 T s / 2 G 2 [ η ( t ) ] · σ ( t ) sin [ η ( t ) ] / R 3 ( t ) d t
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.

3. 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 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.

3.1. 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):
S ( f τ , t , n ) = σ n · W a ( t , n ) · W r [ f τ ] · exp [ j 4 π ( f 0 + f τ ) / c · R n ( t ) ] · exp [ ( 1 ) | n | j π f τ 2 / k r ]
where fτ is the range frequency, f0 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]:
H c o n t r a r y ( f τ ) = exp ( j π f τ 2 / k r )
We can get the range compression result through the following matching filtering process:
s r ( τ , t , n ) = F 1 [ S ( f τ , t , n ) H c o n t r a r y ( f τ ) ]
where F−1[·] represents the inverse FFT(IFFT).
The specific expression is as follows:
s o d d ( τ , t , n ) = σ n · W a ( t , n ) · sin c [ τ 2 R n ( t ) / c ] · exp [ j 4 π R n ( t ) / λ ]
s e v e n ( τ , t , n ) = σ n · W a ( t , n ) · W r [ ( τ 2 R n ( t ) / c ) / 2 ] · exp [ j 4 π R n ( t ) / λ ] · exp [ j π k r / 2 · ( τ 2 R n ( t ) / c ) 2 ]
where sodd(τ,t,n) represents the compression results of odd n, seven(τ,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 sodd(τ,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]:
S r d _ o d d ( τ , f t , n ) = σ n · W a ( f t , n ) · sin c [ τ 2 R n ( f t ) / c ] · exp [ j 4 π R 0 / λ · ( 1 ( λ f t / 2 V ) 2 ) ]
where ft is the azimuth frequency, R0 is the center slant range, V is the equivalent velocity between satellite and target at center time. Rn(ft) 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]:
H a z i m u t h ( f t ) = exp [ j 4 π R r e f / λ · ( 1 ( λ f t / 2 V ) 2 1 ) ]
where Rref represents the reference slant range.
According to the SAR imaging algorithm, Rref 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 Rref for the ambiguous area. According to the range ambiguity theory, it can be expressed as:
R r e f = R m a i n + c / 2 · n / p r f
where Rmain 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:
s i m ( τ , t , n ) = F 1 { [ S r d _ o d d ( τ , f t , n ) + S r d _ e v e n ( τ , f t , n ) ] H a z i m u t h ( f t ) }
where Srd_even(τ,ft,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 Rref for main ambiguous area which shows strongest intensity in the focused image. Considering the mismatched Rref leads to defocusing targets, we just need to test the quality of images acquired by different Rref(n) and select the best focusing one. Generally, the focusing quality of a SAR image can be measured by image entropy [21]:
E i m = x = 1 M y = 1 N I ( x , y ) S log I ( x , y ) S S = x = 1 M y = 1 N I ( x , y )
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(Rx,Ry,Rz) 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:
R x = R T · cos α · cos β R y = R T · cos α · sin β R z = R T · sin α
where RT 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]:
| R s R | = R c e n V s × ( R s R ) R c e n = λ f d c 2 ( R x 2 + R y 2 ) R e q u + R z 2 R p o l = 1
where Rs and Vs represent the position and velocity vector of satellite at the center time, fdc the Doppler center frequency, Requ and Rpol 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:
c e n R = R c e n + n / p r f · c / 2
where Rcen 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.

3.2. 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 σ in background window. The general flowchart is given in the left part of Figure 2. The specific steps are as follows.
• 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:
T = σ · t 1 + μ t 1 = Φ 1 ( 1 P F A )
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:
I t a r > σ · t 1 + μ
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.

3.3. 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:
s i m ( x , y ) = i , j C r e s i m ( i , j ) + p , q C de s i m ( p , q )
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:
s p r ( x , y ) = i , j C r e s i m ( i , j ) + p , q C de s i m ( p , q ) / N
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.

3.4. 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.

4. 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 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:
E n = 1 = 2.8438 E n = 1 = 0.1278
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.
L o n = 120.921 ° E L a t = 48.833 ° N
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. In fact, it can reach almost the same detection rate with few error targets for t1 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 t1 = 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.

5. 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.

Author Contributions

Conceptualization, X.Q. and C.D.; Methodology, X.Q., X.W. and B.L.; Software, X.W.; Validation, X.Q., X.W. and B.H.; Formal Analysis, X.Q.; Investigation, X.Q. and C.D.; Resources, Q.C. and B.L.; Data Curation, Q.C. and X.W.; Writing-Original Draft Preparation, X.W.; Writing-Review & Editing, X.W. and X.Q.; Visualization, C.D.; Supervision, C.D. and B.L.; Project Administration, X.Q.; Funding Acquisition, C.D.

Funding

This work was supported by the National Science Foundation of China (grant 61331017) and the GF-3 High-Resolution Earth Observation System under Project (30-Y20A12-9004-15/16, 03-Y20A11-9001-15/16).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Harger, R.O. Synthetic Aperture Radar Systems: Theory and Design; H36-1; Academic Press: New York, NY, USA, 1970. [Google Scholar]
  2. Wu, Y.; Yu, Z.; Xiao, P.; Li, C. Azimuth ambiguity suppression based on minimum mean square error estimation. In Proceedings of the IEEE Geoscience and Remote Sensing Symposium, Milan, Italy, 26–31 July 2015; pp. 2425–2428. [Google Scholar]
  3. Guo, X.; Gao, Y.; Wang, K.; Liu, X. Suppression of Azimuth Ambiguities of Strong Point-Like Targets for Multichannel SAR Systems. IEEE Geosci. Remote Sens. Lett. 2017, 14, 1046–1050. [Google Scholar] [CrossRef]
  4. Wang, K.; Chen, J.; Yang, W.; Li, Z. Suppression of azimuth ambiguities in spaceborne stripmap SAR using accurate restoration modeling. In Proceedings of the IEEE Geoscience and Remote Sensing Symposium, Milan, Italy, 26–31 July 2015; pp. 2429–2432. [Google Scholar]
  5. Grifflths, H.D.; Mancini, P. Ambiguity Suppression in Sars Using Adaptive Array Techniques. In Proceedings of the International IEEE Geoscience and Remote Sensing Symposium (IGARSS ’91), Espoo, Finland, 3–6 June 1991; pp. 1015–1018. [Google Scholar]
  6. Krieger, G.; Huber, S.; Villano, M.; Younis, M.; Rommel, T.; Dekker, P.L.; de Almeida, F.Q.; Moreira, A. CEBRAS: Cross elevation beam range ambiguity suppression for high-resolution wide-swath and MIMO-SAR imaging. In Proceedings of the IEEE Geoscience and Remote Sensing Symposium, Milan, Italy, 26–31 July 2015; pp. 196–199. [Google Scholar]
  7. Yang, J.; Sun, G.; Wu, Y.; Xing, M. Range ambiguity suppression by azimuth phase coding in multichannel SAR systems. In Proceedings of the IET International Radar Conference, Xi’an, China, 14–16 April 2013; pp. 1–5. [Google Scholar]
  8. Guo, L.; Tan, X.; Dang, H. Range ambiguity suppression for multi-channel SAR system near singular points. In Proceedings of the IEEE Geoscience and Remote Sensing Symposium, Beijing, China, 10–15 July 2016; pp. 2086–2089. [Google Scholar]
  9. Krieger, G.; Gebert, N.; Moreira, A. Multidimensional Waveform Encoding: A New Digital Beamforming Technique for Synthetic Aperture Radar Remote Sensing. IEEE Trans. Geosci. Remote Sens. 2007, 46, 31–46. [Google Scholar] [CrossRef] [Green Version]
  10. Tomiyasu, K. Tutorial review of synthetic-aperture radar (SAR) with applications to imaging of the ocean surface. Proc. IEEE 1978, 66, 563–583. [Google Scholar] [CrossRef]
  11. Tomiyasu, K. Image processing of synthetic aperture radar range ambiguous signals. IEEE Trans. Geosci. Remote Sens. 1994, 32, 1114–1117. [Google Scholar] [CrossRef]
  12. Mittermayer, J.; Martinez, J.M. Analysis of range ambiguity suppression in SAR by up and down chirp modulation for point and distributed targets. In Proceedings of the IEEE Geoscience and Remote Sensing Symposium, Toulouse, France, 21–25 July 2003; pp. 4077–4079. [Google Scholar]
  13. Mo, H.; Zeng, Z. Investigation of multichannel ScanSAR with up and down chirp modulation for range ambiguity suppression. In Proceedings of the IEEE Geoscience and Remote Sensing Symposium, Beijing, China, 10–15 July 2016; pp. 1130–1133. [Google Scholar]
  14. Riché, V.; Méric, S.; Baudais, J.Y.; Pottier, E. Investigations on OFDM Signal for Range Ambiguity Suppression in SAR Configuration. IEEE Trans. Geosci. Remote Sens. 2014, 52, 4191–4197. [Google Scholar] [CrossRef]
  15. Pyne, B.; Ravindra, V.; Saito, H. An improved pulse repetition frequency selection scheme for synthetic aperture radar. In Proceedings of the IEEE Radar Conference, Paris, France, 9–11 September 2015; pp. 257–260. [Google Scholar]
  16. Lin, C.; Huang, P.; Wang, W.; Li, Y.; Xu, J. Unambiguous Signal Reconstruction Approach for SAR Imaging Using Frequency Diverse Array. IEEE Geosci. Remote Sens. Lett. 2017, 14, 1628–1632. [Google Scholar] [CrossRef]
  17. Xu, J.; Liao, G.; So, H.C. Space-Time Adaptive Processing with Vertical Frequency Diverse Array for Range-Ambiguous Clutter Suppression. IEEE Trans. Geosci. Remote Sens. 2016, 54, 5352–5364. [Google Scholar] [CrossRef]
  18. Curlander, J.C.; Mcdonough, R.N. Synthetic Sperture Radar: Systems and Signal Processing; Publishing House of Electronics Industry: Beijing, China, 2006. [Google Scholar]
  19. Cumming, I.G.; Wong, F.H. Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation; Publishing House of Electronics Industry: Beijing, China, 2012. [Google Scholar]
  20. Wei, Z.Q. Synthetic Aperture Radar Satellite; Science Publishing Company: Beijing, China, 2001. [Google Scholar]
  21. Zhang, S.; Liu, Y.; Li, X. Fast Entropy Minimization Based Autofocusing Technique for ISAR Imaging. IEEE Trans. Signal Process. 2015, 63, 3425–3434. [Google Scholar] [CrossRef]
  22. Yuan, X.K. Target location method for synthetic aperture radar satellite. Spacecr. Eng. 1998, 232–238. [Google Scholar]
  23. Hao, C.P. A Two Parameter CFAR Detector in K-Distribution Clutter. J. Electron. Inf. Technol. 2007, 29, 756–759. [Google Scholar]
  24. Ai, J.Q.; Qi, X.Y.; Yu, W.D. Improved two parameter CFAR ship detection algorithm in SAR images. J. Electron. Inf. Technol. 2009, 31, 2881–2885. [Google Scholar]
  25. Xu, J.; Ma, Y.; Li, J.; Peng, Y. Optimizing CFAR-based SAR Target Detection Algorithm for DSP Platform. In Proceedings of the IEEE Sixth International Conference on Instrumentation & Measurement, Computer, Communication and Control, Harbin, China, 21–23 July 2016; pp. 492–495. [Google Scholar]
  26. Leng, X.; Ji, K.; Yang, K.; Zhou, H. A Bilateral CFAR Algorithm for Ship Detection in SAR Images. IEEE Geosci. Remote Sens. Lett. 2015, 12, 1536–1540. [Google Scholar]
  27. Chen, P.; Yang, J.; Wang, J. Marine targets detection using GF-3 SAR data. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Fort Worth, TX, USA, 23–28 July 2017; pp. 1978–1980. [Google Scholar]
  28. Greet, L.D.; Harris, S. Detection and clutter suppression using fusion of conventional CFAR and two parameter CFAR. In Proceedings of the IEEE International Conference on Electronics Computer Technology, Kanyakumari, India, 8–10 April 2011; pp. 216–219. [Google Scholar]
Figure 1. This is a figure which shows the principle of range ambiguity.
Figure 1. This is a figure which shows the principle of range ambiguity.
Sensors 18 01454 g001
Figure 2. Figure which shows the detailed steps of the proposed range ambiguity suppression method.
Figure 2. Figure which shows the detailed steps of the proposed range ambiguity suppression method.
Sensors 18 01454 g002
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.
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.
Sensors 18 01454 g003
Figure 4. GF-3 satellite images in full-polarization strip mode which acquired in the near Argun River of Inner Mongolia. (a) is the HH-polarization mode, (b) is the VV-polarization mode.
Figure 4. GF-3 satellite images in full-polarization strip mode which acquired in the near Argun River of Inner Mongolia. (a) is the HH-polarization mode, (b) is the VV-polarization mode.
Sensors 18 01454 g004
Figure 5. Image results of the ambiguous area: (a) represents the azimuth matched filter Rref with n = −1; (b) represents the azimuth matched filter Rref with n = 1.
Figure 5. Image results of the ambiguous area: (a) represents the azimuth matched filter Rref with n = −1; (b) represents the azimuth matched filter Rref with n = 1.
Sensors 18 01454 g005
Figure 6. (a) magnification details of the main ambiguous area image with n = −1 shown in Figure 5a; (b) shows the corresponding optical image obtained by latitude and longitude information of n = −1.
Figure 6. (a) magnification details of the main ambiguous area image with n = −1 shown in Figure 5a; (b) shows the corresponding optical image obtained by latitude and longitude information of n = −1.
Sensors 18 01454 g006
Figure 7. Detection results of ambiguous area targets based on CFAR. (a) represents the ambiguous area imaging results; (b) represents the detection results.
Figure 7. Detection results of ambiguous area targets based on CFAR. (a) represents the ambiguous area imaging results; (b) represents the detection results.
Sensors 18 01454 g007aSensors 18 01454 g007b
Figure 8. The range ambiguity suppression result of HH-polarization mode.
Figure 8. The range ambiguity suppression result of HH-polarization mode.
Sensors 18 01454 g008
Figure 9. The range ambiguity suppression result of VV-polarization mode.
Figure 9. The range ambiguity suppression result of VV-polarization mode.
Sensors 18 01454 g009
Table 1. GF-3 Satellite parameters.
Table 1. GF-3 Satellite parameters.
Radar ParameterValueRadar ParameterValue
λ0.055517 mPRF1292.0768
Transmit Band40 MHzLook angle38.91°
Sample rate66.667 MHzIncidence angle44.64°
Modulation rate1.6006 × 1012Near range1001.7 km
Satellite velocity(−1677.18, 5525.42, −4885.91)Reference range1015.3 km
Satellite position(−2,870,758.09, 3,815,169.12, 5,287,687.27)Center range1015.3 km
Target position118.4872° E, 49.2979° NTarget radius6,368,250 m
Equatorial radius6,378,140 mPolar radius6,356,755 m
equivalent velocity7097.4 m/sFdc6.508994
Table 2. Parameters used in CFAR.
Table 2. Parameters used in CFAR.
ParameterValue
Target window length2
Protection window length8
Background window length32
t12.5~3.5

Share and Cite

MDPI and ACS Style

Wen, X.; Qiu, X.; Han, B.; Ding, C.; Lei, B.; Chen, Q. A Range Ambiguity Suppression Processing Method for Spaceborne SAR with Up and Down Chirp Modulation. Sensors 2018, 18, 1454. https://doi.org/10.3390/s18051454

AMA Style

Wen X, Qiu X, Han B, Ding C, Lei B, Chen Q. A Range Ambiguity Suppression Processing Method for Spaceborne SAR with Up and Down Chirp Modulation. Sensors. 2018; 18(5):1454. https://doi.org/10.3390/s18051454

Chicago/Turabian Style

Wen, Xuejiao, Xiaolan Qiu, Bing Han, Chibiao Ding, Bin Lei, and Qi Chen. 2018. "A Range Ambiguity Suppression Processing Method for Spaceborne SAR with Up and Down Chirp Modulation" Sensors 18, no. 5: 1454. https://doi.org/10.3390/s18051454

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