Next Article in Journal
Joint Center Estimation Using Single-Frame Optimization: Part 2: Experimentation
Next Article in Special Issue
Design and Implementation of a Novel Polarimetric Active Radar Calibrator for Gaofen-3 SAR
Previous Article in Journal
Self-Calibration Algorithm for a Pressure Sensor with a Real-Time Approach Based on an Artificial Neural Network
Previous Article in Special Issue
Geo-Positioning Accuracy Improvement of Multi-Mode GF-3 Satellite SAR Imagery Based on Error Sources Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Real-Time Imaging Algorithm Based on Sub-Aperture CS-Dechirp for GF3-SAR Data

1
National Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China
2
School of Physics and Optoelectronic Engineering, Xidian University, Xi’an 710071, China
3
College of Geomatics, Xi’an University of Science and Technology, Xi’an 710054, China
*
Author to whom correspondence should be addressed.
Submission received: 30 June 2018 / Revised: 29 July 2018 / Accepted: 2 August 2018 / Published: 5 August 2018
(This article belongs to the Special Issue First Experiences with Chinese Gaofen-3 SAR Sensor)

Abstract

:
Conventional synthetic aperture radar (SAR) imaging algorithms usually require a period of time to process data that is longer than the time it takes to record one synthetic aperture or that corresponding to an adequate azimuth resolution. That is to say, the real-time processing system is idle during the long data recording time and the utilization of computational resources is low. To deal with this problem, a real-time imaging algorithm based on sub-aperture chirp scaling dechirp (CS-dechirp) is proposed in this paper. With CS-dechirp, the sub-aperture data could be processed to form an image with relatively low resolution. Subsequently, a few low-resolution images are generated as longer azimuth data are recorded. At the stage of full-resolution image generation, a coherent combination method for the low-resolution complex-value images is developed. As the low-resolution complex-value images are coherently combined one by one, the resolution is gradually improved and the full-resolution image is finally obtained. The results of a simulation and real data from the GF3-SAR validate the effectiveness of the proposed algorithm.

1. Introduction

Spaceborne synthetic aperture radar (SAR) can perform two-dimensional high-resolution imaging of ground targets at a long distance and in all-weather and all-day conditions [1,2,3], which makes it a key method for real-time information acquisition. SAR has been widely used in battlefield reconnaissance, target identification, resource exploration, disaster detection, and many other important areas [4,5]. Real-time imaging processing is a key technology of spaceborne SAR in earth observation. To achieve real-time imaging, the processing system is required to produce the final image at the same time as data recording ends or with a relatively short delay [6]. Therefore, how to improve the real-time processing capability of spaceborne SAR is a key issue.
In ideal conditions, some classic algorithms such as the range Doppler algorithm (RDA) [7], chirp scaling algorithm (CSA) [8], range migration algorithm (RMA or omega-k algorithm) [9], and polar format algorithm (PFA) [10] can obtain well-focused images, which are usually used in SAR real-time processing. These algorithms usually process data in a time period longer than that required for one synthetic aperture or that corresponding to an adequate azimuth resolution [11,12]. That is to say, the real-time processing system is idle during the long data recording times and the utilization of computational resources is low. This problem is particularly obvious in spaceborne SAR systems due to their long radar range, high azimuth resolution, and thus long synthetic aperture length [13]. In this case, if the imaging processing starts after recording the azimuth data of a full-aperture time, the required data storage and computational load increases substantially [14,15], which undoubtedly increases the requirements for the hardware processing system of a spaceborne SAR and hinders real-time imaging processing.
To deal with this problem, a new sub-aperture approach for real-time SAR processing has been studied by Alberto Moreira [16]. This approach has shown good real-time performance and imaging results. The signal is divided into sub-apertures that are small enough, so the range cell migration in one sub-aperture datum can be ignored. However, the range cell migration is considered in sub-aperture signal stitching. The interpolation of samples is used for range cell migration correction in this approach. This correction is limited to migration within one sub-aperture, which is less than half of the range resolution. In addition, the overlap of the sub-apertures must be greater than 21% for sufficient attenuation of the grating lobes or paired echoes, which increases the computational load. As shown in the simulation, the grating lobes appear in the imaging results. In our opinion, this phenomenon occurs due to the overlapping sub-aperture data. On this basis, this paper proposes a real-time imaging algorithm based on sub-aperture chirp scaling dechirp (CS-dechirp), which can perform imaging processing while the data are recorded and thus without waiting for a full-aperture time. In the proposed method, the low-resolution complex-value image is formed by using sub-aperture data, which can be much shorter than a synthetic aperture length. Subsequently, a few low-resolution images are generated as the longer azimuth data is recorded. At the stage of full-resolution image generation, a coherent combination method for the low-resolution complex-value images is developed. As the low-resolution complex-value images are coherently combined one by one, the resolution is gradually improved and the full-resolution image is finally obtained. Compared with the study by Alberto Moreira [16], the method proposed in this paper does not require sub-aperture data to overlap, and the grating lobes did not appear in the simulation results. The range cell migration can be adequately corrected by the sub-aperture CSA without limitation. Since the processing system produces the image at the same time as recording the data, this pipeline structure is very suitable for the data stream in SAR real-time processing.
This paper is organized as follows. In Section 2, the sub-aperture signal model for spaceborne SAR is established. In Section 3, the flow chart and equation derivation of the sub-aperture CS-dechirp imaging algorithm are detailed. In Section 4, the point targets simulation and the imaging results of the real data from the GF3-SAR are used to validate the effectiveness of the proposed algorithm. In Section 5, the content of this paper is summarized and analyzed.

2. Sub-Aperture Signal Model

The geometry of the spaceborne SAR is depicted in Figure 1. The satellite moves along the trajectory from P k to P k + 1 for the k -th sub-aperture data recording, where t k is the center of the k -th sub-aperture data and t s u b is the slow time in the azimuth of the sub-aperture data. According to Figure 1, the instantaneous slant range between the point target Q ( R B , X ) located in the scene and the ideal moving satellite can be written as:
R ( t s u b ; R B ) = R B 2 + ( v t k + v t s u b X ) 2
where v is the equivalent speed of the spaceborne SAR, R B is the coordinate in the range, and X is the coordinate in the azimuth of the point target Q ( R B , X ) .
Supposing that the transmitted signal is a linear frequency modulation (LFM) signal s t ( t ) = a r ( t ) exp ( j π γ t 2 ) , where γ is the chirp rate of the LFM signal, t is the fast time in the range, and a r ( t ) is the window function of the LFM signal, which is a rectangular window [17], then the received sub-aperture signal of the point target Q ( R B , X ) can be written in the t t s u b domain as:
s ( t , t s u b ; R B ) = a r ( t 2 R ( t s u b ; R B ) c ) a a ( t k + t s u b ) × exp [ j π γ ( t 2 R ( t s u b ; R B ) c ) 2 ] exp [ j 4 π λ R ( t s u b ; R B ) ]
where a a ( · ) is the beam window function in the azimuth, c is the speed of light, and λ is the wavelength of the signal.
From Equations (1) and (2), the difference between the sub-aperture expression of the signal and the traditional full-aperture expression is only a sub-aperture center term, t k , which plays a very important role in image stitching, as detailed in Section 3.2.3.

3. Real-Time Imaging Algorithm Based on Sub-Aperture CS-Dechirp

The real-time imaging algorithm based on sub-aperture CS-dechirp can perform imaging while recording data. The data from each sub-aperture are processed by the CS-dechirp algorithm to obtain a low-resolution complex-value image; subsequently, a high-resolution complex-value image of all received data can be obtained by combining all the sub-aperture images coherently.

3.1. Description of Real-Time Imaging Algorithm

A sketch of the real-time imaging algorithm is shown in Figure 2. The GF3 satellite flies along the trajectory from P1 to Pn, experiencing n times of data recording. During data recording, the corresponding sub-aperture data are processed to form a low-resolution complex-value image, which is then coherently combined with the foregone images to form a new image with higher resolution. As the sub-aperture images are coherently combined one by one, the resolution is gradually improved until full resolution is achieved. This forms the pipeline structure of the sub-aperture data stream, in which the final imaging result can be gradually generated as sub-aperture data are continuously recorded, which is suitable for the GF3-SAR real-time imaging processing.
The flow chart of the real-time imaging algorithm based on sub-aperture CS-dechirp is shown in Figure 3. Data from each sub-aperture recorded by the GF3-SAR are processed by the CS-dechirp algorithm to obtain a low-resolution complex-value image; the high-resolution complex-value image of all received data can then be obtained by stitching together all the sub-aperture images.

3.2. Data Processing Based on CS-Dechirp Algorithm

From Equation (1), the slant range R ( t s u b ; R B ) changes with the slow time t s u b , and the migration amount varies with R B , which shows the space-variance characteristic of the range cell migration. Because the scaling processing in the CSA can adequately compensate the range cell migration with the space-variance characteristic and the CSA is suitable for hardware implementation [18], the real-time imaging algorithm takes the CS-dechirp algorithm as the core to process sub-aperture data.
According to the flow chart shown in Figure 3, the data from each sub-aperture go through the sub-aperture CSA module, azimuth dechirp module, and sub-aperture complex-value images stitching module to obtain a SAR image with higher resolution.

3.2.1. Sub-Aperture CSA

Supposing the sub-aperture Doppler frequency corresponding to t s u b is f s u b , the expression of the k -th sub-aperture echo signal shown in Equation (2) can be written in the t f s u b domain as:
s k ( t , f s u b ; R B ) = a r ( t 2 R ( f s u b ; R B ) c ) a a ( R B λ f s u b 2 v 2 1 ( f s u b / f a M ) 2 ) × exp [ j 2 π v R B f a M 2 f s u b 2 ] exp [ j 2 π f s u b ( X v t k ) ] × exp [ j π γ e ( f s u b ; R B ) ( t 2 R ( f s u b ; R B ) c ) 2 ]
where 1 γ e ( f s u b ; R B ) = 1 γ R B 2 λ sin 2 θ c 2 cos 3 θ , θ is the squint angle and f a M = 2 v / λ . From Equation (3), the difference between the sub-aperture expression of the signal and the traditional full-aperture expression is the second exponent term of t k , which shows the phase of the point target relative to t k in sub-aperture imaging.
The chirp scaling quadratic phase function for the range cell migration correction is:
H 1 ( t , f s u b ; R 0 ) = exp [ j π γ e ( f s u b ; R B ) a ( f s u b ) ( t 2 R ( f s u b ; R 0 ) c ) 2 ]
where a ( f s u b ) = 1 / 1 ( f s u b / f a M ) 2 1 is the CS factor, R ( f s u b ; R 0 ) is the slant range with f s u b as an argument, and R 0 is the reference distance between the scene center and the radar trajectory. Since γ e ( f s u b ; R B ) changes less within the observed scene, R B in γ e ( f s u b ; R B ) can be replaced by R 0 . In the following part, we use γ e ( f s u b ; R 0 ) to substitute γ e ( f s u b ; R B ) .
After multiplying Equations (4) with (3) in the t f s u b domain, the signal is transformed to the f r f s u b domain as:
s k ( f r , f s u b ; R B ) = a r [ f r γ e ( f s u b ; R 0 ) [ 1 + a ( f s u b ) ] ] a a ( R B λ f s u b 2 v 2 1 ( f s u b / f a M ) 2 ) × exp [ j π f r 2 γ e ( f s u b ; R 0 ) [ 1 + a ( f s u b ) ] ] exp [ j 4 π c [ R B + R 0 a ( f s u b ) ] f r ] × exp [ j 2 π v R B f a M 2 f s u b 2 ] exp [ j 2 π f s u b ( X v t k ) ]
where Θ Δ ( f s u b ; R B ) = 4 π c 2 γ e ( f s u b ; R 0 ) a ( f s u b ) [ 1 + a ( f s u b ) ] ( R B R 0 ) 2 is the residual phase due to the operation of the chirp scaling quadratic phase function. The first exponent term in Equation (5) is the modulation phase function in the range frequency domain. R 0 a ( f s u b ) in the second exponent term shows that the migration amount of all point targets in the observed scene becomes the same after the CS operation.
Subsequently, the phase function for the range compressing, secondary range compressing, and range cell migration correction can be constructed as:
H 21 ( f r , f s u b ; R 0 ) = exp [ j π f r 2 γ e ( f s u b ; R 0 ) [ 1 + a ( f s u b ) ] ] exp [ j 4 π R 0 a ( f s u b ) c f r ]
After multiplying Equations (6) with (5), the signal is transformed to the t f s u b domain as follows:
s k ( t , f s u b ; R B ) = sin c a r ( t 2 R B c ) a a ( R B λ f s u b 2 v 2 1 ( f s u b / f a M ) 2 ) × exp [ j 2 π v R B f a M 2 f s u b 2 ] exp [ j Θ Δ ( f s u b ; R B ) ] × exp [ j 2 π f s u b ( X v t k ) ]
Here, the CSA completes the range compressing and range cell migration correction. From Equation (7), the first exponent term is required to process matched filtering, and the second exponent term is the residual phase term that is required to compensate. Based on the above, the residual phase compensation function for Equation (7) can be written as:
H 22 ( t , f s u b ; R 0 ) = exp [ j Θ Δ ( f s u b ; R B ) ]
However, the matched filtering used for the sub-aperture signal of Equation (7) will cause image aliasing, which can be illustrated by the time-frequency distribution in Figure 4. To avoid this aliasing caused by the matched filtering, zero padding in the azimuth must be used for the traditional CSA, which undoubtedly increases the computation burden. Therefore, the azimuth dechirp can be used for azimuth compression. The azimuth dechirp can solve the problem of image aliasing without zero padding, which can be illustrated by the time-frequency distribution in Figure 5.

3.2.2. Azimuth Dechirp

Since the high-order term in Equation (7) leads to the residual phase in the azimuth dechirp, the high-order term must be converted to a quadratic phase function. The conversion function can be written as:
H 3 ( t , f s u b ; R B ) = exp [ j 2 π v R B f a M 2 f s u b 2 j π k a 0 f s u b 2 ]
where k a 0 = 2 v 2 / λ R 0 is the time-domain chirp rate of the signal.
After the compensation of Equation (8) and the conversion of Equation (9), the signal can be written as:
s k ( t , f s u b ; R B ) = sin c a r ( t 2 R B c ) a a ( R B λ f s u b 2 v 2 1 ( f s u b / f a M ) 2 ) × exp [ j π k a 0 f s u b 2 j 2 π ( X v t k ) f s u b ]
The signal of Equation (10) is transformed to the t t s u b domain as:
s k ( t , t s u b ; R B ) = sin c a r ( t 2 R B c ) a a ( t k + t s u b X v ) exp [ j π k a 0 ( t k + t s u b X v ) 2 ]
From Equation (11), there are only the first and quadratic terms in its phase. A quadratic phase function can be constructed to complete the azimuth dechirp as:
H 4 ( t , t s u b ) = exp ( j π k a 0 ( t k + t s u b ) 2 )
After the dechirp processing in Equation (12), the signal of Equation (11) becomes:
s k ( t , t s u b ; R B ) = sin c a r ( t 2 R B c ) a a ( t k + t s u b X v ) × exp [ j 2 π k a 0 X v ( t k + t s u b ) + j π k a 0 ( X v ) 2 ]

3.2.3. Sub-Aperture Complex-Value Image Stitching

The theoretical form of sub-aperture image stitching can be described as follows. For simplicity, Equation (13) can be rewritten as:
s k ( t , t s u b ; R B ) = sin c a r ( t 2 R B c ) a a ( t k + t s u b X v ) × exp [ j 2 π f d t s u b j 2 π f d t k + j π k a 0 ( X v ) 2 ]
where f d = k a 0 X / v .
According to Equation (14), the Fourier transformation (FT) must be performed on the azimuth signal to obtain the sub-aperture focusing results. However, the focusing position f d in the azimuth exceeds the range of (−PRF/2, PRF/2) (PRF, Pulse Repeating Frequency), resulting in azimuth-aliasing after the FT for the signal. Therefore, the FT cannot be directly performed. To deal with this problem, an azimuth compensation function for spectrum shifting is introduced:
H 5 ( t , t s u b ) = exp ( j 2 π f M t s u b )
where f M = k a 0 t k . This function ensures that the focusing positions of point targets are within (−PRF/2, PRF/2). Multiplying Equations (15) with (14), the signal becomes:
s k ( t , t s u b ; R B ) = sin c a r ( t 2 R B c ) a a ( t k + t s u b X v ) × exp [ j 2 π ( f d f M ) t s u b j 2 π f d t k + j π k a 0 ( X v ) 2 ]
Then, the signal is transformed to the t f s u b domain as:
s k ( t , f s u b ; R B ) = sin c a r ( t 2 R B c ) δ a ( f s u b + f d f M ) × exp ( j 2 π f d t k ) exp ( j π k a 0 ( X v ) 2 )
Because the focus position f M f d is smaller than PRF, the single in Equation (17) is non-azimuth-aliasing.
In addition, a coherent combination is required for sub-aperture images. Therefore, the phase coherence and simplicity of the combination also need to be considered when performing the above frequency compensation. Thus, a combination method based on the integer points of the azimuth frequency is proposed.
According to Equation (17), one can note that exp ( j 2 π f d t k ) is a residual phase, which varies with t k or the azimuth sub-aperture index k . To combine the sub-aperture complex-value images coherently, this phase term should be compensated, and the phase compensation function can be written as:
H 6 ( t , f s u b ) = exp ( j 2 π ( f M f s u b ) t k )
After the compensation of Equation (18), the SAR complex-value image of the k -th sub-aperture data can be obtained as:
s k ( t , f s u b ; R B ) = sin c a r ( t 2 R B c ) δ a ( f s u b + f d f M ) exp ( j π k a 0 ( X v ) 2 )

4. Simulation and Real Data Results

In order to illustrate the effectiveness of the proposed algorithm, the imaging results of a point targets simulation and the GF3-SAR real data, which are separately processed by the full-aperture standard CSA and the sub-aperture CS-dechirp algorithm, are used for comparative analysis.

4.1. Point Targets Simulation

The simulation parameters of the spaceborne SAR are shown in Table 1 and the lattice in the scene is 3 (range) × 5 (azimuth). The distance between point targets is 1500 m in the range and 625 m in the azimuth.
The imaging results with the standard full-aperture CSA and the sub-aperture CS-dechirp algorithm are shown in Figure 6. To show the process of sub-aperture data imaging and coherent sub-aperture stitching, the focusing results and corresponding azimuth spectrums of the point target in the red rectangle of Figure 6b are shown in Figure 7 and Figure 8. Figure 7 and Figure 8 show the change in focusing results and azimuth spectrums as the sub-aperture images are combined. In this simulation, the length of the sub-aperture is one-fifth of the full-aperture length. The focusing results of 1–5 sub-apertures are presented in Figure 7a–e. From Figure 7, one can note that the resolution of the point target is gradually improved as the sub-aperture synthesis amount increases.
Figure 8a–e, corresponding to Figure 7a–e, show the change in the azimuth frequency spectrum with sub-aperture combination. From Figure 8a–e, the bandwidths of the azimuth frequency spectrum are gradually increased as the sub-aperture synthesis amount increases. From Figure 7 and Figure 8, it can be seen that the sub-aperture image combination in our method is equivalent to the coherent combination of the azimuth spectrum.
The contours and azimuth profiles of the point targets marked by the red box in Figure 6a,b is show in Figure 9. From Figure 9, the focusing result processed by the sub-aperture CS-dechirp algorithm is consistent with the standard CSA, which proves that the proposed real-time imaging algorithm in this paper is feasible. Figure 10 shows the phase of the point target processed by the sub-aperture CS-dechirp algorithm in the range and azimuth, indicating the good coherence of the sub-aperture signal. Since the high-order phase is compensated by the sub-aperture CS-dechirp method and the sub-aperture overlapping is avoided, the grating lobes do not appear in our method. In other words, our method has better performance in grating lobe suppression than that in Reference [16].

4.2. GF3-SAR Data Results

The parameters of the GF3-SAR with stripmap mode are shown in Table 2. The imaging result based on the standard full-aperture CSA is shown in Figure 11a and that based on the sub-aperture CS-dechirp algorithm is shown in Figure 11b. By enlarging the area marked by the red box in Figure 11a,b, detailed images are shown in Figure 12, Figure 13 and Figure 14. From the detailed images, one can note that the result processed by the sub-aperture CS-dechirp algorithm is consistent with the standard full-aperture CSA.

4.3. Computational Load Analysis

Assuming the final image is of N r × N a pixels, where N r and N a represent the range and azimuth points, respectively, the computational complexities of the standard full-aperture CSA and the proposed method are analyzed as follows. The floating-point operations can be derived as 5 N r N a log 2 N r for FFT in range and 5 N r N a log 2 N r for FFT in azimuth. The complex multiplication needs 6 N r N a operations [19].
For the standard full-aperture CSA, the processing system processes all recorded data. The standard full-aperture CSA contains four FFTs and three complex multiplications. Then, its computational complexity can be achieved by:
S c s a = 10 N r N a log 2 N r + 10 N r N a log 2 N a + 18 N r N a
For the proposed method, the processing system processes data from one sub-aperture each time. The proposed method contains five FFTs and five complex multiplications. Assuming the sub-aperture number in the azimuth is N , the number of azimuth points of sub-aperture data is N a / N . Then, the computational complexity for one sub-aperture data can be achieved by:
S s u b = 10 N r N a N log 2 N r + 15 N r N a N log 2 N a / N + 30 N r N a N = 10 N N r N a log 2 N r + 15 N N r N a log 2 N a + 30 N N r N a 15 N N r N a log 2 N
Thus, the computational complexity of the proposed method is:
S p r o = N · S s u b = 10 N r N a log 2 N r + 15 N r N a log 2 N a + 30 N r N a 15 N r N a log 2 N
The magnitudes of S c s a and S p r o are similar, so the computational load of the proposed method is almost the same as that of the standard full-aperture CSA. However, the proposed method can perform imaging processing while recording data; therefore, the proposed method based on the pipeline structure has a relatively high computational efficiency.

5. Conclusions

The GF3-SAR real-time imaging algorithm based on sub-aperture CS-dechirp can perform imaging processing while recording data, which reduces the amount of data storage and computational load for the GF3-SAR system and has good real-time performance. In order to illustrate the effectiveness of the algorithm proposed in this paper, the imaging results of point targets simulation and GF3-SAR real data, which are separately processed, are used for comparative analysis. The algorithm based on sub-aperture CS-dechirp is very suitable for the real-time imaging processing of the GF3-SAR.

Author Contributions

All authors contributed extensively to the study presented in this manuscript. Conceptualization and formal analysis, G.-C.S.; Writing-Original Draft Preparation, Writing-Review & Editing, Y.L.; Investigation, J.Y.; Resources, L.G.; Supervision and validation, M.X. and S.W.

Funding

This work was supported by the National Key R&D Program of China under grant 2017YFC1405600 and the Fundamental Research Funds for the Central Universities under grant JB180213.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sigurd, H.; Felipe, Q.; Michelangelo, V.; Marwan, Y.; Gerhard, K.; Alberto, M. Tandem-L: A Technical Perspective on Future Spaceborne SAR Sensors for Earth Observation. IEEE Trans. Geosci. Remote Sens. 2018, 50, 4792–4807. [Google Scholar]
  2. Ren, Y.; Li, X.; Gao, G.; Thomas, E. Derivation of Sea Surface Tidal Current From Spaceborne SAR Constellation Data. IEEE Trans. Geosci. Remote Sens. 2017, 55, 3236–3247. [Google Scholar] [CrossRef]
  3. Felipe, Q.; Marwan, Y.; Gerhard, K.; Alberto, M. An Analytical Error Model for Spaceborne SAR Multichannel Azimuth Reconstruction. IEEE Geosci. Remote Sens. Lett. 2018, 15, 853–857. [Google Scholar]
  4. Ni, C.; Tan, H.; Liu, Q. VPX based spaceborne SAR real-time echo simulation. In Proceedings of the IEEE 6th International Symposium on MAPE, Shanghai, China, 28–30 October 2015; pp. 759–763. [Google Scholar]
  5. Moreira, A.; Scheiber, R.; Mittermayer, J.; Spielbauer, R. Real-time implementation of the extended chirp scaling algorithm for air- and spaceborne SAR-processing. In Proceedings of the International Geoscience and Remote Sensing Symposium (IGARSS), Firenze, Italy, 10–14 July 1995; pp. 2286–2288. [Google Scholar]
  6. Gu, C.; Chang, W.; Li, X.; Jia, G.; Luan, X. A New Distortion Correction Method for FMCW SAR Real-Time Imaging. IEEE Geosci. Remote Sens. Lett. 2017, 14, 429–433. [Google Scholar] [CrossRef]
  7. Smith, A. A new approach to range-Doppler SAR processing. Int. J. Remote Sens. 1991, 12, 235–251. [Google Scholar] [CrossRef]
  8. Keith, R.; Runge, H.; Bamler, R.; Cumming, G.; Wong, H. Precision SAR processing using chirp scaling. IEEE Trans. Geosci. Remote Sens. 1994, 32, 786–799. [Google Scholar]
  9. Cafforio, C.; Prati, C.; Rocca, F. SAR data focussing using seismic migration techniques. IEEE Trans. Aerosp. Electron. Syst. 1991, 27, 199–207. [Google Scholar] [CrossRef]
  10. Zhu, D.; Ye, S.; Zhu, Z. A polar format algorithm using chirp scaling for spotlight SAR image formation. IEEE Trans. Aerosp. Electron. Syst. 2008, 44, 1433–1448. [Google Scholar] [CrossRef]
  11. Sun, G.C.; Xing, M.; Xia, X.G.; Wu, Y.; Huang, P.; Wu, Y.; Bao, Z. Multichannel Full-Aperture Azimuth Processing for Beam Steering SAR. IEEE Trans. Geosci. Remote Sens. 2013, 51, 4761–4778. [Google Scholar] [CrossRef]
  12. Sun, G.C.; Xing, M.; Xia, X.G.; Yang, J.; Wu, Y.; Bao, Z. A Unified Focusing Algorithm for Several Modes of SAR Based on FrFT. IEEE Trans. Geosci. Remote Sens. 2013, 51, 3139–3155. [Google Scholar] [CrossRef]
  13. Anghel, A.; Vasile, G.; Ioana, C.; Cacoveanu, R.; Ciochina, S. Micro-Doppler Reconstruction in Spaceborne SAR Images Using Azimuth Time–Frequency Tracking of the Phase History. IEEE Trans. Geosci. Remote Sens. 2016, 13, 604–608. [Google Scholar] [CrossRef]
  14. Zuo, S.S.; Xing, M.; Xia, X.G.; Sun, G.C. Improved Signal Reconstruction Algorithm for Multichannel SAR Based on the Doppler Spectrum Estimation. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2017, 10, 1425–1442. [Google Scholar] [CrossRef]
  15. Sun, G.C.; Xing, M.; Wang, Y.; Wu, Y.; Wu, Y.; Bao, Z. Sliding Spotlight and TOPS SAR Data Processing Without Subaperture. IEEE Geosci. Remote Sens. Lett. 2011, 8, 1036–1040. [Google Scholar] [CrossRef]
  16. Moreira, A. Real-time synthetic aperture radar (SAR) processing with a new subaperture approach. IEEE Trans. Geosci. Remote Sens. 1992, 30, 714–722. [Google Scholar] [CrossRef]
  17. Wang, X.; Wang, R.; Deng, Y.; Wang, W.; Li, N. SAR Signal Recovery and Reconstruction in Staggered Mode with Low Oversampling Factors. IEEE Geosci. Remote Sens. Lett. 2018, 15, 704–708. [Google Scholar] [CrossRef]
  18. Yuan, Y.; Chen, S.; Zhang, S.; Zhao, H. A chirp scaling algorithm for forward-looking linear-array SAR with constant acceleration. IEEE Geosci. Remote Sens. Lett. 2018, 15, 88–91. [Google Scholar] [CrossRef]
  19. Bie, B.; Xing, M.; Xia, X.; Sun, G. A frequency domain backprojection algorithm based on local cartesian coordinate and subregion range migration correction for high-squint SAR mounted on maneuvering platforms. IEEE Trans. Geosci. Remote Sens. 2018, 99, 1–16. [Google Scholar] [CrossRef]
Figure 1. Geometry of the spaceborne synthetic aperture radar (SAR).
Figure 1. Geometry of the spaceborne synthetic aperture radar (SAR).
Sensors 18 02562 g001
Figure 2. Sketch of the sub-aperture real-time imaging algorithm.
Figure 2. Sketch of the sub-aperture real-time imaging algorithm.
Sensors 18 02562 g002
Figure 3. Flow chart of the real-time imaging algorithm.
Figure 3. Flow chart of the real-time imaging algorithm.
Sensors 18 02562 g003
Figure 4. Time-frequency distribution for matched filtering: (a) unprocessed sub-aperture signal; (b) sub-aperture signal after matched filtering.
Figure 4. Time-frequency distribution for matched filtering: (a) unprocessed sub-aperture signal; (b) sub-aperture signal after matched filtering.
Sensors 18 02562 g004
Figure 5. Time-frequency distribution for azimuth dechirp: (a) unprocessed sub-aperture signal; (b) sub-aperture signal after azimuth dechirp.
Figure 5. Time-frequency distribution for azimuth dechirp: (a) unprocessed sub-aperture signal; (b) sub-aperture signal after azimuth dechirp.
Sensors 18 02562 g005
Figure 6. Imaging results of point targets: (a) standard full-aperture chirp scaling algorithm (CSA); (b) sub-aperture chirp scaling (CS)-dechirp algorithm.
Figure 6. Imaging results of point targets: (a) standard full-aperture chirp scaling algorithm (CSA); (b) sub-aperture chirp scaling (CS)-dechirp algorithm.
Sensors 18 02562 g006
Figure 7. Imaging results after sub-aperture synthesis. (a) With one sub-aperture datum; (b) With two sub-aperture data; (c) With three sub-aperture data; (d) With four sub-aperture data; (e) With five sub-aperture data.
Figure 7. Imaging results after sub-aperture synthesis. (a) With one sub-aperture datum; (b) With two sub-aperture data; (c) With three sub-aperture data; (d) With four sub-aperture data; (e) With five sub-aperture data.
Sensors 18 02562 g007aSensors 18 02562 g007b
Figure 8. Frequency spectrum of the azimuth signal after sub-aperture synthesis. (a) With one sub-aperture datum; (b) With two sub-aperture data; (c) With three sub-aperture data; (d) With four sub-aperture data; (e) With five sub-aperture data.
Figure 8. Frequency spectrum of the azimuth signal after sub-aperture synthesis. (a) With one sub-aperture datum; (b) With two sub-aperture data; (c) With three sub-aperture data; (d) With four sub-aperture data; (e) With five sub-aperture data.
Sensors 18 02562 g008aSensors 18 02562 g008b
Figure 9. Results of the standard full-aperture CSA and the proposed method: (a) standard CSA; (b) proposed method; (c) azimuth profiles based on the standard CSA; (d) azimuth profiles based on the proposed method.
Figure 9. Results of the standard full-aperture CSA and the proposed method: (a) standard CSA; (b) proposed method; (c) azimuth profiles based on the standard CSA; (d) azimuth profiles based on the proposed method.
Sensors 18 02562 g009aSensors 18 02562 g009b
Figure 10. Phase of the point target processed by the proposed algorithm: (a) phase in range; (b) phase in azimuth.
Figure 10. Phase of the point target processed by the proposed algorithm: (a) phase in range; (b) phase in azimuth.
Sensors 18 02562 g010
Figure 11. Imaging results of the GF3-SAR stripmap data (from top to bottom is range, and from left to right is azimuth): (a) standard full-aperture CSA; (b) proposed algorithm.
Figure 11. Imaging results of the GF3-SAR stripmap data (from top to bottom is range, and from left to right is azimuth): (a) standard full-aperture CSA; (b) proposed algorithm.
Sensors 18 02562 g011
Figure 12. Results of area A (from top to bottom is range, and from left to right is azimuth): (a) standard full-aperture CSA; (b) proposed algorithm.
Figure 12. Results of area A (from top to bottom is range, and from left to right is azimuth): (a) standard full-aperture CSA; (b) proposed algorithm.
Sensors 18 02562 g012
Figure 13. Results of area B (from top to bottom is range, and from left to right is azimuth): (a) standard full-aperture CSA; (b) proposed algorithm.
Figure 13. Results of area B (from top to bottom is range, and from left to right is azimuth): (a) standard full-aperture CSA; (b) proposed algorithm.
Sensors 18 02562 g013
Figure 14. Results of area C (from top to bottom is range, and from left to right is azimuth): (a) standard full-aperture CSA; (b) proposed algorithm.
Figure 14. Results of area C (from top to bottom is range, and from left to right is azimuth): (a) standard full-aperture CSA; (b) proposed algorithm.
Sensors 18 02562 g014
Table 1. Simulation parameters of the spaceborne SAR.
Table 1. Simulation parameters of the spaceborne SAR.
ParameterValue
Carrier frequency9.63 GHz
Bandwidth50 MHz
Sample frequency60 MHz
Velocity7391 m/s
PRF2738 Hz
Center line distance617 km
Table 2. Parameters of the GF3-SAR in stripmap mode.
Table 2. Parameters of the GF3-SAR in stripmap mode.
ParameterValue
Bandwidth100 MHz
Sample frequency133 MHz
Wavelength0.055 m
Velocity7132 m/s
PRF2580 Hz
Center line distance842 km

Share and Cite

MDPI and ACS Style

Sun, G.-C.; Liu, Y.; Xing, M.; Wang, S.; Guo, L.; Yang, J. A Real-Time Imaging Algorithm Based on Sub-Aperture CS-Dechirp for GF3-SAR Data. Sensors 2018, 18, 2562. https://0-doi-org.brum.beds.ac.uk/10.3390/s18082562

AMA Style

Sun G-C, Liu Y, Xing M, Wang S, Guo L, Yang J. A Real-Time Imaging Algorithm Based on Sub-Aperture CS-Dechirp for GF3-SAR Data. Sensors. 2018; 18(8):2562. https://0-doi-org.brum.beds.ac.uk/10.3390/s18082562

Chicago/Turabian Style

Sun, Guang-Cai, Yanbin Liu, Mengdao Xing, Shiyu Wang, Liang Guo, and Jun Yang. 2018. "A Real-Time Imaging Algorithm Based on Sub-Aperture CS-Dechirp for GF3-SAR Data" Sensors 18, no. 8: 2562. https://0-doi-org.brum.beds.ac.uk/10.3390/s18082562

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