|
1.IntroductionIn 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) array3–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.11–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-SARAs shown in Fig. 1, the radar is mounted on an airplane flying at altitude H with the velocity . Let the -, - and -axes denote the cross-track, along-track, and height directions, respectively. In cross-track, a sparse MIMO array with transmitters and receivers is adopted to form a uniform virtual array, which has virtual elements and total length , where represents the element spacing. The position of the ’th virtual element is , where represents the along-track spatial sample coordinate. The illuminated scenario is located at the nadir area of the platform, and for a point target with reflectivity located at in the 3-D target support band , the instantaneous target-to-’th virtual element distance is Assuming that orthogonal pulses with pulse length , frequency bandwidth , and carrier frequency 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: where denotes the fast-time variable, denotes the along-track synthetic aperture length, and 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 is directly accessible, which is where is the fast-time wavenumber variable, with , , and .Usually, the process of imaging is to estimate the distribution of the target reflectivity from the range compressed signal or its wavenumber domain signal . 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 Imaging3.1.Multichannel Amplitude-Phase ErrorsThe 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, and represent the cross-track and slant range directions, respectively. Assume each virtual element’s position is , where . Using to replace the instantaneous range in Eq. (1) and taking the multichannel amplitude-phase error into consideration, the signal corresponding to the ’th channel can be expressed as where denotes the gain variations due to the amplitude error and denotes the phase error, which can be divided into two parts: the first term representing the different time delay (or shift) components for each channel is referred to as range deviation phase error, and the other part is referred to as residue phase error.For almost all the imaging algorithms, the cross-track focusing is done by coherent accumulation along the array direction, i.e., 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 and . 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 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 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 , center frequency GHz, and the random range deviation corresponding to each channel conforms to the uniform distribution, i.e., , , where is a constant. Figure 3(a) shows the curve where the cross-track ISLR changes with parameter . It is apparent that 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. The residue phase error 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 , and the ideal target response as , then according to the Fourier transform equation,20 the practical signal can be processed after coherent accumulation, which yields where * denotes the discrete convolution. Equation (6) implies that the ideal target impulse response spreads over the axis because of the presence of . Again, using the ISLR evaluation methodology, except that the residue phase error is set to conform to a uniform distribution , . As shown in Fig. 3(b), the target response is deteriorated with the increase of the residue phase error. For example, when distributes randomly between 0 and , the level of ISLR has been creeping up from to , where correct imaging is impossible to achieve.3.2.Virtual Element Position ErrorsThe 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 from the ideal situation , where , , and 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 ’th virtual element is represented as , so the actual target-to-sensor distance is From Eqs. (2) and (7), and using , we can write the range compressed impulse response as Taking the binomial expansion of Eq. (8), which yields where is the minimum range between the point target to the ’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 , results inThe 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 , and represent phase error terms resulting from the virtual element position errors , , and , respectively, which are defined as 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 , , and 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 , , and as discrete, white, Gaussian stochastic process with zero-mean and variance of , , and , respectively. , , and 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 is a compound of linear phase error in the along-track direction and random phase error in the cross-track direction, whereas and 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 denote the ’th order orthogonal expansion term coefficient of , and as its variance. Figure 5(d) implies that to meet a typical ISLR demand for SAR imaging,23 has to be kept below . 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 , which are evaluated with four different point targets. For targets far away from the array, the coefficient variance is large and a coefficient variance corresponds to with the value of 0.5 mm, which is taken as the calibration precision requirement of virtual element position error . 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 of changing with the standard deviation of the virtual element position error component , again evaluated with four point targets located at different positions. It is apparent that the four curves almost overlapped, and a coefficient variance corresponds to value of approximately 0.075 mm. Some calibration or compensation procedures would be necessary if exceeds this limit. Compared with , the beamforming performance along the cross-track is more sensitive to , 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 . When is reduced to 0.5 mm, as shown in Fig. 6(b), the ISLR is reduced to an acceptable level of 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 direction will still spread out energy across the impulse response and give rise to an ISLR of . In Fig. 6(d), when the is set to 0.075 mm, a favorable cross-track impulse response with ISLR of 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 should be taken as a reference for calibrating the virtual element position errors. 4.Calibration of the Array Errors4.1.Calibration of the Multichannel Amplitude-Phase Errors4.1.1.Amplitude error calibration using power normalizationFor 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 , 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 ’th channel is , with the power value of and fast-time samples number , then power normalization can be done through where is the ’th channel data after amplitude error calibration.4.1.2.Phase error calibration using external calibratorsThe 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. 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 of the calibrator in all channels are estimated by cross-correlating the data with the range point spread function : Then the difference of and the ideal distance value , i.e., , is taken as the estimated value of the range deviation error corresponding to the ’th channel: After correcting the range deviations, the energy trajectory of the calibrator locates at the range gate of and manifests as a parallel line with respect to the array. Then the phase values , on the calibrator response peak can be extracted as the estimated residue phase errors, and thus can be corrected throughThe 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 CalibratorsFor 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.25–27 In the downward-looking LA-3D-SAR imaging model, , , and 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 , the position errors are deterministic and constant. For a certain , 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 , 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) is located at in Fig. 8, where and represent the along-track and cross-track looking angle with respect to the cross-track array. At time , let denote the along-track position of the LA-3D-SAR, the transmitting antennas emit orthogonal waveforms , the calibrator relays the signals and later, the signals are received by all receiving antennas simultaneously. By aperture synthesis, a total number of virtual elements are formed and the channel data corresponding to the ’th virtual element are where is the fast-time range variable, is the center wave number, represents additive channel noise which is a Gaussian, white stochastic process and independent of the signal components, is the instantaneous range from the ’th virtual element to the calibrator, with the expression:Let denote the range from the calibrator to the reference element, then taking the Taylor series expansion of ,31 and neglecting the small quadratic and higher-order terms, yields Inserting Eq. (18) into Eq. (16), and transforming the result into vector-matrix form, which yields where where represents vector transpose, represents a diagonal matrix whose diagonal is vector , 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 is where 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 :In Eq. (28), represents the signal subspace spanned by the eigenvectors corresponding to big eigenvalues, whereas 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 denote this eigenvector, then where is a real constant. is the estimated steering vector that contains virtual element position errors. From Eq. (22)–(24), we can estimate the diagonal matrix by solving the following optimization problem:Let and denote the ’th elements of vectors and , respectively, i.e., and , and considering that is a diagonal square matrix, i.e., , then it is easy to get Because the above expression is non-negative, the minimization is achieved when each additional term of in Eq. (31) becomes 0, which means that the estimated value of is , so that for which , all have nonzero values according to the definition of in Eq. (21).Suppose there are a total number of time-divided active calibrators in the scene. For each calibrator, we can obtain an estimated diagonal matrix , let denote the vector formed by the diagonal elements of matrix , then the virtual antenna element position errors can be estimated by solving the matrix equation: where where is the operator for obtaining the angle of a complex quantity. Because contains unknown variables and each calibrator can provide () equations to the linear system, then there should exist at least three calibrators to solve the matrix equation. If , the unknown variables can be determined uniquely by multiplying the inverse of a nonsingular matrix: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., , we say the linear system is overdetermined and thus can be solved by the least square method.32 However, if , 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 Results5.1.Indoor Experiments of the Multichannel Amplitude-Phase Error Calibration5.1.1.Two-dimensional imaging experimentTo 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. Table 12-D imaging experiment parameters.
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). 5.1.2.Three-dimensional imaging experimentA 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. Table 23-D imaging experiment parameters.
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. 5.2.Numerical Simulations of the Virtual Element Position Error CalibrationWe 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 , , and are both set to 1 mm. Three time-divided active calibrators are used and the SNR is set to be 20 dB. Table 3Typical Ka-band parameters for simulation use.a
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. 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. 6.ConclusionsThe 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. AcknowledgmentsThis 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). ReferencesM. Soumekh, Synthetic Aperture Radar Signal Processing with Matlab Algorithms, Wiley, New York
(1999). Google Scholar
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). http://dx.doi.org/10.1109/TGRS.2009.2031430 IGRSD2 0196-2892 Google Scholar
Y. Cao et al.,
“A mathematical model for MIMO imaging,”
Proc. SPIE, 8393 5
(2012). http://dx.doi.org/10.1117/12.918920 Google Scholar
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/rs70302952 Google Scholar
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/587276 Google Scholar
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.094097 Google Scholar
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.095090 Google Scholar
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
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,
1
–4
(2006). Google Scholar
J. F. Nouvel et al.,
“ONERA DRIVE project,”
in Proc. European Conf. on Synthetic Aperture Radar,
1
–3
(2008). Google Scholar
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-2 Google Scholar
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/rs5105304 Google Scholar
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.2351792 Google Scholar
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/612326 Google Scholar
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.2365611 Google Scholar
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
W. G. Carrara et al., Spotlight Synthetic Aperture Radar: Signal Processing Algorithms, Artech House, Boston, Massachusetts
(1995). Google Scholar
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.,
5310
–5314
(2007). Google Scholar
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,
1
–4
(2008). Google Scholar
A. V. Oppenheim, Discrete-Time Signal Processing, 2/E, Pearson Education, India
(2006). Google Scholar
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.2332311 Google Scholar
W. Philips,
“ECG data compression with time-warped polynomials,”
IEEE Trans. Biomed. Eng., 40 1095
–1101
(1993). http://dx.doi.org/10.1109/10.245626 IEBEAX 0018-9294 Google Scholar
I. G. Cumming and F. H. Wong, Digital Signal Processing of Synthetic Aperture Radar Data: Algorithms and Implementation, Artech House, Norwood
(2004). Google Scholar
J. R. Fienup,
“Phase error correction for synthetic-aperture phased-array imaging systems,”
Proc. SPIE, 4123 47
–55
(2000). http://dx.doi.org/10.1117/12.409285 PSISDG 0277-786X Google Scholar
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.,
67
–69
(1993). http://dx.doi.org/10.1109/NAECON.1993.290894 Google Scholar
N. Fistas and A. Manikas,
“A new general global array calibration method,”
in Proc. IEEE ICASSP-94 Conf.,
IV/73
–IV/76
(1994). Google Scholar
K. V. Stavropoulos and A. Manikas,
“Array calibration in the presence of unknown sensor characteristics and mutual coupling,”
in 10th European Signal Processing Conf.,
1
–4
(2000). http://dx.doi.org/10.1109/ICASSP.1994.389873 Google Scholar
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.526899 Google Scholar
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,
307
–310
(1980). Google Scholar
A. Freeman,
“SAR calibration: an overview,”
IEEE Trans. Geosci. Remote Sens., 30 1107
–1121
(1992). http://dx.doi.org/10.1109/36.193786 IGRSD2 0196-2892 Google Scholar
Y. Qixiao and S. Yonghuan, Practical Mathematics Handbook, Science Press, Beijing, China
(2006). Google Scholar
G. Golub and C. V. Loan, Matrix Computations, 3rd ed.Johns Hopkins University Press(1996). Google Scholar
KeySight (Agilent Technologies),
“PNA series network analyzers online help (User Guide, Programming Guide),”
2015). http://na.support.keysight.com/pna/help/latest/help.htm Google Scholar
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,
1
–4
(2010). Google Scholar
Y. Lin et al.,
“Airborne circular SAR imaging: results at P-band,”
in Proc. IEEE Int. Geosci. Remote Sens. Symp.,
5594
–5597
(2012). http://dx.doi.org/10.1109/IGARSS.2012.6352051 Google Scholar
BiographyWeixian 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. |