Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Road Surface Classification Using a Deep Ensemble Network with Sensor Feature Selection
Previous Article in Journal
Magnetoconductivity and Terahertz Response of a HgCdTe Epitaxial Layer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sub-Nyquist SAR Based on Pseudo-Random Time-Space Modulation

1
School of Electronics and Information Engineering, Beihang University, Beijing 100083, China
2
Department of Applied Science and Frontier Technology, Qian Xuesen Laboratory of Space Technology, Beijing 100094, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(12), 4343; https://doi.org/10.3390/s18124343
Submission received: 20 October 2018 / Revised: 4 December 2018 / Accepted: 6 December 2018 / Published: 9 December 2018
(This article belongs to the Section Remote Sensors)

Abstract

:
Sub-Nyquist sampling technology can ease the conflict between high resolution and wide swath in a synthetic aperture radar (SAR) system. However, the existing sub-Nyquist SAR imposes a constraint on the type of the observed scene and can only reconstruct the scene with small sparsity (i.e., number of significant coefficients). The information channel model of microwave imaging radar based on information theory, in which scene, echo, and the mapping relation between the two correspond to information source, sink, and channel, is built, and noisy-channel coding theorem explains the reason for the aforementioned under this model. To allow the wider application of sub-Nyquist SAR, this paper proposes sub-Nyquist SAR based on pseudo-random space-time modulation. This modulation is the spatial and temporal phase modulation to the traditional SAR raw data and can increase the mutual information of information source and sink so that the scenes with large sparsity can be reconstructed. Simulations of scenes with different sparsity, e.g., an ocean with several ships and urban scenes, were run to verify the validity of our proposed method, and the results show that the scenes with large sparsity can be successfully reconstructed.

1. Introduction

Sub-Nyquist sampling technology [1,2,3,4,5,6,7], in which the sampling frequency is smaller than the Nyquist sampling rate [8], was proposed because there is redundancy in any information [9]. This technology can be applied in various fields, e.g., wireless communication [10] and synthetic aperture radar (SAR) [11,12,13,14,15,16,17]. Synthetic aperture radar takes samples of signals reflected from an illuminated scenario/scene to obtain a high-resolution image after processing at all times and in all weather conditions [18,19]. In a traditional SAR system, the pulse repetition frequency (PRF) should satisfy the Nyquist sampling theorem and exceed the Doppler bandwidth to avoid spectrum aliasing, and the Doppler bandwidth is inversely proportional to the azimuth resolution. This means that the better the resolution, the higher the PRF. Additionally, the echo must be completely received within one pulse repetition interval, i.e., the reciprocal of the PRF. This means that the larger the swath, the lower the PRF. Therefore, high resolution requiring high PRF and wide swath requiring low PRF are a contradiction [18,19]. When sub-Nyquist sampling is adopted on the azimuth dimension, the contradiction between high resolution and wide swath can be eased under compressive sensing (CS) theory.
The restricted isometry property (RIP) is a sufficient and necessary condition for successful reconstruction in CS theory [20]. The RIP implies that the sparser the scene is, the more easily the RIP is satisfied. However, approximately certifying the RIP is challenging [21]. In the SAR system, the RIP is usually satisfied by random sub-Nyquist sampling on the azimuth dimension, and this sampling method is achieved by randomly transmitting a smaller number of pulses or randomly selecting from the Nyquist samples [7,15]. Under this sub-Nyquist sampling method, the reconstructed scenes are required to be sparse according to CS theory [20,22], and most references [7,11,15,23] directly assume that the reconstructed scene is sparse when a small number of very strong reflectivity areas exist in a darker background. Because the azimuth reconstructed matrix is decided by the Doppler movement, SAR has the limitation of application on the type of the observed scene only by random sub-Nyquist sampling method. In this paper, we explain the reason for this constraint from the perspective of information theory and relieve this constraint to make the scene more widely applicable. Additionally, greedy algorithms [24,25], l1 relaxation methods [26,27,28,29], and Bayesian-based methods [30,31,32] can be used under the assumption of satisfying the RIP. This paper mainly adopts l1 relaxation methods with small recovery error compared with the other two algorithms [33].
Receiving the raw data from the observed scene is essentially the procedure of information transferring in the microwave imaging radar system [34,35]. Therefore, the information channel model of the microwave imaging radar is built to analyze the transferring information entropy from the perspective of information theory, and the observed scene, the raw data and the mapping relation between the two correspond to information source, sink, and channel. The carrier of transferring information is information channel. Channel capacity is a quantitative index which evaluates the quality of information channel. Under this model, noisy-channel coding theorem indicates that all information is nearly error-free transferred in one coded way, while the transferring information entropy is smaller than the channel capacity, otherwise reliable transferring cannot be achieved [36]. This theorem explains the necessary requirement of the reconstruction and the reason for the constraint on the sparsity of scene.
To relieve the constraint on the type of the observed scene in the existing sub-Nyquist SAR system so that the scene with large sparsity can be reconstructed, this paper proposes sub-Nyquist SAR based on pseudo-random space-time modulation. Under the information channel model of the microwave imaging radar based on information theory, this modulation is the spatial and temporal phase modulation to the traditional SAR raw data and can increase the mutual information of information source and sink so that more information about scene can be achieved. After the pseudo-random space-time modulation, the sub-Nyquist SAR system can also reconstruct the scene with large sparsity. This system lays the foundation to achieve a high-resolution and wide-swath SAR system.
The rest of this paper is organized as follows. Section 2 presents the constraint on the type of the observed scene in the existing sub-Nyquist SAR system. In Section 3, the information channel model to the microwave imaging radar based on information theory is established, and noisy-channel coding theorem is used to explain the reason for the constraint under this model. In Section 4, the pseudo-random space-time modulation is applied to the traditional SAR raw data to relieve this constraint, and the compressive reflector antenna (CRA) is taken as the carrier of this modulation. Section 5 proposes the sub-Nyquist SAR based on pseudo-random space-time modulation, which is described from four aspects: the choice of the sampling method, the echo signal model, the reconstruction method, and the reconstruction performance. The simulation results are presented in Section 6, and Section 7 concludes the paper.

2. SAR Image Reconstruction Based on CS

After removing carrier frequency through quadrature demodulation, the SAR raw data for point targets can be expressed as
s c ( τ , η ) = x i , y i σ i W i ( τ , η ) exp { j π K r [ τ 2 R i ( η ) c ] 2 j 4 π R i ( η ) λ } + n ( τ , η ) ,
where τ and η are the fast time along the range and the slow time along the azimuth, respectively, σ i and W i ( τ , η ) are the backscattering cross-section and the weighting pattern corresponding to the target at ( x i , y i ) , x i and y i denote the azimuth and range coordinates, respectively, K r denotes the frequency modulation rate of the pulse. R i ( η ) represents the distance between the radar and the point target at the slow time η . c is the speed of light, λ is the wavelength, and n ( τ , η ) denotes the system noise.
For moderating the contradiction between high-resolution and wide-swath, in this paper sub-Nyquist sampling is implemented along the azimuthal dimension in this paper. Accordingly, SAR image reconstruction includes three steps: range compression based on matched filtering, range cell migration correction (RCMC), and azimuth reconstruction from sub-Nyquist samples based on CS. After implementing range compression and RCMC to Equation (1), the signal at a certain range cell is represented by:
s c ( τ 0 , η ) = x i , y i σ i W i ( τ 0 , η ) T r × sin c { K r T r ( τ 0 2 R i ( η η c i ) c ) } × exp { j 4 π R i ( η η c i ) λ } + n ( τ 0 , η ) ,
where η c i is the beam center crossing time for the target ( x i , y i ) . Equation (2) can also be expressed as the vector-matrix form
s N × 1 = D N × M σ M × 1 + n N × 1 ,
where N is the number of samples on the azimuthal dimension, and M is the number of resolution cells at the certain range cell in the observed scene, s N × 1 = [ s c ( τ 0 , η 1 ) , s c ( τ 0 , η 2 ) , , s c ( τ 0 , η N ) ] T , σ M × 1 = [ σ 1 , σ 2 , , σ M ] T , D N × M = { D i ( τ 0 , η n ) } n = 1 , i = 1 N , M , D i ( τ 0 , η n ) = W i ( τ 0 , η n ) · T r · exp { j 4 π R i ( η n η c i ) / λ } , and n N × 1 = [ n ( τ 0 , η 1 ) , n ( τ 0 , η 2 ) , , n ( τ 0 , η N ) ] T .
To exactly reconstruct σ M × 1 from the undetermined system (Equation (3)) for which the unknown number M is larger than the known number N , the RIP of the matrix D N × M should be satisfied. Assuming that Λ { 1 , 2 , , M } is an index subset with S elements and 1 S M , the RIP is quantitatively denoted as [20] for all σ Λ S :
( 1 δ S ) σ Λ 2 2 D Λ σ Λ 2 2 ( 1 + δ S ) σ Λ 2 2 ,
where the restricted isometry constant δ S ( 0 , 1 ) is a small quantity, · 2 is l2-norm of vector [37], D Λ denotes the sub-matrix formed from the columns of D indexed by Λ , and σ Λ denotes the sub-vector formed from the rows of σ indexed by Λ . The RIP essentially explains that the distance between two signals D σ and D σ is proportional to the distance between the two signals σ and σ , it can guarantee the exact reconstruction. The sparsity S , the non-zero number of the unrecovered signal σ M × 1 , is introduced to express the RIP definition. The fewer the number of non-zero elements, the sparser the signal is and the more easily the RIP is satisfied. Additionally, the improvement of RIP condition contributes to enlarging the sparsity so that the reconstructed scene is not only limited to the sparse type.
The RIP implies that randomness plays a crucial role in the reconstructed matrix D N × M . Usually D N × M has the two following compositional forms [22]: (1) all variables of the matrix obey one certain random distribution; and (2) the rows of the matrix are randomly selected from an orthogonal basis. The azimuthal signal in the SAR system is decided by the Doppler movement between the radar and the scene so that the elements of D N × M are the Doppler signal, and D N × M under the Nyquist sampling method is an orthogonal basis. Therefore, if sub-Nyquist samples are randomly selected from the Nyquist samples, the scene may be recovered according to the Form 2. Under the assumption of satisfying the RIP, three categories of the reconstructed algorithm can be used. The first is greedy algorithms, e.g., orthogonal matching pursuit (OMP) [24] and thresholding [25]; the second is l1 relaxation methods, e.g., the Dantzig selector [26] and basis pursuit denoising (BPDN) [27,28,29]; and the third is Bayesian-based methods, which include maximum a posteriori (MAP) estimation and Hierarchical Bayesian framework [30,31,32]. The greedy method is faster than the other two methods in term of the recovery time, l1 relaxation method performs better in term of small recovery error, Bayesian method balances between small recovery error and short recovery time [33]. This paper mainly adopts l1 relaxation methods with small recovery error. The optimization equation of Equation (3) without matrix subscript is
σ ^ = arg min σ { s D σ 2 2 + γ σ 1 } ,
where γ is the regularization factor. The first item s D σ 2 2 guarantees the recovery accuracy, the second item σ 1 ensures the sparsity of the recovered scene, and γ balances between the recovery accuracy and the sparsity. The parameter γ , an experience value, is selected by minimizing the reconstructed error for each application.
Two TerraSAR-X images, i.e., ocean containing ships with 0.8% sparsity level, and an urban scene with 25% sparsity level, were selected. The sparsity level is percent of significant coefficients in the signal. Their reflectivity functions were used to simulate raw data in Figure 1, and the simulated parameters are shown in Table 1. The results in Figure 1a,c are reconstructed with Nyquist samples based on a chirp scaling algorithm [38]. The results shown in Figure 1b,d are reconstructed with sub-Nyquist samples obtained according to the abovementioned Form 2, and the reconstructed algorithm is based on the optimization equation of Equation (5), the corresponding parameters γ are about 100 and 1, respectively. From the simulation results, we can see that most details of ships on the ocean in Figure 1b are maintained, but the urban scene in Figure 1d is not successfully reconstructed. The explained reason that the reconstructed results with the same simulated parameters are different may be that the sparsity of urban scene does not satisfy the RIP condition so that it is not successfully reconstructed. It is challenging to quantitatively certify the RIP [21]. The simulated results also show that this sub-Nyquist SAR system has constraint on the scene type. The reason for the unsuccessful reconstruction is explained in Section 3.

3. Information Channel Model of Microwave Imaging Radar Based on Information Theory

Receiving the echo from the observed scene in the microwave imaging radar system is actually the procedure of information transferring [34,35]. Information transferring is a major field of information theory and transferred information can be quantitatively calculated by the concept of entropy. Scenes with different sparsity have different information entropy. The sparser the scene, the smaller the information entropy. To quantitatively analyze this procedure of information transferring, an information channel model of microwave imaging radar model based on information theory was built and the corresponding relationship between the two should be established.
The key factors in the radar imaging system are the scene, the echo, and the mapping relation between the scene and the echo. The mapping relation between the two is influenced by the slant range, the amplitude, phase of the antenna, etc. Analogous to the procedure of information transferring in information theory [39,40], the scene, the mapping relation, and the echo correspond to information source, channel, and sink, respectively, as illustrated in Figure 2. Under this model, if the scene is successfully reconstructed from the echo, it means that the information channel successfully transferred information in a coded way.
In information theory, the necessary condition of successfully transferring information can be described by noisy-channel coding theorem. This theorem indicates that the transferred information entropy must be smaller than the channel capacity so that all information is nearly error-free transferred to the receiver in a coded way, otherwise, reliable transferring cannot be achieved [40]. Channel capacity, a quantitative index of analyzing information channel, denotes the maximum transferring information entropy with arbitrarily small error probability [40]. Under the information channel model of microwave imaging radar, channel capacity is decided by the reconstructed matrix, and it explains the maximum information entropy of scene is obtained from echo. As done for multiple-input multiple-output (MIMO) communication system, the channel capacity in bits/s/Hz is [40].
C n = 1 N log 2 ( 1 + S N R × β n 2 ) ,
where S N R is the signal-to-noise ratio of the data after range compression for that imaging model based on information channel built on the azimuth dimension, and β n is the non-zero singular value of D N × M . Figure 3 illustrates variable curves of channel capacity against S N R at different sampling frequencies. The simulation result shows that channel capacities under the sub-Nyquist sampling method are smaller than those under the Nyquist sampling method, and that smaller sampling frequencies are associated with smaller channel capacities.
Corresponding to the noise-channel coding theorem, the necessary condition of the successful reconstruction is that the channel capacity should be larger than the information entropy of the scene. The reconstructed matrices of Figure 1a,c is the same, so the channel capacity corresponding to Figure 1a,c is also the same. Similarly, the channel capacity corresponding to Figure 1b,d is also the same. However, the information entropy of the scenes in Figure 1a,c is different, for the ocean with several ships is sparser than the urban scene, so that the information entropy of the urban scene is larger than that of the ocean with several ships. The same is true for the scenes in Figure 1b,d. Comparing the result in Figure 1b with that in Figure 1d, the reason for the result in Figure 1d is that the information entropy of the urban scene is larger than the channel capacity under the sub-Nyquist sampling method. Additionally, the simulation results also show that the channel capacity under the Nyquist sampling method is larger than that under the sub-Nyquist sampling method, comparing the result in Figure 1c with that in Figure 1d, as illustrated in Figure 3. Although the scene information entropy is not quantitatively calculated, the above explanation and inference to the results in Figure 1 are logical according to the simulated results.
Therefore, sub-Nyquist samples in the existing sub-Nyquist SAR system only reconstruct scenes with small sparsity for the reconstructed matrix with small channel capacity, and so the existing sub-Nyquist SAR system has great limitation. In the next section, this limitation will be relieved by pseudo-random space-time modulation.

4. Pseudo-Random Space-Time Modulation

For that the reconstructed matrices are decided by the Doppler movement between the moving platform and scene in the existing sub-Nyquist SAR system, the randomness required by the RIP to large value of sparsity may be not enough when the sub-Nyquist sampling rate is certain. Therefore, only changing the movement or sampling method does not improve the reconstruction to large value of sparsity. The reconstructed matrices with good RIP, e.g., random matrix, hold for large values of the sparsity [20,41]. If a scene with large sparsity, e.g., an urban scene, is to be reconstructed with sub-Nyquist samples, it adopts pseudo-random space-time modulation to form the random reconstructed matrix in this paper.

4.1. Pseudo-Random Space-Time Modulation

The RIP implies that randomness plays a major role in the reconstruction. Random measurement waveform or random sub-Nyquist sampling can both yield near-optimal reconstruction with overwhelmingly high probability [42]. The existence of randomness is achieved by pseudo-random modulation regardless of which sub-Nyquist sampling method is adopted. Additionally, the information of slant range is implied in the echo phase [19]. Therefore, this modulation is set to be phase modulation which generates random phase obeying a certain distribution along the azimuthal dimension so that the random reconstructed matrix is formed. This method is essentially the phase modulation, and we named it the pseudo-random space-time modulation in this paper. As illustrated in Figure 4, Information Channel, based on the Doppler movement, uses the modulated scene information as information source, and is echoed as information sink. The procedure involves transferring information in the SAR system.
After this modulation, multiplying the traditional raw data with the complex exponent with variable random phase on space and time generates the new raw data, and the random reconstructed matrix is formed.

4.2. Choice of Pseudo-Random Space-Time Modulation

For that Information Channel C is decided by the Doppler movement so that it is immutable at a certain sub-Nyquist sampling frequency, increasing the mutual information between information source X1 and information sink Y is to improve this modulation. In this modulation, the phase of complex exponent is a stochastic variable. The distribution of stochastic variable has many distributions, e.g., uniform distribution or Gaussian distribution, and different distributions lead to different reconstruction performance. While the ratio of the maximum to minimum non-zero singular values is small, i.e., βmax/βmin→1, it can transfer more scene information, where βmax and βmin are the maximum and minimum non-zero singular values of random modulation matrix, respectively. It can be explained by the restricted isometry constant δs. The RIP can be simplified as (1 − δs) ≤ βmin2βmax2 ≤ (1 + δs), which implies that the smaller restricted isometry constant δs leads to the more exact reconstruction, and so it also verifies that the smaller value of βmax/βmin has the better reconstructed performance from the perspective of restricted isometry constant δs. Additionally, singular value decomposition explains that singular value corresponds to the information implied in the matrix and the importance of information is positively correlated with the magnitude of the singular value. Therefore, the maximum non-zero singular value βmax should be improved to maintain the scene details. During the quantitative analysis of βmax/βmin and βmax, the mutual coherence coefficient u is a widely used variable to be introduced [42], this coefficient explains the maximum similarity between any two columns in D N × M and is quantitatively denoted as:
u = max 1 i j M | D i , D j | D i 2 · D j 2 ,
βmax/βmin and βmax are derived to (see Appendix A)
β max β min 1 ( 1 u · M · ( M 1 ) · ( 1 + 1 M N + 1 ) 1 + M N N · ( M 1 ) ) 1 / 2 ,
and
( m + m · M N N · ( M 1 ) ) 1 / 2 β max ( m + m · u · ( M 1 ) ) 1 / 2 ,
respectively. In Equation (8), the smaller the mutual coherence u, the smaller the ratio βmax/βmin is and the exacter the whole reconstruction is. Equation (9) denotes that the smaller the mutual coherence u is, the smaller βmax is and the less the scene details are maintained. Therefore, it is undesirable for that the distribution of random phase to choose u with a large or small value; rather, it should try to maintain the scene details under guaranteeing the whole error. Equations (8) and (9) indicate the choice of pseudo-random space-time modulation. The sparsity achieved by different modulations can refer to Reference [43].

4.3. Carrier of Pseudo-Random Space-Time Modulation

In the pseudo-random space-time modulation, “space” means that different imaging areas have different random phases, and “time” means that different sampling moments have different random phases. This modulation requires the active source of generating the complex exponent with variable random phase on space and time. Then, the question remains as to what is to be taken as the carrier of this modulation.
In the SAR system, antenna is an active converter and is hoped to generate spatially and temporally variable random phase. In reality, the antenna with a random phase radiation pattern already exists. The Lincoln Laboratory in the Massachusetts Institute of Technology (MIT) proposed CRA that is an antenna generating a random phase radiation pattern varying in space and time [44,45]. The beamforming of CRA is based on the following two-dimensional coding: (a) spatial coding performed by introducing dielectric or metallic scatterers on the surface of the reflector, and (b) temporal coding through the use of temporal multiplexing of transmitting and receiving horn arrays [46]. The radiation pattern of CRA is shown in Figure 5b,d. The phased array antenna in the traditional SAR system [47] is used for comparison with CRA, and usually generates phase varying linearly as the beam pointing is not random, as illustrated in Figure 5a,c.
The phased array antenna has a high-gain mainlobe and low-gain sidelobe to guarantee the SAR system performance, e.g., the azimuth ambiguity-to-signal ratio (AASR), range ambiguity-to-signal ratio (RASR), and noise-equivalent sigma zero ( N E σ 0 ) [19]. The CRA generates an amplitude radiation pattern with multiple lobes. In this paper, pseudo-random space-time modulation can be achieved by the combination of the phased array antenna and CRA. When transmitting pulses, the traditional antenna works so that the one-way gain is guaranteed. When receiving echoes, CRA works. Although the two-way gain decreases compared with the traditional SAR system, increasing the transmitting power or the technology of depressing sidelobe, e.g., beamforming technology [48], can be adopted to improve the system performance.

5. Sub-Nyquist SAR Based on Pseudo-Random Space-Time Modulation

In this section, sub-Nyquist SAR based on pseudo-random space-time modulation is introduced from four perspectives: the choice of sub-Nyquist sampling method; the echo signal model; the reconstructed method; and the reconstruction performance.

5.1. Choice of Sub-Nyquist Sampling Method

Compared with the traditional SAR satisfying the Nyquist sampling theorem, SAR based on pseudo-random space-time modulation works with sub-Nyquist samples. In this paper, sub-Nyquist sampling refers to sampling when the sampling frequency is smaller than the Nyquist sampling frequency in this paper. There are two typical methods for sub-Nyquist sampling method: random sub-Nyquist sampling, e.g., random equivalent sampling (RES) [49,50] and multicoset sampling [51,52], and uniform sub-Nyquist sampling [53,54]. The random sub-Nyquist sampling method involves sampling at an arbitrary sampling interval, and the uniform sub-Nyquist sampling method involves sampling at a certain sampling interval.
Because the observational distance is far in the spaceborne SAR system, the time interval between transmitting and receiving is several times the Nyquist pulse repetition time (PRT). The system design should ensure that the echo is completely received and avoids the nadir echo. In the random sub-Nyquist sampling method, the received echo may overlap with the transmitted pulse or the nadir echo. Unless a SAR system adopts one way to avoid the overlapping, the random sub-Nyquist sampling method is not feasible, as illustrated in Figure 6a. Considering that it should avoid the overlapping between transmitting and receiving and that the pseudo-random space-time modulation already guarantees the existence of randomness required by reconstruction, sub-Nyquist samples can be achieved by the uniform sub-Nyquist sampling method in this paper.
Although the uniform sub-Nyquist sampling on the azimuth dimension is sparse, its azimuthal spectrum span is equal to the entire Doppler bandwidth. This lays a good foundation for image reconstruction without the loss of resolution. Additionally, because sub-Nyquist sampling reduces the SNR in Equation (6), the transmitting power can be increased to guarantee information capacity.

5.2. Echo Signal Model after Pseudo-Random Space-Time Modulation

Traditional echo model, it only has the antenna amplitude weighting, not the antenna phase weighting for which the traditional antenna generates a linear phase, and it has little effect on the echo. Unlike with the traditional antenna, the phase is random, so that the phase weighting must be considered in the echo signal model. Therefore, in order to analyze SAR based on the pseudo-random space-time modulation, the corresponding echo signal model should be set up. Different scene positions corresponding to different beam pointing have different time-varying random phase weighting. During the physical formation process of this modulated system based on CRA, it should distinguish the key factors affecting the echo signal, e.g., platform movement, antenna amplitude and phase, and the transmitting signal. After the pseudo-random space-time modulation, the SAR echo for point targets when transmitting the same signal with Equation (1) is denoted as
p 0 ( τ , η ) = x i , y i σ i W i ( τ , η ) r e c t { ( τ 2 R i ( η ) c ) / T r } · exp { j φ i ( η ) }   exp { j π K r [ τ 2 R i ( η ) c ] 2 } · exp { j 4 π R i ( η ) λ } + n ( τ , η )
where φ i ( η ) is the random phase generated by CRA corresponding to the target at (xi,yi). The different scene position (xi,yi) corresponding to different beam pointing has different time-varying random phase φ i ( η ) .

5.3. Reconstructed Method after Pseudo-Random Space-Time Modulation

As mentioned in Section 2, no matter which algorithm is adopted, the reconstructed algorithm includes three steps: the range compression, RCMC, and the azimuth compression. For the range compression, when the transmitted signal is the linear frequency modulated (LFM) signal, the range signal can be compressed based on matched filtering (MF). Most of the already existing algorithms correct range migration in the azimuth frequency domain. When targets on the same range cell have the same range migration in the azimuth frequency domain, correcting the range migration of one target is equivalent to correcting that of targets on the same range cell so that correcting range migration in the frequency domain is efficient [19]. For azimuth compression in the time or frequency domain, MF is also often adopted when it at least has Nyquist samples. However, in the sub-Nyquist SAR based on the pseudo-random space-time modulation, the echo signal is sampled at a sub-Nyquist frequency, so that the azimuth compression based on MF is no longer applicable. Based on the analysis in Section 2, the azimuth signal with sub-Nyquist samples can be recovered by the optimization equation of Equation (5). The demonstration of the reconstructed process is given as follows:
(1) Range Compression
Because the pseudo-random space-time modulated SAR system transmits a LFM signal, range compression can still adopt MF. Assuming that
H ( f τ ) = r e c t ( f τ | K r | T r ) exp ( j π f τ 2 K r ) ,
is matched filter of the range dimension in the frequency domain, the signal after range compression is
p r c ( τ , η ) = IFFT τ { p 0 ( f τ , η ) H ( f τ ) } x i , y i σ i W i ( τ , η ) T r · sin c { K r T r [ τ 2 R i ( η c i ) c ] } · exp { j φ i ( η ) } · exp { j 4 π R i ( η ) λ } + n ( τ , η )
where P 0 ( f τ , η ) is the Fourier transform of the signal p 0 ( τ , η ) on the range dimension, IFFT τ { · } is the inverse Fourier transform on the range dimension, and sin c ( · ) is the sinc function.
(2) Range Cell Migration Correction (RCMC)
In the traditional SAR system, the azimuthal signal is nearly an LFM signal, and the traditional algorithms, e.g., range Doppler algorithm (RDA) or chirp scaling, efficiently correct range migration in the azimuthal frequency domain [19]. After pseudo-random space-time modulation, the azimuthal signal in the frequency domain does not have the explicit expression. The RCMC is implemented in the azimuth time domain in this paper. The procedure of the RCMC is as follows: (1) firstly, it should choose the imaging area at least guaranteeing to contain the observed area in the raw data. Cell division of the imaging area should be smaller than or equal to the resolution of the raw data. When the cell of the imaging area is divided more precisely, the resolution will not change, and this seriously affects the imaging efficiency. Therefore, the cell division of the imaging area should be slightly smaller than or equal to the resolution of the raw data; (2) it takes the ascending sampling to the raw data along the range dimension and calculates the slant range from each division cell of the imaging area to radar at each sampling time. Then, find the position of each division cell on the ascending sampling range cell at each sampling moment according to the time delay calculated by the slant range and take the data on this position; (3) calculate the sum on the data of division cell on the same nearest slant range. After this, the RCMC is finished. After application of the RCMC, the signal is denoted by
p r c m ( τ , η ) = x i , y i σ i W i ( τ , η ) T r · sin c { K r T r [ τ 2 R i ( η c i ) c ] } · exp { j φ i ( η ) } + n ( τ , η )
(3) Azimuth Compression
At the range sampling moment τ 0 , Equation (13) can be expressed corresponding to each sampling moment η n as
p N × 1 = [ p r c m ( τ 0 , η 1 ) p r c m ( τ 0 , η 2 ) p r c m ( τ 0 , η n ) p r c m ( τ 0 , η N ) ] = φ N × M σ M × 1 + n N × 1 ,
where σ M × 1 = [ σ 1 , σ 2 , , σ M ] T , n N × 1 = [ n ( τ 0 , η 1 ) , n ( τ 0 , η 2 ) , , n ( τ 0 , η N ) ] T , and φ N × M = { W i ( τ 0 , η n ) · T r · sin c { K r T r [ τ 0 2 R i ( η c i ) c ] } · exp [ j φ i ( η n ) ] } n = 1 , i = 1 N , M .
Because the samples are achieved with the sub-Nyquist frequency, the optimization equation of Equation (14) without subscript is as follows:
σ ^ = arg min σ { p φ σ 2 2 + γ σ 1 } ,
where σ ^ is the reconstructed result at the range sampling moment τ 0 . The imaging area is reconstructed at all the range sampling moments after the RCMC.

5.4. Reconstructed Performance

The reconstruction equation without subscript is rewritten as:
p = φ σ + n ,
Usually, the components of noise n are approximated as a Gaussian distribution with zero-mean and variance σ n 2 [30]. The scene backscattering cross-section is affected by the observation incidence angle, wavelength, surface structure, etc., and many effect factors are not totally accounted for by a simple deterministic data model [55]. The prior distribution of σ is approximate from SAR image even although speckle noise exists. The statistics of SAR image have been investigated under the assumption of Gaussian statistics with zero-mean and variance σ x 2 for the backscattering cross-section [56,57]. Under the above assumption about the noise n and the scene σ , it follows that:
{ p n ( n ) = p p | σ ( p | σ ) = 1 2 π σ n exp { p φ σ 2 2 2 σ n 2 } p σ ( σ ) = 1 2 π σ x exp { σ 2 2 2 σ x 2 } ,
The knowledge of the reconstruction performance based on CS is interesting. Since the CS recovery algorithm is an estimation method, classical point target evaluation system (e.g., 3 dB resolution, integral side lobe ratio (ISLR), and peak side lobe ratio (PSLR)) are not available in the SAR system based on CS. For a given scene σ , the reconstruction performance can be evaluated by the mean square error (MSE) E { ( σ ^ σ ) 2 } , where σ ^ is the reconstruction result. In the Bayesian framework, the estimation error is simplified as (see Appendix B)
E { ( σ ^ σ ) 2 } t r a c e [ ( φ H φ σ n 2 + I σ x 2 ) 1 ] ,
where t r a c e ( · ) denotes the trace of a matrix and I is the unit matrix.

6. Validation and Analysis

Sub-Nyquist SAR based on pseudo-random space-time modulation is simulated and compared with the traditional SAR with Nyquist samples. The corresponding simulated parameters are listed in Table 1. Figure 7c shows that the urban scene is successfully reconstructed, and its main details are maintained after pseudo-random space-time modulation compared with the result in Figure 7b. Although the random sub-Nyquist sampling method is adopted in Figure 7b and the uniform sub-Nyquist sampling method in Figure 7c, the average PRFs are both 278 Hz.
According to the optimization reconstruction performance in Equation (18), Figure 8 shows that the reconstruction error decreases after pseudo-random space-time modulation.
The simulated results in Figure 9 verify the decision to indicate the choice of pseudo-random space-time modulation in Section 4. As previously mentioned, the phase of complex exponent generated by CRA is a stochastic variable. The simulation chooses that the phase has a uniform distribution with different variances to compare and analyze. After pseudo-random space-time modulation, results in Figure 9a,b are successfully reconstructed with sub-Nyquist samples compared to Figure 1d. Results of Figure 9a,b are simulated when the phase of pseudo-random space-time modulation obeys a uniform distribution with different variances. Results of Figure 9c,d denote the curves of mutual coherence coefficient, and singular value amplitude under two distributions. For the convenience of description, assuming the mutual coherence coefficient, the maximum singular value and the minimum singular value under two different random distributions are denoted by uu, ug, βmax_u, βmax_g, βmin_u, and βmin_g, respectively. As illustrated in Figure 9, uu is larger than ug in Figure 9c, βmax_u and βmax_u/βmin_u are larger than βmax_g and βmax_g/βmin_g in Figure 9d, respectively, and more road details are maintained in the red ellipsoid of Figure 9a compared with Figure 9b. It is consistent with Equations (8) and (9) that random phase distribution with large mutual coherence leads to many scene details under guaranteeing the whole error. Pseudo-random space-time can improve the imaging performance, and Equations (8) and (9) indicate the choice of the modulation.

7. Conclusions

This paper presents sub-Nyquist SAR based on pseudo-random space-time modulation. To achieve this system, CRA with high-sensing capacity was adopted to be a carrier of pseudo-random space-time modulation. The following two major findings have been obtained:
(1) To explain the reason for the constraint on the scene sparsity in the traditional sub-Nyquist SAR, this paper establishes an information channel model of microwave imaging radar based on information theory. Under this model, noisy-channel coding theorem indicates the necessary condition of the successful reconstruction and explains the reason for this constraint.
(2) To relieve the constraint on the scene sparsity so that the scene with large sparsity can be reconstructed, this paper proposes sub-Nyquist SAR based on pseudo-random space-time modulation. This modulation is the phase modulation of the traditional raw data and it can reduce the mutual information of successive measurements so that more scene information can be obtained.
The proposed method is applied not only to the strip mapping mode but also to other modes such as sliding spotlight and terrain observation with progressive scan (TOPS) modes. Our future research will focus on a wide-swath SAR system based on pseudo-random space-time modulation.

Author Contributions

Data curation, W.C. and Z.Y.; Formal analysis, W.C. and Z.Y.; Methodology, W.C.; Z.Y.; C.L. and P.X.; Project administration, Z.Y.; Software, W.C.; Validation, Z.Y.; Visualization, W.C.; Writing—original draft preparation, W.C., Z.Y., C.L. and P.X.; Writing—review and editing, W.C., Z.Y., C.L. and P.X.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Assuming that A M × M = D N × M H D N × M and m = D i 2 2 ¯ , i { 1 , 2 , , M } . For that A M × M is Hermitian matrix and N < M , A M × M has N non-zero real eigenvalues α 1 , α 2 , , α N in descending numerical order, M N zero eigenvalues α N + 1 = α N + 2 = = α M = 0 . The singular value of D N × M is β i = α i . The trace of A M × M and A M × M 2 are
t r a c e ( A M × M ) = n = 1 N α n = M · m ,
t r a c e ( A M × M 2 ) = n = 1 N α n 2 = i = 1 M j = 1 M | D i , D j | 2 = M · m 2 + 2 · i = 1 M j = i + 1 M | D i , D j | 2 M · m 2 + M · ( M 1 ) · u 2 · m 2
Cauchy-Schwarz inequality shows
n = 1 N α n 2 ( n = 1 N α n ) 2 N = M 2 · m 2 N ,
Combining Equations (A2) and (A3) gives:
M 2 · m 2 N t r a c e ( A M × M 2 ) M · m 2 + M · ( M 1 ) · u 2 · m 2 ,
Assuming that b = t r a c e ( A M × M 2 ) / M m 2 , Equation (A4) is substituted into b :
m · M N N b m · u · M 1 ,
Reference [58] indicates
m + b / M 1 α max m + b · M 1 ,
α k α l b · M · ( 1 k + 1 M l + 1 ) , 1 k l M ,
Equation (A5) is substituted into Equation (A6),
( m + m · M N N · ( M 1 ) ) 1 / 2 β max ( m + m · u · ( M 1 ) ) 1 / 2 ,
When k = 1 and l = N , Equation (A7) is simplified as
α 1 α N b · M · ( 1 + 1 M N + 1 ) ,
It divides by α 1 at two sides of Equation (A9) and applies the range of α 1 ,
1 α N α 1 b · M · ( 1 + 1 M N + 1 ) α 1 b · M · ( 1 + 1 M N + 1 ) m + m · M N N · ( M 1 ) ,
The ratio of the maximum to minimum eigenvalue α 1 / α N is achieved by Equation (A10)
α 1 α N 1 1 b · M · ( 1 + 1 M N + 1 ) m + m · M N N · ( M 1 ) ,
β 1 / β N is further achieved by applying the range of b
β 1 β N 1 ( 1 u · M · ( M 1 ) · ( 1 + 1 M N + 1 ) 1 + M N N · ( M 1 ) ) 1 / 2 ,

Appendix B

Assuming that the vector σ is a random variable vector, the mean square error (MSE) of the estimation σ ^ in a Bayesian framework is denoted by Reference [59]
MSE ( σ ^ ) = m = 1 M E { ( σ ^ m σ m ) 2 } t r a c e ( { E [ 2 ln p p , σ ( p , σ ) σ m 1 σ m 2 ] } 1 ) = t r a c e ( { E [ 2 ln p p / σ ( p / σ ) σ m 1 σ m 2 ] E [ 2 ln p σ ( σ ) σ m 1 σ m 2 ] } 1 )
where σ ^ m and σ m are the mth element of σ ^ and σ , respectively. t r a c e ( · ) denotes the trace of a matrix. p p , σ ( p , σ ) , p p / σ ( p / σ ) and p σ ( σ ) are the joint probability density function of the vector p and σ , the conditional probability density function of the vector p and σ , the prior probability density function of the vector σ , respectively. The data information matrix J D and the prior information matrix J P represent E [ 2 ln p p / σ ( p / σ ) σ m 1 σ m 2 ] and E [ 2 ln p σ ( σ ) σ m 1 σ m 2 ] , respectively.
The following analyzes the Bayesian information matrix J B = J D + J P from these two information matrices:
(1)
the calculation of the data information matrices J D :
J D = E [ 2 ln p p / σ ( p / σ ) σ m 1 σ m 2 ] = E [ 1 2 σ n 2 · p φ σ 2 2 σ m 1 σ m 2 ] = E { 1 2 σ n 2 · σ m 1 σ m 2 [ p H p p H φ σ σ H φ H σ + σ H φ H φ σ ] } = E { 1 2 σ n 2 · σ m 2 [ 2 φ H φ σ 2 φ H p ] } = φ H φ σ n 2
(2)
the calculation of the prior information matrix J P :
J P = E [ 2 ln p σ ( σ ) σ m 1 σ m 2 ] = E [ 1 2 σ x 2 · σ 2 2 σ m 1 σ m 2 ] = I σ x 2
where I is the unit matrix. Substituting Equations (A13) and (A14) into Equation (A15), we obtain
E { ( σ ^ σ ) 2 } t r a c e [ ( φ H φ σ n 2 + I σ x 2 ) 1 ] ,

References

  1. Abo-Zahhad, M.M.; Hussein, A.I.; Mohamed, A.M. Compressive sensing algorithms for signal processing applications: A survey. Int. J. Commun. Netw. Syst. Sci. 2015, 8, 197–216. [Google Scholar]
  2. Matusiak, E.; Eldar, Y.C. Sub-Nyquist sampling of short pulses. IEEE Trans. Signal Process. 2012, 60, 1134–1148. [Google Scholar] [CrossRef]
  3. Xiao, L.; Xia, X. Frequency determination from truly sub-Nyquist samplers based on robust Chinese remainder theorem. Signal Process. 2018, 15, 248–258. [Google Scholar] [CrossRef]
  4. Grigoryan, R. Acquisition of Multi-Band Signals via Compressed Sensing. Ph.D. Thesis, Aalborg University, Aalborg, Denmark, 2015. [Google Scholar]
  5. Lagunas, E. Compressive Sensing Based Candidate Detector and Its Applications to Spectrum Sensing and through-the-Wall Radar Imaging. Ph.D. Thesis, Universitat Politecnica de Catalunya (UPC), Barcelona, Spain, 2014. [Google Scholar]
  6. Liu, C.; Wang, H.; Zhang, J.; He, Z. Wideband Spectrum Sensing Based on Single-Channel Sub-Nyquist Sampling for Cognitive Radio. Sensors 2018, 18, 2222. [Google Scholar] [CrossRef]
  7. Aberman, K.; Eldar, Y.C. Sub-Nyquist SAR via Fourier domain range-Doppler processing. IEEE Trans. Geosci. Remote Sens. 2017, 55, 6228–6244. [Google Scholar] [CrossRef]
  8. Nyquist, H. Certain topics in telegraph transmission theory. Trans. Am. Inst. Electr. Eng. 1928, 47, 617–644. [Google Scholar] [CrossRef]
  9. Shannon, C.E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef]
  10. Martian, A.; Craciunescu, R.; Vulpe, A.; Suciu, G.; Fratu, O. Access to RF white spaces in Romania: Present and future. Wirel. Pers. Commun. 2016, 87, 693–712. [Google Scholar] [CrossRef]
  11. Yang, J.G.; Thompson, J.; Huang, X.T.; Jin, T.; Zhou, Z.M. Segmented reconstruction for compressed sensing SAR imaging. IEEE Trans. Geosci. Remote Sens. 2013, 51, 4214–4225. [Google Scholar] [CrossRef]
  12. Potter, L.C.; Ertin, E.; Parker, J.T.; Cetin, M. Sparsity and compressed sensing in radar imaging. Proc. IEEE 2010, 98, 1006–1020. [Google Scholar] [CrossRef]
  13. Ender, J.H.G. On compressive sensing applied to radar. Signal Process. 2010, 90, 1402–1414. [Google Scholar] [CrossRef]
  14. Herman, M.A.; Strohmer, T. High-Resolution Radar via Compressed Sensing. IEEE Trans. Signal Process. 2009, 57, 2275–2284. [Google Scholar] [CrossRef] [Green Version]
  15. Alonso, M.T.; López-Dekker, P.; Mallorquí, J.J. A novel strategy for radar imaging based on compressive sensing. IEEE Trans. Geosci. Remote Sens. 2011, 48, 4285–4295. [Google Scholar] [CrossRef]
  16. Ni, J.; Zhang, Q.; Luo, Y.; Sun, L. Compressed sensing SAR imaging based on centralized sparse representation. IEEE Sens. J. 2018, 18, 4920–4932. [Google Scholar] [CrossRef]
  17. Bi, H.; Bi, G.; Zhang, B.; Hong, W. Complex-image-based sparse SAR imaging and its equivalence. IEEE Trans. Geosci. Remote Sens. 2018, 56, 5006–5013. [Google Scholar] [CrossRef]
  18. Curlander, J.C.; McDonough, R.N. Synthetic Aperture Radar: Systems and Signal. Processing; Wiley: New York, NY, USA, 1991. [Google Scholar]
  19. Cumming, I.; Wong, F. Digital Signal Processing of Synthetic Aperture Radar Data: Algorithms and Implementation; Artech House: Norwood, MA, USA, 2004. [Google Scholar]
  20. Foucart, S.; Rauhut, H. A Mathematical Introduction to Compressive Sensing; Birkhäuser: Basel, Switzerland, 2013. [Google Scholar]
  21. Weed, J. Approximately certifying the restricted isometry property is hard. IEEE Trans. Inf. Theory 2017, 64, 5488–5497. [Google Scholar] [CrossRef]
  22. Romberg, J.K. Compressive sensing by random convolution. SIAM J. Imaging Sci. 2009, 2, 1098–1128. [Google Scholar] [CrossRef]
  23. Yang, D.; Liao, G.; Zhu, S.; Yang, X.; Zhang, X. SAR imaging with undersampled data via matrix completion. IEEE Geosci. Remote Sens. Lett. 2014, 11, 1539–1543. [Google Scholar] [CrossRef]
  24. Tropp, J.A.; Gilbert, A.C. Signal recovery from partial information via orthogonal matching pursuit. IEEE Trans. Inf. Theory 2007, 53, 4655–4666. [Google Scholar] [CrossRef]
  25. Cerrone, C.; Cerull, R.; Golden, B. Carousel greedy: A generalized greedy algorithm with applications in optimization. Comput. Oper. Res. 2017, 85, 97–112. [Google Scholar] [CrossRef]
  26. Candes, E.; Tao, T. The Dantzig selector: Statistical estimation when p is much larger than n. Ann. Stat. 2007, 35, 2313–2351. [Google Scholar] [CrossRef]
  27. Tibshirani, R. Regression shrinkage and selection via the Lasso. J. R. Stat. Soc. B 1996, 58, 267–288. [Google Scholar] [CrossRef]
  28. Chen, S.S.; Donoho, D.L.; Saunders, M.A. Atomic decomposition by basis pursuit. SIAM J. Sci. Comput. 1998, 20, 33–61. [Google Scholar] [CrossRef]
  29. Chen, X.; Wang, H.; Wang, R. A null space analysis of the L1 synthesis method in frame-based compressed sensing. Appl. Comput. Harmon. Anal. 2014, 37, 492–515. [Google Scholar] [CrossRef]
  30. Ji, S.; Xue, Y.; Carin, L. Bayesian compressive sensing. IEEE Trans. Signal Process. 2008, 56, 2346–2356. [Google Scholar] [CrossRef]
  31. Babacan, S.D.; Molina, R.; Katsaggelos, A.K. Bayesian compressive sensing using Laplace priors. IEEE Trans. Image Process. 2010, 19, 53–63. [Google Scholar] [CrossRef] [PubMed]
  32. Baron, D.; Sarvotham, S.; Baraniuk, R.G. Bayesian compressive sensing via belief propagation. IEEE Trans. Signal Process. 2010, 58, 269–280. [Google Scholar] [CrossRef]
  33. Arjoune, Y.; Kaabouch, N.; Ghazi, H.E.; Tamtaoui, A. Compressive sensing: Performance comparison of sparse recovery algorithms. In Proceedings of the Computing and Communication Workshop and Conference (CCWC), Las Vegas, NV, USA, 9–11 January 2017; pp. 1–7. [Google Scholar]
  34. Zhang, J.; Yang, K.; Liu, F.; Zhang, Y. Information-theoretic characterization and undersampling ratio determination for compressive radar imaging in a simulated environment. Entropy 2015, 17, 5171–5198. [Google Scholar] [CrossRef]
  35. Zhang, J.; Yang, K. Informational analysis for compressive sampling in radar imaging. Sensors 2015, 15, 7136–7155. [Google Scholar] [CrossRef]
  36. Verdú, S. Fifty Years of Shannon Theory. IEEE Trans. Inf. Theory 1998, 44, 2057–2078. [Google Scholar] [CrossRef]
  37. Golub, G.; Loan, C.F.V. Matrix Computations, 3rd ed.; Johns Hopkins University Press: Baltimore, MD, USA, 1996. [Google Scholar]
  38. Raney, R.K.; Runge, H.; Bamler, R.; Cumming, I.G.; Wong, F.H. Precision SAR processing using chirp scaling. IEEE Trans. Geosci. Remote Sens. 1994, 32, 786–799. [Google Scholar] [CrossRef]
  39. Shanmugam, K.S. Digital and Analog Communication Systems; Wiley: New York, NY, USA, 1979. [Google Scholar]
  40. Tse, D.; Viswanath, P. Fundamentals of Wireless Communication; Cambridge Univ. Press: Cambridge, UK, 2005. [Google Scholar]
  41. Donoho, D.L. Compressed sensing. IEEE Trans. Inf. Theory 2006, 52, 1289–1306. [Google Scholar] [CrossRef]
  42. Donoho, D.L.; Huo, X. Uncertainty principles and ideal atomic decomposition. IEEE Trans. Inf. Theory 2001, 47, 2845–2862. [Google Scholar] [CrossRef] [Green Version]
  43. Arjoune, Y.; Kaabouch, N.; el Ghazi, H.; Tamtaoui, A. A performance comparison of measurement matrices in compressive sensing. WILEY Int. J. Commun. Syst. 2018, 31, 1–18. [Google Scholar] [CrossRef]
  44. Lorenzo, J.A.M.; Juesas, J.H.; Blackwell, W. Single-transceiver compressive antenna for high-capacity sensing and imaging applications. In Proceedings of the 2015 9th European Conference on Antennas and Propagation (EuCAP), Lisbon, Portugal, 13–17 April 2015; pp. 1–3. [Google Scholar]
  45. Martinez-Lorenzo, J.A.; Juesas, J.H.; Blackwell, W. A single-transceiver compressive reflector antenna for high-sensing-capacity imaging. IEEE Antennas Wirel. Propag. Lett. 2015, 15, 968–971. [Google Scholar] [CrossRef]
  46. Molaei, A.; Juesas, J.H.; Lorenzo, J.A.M. Compressive Reflector Antenna Phased Array. In Antenna Arrays and Beam-Formation; InTechOpen: Lodon, UK, 2017. [Google Scholar] [Green Version]
  47. Mailloux, R.J. Phased Array Antenna Handbook; Artech House: Boston, MA, USA, 1994. [Google Scholar]
  48. Sohrabi, F.; Yu, W. Hybrid digital and analog beamforming design for large-scale antenna arrays. IEEE J. Sel. Top. Signal Process. 2016, 10, 501–513. [Google Scholar] [CrossRef]
  49. Sonnaillon, M.O.; Urteaga, R.; Bonetto, F.J. High-Frequency Digital Lock-In Amplifier Using Random Sampling. IEEE Trans. Instrum. Meas. 2008, 57, 616–621. [Google Scholar] [CrossRef]
  50. Zhao, Y.; Wang, L.; Wang, H.; Liu, C. Minimum rate sampling and spectrum blind reconstruction in random equivalent sampling. Circuits Syst. Signal Process. 2015, 34, 2667–2680. [Google Scholar] [CrossRef]
  51. Herley, C.; Wong, P.W. Minimum rate sampling and reconstruction of signals with arbitrary frequency support. IEEE Trans. Inf. Theory 1999, 45, 1555–1564. [Google Scholar] [CrossRef]
  52. Ma, Y.; Gao, Y.; Liang, Y.C.; Cui, S. Reliable and efficient sub-Nyquist wideband spectrum sensing in cooperative cognitive radio networks. IEEE J. Sel. Areas Commun. 2016, 34, 2750–2762. [Google Scholar] [CrossRef]
  53. Goodman, N.; Lin, S.C.; Rajakrishna, D.; Stiles, J. Processing of multiple-receiver spaceborne arrays for wide-area SAR. IEEE Trans. Geosci. Remote Sens. 2002, 40, 841–852. [Google Scholar] [CrossRef] [Green Version]
  54. Di Martino, G.; Iodice, A. Coprime Synthetic Aperture Radar (CopSAR): A new acquisition mode for maritime surveillance. IEEE Rans. Geosci. Remote Sens. 2015, 53, 3110–3123. [Google Scholar] [CrossRef]
  55. Kuruoglu, E.E.; Zerubia, J. Modeling SAR images with a generalization of the Rayleigh distribution. IEEE Trans. Image Process. 2004, 13, 527–533. [Google Scholar] [CrossRef] [PubMed]
  56. Sarabandi, K. Derivation of phase statistics from the Mueller matrix. Radio Sci. 1992, 27, 553–560. [Google Scholar] [CrossRef]
  57. Kong, J.A.; Swartz, A.A.; Yueh, H.A.; Novak, L.M.; Shin, R.T. Identification of terrain cover using the optimum polarimetric classifier. J. Electromagnet. Waves Appl. 1987, 2, 171–194. [Google Scholar]
  58. Seber, G. A Matrix Handbook for Statisticians; Wiley-Interscience: Hoboken, NJ, USA, 2008. [Google Scholar]
  59. Hvan Trees, L.; Bell, K.L.; Tian, Z. Detection, Estimation, and Modulation Theory: Part 1: Detection, Estimation, and Filtering Theory, 2nd ed.; Wiley: Hoboken, NJ, USA, 2013. [Google Scholar]
Figure 1. (a) Reconstructed result of ocean containing ships from Nyquist samples based on chirp scaling; (b) Reconstructed result of ocean containing ships from sub-Nyquist samples based on compressive sensing (CS); (c) Reconstructed result of an urban scene from Nyquist samples based on chirp scaling; (d) Reconstructed result of an urban scene from sub-Nyquist samples based on CS.
Figure 1. (a) Reconstructed result of ocean containing ships from Nyquist samples based on chirp scaling; (b) Reconstructed result of ocean containing ships from sub-Nyquist samples based on compressive sensing (CS); (c) Reconstructed result of an urban scene from Nyquist samples based on chirp scaling; (d) Reconstructed result of an urban scene from sub-Nyquist samples based on CS.
Sensors 18 04343 g001
Figure 2. Analogous to the communication system in information theory, the scene, the mapping relation, and the echo correspond to information source, channel, and sink, respectively.
Figure 2. Analogous to the communication system in information theory, the scene, the mapping relation, and the echo correspond to information source, channel, and sink, respectively.
Sensors 18 04343 g002
Figure 3. Channel capacity against signal-to-noise ratio (SNR) at different sampling frequencies. As the sampling frequency decreases, channel capacity decreases.
Figure 3. Channel capacity against signal-to-noise ratio (SNR) at different sampling frequencies. As the sampling frequency decreases, channel capacity decreases.
Sensors 18 04343 g003
Figure 4. The information channel after pseudo-random space-time modulation. Information channel is based on Doppler movement under the sub-Nyquist sampling method. After pseudo-random space-time modulation, the mutual information between information source X1 and information sink Y is increased so that the scene with large value of sparsity is reconstructed.
Figure 4. The information channel after pseudo-random space-time modulation. Information channel is based on Doppler movement under the sub-Nyquist sampling method. After pseudo-random space-time modulation, the mutual information between information source X1 and information sink Y is increased so that the scene with large value of sparsity is reconstructed.
Sensors 18 04343 g004
Figure 5. The amplitude and phase patterns of the phased array antenna and compressive reflector antenna (CRA), compared with the linear phase of the phased array antenna; the phase of the CRA is random.
Figure 5. The amplitude and phase patterns of the phased array antenna and compressive reflector antenna (CRA), compared with the linear phase of the phased array antenna; the phase of the CRA is random.
Sensors 18 04343 g005
Figure 6. Illustration of sub-Nyquist sampling method for sampling duration T1. The dashed line denotes the Nyquist sample. The solid line denotes the sub-Nyquist sample in random and uniform sub-Nyquist sampling method randomly and uniformly selects from dashed line, respectively. (a) random sub-Nyquist sampling: the conflict between transmitted pulse and receiving echo. (b) uniform sub-Nyquist sampling: echo is completely received and has no conflict with the transmitted pulse.
Figure 6. Illustration of sub-Nyquist sampling method for sampling duration T1. The dashed line denotes the Nyquist sample. The solid line denotes the sub-Nyquist sample in random and uniform sub-Nyquist sampling method randomly and uniformly selects from dashed line, respectively. (a) random sub-Nyquist sampling: the conflict between transmitted pulse and receiving echo. (b) uniform sub-Nyquist sampling: echo is completely received and has no conflict with the transmitted pulse.
Sensors 18 04343 g006
Figure 7. (a) Reconstructed result from Nyquist samples; (b) Reconstructed result from sub-Nyquist samples; (c) Reconstructed result from sub-Nyquist samples after the pseudo-random space-time modulation.
Figure 7. (a) Reconstructed result from Nyquist samples; (b) Reconstructed result from sub-Nyquist samples; (c) Reconstructed result from sub-Nyquist samples after the pseudo-random space-time modulation.
Sensors 18 04343 g007
Figure 8. The reconstruction error before and after the pseudo-random time-space modulation.
Figure 8. The reconstruction error before and after the pseudo-random time-space modulation.
Sensors 18 04343 g008
Figure 9. Simulated results of (a,b) reconstructed scene, (c) mutual coherence coefficient, and (d) singular value.
Figure 9. Simulated results of (a,b) reconstructed scene, (c) mutual coherence coefficient, and (d) singular value.
Sensors 18 04343 g009aSensors 18 04343 g009b
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParameterValues of the Nyquist Sampling MethodValues of the Sub-Nyquist Sampling Method
Orbital height (km)693693
Wavelength (m)0.05550.0555
Pulse width (µs)5050
Antenna height (m)8.938.93
Signal bandwidth (MHz)100100
The sampling frequency (MHz)110110
Incidence angles (°)38.73–43.5038.73–43.50
(Average) pulse repetition frequency (PRF) (Hz)1946278

Share and Cite

MDPI and ACS Style

Chen, W.; Li, C.; Yu, Z.; Xiao, P. Sub-Nyquist SAR Based on Pseudo-Random Time-Space Modulation. Sensors 2018, 18, 4343. https://doi.org/10.3390/s18124343

AMA Style

Chen W, Li C, Yu Z, Xiao P. Sub-Nyquist SAR Based on Pseudo-Random Time-Space Modulation. Sensors. 2018; 18(12):4343. https://doi.org/10.3390/s18124343

Chicago/Turabian Style

Chen, Wenjiao, Chunsheng Li, Ze Yu, and Peng Xiao. 2018. "Sub-Nyquist SAR Based on Pseudo-Random Time-Space Modulation" Sensors 18, no. 12: 4343. https://doi.org/10.3390/s18124343

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