International Journal of Antennas and Propagation

Volume 2014 (2014), Article ID 161294, 8 pages

http://dx.doi.org/10.1155/2014/161294

Research Article

Correction of Channel Imbalance for MIMO SAR Using Stepped-Frequency Chirps

Xiulian Luo,1,2 Yunkai Deng,¹ Robert Wang,¹ Lei Guo,1,2 and Mingjiang Wang1,2

¹Institute of Electronics, Chinese Academy of Sciences, Beijing 100190, China

²Graduate University of the Chinese Academy of Sciences, Beijing 100080, China

Received 20 December 2013; Accepted 19 February 2014; Published 3 April 2014

Academic Editor: Wei Xu

Copyright © 2014 Xiulian Luo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

To simultaneously achieve two-dimensional high resolution and wide swath in synthetic aperture radar (SAR), azimuth MIMO structure combined with stepped-frequency chirp signals was developed via splitting the antenna into subapertures. During transmitting each subaperture transmits a chirp pulse at a different carrier frequency, while during receiving every subaperture receives the scattered pulses at the same time. Separating the scattered pulses received by each subaperture and downlinking them to the ground yield different signal paths. Due to the dedicated network in the SAR system, the channel imbalance is inevitable. To correct the channel imbalance, this paper presents an external calibration method, where the channel characteristics are estimated from the peak value of a strong point target for each channel. Simulation and real raw data experiments are performed to validate the proposed method.

1. Introduction

Resolution and swath width are two key specifications for the spaceborne stripmap SAR. However, high azimuth resolution and wide range swath imaging pose contradicting requirement on conventional SAR system design. In order to overcome this inherent contradiction, the displaced phase center antenna (DPCA) techniques [1–3] and the multiple-input and multiple-output (MIMO) antenna techniques [4–7] are introduced to acquire additional spatial sampled information. On the other hand, the range resolution, determined by the pulse bandwidth, has to be increased to match the azimuth resolution. However, transmitting and receiving wideband signals burden the active electronically steered array (AESA) antenna and the echo sampled equipment [8, 9], thus confining the range resolution. Fortunately, a wideband chirp signal can be synthesized by a group of narrowband subchirps centered at stepped frequencies, which requires the synthetic bandwidth technique [8–10].

As the synthetic bandwidth technique will reduce the number of equivalent phase centers (EPCs), the additional sampled information in DPCA SAR is actually used to increase the range bandwidth other than the equivalent PRF [11]. Therefore, this paper focuses on the MIMO SAR using stepped-frequency chirps (STFC) [4, 5], where both the equivalent PRF and the range bandwidth can be times increased. However, as this MIMO system yields different signal paths, the channel imbalance is inevitable. For each subband, the collected data is from a DPCA system, where the channel imbalance can be regarded as the inner-band imbalance which results in azimuth ambiguities [12]. The imbalance among different subbands can be regarded as the interband imbalance which leads to degradation of the range performance after synthetic bandwidth [13].

Therefore, in order to acquire high quality images, effective correction of the channel imbalance is indispensable. The available correction methods can be classified as two categories: methods on internal calibration [14] and methods on raw data [15, 16]. However, all these correction methods only either handled the inner-band imbalance in the DPCA system with single subband [15, 16] or concerned the multiple subbands system with single subaperture [14]. In [14], Deng et al. gave the internal calibration approach, where the imbalance information is obtained from the internal calibration data acquired from the calibration subsystem. This method can be extended to the MIMO system using STFC. Nevertheless, internal calibration subsystems are required, which dramatically increases the complexity of the SAR systems. In addition, it cannot correct imbalance introduced by the antennas. The subspace projection method proposed in [16] and the azimuth cross-correlation method presented in [15] do not require additional subsystem. However, they cannot correct the interband imbalance, as the phase difference caused by different carrier frequencies is range dependent and hard to compensate.

Usually, the channel imbalance information (phase and amplitude) does not change within a certain illumination time [15]. Moreover, many artificial corner reflectors have been distributed here and there for radiometric calibration [17]. Therefore, this paper proposes an external calibration method where the channel imbalance information is estimated from the peak value of a corner reflector or a strong point target in small complex images. Afterwards, the estimated channel errors are used to correct the channel imbalance when imaging for much larger data.

Succeeding sections are organized as follows. In Section 2, the MIMO SAR using STFC is conceptually designed and the corresponding image formation processing is briefly summarized. In Section 3, the external correction method is presented, followed by the simulation and real raw data experiments in Section 4 to validate the proposed method. Finally, the conclusion is drawn in Section 5.

2. MIMO SAR Using Stepped-Frequency Chirps

2.1. Conceptual Design

The MIMO SAR simultaneously transmits a set of stepped-frequency chirps, regarded as subbands, via multiple subapertures. During receiving, every subaperture receives all the subbands in one PRI. The schematic diagram and the corresponding timeline of transmitted pulses and received echoes are illustrated in Figure 1, where the circles represent the location of the equivalent phase centers (EPCs) of the subbands, and the number inside the circles represents the sequence number of the subbands.

Figure 1: Schematic diagram of the MIMO SAR using STFC.

In order to separate two frequency-adjacent subbands and absent gaps from the connection of subbands at the same time, the subband bandwidth is required to be equal to the frequency step . The relationship of the time frequency for the transmitted pulses is illustrated in Figure 2, where () is the carrier frequency of the th subband.

Figure 2: Relationship of time frequency of transmitted pulses.

The baseband echo of the th subband for each subaperture is acquired as follows. Demodulate the radio frequency (RF) pulses with the carrier frequency . Then filter the demodulated signal with an analog low pass filter (LPF) whose bandwidth is greater than . Afterwards, convert the analog signal to a digital signal via an AD convertor with a sampling rate greater than the cutoff frequency of the LPF. As illustrated in Figure 3, due to the imperfect filtering and the Gibbs phenomenon of the signal spectrum, some spectral component of one subband is leaked into the adjacent subband, thus resulting in subband interference. In order to reduce the interference, a digital rectangular window can be applied to the sampled data before imaging. Then, the interference only results from the spectral component outside the spectrum bandwidth, as illustrated by the blue-shaded region in Figure 3. As stated in [13], the influence of the subband interference under this situation can be ignored because it results in −47 dB ambiguities in range. In fact, when the subband bandwidth is larger, the ambiguities in range will be reduced further. This is because the relative amplitude of the ambiguities is determined by relating the leakage energy to the whole energy [13]. Moreover, the leakage energy is inversely proportional to the bandwidth-time product.

Figure 3: Separation of two frequency-adjacent subbands.

Therefore, in this paper, for each subaperture, adjacent subbands are assumed to be completely separated; that is, the subband interference is not considered.

2.2. Image Formation Processing

When the uniform sampling condition

with the system pulse repetition frequency, the length of the subaperture, and the velocity of the sensor, is not satisfied, multichannel reconstruction [1] is required. The reference position during the multichannel reconstruction should be the center of the whole antenna to make the EPCs to be synthesized coincide spatially. Afterwards, the frequency domain bandwidth synthesis method [9, 10] is applied to acquire a large bandwidth pulse, followed by the conventional imaging. It should be noted that the synthesized signal should be firstly converted to the original echo domain if the CS imaging algorithm is adopted. The block diagram of the image formation processing for this MIMO SAR system is presented in Figure 4. Assume that , and then the variation of the distribution of EPCs introduced by multichannel reconstruction and bandwidth synthesis is shown in Figure 5 (ignore the first row and the dashed circles, which will be used in Section 4).

Figure 4: Block diagram of the image formation processing.

Figure 5: Distribution of the EPCs.

3. Correction of the Channel Imbalance

In this section, the correction of both amplitude and phase imbalance among channels will be presented. As this correction method is based on a focused strong point target in small complex images, separating imaging for each subband and each subaperture is required.

The base-band echo for the th subband (transmitted by the th subaperture) received by the th subaperture is

with the range time, the azimuth time, the chirp rate, the light velocity, the pulse duration, the amplitude characteristic, the phase characteristic, the system noise and the echo of other targets, and the instantaneous distance in propagation between the transmitter and the receiver for the corner reflector (point target). Applying the Taylor series expansion and ignoring the quadratic phase termquadratic term error during the equivalence of the two-way slant ranges, can be approximated as

where

is the azimuth time interval between the corresponding EPC and the center of the whole antenna, and

represents the two-way slant range from the center of the whole antenna to the point target. In (5), and indicate the slant range and azimuth coordinates of the point target, respectively.

First, the imaging for the echoes is independently performed. Due to lower PRF than Nyquist sampling rate, strong azimuth ambiguities of the point target appear in the images. However, the ambiguities are far away from the real point target, thus scarcely impacting the estimation. The imaging result can be approximated as [18]

where is the azimuth chirp rate, and are the gain resulting from the imaging of the “real” part and the ambiguous part of the point target, respectively, and are the “real” and ambiguous response of other targets and noise, and and are the amplitude of the azimuth and range impulse responses, two sinc-like functions [18]. According to (3), can be rewritten as

where is also a sinc-like function.

Considering that the peak of the corner reflector is generally much greater than other terms at position , only the “real” response of the corner reflector is taken into account when deriving the correction method, which is

Second, the peak position and value of the corner reflector are estimated from these small complex images. In order to make the estimation more accurate, the 2D-interpolation over the peak point and its surrounding points is performed. Consequently, the peak position of the corner reflector is acquired to be . The amplitude of the peak is , and the phase of the peak is

In fact, in order to reduce the influence of the noise, and are obtained by averaging the values of the points within the two-dimensional 1 dB main lobe. The first channel (subband 1 and subaperture 1) is set as the reference channel. Then the estimated amplitude error matrix is

Third, compensate the inherent phase caused by the offset of the carrier frequencies from the center frequency of the synthesized wideband signal, . After this, the phase of the peak value is changed as

where

Therefore, the estimated phase error matrix can be expressed as

Finally, correct the channel imbalance according to (10) and (13). The correction matrix is

Therefore, the correction is implemented as

It should be noted that the data size of in (2) can be much smaller than that in (15) to improve the processing efficiency. After correction, the balanced data is used to form the final image with high resolution. The whole processing is shown in Figure 6, where the image formation processing strategy is presented in Figure 4.

Figure 6: Block diagram of the channel imbalance correction.

4. Experimental Results

4.1. Simulation Experiment

In the simulation, the 2D imaging of a plane scene containing a letter “A” and a 4-time (6 dB) brighter point target is performed, and this point target is inside the letter “A.” The simulation parameters are specified in Table 1.

Table 1: System parameters.

Parameter
Symbol
Value

Center frequency (GHz)

9.685

Slant range of scene center (km)

30

Sensor velocity (m/s)

215

Subaperture length (m)

2.5

System PRF (Hz)
PRF
140

Transmit pulse length (us)

10

Number of subbands/subapertures

2

Transmitted pulse bandwidth (MHz)

60

Sampling rate in range (MHz)

72

SNR (dB)
SNR
6

The accurate amplitude and phase errors are given in the 2nd and 3rd columns of Table 2, respectively. The first channel is set as the reference channel. According to the estimation method proposed in the last section, the amplitude and phase errors are estimated and then presented in the 4th and 5th columns of Table 2, respectively. The imaging results without and with correction of the channel imbalance are shown in Figures 7(a) and 7(b), respectively. The range comparison of the point target and its surroundings is shown in Figure 7(c). From Figure 7, one can observe that the ambiguities in azimuth (caused by the inner-band imbalance) and the raising of the side-lobes in range (stemming from the interband imbalance) are removed by the correction of channel imbalance.

Table 2: Channel errors (CE).

Channel (, )
Exact CE
Estimated CE for simulated data
Estimated CE for constructed raw data

(°)

(°)

(°)

1, 1
1
0
1
0
1
0

1, 2
1.3
25
1.296
24.308
1.253
28.201

2, 1
1.5
30
1.500
29.649
1.468
33.325

2, 2
1.4
45
1.401
44.426
1.352
42.720

(a)

(b)

(c)

Figure 7: Simulation results for OFDM-MIMO SAR. (a) Imaging result without correction of channel imbalance. (b) Imaging result with correction of channel imbalance. (c) Range comparison of the amplification of the point target in (a) and (b).

From the simulation experiment, it is concluded that the proposed correction (compensation) method can correct both the inner-band and the interband imbalance for the simulated MIMO SAR data.

4.2. Real Raw Data Experiment

For the real raw data, the imaging scene is much more complicated than the simulated data, thus possibly affecting the estimation accuracy. As the real raw data collected by a MIMO SAR using STFC is not available, in this subsection, the proposed correction method is assessed with an X-band airborne SAR using a single aperture and two STFCs which are transmitted in two consecutive PRIs. In this SAR, both the subband bandwidth and the frequency step are 30 MHz, the original PRF is 1600 Hz, and the sensor velocity is 215 m/s. The distribution of the original EPCs is illustrated by the first row in Figure 5.

To validate the proposed correction method on the real raw data, the echo data for the MIMO system have to be firstly constructed. As discussed before, in the MIMO system using STFCs, the echo for each subband can be regarded from a DPCA system [1]. Further, in the DPCA system, the data of each channel is equivalent to the counterpart single channel data processed by an azimuth time delay function and then times subsampling [2]. Herein, the counterpart single channel data is the original subband data, and the delay function in Doppler domain for the th subband received by the th subaperture is set to

with the Doppler frequency, , , , and . Then, the distribution of the EPCs for the four channels can be illustrated by the 2nd, 3rd, 4th, and 5th rows in Figure 5 where the dashed circles represent the discarded pulses to realize a two-time subsampling. With consideration of the channel errors specified in the 2nd and 3rd columns of Table 2, the four echo data of the MIMO system using STFCs are well constructed.

According to the error estimation process presented in Figure 6, the measured channel errors are listed in the 6th and 7th columns of Table 2. The final imaging results without and with correction of the channel imbalance are shown in Figures 8(a) and 8(b), respectively. The corner reflector used for external calibration is marked by ellipse in the two figures. After interpolation and amplification, the range comparison of the corner reflector in Figures 8(a) and 8(b) is shown in Figure 8(c), where one can see that the raising of the side lobes is eliminated by the compensation. Additionally, in Figure 8(a) the peak amplitude of the ambiguity marked by dashed circle is −15.89 dB while in Figure 8(b) the ambiguity is not distinguishable from the background.

(a)

(b)

(c)

Figure 8: Real raw data results for OFDM-MIMO SAR. (a) Without correction. (b) With correction. (c) Range comparison of the amplification of the corner reflector.

Therefore, for the constructed MIMO raw data, the presented correction method is also effective in removing the channel imbalance.

5. Conclusion

MIMO SAR using stepped frequency chirps is of great potential for future high-resolution wide-swath SAR missions. Channel imbalance is a key problem for this novel SAR mode, as it relates to the imaging performance. This paper proposed a simple and effective method to correct this channel imbalance, which was validated by simulation and real raw data experiments. Other problems, such as the system design and the orbit model, are not discussed in this paper. Nevertheless, they are also indispensable, especially when this novel mode is applied to spaceborne missions.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

N. Gebert, G. Krieger, and M. A. Moreira, “Digital beamforming on receive: techniques and optimization strategies for high-resolution wide-swath SAR imaging,” IEEE Transactions on Aerospace and Electronic Systems, vol. 45, no. 2, pp. 564–592, 2009.

W. Jing, M. Xing, C.-W. Qiu, Z. Bao, and T.-S. Yeo, “Unambiguous reconstruction and high-resolution imaging for multiple-channel SAR and Airborne experiment results,” IEEE Geoscience and Remote Sensing Letters, vol. 6, no. 1, pp. 102–106, 2009.

A. Currie and M. A. Brown, “Wide-swath SAR,” IEE Proceedings, Part F: Radar and Signal Processing, vol. 139, no. 2, pp. 122–135, 1992.

W.-Q. Wang, “Space-time coding MIMO-OFDM SAR for high-resolution imaging,” IEEE Transactions on Geoscience and Remote Sensing, vol. 49, no. 8, pp. 3094–3104, 2011.

Y. Song and R. Yang, “High resolution, wide swath SAR using sub-aperture sub-band technique,” in Proceedings of the CIE International Conference on Radar (ICR '06), Shanghai, China, October 2006.

G. Krieger, “MIMO-SAR: opportunities and pitfalls,” IEEE Transactions on Geoscience and Remote Sensing, no. 99, pp. 1–18, 2013.

J. H. G. Ender and J. Klare, “System architectures and algorithms for radar imaging by MIMO-SAR,” in Proceedings of the IEEE Radar Conference (RADAR '09), Pasadena, Calif, USA, May 2009.

A. W. Doerry, SAR Processing with Stepped Chirps and Phased Array Antennas, Sandia Nat. Lab., Albuquerque, NM, USA, 2006.

P. Berens, “SAR with ultra-high range resolution using synthetic bandwidth,” in Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS '99), pp. 1752–1754, Hamburg, Germany, July 1999.

W. Nel, J. Tait, R. Lord, and A. Wilkinson, “The use of a frequency domain stepped frequency technique to obtain high range resolution on the CSIR X-Band SAR system,” in Proceedings of the 6th IEEE AFRICON Conference in Africa, pp. 327–332, October 2002.

J. H. G. Ender and A. R. Brenner, “PAMIR—a wideband phased array SAR/MTI system,” IEE Proceedings: Radar, Sonar and Navigation, vol. 150, no. 3, pp. 165–172, 2003.

C.-G. Gao, Y.-K. Deng, and J. Feng, “Theoretical analysis on the mismatch influence of displaced phase center multiple-beam SAR systems,” Journal of Electronics & Information Technology, vol. 33, no. 8, pp. 1828–1832, 2011.

X. Luo, R. Wang, Y. Deng, and W. Xu, “Influences of channel errors and interference on the OFDM-MIMO SAR,” in Proceedings of the IEEE Radar Conference, pp. 1–5, Ottawa, Canada, 2013.

Y. Deng, H. Zheng, R. Wang, J. Feng, and Y. Liu, “Internal calibration for stepped-frequency chirp SAR imaging,” IEEE Geoscience and Remote Sensing Letters, vol. 8, no. 6, pp. 1105–1109, 2011.

J. Feng, C. Gao, Y. Zhang, and R. Wang, “Phase mismatch calibration of the multichannel SAR based on azimuth cross-correlation,” IEEE Geoscience and Remote Sensing Letters, vol. 10, no. 4, pp. 903–907, 2013.

L. Zhang, M.-D. Xing, C.-W. Qiu, and Z. Bao, “Adaptive two-step calibration for high-resolution and wide-swath SAR imaging,” IET Radar, Sonar and Navigation, vol. 4, no. 4, pp. 548–559, 2010.

J. C. Curlander and R. N. Mcdonough, Synthetic Aperture Radar: Systems and Signal Processing, John Wiley & Sons, New York, NY, USA, 1991.

I. G. Cumming and F. H. Wong, Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation, Artech House, Norwood, Mass, USA, 2005.