Open AccessArticle

Channel Phase Calibration for High-Resolution and Wide-Swath SAR Imaging with Doppler Spectrum Sharpness Optimization

School of Electronic and Communication Engineering, Guangzhou University, Guangzhou 510006, China

Beijing Institute of Tracking Telementry & Telecommunication, Beijing 100094, China

Beijing Institute of Remote Sensing Satellite, Beijing 100854, China

Author to whom correspondence should be addressed.

Sensors 2022, 22(5), 1781; https://doi.org/10.3390/s22051781

Submission received: 21 December 2021 / Revised: 30 January 2022 / Accepted: 8 February 2022 / Published: 24 February 2022

(This article belongs to the Special Issue Signal Processing in Radar Systems)

Download

Browse Figures

Versions Notes

Abstract

Channel phase calibration is a crucial issue in high resolution and wide swath (HRWS) imagery with azimuth multi-channel synthetic aperture radar (SAR) systems. Precise phase calibration is definitely required in reconstructing the full Doppler spectrum for precise HRWS imagery without high-level ambiguities. In this paper, we propose a novel calibration for HRWS SAR imagery by optimizing the reconstructed unambiguous Doppler spectrum. The sharpness of the reconstructed Doppler spectrum is applied as the metric to measure the unambiguity quality, which is maximized to retrieve the element phase error caused by channel imbalance. Real data experiments demonstrate the performance of the proposed calibration for ambiguity suppression in HRWS SAR imagery.

Keywords:

synthetic aperture radar (SAR); high resolution and wide swath (HRSW); phase calibration; maximum sharpness optimization

1. Introduction

Synthetic aperture radar (SAR) techniques require generating imagery with simultaneous high resolution and wide swath (HRWS) capability. However, in order to receive unambiguous echoes in range or/and Doppler, conventional SAR observed swath should be restricted by the fundamental trade-off between Doppler ambiguity and spatial resolution for a single-channel radar system [1]. With the assistance of an azimuth multi-channel system configuration, this basic restriction can be overcome by digital beamforming (DBF) techniques [2,3,4] on multiple receiver channel data. The received signal of each channel is simultaneously recorded and digitalized. By a Doppler spectrum reconstruction with DBF processing, the unambiguous and full Doppler spectrum can be retrieved to generate the HRWS SAR imagery, which systematically breaks the fundamental trade-off of azimuth resolution and range swath in a single channel SAR system. In general, these multi-channel HRWS SAR systems [5,6,7] can be sorted into multiple antenna systems and distributed systems. Several novel DBF algorithms have also been proposed for resolving Doppler ambiguity in the spatial domain [8,9,10,11,12,13,14,15] for these systems. Azimuth multi-channel SAR with DBF is naturally an effective approach to obtain unambiguous HRWS imagery in practice. However, unavoidable channel errors would significantly degrade the performance of ambiguity suppression, especially for an observing scene with strong targets, such as a strong ship under ocean circumstances. Channel mismatch calibration is crucial for HRWS imagery in azimuth multi-channel SAR systems.

Due to unexpected systematic and circumstance factors, such as mismatched transmitting and receiving channel characteristics. Channel imbalance brings additional phase errors between multiple channels; these errors degrade the DBF Doppler spectrum reconstruction performance dramatically. Before applying the receiving beamforming technique for Doppler spectrum reconstruction, the multi-channel SAR system should be well calibrated by compensating for the phase and amplitude imbalance between channels. In particular, the phase errors cause unexpected spectrum folding and mismatch in the DBF process, leaving serious ambiguities in the SAR imagery. Their presence will deteriorate the azimuth ambiguity-to-signal ratio (AASR) [16] of the HRWS imagery. In general, the imbalance phase errors are expected to be well compensated by calibration techniques. Recently, a series of innovative calibration algorithms [17,18,19,20,21,22,23] have been developed, which are roughly classified into two sorts, including inner calibration and data-driven calibration. In the inner calibration methods, the azimuth channel imbalance is expected to be compensated by an extra antenna measurement system [18,23], which would limit its real applications due to increased system complexity. On the other hand, data-driven calibration algorithms try to estimate the channel imbalance from data itself in an adaptive manner without any increase in system complexity. In this paper, we focus only on data-driven calibration for HRWS SAR imagery. Among the data-driven calibration algorithms, subspace projection algorithms [16,17,20,21] have been well investigated in many works. They generally estimate the channel phase errors for multi-channel SAR HRWS imagery with the same basic principle. The covariance matrix in the range-Doppler domain was expected to be estimated from the data itself. Then, eigenvalue decomposition on the covariance matrix was followed to separate the signal subspace and noise subspace. The signal subspace was spanned by the signal steering vectors with phase errors, which are orthogonal to the noise subspace. By minimizing the projection of the estimated steering vectors on the noise subspace, the imbalance phase errors are obtained. Typically, in reference [16], an auto-calibration algorithm developed on the virtual-source calibration was presented to estimate the channel errors for the distributed multi-receive radar configuration. The multi-channel data in the azimuth time domain is transferred into the Doppler domain via Fourier transform, and signal elements corresponding to different directions within each Doppler bin are regarded as the correction virtual sources. In subspace projection kind algorithms, by an iterative optimization projection into the noise subspaces, the element phase errors are expected to be obtained. Subspace projection-based channel calibrations are very popular in HRWS imaging with multi-channel SAR systems. In these works, the imbalance phase errors are usually assumed to be spatial-invariant. However, the extension to the azimuth- and range-dependent imbalance calibration is not hard by introducing the local invariant assumption and using local covariance matrix decomposition processing in the range-Doppler domain. It is necessary to emphasize that, in the subspace channel calibrations, an additional spatial degree of freedom is required, which means the channel number should be larger than that of the ambiguity. This requirement constrains both the azimuth multi-channel system design and HRWS maximum performance, which motivates our work in this paper.

In this paper, we developed a sharpness optimization-based phase calibration for HRWS SAR imagery with azimuth multi-channel configuration. In the algorithm, the imbalance phase error is assumed to be spatial-invariant. Based on the fact that the reconstruction of the full Doppler spectrum depends on the channel phase errors, the signal energy of the reconstructed Doppler spectrum will not separate into different ambiguous bins presenting sharp formation in the full Doppler domain. In the azimuth phase calibration, the estimation of the channel phase errors turns into maximizing the sharpness of the reconstruction full Doppler spectrum by correcting phase errors within ambiguity resolving, which simplifies the problem of calibration due to no need for additional spatial freedom degrees. In contrast to the former approaches, this method proposed in this paper calibrates phase imbalance in a robust manner by using all range bins as calibration sources. In addition, we can obtain the global optimal solution efficiently by accelerating the optimization with some useful skills. Detailed analysis of its performance is presented by utilizing a real SAR dataset.

This paper is organized as follows. In Section 2, the Doppler ambiguity signal model and the inverse filtering for ambiguity removal are constructed. Section 3 extends to the proposed channel phase calibration. Section 4 evaluates the comparison performance of the approach. Finally, some conclusions are drawn.

2. HRWS Imaging with Azimuth Multi-Channel SAR Configuration

In this section, we introduce the azimuth multi-channel SAR for HRWS imagery generation and the Doppler ambiguity resolving with inverse-filtering processing.

2.1. Azimuth Multi-Channel SAR Mathematical Model

The azimuth multi-channel SAR geometry is shown in Figure 1. The reference equivalent phase center (EPC) is at the origin of the Cartesian coordinate. The platform moves along the x-axis at velocity

v_{s}

. Usually, the system uses a reference antenna to transmit a wide-beam signal, and all antennas receive returns simultaneously. The data can be converted into the equivalent self-transmitting and self-receiving data by compensating a constant phase with respect to the reference channel [9]. In the ideal model, the channels are supposed to be calibrated and the baseline between equivalent phase centers is constant during the coherent processing interval (CPI). Assume that m = 1, 2, …, M and M is the number of equivalent phase centers. The

m t h

equivalent phase center is

(x_{m}, 0, 0)

at the time

t = 0

and

(x_{m} + v_{s} t, 0, H)

at time

t

. Referring to [17], the echo received at the

m t h

EPC at time

t

can be written as follows:

\begin{matrix} s_{c, m} (τ, t) = & \iint σ (x, y) \cdot h [τ - \frac{2 r_{c, m} (x, y, t)}{c}] \cdot g (t - \frac{x - x_{m}}{v_{s}}) \\ \cdot \exp [- j 4 π \frac{r_{c, m} (x, y, t)}{λ}] d x d y \end{matrix}

(1)

where

R (x, y, t) = \sqrt{{(x - x_{m} - v_{s} t)}^{2} + R_{b}^{2}}

(2)

h (τ) = a (τ) \cdot \exp (- j π α τ^{2})

(3)

and

α

is the chirp rate;

t

and

τ

denote the azimuth slow-time and range fast-time respectively;

R_{b}^{2} = \sqrt{y^{2} + H^{2}}

;

c

is the propagation velocity;

σ (x, y)

is the clutter complex reflectivity of scatter

(x, y)

on ground;

g (t)

represents the antenna pattern and other slow time-variant characters.

After transforming (1) into the two-dimension (2D) frequency domain via a 2D Fourier transform (FT), we can obtain the following equation by using the concept of instantaneous wavenumber:

S_{c, m} (K_{r}, f_{d}) = s (K_{r}, f_{d}) \cdot \exp (- j 2 π \frac{x_{m}}{v_{s}} f_{d})

(4)

where

\begin{array}{l} s (K_{r}, f_{d}) = \iint σ (x, y) \cdot H (K_{r}) \cdot G (f_{d}) \\ \cdot \exp (- j R_{b} \sqrt{K_{r}^{2} - {(2 π f_{d})}^{2} / v_{s}^{2}}) \cdot \exp (- j 2 π f_{d} \frac{x}{v_{s}}) d x d y \end{array}

(5)

and

H (K_{r}) = F T [h (τ)]

. As a result, there is only one linear phase difference

\exp (- j 2 π \frac{x_{m}}{v_{s}} f_{d})

between the

m t h

and the reference EPC for the same Doppler bin output. In the Doppler domain, for a single Doppler bin containing all clutter echoes from the same azimuth angle

θ

. The relation between the Doppler frequency

f_{d}

and the corresponding azimuth angle

θ

is given by

f_{d} = \frac{2 v_{s} \sin θ}{λ}

(6)

This linear relationship between

f_{d}

and

\sin θ

can be represented in the spatial-time spectrum, as Figure 2 shows.

In the HRWS SAR imaging, the pulse repetition frequency (PRF)

f_{r}

is usually much lower than the instantaneous Doppler bandwidth. Thus, the Doppler spectrum is ambiguous, as shown in Figure 2b, where the dots represent signal components from different directions within the identical Doppler bin. We rewrite (4) as follows:

S_{c, m} (K_{r}, f_{d}) = \sum_{i = - N + 1}^{N} s (K_{r}, f_{d} + i \cdot f_{r}) \cdot \exp [- j 2 π \frac{x_{m}}{v_{s}} (f_{d} + i \cdot f_{r})]

(7)

where the total summation term

2 N

representing the Doppler ambiguity number, should not be more than the number of EPC, which is a basic constraint for all Doppler resolving algorithms, and

- f_{r} / 2 \leq f_{d} \leq + f_{r} / 2

. For expression simplicity, we define the ambiguity number as an even number, which does not affect handling ambiguity resolving for odd cases. Therefore, the array steering matrix for the Doppler bin

f_{d}

is given by

A (f_{d}) = [\begin{matrix} a_{- N + 1}, & \cdot \cdot \cdot & a_{i}, & \cdot \cdot \cdot & a_{N} \end{matrix}]

(8)

The array steering vector for the Doppler component

f_{d} + i \cdot f_{r}

a_{i} = {[\begin{matrix} 1, & \cdot \cdot \cdot & \exp (- j 2 π \frac{x_{m}}{v_{s}} (f_{d} + i \cdot f_{r})), & \cdot \cdot \cdot & \exp (- j 2 π \frac{x_{M}}{v_{s}} (f_{d} + i \cdot f_{r})) \end{matrix}]}^{T}

(9)

where

{[\cdot]}^{T}

denotes the vector transpose operator. The array output of the Doppler bin can be written as follows:

\begin{matrix} S_{c} (f_{d}) = {[\begin{matrix} S_{c, 1} (K_{r}, f_{d}), & \cdot \cdot \cdot & S_{c, m} (K_{r}, f_{d}), & \cdot \cdot \cdot & S_{c, M} (K_{r}, f_{d}) \end{matrix}]}^{T} \\ S_{c} (f_{d}) = A (f_{d}) s (f_{d}) + e (f_{d}) \end{matrix}

(10)

where

s (f_{d}) = {[\begin{matrix} s (K_{r}, f_{d} + (- N + 1) \cdot f_{r}), & \dots & s (K_{r}, f_{d} + i \cdot f_{r}), & \dots & s (K_{r}, f_{d} + N \cdot f_{r}) \end{matrix}]}^{T}

is the signal vector, and

e (f_{d})

denotes the additive noise vector, which is assumed to be independent and white noise. The problem now is to unfold the spatial-time spectrum and reconstruct the full Doppler bandwidth spectrum.

2.2. Unambiguous Doppler Spectrum Reconstruction with Inverse-Filtering

In this subsection, we introduce non-adaptive beamforming for Doppler ambiguity resolving [10], which is seamless with the following channel calibration we proposed in this paper. The Doppler spectrum of (7) is divided into many short bins, and each bin corresponds to

2 N

Doppler components from corresponding azimuth angles. The non-adaptive beamforming approach first extracts these spectrum components from each Doppler bin by means of spatial filtering. By rearranging the extracted Doppler components in order, the full Doppler spectrum can be reconstructed. Considering no phase error involved, the spatial filtering process for Doppler ambiguity resolving is given by the following equations [16].

\begin{array}{l} \hat{s} (f_{d}) = W^{H} (f_{d}) \cdot S_{c} (f_{d}) \\ \approx W^{H} (f_{d}) A (f_{d}) s (f_{d}) \end{array}

(11)

where

\hat{s} (f_{d})

is the reconstructed Doppler spectrum,

{[\cdot]}^{H}

denotes the conjugation transpose operator, and the spatial filtering for the

2 N

Doppler components is given by

\begin{array}{l} W^{H} (f_{d}) = H \cdot A^{- 1} (f_{d}) \\ = {[w (f_{d} + (- N + 1) \cdot f_{r}), \dots, w (f_{d} + i \cdot f_{r}), w (f_{d} + N \cdot f_{r})]}_{M \times 2 N}^{H} \end{array}

(12)

H = diag {1, 1, \dots, 1}_{2 N \times 2 N}

(13)

It can be noted that beamforming filtering requires the inverse matrix of the steering matrix

A (f_{d})

, which means the channel number

M

should be not smaller than the number of ambiguity Doppler components

2 N

. After all the spectrum components are extracted, we rearrange them and obtain the unambiguous full spectrum. Then, conventional SAR imaging approaches can be used to focus the wide-swath image with high resolution. It is necessary to emphasize that the adaptive beamforming technique [18] usually provides better performance under non-uniform sampling than the non-adaptive beamforming Doppler ambiguity resolving in (11), in the circumstances of no imbalance phase errors involved. The adaptive beamforming filter can automatically match the interference subspaces even with the phase imbalance to generate optimal matches by using the coherence matrix decomposition. Herein, we use the non-adaptive beamforming filter to implement the phase error estimation in the following section. It is noted that the non-adaptive filter is fixed and independent of the imbalance phase errors. Non-adaptive filtering is sensitive to the phase imbalance, which can provide ideal Doppler spectrum reconstruction only when the phase imbalance is accurately corrected. This sensitivity would be helpful to develop the sharpness optimization to retrieval the channel phase errors. In real applications, after phase calibration, adaptive beamforming can be adapted to achieve optimal performance.

Now, we introduce the imbalance phase error into the signal model to pave the way for developing the following calibration optimization. As this work only focuses on the imbalance phase errors between azimuth multi-channels, by introducing the phase error matrix

Γ_{e} (ε)

, the multi-channel signal model in (10) can be rewritten by the following equation:

{\tilde{S}}_{c} (f_{d}; ε) = Γ_{e} (ε) A (f_{d}) s (f_{d}) + e (f_{d})

(14)

where the phase error matrix is a diagonal matrix represented as

\begin{array}{l} Γ_{e} (ε) = diag [\exp (j ε)] \\ = diag [\exp (j ε_{1}), \exp (j ε_{2}), \dots, \exp (j ε_{M})] \end{array}

(15)

where

Γ_{e} (ε)

denotes the imbalance phase error matrix with

ε = [ε_{i}], (i = 1, 2, \dots, M)

as the phase error vector.

Taking the imbalance phase errors into Doppler ambiguity resolving with the non-adaptive filtering in (11), the reconstructed Doppler signal depends on the phase error given by

\begin{array}{l} \hat{s} (f_{d}; ε) = W^{H} (f_{d}) \cdot {\hat{Γ}}_{e}^{- 1} \cdot {\tilde{S}}_{c} (f_{d}; ε) \\ \approx W^{H} (f_{d}) \cdot {\hat{Γ}}_{e}^{- 1} \cdot Γ_{e} \cdot A (f_{d}) s (f_{d}) \end{array}

(16)

It is notable from (14)–(16) that, in the presence of imbalance phase errors, Doppler ambiguity removal with spatial filtering would fail. The serious mismatch between the ideal and real steering vectors will induce not only signal loss but also significant ambiguity residue. In the next section, the phase calibration to estimate

Γ_{e} (ε)

will be introduced in detail.

3. Channel Phase Calibration with Maximum Sharpness Optimization

By optimizing the pattern design of the SAR system, the unambiguous Doppler spectrum has relative focused intensity around the Doppler centroid, while the intensity decreases away from the centroid. Due to the presence of channel phase errors, the distribution of the multi-channel reconstructed azimuth Doppler spectrum is unfocused due to the Doppler ambiguous components folding and smearing. Compared with the ideal multi-channel reconstructed Doppler spectrum, the reconstructed Doppler distribution contaminated by channel phase errors will be spread flat. Therefore, channel phase errors will lead to the broadening of the reconstructed Doppler spectrum after ambiguity removal filtering. This intrinsic phenomenon paves the way for estimating the phase errors by measuring the reconstructed Doppler spectrum distribution. At first, the reconstructed Doppler spectrum intensity is defined as follows.

\begin{array}{l} I (f_{d}; ε) = {| \hat{s} (f_{d}; ε) |}^{2} \\ = \hat{s} (f_{d}; ε) \cdot {\hat{s}}^{*} (f_{d}; ε) \end{array}

(17)

To sum up all reconstructed Doppler components, the sharpness function [24,25] of the reconstructed Doppler spectrum intensity is introduced as the metric of the reconstruction with the phase error as its variable.

P (ε) = \sum_{f_{d}} I {(f_{d}; ε)}^{2}

(18)

Herein, the sharpness optimization for phase calibration can be exploited by maximizing the sharpness of the reconstructed Doppler spectrum. The maximum sharpness optimization is given by

\hat{ε} = \arg \max P (ε)

(19)

Herein, we have to emphasize that some other metric function can also be introduced to represent the intensity distribution property of the reconstructed Doppler spectrum, such as entropy in [26] and contrast in [27]. The reason we select sharpness mainly lies in its convenience in deriving a gradient-based solver. Newton’s solver is convenient for obtaining an optimal estimation with high efficiency [28]. This solver is implemented in an iterative manner, and we assume in the (k + 1)th iteration, the phase error vector is updated with the gradient vector, and the Hessian matrix is given by

\begin{array}{l} \hat{ε} (k + 1) = \hat{ε} (k) + Δ ε_{k} \\ = \hat{ε} (k) + [- H^{- 1} (ε) \cdot \frac{\partial P (ε)}{\partial ε}] \end{array}

(20)

where

H^{- 1} (ε)

denotes the Hessian inverse of the optimization function with respect to the phase error vector. We also introduce a variable

{\hat{Γ}}_{e} = diag {\exp [j \hat{ε} (k)]}

to denote the estimated phase error matrix in the last iteration. In addition, the phase compensation process applies

{\hat{Γ}}_{e}^{- 1}

on the multi-channel signal via

{\hat{Γ}}_{e}^{- 1} \cdot {\tilde{S}}_{c} (f_{d}; ε)

. The first-order gradient is expressed by the following derivations:

\frac{\partial P (ε)}{\partial ε} = \sum_{k, r} 2 \cdot I (f_{d}; ε) \cdot \frac{\partial I (f_{d}; ε)}{\partial ε}

(21)

\frac{\partial I (f_{d}; ε)}{\partial ε} = \frac{\partial {| \hat{s} (f_{d}; ε) |}^{2}}{\partial ε} = 2 \cdot Re {{\hat{s}}^{*} (f_{d}; ε) \cdot \frac{\partial \hat{s} (f_{d}; ε)}{\partial ε}}

(22)

where

\frac{\partial \hat{s} (f_{d}; ε)}{\partial ε} = {[W^{H} (f_{d}) \cdot {\hat{Γ}}_{e}^{- 1} \cdot {\tilde{S}}_{c} (f_{d}; ε)]}^{*}

(23)

and the Hessian matrix of the sharpness optimization Hessian can be derived by the following expressions:

\begin{array}{l} H (ε) = \frac{\partial^{2} P (ε)}{\partial ε_{i} \partial ε_{j}} = \sum_{f_{d}} 2 \cdot I (f_{d}; ε) \frac{\partial I (f_{d}; ε)}{\partial ε_{i}} \\ = \sum_{f_{d}} 2 \cdot [\frac{\partial I}{\partial ε_{j}} \frac{\partial I}{\partial ε_{i}} + I (f_{d}; ε) \cdot \frac{\partial^{2} I}{\partial ε_{i} \partial ε_{j}}] \end{array}

(24)

\begin{array}{l} \frac{\partial^{2} I (f_{d}; ε)}{\partial ε_{i} \partial ε_{j}} = \frac{\partial^{2} {| \hat{s} (f_{d}; ε) |}^{2}}{\partial ε_{i} \partial ε_{j}} \\ = 2 \cdot Re {{(\frac{\partial \hat{s}}{\partial ε_{i}})}^{*} \frac{\partial \hat{s}}{\partial ε_{i}} + {\hat{s}}^{*} (f_{d}; ε) \frac{\partial^{2} \hat{s}}{\partial ε_{i} \partial ε_{j}}} \end{array}

(25)

where

\frac{\partial^{2} \hat{s}}{\partial ε_{i} \partial ε_{j}} = - W^{H} (f_{d}) \cdot {\hat{Γ}}_{e}^{- 1} \cdot {\tilde{S}}_{c} (f_{d}; ε)

(26)

In general, only with several iterations can the maximum of sharpness optimization be achieved. The convergence can be controlled by setting a maximum iteration number or judging the estimate difference between two sequential iterations. One can note that, in solving the optimization in (20), we need to calculate the first-order and second-order partial derivatives in each iteration. The Hessian matrix inverse determining the search step is calculated for every iteration, which involves the majority of the expensive computational load. In order to accelerate the solver to the sharpness optimization for the phase error matrix, we also give an accelerating scheme. By considering the termination condition to ensure the convergence property, it is possible to achieve a balance between the searching step accuracy and optimization efficiency. According to the principle of the BFGS (Broyden, Fletcher, Goldfard and Shanno) [29] algorithm, the Hessian matrix can be calculated approximately by an updated formula. In this manner, only the first initial direct calculation of Hessian and its inverse is calculated, and in the following iteration, Hessian can be obtained by an updating manner. This update can be given by the following equations:

\begin{array}{l} H_{k + 1} & = H_{k} - \frac{H_{k} s_{k} s_{k}^{T} H_{k}}{s_{k}^{T} H_{k} s_{k}} + \frac{y_{k} y_{k}^{T}}{y_{k}^{T} y_{k}} \\ s_{k} & = {\hat{ε}}_{k + 1} - {\hat{ε}}_{k} \\ y_{k} & = \frac{\partial P (ε; k + 1)}{\partial ε} - \frac{\partial P (ε; k)}{\partial ε} \end{array}

(27)

In terms of clarity, the complete calibration flowchart with the proposed maximum sharpness optimization is shown in Figure 3. One can note that the full processing flow contains three steps. The first pre-processing performs the range and azimuth FFTs on the raw multi-channel data to transfer it into the two-dimensional frequency domain. The second step calibrates the imbalance phase together with reconstructing the full Doppler spectrum without ambiguity. Finally, SAR imaging processing with the conventional range-Doppler algorithm [30] is accomplished on the reconstructed full Doppler spectrum to achieve HRWS imagery. Of course, one can also calibrate raw data with the obtained phase errors and apply adaptive beamforming for ambiguity resolving, which would take advantage of the adaptive beamformer to achieve improved image performance. The alternation of beamformers is not the key point of this paper, and in the following experiment, we perform experiments with the given flowchart. In the above deviations, the imbalance phase errors are assumed to be spatial-invariant. However, it should be noted that the extension to the range-variant imbalance calibration is straightforward by dividing the data into several range blocks and performing independent block calibration with the proposed algorithm to handle the range-dependence of array imbalance. The complete HRWS SAR processing flowchart is shown in Figure 3. The flowchart includes three import parts, namely data preprocessing, calibration phase optimization and SAR imaging processing. It is clear that the proposed algorithm is contained in the second part.

4. Experiments and Performance Analysis

In the following experiment, azimuth multi-channel SAR data is generated by using the real space-borne SAR of Sentinel-1 [30]. The reason for using this dataset to perform the comparative experiment is that the image contains different scenes of land with different strong buildings and open ocean areas with several strong ship targets. These ship targets provide necessary convenience for us to calculate the azimuth ambiguity suppression ratio quantitatively. The original single channel SAR data is sampled with PRF up to 1676 Hz, which is unambiguous in azimuth. In order to generate ambiguous multi-channel data, the data is re-sampled in azimuth time and down-sampled into 419 Hz to simulate a 4-channel uniform linear array data. According to the above parameters, the array baseline is about 8.427 m between two adjacent channels, and the equivalent system PRF is about 419 Hz. It should be emphasized that, in the simulated multi-channel data Doppler ambiguity number is three with the multi-channel system spatial degree of freedom is four. By the redundancy of the channel degree, we can perform the subspace projection algorithm [14] as the comparing rival for channel equalization. In order to verify the effectiveness of the proposed algorithm, imbalance phase errors between channels are randomly generated with the scope of

[- π, π]

and added into the multi-channel SAR data. Some important parameters of the synthetic multi-channel SAR data are listed in Table 1.

In the experiment, the multi-channel SAR imaging results without calibration and with phase calibration from both subspace projection algorithm and maximum Doppler spectrum sharpness algorithm are given in Figure 4a–c. As one can be seen from Figure 4a, because of the random phase errors present between channels, serious ambiguities and blurring are left in the imaging results. In the enlarged sub-images, the ambiguous components of ship targets are very obvious, which would be regarded as real ship in detection processing. On the other hand, multi-channel results from both calibrations perform well in ambiguity suppression with identical spatial filtering in the first experiment with a high SNR level. The ambiguity components present in Figure 4a are well suppressed in vision. This experiment demonstrates that when imbalance phase errors occur, calibration processing is necessary in order to achieve an optimal ambiguity suppression performance.

In order to further show the blur suppression performance of the proposed calibration algorithm, the multi-channel SAR images of two prominent targets on ocean in Figure 4 is enlarged in Figure 5a–c correspondingly. It can be seen that the ambiguity components of ship target 1 and target 2 are greatly suppressed after the calibrations by comparing with original results in Figure 5a.

In order to present the azimuth ambiguity suppression performance of different calibration schemes in the above experiment, the azimuth profile containing two ships defined as target 1 and target 2 is shown in Figure 6. We can note that the ambiguity components of both targets are greatly suppressed. The relative ambiguity amplitudes of targets 1 and 2 are suppressed by the proposed calibration from −9 dB and −27 dB down to about −35 dB and −50 dB, providing an improvement in the azimuth ambiguity suppression ratio up to about 28 dB. Similarly, successful calibration is also obtained by the subspace projection. It is worth noting that the relative ambiguity amplitude of target 1 is suppressed down to −32 dB, which is higher than that from the proposed calibration, reflecting the improved performance achieved by our work. Additionally, it is notable that the subspace projection calibration is less robust in the iteration number than the sharpness optimization. In real applications, the iteration number should be well set to ensure considerably optimal calibration.

In this section, we would like to further confirm the robustness of the proposed sharpness optimization calibration to strong noise. A channel calibration experiment is carried out by adding different noise into the multi-channel data with phase imbalance. Datasets with different signal-to-noise ratios (SNR) from −15 dB~20 dB are generated by adding complex white Gaussian noise to the original multi-channel data implemented by the AWGN function of MATLAB 2020a software. The experimental results are shown in Figure 7a–h. As you can see from the picture, in the cases of low SNR, such as SNR down to 0 dB, the AASR of the typical target 1 from subspace projection algorithm is about −15 dB. On the other hand, with an increase in SNR setting, such as SNR up to 5 dB, the calibration performance of the subspace projection algorithm is gradually guaranteed, as the AASR usually reaches up to −20 dB. When SNR = 20 dB, the AASR from the subspace algorithm is optimal, down to about −35 dB. Clearly, the noise effect on the subspace projection algorithm for calibration is demonstrated. We can conclude that a considerably high SNR is the fundamental condition for the success of the subspace projection algorithm. On the other hand, the sharpness optimization calibration algorithm is generally robust to strong noise. This is because it makes use of all signals to solve the phase error estimation. Even in the scenarios under low SNR, ambiguity suppression with Doppler spectrum sharpness optimization maintains its optimal performance. In the experiment, with varying SNRs, the AASRs obtained by the proposed algorithm remained stable and optimally maintained down to −25 dB. For clarity of the comparison, we provide the AASRs with varying SNR from −15 dB to 20 dB, from which we can note the advantage of the proposed calibration in handling strong noise over subspace calibration. To illustrate the excellent performance of the proposed algorithm, the AASRs of different algorithms with different SNRs are shown in Figure 8. In the experiment, with varying SNRs, the AASRs obtained by the proposed algorithm remained stable and optimally maintained down to −25 dB. For clarity of the comparison, we provide the AASRs with varying SNR from −15 dB to 20 dB, from which we can note the advantage of the proposed calibration in handling strong noise over subspace calibration. From these comparing results, one can conclude that the proposed algorithm has robustness in front of strong noise, which is suitable for calibrating a wide range of multi-channel SAR imagery applications.

5. Conclusions

Imbalance calibration is an important task for HRWS unambiguous SAR imaging with multi-channel SAR. In this paper, we proposed a phase calibration method with maximum sharpening of the reconstructed Doppler spectrum. The sharpness maximization optimization was developed to bridge the imbalance phase estimation and solve the optimization with a gradient-based algorithm. Without requirement on additional spatial freedom degrees, the proposed calibration provides reliable phase correction performance in both precision and efficiency aspects. Synthetic experiments are used to confirm the validity of the calibration by comparing conventional algorithms, which confirms the superiorities of the proposed calibration. We have to emphasize that, range- and azimuth-variant channel imbalance would be present in real HRWS SAR applications, which induces many difficulties for channel calibration with high efficiency and robustness. These interesting problems are still open in HRWS SAR research and left to future work.

Author Contributions

Conceptualization, M.Z.; methodology, M.Z. and Z.W.; software, S.H. and Z.Z.; validation, S.H.; formal analysis, S.H.; writing—original draft preparation, M.Z. and Z.Z.; writing—review and editing, M.Z. and Z.Z.; supervision, M.Z.; project administration, S.H.; funding acquisition, S.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Scientific Research Project of Guangzhou University (No.2900603340).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank all the anonymous reviewers and editors for their useful comments and suggestions that have permitted improvement in the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

Cantafio, L.J. Space-Based Radar Handbook; Artech House: Norwood, MA, USA, 1989. [Google Scholar]
Gebert, N.; Krieger, G.; Moreira, A. Digital beamforming on receive: Techniques and optimization strategies for high-resolution wide-swath SAR imaging. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 564–592. [Google Scholar] [CrossRef] [Green Version]
Krieger, G.; Gebert, N.; Moreira, A. Multidimensional waveform encoding: A new digital beamforming technique for synthetic aperture radar remote sensing. IEEE Trans. Geosci. Remote Sens. 2008, 146, 31–46. [Google Scholar] [CrossRef] [Green Version]
Younis, M.; Wiesbeck, W. Digital beamforming in SAR systems. IEEE Trans. Geosci. Remote Sens. 2003, 41, 1735–1739. [Google Scholar] [CrossRef]
Krieger, G.; Gebert, N.; Younis, M.; Patyuchenko, A.; Moreira, A. Advanced concepts for ultra-wide-swath SAR imaging. In Proceedings of the 7th European Conference on Synthetic Aperture Radar, Friedrichshafen, Germany, 2–5 June 2008; pp. 1–4. [Google Scholar]
Currie, A.; Brown, M.A. Wide-swath SAR. Proc. Inst. Electr. Eng. Feb. 1992, 139, 122–135. [Google Scholar] [CrossRef]
Moreira, A.; Krieger, G.; Hajnsek, I.; Papathanassiou, K.; Younis, M.; Lopez-Dekker, P.; Huber, S.; Villano, M.; Pardini, M.; Eineder, M.; et al. Tandem-L: A highly innovative bistatic SAR mission for global observation of dynamic processes on the Earth’s surface. IEEE Geosci. Remote Sens. Mag. 2015, 3, 8–23. [Google Scholar] [CrossRef]
Sikaneta, I.; Gierull, C.H.; Cerutti-Maori, D. Optimum signal processing for multichannel SAR: With application to high-resolution wide-swath imaging. IEEE Trans. Geosci. Remote Sens. 2014, 52, 6095–6109. [Google Scholar] [CrossRef]
Griffiths, H.D.; Mancini, P. Ambiguity suppression in SAR using adaptive array techniques. In Proceedings of the 11th Annual International Geoscience and Remote Sensing Symposium (IGARSS ’91), Espoo, Finland, 3–6 June 1991; pp. 1015–1018. [Google Scholar]
Krieger, G.; Gebert, N.; Moreira, A. Unambiguous SAR signal reconstruction from nonuniform displaced phase center sampling. IEEE Geosci. Remote Sens. Lett. 2004, 1, 260–264. [Google Scholar] [CrossRef] [Green Version]
Goodman, N.A.; Lin, S.C.; Rajakrishna, D.; Stiles, J.M. Processing of multiple receiver, spaceborne arrays for wide-area SAR. IEEE Trans. Geosci. Remote Sens. 2002, 40, 841–852. [Google Scholar] [CrossRef]
Li, Z.; Wang, H.; Bao, Z. Generation of wide-swath and high-resolution SAR Images from multichannel small spaceborne SAR system. IEEE Geosci. Remote Sens. Lett. 2005, 2, 82–86. [Google Scholar] [CrossRef]
Yang, T.; Li, Z.; Suo, Z.; Liu, Y.; Bao, Z. Performance analysis for multi-channel HRWS SAR systems based on STAP approach. IEEE Geosci. Remote Sens. Lett. 2013, 10, 1409–1413. [Google Scholar] [CrossRef]
Li, Z.; Wang, H.; Bao, Z.; Liao, G. Performance improvement for constellation SAR using signal processing techniques. IEEE Trans. Aerosp. Electron. Syst. 2006, 42, 436–452. [Google Scholar]
Liu, Y.; Li, Z.; Suo, Z.; Bao, Z. Azimuth resolution improvement of spaceborne SAR images with nearly non-overlapped Doppler bandwidth. In Proceedings of the IET International Radar Conference 2013, Xi’an, China, 14–16 April 2013; pp. 1–4. [Google Scholar]
Yang, T.; Li, Z.; Liu, Y.; Bao, Z. Channel error estimation methods for multichannel SAR systems in azimuth. IEEE Geosci. Remote Sens. Lett. 2013, 10, 548–552. [Google Scholar] [CrossRef]
Zhang, L.; Xing, M.D.; Qiu, C.W.; Bao, Z. Adaptive two-step calibration for high-resolution and wide-swath SAR imaging. IET Radar Sonar Navig. 2010, 4, 548–559. [Google Scholar] [CrossRef]
Laskowski, P.; Bordoni, F.; Younis, M. Error analysis and calibration techniques for multichannel SAR instruments. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Melbourne, VIC, Australia, 21–26 July 2013; pp. 4503–4506. [Google Scholar]
Kim, J.; Younis, M.; Becker, D.; Wiesbeck, W. Experimental performance analysis of digital beam forming on synthetic aperture radar. In Proceedings of the 7th European Conference on Synthetic Aperture Radar, CD-ROM, Friedrichshafen, Germany, 2–5 June 2008. [Google Scholar]
Guo, X.; Gao, Y.; Wang, K.; Liu, X. Improved Channel Error Calibration Algorithm for Azimuth Multichannel SAR Systems. IEEE Geosci. Remote Sens. Lett. 2016, 13, 1022–1026. [Google Scholar] [CrossRef]
Aifei, L.; Guisheng, L.; Lun, M.; Qing, X. An array error estimation method for constellation SAR systems. IEEE Geosci. Remote Sens. Lett. 2010, 7, 731–735. [Google Scholar]
Martino, G.D.; Iodice, A.; Riccio, D.; Ruello, G. Filtering of Azimuth ambiguity in stripmap synthetic aperture radar images. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2014, 7, 3967–3978. [Google Scholar] [CrossRef] [Green Version]
Wang, W.; Wang, R.; Deng, Y.; Xu, W.; Guo, L.; Hou, L. Azimuth ambiguity suppression with an improved reconstruction method based on antenna pattern for multichannel synthetic aperture radar systems. IET Radar Sonar Navigat. 2015, 9, 492–500. [Google Scholar] [CrossRef]
Morrison, R.L.; Do, M.N.; Munson, D.C. SAR image autofocus by sharpness optimization: A theoretical study. IEEE Trans. Image Process. 2007, 16, 2309–2321. [Google Scholar] [CrossRef]
Fienup, J.R. Synthetic-aperture radar autofocus by maximizing sharpness. Optic. Lett. 2000, 25, 221–223. [Google Scholar] [CrossRef]
Thomas, K.J.; Alaa, K.A. Monotonic iterative algorithm for minimum-entropy autofocus. In Proceedings of the Adaptive Sensor Array Processing Workshop, Lexington, MA, USA, 6–7 June 2006; pp. 645–648. [Google Scholar]
Berizzi, F.; Corsini, G. Autofocusing of inverse synthetic aperture radar images using contrast optimization. IEEE Trans. AES 1996, 32, 1185–1191. [Google Scholar] [CrossRef]
Larry, A. Minimization of functions having Lipschitz continuous first partial derivatives. Pac. J. Math 1996, 16, 1–3. [Google Scholar]
Fletcher, R. Practical Methods of Optimization, 2nd ed.; John Wiley & Sons: New York, NY, USA, 1987. [Google Scholar]
Cumming, I.G.; Wong, F.H. Digital Signal Processing of Synthetic Aperture Radar Data: Algorithms and Implementation; Artech House: Norwood, MA, USA, 2004. [Google Scholar]

Figure 1. Geometry of azimuth multi-channel SAR.

Figure 2. Spatial-time spectrum (a) without Doppler ambiguity and (b) with Doppler ambiguity.

Figure 3. The flowchart of HRWS SAR imaging with phase calibration.

Figure 4. Multi-channel SAR image results with different calibrations (horizontal azimuth, vertical range). (a) Result without calibration. (b) Result with subspace projection calibration. (c) Result with the proposed calibration.

Figure 5. Enlarged target image results with different calibrations (horizontal azimuth, vertical range). (a) Result without calibration. (b) Result with subspace projection calibration. (c) Result with the proposed calibration.

Figure 6. Azimuth profiles of ship targets from different calibrations.

Figure 7. Azimuth ambiguity suppression of different calibrations under SNRs from −15 dB to 20 dB. (a) Results under SNR = −15 dB, (b) Results under SNR = −10 dB, (c) Results under SNR = −5 dB, (d) Results under SNR = 0 dB, (e) Results under SNR = 5 dB, (f) Results under SNR = 10 dB, (g) Results under SNR = 15 dB, (h) Results under SNR = 20 dB.

Figure 8. AASRs with different SNRs.

Table 1. System parameters.

Wavelength	Pulse Width	Sample Rate	Bandwidth	Azimuth Points	Range Points	PRF	Equivalent Velocity
0.056 m	41.750 μs	32.3 MHz	30.1 MHz	2048	2048	1257 Hz	7062 m/s

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Zhang, M.; Huan, S.; Zhao, Z.; Wang, Z. Channel Phase Calibration for High-Resolution and Wide-Swath SAR Imaging with Doppler Spectrum Sharpness Optimization. Sensors 2022, 22, 1781. https://doi.org/10.3390/s22051781

AMA Style

Zhang M, Huan S, Zhao Z, Wang Z. Channel Phase Calibration for High-Resolution and Wide-Swath SAR Imaging with Doppler Spectrum Sharpness Optimization. Sensors. 2022; 22(5):1781. https://doi.org/10.3390/s22051781

Chicago/Turabian Style

Zhang, Man, Sha Huan, Zeya Zhao, and Zhibin Wang. 2022. "Channel Phase Calibration for High-Resolution and Wide-Swath SAR Imaging with Doppler Spectrum Sharpness Optimization" Sensors 22, no. 5: 1781. https://doi.org/10.3390/s22051781

APA Style

Zhang, M., Huan, S., Zhao, Z., & Wang, Z. (2022). Channel Phase Calibration for High-Resolution and Wide-Swath SAR Imaging with Doppler Spectrum Sharpness Optimization. Sensors, 22(5), 1781. https://doi.org/10.3390/s22051781

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Channel Phase Calibration for High-Resolution and Wide-Swath SAR Imaging with Doppler Spectrum Sharpness Optimization

Abstract

1. Introduction

2. HRWS Imaging with Azimuth Multi-Channel SAR Configuration

2.1. Azimuth Multi-Channel SAR Mathematical Model

2.2. Unambiguous Doppler Spectrum Reconstruction with Inverse-Filtering

3. Channel Phase Calibration with Maximum Sharpness Optimization

4. Experiments and Performance Analysis

5. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI