Next Article in Journal
Effectiveness of Semi-Supervised Learning and Multi-Source Data in Detailed Urban Landuse Mapping with a Few Labeled Samples
Next Article in Special Issue
Time Series Surface Deformation of Changbaishan Volcano Based on Sentinel-1B SAR Data and Its Geological Significance
Previous Article in Journal
Performance Analysis of Zero-Difference GPS L1/L2/L5 and Galileo E1/E5a/E5b/E6 Point Positioning Using CNES Uncombined Bias Products
Previous Article in Special Issue
Entropy-Based Non-Local Means Filter for Single-Look SAR Speckle Reduction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel High-Squint Spotlight SAR Raw Data Simulation Scheme in 2-D Frequency Domain

1
College of Computer and Information Engineering, Henan University, Kaifeng 475004, China
2
Henan Engineering Research Center of Intelligent Technology and Application, Henan University, Kaifeng 475004, China
3
Henan Key Laboratory of Big Data Analysis and Processing, Henan University, Kaifeng 475004, China
4
School of Software, Henan University, Kaifeng 475004, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(3), 651; https://doi.org/10.3390/rs14030651
Submission received: 27 December 2021 / Revised: 26 January 2022 / Accepted: 28 January 2022 / Published: 29 January 2022
(This article belongs to the Special Issue Radar Imaging Theory, Techniques, and Applications II)

Abstract

:
Raw data simulation is the front-end work of synthetic aperture radar (SAR), which is of great significance. For high-squint spotlight SAR, the frequency domain simulation algorithm is invalid because of the range-azimuth coupling effect. In order to realize high-squint spotlight SAR raw data simulation in the frequency domain, an algorithm based on coordinate transformation and non-uniform fast Fourier transform (NUFFT) is proposed. This algorithm generates broadside raw data using a two-dimensional (2-D) frequency simulation algorithm; then, coordinate transformation is used by analyzing the characteristics of broadside and high-squint spotlight SAR. After coordinate transformation, NUFFT is carried out to realize the coupling relation in the 2-D frequency domain. Since the coordinate transformation ignores the influence of range walk, the range walk is compensated after NUFFT. As a result, compared with the traditional squint spotlight SAR frequency domain simulation algorithm, the proposed algorithm can improve the accuracy of point and distributed target imaging results, and the efficiency of the proposed algorithm can be significantly improved in contrast the traditional time domain algorithm.

Graphical Abstract

1. Introduction

Synthetic aperture radar (SAR) is an active imaging technology, based on microwaves, which has the ability to observe the Earth, all day and during all weathers. SAR can be operated in a number of different modes, such as stripmap, spotlight, scan and terrain observation by progressive scan (TOPS) [1,2,3,4,5,6,7]. Each mode has its own strengths and weaknesses. For the spotlight SAR mode, it steers the beam over the full acquisition time to illuminate the same area (a spot on the ground), so as to improve the azimuth resolution. Due to the mobility of SAR platforms, a high or low squint angle always exists [8,9]. Therefore, squint spotlight SAR has been widely used in high-resolution Earth-observation missions.
With the further development of SAR technology, raw data simulation has become an important research direction. Raw data simulation is helpful to verify the performance of a SAR system, and the system parameters can be optimized after simulation [10,11,12,13,14,15,16]. The time domain algorithm and the frequency domain algorithm are two ways of performing SAR raw data simulations. The high computational complexity of the time domain algorithm makes it difficult to simulate large scenes. Some researchers have proposed several time domain raw data simulation acceleration algorithms. Fu et al. used a graphics processing unit (GPU) to accelerate the time domain raw data simulation of high-squint spotlight SAR [17], and Zhang et al. also used GPU to accelerate raw data simulation of multiple-mode SAR [18,19,20]. However, these algorithms have high hardware requirements and the effect of acceleration is limited [21]. Therefore, it is necessary to research the frequency domain algorithm. In [22,23,24,25,26,27,28,29,30], Franceschetti et al. described the methods of frequency domain raw data simulation in particular. Such a simulator first builds a geometric model, and then, the scattering coefficient is calculated based on an asymptotic evaluation of SAR unit response. Finally, the data are processed in the two-dimensional (2-D) frequency domain to obtain the final raw data. The proposal of this direction is a significant milestone in the development of SAR raw data simulations. Subsequently, Franceschetti et al. proposed a method of spotlight SAR raw data simulation in the frequency domain. This method can simulate the range migration and the space-variant feature, andsimulation results verified the correctness of the algorithm. However, all the above methods mainly focus on the broadside mode. Wang et al. proposed a method to solve squint spotlight SAR raw data simulations [31]. It uses optical principles to efficiently calculate the interpolation results. Diao and Xu presented an efficient method for squint mode SAR raw data simulation, based on scaled Fourier transform (SCFT) in [32]. These methods, however, are only applicable to the case of small-squint spotlight SAR. In recent years, Li et al. proposed a method based on sub-aperture Keystone transform for squint SAR with a curved trajectory [33]. Although this method can accurately simulate the original raw data of extended scenes for high-squint SAR mounted on a maneuvering platform, its efficiency is much lower than 2-D frequency domain raw data simulation algorithms.
In order to precisely and efficiently calculate the high-squint raw data in the 2-D frequency domain, this work developed an algorithm based on coordinate transformation and non-uniform fast Fourier transform (NUFFT). First, due to the simplicity of the broadside SAR geometric model, calculating the broadside SAR raw data in the 2-D frequency domain was uncomplicated. The key to this algorithm was how to convert the raw data of broadside SAR to be equivalent to that of a high-squint SAR, and this problem can be solved by coordinate transformation. Then, in the process of coordinate transformation, interpolation was used to simulate the coupling of the azimuth and range. Because of the complexity of interpolation operations in data processing, we used the NUFFT instead of interpolation and fast Fourier transform (FFT). Finally, we simulated the linear range walk of high-squint SAR in the range frequency domain. This algorithm only required complex multiplications, and it was suitable for the characteristics of high-squint spotlight SAR raw data.
This article is organized as follows. In Section 2, the spotlight SAR geometric model and the broadside SAR raw data simulations in the frequency domain are introduced; more importantly, the coordinate transformation method and the application of NUFFT in this algorithm are described. In Section 3, some simulation and imaging results are shown. Moreover, the computational complexity of the algorithm is analyzed. In Section 4, we address relevant discussions and analyze the experimental results. Finally, some conclusions are reported in Section 5.

2. Models and Methods

2.1. Models

The geometries of broadside and high-squint spotlight SAR are shown in Figure 1, where the origin of the coordinate system is located at scene center O, and the X, Y and Z axes represent the azimuth, range, and geocentric directions, respectively. R c is the mid-range coordinate. R s t a r t and R e n d are the target-to-antenna distances at the beginning and end of SAR operations, r is the distance from the sensor to the center of the scene on the Zero-Doppler plane. As can be seen in Figure 1, the difference between the broadside and the high-squint is the look angle of the platform. However, the range-azimuth coupling effect is inegligible due to the squint angle. The signal model of the broadside and high-squint SAR are introduced in what follows.

2.1.1. Broadside SAR

In reference to Figure 1a, the slant range distance between the radar and the target (x, y) in the case of broadside SAR is:
R ( t a ) = ( v t a x ) 2 + ( R c y ) 2 2 ( R c y ) ( v t a x )
where v is the SAR sensor’s velocity and t a is the azimuth slow time. Then, the demodulated raw data can be expressed as:
S ( t , t a ) = r e c t ( t 2 R ( t a ) c T r ) · e x p ( j 4 π R ( t a ) λ ) · e x p ( j π K r ( t 2 R ( t a ) c ) 2 )
where t is the fast time, λ is the wavelength, Kr is the frequency modulation rate of the transmitted signal, Tr is the signal pulse width, and c is the propagation velocity of light.
The 2-D frequency spectrum of the scene raw signal is:
S ( f r , f a ) = d x d y σ ( x , y ) · H 0 ( f r , f a ) = W ( f r ) · e x p ( j π f r 2 K r ) · e x p ( 2 j π x v f a ) × e x p ( 4 j π y ( f c + f r c ) 2 ( f a 2 v ) 2 )
where σ ( x , y ) is the scattering coefficient of the scene and W ( · ) is the envelopment of the transmitted signal window function after Fourier transform. f r is the range frequency, f a is the azimuth frequency, and f c is the carried frequency.
A flow chart of the frequency domain raw data simulation algorithm is shown in Figure 2. First, 2-D Fourier transform is performed according to scattering coefficient σ ( x , y ) , and then the raw data are calculated based on Equation (3). Finally, we determine the inverse Fourier transform using the data.

2.1.2. High-Squint SAR

With reference to Figure 1b, the slant range distance between the radar and the target ( x , y ) in the case of high-squint SAR is:
R ( t a ) = ( v t a x ) 2 + ( R c y ) 2 2 ( R c y ) ( v t a x ) s i n θ
where θ is the center squint angle. In order to better analyze the distance characteristics in the case of squint, Equation (4) is expanded using the Taylor series at the central illuminated time:
R ( t a ) = R ( t a ) v · s i n θ · ( t a t 0 ) + 1 2 · v 2 c o s 2 θ R ( t 0 ) ( t a t 0 ) 2 + 1 2 · v 3 s i n θ c o s 2 θ R 2 ( t 0 ) ( t a t 0 ) 3 + = R ( t a ) + X 1 + X 2 + X 3
Equation (5) consists of range walk, X 1 = v · s i n θ · ( t a t 0 ) , range curvature, X 2 = 1 2 · v 2 c o s 2 θ R ( t 0 ) ( t a t 0 ) 2 , and higher order terms, X 3 = 1 2 · v 3 s i n θ c o s 2 θ R 2 ( t 0 ) ( t a t 0 ) 3 + . The range walk is larger than the range curvature, and the range curvature is larger than the higher order terms. However, when the squint angle exists, the higher order terms cannot be ignored in the frequency domain raw data simulation. Therefore, the error increases rapidly with an increase in squint angle. Assuming that Δ ψ / y represents the error, we can obtain:
Δ ψ / y = ( R ( t a ) R ( t a ) ) / y
The frequency domain raw data simulation algorithm is suitable for large scene simulations in the case of limited time. However, in Figure 1a, there is a high-squint angle in the SAR system; according to Equations (2) and (5), the transfer function of the SAR system is related to the slant range distance. Moreover, Equation (6) reveals that when the squint angle increases, the error will increase exponentially; thus, the raw data have spatial variability, and cannot be directly calculated in a 2-D frequency. Therefore, the frequency domain raw data simulation algorithm is invalid when a squint angle exists.

2.2. Methods

In order to solve the problem that the frequency domain raw data simulation algorithm is invalid in the case of high-squint SAR, an algorithm based on coordinate transformation and NUFFT was proposed. Due to the differences in the geometric model between broadside and high-squint SAR, we first generated broadside raw data using 2-D frequency simulation algorithm. To simulate the characteristics of high-squint SAR, a new approach based on coordinate transformation and NUFFT is presented in this section.

2.2.1. Coordinate Transformation

Figure 3a represents the geometric diagram between the platform and target, A B denotes the motion trajectory array of the platform, and d is the distance between two adjacent array elements. However, Figure 3a shows a geometric diagram of broadside SAR. To get the raw data of the high-squint SAR, rotating the coordinate series is a new way to address the problem of frequency domain raw data simulation. In an equivalent geometric diagram, depicted in Figure 3b, A B is the motion trajectory array of the platform after coordinate transformation, and θ 0 is the projection angle. As such, the distance between two adjacent array elements is changed to d / c o s θ 0 . To explore the difference between motion trajectory arrays A B and A B , the derivation process is shown as follows.
Assuming that the incident angle is θ 1 , the signal received by different coordinate series can be expressed as:
S A B ( t , n ) = s ( t ) · e x p ( j φ 0 + j 4 π λ n d sin ( θ 1 θ 0 )
S A B ( t , n ) = s ( t ) · e x p ( j φ 0 + j 4 π λ n d cos θ 0 s i n θ 1 )
n [ ( N 1 ) 2 , ( N 1 ) 2 ]
where s ( t ) is the amplitude of the signal and φ 0 is the phase of the received signal. In Equation (9), N is the serial number of samples.
Assuming that θ 1 θ 0 , the following approximation can be used in this case:
s i n θ 1 s i n θ 0 sin ( θ 1 θ 0 ) cos θ 0
According to Equation (10), the signal in Equation (8) can be rewritten as:
S A B = S ( t ) · e x p ( j φ 0 + j 4 π λ n d cos θ 0 s i n θ 1 j 4 π λ n d cos θ 0 s i n θ 0 )
where j 4 π λ n d cos θ 0 s i n θ 0 is the order of the range walk, and this order is considered in Section 2.2.3. Then we can obtain:
S A B ( t , k ) = S ( t ) · e x p ( j φ 0 ) δ ( k 2 N d sin ( θ 1 θ 0 ) λ )
S A B ( t , k ) = S ( t ) · e x p ( j φ 0 ) δ ( k 2 d cos θ 0 N ( sin θ 1 sin θ 0 ) λ )
Equations (12) and (13) are obtained by taking the Fourier transform with respect to the azimuth from Equations (7) and (12). Additionally, the Doppler frequency is changed as follows:
f A B = d · P R F λ sin ( θ 1 θ 0 )
f A B = d · P R F cos θ 0 · λ ( sin θ 1 sin θ 0 )
where P R F is the pulse repetition frequency, f A B is the azimuth frequency without the range walk order of motion trajectory array A B , and f A B is the azimuth frequency of the motion trajectory array A B . According to Equations (12)–(15), we have:
S A B ( t , k ) = S ( t ) · e x p ( j φ 0 ) δ ( k 2 N f A B P R F )
S A B ( t , k ) = S ( t ) · e x p ( j φ 0 ) δ ( k 2 N f A B P R F )
Observing the equations above, it is obvious that the differences in Equations (16) and (17) are only between f A B and f A B . In order to eliminate the difference from motion trajectory arrays A B and A B , we only need to consider the difference between f A B and f A B . Let
f A B = f A B sin θ 1 sin θ 0 cos θ 0 sin ( θ 1 θ 0 ) = d · P R F λ sin ( θ 1 θ 0 )
According to Equation (17), we obtain:
sin θ 1 = λ f A B v + sin θ 0
Then, sin ( θ 1 θ 0 ) can be expressed as:
sin ( θ 1 θ 0 ) = sin θ 1 sin θ 0 cos θ 1 cos θ 0 = cos θ 0 ( λ f A B v + sin θ 0 ) sin θ 0 1 ( λ f A B v + sin θ 0 ) 2
Substituting (20) into (18), f A B can be rewritten as:
f A B = d · P R F λ ( cos θ 0 ( λ f A B v + sin θ 0 ) sin θ 1 ( λ f A B v + sin θ 0 ) 2 )
The relevance between f A B and f A B is shown in Equation (21); after simplifying this equation and the equivalent substitution, we can obtain the following equation according to the “go-stop-go” assumption.
f a cos θ = ( f a + 2 v sin θ 0 c ( f c + f r ) tan θ 0 ( 2 v c ( f c + f r ) ) 2 ( f a + 2 v sin θ 0 c ( f c + f r ) ) 2 = ( f a + f d c ) tan θ ( 2 v c ( f c + f r ) ) 2 ( f a + f d c ) 2 )
with
f d c = ( 2 v sin θ 0 / c ) ( f c + f r )
Equation (22) is the frequency after coordinate transformation. Substituting (22) into (3), we can obtain:
S ( f r , f a ) = W r ( f r ) × exp ( j π f r 2 γ ) exp ( 2 j π ( f a + f d c ) ( R 0 sin θ 0 + v t c ) v ) × e x p ( 4 j π R 0 cos θ 0 ( f c + f r c ) 2 ( f a + f d c 2 v cos θ 0 ) 2 )
The above equations describe the relevance of the azimuth frequency before and after coordinate transformation. Equation (24) is the 2-D frequency spectrum of the scene raw signal after coordinate transformation. The process from Equation (3) to Equation (24) requires interpolation, butmotion trajectory array A B is non-uniform. Among the interpolation methods, the great majority of methods are not suitable for this situation. Those that are generally used may reduce computational efficiency and even cause deviations. Therefore, it is important to use a suitable method.

2.2.2. NUFFT

NUFFT is not a standard fast Fourier transform algorithm, but it is suitable for interpolation calculations of non-uniform data. According to Section 2.2.1, the array is non-uniform after coordinate transformation. In this case, NUFFT can not, only complete interpolation, but also complete azimuth Fourier transform. This operation is convenient to complete subsequent operations in the azimuth time domain, so the algorithm structure is optimized, and the efficiency is improved [34,35,36].
There are three kinds of NUFFTof interest: type-1 is non-equispaced data (NED) NUFFT, in which the spatial samples are located irregularly, and the spectral samples regular. Type-2 is non-equispaced result (NER) NUFFT, where the spatial samples are located regularly, and the spectral samples irregularly. In type-3, both the spatial and spectral samples are located irregularly. According to the process of coordinate transformation, NED NUFFT is suitable for this approach.
The NED NUFFT of samples { Z l } l = 1 M , located at non-equispaced points, { x l } l = 1 M , and evaluated on an equispaced grid, is defined as:
Z ^ k = l = 1 M Z l e x p ( j 2 π x l k N ) ,   k = N 2 , ,   N 2 1
In Equation (25), if M = N ,   x l = l , the NUFFT will become a uniform Fourier transform. According to the above process, the NED NUFFT uses a Fourier series, which is truncated, to approximate signal e x p ( j ξ x ) with Gaussian kernel e x p ( b x 2 ) .
φ ( x ) = e x p ( b x 2 ) e x p ( j ξ x )
where b is a coefficient of the Gaussian kernel. After derivation, the following equation is obtained:
Z ^ k = ( 2 π ) 1 / 2 ω ( 2 π k ξ N ) | m | K l = 1 M Z l ^ ( ξ x m ) e x p ( j 2 π m k / ξ N ) ,   k = N 2 , ,   N 2 1
where ω ( · ) is the window function, m is an independent variable regarding signal, and ξ can be interpreted as an oversampling factor. Generally, ξ = 2 . ^ represents the interpolation function.
According to Equations (26) and (27),functions of e x p ( j ξ x ) can be approximated by the finite order of φ ( x ) . The steps of the NUFFT algorithm according to [37] are:
(1)
The window function is used to process the raw data and make the data relatively smooth;
(2)
The oversampling technique is used to calculate the Fourier transform;
(3)
The result from the second step and the window function are used to perform convolution processing.
We use the above steps of coordinating transformation and NUFFT, to process the data, and then the signal can be expressed as:
S ( f r , t a ) = W r ( f r ) exp ( j π f r 2 K r ) × exp ( j 4 π c ( f c + f r ) ( R ( t a ) + v s i n θ 0 t a ) )
From the derivation of the coordinate transformations, the order of the range walk is not considered in Equation (11). Thus, Equation (28) requires compensation of the range walk.

2.2.3. Compensation of Range Walk

The range walk is an essential characteristic in squint spotlight SAR. After coordinate transformation and NUFFT, the characteristics of squint spotlight SAR have been simulated, except for range walk. Therefore, range walk compensation will be carried out.
According to Equation (11), in Section 2.2.1, and Equation (28), linear range walk is not considered. The high-squint raw data were not completed. Therefore, it is important to compensate for the range walk of the signal at this time. Assuming that the compensation function is expressed as:
H 1 ( f r , t a ) = exp ( j 4 π v s i n θ 0 c ( f c + f r ) · t a )
Multiplying Equation (29) with Equation (28), the Fourier transform is then taken with respect to the fast time. The raw data in case of high-squint spotlight SAR are completed.
To summarize, the raw data simulation in the 2-D frequency algorithm based on coordinate transformation and NUFFT can accurately calculate the data in the case of high-squint. Comparing to the time domain algorithm, the calculated efficiency has been significantly improved. A flow chart of the 2-D frequency domain raw data simulation algorithm is shown in Figure 4.

3. Experiment Results

To verify the performance of the proposed algorithm, the point target and distributed target simulation experiments are carried out. The comparison experiments with the traditional time domain algorithm and traditional squint frequency algorithm [32] are shown in this section. Moreover, the computational complexity of the proposed algorithm is shown.

3.1. Simulation Result

The main simulation parameters are shown in Table 1.
In this experiment, a point target is arranged in the scene. Figure 5a,b shows the 2-D frequency spectrums before and after coordinate transformation, and Figure 5c shows the final 2-D frequency spectrum after compensation of the range walk. It is evident that the space-variant feature, which is observed in the red circles, is well simulated after NUFFT is applied. The simulated imaging results are shown in Figure 6. The point target is well focused in the case of a high squint.
It is widely known that the time domain algorithm has the highest precision, but it has a large computational complexity. Therefore, evaluating the performance of frequency domain algorithms is usually compared with the time domain algorithm. Three different raw data simulation algorithms are used in the same conditions, and different simulation results are imaged via a back-projection algorithm in Figure 6. The peak side lobe ratio (PSLR) and the integrated side lobe ratio (ISLR) are investigated to evaluate the focus quality, and the results are summarized in Table 2. In addition, the simulation times of different algorithms are shown in Table 3.
As can be seen from Figure 6, Table 2 and Table 3, for a the perspective of simulation quality, the time domain algorithm is not affected by the squint angle and always has high precision. The point target of the traditional squint frequency domain algorithm is well focused under a squint angle of 30°, but the values of PSLR and ISLR are slightly higher than those of the time domain algorithm. Under squint angles of 45° and 60°, the imaging results of the traditional squint frequency domain algorithm are poor, and the values of PSLR and ISLR are significantly higher than those of the time domain algorithm and the proposed algorithm. Therefore, the traditional frequency domain algorithm is not suitable for simulations in cases of high squint. However, the quality of the proposed algorithm is high whatever the squint angle is 30°, 45°, or 60°. Although the proposed algorithm is not as accurate as the time domain algorithm, its imaging result of the point target is well focused. The values of PSLR and ISLR are almost the same as those of the time domain algorithm. In particular, when the squint angle increases, the PSLR and ISLR increase slightly.
From the perspective of simulation efficiency, obviously, the computational efficiency of the frequency domain algorithm is higher than that of the time domain algorithm. The time domain algorithm calculates the raw data of each target in a scene, and this is the reason why the time domain algorithm has a high accuracy and low efficiency. The traditional squint frequency domain algorithm has almost an equivalent efficiency to the proposed algorithm; however, with the increase in squint angle, the computing time of the proposed algorithm is shorter than the computing time of the traditional squint frequency domain algorithm, because the NUFFT of the proposed algorithm saves more time than the interpolation operation of the traditional frequency domain algorithm when the squint angle increases. The computational complexity will be analyzed in Section 3.2.
The phase error of the proposed algorithm compared to the time domain algorithm increases as the squint angle increases, and the phase error is at its maximum under a squint of 60°. However, the phase errors of the azimuth and the range are less than 20°, the phase errors are shown in Figure 7, the imaging result will not be affected. Therefore, the accuracy of proposed algorithm is high.
The input image and imaging results of the distributed target simulation after geometric correction are shown in Figure 8, and the simulation parameters of the distributed target are shown in Table 4. From the results of the point target and the distributed target, the proposed algorithm can simulate the original image better than the traditional frequency domain simulation algorithm, and the imaging results of the proposed method are almost the same as those of the input image. More importantly, the proposed algorithm only takes 21 s to complete a distributed target simulation. Compared with the time domain algorithm, the computational efficiency is improved by an order of magnitude.

3.2. Performance Analysis

As noted in the previous description, the entire proposed algorithm primarily includes complex multiplication and additions. The computational complexity is expressed through complex multiplication. Assuming that the size of the raw data is N a × N r , N a is the number of azimuth samples and N r is the number of range samples. For simplicity, the scene size is also set to N a × N r . According to Equation (2), the computational complexity of the time domain algorithm can be expressed as:
N 1 ( N a · N r ) 2
According to Figure 2, compared with the broadside frequency domain algorithm, the traditional squint frequency domain algorithm requires interpolation operations. Assuming that the computational complexity of interpolation is N i , the computational complexity of traditional squint frequency domain algorithm can be expressed as:
N 2 N i + 3 N a · N r + 2 N a · N r log 2 N a + 2 N a · N r log 2 N r
According to Section 2, assuming that the computational complexity of NUFFT is N n , the computational complexity of the proposed algorithm is expressed as:
N 3 N n + 3 N a · N r + 2 N a · N r log 2 N a + N a · N r log 2 N r
It can be seen from Equations (30)–(32) that the proposed algorithm has the lowest complexity, followed by the traditional frequency domain algorithm; the time domain algorithm is the most complicated. Therefore, if raw data simulation seeks a high accuracy and there is no limitation to efficiency, the time domain algorithm should be selected. However, if a mission needs to simulate an image quickly and accurately in the case of high squint, the proposed algorithm is undoubtedly the best choice.

4. Discussion

From the above experiment results, we discuss the proposed method from the following three aspects: hardware requirement, accuracy, and efficiency.
First, the proposed algorithm can quickly simulate raw data without high-performance hardware. The simulation platform and hardware information are listed in Table 5. It can be seen that the hardware requirements of the proposed algorithm are not strict. However, other time domain acceleration algorithms proposed by some researchers require a high-performance central processing unit (CPU) and graphics processing unit (GPU), for example, the algorithms in [17,18,19,20,21] require a powerful CPU and GPU for acceleration.
Second, in this paper, a back-projection algorithm is used to image the generated raw data. A back-projection algorithm is an accurate imaging algorithm, which is not affected by any factors. This imaging algorithm can accurately correspond to the information of raw data; therefore, in order to evaluate the accuracy of the generated raw data, the back-projection algorithm is used. We take the results from the BP algorithm as an evaluation standard to determine the accuracy of the raw data simulation. The accuracy of existing frequency domain raw data simulation algorithms is not satisfactory, especially in the case of high squint. The proposed algorithm is different from the traditional frequency domain simulation algorithm; we considered another approach, using the method of coordinate transformation. On the one hand, the proposed algorithm retains the high efficiency of the traditional frequency domain raw data simulation algorithm. On the other hand, through coordinate transformation and NUFFT, the accuracy is higher than that of the traditional frequency domain raw data simulation algorithm in the case of high squint.
Finally, in terms of efficiency, the traditional time domain algorithm has a very large amount of needed computation, and it calculates the raw data of each point in turn. However, the proposed algorithm in this paper adopts a unified processing method for all points through a transfer function in the frequency domain, and this computing mode is efficient. In addition, the computational complexity is shown in Section 3.2, and, compared with the time domain algorithm, the proposed algorithm has great advantages in terms of the amount of computation. In order to clearly show the operation time of the traditional squint frequency algorithm and the proposed algorithm, we created line charts, as shown in Figure 9, according to Table 3, and these line charts can more clearly show the changes in the trends of the computational efficiency. Due to the computational time of the time domain algorithm being much higher than that of the frequency domain algorithm, we do not show the time of the time domain algorithm on the line chart. As can be seen from Figure 9, the traditional squint frequency domain algorithm has almost the same efficiency as the proposed algorithm. With an increase in squint angle, the computing time of the proposed algorithm is shorter than that of the computing time of the traditional squint frequency domain algorithm. The proposed algorithm is more suitable for raw data simulation than the traditional squint frequency domain algorithm in the case of large squint.
In future work, we will continue to focus on how to improve the accuracy of the proposed algorithm and maintain the original computational efficiency. By referring to the imaging algorithm, we will continue to seek an algorithm that can improve the accuracy of raw data simulations. In addition, the recent deep automated image quality analysis method is also a potential direction.

5. Conclusions

High-squint SAR suffers from the range-azimuth coupling effect, and it is the reason why the frequency domain simulation algorithm is invalid when a squint angle exists. We analyzed the geometries of broadside and high-squint spotlight SAR. According to the different geometries, a frequency domain simulation algorithm based on coordinate transformation and NUFFT is proposed. NUFFT is used to solve problems after coordinate transformations. Comparative experiments show that the proposed algorithm simulates the characteristics of high-squint SAR, as well as point and distributed target imaging results and is well focused under a high squint angle. More importantly, the accuracy of the proposed algorithm is higher than that of the traditional algorithm. In addition, the computational efficiency of the proposed algorithm can be significantly improved, in contrast to the traditional time domain algorithm. In the future, relevant research will continue. We will focus on methods to reduce any artifacts in high-squint spotlight SAR simulated data generation; for example, recent deep automated image quality analysis methods [38] of the proposed algorithm show great potential, in the simulation of high-squint spotlight SAR raw data simulations.

Author Contributions

Analysis, Z.G. and Z.F.; methodology, N.L., Z.G. and Z.F.; validation, N.L., Z.F. and J.C.; resources, N.L., J.C., Z.G. and L.W.; writing—original draft preparation, Z.G. and Z.F.; writing—review and editing, N.L., L.W., Z.G. and J.C.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under grants 61871175 and 41901302, and in part by the Foundation of Key Laboratory of Radar Imaging and Microwave Photonics, Ministry of Education under grant RIMP2020003.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The study did not report any data.

Acknowledgments

The authors would like to thank the anonymous reviewers for their valuable and detailed comments, that are crucial in improving the quality of this paper. The authors would also like to thank all the teachers and students of the “SAR information processing” team for their help in this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, D.; Rodriguez-Cassola, M.; Prats-Iraola, P.; Wu, M.; Moreira, A. Reverse Backprojection Algorithm for the Accurate Generation of SAR Raw Data of Natural Scenes. IEEE Geosci. Remote Sens. Lett. 2017, 14, 2072–2076. [Google Scholar] [CrossRef] [Green Version]
  2. Zhang, Y.; Sun, J.; Lei, P.; Hong, W. SAR-Based Paired Echo Focusing and Suppression of Vibrating Targets. IEEE Trans. Geosci. Remote Sens. 2014, 52, 7593–7605. [Google Scholar] [CrossRef]
  3. Niu, S.; Qiu, X.; Lei, B.; Ding, C.; Fu, K. Parameter Extraction Based on Deep Neural Network for SAR Target Simulation. IEEE Trans. Geosci. Remote Sens. 2020, 58, 4901–4914. [Google Scholar] [CrossRef]
  4. Tong, X.; Bao, M.; Sun, G.; Han, L.; Zhang, Y.; Xing, M. Refocusing of Moving Ships in Squint SAR Images Based on Spectrum Orthogonalization. Remote Sens. 2021, 13, 2807. [Google Scholar] [CrossRef]
  5. Huai, Y.; Liang, Y.; Ding, J.; Xing, M.; Zeng, L.; Li, Z. An Inverse Extended Omega-K Algorithm for SAR Raw Data Simulation With Trajectory Deviations. IEEE Geosci. Remote Sens. Lett. 2016, 13, 826–830. [Google Scholar] [CrossRef]
  6. Li, N.; Zhang, H.; Zhao, J.; Wu, L.; Guo, Z. An Azimuth Signal-Reconstruction Method Based on Two-Step Projection Technology for Spaceborne Azimuth Multi-Channel High-Resolution and Wide-Swath SAR. Remote Sens. 2021, 13, 4988. [Google Scholar] [CrossRef]
  7. Li, N.; Shen, Q.; Wang, L.; Wang, Q.; Guo, Z.; Zhao, J. Optimal time selection for ISAR imaging of ship targets based on time-frequency analysis of multiple scatterers. IEEE Geosci. Remote Sens. Lett. 2021, 19, 1–5. [Google Scholar] [CrossRef]
  8. Xing, M.; Wu, Y.; Zhang, Y.; Sun, G.; Bao, Z. Azimuth Resampling Processing for Highly Squinted Synthetic Aperture Radar Imaging with Several Modes. IEEE Trans. Geosci. Remote Sens. 2014, 52, 4339–4352. [Google Scholar] [CrossRef]
  9. Fan, W.; Zhang, M.; Li, J.; Wei, P. Modified Range-Doppler Algorithm for High Squint SAR Echo Processing. IEEE Geosci. Remote Sens. Lett. 2019, 16, 422–426. [Google Scholar] [CrossRef]
  10. Chiang, C.-Y.; Chen, K.-S.; Yang, Y.; Zhang, Y.; Zhang, T. SAR Image Simulation of Complex Target including Multiple Scattering. Remote Sens. 2021, 13, 4854. [Google Scholar] [CrossRef]
  11. Jiang, W.; Wang, Y.; Li, Z.; Liu, G.; Zhang, M. Simulation of a Wideband Radar Echo of a Target on a Dynamic Sea Surface. Remote Sens. 2021, 13, 3186. [Google Scholar] [CrossRef]
  12. Ku, C.; Chen, K.; Chang, P.; Chang, Y. Imaging Simulation for Synthetic Aperture Radar: A Full-Wave Approach. Remote Sens. 2018, 10, 1404. [Google Scholar] [CrossRef] [Green Version]
  13. Xiao, P.; Liu, M.; Guo, W.; Chen, W. Reconstruction of Synthetic Aperture Radar Raw Data under Analog-To-Digital Converter Saturation Distortion for Large Dynamic Range Scenes. Remote Sens. 2019, 11, 1043. [Google Scholar] [CrossRef] [Green Version]
  14. Qian, Y.; Zhu, D. Image Formation of Azimuth Periodically Gapped SAR Raw Data with Complex Deconvolution. Remote Sens. 2019, 11, 2698. [Google Scholar] [CrossRef] [Green Version]
  15. Ma, B.; Yu, B.; Liu, C.; Wang, Y. A 2-D Fourier Domain Algorithm for StripMap FMCW-SAR Raw Signal Simulation. J. Electron. Inf. Technol. 2011, 33, 375–380. [Google Scholar] [CrossRef]
  16. Wei, L.; Li, S.; Wu, Y.; Xiang, M. Performance Comparison of Algorithms for SAR Raw Signal Generation. J. Electron. Inf. Technol. 2005, 27, 262–265. [Google Scholar]
  17. Fu, Z.; Bai, L.; Guo, Z.; Min, L.; Li, N. Research on Acceleration Algorithm for Raw Data Simulation of High Resolution Squint Spotlight SAR. In Proceedings of the 2021 IEEE International Geoscience and Remote Sensing Symposium (IGARSS 2021), Brussels, Belgium, 11–16 July 2021; pp. 3281–3284. [Google Scholar] [CrossRef]
  18. Zhang, F.; Hu, C.; Li, W.; Hu, W.; Wang, P.; Li, H. A Deep Collaborative Computing Based SAR Raw Data Simulation on Multiple CPU/GPU Platform. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2017, 10, 387–399. [Google Scholar] [CrossRef]
  19. Zhang, F.; Yao, X.; Tang, H.; Yin, Q.; Hu, Y.; Lei, B. Multiple Mode SAR Raw Data Simulation and Parallel Acceleration for Gaofen-3 Mission. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 2115–2126. [Google Scholar] [CrossRef]
  20. Zhang, F.; Hu, C.; Li, W.; Hu, W.; Wang, P. Accelerating Time-Domain SAR Raw Data Simulation for Large Areas Using Multi-GPUs. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 3956–3966. [Google Scholar] [CrossRef]
  21. Jing, G.; Zhang, Y.; Li, Z.; Sun, G.; Xing, M.; Bao, Z. Efficient Realiazation of SAR Echo Simulation Based on GPU. J. Syst. Eng. Electron. 2016, 38, 2493–2498. [Google Scholar] [CrossRef]
  22. Franceschetti, G.; Schirinzi, G. A SAR Processor Based on Two Dimensional FFT Codes. IEEE Trans. Aerosp. Electron. Syst. 1990, 26, 356–366. [Google Scholar] [CrossRef]
  23. Franceschetti, G.; Migliaccio, M.; Riccio, D.; Schirinzi, G. SARAS: A Synthetic Aperture Radar (SAR) Raw Signal Simulator. IEEE Trans. Geosci. Remote Sens. 1992, 30, 110–123. [Google Scholar] [CrossRef]
  24. Franceschetti, G.; Migliaccio, M.; Riccio, D. On ocean SAR raw signal simulation. IEEE Trans. Geosci. Remote Sens. 1998, 36, 84–100. [Google Scholar] [CrossRef]
  25. Franceschetti, G.; Lanari, R.; Marzouk, E. A New Two-Dimensional Squint Mode SAR Processor. IEEE Trans. Aerosp. Electron. Syst. 1996, 32, 854–863. [Google Scholar] [CrossRef]
  26. Franceschetti, G.; Lodice, A.; Riccio, D.; Ruello, G.; Siviero, R. SAR Raw Signal Simulation of Oil Slicks in Ocean Environments. IEEE Trans. Geosci. Remote Sens. 2002, 40, 1935–1949. [Google Scholar] [CrossRef]
  27. Franceschetti, G.; Lodice, A.; Riccio, D.; Ruello, G. SAR Raw Singal Simulation for Urban Structures. IEEE Trans. Geosci. Remote Sens. 2003, 41, 1986–1995. [Google Scholar] [CrossRef]
  28. Franceschetti, G.; Lodice, A.; Perna, S.; Riccio, D. Efficient Simulation of Airborne SAR Raw Data of Extended Scenes. IEEE Trans. Geosci. Remote Sens. 2006, 44, 2851–2860. [Google Scholar] [CrossRef]
  29. Franceschetti, G.; Guida, R.; Lodice, A.; Riccio, D.; Ruello, G. Efficient Simulation of Hybrid Stripmap/Spotlight SAR Raw Signals from Extended Scenes. IEEE Trans. Geosci. Remote Sens. 2006, 42, 2385–2396. [Google Scholar] [CrossRef]
  30. Franceschetti, G.; Lodice, A.; Perna, S.; Riccio, D. SAR Sensor Trajectory Deviations: Fourier Domain Formulation and Extended Scene Simulation of Raw Signal. IEEE Trans. Geosci. Remote Sens. 2006, 44, 2323–2334. [Google Scholar] [CrossRef]
  31. Wang, Y.; Zhang, Z.; Deng, Y. Squint Spotlight SAR Raw Signal Simulation in the Frequency Domain Using Optical Principles. IEEE Trans. Geosci. Remote Sens. 2008, 46, 2208–2215. [Google Scholar] [CrossRef]
  32. Diao, G.; Xu, X. Fast Algorithm for SAR Raw Data Simulation with Large Squint Angles. J. Electron. Inf. Technol. 2011, 33, 684–689. [Google Scholar] [CrossRef]
  33. Li, G.; Ma, Y.; Hou, J.; Xu, G. Sub-aperture Keystone Transform Based Echo Simulation Method for High-squint SAR with a Curve Trajectory. J. Electron. Inf. Technol. 2020, 42, 2261–2268. [Google Scholar] [CrossRef]
  34. Breglia, A.; Capozzoli, A.; Curcio, C.; Liseno, A. NUFFT-Based Interpolation in Backprojection Algorithms. IEEE Geosci. Remote Sens. Lett. 2021, 18, 2117–2121. [Google Scholar] [CrossRef]
  35. Hu, R.; Li, X.; Yeo, T.; Yang, Y.; Chi, C.; Zuo, F.; Hu, X.; Pi, Y. Refocusing and Zoom-In Polar Format Algorithm for Curvilinear Spotlight SAR Imaging on Arbitrary Region of Interest. IEEE Trans. Geosci. Remote Sens. 2019, 57, 7995–8010. [Google Scholar] [CrossRef]
  36. Zhao, S.; Deng, Y.; Wang, R. Imaging for High-Resolution Wide-Swath Spaceborne SAR Using Cubic Filtering and NUFFT Based on Circular Orbit Approximation. IEEE Trans. Geosci. Remote Sens. 2016, 55, 787–800. [Google Scholar] [CrossRef]
  37. Capozzoli, A.; Curio, C.; Liseno, A. Optimized Nonuniform FFTs and Their Application to Array Factor Computation. IEEE T Antenn Propag. 2019, 67, 3924–3938. [Google Scholar] [CrossRef]
  38. Mukherjee, S.; Valenzise, G.; Cheng, I. Potential of Deep Features for Opinion-Unaware, Distortion-Unaware, No-Reference Image Quality. Lect. Notes Comput. Sci. 2020, 2015, 87–95. [Google Scholar] [CrossRef]
Figure 1. Geometric models: (a) geometry of broadside spotlight SAR; (b) geometry of high-squint spotlight SAR.
Figure 1. Geometric models: (a) geometry of broadside spotlight SAR; (b) geometry of high-squint spotlight SAR.
Remotesensing 14 00651 g001
Figure 2. Flow chart of the frequency domain raw data simulation algorithm.
Figure 2. Flow chart of the frequency domain raw data simulation algorithm.
Remotesensing 14 00651 g002
Figure 3. Geometric diagram: (a) geometric diagram of broadside SAR; (b) geometric diagram after coordinate transformation.
Figure 3. Geometric diagram: (a) geometric diagram of broadside SAR; (b) geometric diagram after coordinate transformation.
Remotesensing 14 00651 g003
Figure 4. Flow chart of the 2-D frequency domain raw data simulation algorithm.
Figure 4. Flow chart of the 2-D frequency domain raw data simulation algorithm.
Remotesensing 14 00651 g004
Figure 5. The 2-D frequency spectrum: (a) before coordinate transformation and NUFFT; (b) after coordinate transformation and NUFFT; (c) after compensation of the range walk.
Figure 5. The 2-D frequency spectrum: (a) before coordinate transformation and NUFFT; (b) after coordinate transformation and NUFFT; (c) after compensation of the range walk.
Remotesensing 14 00651 g005
Figure 6. Imaging results of comparative experiment: (a) squint angle is 60°; (b) squint angle is 45°; (c) squint angle is 30°.
Figure 6. Imaging results of comparative experiment: (a) squint angle is 60°; (b) squint angle is 45°; (c) squint angle is 30°.
Remotesensing 14 00651 g006
Figure 7. The phase errors of the proposed method under the squint of 60°: (a) phase errors of the azimuth; (b) phase errors of the range.
Figure 7. The phase errors of the proposed method under the squint of 60°: (a) phase errors of the azimuth; (b) phase errors of the range.
Remotesensing 14 00651 g007
Figure 8. Imaging results of distributed target: (a) input image; (b) proposed algorithm; (c) traditional squint frequency domain algorithm.
Figure 8. Imaging results of distributed target: (a) input image; (b) proposed algorithm; (c) traditional squint frequency domain algorithm.
Remotesensing 14 00651 g008
Figure 9. Line chart of simulation times: (a) under a squint of 30°; (b) under a squint of 45°; (c) under a squint of 60°.
Figure 9. Line chart of simulation times: (a) under a squint of 30°; (b) under a squint of 45°; (c) under a squint of 60°.
Remotesensing 14 00651 g009
Table 1. Simulation parameters of the point target.
Table 1. Simulation parameters of the point target.
ParameterValue
Signal pulse width30 us
Pulse repetition frequency1500 Hz
Signal bandwidth50 MHz
Height10 km
Velocity200 m/s
Squint angle30°/45°/60°
Slant range of image center60 km
Center frequency9.65 GHz
Table 2. Comparison of evaluation results for simulated point target.
Table 2. Comparison of evaluation results for simulated point target.
Squint AngleAlgorithmsPSLR (dB)
Azimuth/Range
ISLR (dB)
Azimuth/Range
30°Time domain algorithm−13.27/−13.26−11.34/−11.39
Traditional squint algorithm−12.81/−12.63−11.08/−11.13
Proposed algorithm−13.16/−13.21−11.32/−11.28
45°Time domain algorithm−13.23/−13.29−11.70/−11.28
Traditional squint algorithm−10.34/−10.89−10.59/−10.38
Proposed algorithm−13.07/−13.13−11.46/−11.16
60°Time domain algorithm−13.25/−13.21−11.59/−11.32
Traditional squint algorithm−7.62/−6.82−5.89/−6.38
Proposed algorithm−12.89/−12.93−11.02/−10.87
Table 3. Simulation time comparisons.
Table 3. Simulation time comparisons.
Squint AngleAlgorithmsSimulation Times (s)
1 Target10 Targets100 Targets
30°Time domain algorithm3.2733.55329.68
Traditional squint frequency algorithm0.370.692.87
Proposed algorithm0.410.662.91
45°Time domain algorithm3.6735.75332.13
Traditional squint frequency algorithm0.430.772.93
Proposed algorithm0.390.682.88
60°Time domain algorithm3.8737.53339.53
Traditional squint frequency algorithm0.450.733.05
Proposed algorithm0.420.612.91
Table 4. Simulation parameters of the distributed target.
Table 4. Simulation parameters of the distributed target.
ParameterValue
Signal pulse width30 us
Pulse repetition frequency1300 Hz
Signal bandwidth60 MHz
Height5 km
Velocity400 m/s
Squint angle60°
Slant range of image center60 km
Center frequency9.65 GHz
Table 5. Simulation platform and hardware information.
Table 5. Simulation platform and hardware information.
PlatformInformation
CPUInter Xeon Gold 6126
Number of Core: 48
Clock Speed: 2.6 GHz
Memory128 GB
Computing softwareMATLAB 9.4 (R2018a)
Operating systemWindows 10
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, Z.; Fu, Z.; Chang, J.; Wu, L.; Li, N. A Novel High-Squint Spotlight SAR Raw Data Simulation Scheme in 2-D Frequency Domain. Remote Sens. 2022, 14, 651. https://doi.org/10.3390/rs14030651

AMA Style

Guo Z, Fu Z, Chang J, Wu L, Li N. A Novel High-Squint Spotlight SAR Raw Data Simulation Scheme in 2-D Frequency Domain. Remote Sensing. 2022; 14(3):651. https://doi.org/10.3390/rs14030651

Chicago/Turabian Style

Guo, Zhengwei, Zewen Fu, Jike Chang, Lin Wu, and Ning Li. 2022. "A Novel High-Squint Spotlight SAR Raw Data Simulation Scheme in 2-D Frequency Domain" Remote Sensing 14, no. 3: 651. https://doi.org/10.3390/rs14030651

APA Style

Guo, Z., Fu, Z., Chang, J., Wu, L., & Li, N. (2022). A Novel High-Squint Spotlight SAR Raw Data Simulation Scheme in 2-D Frequency Domain. Remote Sensing, 14(3), 651. https://doi.org/10.3390/rs14030651

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