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 (), while the nm-th receiving sub-aperture, in the n-th row and the m-th column, is represented as (, ). 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
, the second one marked with “×” corresponds to transmitter
, and the third one marked with “△” corresponds to transmitter
. 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 (
) transmit waveform (
) is given by
where l denotes the sequence number of the azimuth pulse. Apparently, l can be expressed as
, where
denotes the pulse repetition frequency (PRF) and
denotes the azimuth slow time. Thus, substituting
into (1), the APC modulation phase can be rewritten as
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
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
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
can be written as
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
can be expressed as
where τ denotes the fast time, t
a denotes the slow time, and
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
can be expressed as
where
is the range window function,
is the antenna azimuth pattern modulation, e is the chirp rate, and f
c is the fundamental carrier frequency.
For a point scatterer in the terrain, the echoes received by
in the
n-th row and the
m-th column can be expressed as
where σ denotes the backscattering coefficient of the scatterer,
v denotes the velocity of the radar planform,
Xnm denotes the azimuth coordinate of
. 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 to the scatterer, and then back to the receiver
.
After performing an azimuth Fourier transform, the signal in (8) can be expressed in the Doppler domain as
where
is the Doppler spectrum of
,
is the azimuth Doppler frequency, and
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
.
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
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
and azimuth squint instantaneous angle
can be described by
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
Substituting (11) into (9), the signal in the Doppler domain can be rewritten as
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 (
) Doppler sub-band can be described by
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
where
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
where
denotes the extracted echo,
denotes the weight vector, and
denotes the received signal vector.
is a 1 × K-dimensional vector constructed by the signal from the
n-th row receiving sub-apertures, expressed as
where the exact expression of
has been shown in (12).
The weight vector can be obtained by solving the following equations
where superscript −1 denotes the matrix inverse and
is a unit vector, and
,
, which means in vector
only one element is 1 and not 0. In (17),
is an M × K-dimensional matrix, and
is the azimuth sub-aperture array steering vector which can be constructed as
where superscript
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
where
denotes the azimuth coordinate of the
(
).
As a measure for the variation of SNR caused by the DBF network, the SNR scaling factor
[
38] can be obtained by
where SNR
in and SNR
out denote the SNRs before and after DBF, respectively. The operator
represents the mean squared value operator, and
is the k-th column in the matrix
in (16) in the Doppler frequency domain.
When the DPCA condition is satisfied, the optimum is obtained. However, will worsen with a rising mean squared value of 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
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.
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.