2 May 2016 Array error calibration methods in downward-looking linear-array three-dimensional synthetic aperture radar
Author Affiliations +
J. of Applied Remote Sensing, 10(2), 025010 (2016). doi:10.1117/1.JRS.10.025010
Abstract
In order to achieve high-precision three-dimensional (3-D) imaging with an airborne downward-looking linear-array 3-D synthetic aperture radar (LA-3D-SAR), a uniform virtual antenna array can be obtained by aperture synthesis of the cross-track sparse multiple-input-multiple-output array. However, the actual 3-D imaging quality is unavoidably degraded by array errors such as the multichannel amplitude-phase errors due to the nonideal antenna characteristics, and the virtual element position errors due to vibrations and motion measurement deviations. We investigate the effects of these errors on the forms and the degrees of image quality degradation and consider the use of corresponding calibration methods to eliminate the effects of errors. For the multichannel amplitude-phase errors, the target response is subject to an integrated sidelobe level increase introduced by the phase error, which can be calibrated based on external (parallel or point target) calibrators, as proposed in the paper. For the virtual element position errors, they mainly the result of contrast degradation and noise in the image along the cross-track direction and have little impact on the range and along-track directions. The imaging performance is more sensitive to the error component in the height direction as compared to other components, the precision requirement of which should be established as the calibration reference. A calibration method based on time-divided active calibrators is proposed to estimate and correct the virtual element position errors. Both numerical simulations and real data experiments have shown the validity of the analyses as well as the effectiveness of the proposed calibration methods.
Tan, Huang, Han, Liu, and Peng: Array error calibration methods in downward-looking linear-array three-dimensional synthetic aperture radar

1.

Introduction

In the last decade, an increased interest has been devoted to three-dimensional synthetic aperture radar (3D-SAR),1 because of its capability to generate 3-D radar images by retrieving the reflectivity and real 3-D localization information of the scattering objects. Among them, the downward-looking linear array SAR (LA-3D-SAR)2 with a linear antenna array [e.g., a sparse multiple-input-multiple-output (MIMO) array34.5.6.7.8 is usually adopted to reduce the complexity of the linear array antenna, where multiple antennas are used for transmitting while others are for receiving.] deployed under the wings of airplanes as the most feasible configuration for topographical survey application. With the rectilinear motion and the downward-looking working mode, LA-3D-SAR synthesizes a 2-D plane array to achieve 3-D imaging of illuminated areas and overcomes the layover and shadowing problems existed in 2-D SAR. Until now, two LA-3D-SAR systems have been developed,9,10 and some 3-D image formation algorithms have been proposed.1112.13.14.15.16

As a multichannel coherent system, LA-3D-SAR is particularly susceptible to channel inconsistency, which will cause geometric distortion, loss of resolution, decrease in image contrast, spurious targets, and a reduction in signal-to-noise ratio (SNR).17 However, under realistic operational conditions, the antenna elements are subject to array errors such as multichannel amplitude-phase errors and antenna position errors. These errors come from the nonideal antenna characteristics, vibrations, and motion measurement errors, which may cause channel inconsistency that is large enough to deteriorate the 3-D imaging quality, especially for the LA-3D-SAR working in a high frequency band (e.g., X, Ka band). Thus, the study on the high-precision amplitude-phase consistency of LA-3D-SAR is anticipated across all MIMO channels.

In Refs. 18 and 19, the error effects and calibration methods for channel inconsistency based on the deterministic patterns of the wings’ oscillation and the flight attitude variations are investigated. In this paper, we focus on two other kinds of array errors: multichannel amplitude-phase errors and random virtual element position errors, neither of which has received sufficient attention in the current research. According to their mechanisms, different error models are established, and their effects on the LA-3D-SAR images are studied. A close examination of the effects of the errors with regard to their characteristics, the levels of severity, and the relevant calibration methods is presented to specify the tolerable levels, which is important for the system design and optimization. Specific calibration methods are proposed to deal with the array errors, which are verified by either the ground-based experiment data or some numerical simulations.

The remaining sections are organized as follows. The imaging principle of LA-3D-SAR is first reviewed in Sec. 2, where the geometric and signal models are given. In Sec. 3, with different error models, the error effects of the multichannel amplitude-phase errors and the virtual element position errors are analyzed. In Sec. 4, different calibration methods based on external calibrators are proposed to deal with the errors, respectively. The ground-based experiment and numerical simulation verification are presented in Sec. 5. This paper is concluded in Sec. 6.

2.

Imaging Principle of Downward-Looking LA-3D-SAR

As shown in Fig. 1, the radar is mounted on an airplane flying at altitude H with the velocity Va. Let the X-, Y- and Z-axes denote the cross-track, along-track, and height directions, respectively. In cross-track, a sparse MIMO array with MT transmitters and MR receivers is adopted to form a uniform virtual array, which has M=MTMR virtual elements and total length L=(M1)d, where d represents the element spacing. The position of the m’th virtual element is (xm,yn,H), where yn represents the along-track spatial sample coordinate. The illuminated scenario is located at the nadir area of the platform, and for a point target P with reflectivity σp located at (xp,yp,zp) in the 3-D target support band Ω, the instantaneous target-to-m’th virtual element distance is

(1)

R(xm,yn;P)=(xmxp)2+(ynyp)2+(Hzp)2.
Assuming that orthogonal pulses with pulse length Tp, frequency bandwidth B, and carrier frequency fc are transmitted, the echo signal data corresponding to each virtual element channel will transform to the following form after carrier frequency down-conversion, orthogonal demodulation,6 and range compression:

(2)

src(t,xm,yn)=Ωrect[xmL]rect[ynypLsyn]σpej4πfccRsinc[B(t2Rc)]dxpdypdzp,
where t denotes the fast-time variable, Lsyn denotes the along-track synthetic aperture length, and c is the speed of light. For notational simplicity, we use rectangular window to represent the antenna patterns and signal weighing functions in Eq. (2). In some continuous wave radar, the wavenumber domain form of the range compressed signal src(t,xm,yn) is directly accessible, which is

(3)

Src(k,xm,yn)=Ωrect[kkcB]rect[xmL]rect[ynypLsyn]×σpexp(jk2R)dxpdypdzp,
where k[kmin,kmax] is the fast-time wavenumber variable, with kmin=2π(fcB/2)/c, kmax=2π(fc+B/2)/c, and kc=2πfc/c.

Fig. 1

Imaging geometry of downward-looking LA-3D-SAR.

JARS_10_2_025010_f001.png

Usually, the process of imaging is to estimate the distribution of the target reflectivity σp from the range compressed signal src(t,xm,yn) or its wavenumber domain signal Src(k,xm,yn). For example, in Ref. 9, after range compression, the compression in along-track is done using the synthetic aperture principle and the compression in cross-track is done using a digital beamforming. After that, true 3-D resolving ability is achieved so that the LA-3D-SAR can overcome the layover and shadowing effects in the conventional 2-D SARs. However, this is not the case when array errors, such as the multichannel amplitude-phase errors or virtual element position errors are present, where the signal coherency in the LA-3D-SAR is deteriorated and the final image may degrade drastically. Thus, it is essential to carry on a detailed investigation of the error effects and the related calibration methods.

3.

Array Errors and the Effects on LA-3D-SAR Imaging

3.1.

Multichannel Amplitude-Phase Errors

The multiple transmit and receive channels consist of active components (e.g., high frequency power amplifiers, mixers, intermediate frequency amplifiers and so on), whose amplitude and phase characteristics are slowly varying with the time. In addition, due to the nonideal electrical characteristic of antennas in the array, each channel may have a random deviation in the gain (or amplitude) and phase-shift from the ideal situation, which results in the multichannel amplitude-phase errors. This type of error is independent of the along-track direction and mainly affects the imaging along the array (i.e., the cross-track direction), so that it is appropriate to simplify the analysis using only an equivalent 2-D geometric and signal model.

As shown in Fig. 2, X and Y represent the cross-track and slant range directions, respectively. Assume each virtual element’s position is (xm,0), where m=0,1,2M1. Using R=yp2+(xmxp)2 to replace the instantaneous range in Eq. (1) and taking the multichannel amplitude-phase error into consideration, the signal corresponding to the m’th channel can be expressed as

(4)

sm(k)=AmΩσP×exp{jk·2yp2+(xmxp)2+jΔΦm}dΩ=AmΩσP×exp{jk·2yp2+(xmxp)2+j(k·Δrm+ϕm,c)}dΩ,
where Am denotes the gain variations due to the amplitude error and ΔΦm=k·Δrm+ϕm,c denotes the phase error, which can be divided into two parts: the first term k·Δrm representing the different time delay (or shift) components for each channel is referred to as range deviation phase error, and the other part ϕm,c is referred to as residue phase error.

Fig. 2

Equivalent 2-D imaging geometry.

JARS_10_2_025010_f002.png

For almost all the imaging algorithms, the cross-track focusing is done by coherent accumulation along the array direction, i.e.,

(5)

s=m=0M1AmejΦejΔΦm=ejΦm=0M1AmejΔΦm,
where Φ is the coherent phase that the algorithm compensates to. From Eq. (5), it is apparent that the target response does not achieve coherent accumulation because of the presence of multichannel amplitude-phase error Am and ΔΦm. According to the geometric and signal model, each antenna element in the array almost illuminates the same target area, so that it is reasonable to eliminate the target response deterioration introduced by the amplitude error Am by a simple power normalization of the channel data.

Compared to the amplitude error, the target response is more sensitive to the phase error. The range deviation phase error k·Δrm affects the imaging in a similar way as the range cell migration effects in conventional SAR imaging. We can use the integrated sidelobe ratio (ISLR) along the array as the figure of merit to evaluate the severity of the degradation induced by this error. Assume the channel number is M=101, center frequency fc=31 GHz, and the random range deviation corresponding to each channel conforms to the uniform distribution, i.e., ΔrmU[0,C], m=0,1,2M1, where C is a constant. Figure 3(a) shows the curve where the cross-track ISLR changes with parameter C. It is apparent that ISLR=10.43  dB for the error-free case, which accords with the ideal SAR impulse response, and the ISLR gets worse with the increase of range deviation error, e.g., when the deviation is beyond 0.3 wavelength, the target energy has spread out across the array totally, so that the ISLR measurement becomes meaningless.

Fig. 3

Variation of cross-track ISLR with (a) the range deviation error and (b) the residue phase error.

JARS_10_2_025010_f003.png

The residue phase error ϕm,c in Eq. (4) manifests itself as a constant (within one channel) phase error attached to each different array channel. Let us denote the signal along the array direction as s(xm), and the ideal target response as psfx(x), then according to the Fourier transform equation,20 the practical signal s(xm)·exp{jϕm,c} can be processed after coherent accumulation, which yields

(6)

s(x)=psfx(x)*IDFT{exp(jϕm,c)},
where * denotes the discrete convolution. Equation (6) implies that the ideal target impulse response spreads over the X axis because of the presence of IDFT{exp(jϕm,c)}. Again, using the ISLR evaluation methodology, except that the residue phase error ϕm,c is set to conform to a uniform distribution ϕm,cU[0,C·π], m=0,1,2M1. As shown in Fig. 3(b), the target response is deteriorated with the increase of the residue phase error. For example, when ϕm,c distributes randomly between 0 and π, the level of ISLR has been creeping up from 3 to 2  dB, where correct imaging is impossible to achieve.

3.2.

Virtual Element Position Errors

The virtual element position error is usually caused by wind deformation, gravity deformation, platform vibration, installation or motion measurement errors when a conformally designed antenna array is utilized in the airplane platform. It is usually expressed as the virtual element position deviation δm=Δxm·ix+Δym·iy+Δzm·iz from the ideal situation (xm,yn,H), where ix, iy, and iz represent the unit vector in cross-track, along-track, and height direction, respectively. An intuitive visualization of the virtual element position errors is provided in Fig. 4, where the actual position of the m’th virtual element is represented as (xm+Δxm,yn+Δym,H+Δzm), so the actual target-to-sensor distance is

(7)

R=(xm+Δxmxp)2+(yn+Δymyp)2+(H+Δzmzp)2.

Fig. 4

Virtual element position errors.

JARS_10_2_025010_f004.png

From Eqs. (2) and (7), and using r=(c/2)t, we can write the range compressed impulse response as

(8)

sr(r,xm,yn)=ejkc2Rpsfr(rR).
Taking the binomial expansion of Eq. (8), which yields

(9)

RRm0+(ynyp)22Rm0+Δxm2+Δym2+Δzm22Rm0+(xmxp)ΔxmRm0+(ynyp)ΔymRm0+(Hzp)ΔzmRm0,
where Rm0=(xmxp)2+(Hzp)2 is the minimum range between the point target to the m’th virtual element. The third term in Eq. (9) is so small that can be neglected, then inserting the remaining terms into Eq. (8), and neglecting the range curvature migration term in psfr(·), results in

(10)

sr(r,xm,yn)=ejkc2Rm0ejkc(ynyp)2Rm0ej(Φx+Φy+Φz)psfr(rRm0).

The first exponential term in Eq. (10) contains information for cross-track beamforming, the second exponential term is the modulation term for along-track compression, and Φx, Φy and Φz represent phase error terms resulting from the virtual element position errors Δxm, Δym, and Δzm, respectively, which are defined as

(11)

{Φx=2kc(xmxp)Δxm/Rm0,Φy=2kc(ynyp)Δym/Rm0,Φz=2kc(Hzp)Δzm/Rm0.

As detailed in Ref. 5, for the range compressed data, a 3-D SAR image can be obtained by performing the along-track focusing and cross-track beamforming. However, the 3-D imaging quality is practically degraded by the phase error terms Φx, Φy, and Φz presented in Eq. (10) with high probability. In this paper, we relate the errors’ effects to the image degradation by deriving appropriate statistical measures. We first model the virtual element position errors Δxm, Δym, and Δzm as discrete, white, Gaussian stochastic process with zero-mean and variance of σx2, σy2, and σz2, respectively. Δxm, Δym, and Δzm are random variables over the cross-track virtual element positions and slowly time-varying along the along-track direction. From Eq. (11), we can see that Φy is a compound of linear phase error in the along-track direction and random phase error in the cross-track direction, whereas Φx and Φz are both random phase errors in the cross-track dimension.

In our previous research,21 we decompose the phase errors into different components and discuss their respective effects on imaging. Usually, a constant phase error component will have little effect on imaging, a linear phase error component will cause target displacement, a second-order component representing quadratic phase error will cause mainlobe broadening and sidelobe increase, and higher-order components manifesting as oscillation functions will cause paired echoes, which increase the integrated sidelobe level. We propose to use the discrete Legendre orthogonal transform22 to guarantee the mutual orthogonality of the expansion terms, and then appropriate statistical measures relating the effects of random phase errors to the degradation of the final image can be derived. From the analyses, it can be concluded that the virtual element position errors will contribute to a high cross-track integrated sidelobe level, which results in noise and contrast degradation in the image, and has little impact along the range and along-track directions.

Figure 5 shows the relationship between the error effects and the orthogonal expansion coefficients as well as the virtual element position error statistics for a typical Ka-band LA-3D-SAR parameters set. Let ak denote the k’th order orthogonal expansion term coefficient of Φx, and σak2 as its variance. Figure 5(d) implies that to meet a typical 17  dB ISLR demand for SAR imaging,23 σak2 has to be kept below 104. Otherwise, the resultant image will suffer from intolerable noise and contrast degradation in the image. Figure 5(e) shows the curves of the coefficient variance varying with σx, which are evaluated with four different point targets. For targets far away from the array, the coefficient variance is large and a 104 coefficient variance corresponds to σx with the value of 0.5 mm, which is taken as the calibration precision requirement of virtual element position error Δxm. We can also derive from Figs. 5(a)5(e) that, for virtual element position errors with millimeter-scale levels, the sidelobe increase due to paired echoes dominates the error effects, while the mainlobe is almost intact due to small low-order phase error components. Figure 5(f) gives the curves of the coefficient variance σbk2 of Φz changing with the standard deviation σz of the virtual element position error component Δzm, again evaluated with four point targets located at different positions. It is apparent that the four curves almost overlapped, and a 104 coefficient variance corresponds to σz value of approximately 0.075 mm. Some calibration or compensation procedures would be necessary if Δzm exceeds this limit.

Fig. 5

Curves showing: (a) the cross-track displacement as a function of the first-order coefficient variance; (b) and (c) cross-track mainlobe broadening and ISLR as functions of the second-order coefficient variance; (d) cross-track ISLR as a function of the high-order coefficient variance; (e) and (f) the dependence of coefficient variance on the virtual element position errors. Ka-band is assumed where the parameters are the same as those given in Table 1 in Ref. 13.

JARS_10_2_025010_f005.png

Compared with Δym, the beamforming performance along the cross-track is more sensitive to Δzm, which is illustrated in Fig. 6 with a Ka-band point analysis example. In Fig. 6(a), though the mainlobe of the cross-track impulse response is intact as the error-free case, the impulse response is corrupted by high sidelobes (an ISLR above 0 dB) resulting from virtual element position errors when σx=5  mm. When σx is reduced to 0.5 mm, as shown in Fig. 6(b), the ISLR is reduced to an acceptable level of 17  dB correspondingly, which satisfies the imaging requirement. However, as shown in Fig. 6(c), the same 0.5 mm level of errors as applied in the Z direction will still spread out energy across the impulse response and give rise to an ISLR of 1.5  dB. In Fig. 6(d), when the σz is set to 0.075 mm, a favorable cross-track impulse response with ISLR of 17  dB is obtained. Thus, in order to mitigate the impact of the rising sidelobes, calibration measures should be taken to keep the error below a tolerable level, and the precision requirement of Δzm should be taken as a reference for calibrating the virtual element position errors.

Fig. 6

The cross-track impulse response when the virtual element position errors Δyn or Δzn exists.

JARS_10_2_025010_f006.png

4.

Calibration of the Array Errors

4.1.

Calibration of the Multichannel Amplitude-Phase Errors

4.1.1.

Amplitude error calibration using power normalization

For LA-3D-SAR imaging, each antenna element in the array almost illuminates the same target area, which implies that the amplitude variations among different channel data are caused by the random amplitude error Am, thus can be compensated by a simple power (sum of the magnitude squared) normalization of the channel data, i.e., assuming the signal corresponding to the m’th channel is sm(n), with the power value of Em and fast-time samples number N, then power normalization can be done through

(12)

s˜m(n)=sm(n)Em,m=0,1,,M1,n=1,2,,N,
where s˜m(n) is the m’th channel data after amplitude error calibration.

4.1.2.

Phase error calibration using external calibrators

The multichannel phase errors can be calibrated on a static platform indoors before the array is applied on the airplane platform. Due to the slowly-varying characteristic of the phase errors, the calibration parameters can be stored onboard for the SAR processor to look up and compensate within a certain period of time.

Figure 7 shows the calibration geometries using two different external calibrators. Figure 7(a) uses an external calibrator, which is placed in parallel to the array direction, such that the distance (range) histories and phase histories are expected to be the same for all channel data. Figure 7(b) uses a point target (e.g., a corner reflector) with strong reflectivity as the external calibrator, whereby the phase error calibration parameters can be estimated by virtue of the triangular geometry relationship of the array and the point target calibrator.

Fig. 7

Calibration geometry: (a) with a parallel external calibrator and (b) with a point target calibrator.

JARS_10_2_025010_f007.png

More specifically, take the calibration procedure in Fig. 7(a) as an example. First, the range compressed data (in 2-D time or spatial domain) are upsampled along the fast-time direction, then the distance histories R^m of the calibrator in all channels are estimated by cross-correlating the data with the range point spread function psfr(t):

(13)

{μm(τ)=sm(t)psfr(tτ)dtτ^m=argmaxτμm(τ)R^m=c2×τ^m.
Then the difference of R^m and the ideal distance value R0, i.e., Δr^m=R^mR0, is taken as the estimated value of the range deviation error corresponding to the m’th channel:

(14)

Sm(k)=Sm(k)×exp{j2kΔr^m}.
After correcting the range deviations, the energy trajectory of the calibrator locates at the range gate of R0 and manifests as a parallel line with respect to the array. Then the phase values ϕ^m,c, m=0,1,,M1 on the calibrator response peak can be extracted as the estimated residue phase errors, and thus can be corrected through

(15)

Sm(k)=Sm(k)×exp{jϕ^mc}.

The calibration method in Fig. 7(b) is similar, except that the ideal energy trajectory exhibits the form of a hyperbolic curve rather than a straight line, and the range deviation phase error values are estimated by virtue of the triangular geometry, whereas the residue phase error values are extracted along the curved energy trajectory. Ground-based experiment examples given in the latter part of the paper have confirmed the feasibility of the calibration method.

4.2.

Calibration of the Virtual Element Position Errors Using Time-Divided Active Calibrators

For the virtual element position errors, one way of calibration is to utilize distributed motion measurement units to retrieve the exact position of each antenna and then correct the error effects by postprocessing.24 However, considering the big amount of elements and the limited precision of the state-of-the-art measurement device, this approach will be too complicated and expensive.

Another feasible approach is to use the array calibration methods in array signal processing theory.2526.27 In the downward-looking LA-3D-SAR imaging model, Δxm, Δym, and Δzm are random variables over the cross-track virtual element positions and slowly time-varying over the along-track integration time, which means that for a set of adjacent along-track spatial samples yn, the position errors are deterministic and constant. For a certain yn, the cross-track beamforming is analogous to the direction of arrival (DOA) estimation problem28 in array signal processing theory.

In a typical DOA radar, the array calibration methods usually use active calibrators (or signal sources) that have a priori position and angle (DOA) information relative to the array. By reversely deriving the DOA estimation algorithm (e.g., the MUSIC algorithm29), the actual antenna phase centers of the array elements can be retrieved. For a fixed array, the steering vector is determined by the element positions and the calibrator’s DOA, so that the element position errors are also involved in the steering vector, which can be retrieved by some optimization methods. The virtual element position errors in LA-3D-SAR can be calibrated in a similar manner, except that the data channels are composed of virtual elements formed by aperture synthesis of a sparse MIMO array, rather than of real receiving elements. A practical solution is to use active calibrators that have relay functions, such as active transponders.30 Such a device is analogous to a corner reflector, in which the calibrator’s receiving antenna receives the SAR signal and then feeds it back to the radar sensor after gain amplification. The transponders adopt transparent relay technology such that the signal phase is preserved and the signal characteristics are not deteriorated. The calibrators can be arranged in the illuminated scene, where the angle information (DOA) with respect to the array in the airborne platform can be accurately measured by ground and airborne position measurement devices. Using this information, the virtual element position errors can be calibrated by estimating the steering vectors of the array.

The following is a detailed description of our proposed calibration method, which is based on time-divided active calibrators. Here, “time-divided” means that the active calibrators are activated time-dividedly one after another, at different along-track times during the calibration process. As shown in Fig. 8, set the central virtual element of the cross-track array as the reference and assuming the along-track time η=0, so that the reference element is located at the coordinate origin (0, 0, 0) and is assumed to be error-free. Assume a calibrator (transponder) Pi is located at (xi,yi,zi) in Fig. 8, where αi and βi represent the along-track and cross-track looking angle with respect to the cross-track array. At time η, let yn denote the along-track position of the LA-3D-SAR, the MT transmitting antennas emit orthogonal waveforms p0(t),p1(t),,pMT1(t), the calibrator Pi relays the signals and later, the signals are received by all MR receiving antennas simultaneously. By aperture synthesis, a total number of M virtual elements are formed and the channel data corresponding to the m’th virtual element are

(16)

sm(r)=ejkc2Rpsfr(rR)+Nm(r),
where r is the fast-time range variable, kc is the center wave number, Nm(r) represents additive channel noise which is a Gaussian, white stochastic process and independent of the signal components, R is the instantaneous range from the m’th virtual element to the calibrator, with the expression:

(17)

R=(xm+Δxmxi)2+(yn+Δymyi)2+(H+Δzmzi)2.

Fig. 8

Calibration of virtual element position error using time-divided active calibrators.

JARS_10_2_025010_f008.png

Let Ri0=(ynyi)2+xi2+(Hzi)2 denote the range from the calibrator to the reference element, then taking the Taylor series expansion of R,31 and neglecting the small quadratic and higher-order terms, yields

(18)

RRi022(ynyi)Δym2(xm+Δxm)xi+2(Hzi)ΔzmRi0(ynyi)(ynyi)2+xi2(ynyi)2+xi2Ri0Δymxi(ynyi)2+xi2(ynyi)2+xi2Ri0(xm+Δxm)+(Hzi)Ri0Δzm=Ri0Δymsin(αi)sin(βi)(xm+Δxm)cos(αi)sin(βi)+Δzmcos(βi).

Inserting Eq. (18) into Eq. (16), and transforming the result into vector-matrix form, which yields

(19)

r=a˜(α,β)s(t)+N,
where

(20)

r=[s0(r)s1(r)sM1(r)]T,

(21)

a(αi,βi)=[exp{j2kcx0cos(αi)sin(βi)},,exp{j2kcxM1cos(αi)sin(βi)}]T,

(22)

a˜(αi,βi)=Δia(αi,βi),

(23)

Δi=diag{[exp{jφ0(αi,βi)}exp{jφM1(αi,βi)}]T},

(24)

φm(αi,βi)=2kc[Δymsin(αi)sin(βi)+Δxmcos(αi)sin(βi)Δzmcos(βi)],

(25)

s(r)=ejkc2R0ipsfr(rR),

(26)

N=[N0(t)N1(t)NM1(t)]T,
where [·]T represents vector transpose, diag{v} represents a diagonal matrix whose diagonal is vector v, Eq. (21) is the error-free steering vector of the array, whereas Eq. (22) is the steering vector when virtual element position errors exist.

The covariance matrix of the array signal r is

(27)

Rrr=E[rrH]=a˜E[SSH]a˜H+σ2I=a˜RSa˜H+σ2I,
where [·]H represents Hermitian transpose. Because the signal and noise are independent, the covariance matrix can be decomposed into two uncorrelated parts corresponding to the signal and noise. Let us take the eigenvalue decomposition of matrix Rrr:

(28)

Rrr=USSUSH+UNNUNH.

In Eq. (28), US represents the signal subspace spanned by the eigenvectors corresponding to big eigenvalues, whereas UN represents the noise subspace spanned by the eigenvectors corresponding to small eigenvalues. Assuming that there is only one active calibrator in the scene, according to the MUSIC algorithm, the eigenvector of the covariance matrix of the array output corresponding to the largest eigenvalue is proportional to the steering vector. Let q denote this eigenvector, then

(29)

a˜^(αi,βi)=c0q,
where c0 is a real constant. a˜^(αi,βi) is the estimated steering vector that contains virtual element position errors. From Eq. (22)–(24), we can estimate the diagonal matrix Δi by solving the following optimization problem:

(30)

Δ^i=argminΔi[a˜^(αi,βi)Δia(αi,βi)]H[a˜^(αi,βi)Δia(αi,βi)]s.t.the matrixΔiwithMrows andMcolumns is diagonal.

Let bm and am denote the m’th elements of vectors a˜^(αi,βi) and a(αi,βi), respectively, i.e., a˜^(αi,βi)=[b1,bm,bM]T and a(αi,βi)=[a1,am,aM]T, and considering that Δi is a diagonal square matrix, i.e., Δi=diag(x1,xm,xM), then it is easy to get

(31)

[a˜^(αi,βi)Δia(αi,βi)]H[a˜^(αi,βi)Δia(αi,βi)]=|b1x1a1|2++|bmxmam|2++|bMxMaM|20.

Because the above expression is non-negative, the minimization is achieved when each additional term of |bmxmam|2 in Eq. (31) becomes 0, which means that the estimated value of xm is bm/am, so that

(32)

Δ^i=diag(b1a1,,bmam,,bMaM)
for which am, m=1,2,M all have nonzero values according to the definition of a(αi,βi) in Eq. (21).

Suppose there are a total number of K time-divided active calibrators in the scene. For each calibrator, we can obtain an estimated diagonal matrix Δ^i, let gi denote the vector formed by the diagonal elements of matrix Δ^i, then the virtual antenna element position errors can be estimated by solving the matrix equation:

(33)

g=Tδ,
where

(34)

g=[Arg{g1}Arg{g2}Arg{gK}]K×MT,

(35)

δ=[Δx0Δx1ΔxM1Δy0Δy1ΔyM1Δz0Δz1ΔzM1]3×M,

(36)

T=2kc[sin(α1)sin(β1)cos(α1)sin(β1)cos(β1)sin(α2)sin(β2)cos(α2)sin(β2)cos(β2)sin(αK)sin(βK)cos(αK)sin(βK)cos(βK)]K×3,
where Arg{·} is the operator for obtaining the angle of a complex quantity. Because δ contains 3(M1) unknown variables and each calibrator can provide (M1) equations to the linear system, then there should exist at least three calibrators to solve the matrix equation. If K=3, the unknown variables can be determined uniquely by multiplying the inverse of a nonsingular matrix:

(37)

δ=T1g.

In such a case, the DOA angular interval among the calibrators should not be too close in order to avoid the matrix singularity. If more than three calibrators are provided, i.e., K>3, we say the linear system is overdetermined and thus can be solved by the least square method.32 However, if K<3, then the linear system is underdetermined, which means the unknowns cannot be determined uniquely. So, to complete the calibration mission in practice, more than three calibrators with a dispersed distribution in the scene are preferable.

5.

Experimental Verification Results

5.1.

Indoor Experiments of the Multichannel Amplitude-Phase Error Calibration

5.1.1.

Two-dimensional imaging experiment

To verify the multichannel amplitude-phase error calibration methods, we utilize a Ka-band horn antenna array horizontally placed on a static indoor platform, which is shown in Fig. 9. An Agilent E8363C vector network analyzer (VNA) is used as the signal transceiver.33 With the help of a microwave multiplexer to switch the transmit and receive channel from one antenna to another on a pulse to pulse basis, a linear array with 60 virtual elements and an 8-mm inner-element spacing is synthesized. The test target scenario is composed of three metallic spheres, and before acquiring the backscattered data from them, a metal sheet functioning as a calibrator is placed parallel to the array. The experiment parameters are listed in Table 1.

Fig. 9

2-D imaging experiment setup for calibrating the multichannel amplitude-phase errors.

JARS_10_2_025010_f009.png

Table 1

2-D imaging experiment parameters.

ParametersValue
Center frequency31 GHz
Frequency bandwidth6 GHz
Frequency points1601
Array aperture length0.472 m
Inner-element spacing8 mm
Virtual element number60

The calibration procedure is implemented following Eqs. (12)–(15), then the calibration parameters are applied, and the results are shown in Fig. 10. Figure 10(b) shows the degraded imaging result without the multichannel amplitude-phase error correction. It is apparent that the image is defocused with blurred and ghost targets, which are caused by channel inconsistency. After the calibration parameters are applied, a clear focused image shown in Fig. 10(c) was obtained, within which the geometric locations of all three target spheres accord with those of the true scenario, as shown in the photo of Fig. 10(a).

Fig. 10

2-D imaging experiment: (a) target scene with three metallic spheres; (b) 2-D imaging result without any calibration; and (c) 2-D imaging result after calibration.

JARS_10_2_025010_f010.png

5.1.2.

Three-dimensional imaging experiment

A 3-D imaging experiment with the same Ka-band array as in Fig. 9 is also conducted to verify the calibration method. As seen in Fig. 11, the array is vertically fixed on a horizontal rail of a ground-based platform.34 When the VNA is programmed to transmit and receive signal pulses, a synthetic aperture is realized by scanning the array along the horizontal rail and a virtual array aperture is realized by aperture synthesis of the antennas. The geometry assembles a downward-looking LA-3D-SAR in that it is simply a rotation version of the coordinates introduced in Fig. 1. In the current experiment, the array antennas synthesize a linear array with 65 virtual elements and an 8-mm inner-element spacing. The target scenario comprises five metallic spheres in front of the array and the experiment parameters are listed in Table 2.

Fig. 11

3-D imaging experiment setup for calibrating the multichannel amplitude-phase errors.

JARS_10_2_025010_f011.png

Table 2

3-D imaging experiment parameters.

ParametersValue
Center frequency35 GHz
Frequency bandwidth4 GHz
Frequency points801
Inner-element spacing8 mm
Virtual element number65
Movement interval on rail8 mm
Samples number on rail101

The phase error is calibrated by setting one point scatterer (the central metallic sphere in this experiment) as the reference target and calculating the difference between the ideal phase history along the antenna array and the real phase extracted from the 2-D compressed image corresponding to each channel. The ideal phase is obtained with the help of the relative geometric relationship between the array and the reference target and the extraction of the real phase can be performed on the peak location in the 2-D image. Then the phase error can be compensated and hence, an exact 3-D image could be expected.

In Fig. 12, we show the reconstructed 3-D images of the measurement data with and without the calibration steps, thus the importance of amplitude-phase error compensation is demonstrated. In Fig. 12(b), a 2-D image corresponding to the 40th data channel was displayed with which we can clearly see focused five spheres. The 2-D images of the other channels are similar. Figure 12(c) shows the defocused 3-D image obtained without any error compensation, whereas Fig. 12(d) gives the 3-D result after calibration. In comparison with Fig. 12(c), the five targets in Fig. 12(d) are focused in all three dimensions with higher precision. And the locations of the focused targets are in exact accordance with the target scene setup in Fig. 12(a). However, a slight defocusing still exists and this may result from the slight vibration of the array during scanning along the horizontal rail.

Fig. 12

3-D imaging experiment: (a) target scene set up with five metallic spheres; (b) 2-D imaging result corresponding to the 40th data channel; (c) 3D imaging result without any calibration; and (d) 3-D imaging result after calibration.

JARS_10_2_025010_f012.png

5.2.

Numerical Simulations of the Virtual Element Position Error Calibration

We test the calibration method with a simulated 3-D distributed target scene with the typical Ka-band LA-3D-SAR parameters listed in Table 3. The target scene comprises scatterers, whose locations and reflectivities come from a 2-D CSAR image and its digital elevation model (DEM) data,35 as shown in Fig. 13. The virtual element position errors are artificially added to the SAR raw data, where σx, σy, and σz are both set to 1 mm. Three time-divided active calibrators are used and the SNR is set to be 20 dB.

Table 3

Typical Ka-band parameters for simulation use.a

ParametersValueParametersValue
Airplane height H1500 mCenter frequency fc36.5 GHz
Array length L4 mElement number N268
Δxn, Δyn, and Δzn1 mm standard deviation

a

The parameters are the same as those given in Table 1 in Ref. 13.

Fig. 13

3-D distributed target scene simulation input: DEM data and the corresponding 2-D reflectivity image.

JARS_10_2_025010_f013.png

Figure 14 shows the reconstructed 3-D image of the simulated target scene. The image with virtual element position errors is filled with strips of energy, which represent the cross-track degradation due to a high integrated sidelobe level, making the targets in the distributed scenario indiscernible. If we apply the proposed calibration method based on time-divided active calibrators, the position errors can be calibrated within a very high accuracy of 0.02 mm, which is shown in Fig. 15, so that the reconstructed 3-D image with calibration exhibits a good performance in the same manner as the error-free case.

Fig. 14

Reconstructed 3-D images of the simulated target scene (displayed at the isosurface value of 9 dB).

JARS_10_2_025010_f014.png

Fig. 15

The estimated values of the virtual element position errors in (a) cross-track direction, (b) along-track direction, and (c) height direction as compared to their respective true values.

JARS_10_2_025010_f015.png

Subsequently, Monte Carlo simulations with 200 iterations are conducted to demonstrate the root mean square error (RMSE) of the estimator versus the SNR, as well as the number of calibrators. The results are given in Fig. 16. In Fig. 16(a), it can be seen that the estimator’s accuracy improves with the increase of SNR, and the SNR with a value of 20 dB can keep the RMSE well below 0.02 mm, satisfying the calibration requirement of the virtual element position error (with a standard deviation of 0.075 mm). In Fig. 16(b), the RMSE values corresponding to different numbers of randomly distributed calibrators under a 20-dB SNR level are given. Though the trend of RMSE improvement with the increase of the calibrator’s number is obvious, we can clearly see that using only three calibrators is quite enough for keeping the accuracy below the calibration requirement mentioned before. Consequently, using three calibrators is the cost-effective solution for the time being.

Fig. 16

RMSE versus (a) SNR and (b) number of calibrators (SNR=20  dB).

JARS_10_2_025010_f016.png

6.

Conclusions

The LA-3D-SAR is a coherent system and is particularly susceptible to array errors, such as the multichannel amplitude-phase errors and virtual element position errors. These errors are traceable to the nonideal electrical characteristics of antennas, vibration due to flight disturbance, and motion measurement errors. In this paper, the impact of the multichannel amplitude-phase errors and virtual element position errors on 3-D imaging quality with regard to their characteristics and level of severity is investigated and it has been shown that a small amount of error may cause massive degradations within the final 3-D SAR image. For the multichannel amplitude-phase errors, which increase the integrated sidelobe level in the array direction, the target response is more sensitive to the phase error than to the amplitude error. The amplitude error can be corrected by a simple power normalization while the phase error can be calibrated using external (parallel or point target) calibrators, as proposed and experimentally verified in the paper. For the virtual element position errors, it has been shown that a high integrated sidelobe level, which generates contrast degradation and noise can be caused in the cross-track direction, but with little impact on the range and along-track directions. By contrast, the cross-track beamforming performance is more sensitive to the error component in the height direction, and for typical Ka-band LA-3D-SAR parameters, the tolerable level of the error standard deviation is found to be below the submillimeter level. To mitigate the error effects, calibration and compensation measures based on time-divided active calibrators are proposed and simulated in the paper. It is worth noting that in designing an LA-3D-SAR system, as well as in processing the LA-3D-SAR data, the analyses and investigations in this paper will provide a useful reference to specify the system performance, calibration requirements, and how to eliminate the error effects and optimize the imaging result.

Acknowledgments

This work is supported by the National Natural Science Foundation for Young Scientists of China (Grant Nos. 61201433 and 61401434), the National Natural Science Foundation Project (Grant No. 61372186), Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region (Grant No. NJYT-14-B09).

References

1. 

M. Soumekh, Synthetic Aperture Radar Signal Processing with Matlab Algorithms, Wiley, New York (1999).Google Scholar

2. 

J. Shi et al., “APC trajectory design for ‘one-active’ linear-array three-dimensional imaging SAR,” IEEE Trans. Geosci. Remote Sens. 48, 1470–1486 (2010).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2009.2031430Google Scholar

3. 

Y. Cao et al., “A mathematical model for MIMO imaging,” Proc. SPIE 8393, 5 (2012).http://dx.doi.org/10.1117/12.918920Google Scholar

4. 

P. Huang and W. Xu, “ASTC-MIMO-TOPS mode with digital beam-forming in elevation for high-resolution wide-swath imaging,” Remote Sens. 7, 2952–2970 (2015).http://dx.doi.org/10.3390/rs70302952Google Scholar

5. 

W. Q. Wang, “Virtual antenna array analysis for MIMO synthetic aperture radars,” Int. J. Antennas Propag. 2010(6), 1–10 (2012).http://dx.doi.org/10.1155/2012/587276Google Scholar

6. 

H. Z. Shao and W. Q. Wang, “Multiple-input multiple-output synthetic aperture radar antenna arrangement for wide-swath remote sensing,” J. Appl. Remote Sens. 9, 094097 (2015).http://dx.doi.org/10.1117/1.JRS.9.094097Google Scholar

7. 

J. C. Yang, W. M. Su and H. Gu, “Two-dimensional superresolution multiple-input multiple-output and inverse synthetic aperture radar imaging based on spatial frequency ambiguity resolving,” J. Appl. Remote Sens. 9(1), 095090 (2015).http://dx.doi.org/10.1117/1.JRS.9.095090Google Scholar

8. 

X. Yang et al., “Amplitude and phase errors correction for array 3D SAR system based on single prominent point like target echo data,” J. Radars 4, 409–418 (2014).Google Scholar

9. 

J. Klare, A. R. Brenner and J. H. Ender, “A new airborne radar for 3D imaging-image formation using the ARTINO principle,” in Proc. 6th European Conf. on Synthetic Aperture Radar, pp. 1–4, Aachen, Germany (2006).Google Scholar

10. 

J. F. Nouvel et al., “ONERA DRIVE project,” in Proc. European Conf. on Synthetic Aperture Radar, pp. 1–3 (2008).Google Scholar

11. 

W. Tan et al., “Synthetic aperture radar tomography sampling criteria and three-dimensional range migration algorithm with elevation digital spotlighting,” Sci. China Ser. F: Inf. Sci. 52(1), 100–114 (2009).http://dx.doi.org/10.1007/s11432-009-0003-2Google Scholar

12. 

X. Peng et al., “Airborne downward looking sparse linear array 3-D SAR heterogeneous parallel simulation,” Remote Sens. 5, 5304–5329 (2013).http://dx.doi.org/10.3390/rs5105304Google Scholar

13. 

K. Han et al., “Efficient pseudopolar format algorithm for down-looking linear-array SAR 3-D imaging,” IEEE Geosci. Remote Sens. Lett. 12, 572–576 (2015).http://dx.doi.org/10.1109/LGRS.2014.2351792Google Scholar

14. 

X. Ren, L. Chen and J. Yang, “3D imaging algorithm for down-looking MIMO array SAR based on Bayesian compressive sensing,” Int. J. Antennas Propag. 2014, 1–9 (2014).http://dx.doi.org/10.1155/2014/612326Google Scholar

15. 

S. Zhang, Y. Zhu and G. Kuang, “Imaging of downward-looking linear array three-dimensional SAR based on FFT-MUSIC,” IEEE Geosci. Remote Sens. Lett. 12(4), 885–889 (2015).http://dx.doi.org/10.1109/LGRS.2014.2365611Google Scholar

16. 

J. Shi et al., “An optimal DEM reconstruction method for linear array synthetic aperture radar based on variational model,” J. Radars 4(1), 20–28 (2015).Google Scholar

17. 

W. G. Carrara et al., Spotlight Synthetic Aperture Radar: Signal Processing Algorithms, Artech House, Boston, Massachusetts (1995).Google Scholar

18. 

J. Klare et al., “Image quality analysis of the vibrating sparse MIMO antenna array of the airborne 3D imaging radar ARTINO,” in Proc. IEEE Int. Geosci. Remote Sens. Symp., pp. 5310–5314 (2007).Google Scholar

19. 

J. Klare, A. Brenner and J. Ender, “Impact of platform attitude disturbances on the 3D imaging quality of the UAV ARTINO,” in Proc. European Conf. on Synthetic Aperture Radar, pp. 1–4 (2008).Google Scholar

20. 

A. V. Oppenheim, Discrete-Time Signal Processing, 2/E, Pearson Education, India (2006).Google Scholar

21. 

K. Han et al., “Statistical analysis of the effects of virtual element position errors on airborne down-looking LASAR 3-D imaging,” IEEE Geosci. Remote Sens. Lett. 12, 1–5 (2015).http://dx.doi.org/10.1109/LGRS.2014.2332311Google Scholar

22. 

W. Philips, “ECG data compression with time-warped polynomials,” IEEE Trans. Biomed. Eng. 40, 1095–1101 (1993).IEBEAX0018-9294http://dx.doi.org/10.1109/10.245626Google Scholar

23. 

I. G. Cumming and F. H. Wong, Digital Signal Processing of Synthetic Aperture Radar Data: Algorithms and Implementation, Artech House, Norwood (2004).Google Scholar

24. 

J. R. Fienup, “Phase error correction for synthetic-aperture phased-array imaging systems,” Proc. SPIE 4123, 47–55 (2000).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.409285Google Scholar

25. 

M. Zhang and Z.-D. Zhu, “DOA estimation with sensor gain, phase and position perturbations,” in Proc. of the IEEE 1993 National Aerospace and Electronics Conf., pp. 67–69 (1993).http://dx.doi.org/10.1109/NAECON.1993.290894Google Scholar

26. 

N. Fistas and A. Manikas, “A new general global array calibration method,” in Proc. IEEE ICASSP-94 Conf., pp. IV/73–IV/76 (1994).Google Scholar

27. 

K. V. Stavropoulos and A. Manikas, “Array calibration in the presence of unknown sensor characteristics and mutual coupling,” in 10th European Signal Processing Conf., pp. 1–4 (2000).http://dx.doi.org/10.1109/ICASSP.1994.389873Google Scholar

28. 

H. Krim and M. Viberg, “Two decades of array signal processing research: the parametric approach,” IEEE Trans. Signal Process. Mag. 13, 67–94 (1996).http://dx.doi.org/10.1109/79.526899Google Scholar

29. 

G. Bienvenu and L. Kopp, “Adaptivity to background noise spatial coherence for high resolution passive methods,” in Int. Conf. on Acoustics, Speech, and Signal Processing, Vol. 5, pp. 307–310 (1980).Google Scholar

30. 

A. Freeman, “SAR calibration: an overview,” IEEE Trans. Geosci. Remote Sens. 30, 1107–1121 (1992).IGRSD20196-2892http://dx.doi.org/10.1109/36.193786Google Scholar

31. 

Y. Qixiao and S. Yonghuan, Practical Mathematics Handbook, Science Press, Beijing, China (2006).Google Scholar

32. 

G. Golub and C. V. Loan, Matrix Computations, 3rd ed., Johns Hopkins University Press (1996).Google Scholar

33. 

KeySight (Agilent Technologies), “PNA series network analyzers online help (User Guide, Programming Guide), a.10.40 ed.,” http://na.support.keysight.com/pna/help/latest/help.htm (2015).Google Scholar

34. 

W. Hong et al., “Development and experiments of ground-based SAR in IECAS for advanced SAR imaging technique validation,” in Proc. European Conf. on Synthetic Aperture Radar, pp. 1–4, Aachen, Germany (2010).Google Scholar

35. 

Y. Lin et al., “Airborne circular SAR imaging: results at P-band,” in Proc. IEEE Int. Geosci. Remote Sens. Symp., pp. 5594–5597 (2012).http://dx.doi.org/10.1109/IGARSS.2012.6352051Google Scholar

Biography

Weixian Tan is an associate professor at the Inner Mongolia University of Technology, Hohhot, China. He received his PhD in signal and information engineering from the Institute of Electronics, Chinese Academy of Sciences, in 2009. He is the author of more than 50 journal and conference papers. His current research interests include microwave imaging, holography, and synthetic aperture radar systems.

Pingping Huang received his PhD from the Institute of Electronics, Chinese Academy of Sciences, Beijing, China, in 2010. He is currently with the College of Information Engineering, Inner Mongolia University of Technology, Hohhot, China. His current major research interests are signal processing, digital beamforming, polarimetric interferometry, and spaceborne synthetic aperture radar system design.

Kuoye Han received his PhD from the Institute of Electronics, Chinese Academy of Sciences, Beijing, China, in 2015. He is currently with the Information Science Academy, China Electronics Technology Group Corporation. His current major research interests are phased-array radar, synthetic aperture radar system design, digital beamforming, and microwave imaging.

Qi Liu is a senior engineer at China Institute of Water Resources and Hydropower Research, Beijing, China. She received her PhD in signal and information engineering from the Institute of Electronics, Chinese Academy of Sciences, in 2009. Her current research interests include microwave imaging, water information, and rain radar.

Xueming Peng received his bachelor’s degree in electrical engineering and automation from Anhui University, Hefei, China, in 2009, and his PhD from the University of Chinese Academy of Sciences, Beijing, China, in July 2014. From August 2014 to May 2015, he was with the No. 802 Research Institute, Shanghai Academy of Spacelight Technology (also known as Shanghai Radio Equipment Institute). He is currently a researcher with the Global R&D Centre (Shanghai), Carestream Health, Inc., Shanghai, China. His main research interests are radar signal processing, heterogeneous parallel computing and real-time operating systems, and 3-D polygon mesh processing.

WeiXian Tan, Pingping Huang, Kuoye Han, Qi Liu, Xueming Peng, "Array error calibration methods in downward-looking linear-array three-dimensional synthetic aperture radar," Journal of Applied Remote Sensing 10(2), 025010 (2 May 2016). http://dx.doi.org/10.1117/1.JRS.10.025010
Submission: Received ; Accepted
JOURNAL ARTICLE
20 PAGES


SHARE
KEYWORDS
Calibration

Error analysis

Antennas

3D image processing

Synthetic aperture radar

Chemical elements

3D acquisition

Back to Top