# Channel phase error compensation for MIMO-SAR.

1. IntroductionSynthetic aperture radar (SAR) is a powerful microwave instrument for remote sensing [1]. 2D even 3D [2] image of the terrestrial surface can be obtained independent of weather and sunlight illumination. High-resolution as well as wide swath (HRWS) is one of the main goals of the system engineers [3, 4].

For the sake of high azimuth resolution, people have developed spotlight mode [5, 6] in which the antenna points toward the same region during the whole data acquisition. While obtaining large imaging swath, ScanSAR mode [7] is proposed in which the antenna directivity is fixed, swapping between different strips. Limited by the minimum-antenna-area constraint [1], conventional spaceborne SAR is hard to achieve high-resolution in azimuth meanwhile wide imaging swath. High azimuth resolution requires large Doppler bandwidth, which means a high pulse repetition frequency (PRF) to comply with the Nyquist theory, whereas, wide imaging swath demands wide range beam width in which the PRF must be low to guarantee the echo's completeness. Modes such as sliding-spotlight and TOPS [8, 9] are the trade-off between the two demands.

In range direction, SAR transmits chirp pulse. The resolution is proportional to the signal bandwidth, which is defined as

[[sigma].sub.r] = c/2B, (1)

where [[sigma].sub.r] is the range resolution, c is the velocity of light, and B is the signal bandwidth. For 0.1m resolution in slant range, 1.5GHz bandwidth is necessary. This implies at least 1.8GHz sampling rate (with oversampling rate 20%). The high quality, that is, radiation hardened, A/D convertor is quite difficult to manufacture.

Multi-input multioutput (MIMO) technique provides an opportunity to break through these constraints [10, 11]. The antenna being divided into several apertures, each aperture transmits pulses at different center frequencies simultaneously and all apertures receive the echoes. In this way, numerous effective phase centers are formed. These effective phase centers may be overlapped in spacetime or in frequency. Then using the signal processing technique signals with large bandwidth (both in range and in Doppler) as well as wide swath are obtained [12, 13].

Channel imbalance is inevitable in multichannel system due to the effects of the qualities of modules, the construction of the system, and the atmosphere. The imbalance presents in echo data in terms of amplitude error and phase error. These errors impair the signal synthesis in which the final SAR image will be ambiguous even unrecognized.

Channel imbalance can be removed by means of system internal calibration [1] or analyzing through raw data. Methods based on internal calibration are precise but cannot calibrate the phase errors caused by the antenna, whereas literature belonging to the latter is rare to see. In this paper, an algorithm based on raw data to handle channel phase imbalance in MIMO-SAR system is proposed. The method models phase imbalance as two parts: the transmit phase error (TXE) and the receive phase error (RXE). The transmit phase error is estimated by correlation operation, whereas the receive phase error is removed using cost function at the azimuth processing stage.

This paper is organized as follows. In Section 2, the signal mode is introduced. The proposed method is presented in Section 3 with simulations and discussions in Section 4. Section 5 will give the conclusion.

2. Signal Mode

Figure 1 is a diagram of a MIMO-SAR system with five channels in azimuth. Antennas of each channel and the corresponding phase-centers are marked with rectangles and triangles, respectively. Each channel transmits pulse simultaneously and receives the echoes transmitted by not only itself but also the others. Suppose channel m transmits pulse at center frequency [f.sub.c,m]. By compensating a known constant phase between channels, the received echo-signal can be converted into the equivalent self-transmit and self-receive signal [14]. In Figure 1(b), the effective phase-centers are marked with circles which are located halfway in between the transmitting and the respective receiving phase-centers. The effective phase-centers of different center frequencies are marked with distinct colors. Theoretically, if a MIMO-SAR transmits pulse with M channels and receives with N, the total effective phase-centers is Mx N.

The signal transmitted by channel m is

[s.sub.t,m]([tau], [eta]) = rect ([tau]/T) exp {j2[pi] ([f.sub.c,m][tau] +1/2 [K.sub.r][[tau].sup.2])} , (2)

where [tau] is the fast time, [eta] is the slow time, T is the pulse length, and [K.sub.r] is the chirp rate. Then the signal after mixing with center frequency [f.sub.c,m] received by channel n is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)

where [sigma] is the backscattering coefficient of the target, c is the velocity of light, and [R.sub.n,m]([eta]) is the slant range

[R.sub.n,m] ([eta]) = [square root of ([R.sup.2.sub.0] + [([V.sub.r][eta] + (n - 1) x d/2 + (m - 1) x d/2).sup.2])], (4)

where [R.sub.0] is the nearest range, [V.sub.r] is the platform velocity, d is the antenna length in azimuth of a single channel, and m and n ranges from 1 to M. In (4), the reference channel is channel 1; that is, the reference slant range is [R.sub.1,1]([eta]) = [square root of ([R.sup.2.sub.0] + [V.sup.2.sub.r][[eta].sup.2])].

Consider the signals with a specific center frequency. They can be treated as the azimuth multichannel SAR signals. Usually, all channels have the same the antenna lengths. According to the slow time delay, after compensating the known constant phase, the received signals have

[s.sub.r,n,m] ([tau], [eta]) = [s.sub.r,1,m] ([tau], [eta] + (n - 1) x d/2[V.sub.r]). (5)

If the radar platform velocity and the PRF have a relationship

PRF =Md/2[V.sub.r], (6)

where M is the number of receive channels and d is the antenna length of a single channel, the effective phase-centers are distributed uniformly in azimuth. However, (6) is hard to meet strictly in practical SAR systems because of the timing restriction, and therefore the effective phase-centers are nonuniform distribution as the upper part of Figure 1(b) shows. Using the signal processing technique such as system filters [14] or STAP [15] (the matrices are identical indeed), the data can be reconstructed as uniformly sampled. The lower part of Figure 2(b) is the effective phase-centers after reconstruction where the dashed line denotes the zero delay. It is interesting that the time delays are different for each [f.sub.c,m] because (4) associates with m.

From Figure 1(b), it can be seen that the effective phasecenters of different [f.sub.c,m] are overlapped after reconstruction. If the pulse bandwidth of all transmit channels is identical and [f.sub.c,m] steps following

[f.sub.c,m] - [f.sub.c,m-1] [less than or equal to] B, (7)

where B is the chirp bandwidth, then a chirp signal with large bandwidth can be obtained using the subband synthesis technique with the overlapped samples [16]. Subband synthesis can be realized either in the range time domain or in the range frequency domain. Besides, synthesis can be before range compression or after compression.

In summary, the MIMO-SAR processing steps followed by this paper are

(1) samples reconstruction in azimuth;

(2) imaging each subband;

(3) subband synthesis in range.

3. Phase Error and the Compensation

Channel imbalance is inevitable in multichannel system, not expecting MIMO-SAR. It presents in echo data in terms of amplitude error and phase error. We only focus on phase error in this paper. Denote the transmission characteristic of each transmit channel as [DELTA][[psi].sub.1], [DELTA][[psi].sub.2], ... , [DELTA][[psi].sub.M] and each receive channel as [DELTA][[phi].sub.1,1], [DELTA][[phi].sub.1,2], ... , [DELTA][[phi].sub.1,M], [DELTA][[phi].sub.2,1], [DELTA][[phi].sub.2,2] , ... , [DELTA][[phi].sub.2,M], ... , [DELTA][[phi].sub.M,1], [DELTA][[phi].sub.M,2], ... , [DELTA][[phi].sub.M,M], respectively. Here, the substrip of [DELTA][psi] indicates the transmit channel number, and [DELTA][[phi].sub.n,m] indicates the characteristic of channel n receiving the signals with center frequency [f.sub.c,m]. The actual signals received are

[s.sub.r,n,m]([tau], [eta]) = [DELTA][[psi].sub.m]x [DELTA][[phi].sub.n,m]x [s.sub.r,1,m] ([tau], [eta] +(n - 1) x d/2[V.sub.r]). (8)

3.1. The Receive Phase Errors. Consider the signals with center frequency [f.sub.c,m] and suppose channel 1 is the reference channel. Define the phase error

[DELTA][[epsilon].sub.r,n,m] =[DELTA][[phi].sub.n,m]/[DELTA][[phi].sub.1,m], n = 2, ... , M (9)

the receive phase error (RXE). RXE is range invariant and if arranged along azimuth, it periodically changes over a cycle of M samples. It can be treated as some form of phase error introduced by motion errors in airborne SAR system; therefore RXE brings blurring on subband imaging that the images are ambiguous and defocused in the azimuth.

The compensation of RXE can be treated as a special case of autofocus. There are numerous autofocus methods; however not all are suitable for RXE removing. Because the azimuth reconstruction will disturb the phase error, the phase estimated by autofocus is usually not the original one. We resort to the method using cost function [17] in which the phase error is trial rather than directly derived from the imagery.

Generally, the value of cost function relates to the extent of an image into focus. The broaden mainlobe of a point target and the ghosts along azimuth are the two main characteristics blurred images caused by RXE exhibited [18]. Although defocused, the ghosts share the same shape with the actual targets. According to these two characteristics, we develop a new cost function called image self-correlation in azimuth (ISCA):

[CF.sub.ISCA] (k) = [E.sub.m]{I (p, q) x I*(p, q + r)}/[E.sub.m]{I (p, q) x I*(p, q)}, (10)

where I(p, q) is a grey-scale SAR image, p and q are the range coordinate and the azimuth one, respectively, r is the azimuth offset, * denotes conjugation, and [E.sub.m] implies average in range. ISCA is normalized by the denominator in which the maximum value is 1 which happens only when k = 0.

ISCA describes the RXE in two ways. First, the ISCA of a focused image descends faster around k = 0 (zero peak) than that of a defocused one; that is, the zero peak is sharper, because the mainlobe is wider on the defocused image. Second, besides the zero peak, the curve of self-correlation may have other couples of peaks (nonzero peak). When the image is fully focused, the ghosts are suppressed that these nonzero peaks become flat.

Thus, the estimation of RXE can be described as an optimization procedure

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)

where [DELTA][[??].sub.r,m] is the optimal solution, [DELTA][[epsilon].sub.r,m] = [[DELTA][[epsilon].sub.r,2,m], [DELTA][[epsilon].sub.r,3,m], ... , [DELTA][[epsilon].sub.r,M,m], and operator ext denotes judging via the width of the zero peak or with the height of nonzero peaks of ISCA.

As the optimization in (11) has no closed-form solution, an iterative method called coordinate descent algorithm [19] is adopted in which the elements of [DELTA][[epsilon].sub.r,m] are obtained in sequence. The iterative processes as well as the computational burden are distinctively different according to the imaging algorithm applied.

The iteration flow, when frequency domain imaging algorithm such as the familiar chirp scaling (CSA) [20] is chosen, is shown on the left of Figure 2. Because some time-consuming progresses are involved in iteration, the time consumption of the whole iteration will be quite large especially when the number of channels is big. This is the inherent limitation of cost function autofocusing.

Situations will turn better when time domain imaging algorithm such as back projection imaging algorithm (BPA) [21] is used. Integrating along the target track precisely, BPA can be used in almost all SAR working modes. Resembling the idea of fast factorized BP (FFBP), the data of each single channel are imaged first and then integrated as the right of Figure 2 shows. As only multichannel integration is involved in iteration, less time is required. After each single channel is imaged for the first time, we can select some range bins with highest contrast as a new data set for the iteration ahead. This will speed up the iteration further.

3.2. The Transmit Phase Errors. When RXE is compensated, we obtain M frame focused images with different center frequencies and (8) becomes

[s'.sub.r,n,m]([tau], [eta]) = [DELTA][[psi].sub.m]x [DELTA][[phi].sub.1,m]x [s.sub.r,1,m]([tau], [eta] +(n - 1) x d/2[V.sub.r]). (12)

Similarly, suppose channel 1 is the reference channel. Define the phase error

[DELTA][[epsilon].sub.t,m] = [DELTA][[psi].sub.m][DELTA][[phi].sub.1,m]/[DELTA][[psi].sub.1][DELTA][[phi].sub.1,m], m = 2, ... , M (13)

the transmit phase error (TXE). TXE will reduce the performance of subband synthesis in which the peak sidelobe level ratio (PSLR) and the integrated sidelobe level ratio (ISLR) are worsened [22].

On the focused images, the target is compressed as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (14)

where A is a irrespective coefficient, [p.sub.r] and [p.sub.a] are the amplitude of impulse response; that is, the sinc function, [f.sub.dc] is the Doppler center. [R.sub.n,m] ([eta]) is knownas (4); so the estimation of TXE is simple compared to RXE in which

[DELTA][[epsilon].sub.t,m] = [s.sub.c,n,m]([tau], [eta])/[s.sub.c,n,1]([tau], [eta]) x exp {j 4[pi] ([R.sub.n,m]([eta]) [f.sub.c,m]- [R.sub.n,1]([eta]) [f.sub.c,1])/c} . (15)

To make the estimation more precisely and fast, not the whole image but parts of the image with high signal to cluster ratio (SCR) are used. A rectangular sliding window of appropriate size is traveled throughout the image, with the variance of the pixels in the window calculated, those high SCR targets are likely to exist in the windows with high variance [23].

4. Simulation

In this Section, point target simulations are introduced to analyze the proposed algorithm. The simulation parameters are listed in Table 1. The system has 6 channels in azimuth with center frequency ranged from 9.35GHz to 9.85GHz stepped by 100MHz, in which the overlapped ratio of subband in range frequency is zero. The full antenna length is 6min azimuth which means each aperture can share one out of six.

First, impacts of the phase error on azimuth reconstruction as well as range matched filtering are confirmed. Figure 3 gives the results of azimuth reconstruction. The RXEs added are 0.030[pi], -0.055[pi], 0.041[pi], -0.035[pi], and -0.018[pi], respectively. It is clear that the energy is confused across the spectrum as the existence of RXE. And the most important influence is that the phase no longer holds the form of quadratic which means the azimuth compression will fail. Figure 4 gives the results of subband synthesis. The TXEs added are 0.50[pi], -0.55[pi], -0.21[pi], 0.09[pi], and -0.47[pi], respectively. As mentioned in Section 3.2, the sidelobe indicators are worsened that the shape is distorted and no longer symmetrical, whereas the mainlobe remains the same; that is, the range resolution is not declined. Comparing the order of RXE with that of TXE, it can be found that azimuth reconstruction is more sensitive to the phase error.

Second, the validation of the proposed algorithm is verified. The system SNR is set to 20 dB. As only point target is involved in the simulations, the computation cost of frequency domain imaging algorithm is acceptable. Thus the chirp scaling algorithm is used. Figure 5 gives the imaging results. The ghost targets are clearly seen on Figure 5(a) and Figure 6 in which the real target is defocused without doubt. After the compensation, the ambiguity is suppressed; see Figure 5(b). The estimation results are given in Table 2.

At last, the performance of the proposed algorithm is investigated using Monte Carlo experiment. The TXE interval as well as that of RXE is set to [-0.5[pi], 0.5[pi]] with uniform distribution. The system SNR is set from 0 dB to 30 dB with step of 2 dB. For each SNR, the experiment is carried out 100 times. The experiment results are given by Figure 7, where the averaged root-mean-square error (ARMSE) is defined by

ARMSE = 1/M-1 [M.summation over (m=2)] (1/N [square root of ([N.summation over (n=1)][([[??].sub.m,n]- [[epsilon].sub.m]).sup.2])]), (16)

where M is the total channels, N is the simulation times, [epsilon] is the error set and [??] is the error estimated, and substrip m and n indicate the channel and the sample, respectively. From Figure 6, it is shown that ARMSE of RXE estimation is almost independent on system SNR. This is because the cost function is obtained from SAR image where the energy is focused, and noise is greatly suppressed by the defined cost function ISCA (10). However, ARMSE of TXE estimation is dependent on system SNR which means that (16) has room for improvement. The ARMSE order of TXE is bigger than that of RXE because the estimation of RXE is prior to that of TXE; that is, the estimation error propagates.

5. Conclusion

MIMO technique offers an opportunity to map wider image swaths with improved spatial resolution in SAR system. Aiming at the inevitable channel imbalance of multichannel system, this paper has described a new algorithm for compensating the phase errors in MIMO-SAR systems based on raw data. Compared with the internal calibration, not only the transmitter and the receiver phase imbalances but also the antenna phase imbalances can be removed. The phase error is classified as transmit phase error or receive phase error which is compensated separately. Simulations confirm the contribution of different kinds of error to image quality decline. Meanwhile the performance with different SNR is also investigated. Algorithm considering gain errors will be our future research focus.

Conflict of Interests

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

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

References

[1] J. C. Curlander and R. N. McDonough, Synthetic Aperture Radar-Systems and Signal Processing, John Wiley & Sons, Hoboken, NJ, USA, 1991.

[2] J. Klare, A. Brenner, and J. Ender, "A new airborne radar for 3D imaging-Image formation using the ARTINO principle," in Proceedings of the 6th European Conference on Synthetic Aperture Radar (EUSAR '06), pp. 16-18, 2006.

[3] A. Currie and M. A. Brown, "Wide-swath SAR," IEE Proceedings F: Radar and Signal Processing, vol. 139, no. 2, pp. 122-135, 1992.

[4] J.-H. Kim, M. Younis, P. Prats-Iraola,M. Gabele, and G. Krieger, "First spaceborne demonstration of digital beamforming for azimuth ambiguity suppression," IEEE Transactions on Geoscience and Remote Sensing, vol. 51, no. 1, pp. 579-590, 2013.

[5] W. G. Carrara, R. S. Goodman, and R. M. Majewski, Spotlight Synthetic Aperture Radar, ArtechHouse, Norwood, Mass,USA, 1995.

[6] J. Mittermayer, A. Moreira, and O. Loffeld, "Spotlight SAR data processing using the frequency scaling algorithm," IEEE Transactions on Geoscience and Remote Sensing, vol. 37, no. 5, pp. 2198-2214, 1999.

[7] R. K.Moore, J. P. Claassen, and Y.H. Lin, "Scanning spaceborne synthetic aperture radar with integrated radiometer," IEEE Transactions on Aerospace and Electronic Systems, vol. 17, no. 3, pp. 410-421, 1981.

[8] P. Prats, R. Scheiber, J. Mittermayer, A. Meta, and A. Moreira, "Processing of sliding spotlight and TOPS SAR data using baseband azimuth scaling," IEEE Transactions on Geoscience and Remote Sensing, vol. 48, no. 2, pp. 770-780, 2010.

[9] W. Xu, P. Huang, R. Wang, and Y. Deng, "Processing of multichannel sliding spotlight and TOPS synthetic aperture radar data," IEEE Transactions on Geoscience and Remote Sensing, vol. 51, no. 8, pp. 4417-4429, 2013.

[10] W.-Q.Wang, "Space-time coding MIMO-OFDMSAR for high-resolution imaging," IEEE Transactions on Geoscience and Remote Sensing, vol. 49, no. 8, pp. 3094-3104, 2011.

[11] W.-Q. Wang, "Virtual antenna array analysis for MIMO synthetic aperture radars," International Journal of Antennas and Propagation, vol. 2012, Article ID 587276, 10 pages, 2012.

[12] S. Yuepeng and Y. Ruliang, "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.

[13] J.H. G. Ender and J. Klare, "System architectures and algorithms for radar imaging by MIMO-SAR," in Proceedings of the IEEE Radar Conference, Pasadena, Calif, USA, May 2009.

[14] G. Krieger, N. Gebert, and A. Moreira, "Unambiguous SAR signal reconstruction from nonuniform displaced phase center sampling," IEEE Geoscience and Remote Sensing Letters, vol. 1, no. 4, pp. 260-264, 2004.

[15] W. Jing, M. Xing, C.-W. Qiu, Z. Bao, and T.-S. Yeo, "Unambiguous reconstruction and high-resolution imaging for multiplechannel SAR and Airborne experiment results," IEEE Geoscience and Remote Sensing Letters, vol. 6, no. 1, pp. 102-106, 2009.

[16] R. T. Lord and M. R. Inggs, "High resolution SAR processing using stepped-frequencies," in Proceedings of the IEEE International Geoscience and Remote Sensing Symposium(IGARSS '97), pp. 490-492, August 1997.

[17] F. Berizzi and G. Corsim, "Autofocusing of inverse synthetic aperture radar images using contrast optimization," IEEE Transactions on Aerospace and Electronic Systems, vol. 32, no. 3, pp. 1185-1191, 1996.

[18] C. G. Gao, Study on new techniques for high-resolution and wide-swath synthetic aperture radar [Ph.D. thesis], University of Chinese Academy of Sciences, 2012.

[19] J. K. Thomas, "Monotonic iterative algorithm for minimum-entropy autofocus," in Proceedings of the Adaptive Sensor Array Processing Workshop (ASAP '06), June 2006.

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

[21] M. D. Desai and W. K. Jenkins, "Convolution backprojection image reconstruction for spotlight mode synthetic aperture radar," IEEE Transactions on Image Processing, vol. 1, no. 4, pp. 505-517, 1992.

[22] Q. Chen, Study on new techniques for high-resolution and wide-swath spaceborne synthetic aperture radar [Ph.D. thesis], University of Chinese Academy of Sciences, 2013.

[23] R. C. Gonzalez and R. E. Woods, Digital Image Processing, Prentice Hall, 3rd edition, 2007.

Lei Zhang, (1,2) Yunkai Deng, (1) and Robert Wang (1)

(1) Department of Spaceborne Microwave Remote Sensing System, Institute of Electronics, Chinese Academy of Sciences, Beijing 100190, China

(2) University of Chinese Academy of Sciences, Chinese Academy of Sciences, Beijing 100190, China

Correspondence should be addressed to Lei Zhang; 314forever@163.com

Received 27 February 2014; Accepted 2 June 2014; Published 22 June 2014

Academic Editor:Wen-QinWang

Table 1: Simulation parameters. Parameter Value Channels 6 Center frequency 9.35 GHz~9.85 GHz Chirp bandwidth 100 MHz Pulse length 20 us Platform velocity 7000 m/s PRF 900 Hz Antenna length (full) 6m Range 800 km Table 2: The phase errors added and estimated. TXE Channel 2 3 4 5 Error 0.50[pi] -0.55[pi] -0.21[pi] 0.09[pi] Estimation 0.512[pi] -0.601[pi] -0.213[pi] 0.078[pi] TXE RXE Channel 6 2 3 4 Error -0.47[pi] 0.030[pi] -0.055[pi] 0.041[pi] Estimation -0.455[pi] 0.0321[pi] -0.0564[pi] 0.0440[pi] RXE Channel 5 6 Error -0.035[pi] -0.018[pi] Estimation -0.0372[pi] -0.0203[pi]

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Zhang, Lei; Deng, Yunkai; Wang, Robert |

Publication: | International Journal of Antennas and Propagation |

Date: | Jan 1, 2014 |

Words: | 3972 |

Previous Article: | An analytical study of the impact of arbitrary receive antenna correlation in TAS/MRC. |

Next Article: | Miniaturized antenna solution based on lossy planar microresonators to conjointly control radiation and selectivity. |

Topics: |