Next Article in Journal
Federation of Internet of Things Testbeds for the Realization of a Semantically-Enabled Multi-Domain Data Marketplace
Next Article in Special Issue
Synthetic Aperture Radar Processing Approach for Simultaneous Target Detection and Image Formation
Previous Article in Journal
Gradient Projection with Approximate L0 Norm Minimization for Sparse Reconstruction in Compressed Sensing
Previous Article in Special Issue
Wide-Angle, Ultra-Wideband ISAR Imaging of Vehicles and Drones
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel MIMO–SAR Solution Based on Azimuth Phase Coding Waveforms and Digital Beamforming

1
School of Computer and Information, Hefei University of Technology, Hefei 230009, China
2
East China Research Institute of Electronic Engineering, Hefei 230088, China
3
Institute of Electronic Engineering, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(10), 3374; https://doi.org/10.3390/s18103374
Submission received: 24 July 2018 / Revised: 5 October 2018 / Accepted: 6 October 2018 / Published: 9 October 2018
(This article belongs to the Special Issue Automatic Target Recognition of High Resolution SAR/ISAR Images)

Abstract

:
In multiple-input multiple-output synthetic aperture radar (MIMO–SAR) signal processing, a reliable separation of multiple transmitted waveforms is one of the most important and challenging issues, for the unseparated signal will degrade the performance of most MIMO–SAR applications. As a solution to this problem, a novel APC–MIMO–SAR system is proposed based on the azimuth phase coding (APC) technique to transmit multiple waveforms simultaneously. Although the echo aliasing occurs in the time domain and Doppler domain, the echoes can be separated well without performance degradation by implementing the azimuth digital beamforming (DBF) technique, comparing to the performance of the orthogonal waveforms. The proposed MIMO–SAR solution based on the APC waveforms indicates the feasibility and the spatial diversity of the MIMO–SAR system. It forms a longer baseline in elevation, which gives the potential to expand the application of MIMO–SAR in elevation, such as improving the performance of multibaseline InSAR and three-dimensional SAR imaging. Simulated results on both a point target and distributed targets validate the effectiveness of the echo separation and reconstruction method with the azimuth DBF. The feasibility and advantage of the proposed MIMO–SAR solution based on the APC waveforms are demonstrated by comparing with the imaging result of the up- and down-chirp waveforms.

1. Introduction

Synthetic aperture radar (SAR) is a powerful remote sensing technique independent of weather and sunlight illumination. With multiple transmitters and multiple receivers employed, multiple-input multiple-output SAR (MIMO–SAR) enables not only the acquisition of additional phase centers and long baselines for high-resolution wide-swath (HRWS) SAR imaging [1,2,3], but also SAR applications like multibaseline interferometry or three-dimensional imaging [4,5,6,7,8]. Moreover, it enables the possibility to utilize multiple SAR observing modes simultaneously in one MIMO–SAR system.
Design of the transmitted waveforms from multiple transmitters is one of the most important and challenging issues in MIMO–SAR implementation. To this end, time-division multiplexing (TDM) waveforms [9,10,11], frequency-division multiplexing (FDM) waveforms [12,13] and code-division multiplexing (CDM) waveforms [14,15,16,17,18] have been proposed in recent decades. TDM uses a time filter, FDM uses a bandpass filter, while the CDM usually uses a matched filter for the reliable separation of radar echoes. In those waveforms, the CDM waveforms, especially the orthogonal waveforms [14,15,16,17,18], are widely discussed for their capacity to obtain a high-resolution wide-swath SAR image. However, those transmitted orthogonal signals share the same frequency band. Thus, the signal energies of all unmatched waveforms are present simultaneously in the focused signal, which will degrade the performance of most MIMO–SAR systems [19,20,21]. For this reason, these mutually orthogonal waveforms do not suit the senior implementation of the MIMO–SAR system very well.
As a solution to this challenge, the azimuth phase coding (APC) waveforms have been presented in the literature, References [22,23], where the Doppler bandpass filter can enable reliable separation of the echoes corresponding to each transmitted waveform. However, the main drawback of the APC technique in [22,23] is that the pulse repetition frequency (PRF) needs to be K-times as large as the Doppler bandwidth, where K is the number of the simultaneously transmitting waveforms. That drawback will lead to a significant reduction of the swath coverage and will limit the application in spaceborne SAR. Thus, finding a way to reduce the PRF is necessary and beneficial.
By transmitting a set of short-term shift-orthogonal waveforms, the radar echoes can be separated by the digital beamforming (DBF) technique [24], which makes the reduction of the PRF possible. Krieger et al. introduced a set of short-term shift-orthogonal waveforms and the DBF technique in elevation in [25], which put up an effective and inspiring approach to the echo separation for MIMO–SAR.
In this paper, a novel approach to apply the APC waveforms in combination with the DBF technique in the azimuth is proposed as a new MIMO–SAR solution, named APC–MIMO–SAR. The simultaneously transmitted APC waveforms are modulated to carry a set of phases, which makes the Doppler of the received signal distinguishable. By utilizing an antenna with several additional azimuth-displaced phase centers [26], the echoes can be well separated by the DBF technique. Additionally, the spatial diversity of the MIMO–SAR can be restored. In this case, the PRF of this MIMO–SAR system just needs to be slightly larger than the Doppler bandwidth [27,28]. The proposed APC waveforms, together with the azimuth DBF, can be used to exploit the potentials of MIMO–SAR without performance degradation.
This paper is organized as follows. In Section 2, the proposed MIMO–SAR solution with the APC waveforms is introduced in detail. Section 3 proposes an echo separation method by the azimuth DBF processing. Section 4 discusses the signal-to-noise ratio and the requirement of the antenna array size. Section 5 presents the simulations on both the point target and distributed targets to show the advantage of the APC–MIMO–SAR and the effectiveness of the proposed solution. Section 6 draws the conclusions.

2. Methods

This section describes the aim, architecture and advantage of the proposed APC–MIMO–SAR solution. The architecture of the MIMO–SAR system is introduced first, and then the modulated APC waveform adopted in this paper is analyzed. Finally, we discuss new problems in the signal model and Doppler spectra.

2.1. Architecture for Proposed MIMO–SAR System

A novel instrument architecture for transmitting and receiving in the MIMO–SAR has been discussed in detail in this section. This architecture can provide a longer baseline in elevation, which is of benefit to expand the application of MIMO–SAR.
Assuming that the whole aperture is divided into N × M sub-apertures, N in elevation and M in the azimuth. K sub-apertures in the first column are chosen to be the transmitters, and all the sub-apertures are chosen to be the receivers. Thus, there are K transmitters and N × M receivers, apparently K ≤ N. The N × M receiving sub-apertures can receive the echoes simultaneously. The k-th transmitting sub-aperture, in the k-th row and the first column, is represented as Tx k ( k   =   1 ,   2 ,     , K ), while the nm-th receiving sub-aperture, in the n-th row and the m-th column, is represented as Rx n m ( n   =   1 ,   2 ,   ,   N , m   = 1 ,   2 ,   ,   M ). For separating the echoes in the following processing, plenty of spatial degrees of freedom are needed, which means the azimuth receivers’ number M should be not smaller than the transmitters’ number N (N ≤ M).
Since the architecture of the proposed MIMO–SAR system is diversiform, one possible architecture of the transmitting and receiving aperture is shown in Figure 1. For convenience, Figure 1 gives a 3 × 3 sub-aperture.
The proposed MIMO–SAR system architecture can provide additional displaced phase centers [26] both in elevation and azimuth, comparing to a SIMO (single-input multi-output) SAR system. Figure 2 shows the sketch of the elevation- and azimuth-displaced phase centers. In Figure 2, the antenna aperture and the displaced phase centers of two kinds of SIMO–SAR systems are listed on the left, comparing to those of the proposed MIMO–SAR system on the right. Though the displaced phase centers of the proposed MIMO–SAR system are in the same plane, we disassemble them into three layers just for a clearer expression. The first layer marked with “•” corresponds to transmitter Tx 1 , the second one marked with “×” corresponds to transmitter Tx 2 , and the third one marked with “△” corresponds to transmitter Tx 3 . The dotted lines connecting different layers indicate that those phase centers are in the same position on the original plane.
From Figure 3, it is easy to conclude that the long baseline in elevation is formed, doubling the length of the baseline in elevation in the SIMO–SAR system. The longer baseline in elevation provides the additional elevation-displaced phase centers, which can be used to obtain additional degrees of distance freedom. This benefits the MIMO–SAR processing in two aspects. One is to separate echoes and suppress range ambiguities [29,30,31], and the other is to extend the application into multibaseline InSAR [32] and three-dimensional SAR imaging [33]. The additional azimuth-displaced phase centers provide more spatial degrees of freedom, which will make the separation of the echoes’ overlapped spectra via the DBF processing in the azimuth possible [34,35]. This is the new aspect in our paper, which can be used to exploit the potentials of MIMO–SAR without performance degradation.

2.2. Modulated APC Waveforms

The APC technique is conceived for single transmit antenna SAR systems to suppress range ambiguities [29,30,31], and it is applied to MIMO–SAR for the design of APC waveforms proposed by Cristallini et al. [22] and Meng et al. [23]. This generation of APC waveforms is proposed to modulate an individual azimuth phase on the original transmit signal for each transmitting aperture, and to ensure the Doppler spectra of echoes from different transmitting apertures occupy different Doppler bands without overlapping each other. Both the conventional chirp signal and the orthogonal waveforms can be used as the original signal/signals to generate a class of the APC waveforms. In the receivers, the Doppler bandpass filter can be used to separate the echoes reliably from different transmit antennas in the Doppler domain. However, the PRF needs to be at least K-times as large as the Doppler bandwidth to get an unaliasing Doppler spectrum, when transmitting K kinds of modulated APC waveforms simultaneously. The overclaim on PRF in [22] and [23] significantly narrows the swath coverage and limits the waveforms’ application, especially in spaceborne MIMO–SAR. To solve this problem, a novel approach to apply the APC waveforms mentioned above is proposed in this section, named APC–MIMO–SAR. The APC–MIMO–SAR receives the echoes by the azimuth sub-apertures to get additional azimuth-displaced phase centers, which makes the reduction of the PRF possible.
Figure 3 shows the scheme of the generation of the modulated APC waveforms. An APC–MIMO–SAR system with K transmitting sub-apertures can transmit K kinds of modulated APC waveforms simultaneously. The APC modulation phase for the k-th ( k   = 1 ,   2 ,   ,   K ) transmit waveform ( Tx k ) is given by
φ mod , k ( l ) = exp ( j π K ( l + k 1 ) 2 )  
where l denotes the sequence number of the azimuth pulse. Apparently, l can be expressed as l = t a f P R F , where f P R F denotes the pulse repetition frequency (PRF) and t a denotes the azimuth slow time. Thus, substituting l = t a f P R F into (1), the APC modulation phase can be rewritten as
φ mod , k ( t a ) = exp ( j π K ( t a f P R F + k 1 ) 2 )  
There are N × M receiving sub-apertures which can receive echoes simultaneously, as Figure 4 shows. In order to separate the echoes, the receiver needs to do a demodulating processing after receiving. The APC demodulation phases for each of the receiving sub-apertures are defined as
φ d e ( t a )   =   exp ( j π K t a 2 f P R F 2 )  
Obviously, the APC demodulation phase is a function which does not vary with the sequence number Rxnm of the receiving sub-aperture.
From Figure 3 and Figure 4, one can notice the convenience of the transmitting and receiving of the APC waveforms. The APC modulation and demodulation can be done by adding a multiplier to the transmitter and receiver, instead of redesigning the signal generation module.
After the APC modulation and demodulation, the residual phase of the k-th APC waveform can be expressed as
φ r e s , k ( t a ) = φ mod , k ( t a ) φ d e ( t a ) = exp ( j 2 π ( k 1 ) f P R F K t a ) exp ( j π ( k 1 ) 2 K )
From (4), one can notice that the residual phase is in connection with the APC waveform order k. The first exponential term is a linear phase of the azimuth slow time ta, and the second exponential term is a constant phase. The first exponential term will be expressed as an additional Doppler shift in the azimuth Doppler domain, and the additional Doppler shift frequency Δ f d , k can be written as
Δ f d , k = ( k 1 ) f PRF K  
From (5), it is explicit that there is no additional Doppler shift for the first APC waveform (when k = 1). The APC demodulation phase in (3) is constructed based on this strategy. Though the echoes alias in the time domain, they can be separated in the Doppler domain, which will be discussed in Section 3.

2.3. Signal Model and Doppler Spectra

The analytic transmitted signal of Tx k can be expressed as
s k ( τ , t a ) = s ( τ , t a ) φ mod , k ( t a )  
where τ denotes the fast time, ta denotes the slow time, and s ( τ , t a ) is defined as the original signal without the APC modulation phase, which is similar to the echo signal of a monostatic SAR system. For a narrowband transmitted linear frequency modulated (LFM) pulse signal, the modulation form of s ( τ , t a ) can be expressed as
s ( τ , t a ) = w r ( τ ) w a ( t a ) exp ( j π γ τ 2 ) exp ( j 2 π f c τ )  
where w r ( ) is the range window function, w a ( ) is the antenna azimuth pattern modulation, e is the chirp rate, and fc is the fundamental carrier frequency.
For a point scatterer in the terrain, the echoes received by Rx n m in the n-th row and the m-th column can be expressed as
s n m ( τ , t a , X n m ) = k = 1 K σ s ( τ t n m , t a + X n m / 2 v ) φ r e s , k ( t a + X n m / 2 v )  
where σ denotes the backscattering coefficient of the scatterer, v denotes the velocity of the radar planform, Xnm denotes the azimuth coordinate of Rx n m . Here, tnm refers to the propagation time between the radar and the scatterer, which can be calculated by the round-trip slant range from the transmitter k Tx k to the scatterer, and then back to the receiver Rx n m .
After performing an azimuth Fourier transform, the signal in (8) can be expressed in the Doppler domain as
S n m ( τ , f a , X n m ) = k = 1 K S ( τ t n m , f a Δ f d , k ) exp ( j 2 π ( f a Δ f d , k ) X n m / 2 v )  
where S ( τ , f a ) is the Doppler spectrum of s ( τ , t a ) , f a is the azimuth Doppler frequency, and Δ f d , k is the additional Doppler shift frequency shown in (5). Thus, the Doppler spectrum of the k-th APC waveform can be regarded as its corresponding original signal with a Doppler shift by Δ f d , k .
In order to obtain a relatively wide swath, the PRF is set to be slightly higher than the Doppler bandwidth, which will cause the shifted Doppler spectra of the APC waveform echoes to span the neighboring PRF. Thus, after the PRF sampling, the Doppler spectra aliasing will occur, as Figure 5 shows.
When the number of the APC waveforms is 3 (K = 3), the Doppler spectrum of the received signal S n m ( τ , t a , X n m ) is shown in Figure 5. Echo 1, Echo 2 and Echo 3 denote the echoes of the first, second and third APC waveforms, respectively. Figure 5a shows the Doppler spectra before the PRF sampling, where the spectrum spans the neighboring PRF. Figure 5b shows the Doppler spectra after the PRF sampling, where the spectrum aliasing occurs. The aliasing causes the energy accumulation to degrade during the imaging processing, which may make it impossible to get an excellent focused image. Therefore, the Doppler spectrum of each APC waveform echo needs to be separated completely before the imaging processing for the APC–MIMO–SAR.

3. Echo Separation and Reconstruction

As the analysis mentioned above, some new problems occur in the proposed APC–MIMO–SAR. The echo signal received simultaneously aliases in the time domain while the one sampled by PRF aliases in the Doppler domain. To deal with these problems, we use the APC technology and the azimuth DBF processing to separate the waveforms, suppress the Doppler ambiguity and reconstruct the echo signal. The approach is discussed in detail below.
Firstly, in a SAR system, it is explicit that the angle–Doppler relation of azimuth Doppler frequency f a and azimuth squint instantaneous angle θ can be described by
sin θ = λ 2 v f a  
where λ is wavelength.
Figure 6 shows the sketch of the angle–Doppler relation of the received signal of the APC waveforms after being demodulated in the azimuth. The sketch provides a concise visualization of the time–frequency relation of the received signal.
From Figure 6, it becomes clear that different APC echo signals present different additional Doppler shift frequencies due to the APC modulation and demodulation. This provides the possibility to separate the time domain-aliased APC waveforms in the azimuth Doppler domain. The angle–Doppler relation of the echo k can be deduced as
sin θ k ( f a ) = λ 2 v ( f a Δ f d , k )  
Substituting (11) into (9), the signal in the Doppler domain can be rewritten as
S n m ( τ , f a , X n m ) = k = 1 K S ( τ t n m , 2 v sin θ k / λ ) exp ( j 2 π X n m sin θ k / λ )  
In the case that the PRF is slightly larger than the Doppler bandwidth, it is also shown in Figure 6 that the Doppler spectra of Echo 2 and Echo 3 alias due to the PRF sampling. Thus, the signal of each Doppler bin can be regarded as a sum of echoes of K APC waveforms where different echoes correspond to different azimuth squint angles. Furthermore, each azimuth squint angle corresponds to K APC waveforms, which makes the direct application of the azimuth DBF processing impossible.
To deal with this problem in the APC–MIMO–SAR, a spatial–temporal filtering can be used to extract the signal from each angle. Firstly, the whole Doppler band has been divided uniformly into K sub-bands as the vertical dotted lines shown in Figure 6. The i-th ( i = 1 , 2 , , K ) Doppler sub-band can be described by
( 1 2 + i 1 K ) f PRF f a ( 1 2 + i K ) f PRF  
where i denotes the order number of the Doppler sub-band. In each Doppler sub-band, the angle–Doppler relation of each APC waveform has one-to-one correspondence. In this case, the angle–Doppler relation of the echo k in the i-th Doppler sub-band can be revised to
sin θ i , k ( f a ) = λ 2 v ( f a Δ f d , k + M i , k f PRF )  
where M i , k = { 0 , k i 1 , k > i is the Doppler ambiguity number.
Secondly, the spatial–temporal filtering is designed with the azimuth DBF technique [36,37]. The filtering is a weighting operation of MIMO–SAR echo for every Doppler bin of a Doppler sub-band. The wanted echo k in the i-th Doppler sub-band without the Doppler ambiguity can be extracted from each angle as Figure 7 shows, and be written as
S u n a m b , i , k ( τ , f a Δ f d , k + M i , k f PRF ) = S ( τ , f a , X n m ) w ( f a )  
where S u n a m b , i , k ( τ , f a Δ f d , k + M i , k f PRF ) denotes the extracted echo, w ( f a ) denotes the weight vector, and S ( τ , f a , X n m ) denotes the received signal vector. S ( τ , f a , X n m ) is a 1 × K-dimensional vector constructed by the signal from the n-th row receiving sub-apertures, expressed as
S ( τ , f a , X n m ) = [ S n 1 ( τ , f a , X n 1 ) , , S n k ( τ , f a , X n k ) , , S n K ( τ , f a , X n K ) ]  
where the exact expression of S n k ( τ , f a , X n k ) has been shown in (12).
The weight vector can be obtained by solving the following equations
w ( f a ) = A 1 ( f a ) e k  
where superscript −1 denotes the matrix inverse and e k = [ e 1 , , e q , , e M ] T is a unit vector, and e q = k = 1 , e q k = 0 , which means in vector e k only one element is 1 and not 0. In (17), A ( f a ) = [ a 1 ( f a ) , , a k ( f a ) , , a K ( f a ) ] is an M × K-dimensional matrix, and a k ( f a ) is the azimuth sub-aperture array steering vector which can be constructed as
a k ( f a ) = [ exp ( j 2 π X n 1 sin θ k / λ ) , , exp ( j 2 π X n M sin θ k / λ ) ] T  
where superscript T denotes the matrix transposition.
As shown in Figure 7a,b, for every Doppler bin in the second Doppler sub-band, the wanted echo can be extracted by steering the formed beam center to the corresponding angle. It is clear that the echo of K APC waveforms can be separated completely via performing the same process on each Doppler sub-band.
Then, the extracted echo signals from different Doppler sub-bands should be rearranged to reconstruct the whole Doppler band signal, as Figure 7c shows, after compensating the residual phase of the k-th APC waveform in (4).
At last, a 2-D focused APC–MIMO–SAR image of the wanted echo k is obtained via the range matching filtering and azimuth focusing processing, which can be applied to the subsequent applications.
For clarity, a useful flowchart of the proposed APC–MIMO–SAR solution is shown in Figure 8. The main steps can be summarized as the APC technique module, the DBF processing module and imaging processing module.

4. Discussion

In fact, the suggested approach, performing each Doppler frequency with a null-steering in the azimuth to reconstruct the unambiguous SAR signal, is closely connected to what has already been proposed in [35]. However, according to [35], any deviation from the multichannel displaced phase center aperture (DPCA) system or DPCA condition will deteriorate the performance of the null-steering and raise the noise level. In this section, the performance deterioration in the form of signal-to-noise ratio (SNR) needs to be discussed first in case of an unsatisfied DPCA condition.
Analogically, the DPCA condition in the suggested approach is
X n m X n 1 = 2 N v f PRF ( m 1 N + k m ) ,   k m = [ 0 , 1 , 2 , ]  
where X n m denotes the azimuth coordinate of the Rx n m ( m = 2 , 3 , , M ).
As a measure for the variation of SNR caused by the DBF network, the SNR scaling factor Φ D B F [38] can be obtained by
Φ D B F = S N R i n / S N R o u t ( S N R i n / S N R o u t ) DPCA = N k = 1 N E [ | a k ( f a ) | 2 ]  
where SNRin and SNRout denote the SNRs before and after DBF, respectively. The operator E [ · ] represents the mean squared value operator, and a k ( f a ) is the k-th column in the matrix A in (16) in the Doppler frequency domain.
When the DPCA condition is satisfied, the optimum Φ D B F is obtained. However, Φ D B F will worsen with a rising mean squared value of a k ( f a ) due to the increased DPCA deviation.
To ensure a reliable echo separation by DBF, the formed Rx beam should be narrow enough, and thus the azimuth length of each receiving sub-aperture should exceed
L a 2 v K / PRF  
It should be emphasized that the sufficient receiving sub-apertures needed by echo separation require a longer azimuth antenna aperture, which may limit their application to the acquisition of additional phase centers and longer baselines in the azimuth.

5. Simulation Results

In this section, simulations on a point target and distributed targets are carried out to verify the validity of the proposed APC–MIMO–SAR solution based on the APC waveforms and echo separation method by azimuth DBF processing. The main parameters of a MIMO–SAR system are given in Table 1.

5.1. Simulation Results on a Point Target

In this part, a simulation is performed with a point target scene to show the effectiveness of the proposed Doppler spectra separation processing, using the azimuth digital beamforming (DBF) technique.
From the experimental parameters in Table 1, one can notice that the PRF is set to be slightly larger than the Doppler bandwidth, which will cause the Doppler spectra aliasing occur after the sampling with PRF. Figure 9 shows the Doppler spectra of echoes before and after echo separation by azimuth DBF, in which the echoes are received by one of the four sub-apertures, for examplee Rx11 For a better presentation in Figure 9, we assume the antenna pattern in the azimuth is in the rectangular shape. It is clearly observed in Figure 9a that the echoes of the two APC waveforms overlap each other because of the APC shift effect and the PRF sampling. After the Doppler sub-band-dependent azimuth DBF processing, the echoes are well separated, as shown in Figure 9b,c, where the Echo 2 shifts in the Doppler domain observably.
Figure 10 presents the azimuth imaging results of the echoes received by Rx11. Figure 10a shows the focusing effect of the point target before the azimuth DBF processing, while Figure 10b shows the one after the processing. It is obvious that the focusing effect in Figure 10a is worse than that in Figure 10b. By comparing Figure 10a with Figure 10b, one can conclude that the focusing effect of the point target has been improved dramatically after the azimuth DBF processing.

5.2. Simulation Results on Distributed Targets

In this part, a simulation on distributed targets with the main parameters given in Table 1 is performed to show the advantage and the effectiveness of the proposed APC–MIMO–SAR.
Figure 11a shows a picture of the reference terrain scene that provides the input for the scene simulation. Figure 11b,c compare the focused images of the MIMO–SAR obtained via different transmit signals. Figure 11b shows the imaging result obtained via the conventional matched filter processing, when the MIMO–SAR sends an up-chirp signal and a down-chirp waveform simultaneously. From Figure 11b, one can notice that the leaked signal energy from the orthogonal waveform (down-chirp in this case) causes the degradation of the focusing performance obviously. Figure 11c shows the imaging result obtained via the DBF processing in the azimuth and the conventional matched filter processing, when the MIMO–SAR transmits the APC waveforms simultaneously. Figure 11c clearly demonstrates the good suppression of the mutual interference of the radar echoes by the proposed APC–MIMO–SAR solution. By comparing Figure 11b with Figure 11c, it is apparent that the quality of the whole image has been improved, and the focusing effect is better.
For a detailed imaging quality comparison, the image entropy and the contrast of the imaging results shown in Figure 11b,c are listed in Table 2. From Table 2, it is clear that the image entropy has decreased and the image contrast has increased in the proposed MIMO–SAR solution.
As the phase preservation of echo separation and reconstruction is very important in the multibaseline SAR system, the interferometric phase of the distributed targets has been analyzed in this section. Figure 12a shows the interference phase of the reference image (Figure 11a) and the complex image obtained by the up- and down-chirp waveforms (Figure 11b). The interference phase of the middle azimuth cell of Figure 12a is clearly shown in Figure 12c, which is up to 2 radians. From Figure 12a,c, it can be seen that the MIMO–SAR modulated by up- and down-chirp waveforms has poor phase-preserving performance. Figure 12b shows the interference phase of Figure 11a and the complex image obtained by the proposed APC–MIMO–SAR (Figure 11c). The interference phase of the middle azimuth cell of Figure 12b is close to 0 radians, as Figure 12d shows. From Figure 12b,d, one can see that APC–MIMO–SAR has better phase-preserving performance.
Based on the comparison of amplitude and phase shown in Figure 11 and Figure 12, the conclusion can be drawn that the proposed APC–MIMO–SAR imaging results are much closer to the ideal ones in both amplitude and phase, in contrast to MIMO–SAR transmitting up- and down-chirp waveforms. From this perspective, the proposed APC–MIMO–SAR has broader application prospects.

6. Conclusions

A novel MIMO–SAR solution, based on the principle of the APC technique combined with azimuth DBF, has been proposed in this paper. The proposed APC–MIMO–SAR provides a longer baseline in elevation, which will contribute more to the processing and application of multibaseline InSAR and three-dimensional SAR imaging. Applying the DBF processing in the azimuth as a spatial–temporal filtering to separate the aliasing echoes sustains a lesser PRF sampling, which will be beneficial to obtain a wider swath. Thus, the proposed APC–MIMO–SAR can be used to exploit the potentials of MIMO–SAR without performance degradation.
Furthermore, the good adaptability featured by APC waveforms allows them to combine with other waveforms easily, such as the multidimensional encoding waveforms or the orthogonal frequency division multiplexing (OFDM) waveforms. However, it is worth noting that the echo separation by azimuth DBF processing requires sufficient spatial degrees of freedom, which may lead to a longer azimuth aperture.

Author Contributions

Conceptualization, F.Z.; Formal analysis, F.Z. and J.Z.; Investigation, J.A.; Methodology, F.Z.; Supervision, M.X.; Validation, F.Z. and Z.D.; Writing—original draft, F.Z.; Writing—review & editing, F.Z.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 61701156, the Postdoctoral Science Foundation of China, grant number 2016M592045, and the Anhui Provincial Natural Science Foundation, grant number 1608085QF142. Ai’s work was supported by the National Natural Science Foundation of China, grant number 61701157 and the Anhui Provincial Natural Science Foundation, grant number 1808085QF206.

Acknowledgments

The authors will thank all the reviewers and editors for their great help and useful suggestions. Thank Sishi Wu for her help for polishing the English writing.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Maori, D.; Sikaneta, I.; Klare, J.; Gierull, C. MIMO SAR Processing for Multichannel High-Resolution Wide-Swath Radars. IEEE Trans. Geosci. Remote Sens. 2014, 52, 5034–5055. [Google Scholar] [CrossRef]
  2. Kriege, G.; Gebert, N.; Moreira, A. Multidimensional waveform encoding: A new digital beamforming technique for synthetic aperture radar remote sensing. IEEE Trans. Geosci. Remote Sens. 2008, 46, 31–46. [Google Scholar] [CrossRef] [Green Version]
  3. Jing, G.; Sun, G.; Xia, X.; Xing, M.; Bao, Z. A Novel Two-Step Approach of Error Estimation for Stepped-Frequency MIMO–SAR. IEEE Geosci. Remote Sens. Lett. 2017, 14, 2290–2294. [Google Scholar] [CrossRef]
  4. Fornaro, G.; Serafino, F.; Soldovieri, F. Three-dimensional focusing with multipass SAR data. IEEE Trans. Geosci. Remote Sens. 2003, 41, 507–517. [Google Scholar] [CrossRef]
  5. Krieger, G.; Hajnsek, I.; Papathanassiou, K.; Younis, M.; Moreira, A. Interferometric synthetic aperture radar (SAR) missions employing formation flying. Proc. IEEE 2010, 98, 816–843. [Google Scholar] [CrossRef] [Green Version]
  6. Klare, J.; Maori, D.C.; Brenner, A.; Ender, J. Image quality analysis of the vibrating sparse MIMO antenna array of the airborne 3D imaging radar ARTINO. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Barcelona, Spain, 23–28 July 2007. [Google Scholar]
  7. Miralles, E.; Multerer, T.; Ganis, A.; Schoenlinner, B.; Prechtel, U.; Meusling, A.; Mietzner, J.; Weckerle, C.; Esteban, H.; Vossiek, M.; et al. Multifunctional and compact 3D FMCW MIMO radar system with rectangular array for medium-range applications. IEEE Aerosp. Electron. Syst. Mag. 2018, 33, 46–54. [Google Scholar] [CrossRef]
  8. Xu, J.; Dai, X.Z.; Xia, X.G.; Wang, L.B.; Yu, J.; Peng, Y.N. Optimizations of multisite radar system with MIMO Radars for target detection. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 2329–2343. [Google Scholar] [CrossRef]
  9. Monte, L.L.; Himed, B.; Corigliano, T.; Baker, C. Performance analysis of time division and code division waveforms in co-located MIMO. In Proceedings of the 2015 IEEE Radar Conference, Arlington, VA, USA, 10–15 May 2015. [Google Scholar]
  10. Boute, R. On the equivalence of time-division and frequency-division multiplexing. IEEE Trans. Commun. 1985, 33, 97–99. [Google Scholar] [CrossRef]
  11. Lutz, S.; Walter, T. Lens based 77 GHz TDM MIMO radar sensor for angular estimation in multitarget environments. In Proceedings of the 2013 European Radar Conference, Nuremberg, Germany, 9–11 October 2013. [Google Scholar]
  12. Mishra, K.V.; Eldar, Y.C.; Shoshan, E.; Namer, M.; Meltsin, M. A cognitive sub-Nyquist MIMO radar prototype. Arxiv, 2018; arXiv:1807.09126. [Google Scholar]
  13. Xu, J.; Liao, G.; Zhu, S.; Huang, L.; So, H.C. Joint range and angle estimation using MIMO radar with frequency diverse array. IEEE Trans. Signal Process. 2015, 63, 3396–3410. [Google Scholar] [CrossRef]
  14. Deng, H. Polyphase code design for orthogonal netted radar systems. IEEE Trans. Signal Process. 2004, 52, 3126–3135. [Google Scholar] [CrossRef]
  15. Zou, B.; Dong, Z.; Liang, D.N. Design and performance analysis of orthogonal coding signal in MIMO–SAR. Sci. Chin. Inf. Sci. 2011, 54, 1723–1737. [Google Scholar] [CrossRef]
  16. Xia, X.; Zhang, T.; Kong, L. MIMO OFDM radar IRCI free range reconstruction with sufficient cyclic prefix. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 2279–2293. [Google Scholar] [CrossRef]
  17. Wang, W.Q. Space-time coding MIMO–OFDM SAR for high-resolution imaging. IEEE Trans. Geosci. Remote Sens. 2011, 49, 3094–3104. [Google Scholar] [CrossRef]
  18. Kim, J.; Younis, M.; Moreira, A.; Wiesbeck, W. A novel OFDM chirp waveform scheme for use of multiple transmitters in SAR. IEEE Geosci. Remote Sens. Lett. 2013, 10, 568–572. [Google Scholar] [CrossRef] [Green Version]
  19. Meng, C.; Xu, J.; Liu, F.; Long, T.; Yang, J.; Peng, Y. Multiple input and multiple output synthetic aperture radar multiple waveform separation based on oblique projection in same frequency coverage. IET Radar Sonar Navig. 2015, 9, 1088–1096. [Google Scholar] [CrossRef]
  20. Wang, J.; Chen, L.; Liang, X.; Ding, C.; Li, K. Implementation of the OFDM Chirp Waveform on MIMO SAR Systems. IEEE Trans. Geosci. Remote Sens. 2015, 53, 5218–5228. [Google Scholar] [CrossRef]
  21. Huang, Y.; Liao, G.; Xu, J.; Li, J.; Yang, D. GMTI and parameter estimation for MIMO SAR system via fast interferometry RPCA method. IEEE Trans. Geosci. Remote Sens. 2018, 56, 1774–1787. [Google Scholar] [CrossRef]
  22. Christallini, D.; Patina, D.; Lombardo, P. Exploiting MIMO SAR potentialities with efficient cross-track constellation configurations for improved range resolution. IEEE Trans. Geosci. Remote Sens. 2011, 49, 38–52. [Google Scholar] [CrossRef]
  23. Meng, C.; Xu, J.; Xia, X.; Long, T.; Mao, E.; Yang, J.; Peng, Y. MIMO-SAR waveform separation based on inter-pulse phase modulation and range-Doppler decouple filtering. IET Electron. Lett. 2013, 49, 420–421. [Google Scholar] [CrossRef]
  24. Younis, M.; Fischer, C.; Wiesbeck, W. Digital beamforming in SAR systems. IEEE Trans. Geosci. Remote Sens. 2003, 41, 1735–1739. [Google Scholar] [CrossRef]
  25. Krieger, G. MIMO-SAR: Opportunities and pitfalls. IEEE Trans. Geosci. Remote Sens. 2014, 52, 2628–2645. [Google Scholar] [CrossRef]
  26. Bellettini, A.; Pinto, M.A. Theoretical Accuracy of Synthetic Aperture Sonar Micronavigation Using a Displaced Phase Center Antenna. IEEE J. Ocean. Eng. 2002, 27, 780–789. [Google Scholar] [CrossRef]
  27. Prats, P.; Scheiber, R.; Mittermayer, J.; Meta, A.; Moreira, A. Processing of Sliding Spotlight and TOPS SAR Data Using Baseband Azimuth Scaling. IEEE Trans. Geosci. Remote Sens. 2010, 48, 770–780. [Google Scholar] [CrossRef] [Green Version]
  28. Sun, G.; Xing, M.; Wang, Y.; Wu, Y.; Wu, Y.; Bao, Z. Sliding spotlight and TOPS SAR data processing without subaperture. IEEE Geosci. Remote Sens. Lett. 2011, 8, 1036–1040. [Google Scholar] [CrossRef]
  29. Krieger, G.; Moreira, A. Potentials of digital beamforming in bi- and multistatic SAR. In Proceedings of the 2003 IEEE International Geoscience and Remote Sensing Symposium, Toulouse, France, 21–25 July 2003. [Google Scholar]
  30. Callaghan, G.D.; Longstaff, I.D. Wide swath space-borne SAR using a quad-element array. IEEE Proc. Radar Sonar Navig. 1999, 146, 159–165. [Google Scholar] [CrossRef]
  31. Bordoni, F.; Younis, M.; Krieger, G. Ambiguity suppression by azimuth phase coding in multichannel SAR systems. IEEE Trans. Geosci. Remote Sens. 2012, 50, 617–629. [Google Scholar] [CrossRef] [Green Version]
  32. Sauer, S.; Ferro-Famil, L.; Reigber, A.; Pottier, E. Polarimetric Dual-Baseline InSAR Building Height Estimation at L-Band. IEEE Geosci. Remote Sens. Lett. 2009, 6, 408–412. [Google Scholar] [CrossRef]
  33. Sauer, S.; Ferro-Famil, L.; Reigber, A.; Pottier, E. Three-Dimensional Imaging and Scattering Mechanism Estimation Over Urban Scenes Using Dual-Baseline Polarimetric InSAR Observations at L-Band. IEEE Trans. Geosci. Remote Sens. 2011, 49, 4616–4629. [Google Scholar] [CrossRef] [Green Version]
  34. Currie, A.; Brown, M.A. Wide-swath SAR. IET Proc. Radar Signal Process. 1992, 139, 122–135. [Google Scholar] [CrossRef]
  35. Krieger, G.; Gebert, N.; Moreira, A. Unambiguous SAR signal reconstruction from nonuniform displaced phase center sampling. IEEE Geosci. Remote Sens. Lett. 2004, 1, 260–264. [Google Scholar] [CrossRef]
  36. Li, Z.F.; Wang, H.Y.; Su, T.; Bao, Z. Generation of wide-swath and high-resolution SAR images from multichannel small spaceborne SAR systems. IEEE Geosci. Remote Sens. Lett. 2005, 2, 82–86. [Google Scholar] [CrossRef]
  37. Sun, G.C.; Xing, M.; Xia, X.G.; Wu, Y.; Huang, P.; Wu, Y.; Bao, Z. Multichannel Full-Aperture Azimuth Processing for Beam Steering SAR. IEEE Trans. Geosci. Remote Sens. 2013, 51, 4761–4778. [Google Scholar] [CrossRef]
  38. Gebert, N.; Krieger, G.; Moreira, A. Digital beamforming on receive: Techniques and optimization strategies for high-resolution Wide-Swath SAR Imaging. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 564–592. [Google Scholar] [CrossRef]
Figure 1. One of the possible architectures of the transmitting and receiving aperture (K = 3, N = 3, M = 3).
Figure 1. One of the possible architectures of the transmitting and receiving aperture (K = 3, N = 3, M = 3).
Sensors 18 03374 g001
Figure 2. Sketch of the elevation- and azimuth-displaced phase centers.
Figure 2. Sketch of the elevation- and azimuth-displaced phase centers.
Sensors 18 03374 g002
Figure 3. Sketch of the modulated azimuth phase coding (APC) waveforms.
Figure 3. Sketch of the modulated azimuth phase coding (APC) waveforms.
Sensors 18 03374 g003
Figure 4. Sketch of the demodulated APC waveforms.
Figure 4. Sketch of the demodulated APC waveforms.
Sensors 18 03374 g004
Figure 5. Doppler spectra of echoes of APC waveforms (K = 3) before (a) and after (b) sampling with f P R F .
Figure 5. Doppler spectra of echoes of APC waveforms (K = 3) before (a) and after (b) sampling with f P R F .
Sensors 18 03374 g005
Figure 6. Sketch of the angle–Doppler relation (K = 3).
Figure 6. Sketch of the angle–Doppler relation (K = 3).
Sensors 18 03374 g006
Figure 7. Azimuth digital beamforming for echo (a) spatial–temporal filtering, (b) separation and (c) reconstruction.
Figure 7. Azimuth digital beamforming for echo (a) spatial–temporal filtering, (b) separation and (c) reconstruction.
Sensors 18 03374 g007
Figure 8. Flow chart of the proposed APC–MIMO–SAR solution.
Figure 8. Flow chart of the proposed APC–MIMO–SAR solution.
Sensors 18 03374 g008
Figure 9. Doppler spectra of echoes of APC waveforms. (a) The overlapped Doppler; (b) the separated Doppler for Echo 1; (c) the separated Doppler for Echo 2.
Figure 9. Doppler spectra of echoes of APC waveforms. (a) The overlapped Doppler; (b) the separated Doppler for Echo 1; (c) the separated Doppler for Echo 2.
Sensors 18 03374 g009
Figure 10. Azimuth imaging results (a) before and (b) after DBF.
Figure 10. Azimuth imaging results (a) before and (b) after DBF.
Sensors 18 03374 g010
Figure 11. Simulation results of MIMO–SAR imaging varying with the transmitted waveforms. (a) Reference terrain scene. (b) Imaging result of up- and down-chirp waveforms. (c) Imaging result of the proposed APC waveforms.
Figure 11. Simulation results of MIMO–SAR imaging varying with the transmitted waveforms. (a) Reference terrain scene. (b) Imaging result of up- and down-chirp waveforms. (c) Imaging result of the proposed APC waveforms.
Sensors 18 03374 g011aSensors 18 03374 g011b
Figure 12. Simulation results of the interferometric phase between the reference terrain scene and (a) up- and down-chirp waveform MIMO–SAR and (b) the proposed APC–MIMO–SAR. (c) Middle azimuth cell of (b). (d) Middle azimuth cell of (d).
Figure 12. Simulation results of the interferometric phase between the reference terrain scene and (a) up- and down-chirp waveform MIMO–SAR and (b) the proposed APC–MIMO–SAR. (c) Middle azimuth cell of (b). (d) Middle azimuth cell of (d).
Sensors 18 03374 g012
Table 1. Main system parameters of simulations.
Table 1. Main system parameters of simulations.
ParametersQuantityParametersQuantity
Wavelength0.03 mSignal bandwidth100 MHz
Center frequency10 GHzAzimuth Doppler3750 Hz
Sensor height600 kmPRF4500 Hz
Platform velocity7500 m/sIncident angle45°
Transmitting aperture number2APC waveform number2
Receiving sub-aperture number4Azimuth resolution2 m
Table 2. The image entropy and the contrast of the imaging results.
Table 2. The image entropy and the contrast of the imaging results.
ParametersEntropyContrast
Figure 11b7.107336.6054
Figure 11c6.8526115.3035

Share and Cite

MDPI and ACS Style

Zhou, F.; Ai, J.; Dong, Z.; Zhang, J.; Xing, M. A Novel MIMO–SAR Solution Based on Azimuth Phase Coding Waveforms and Digital Beamforming. Sensors 2018, 18, 3374. https://doi.org/10.3390/s18103374

AMA Style

Zhou F, Ai J, Dong Z, Zhang J, Xing M. A Novel MIMO–SAR Solution Based on Azimuth Phase Coding Waveforms and Digital Beamforming. Sensors. 2018; 18(10):3374. https://doi.org/10.3390/s18103374

Chicago/Turabian Style

Zhou, Fang, Jiaqiu Ai, Zhangyu Dong, Jiajia Zhang, and Mengdao Xing. 2018. "A Novel MIMO–SAR Solution Based on Azimuth Phase Coding Waveforms and Digital Beamforming" Sensors 18, no. 10: 3374. https://doi.org/10.3390/s18103374

APA Style

Zhou, F., Ai, J., Dong, Z., Zhang, J., & Xing, M. (2018). A Novel MIMO–SAR Solution Based on Azimuth Phase Coding Waveforms and Digital Beamforming. Sensors, 18(10), 3374. https://doi.org/10.3390/s18103374

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