Content area
The performance of distributed small satellite synthetic aperture radar systems degrades significantly due to the unavoidable array errors, including gain, phase, and position errors, in real operating scenarios. In the conventional method proposed in (IEEE T Aero. Elec. Sys. 42:436-451, 2006), the spectrum components within one Doppler bin are considered as calibration sources. However, it is found in this article that the gain error estimation and the position error estimation in the conventional method can interact with each other. The conventional method may converge to suboptimal solutions in large position errors since it requires the joint iteration between gain-phase error estimation and position error estimation. In addition, it is also found that phase errors can be estimated well regardless of position errors when the zero Doppler bin is chosen. In this article, we propose a method obtained by modifying the conventional one, based on these two observations. In this modified method, gain errors are firstly estimated and compensated, which eliminates the interaction between gain error estimation and position error estimation. Then, by using the zero Doppler bin data, the phase error estimation can be performed well independent of position errors. Finally, position errors are estimated based on the Taylor-series expansion. Meanwhile, the joint iteration between gain-phase error estimation and position error estimation is not required. Therefore, the problem of suboptimal convergence, which occurs in the conventional method, can be avoided with low computational method. The modified method has merits of faster convergence and lower estimation error compared to the conventional one. Theoretical analysis and computer simulation results verified the effectiveness of the modified method.[PUBLICATION ABSTRACT]
Xu et al. EURASIP Journal on Advances in Signal Processing 2013, 2013:31 http://asp.eurasipjournals.com/content/2013/1/31
RESEARCH Open Access
Self-calibration method without joint iteration for distributed small satellite SAR systems
Qing Xu*, Guisheng Liao, Aifei Liu and Juan Zhang
Abstract
The performance of distributed small satellite synthetic aperture radar systems degrades significantly due to the unavoidable array errors, including gain, phase, and position errors, in real operating scenarios. In the conventional method proposed in (IEEE T Aero. Elec. Sys. 42:436451, 2006), the spectrum components within one Doppler bin are considered as calibration sources. However, it is found in this article that the gain error estimation and the position error estimation in the conventional method can interact with each other. The conventional method may converge to suboptimal solutions in large position errors since it requires the joint iteration between gain-phase error estimation and position error estimation. In addition, it is also found that phase errors can be estimated well regardless of position errors when the zero Doppler bin is chosen. In this article, we propose a method obtained by modifying the conventional one, based on these two observations. In this modified method, gain errors are firstly estimated and compensated, which eliminates the interaction between gain error estimation and position error estimation. Then, by using the zero Doppler bin data, the phase error estimation can be performed well independent of position errors. Finally, position errors are estimated based on the Taylor-series expansion. Meanwhile, the joint iteration between gain-phase error estimation and position error estimation is not required. Therefore, the problem of suboptimal convergence, which occurs in the conventional method, can be avoided with low computational method. The modified method has merits of faster convergence and lower estimation error compared to the conventional one. Theoretical analysis and computer simulation results verified the effectivenessof the modified method.
Keywords: Array signal processing, Distributed small satellite synthetic aperture radar (DSS-SAR) systems, Error estimation
1. Introduction
With the development of spaceborne synthetic aperture radar (SAR) systems, the functions, such as SAR image, ground moving target indication (GMTI), and SAR interferometry (InSAR), have been well performed [1-4]. In the conventional spaceborne SAR systems, large antennas are required due to the minimum antenna area constraint [5].However, it leads to the failure for the systems to obtain an image of wide area since the illumination area is inversely related to the aperture size of the antennas [6].Distributed small satellite synthetic aperture radar (DSSSAR) systems [6-12] were developed to deal with the problem and have received considerable attention in recent years. In DSS-SAR systems, several small satellites
move in a special orbital configuration and function as a single virtual satellite. A small antenna covering wide area is placed on each satellite with the total antenna area constituted by all the small antennas satisfying the minimum antenna area requirement. Due to the flying formation of DSS-SAR systems, along-track baselines and across-track baselines may exist synchronously. And the across-track baseline is needed for terrain height estimation [7] while the along-track baseline is suitable for the function of SAR and GMTI [8]. In [9], the SAR train configuration is analyzed and many scholars focus their attention on SAR imaging or GMTI based on this configuration [10,11]. However, for the echo received by each small satellite, range or azimuth (Doppler) ambiguities will occur due to the use of small antennas. In order to image a wide swath unambiguously, the echoes of small antennas should be combined coherently in the DSS-SAR system.
* Correspondence: mailto:[email protected]
Web End [email protected] National Laboratory of Radar Signal Processing, Xidian University, Xian 710071, China
2013 Xu et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0
Web End =http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Xu et al. EURASIP Journal on Advances in Signal Processing 2013, 2013:31 Page 2 of 10 http://asp.eurasipjournals.com/content/2013/1/31
In [10,11], the low pulse repetition frequency (PRF) is chosen to avoid range ambiguity which results in Doppler ambiguity, and the approaches are also given to suppress Doppler ambiguity. However, their excellent performance is critically dependent on the knowledge of the array manifold (parameterized by many parameters, such as angle of arrivals, and antennas position information). The array manifold is used in coherent combination of small antennas echoes to suppress Doppler ambiguity [10]. In practice, there always exist various perturbations in the array manifold which are always called as array errors, such as gain, phase, and position errors. Since the array manifold cannot be exactly obtained because of the existence of array errors, the performance of Doppler ambiguity suppression [11,12] can significantly be degraded. Therefore, it is necessary to estimate and calibrate array errors prior to carry out Doppler ambiguity suppression. In [12], the array error calibration method has also been discussed.In array error estimation of DSS-SAR systems, the existing array calibration methods [13-18] can be applied to the systems only when proper calibration sources are chosen. In [12], spectrum components within one Doppler bin are used as calibration sources with known directions which are called virtual calibration sources. Since the number of spectrum components is more than one due to Doppler ambiguity, it is possible to apply array calibration methods with more than one calibration source [15-18] to DSS-SAR systems. In [12], a two-step iterative auto-calibration method is presented to estimate gain-phase and position errors. In the first step, assuming that position errors are known, gain-phase errors are estimated by using the method in [16]. In the second step, based on the gain-phase error estimated in the former step, position errors can be obtained by the least squared method [17]. These two steps should be iterated alternatively to obtain final solutions. For convenience, the array error estimation method in [12] is named as the conventional method (the comparison in this article is limited in terms of array error estimation method).
In the conventional method, two kinds of errors, the gain-phase error and the position error, are respectively estimated under the assumption that other kinds of errors are known. The inherent relationship among gain, phase, and position errors is not considered and analyzed in the conventional method. In this article, by studying of the conventional method, the following two aspects are observed. First, gain error estimation and position error estimation can affect each other, which will influence the convergence rate. And the conventional method may even suffer from suboptimal convergence in large position errors. Second, if spectrum components within the zero Doppler bin are used as calibration sources to estimate the errors, phase error estimation can be performed independent of position
errors. Based on the above two aspects, a modified array error estimation method is proposed in this article. First, in order to eliminate the interaction between the gain and position error estimations, gain errors are first estimated and compensated. Then, phase errors are estimated by using the spectrum components within the zero Doppler bin as calibration sources. Finally, position errors are estimated based on Taylor-series expansion. However, since Taylor-series expansion causes approximation errors, position error estimation should be iterated in order to obtain higher estimation accuracy. In comparison with the conventional method, the modified method can avoid the iteration between gain-phase error estimation and position error estimation, which guarantees that it can converge to optimal solutions with lower computational load and fast convergence speed. Simulation results verify that the modified method performs better than the conventional one.
The remainder of the article is arranged as follows. In Section 2, the modified method for DSS-SAR systems is described in detail. The performance of the modified method is verified by using Monte Carlo simulations in Section 3. The article ends with some conclusions given in Section 4.
2. Modified method for DSS-SAR systems
2.1. Signal model
In this article, the signal model of DSS-SAR systems will be introduced briefly as follows [12]. The SAR train configuration in which all the satellites are arranged along the X-axis as an array is given in Figure 1. The (X, Y, Z) direction is referred to as (along-track, cross-track, and radial). The array operates in the side-looking strip mode and all satellites have an identical along-track velocity, denoted by vs. , , and represent azimuth angle, incidence angle, and cone angle, respectively. Azimuth angle (, fd) can be expressed as (, fd) = o + (, fd), where
Z (up)
Y (left)
O
v
s
X (front)
Figure 1 Configuration of the DSS-SAR array.
Xu et al. EURASIP Journal on Advances in Signal Processing 2013, 2013:31 Page 3 of 10 http://asp.eurasipjournals.com/content/2013/1/31
o and (, fd) are the azimuth angle and the offset from the beam center, respectively. denotes the fast time corresponding to the range bin and fd denotes the
Doppler frequency which corresponds to the slow time.
True coordinates, measured coordinates, and corresponding position errors of the mth satellite are denoted by (xm, ym,
zm), (xmo, ymo, zmo), and (xm, ym, zm), respectively. The relationship among them is
xm; ym; zm xmo; ymo; zmo
xm; ym; zm
;
1 where ymo = 0, zmo = 0. The first satellite is taken as reference, i.e., (x1, y1, z1) = (0, 0, 0) is imposed. Then the clutter echo in the range-Doppler domain received by the mth satellite after demodulation in the strip mode and broadside geometry can be written as [12]
Sm ; fd gme jm
() is a linear frequency modulated signal. G(fd) denotes the Fourier transform of the function which represents the antenna pattern and other time-variant characters (identical to all receiving antennas). rm(x, y, z, fd) is the
slant range from the mth satellite to the ground cell,
rm x; y; z; fd
y2 z2
p
2
f 2d
p
y2 z2
8v2s
6
ymcos ; fd
sin
zmcos
:
Since ymcos(, fd)sin() and zmcos() can be neglected in rm(x, y, z, fd), it can be obtained that r1(x, y, z, fd) rM(x, y, z, fd) r(x, y, z, fd). So, as H1(, fd) HM(, fd) H(, fd).
Based on the previous analysis, (2) can be rewritten as
Sm ; fd gme j
0m
XIi Ie j4 dm ;fdifr Hm ; fd ifr nm ; fd ;
XIi Iej4d0mi Hi nm ; fd 7
where
0m m
4
ym sin
zm cos
2
where gm and 1 are the gain and phase error of the mth array element relative to the first array element (as such, g1 = 1 and 1 = 0). denotes the wavelength of the carrier.
nm(, fd) is the additive white Gaussian noise. 2I + 1 is the number of spectrum components within one range-Doppler bin. The spectrum components have different azimuth angles (, fd + ifr) and the same incidence angles () which can be expressed as follows:
fd ifr
2vs
sin ; fd ifr
;
d0mi xmo xm
sini sin
2vs fd ifr
sin
; i I; . . . ; I: 3 In the following, for simplicity, ( )i is used to denote the variable associated with (fd + ifr), such as i denotes (, fd + ifr). dmi = xm sini sin() + ym cosi sin() + zm cos(). In the side-looking strip mode, o = 0 and then i = i. Moreover, i is very small (|i| 0.43 for a small antenna with the azimuth length of 2 m at X-band) in practice. Thus, cosi = cosi 1 can be obtained with small angle approximation. Since ymo = zmo = 0,
dim xmo xm
sinisin
ymo ym
cosisin
zmo zm
cos
xmo xm
sini sin
ym sin
zm cos
4
Hm(, fd) is the complex envelope of the clutter echo in the range-Doppler domain,
Hm ; fd x; y
h
xmo xm
:
It is found that the contributions of position errors to array outputs can be separated into two parts. One part associated with cross-track and radial position errors is fixed for all clutter echoes and can be regarded as array phase errors. The other part related to along-track position errors changes with spectrum components. For simplicity, the total phase error 0m is called the phase error and the along-track position error is named as the position error in the following. The unknown errors considered here are the position error xm (m = 2,. . .,M) and the gain-phase error gmej0m (m = 2,. . .,M). Using vector notation, Equation (7) can be rewritten as follows
S AH n 8
where S = [S1,. . ., SM]T, diag 1; . . . ; gMej0
M
n o
, A =
[aI,. . .,aI], ai 1; . . . ; ej4
h iT
d0iM
2rm x; y; z; fd
c
G fd
e j0 x;y;z;fddxdy
; H = [HI,. . .,HI]T, ( )T
denotes the transpose operation.
The covariance matrix of S is denoted as RSS. Denote
eigenvalues and corresponding eigenvectors of the covariance matrix RSS with m (listed in descending order)
and um (m = 1,. . ., M). Each column of A is orthogonal to the matrix U = [u2I+2,. . .,uM]. And the orthogonality can be used to estimate the gain-phase and position errors.
5
where (x, y) is the complex surface scattering coeffi
cient. 0 x; y; z; fd 4
y2 z2
p
p
2vs2
2 2fd xvs . h
y2z
Xu et al. EURASIP Journal on Advances in Signal Processing 2013, 2013:31 Page 4 of 10 http://asp.eurasipjournals.com/content/2013/1/31
2.2. The conventional method
According to [12], the cost function to estimate gain-phase errors is
Jp
XIi I jjUHaijj2; 9
where ( )H represents the conjugate transpose operation. The cost function is given based on the MUSIC approach. If the true and ai are obtained, Jp should be zero. Minimizing Jp, the estimation of and ai can be achieved. Thus, ( = diag{}), gain-phase errors, can be obtained by
^ Q 1w=wT Q 1w; 10
where
w 1; 0; . . . ; 0
T; 11
Q
XIi IDi
H
where E = [eIT. . .eIT]T, F = [fIT. . .fIT]T, ei = UHaoi, fi = UHai. The position error X is the least squared solution of (19)
X real FFH 1FH
E
n o
20
where real () takes the real part of the matrix. The phase error introduced by Taylor-series expansion above is within 4 at X-band under the assumption that xm
20 cm. And higher accuracy can be achieved after several iterations. However, from computer simulations, it is found that if xm > 30 cm, the estimate difference of position errors will be larger than 0.008 m. The conventional method can be summarized as joint iteration between the following two steps.
Step 1 Choose some Doppler bin to obtain gain-
phase errors based on (10).
Step 2 Estimate position errors based on (20) and
compensate xmo by xm.
2.3. Formulation of the modified method
Lets consider the gain error estimation of the conventional method.
Based on (12) and (13), we obtain
Q
XIi ICi
H
UUHDi; 12
13
To estimate position errors, ej4d0im (m = 1,. . .,M) is expanded using Taylor-series based on (3):
e
j 4
d0ime
Di diag ai
:
xm sini sin
e
j 4
xmo sini sin
1 j
4
xm
j 2
vs fd ifr
xmo 1 j
Ci
UUH
( )
UUH
21
Ci
2vs fd ifr
XIi ICi
H
14
Then, based on (3) and (14), ai can be rewritten as
ai aio aiX 15
where
aio 1; ej2
vs fdifr
where Ci
5, denotes dot product
(i.e., element-wise multiplication). With the definition of
Zi = (Ci)H Ci and Z
XIi IZi
2
4
Di11 DiMM Di11 DiMM
3
xMo
h iT
x2o ; . . . ; ej2vs fdifr
16
XIi ICi
H
Ci , (21)
ai diag
0; j 2
vs fd ifr
xMo
17
X 1; x2; . . . ; xM
T 18
Therefore, the cost function (9) can be modified as
Jc jjE FXjj2 19
x2o ; . . . ; j 2
vs fd ifr
can be rewritten as
Q Z UUH
ej2
vs fdifr
ej2
vs fdifr
22
The element in the kth row and the lth column of the matrix Zi and Z can be obtained, respectively, as
Zikl Dikk
Di
ll
ej4
vs fdifr
xloxl xko xk
: 23
Xu et al. EURASIP Journal on Advances in Signal Processing 2013, 2013:31 Page 5 of 10 http://asp.eurasipjournals.com/content/2013/1/31
XIi Iej 4vs ifr xlo xko
xl xk
ai
kk
gkej
0k kej
k
Zkl
XIi IZikl ej 4vs fd xlo xl xko xk
^
gk gk
ej
^
0k 0k
kej
k
^
gk ej ^
0k k ^
gk gk ^
gk ej k
0k
^
gk ej ^
0k ^
k ej ^
k
30
e
j 4
vs fd xlo xl xko xk
where ^
gk and gk denote estimated values and differences of gain errors, ^
0k and 0k represent estimated values and differences of phase errors, respectively. Then^
k k^
gk gk
=^
gk and ^
!
k k 0k in
stead of k and k will be obtained if (29) is used to estimate position errors. Thus, gk will affect the position error estimation accuracy. Considering the fact that position errors have an effect on the gain error estimation during the first step while gain errors affect position error estimation during the second step, estimating differences of gain and position errors using (27) and(20) iteratively in the conventional method could converge to suboptimal solutions in large position errors.
So, gain errors should be estimated first and then compensated. Here, the method in [18] is used to estimate gain errors
^
gm
1 P
Xfd
1
XIi1cos 4vs ifr xlo xko xl xk
24
Based on (24), Z is a function of the position error X
Z f X
: 25
And Q is also parameterized by the position error X
Q f X
: 26
In [12], Equation (10) is used to estimate gain errors, which means
^g abs
Q 1 w
s
RSS fd
mm 2n
RSS fd
11 2n
.wT Q 1w ; 27
where g = [1, g2,. . .,gM]T, abs () denotes the complex modulus of the elements of a matrix. The gain error estimate, , will be inaccurate because of the existence of the position error X (derived by Q). Moreover, the estimation differences of gain errors will be larger with the increase of position errors.
However, the estimate differences of gain errors can also influence the estimation accuracy of position errors during the second step in the conventional method. It is analyzed as follows. ai can be rewritten as
ai ej 28
where = [1, . . ., M]T and = [1, . . ., M]T denote the amplitude and phase of ai, respectively.
Since (15) is obtained through the Taylor-series expansion, m (m = 2,. . .,M) approaches one but not exactly equals one. Considering the ith virtual calibration sources, the following formula can be obtained
Jic UHai
2
UHai
2
m 1; . . . ; M
;
31
where P is the number of Doppler bins, n2 can be obtained by averaging 2I+2 to M.
After estimating gain errors, phase and position error estimations will be considered in the following. The conventional method requires joint iteration between phase and position error estimations. Through the following study, the joint iteration is avoided in the modified method.
In the certain Doppler bin fd = 0, there exist 2I + 1 spectrum components which are ifr, i = I, . . ., I, due to the existence of Doppler ambiguity. Based on (13), we obtain
D i Di
:
32
Hence, C i = (Ci) is obtained and so as the following equation
Z i C i
H
C i
29
Based on the solution obtained by (10), the element in the kth row and the kth column of ai (which is a diagonal matrix) can be expressed as
Ci
H
Ci
Ci
Ci
Zi
:
H
33
Xu et al. EURASIP Journal on Advances in Signal Processing 2013, 2013:31 Page 6 of 10 http://asp.eurasipjournals.com/content/2013/1/31
Thus, Z can be obtained as
Z
Xii12 real Zi
real Z0
34
Based on (21) and (34), the phase of Q behaves independently of position errors. Hence, using the data of the zero Doppler bin, phase errors can directly be obtained without the influence of existing position errors by the following estimator
^
0 angle
Q 1 w
Xii Ireal Zi
:
.wT Q 1w 35
where ^
0 01; 02; . . . ; 0M
T , angle() returns the phase
angles of a matrix. This result is consistent with the phenomenon mentioned in [12] that when fd = 0 is chosen, the estimation accuracy of phase errors can be satisfactory with position errors assumed to be zero. Since 0k is induced without the effect of position errors, we do not discuss the influence of 0k on the estimation of the position error, X, described by (30).
Based on the analysis above, the modified method can be summarized as follows.
Step 1 Estimate gain errors based on (31) and
compensate the received data using estimated gain errors.
Step 2 Choose the zero Doppler bin to obtain phase
error estimate ^
0 based on (35).
Step 3 Use
^
XMm2RMSE 0mf g ,
1 M 1
0 to reconstruct diag ej^ n o
and
XMm2RMSE xmf g, respectively. L and L0 are set to be 10 and 3, respectively. The Doppler frequency fd is fixed to be zero.
estimate position errors based on (20).
Through the analysis above, phase error estimation method can perform robustly with the existence of position errors. Since the Taylor-series expansion is used in estimating position errors, Step 3 should be applied repeatedly in order to obtain higher estimation accuracy with xmo compensated by xm (which is estimated in the former iteration) in each iteration.
In the modified method, gain errors are compensated first, hence the relationship between gain and position error estimations is eliminated, which guarantees the estimation accuracy.
The change on the iterative fashion makes it possible for the modified method to work with less computational load. Let L denote the iteration number for the conventional method, the computational load is 3LM3 which mainly arises from the matrix eigen-decomposition (equals M3) and the matrix inversion (equals M3). However, the computational load of the modified method is 2 M3 + L0M3 in which L0 denotes the iteration number of Step 3. From the computer simulations, we find that the estimates of the gain, phase, and position errors through the conventional
method do not converge to the final solutions even when L equals to be 10 and the zero Doppler bin is chosen. And for the modified method, L0 should be more than two to satisfy the estimation accuracy. For example, based on the assumptions that L = 10 and L0 = 3, the
computational load of the modified algorithm is 5 M3,
which is lower than that of the conventional method(30 M3).
3. Simulation experiment
To verify the effectiveness of the modified method, computer simulations are given in this section. Experimental results are obtained based on 200 trials. The simulated DSS-SAR system in the side-looking mode is composed by seven satellites and all the satellites distribute uniformly with the antenna spacing l = vsTr/M (Tr is the
pulse repetition period) in the along-track direction [9,12]. The other parameters are listed in Table 1. The gain, phase, and position errors for each of the M satellites, i.e., gm, 0m, and xm, are defined as uniform random variables in [1 , 1 +], [, ], and [p, p], respectively. They are generated and maintained invariant in each trial. The average root mean squared error (ARMSE) of gain, phase, and position errors are
defined as 1
M 1
XMm2RMSE gmf g ,1 M 1
3.1. Effect of SNR
Suppose p = 0.25 l and = 0.2. The ARMSE curves of gain, phase, and position error estimates versus signal-to -noise ratio (SNR) are shown in Figure 2ac, respectively. From Figure 2, it is observed that both the modified method and the conventional method perform better as the SNR increases. From Figure 2a,c, it can be seen that gain and position error estimation performances of the conventional method are worse than those of the modified method. It is because gain
Table 1 Simulation parameters of DSS-SAR systems Orbit altitude 750 km
Incidence angle 45
Band X-band
Antenna size 2 m 1 m
Bandwidth 100 MHz
PRF 1496 Hz
Velocity 7481.5 m/s
Xu et al. EURASIP Journal on Advances in Signal Processing 2013, 2013:31 Page 7 of 10 http://asp.eurasipjournals.com/content/2013/1/31
(a)
0.35
the conventional method
the modified method
0.3
ARMSE of the gain error estimates
0.2
and position error estimations cannot be separated in the conventional method and they interfere with each other. However, in the modified method, because gain errors are corrected before estimating position errors, position error estimates can converge to optimal solutions. The position error estimate ARMSE is lower than 0.01 m when the SNR is larger than 15 dB. In addition, it is shown in Figure 2b that phase error estimations can work well for both methods. This is consistent with the analysis in modified method for DSS-SAR system, which shows that position errors have no influence on the phase error estimation when the zero Doppler bin is chosen. However, the computational load of the modified algorithm is lower than the conventional method.
0.25
0.15
0.1
0.05
0 0 5 10 15 20 25 30
SNR(dB)
(b)
7
(a)
ARMSE of the phase error estimates (degree)
the conventional method
the modified method
0.57
ARMSE of the phase error estimates (degree)
ARMSE of the position error estimates (degree)
6
0.565
before gain error calibrated
after gain error calibrated
5
0.56
4
0.555
3
0.55
2
0.545
1
0.54
0.535
0 0 5 10 15 20 25 30
SNR(dB)
0.53 0 0.1 0.2 0.3 0.4
(c)
0.14
(b)
the conventional method
the modified method
AESME of the position error estimates (m)
0.12
0.35
0.3
0.1
0.25
before gain error calibrated
after gain error calibrated
0.08
0.2
0.06
0.15
0.04
0.1
0.02
0.05
0 0 5 10 15 20 25 30
SNR(dB)
0 0 0.1 0.2 0.3 0.4
Figure 2 ARMSE of gain, phase and position error estimates versus SNR. (a) ARMSE of gain error estimates versus SNR; (b) ARMSE of phase error estimates versus SNR; (c) ARMSE of position error estimates versus SNR.
Figure 3 ARMSE of phase and position error estimates versus gain error by the modified method. (a) ARMSE of phase error estimates versus gain error; (b) ARMSE of position error estimates versus gain error.
Xu et al. EURASIP Journal on Advances in Signal Processing 2013, 2013:31 Page 8 of 10 http://asp.eurasipjournals.com/content/2013/1/31
(a)
(a)
8
0.31
7
the conventional method
the modified method
0.3
the conventional method
ARMSE of the gain error estimates
ARMSE of the gain error estimaties
6
0.4
0.29
5
0.2
0.28
4
0 0 0.07 0.14 0.21 0.28
3
0.27
2
0.26
1
0.25
0 0 0.07 0.14 0.21 0.28 0.35 0.42 0.49
2 4 6 8 10
number of iterations
p ( antenna spacing l units)
(b)
0.547
(b)
60
ARMSE of the phase error estimates (degree)
0.5465
the conventional method
the modified method
40 0.55
ARMSE of the phase erroe of estimates (degree)
the conventional method
50
0.6
0.546
0.5455
30 0.5 0 0.07 0.14 0.21 0.28
20
0.545
10
0.5445
2 4 6 8 10
number of iterations
0 0 0.07 0.14 0.21 0.28 0.35 0.42 0.49
(c)
p ( antenna spacing l units)
0.12
(c)
ARMSE of the position error estimates (m)
0.1
12
the conventional method
the modified method
0.08
ARMSE of the position error estimates (m)
10
the conventional method
the modified method
0.06
8
0.2
0.1
0.04
6
0 0 0.07 0.14 0.21 0.28
0.02
4
0
2 4 6 8 10
2
number of iterations
0 0 0.07 0.14 0.21 0.28 0.35 0.42 0.49
Figure 5 ARMSE of gain, phase and position error estimates versus number of iterations. (a) ARMSE of gain error estimates versus number of iterations; (b) ARMSE of phase error estimates versus number of iterations; (c) ARMSE of position error estimates versus number of iterations.
p ( antenna spacing l units)
Figure 4 ARMSE of gain, phase and position error estimates versue position error. (a) ARMSE of gain error estimates versus position error; (b) ARMSE of phase error estimates versus position error; (c) ARMSE of position error estimates versus position error.
Xu et al. EURASIP Journal on Advances in Signal Processing 2013, 2013:31 Page 9 of 10 http://asp.eurasipjournals.com/content/2013/1/31
3.2. Effect of gain errors
In the simulation, changes from 0 to 0.4. p = 0.25 l. SNR is 20 dB. The other parameters in the simulation are the same as the former. The effect of gain estimation errors on the following phase and position error estimations in the modified method is given in Figure 3, where the curve gained without calibrating gain errors is denoted as before gain error calibrated while the other curve named after gain error calibrated is obtained by calibrating gain errors first. It can be seen that if gain errors cannot be well calibrated, the position error estimation difference becomes larger with the increment of , while the change of phase error estimation can be ignored. By using the modified method, both of phase and position error estimations can work well. It is clear that gain errors have serious influence on the following position error estimation and it is necessary to estimate gain errors first.
3.3. Effect of position errors
Fixing SNR to 20 dB. = 0.2. ARMSE curves of gain, phase, and position error estimates versus position error are shown in Figure 4ac, respectively. Local simulation results (small figures embedded in Figure 4) are also given for better explanation. Being consistent with the analysis in modified method for DSS-SAR system, gain and position error estimations interfere with each other, which make gain and position error estimations converge to suboptimal solutions in the conventional method. So, the estimation accuracy of the conventional method cannot satisfy the requirement. The gain error estimate ARMSE is larger than 0.1 while the position error estimate ARMSE is close to 0.03 when p = 0.07 l. In addition, phase error estimates obtained by both of the conventional method and the modified method are approximately equal to each other when p 0.28 l. This is consistent with the previous analysis. By using the data of the zero Doppler bin, phase errors can be estimated with existing position errors. However, when p 0.35 l, the conventional method fails to work since gain and position error estimates deviate badly from the true values which has a strong impact on the phase error estimation.
3.4. Effect of iterations
Fix = 0.2, p = 0.25 l and suppose SNR to be 20 dB. The gain, phase, and position error estimation performance of the conventional method versus the number of iterations is shown in Figure 5ac, respectively. It is observed from Figure 5a,c that the conventional method converges slowly and even to suboptimal solutions since gain and position error estimations can interact with each other. However, it is noticed from Figure 5b that the phase error estimation can work well as we analyzed before.
For the modified method, only the position error estimate ARMSE versus the number of iterations is shown in Figure 5c since only the position error estimation is iterated in the modified method. It is noticed from Figure 5c that the position error estimation converges to the optimal solution with less than three iterations. To sum up, the computational load of the modified method is less than that of the conventional method. Meanwhile better performance can be obtained by the modified method.
4. Conclusions
In this study, we focus on gain, phase, and position error estimations of DSS-SAR systems. Based on the conventional method, a modified array error estimation method is proposed here. In the conventional method, the estimation of gain and position errors may converge to suboptimal solutions especially when position errors are large and the joint iteration between the gain-phase and position error estimations are needed. That is because of the interaction between the estimations of gain and position errors. In the modified method, gain errors are first estimated and compensated before the other errors estimation, which guarantees the position error estimation accuracy to be higher than that in the conventional method. Meanwhile, by using the zero Doppler bin data, the phase error estimation can behave independently of position errors. Then the joint iteration strategy between the gain-phase and position error estimations in the conventional method can be avoided, which makes the modified method perform stably with low computational load. Theoretical analysis and simulation results demonstrate the effectiveness of the modified method.
Competing interestsThe authors declare that they have no competing interests.
AcknowledgmentThe authors would like to thank the anonymous reviewers for their constructive comments, which led to significant improvements in this study. Special thanks to Dr. Yinghua Wang and Xuefeng Zheng for improving the English usage and also for their valuable suggestions. This study was sponsored in part by the National Basic Research Program of China (973 Program) under Grant no. 2011CB707001, in part by the Natural Science Foundation of China under Grant no. 61101242, in part by the Fundamental Research Funds for the Central Universities under Grant no. K50511020009, and in part by the Program for Changjiang Scholars and Innovative Research Team in University under Grant no. IRT0954.
Received: 19 June 2012 Accepted: 7 December 2012 Published: 22 February 2013
References1. D Zhen, C Bin, L Diannong, Detection of ground moving targets for two-channel spaceborne SAR-ATI. EURASIP J. Adv. Signal Process., Article ID 230785 (2010). doi:http://dx.doi.org/10.1155/2010/230785
Web End =10.1155/2010/230785
2. F Lombardini, M Pardini, G Fornaro et al, Linear and adaptive spaceborne three-dimensional SAR tomography: a comparison on real data. IET Radar Sonar Navigat. 3(4), 424436 (2009)
3. P Wolfgang, M David, The Terra SAR-X satellite. IEEE Trans. Geosci. Remote Sens. 48(2), 615622 (2010)
Xu et al. EURASIP Journal on Advances in Signal Processing 2013, 2013:31 Page 10 of 10 http://asp.eurasipjournals.com/content/2013/1/31
4. S Chiu, MV Dragosevic, Moving target indication via RADARSAT-2 multichannel synthetic aperture radar processing. EURASIP. J. Adv. Signal Process., (2010). doi:http://dx.doi.org/10.1155/2010/740130
Web End =10.1155/2010/740130
5. K Tomiyasu, Image processing of synthetic aperture radar range ambiguous signals. IEEE Trans. Geosci. Remote Sens. 32(5), 11141117 (1994)
6. NA Goodman, SC Lin, D Rajakrishna, JM Stiles, Processing of multiple receiver, spaceborne arrays for wide-area SAR. IEEE Trans. Geosci. Remote Sens. 40(4), 841852 (2002)
7. S Li, HP Xu, LQ Zhang, An advanced DSS-SAR InSAR terrain height estimation approach based on baseline decoupling. Prog. Electromagn. Res. 119, 207224 (2011)
8. CH Gierull, D Cerutti-Maori, J. Ender, Ground moving target indication with tandem satellite constellations. IEEE Geosci. Remote. Sens. Lett. 5(4), 710714 (2009)
9. JP Aguttes, The SAR train concept: required antenna area distributed over N smaller satellites, increase of performance by N, in Proceedings of the IGARSS03, Toulouse, France (, 2003), pp. 542544
10. ZF Li, HY Wang, T Su, Z Bao, Generation of wide-swath and high-resolution SAR images from multichannel small spaceborne SAR Systems. IEEE Geosci. Remote. Sens. Lett. 2(1), 8286 (2005)
11. Z Lei, Q Cheng-Wei, X Mengdao, B Zheng, SAR imaging and Doppler ambiguity removal with distributed microsatellite arrays. Int. J. Remote. Sens. 31(24), 64416458 (2010)
12. ZF Li, Z Bao, HY Wang, GS Liao, Performance improvement for constellation SAR using signal processing techniques. IEEE Trans. Aerosp. Electron. Syst. 42(2), 436451 (2006)
13. AJ Weiss, B Friedlander, DOA and steering vector estimation using a partially calibrated array. IEEE Trans. Aerosp. Electron. Syst. 32(3), 10471057 (1996)
14. KV Stavropoulos, A Manikas, Array calibration in the presence of unknown sensor characteristics and mutual coupling, in Proceedings of the EUSIPCO, Toulouse, France (, 2000), pp. 14171420
15. L Aifei, L Guisheng, Z Cao et al, An eigenstructure method for estimating DOA and sensor gain-phase errors. IEEE Trans. Signal Process.59(12), 59445956 (2011)
16. B Friedlander, AJ Weiss, Direction finding in the presence of mutual coupling. IEEE Trans. Antennas Propagat. 39(3), 273284 (1991)
17. AJ Weiss, B Friedlander, Array shape calibration using sources in unknown locations-a maximum likelihood approach. IEEE Trans. Acoust. Speech Signal Process. 37(12), 19581966 (1989)
18. MP Wylie, S Roy, RF Schmitt, Self-calibration of linear equi-spaced (LES) arrays, in Proceedings of the ICASSP-93, Minneapolis, USA (, 1993), pp. 281284
doi:10.1186/1687-6180-2013-31Cite this article as: Xu et al.: Self-calibration method without joint iteration for distributed small satellite SAR systems. EURASIP Journal on Advances in Signal Processing 2013 2013:31.
Submit your manuscript to a journal and benet from:
7 Convenient online submission7 Rigorous peer review7 Immediate publication on acceptance7 Open access: articles freely available online 7 High visibility within the eld7 Retaining the copyright to your article
Submit your next manuscript at 7 springeropen.com
The Author(s) 2013