Open AccessArticle

Adaptive Microwave Staring Correlated Imaging for Targets Appearing in Discrete Clusters

Chao Tian

Zheng Jiang

Weidong Chen

and

Dongjin Wang

Key Laboratory of Electromagnetic Space Information, Chinese Academy of Sciences, University of Science and Technology of China, Hefei 230027, China

Author to whom correspondence should be addressed.

Sensors 2017, 17(10), 2409; https://doi.org/10.3390/s17102409

Submission received: 31 August 2017 / Revised: 17 October 2017 / Accepted: 18 October 2017 / Published: 21 October 2017

(This article belongs to the Section Physical Sensors)

Download

Browse Figures

Versions Notes

Abstract

Microwave staring correlated imaging (MSCI) can achieve ultra-high resolution in real aperture staring radar imaging using the correlated imaging process (CIP) under all-weather and all-day circumstances. The CIP must combine the received echo signal with the temporal-spatial stochastic radiation field. However, a precondition of the CIP is that the continuous imaging region must be discretized to a fine grid, and the measurement matrix should be accurately computed, which makes the imaging process highly complex when the MSCI system observes a wide area. This paper proposes an adaptive imaging approach for the targets in discrete clusters to reduce the complexity of the CIP. The approach is divided into two main stages. First, as discrete clustered targets are distributed in different range strips in the imaging region, the transmitters of the MSCI emit narrow-pulse waveforms to separate the echoes of the targets in different strips in the time domain; using spectral entropy, a modified method robust against noise is put forward to detect the echoes of the discrete clustered targets, based on which the strips with targets can be adaptively located. Second, in a strip with targets, the matched filter reconstruction algorithm is used to locate the regions with targets, and only the regions of interest are discretized to a fine grid; sparse recovery is used, and the band exclusion is used to maintain the non-correlation of the dictionary. Simulation results are presented to demonstrate that the proposed approach can accurately and adaptively locate the regions with targets and obtain high-quality reconstructed images.

Keywords:

microwave staring correlated imaging; adaptive imaging; discrete cluster; target strip; spectral entropy; band exclusion

1. Introduction

Microwave staring correlated imaging (MSCI), also called radar coincidence imaging (RCI), was recently proposed as a new microwave imaging method [1,2,3,4]. The essence of MSCI is to construct temporal-spatial stochastic radiation field in the imaging region, which is typically realized by a multi-transmitter configuration emitting independent stochastic waveforms. The stochastic waveforms can be transimitted directly from noise-generating microwave sources, or generated by the transmitters using random phase and frequency modulation [2,5,6]. In comparison with the conventional multiple-input multiple-output (MIMO) radar, which transmits orthogonal waveforms with each transmitter, MSCI focuses on stochastic radiation fields in the imaging plane. Achieving channel separation via waveform orthogonality is not required throughout the MSCI procedure, unlike MIMO radar [3]. With random radiation from multiphase centers, MSCI can increase the resolution of the targets within the same beam coverage under all-weather and all-day circumstances using the correlated imaging process (CIP), which requires the received signal to be combined with the temporal-spatial stochastic radiation field.

MSCI has attracted increasing attention in related research fields. In [7], MSCI was extended to 3D imaging and proven to be able to provide target 3D images using a single pulse without the requirement for motion parameter knowledge. Because there is generally a gain-phase error, which is difficult to have accurate knowledge of, a sparse auto-calibration method was proposed in [8] to perform gain-phase error calibration in the sparsity-driven RCI. To improve the independence degree of the radiation field, Bo Liu et al. presented a new type of ideal independent radiation field: the orthogonal radiation field for MSCI [9]. In [10], a novel adaptive clustered sparse Bayesian learning algorithm was proposed for extended target imaging. The algorithm exploits the continuity structure of the extended target. In [3], a typical MSCI system was designed, and an outfield experiment was performed to verify the effectiveness of the MSCI scheme.

Before the CIP, the MSCI must discretize the continuous target region to a fine grid and accurately compute the measurement matrix [1,2,3,4,7,8,9,10]. The MSCI system typically observes a wide area [1,2,11]; thus, using a uniform fine grid in the entire imaging region will generate an excessive number of grids, which makes the imaging process highly complex and in turn limits the real-time imaging performance. When the targets appear in discrete clusters, such as boats in a calm bay, vehicles in an open-air parking lot, the majority of the observation region does not exhibit scatterers, and, thus, grid subdivision will not benefit the imaging process. To reduce the complexity of the imaging process, the correlation method was employed to first estimate the subareas where the targets exist [11]. On this basis, this paper focuses on adaptive MSCI for targets in discrete clusters and uses narrow pulses of a stochastic signal to first locate target-containing strips in the range direction, which further reduces the complexity of the CIP.

The transmitters of MSCI typically emit pulse waveforms [2,3]. In this paper, we use narrower pulse waveforms to detect the imaging region. Discrete clustered targets are distributed in different range strips in the imaging region, and these strips are illuminated in different time periods when the pulse of the waveform is narrow. Hence, the system separates the echoes of the targets in different strips in the time domain. Here, a modified method that is robust against noise based on the spectral entropy [12] is put forward to detect the echoes of the discrete clustered targets; this method enables the strips with targets to be located in an adaptive manner.

In a strip with targets, the matched filter (MF) reconstruction algorithm [13] is used to locate the regions with targets, and only the regions of interest are discretized to a fine grid. Grid refinement generates more grids than the measurements, and the target scatterers are often sparsely distributed [8]; thus, sparse recovery is used. Because the scatterers are distributed in a continuous scene, the scattering centers are generally located off the grids, regardless of the fineness of the grid. Therefore, a gridding error is generated and is inversely proportional to the grid size. For fine grids, the increasing coherence of the dictionary makes it difficult to ensure the sparse reconstruction, which relies on a low correlation of the columns of the measurement matrix. The band exclusion (BE) [14] is used to solve this problem.

The remainder of this paper is outlined as follows. In Section 2, the MSCI mode is presented. Section 3 presents the adaptive imaging method in detail. In Section 4, the effectiveness of the proposed method is verified using numerical examples. Finally, Section 5 concludes the paper.

2. MSCI Model

MSCI can be realized using a multi-transmitter configuration to transmit time-independent and group-orthogonal waveforms [1,2,4]. As shown in Figure 1, the MSCI system is generally composed of N transmitters and one receiver, whose position vectors are denoted by

{\vec{r}}_{i}^{t}

and

{\vec{r}}^{r}

, respectively. The scatterers are distributed in the imaging region S;

{\vec{r}}_{l}

represents the position vector of the l-th scatterer. S is in the XY plane, which is H below the origin.

R_{0}

is the distance from the center of the imaging region to the origin.

θ

is the angle between antenna aperture D and the XY plane.

The correct model for MSCI is Maxwell’s equations, but the simpler scalar wave equation and first-order Born approximation are commonly used; thus, the echo

E^{s c a} (t)

is modeled as [1,2,15]

E^{s c a} (t) = \int_{S} \frac{1}{4 π |\vec{r} - {\vec{r}}^{r}|} \sum_{i} \frac{S t_{i} (t - \frac{|\vec{r} - {\vec{r}}_{i}^{t}|}{c} - \frac{|\vec{r} - {\vec{r}}^{r}|}{c})}{4 π |\vec{r} - {\vec{r}}_{i}^{t}|} σ (\vec{r}) d \vec{r},

(1)

where c and

S t_{i} (t)

are the speed of light, and excitation signal of the i-th transmitter, respectively. The target information in the imaging area is described by the scattering coefficients

σ (\vec{r})

. The modified radiation field is defined as

E_{r a d} (\vec{r}, t) = \frac{1}{4 π |\vec{r} - {\vec{r}}^{r}|} \sum_{i} \frac{S t_{i} (t - |\vec{r} - {\vec{r}}_{i}^{t}| / c - |\vec{r} - {\vec{r}}^{r}| / c)}{4 π |\vec{r} - {\vec{r}}_{i}^{t}|}

(2)

and

E^{s c a} (t)

is considered to be contaminated by noise

n (t)

; thus, the received signal is

E_{r} (t) = \int_{S} E_{r a d} (\vec{r}, t) σ (\overset{⇀}{r}) d \overset{⇀}{r} + n (t)

(3)

Let

t_{k}, k = 1, \dots, K

be the sample times. Set

E_{r} = (E_{r} (t_{k})) \in C^{K}

as the data vector. Then, the imaging region is discretized as

I_{L} = [{\vec{r}}_{1}, {\vec{r}}_{2}, \dots, {\vec{r}}_{l}, \dots, {\vec{r}}_{L}], l = 1, 2, \dots, L .

(4)

The scattering coefficient vector

σ

is expressed as

σ = {[σ_{1}, σ_{2}, \dots, σ_{l}, \dots, σ_{L}]}^{T}, l = 1, 2, \dots, L,

where the components of

σ

are equal to the scattering coefficients when the grid points are the nearest grid points to positions

\{{\vec{r}}_{l}\}

and zero otherwise. The dynamic range is introduced to describe the difference in scattering coefficients which is defined as

σ_{\max} /σ_{\min}

, where

σ_{\max} = \max_{l} |σ_{l}|, σ_{\min} = \min_{l} |σ_{l}| .

(5)

Let the measurement matrix be

E_{r a d} = [\begin{matrix} E_{r a d} ({\vec{r}}_{1}, t_{1}) & E_{r a d} ({\vec{r}}_{2}, t_{1}) & \dots & E_{r a d} ({\vec{r}}_{L}, t_{1}) \\ E_{r a d} ({\vec{r}}_{1}, t_{2}) & E_{r a d} ({\vec{r}}_{2}, t_{2}) & \dots & E_{r a d} ({\vec{r}}_{L}, t_{2}) \\ \dots & \dots & \dots & \dots \\ E_{r a d} ({\vec{r}}_{1}, t_{K}) & E_{r a d} ({\vec{r}}_{2}, t_{K}) & \dots & E_{r a d} ({\vec{r}}_{L}, t_{K}) \end{matrix}] .

(6)

Then, the mathematic model of MSCI can be described in the form

E_{r} = E_{r a d} σ + d + n,

(7)

where

n = (n (t_{k}))

is the external noise and the gridding error is given by

d = E_{r} - E_{r a d} σ - n

. To estimate the scattering coefficients, the MSCI solves the inverse problem, which recovers

σ

from the measurement vector

E_{r}

and dictionary

E_{r a d}

3. Adaptive Imaging Method for Targets Appearing in Discrete Clusters

According to Equation (7), MSCI must discretize the continuous target space to a fine grid and accurately compute the measurement matrix. The MSCI system commonly observes a wide area; thus, using a uniform fine grid in the entire imaging area will generate an excessive number of grid-cells, which will increase the computational burden. This section will introduce an adaptive imaging method in detail, which significantly reduces the computational burden when the targets appear in discrete clusters.

The transmitters of MSCI typically emit pulse waveforms [2,3], and the waveform from the i-th transmitter is

S t_{i} (t) = \sum_{q = 0}^{Q - 1} r e c t (\frac{t - q T}{T_{p}}) \cdot s t_{i q} (t),

(8)

where

T_{p}

, T and Q are the pulse width, pulse repetition period, and pulse number, respectively. In this paper, we use narrower pulse waveforms to detect the imaging region. When the pulse width

T_{p}

is narrow, the modified radiation field

E_{r a d} (\vec{r}, t)

is distributed only in a part of the imaging region at a given time, whose width

δ y

is associated with the pulse width,

δ y \approx \frac{1}{2} c T_{p} / \sin θ

Discrete clustered targets are distributed in different range strips in the imaging region, and the target strip is defined as a strip with targets, as shown in Figure 2. The spacing between two target strips is larger than

δ y

. These target strips are illuminated in different time periods when the pulse width of the waveform is narrow; thus, the system separates the echoes of the targets in the time domain. Figure 2 shows only the received signal for

q = 0

t_{s, g}^{0}, t_{e, g}^{0} (g = 1, 2, \dots, G)

are the endpoints of the echoes for the targets in different target strips. The endpoints are located using the method introduced in Section 3.1, and the estimates are

t_{s, g}, t_{e, g} (g = 1, 2, \dots, G)

. We define

E_{r, g} = {[E_{r, g}^{0}, E_{r, g}^{1}, \dots, E_{r, g}^{Q - 1}]}^{T},

(9)

where

E_{r, g}^{q} = [E_{r} (t_{s, g} + \frac{T_{p}}{2} + q T), \dots, E_{r} (t_{k} + q T), \dots, E_{r} (t_{e, g} - \frac{T_{P}}{2} + q T)]

, and

t_{k}

is the sample time that satisfies

t_{s, g} + \frac{T_{P}}{2} ⩽ t_{k} ⩽ t_{e, g} - \frac{T_{P}}{2}

. The noise component of

E_{r, g}

n_{g}

. By combining the detected endpoints

t_{s, g}, t_{e, g} (g = 1, 2, \dots, G)

and the imaging geometry, the target strips can be searched in an adaptive manner, as shown in Figure 2. In the g-th target strip, the dictionary

E_{r a d, g}

that corresponds to

E_{r, g}

can be calculated according to Equation (6).

In the same target strip, the echoes of the targets overlap in the time domain; thus, the CIP combining

E_{r a d, g}

and the corresponding received signal vector

E_{r, g}

is used to locate the regions with targets and achieve the target recovery results.

3.1. Location of the Target Strips

In the received signal

E_{r} (t)

, the echo

E^{s c a} (t)

is highly contaminated by noise

n (t)

at low signal-to-noise ratios (SNRs). However, the modified method based on spectral entropy proposed in this paper can locate the endpoints of the echoes to be free of non-echo regions using the difference between the spectrum of the echoes and that of noise.

The received signal for

q = 0

is divided into a sequence of frames with frame length

N_{0}

, and the signals in frame m are

s_{m 0}, s_{m 1}, \dots, s_{m (N_{0} - 1)}

. These signals are converted from the time domain to the frequency domain via fast Fourier transform (FFT), and

N_{F F T}

is the total number of frequency components in the FFT; thus,

S_{m} (f_{0}), S_{m} (f_{1}), \dots, S_{m} (f_{(N_{F F T} - 1)})

are obtained. We define

e s_{m} (f_{k}) = {|S_{m} (f_{k})|}^{γ} .

(10)

To improve the discriminability of the spectrum between echo and non-echo signals,

γ

should be greater than 1. Then, the probability distribution function (pdf) for the spectrum is estimated by normalizing the frequency components:

P_{m i} = e s_{m} (f_{i}) / \sum_{k = 0}^{N_{F F T} - 1} e s_{m} (f_{k}), i = 0, \dots, N_{F F T} - 1 .

(11)

The negative entropy [12]

H_{m}

of frame m is defined as follows:

H_{m} = \sum_{i = 0}^{N_{0} - 1} P_{m i} \log P_{m i} .

(12)

Here,

H_{m}

is adjusted by subtracting

C_{H}

, the average negative entropy of the first 80 frames, and then the entropy features (EFs) to locate the endpoints is obtained as follows:

E F_{m} = {(H_{m} - C_{H})}^{2} .

(13)

We define the contrast of EF as

E F_C = E F_{\max} /E, F_{m e a n},

(14)

where

E F_{\max}

is the max EF, and

E F_{m e a n}

is the average EF of the first 80 frames. A larger

E F_C

is favorable for finding the endpoints; thus, in the practical imaging process, we set

γ = \arg \max_{γ} E F_C (γ) .

(15)

The estimated endpoints

t_{s, g}, t_{e, g} (g = 1, 2, \dots, G)

can be obtained using the algorithm presented in Section 3 of Reference [16]. Then, with the endpoints and configuration parameters of the imaging system, the target strips can be adaptively obtained. The location of the target strips is further analyzed in Section 4.1.

3.2. Location of the Regions with Targets in a Target Strip

According to [2], the imaging performance of MSCI depends on the self-correlation of

E_{r a d} (\vec{r}, t)

, i.e.,

R_{I} (\vec{r}, {\vec{r}}^{'}) = \frac{|\int E_{r a d} (\vec{r}, t) E_{r a d}^{*} ({\vec{r}}^{'}, t) d t|}{\sqrt{\int E_{r a d} (\vec{r}, t) E_{r a d}^{*} (\vec{r}, t) d t} \sqrt{\int E_{r a d} ({\vec{r}}^{'}, t) E_{r a d}^{*} ({\vec{r}}^{'}, t) d t}} .

(16)

Figure 3 shows the coherence pattern

R_{I} (\vec{r}, {\vec{r}}_{0}),

where

{\vec{r}}_{0}

is the center of the imaging region. The full width at half maximum in the x- and y-directions of the coherence pattern are

Δ x_{0}

and

Δ y_{0}

, respectively. The scatterers in the main lobe of the coherence pattern cannot be easily resolved because of the high self-correlation of

E_{r a d} (\vec{r}, t)

. Here, we select

Δ x_{0}

as the initial grid spacing.

In the case of a coarse grid spacing

Δ x = Δ y = Δ x_{0}

, the g-th target strip is discretized as

I_{L, g}^{(0)} = [{\vec{r}}_{g, 1}^{(0)}, {\vec{r}}_{g, 2}^{(0)}, \dots, {\vec{r}}_{g, l}^{(0)}, \dots, {\vec{r}}_{g, L_{g, 0}}^{(0)}], l = 1, 2, \dots, L_{g, 0} .

(17)

Then, we can calculate

E_{r a d, g}^{(0)}

, the dictionary corresponding to

E_{r, g}

for the g-th (

g = 1, 2, \dots, G

) target strip, according to Equation (6). In this case, the gridding error

d^{(0)} = E_{r, g} - E_{r a d, g}^{(0)} σ_{g}^{(0)} - n_{g}

is large, so the MF reconstruction is introduced. The MF reconstruction algorithm has a low resolution but can adapt to the low-SNR and the large-gridding-error environments [13,17] to accurately locate the regions with targets. The scatter coefficients are estimated using the MF reconstruction algorithm as follows:

{\hat{σ}}_{g}^{(0)} = {(E_{r a d, g}^{(0)})}^{H} E_{r, g},

(18)

where

{(E_{r a d, g}^{(0)})}^{H}

is the conjugate transpose of the measurement matrix

E_{r a d, g}^{(0)}

We set the threshold value

ζ

. If

{\hat{σ}}_{g, l}^{(0)} > ζ {\hat{σ}}_{g, \max}^{(0)}

, there may be targets near the position

{\vec{r}}_{g, l}^{(0)}

, and denser grids will be built around it. To reduce the gridding error, we consider the fractional grid

Z / F = \{j / F : j \in Z\}

with a large integer

F \in N

, which is the refinement factor [14].

3.3. Target Reconstruction

3.3.1. Grid Refinement

We reduce the grid spacing in the regions with targets,

Δ x = Δ x / 2, Δ y = Δ y / 2

; the grid refinement method is shown in Figure 4.

Suppose that the measurement matrix

E_{r a d, g}^{(n - 1)}

has been obtained from the (n−1)-th discretization; then, the n-th discretization increases the number of grids by

L_{g, n}

whose positions are

{\vec{r}}_{g, l}^{(n)}, l = 1, \dots, L_{g, n}

. Hence,

E_{r a d, g}^{(n)} = [\begin{matrix} E_{r a d, g}^{(n - 1)} & E_{d, g}^{(n)} \end{matrix}],

(19)

where

E_{d, g}^{(n)} = [\begin{matrix} β_{1}^{(n)} & \dots & \begin{matrix} β_{l}^{(n)} & \dots \end{matrix} & β_{L_{g, n}}^{(n)} \end{matrix}]

and

β_{l}^{(n)} = {[{(β_{l}^{(n)})}^{0}, {(β_{l}^{(n)})}^{1}, \dots, {(β_{l}^{(n)})}^{Q - 1}]}^{T}

{(β_{l}^{(n)})}^{q} = [\begin{matrix} E_{r a d} ({\vec{r}}_{g, l}^{(n)}, t_{s, g} + \frac{T_{p}}{2} + q T) & \dots & \begin{matrix} E_{r a d} ({\vec{r}}_{g, l}^{(n)}, t_{k} + q T) & \dots \end{matrix} & E_{r a d} ({\vec{r}}_{g, l}^{(n)}, t_{e, g} - \frac{T_{p}}{2} + q T) \end{matrix}]

3.3.2. Target Reconstruction Method

The grid refinement generates more grids than the measurements, and the scatterers of targets are often sparsely distributed; thus, sparse recovery is used. The orthogonal matching pursuit (OMP) is a suitable choice for target reconstruction because of its low computational complexity and accurate results in practical applications [18,19]. However, the off-grid scatterer will spill non-zero values into all grid-cells; in this case of denser grids, the increasing coherence of dictionary makes it difficult to ensure the OMP reconstruction, which relies on a low correlation of the columns of the measurement matrix.

Here, we introduce BE [14]. The discrete form of Equation (16) can be denoted by the following pairwise coherence:

μ (m, l) = \frac{|〈φ_{m}, φ_{l}〉|}{∥φ_{m}∥ ∥φ_{l}∥},

(20)

where

φ_{m}

and

φ_{l}

are the m-th and l-th column of

E_{r a d, g}^{(n)}

, respectively. Let

η > 0

and define the

η - coherence

band of index l to be the set

B_{η} (l) = \{m | μ (m, l) > η\}

(21)

and the

η - coherence

band of the index set S is the set

B_{η} (S) = \cup_{l \in S} B_{η} (l)

. The double

η - coherence

band is defined as follows:

B_{η}^{(2)} (S) \equiv B_{η} (B_{η} (S)) .

(22)

When the scatterers in a cluster are sufficiently separated with respect to the size of the region represented by the double

η - coherence

band, the imaging performance will be significantly improved by inbedding BE into the OMP. Considering the double coherence band generated by grid refinement, the band-excluded orthogonal matching pursuit (BOMP) is stated in Algorithm 1. If

η = 1

is satisfactory, the BOMP is converted to the OMP; thus, the OMP is a special form of the BOMP, and the BOMP can adapt to a wider range of scenarios.

Algorithm 1: BOMP Algorithm.

Input: $E_{r a d, g}^{(n)}, E_{r, g}, η > 0$ .
Initialization: $x^{0} = 0$ , $r^{0} = E_{r, g}$ , and $P^{0} = \emptyset$ .
Iteration: For $j = 1, \dots, s$
Step 1: $i_{\max} = \arg \max_{i} |〈r^{j - 1}, φ_{i}〉|, i \notin B_{η}^{(2)} (P^{j - 1})$ ;
Step 2: $P^{j} = P^{j - 1} \cup \{i_{\max}\}$ ;
Step 3: $x^{j} = \arg \min_{z} {∥E_{r a d, g}^{(n)} z - E_{r, g}∥}_{2}$ subject to $supp (z) \in P^{j}$ ;
Step 4: $r^{j} = E_{r, g} - E_{r a d, g}^{(n)} x^{j}$ .
Output: $x^{s}$ .

To achieve the desired resolution, set

n = n + 1

and repeat from the grid refinement until

{∥E_{r a d, g}^{(n)} x^{s} - E_{r, g}∥}_{2} ⩽ ε

n = N_{\max}

; then, the reconstructed image is obtained for the g-th target strip, where

g = 1, 2, \dots, G

Remark 1.

The BOMP outperforms OMP in terms of maintaining the target profile in the recovery image. The condition for using the BOMP is that the scatterers in a cluster are sufficiently separated with respect to the size of the region represented by the double η-coherence band, as shown in Figure 5a. In Figure 5b, the scatterer

B_{2}

is in the double η-coherence band of the scatterer

B_{1}

, so the scatterer

B_{2}

will be excluded when performing Step 1 of the BOMP algorithm, which results in target reconstruction failure.

Remark 2.

For the n-th discretization, the computational cost of the OMP algorithm is approximately

O (s {K^{'}}_{g} {L^{'}}_{g, n})

[19], where

K_{g}^{'}

L_{g, n}^{'}

are the numbers of the rows and columns of

E_{r a d, g}^{(n)}

, respectively. The cost of BOMP increases to

O (s (N_{η b} + 2) {K^{'}}_{g} {L^{'}}_{g, n})

due to the introduction of BE, where

N_{η b}

is the average of the columns in the η-coherence band of the object support. After another grid subdivision,

L_{g, n + 1}^{'} = 4 L_{g, n}^{'}

N_{η b}

also becomes four times the original value, so the total cost increases to approximately 16-fold.

4. Simulation Results and Discussion

Several simulation results are presented in this section to evaluate the effectiveness of the proposed approach. An X-band MSCI radar system with a carrier frequency of

f_{c} = 9.3

GHz is considered. The transceiver configuration is shown in Figure 1. The size of D is

4 \times 1.5

m, and 36 transmitters are uniformly distributed in D with an inter-element spacing of

d 1 = d 2 = 0.5

m. The vertical distance H and angle

θ

are

600

π /6

rad, respectively.

The transmitters emit independent frequency-hopping waveforms with a bandwidth of 500 MHz:

S t_{i} (t) = \sum_{q = 0}^{Q - 1} r e c t (\frac{t - q T}{T_{p}}) \cdot \exp [j 2 π (f_{c} + f_{i q}) (t - q T) + ϕ_{i q}],

(23)

where

T_{p} = 20

ns and

T = 10 μ

Q = 200

and

f_{i q}, ϕ_{i q}

are the baseband frequency and initial phase of the q-th pulse of the i-th transmitter, respectively. The baseband frequency stochastically hops between 300 MHz and 800 MHz with q, and the initial phase is randomly set in the range of 0 to 2

π

. The received signal is mixed with zero-mean Gaussian noise, and the sampling frequency is 4 GHz. We set

ε = 10^{- 3}

N_{\max} = 4

and

η = 0.8

4.1. Location of the Target Strips

Using the target model in Figure 6, the received signal for an SNR of 5 dB can be obtained as shown in Figure 7; the echoes are highly contaminated by noise at low SNRs. However, the endpoints of the echoes can be located to be free of non-echo regions using the difference in their frequency spectra.

The received signal is divided into a sequence of frames with frame length

N_{0} = 80

, and the total number of frequency components in the FFT is

N_{F F T} = 256

. The negative spectral entropy for different

γ

is shown in Figure 8. The spectrum amplitude of the non-echo segments is flat and lower than that of the echo segments; thus, the difference in

e s_{m} (f) = {|S_{m} (f)|}^{γ}

between echo and non-echo segments increases when

γ

is larger than 1, which improves the discriminability of the negative entropy, as shown in Figure 8a,b. However, when

γ

increases further, the value of the negative spectral entropy in the non-echo regions fluctuates considerably, which is intractable for the endpoint detection, as shown in Figure 8c,d.

This paper locates the endpoints of the echoes according to the feature EF, and

E F_C

is used to measure the discriminability of EF between echo and non-echo segments. The performance curve of

E F_C

with respect to

γ

is obtained by increasing

γ

from 1 to 7, as shown in Figure 9. When

1 ⩽ γ < 2.5,

the difference in EF between echo and non-echo segments increases with increases in

γ

, and the peak value is obtained when

γ = 2.5

. However, when

γ

increases further, the negative spectral entropy of the non-echo segment fluctuates considerably, which causes a decrease in

E F_C

, consistent with the phenomenon shown in Figure 8.

Figure 10 shows the negative spectral entropy and EF when

γ = 2.5

. The negative spectral entropy of the non-echo segment is stable, and there is a clear difference in EF between the echo and non-echo segments, which is favorable for detecting the boundaries of the echo segments.

Then, the detected beginning and ending boundaries of the echo segments can be obtained using the algorithm in Section 3 in [16], which are shown in Figure 11 and indicate the highly successful endpoint detection.

The distribution of the normalized modified radiation field at a given time is calculated according to the geometry of the imaging system. Figure 12 shows the normalized modified radiation field in the imaging region at

t = t_{s, 1}

t = t_{e, 1}

t = t_{s, 2}

and

t = t_{e, 2}

. The modified radiation field begins to contact certain clustered targets when

t = t_{s, 1}

, and it leaves the targets when

t = t_{e, 1}

; thus, we can determine that Strip 1 may cover certain targets. In other words, Strip 1 is a target strip. Similarly, we know that Strip 2 is also a target strip.

4.2. Target Reconstruction

For the first target strip, the imaging result obtained using the BOMP algorithm in the coarse grids (

Δ x_{0} = 1

m) is given in Figure 13a, which fails because the gridding error is as large as the echo data. However, in Figure 13b, the MF reconstruction algorithm is robust against gridding error; thus, the target regions can be estimated in the red boxes. Because

ζ

should be less than the reciprocal of the dynamic range

σ_{\max} /σ_{\min} = 6

, we set

ζ = 0.15

when estimating the target regions in this paper. Then, a denser grid is built around the estimated locations with a finer resolution. Figure 14a,b show the poor imaging result obtained by the OMP because, in practice, the scatterers are not located precisely on the pre-discretized grids, and the mutual coherence increases with F when the nearby columns of the measurement matrix become highly correlated. As expected, better recovery results are obtained by the BOMP algorithm in Figure 14c,d, where the object locations are accurately estimated, and the estimated scattering coefficients are more accurate than the OMP result.

In the image reconstructed by the OMP in Figure 14a, false scatterers appear near the reconstructed strong scatterers. Due to the introduction of BE, the false scatterers are eliminated, and a more accurate reconstructed image can be obtained using the BOMP, as shown in Figure 14c. Similarly, in the second target strip, the target regions can be estimated in the red boxes shown in Figure 15a, and the BOMP yields a better recovery result than the OMP for the targets, as shown in Figure 15b,c.

For the target model in Figure 6, we can determine that the targets are distributed in the red boxes in Figure 13b and Figure 15a after locating the target strips and target reconstruction by the MF in the coarse grids. Then, only the regions in the red boxes, which represent just a tiny fraction of the imaging region, are discretized to a fine grid. Thus, the proposed method can significantly reduce the column number of the measurement matrix in Equation (6) and is more computationally efficient than the method using a uniform fine grid for the entire imaging region.

4.3. Performance Simulations

We randomly set 100 discrete clustered target models with a dynamic range of 6, and the spacing of the scatterers in a cluster is greater than the size of the green ellipses in Figure 14a,c. In the target models, all scatterers are distributed in Strips 1 and 2, as shown in Figure 12.

An endpoint detection is counted as a success if

0 ⩽ t_{s, g}^{0} - t_{s, g} ⩽ \frac{T_{p}}{2}

and

0 ⩽ t_{s, g}^{0} - t_{s, g} ⩽ \frac{T_{p}}{2} (g = 1, 2)

. For the 100 target models, the success rates of endpoint detection with respect to different levels of white noise are shown in Table 1. The success rates when

γ = \arg \max_{γ} E F_C (γ)

are higher than those when

γ

is fixed at 1 or 5.

To further evaluate the imaging performance, we perform grid subdivision four times. The refinement factor F increases from 2 to 16, and the performance curve of the probability of success with respect to F is obtained, as shown in Figure 16b. It can be seen that success rate rises with F, and the BOMP outperforms the OMP because of the introduction of BE. Figure 16a shows that the relative gridding error

{∥d∥}_{2} / {∥E_{r}∥}_{2}

is approximately inversely proportional to the refinement factor.

5. Conclusions

MSCI is based on the temporal-spatial stochastic radiation field and can achieve ultra-high resolution in real aperture staring radar imaging under all-weather and all-day circumstances. This paper proposes an adaptive MSCI for targets that appear in discrete clusters. The adaptive MSCI can locate the target strips in an adaptive manner and adaptively estimate the regions with targets in a target strip; thus, a uniform fine grid in the entire imaging area is avoided, and the complexity of the imaging process is reduced. Then, the applicability of MSCI is significantly improved. With several numerical simulations, the proposed approach is proven to be effective and its performance is analyzed. The adaptive imaging method for lower SNR conditions is worthy of further consideration.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (No. 61431016). The authors would like to thank the editors and reviewers for their insightful comments.

Author Contributions

Chao Tian designed the study, analyzed the results and wrote the paper; Zheng Jiang designed and performed the experiments; Weidong Chen and Dongjin Wang supervised its analysis, edited the manuscript, and provided their valuable suggestions to improve this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

Guo, Y.; He, X.; Wang, D. A novel super-resolution imaging method based on stochastic radiation radar array. Meas. Sci. Technol. 2013, 24, 074013. [Google Scholar] [CrossRef]
Li, D.; Li, X.; Qin, Y.; Cheng, Y.; Wang, H. Radar coincidence imaging: An instantaneous imaging technique with stochastic signals. IEEE Trans. Geosci. Remote Sens. 2014, 52, 2261–2277. [Google Scholar]
Cheng, Y.; Zhou, X.; Xu, X.; Qin, Y.; Wang, H. Radar coincidence imaging with stochastic frequency modulated array. IEEE J. Sel. Top. Signal Process. 2017, 11, 414–427. [Google Scholar] [CrossRef]
He, X.; Liu, B.; Chai, S.; Wang, D. A novel approach of high spatial-resolution microwave staring imaging. In Proceedings of the Asia-Pacific Conference on Synthetic Aperture Radar (APSAR), Tsukuba, Spain, 23–27 September 2013; pp. 75–78. [Google Scholar]
Axelsson, S.R.J. Noise radar using random phase and frequency modulation. IEEE Trans. Geosci. Remote Sens. 2004, 42, 2370–2384. [Google Scholar] [CrossRef]
Liu, G.; Shi, X.; Lu, J.; Yang, G.; Song, Y. Design of noise FM-CW radar and its implementation. IEE Proc. F (Radar Signal Proc.) 1991, 138, 420–426. [Google Scholar]
Li, D.; Li, X.; Cheng, Y.; Qin, Y.; Wang, H. Three dimensional radar coincidence imaging. Prog. Electromagn. Res. M 2013, 33, 223–238. [Google Scholar] [CrossRef]
Zhou, X.; Wang, H.; Cheng, Y.; Qin, Y. Sparse auto-calibration for radar coincidence imaging with gain-phase errors. Sensors 2015, 15, 27611–27624. [Google Scholar] [CrossRef] [PubMed]
Liu, B.; Wang, D. Orthogonal radiation field construction for microwave staring correlated imaging. Prog. Electromagn. Res. M 2017, 57, 139–149. [Google Scholar] [CrossRef]
Zhou, X.; Wang, H.; Cheng, Y.; Qin, Y. Radar coincidence imaging by exploiting the continuity of extended target. IET Radar Sonar Navig. 2016, 11, 60–69. [Google Scholar] [CrossRef]
Li, D.; Li, X.; Cheng, Y.; Qin, Y.; Wang, H. Radar coincidence imaging in the presence of target-motion-induced error. J. Electron. Imaging 2014, 23, 023014. [Google Scholar] [CrossRef]
Shen, J.; Hung, J.; Lee, L. Robust entropy-based endpoint detection for speech recognition in noisy environments. In Proceedings of the 5th International Conference on Spoken Language Processing, Sydney, Australia, 30 November–4 December 1998; Volume 98, pp. 232–235. [Google Scholar]
Kay, S.M. Fundamentals of Statistical Signal Processing: Detection Theory; Prentice Hall PTR: Upper Saddle River, USA, 1998; Volume 2. [Google Scholar]
Fannjiang, A.; Liao, W. Coherence pattern–guided compressive sensing with unresolved grids. SIAM J. Imaging Sci. 2012, 5, 179–202. [Google Scholar] [CrossRef]
Nolan, C.J.; Cheney, M. Synthetic aperture inversion. Inverse Probl. 2002, 18, 221. [Google Scholar] [CrossRef]
Huang, L.; Yang, C. A novel approach to robust speech endpoint detection in car environments. In Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing Proceedings, Istanbul, Turkey, 5–9 June 2000; Volume 3, pp. 1751–1754. [Google Scholar]
Yurduseven, O.; Imani, M.F.; Odabasi, H.; Gollub, J.; Lipworth, G.; Rose, A.; Smith, D.R. Resolution of the frequency diverse metamaterial aperture imager. Prog. Electromagn. Res. 2015, 150, 97–107. [Google Scholar] [CrossRef]
Tropp, J.A. Greed is good: Algorithmic results for sparse approximation. IEEE Trans. Inf. Theory 2004, 50, 2231–2242. [Google Scholar] [CrossRef]
Tropp, J.A.; Gilbert, A.C. Signal recovery from random measurements via orthogonal matching pursuit. IEEE Trans. Inf. Theory 2007, 53, 4655–4666. [Google Scholar] [CrossRef]

Figure 1. Imaging geometry of the microwave staring correlated imaging system.

Figure 2. Received signal of the clustered targets in a narrow pulse MSCI. G is the number of target strips.

Figure 3. Coherence pattern. (a) three-dimensional view; and (b) plane view.

Figure 4. Grid refinement method.

Figure 5. Coherence band and the scatterers. (a) the scatterer

A_{2}

is outside the double

η

-coherence band of

A_{1}

; (b) the scatterer

B_{2}

is in the double

η

-coherence band of

B_{1}

Figure 5. Coherence band and the scatterers. (a) the scatterer

A_{2}

is outside the double

η

-coherence band of

A_{1}

; (b) the scatterer

B_{2}

is in the double

η

-coherence band of

B_{1}

Figure 6. Target model in the simulation. The dynamic range is 6.

Figure 7. Received signal for a signal-to-noise ratio of 5 dB.

Figure 8. Negative spectral entropy for

γ = 1

(a);

γ = 3

(b);

γ = 5

(c);

γ = 7

(d).

Figure 8. Negative spectral entropy for

γ = 1

(a);

γ = 3

(b);

γ = 5

(c);

γ = 7

(d).

Figure 9. Contrast of entropy feature versus

γ

Figure 9. Contrast of entropy feature versus

γ

Figure 10. Negative spectral entropy (a) and entropy feature (b) for

γ = 2.5

Figure 10. Negative spectral entropy (a) and entropy feature (b) for

γ = 2.5

Figure 11. Detected boundaries, where

t_{s, 1} = 4.5070 \times 10^{- 6}

t_{e, 1} = 4.5495 \times 10^{- 6}

t_{s, 2} = 4.6265 \times 10^{- 6}

s and

t_{e, 2} = 4.6658 \times 10^{- 6}

s. The solid and dashed lines denote the true boundaries and detected boundaries, respectively.

Figure 11. Detected boundaries, where

t_{s, 1} = 4.5070 \times 10^{- 6}

t_{e, 1} = 4.5495 \times 10^{- 6}

t_{s, 2} = 4.6265 \times 10^{- 6}

s and

t_{e, 2} = 4.6658 \times 10^{- 6}

s. The solid and dashed lines denote the true boundaries and detected boundaries, respectively.

Figure 12. Distribution of the normalized modified radiation field and estimated target strips.

Figure 13. Imaging results using the band-excluded orthogonal matching pursuit algorithm (a) and the matched filter reconstruction algorithm (b) when F = 1 for the targets in the first target strip. The centers of the green ellipses denote the true positions of the scatterers.

Figure 14. Imaging results by the OMP (a,b) and BOMP (c,d) when F = 4 for the targets in the first target strip. The centers of the green ellipses denote the true positions of the scatterers, and the regions represented by the double 0.8-coherence bands of the object supports are covered by the ellipses. The semi major axis and semi minor axis of each ellipse are 0.8 m and 0.5 m, respectively.

Figure 15. Imaging results of the MF reconstruction algorithm when F = 1 (a) and those of OMP (b) and BOMP (c) when F = 4 for the targets in the second target strip.

Figure 16. Results of the imaging performance experiments. (a) the relative gridding error is approximately inversely proportional to the refinement factor; and (b) success probability for target reconstruction by the BOMP and OMP versus the refinement factor F. A reconstruction is considered successful if every reconstructed object is in the double 0.8-coherence band of the object support. The SNR is fixed at 5 dB.

Table 1. Success rates of endpoint detection for different noise levels.

$γ$	Success Rate (%)
$γ$	0 dB	3 dB	5 dB
1	75	96	98
5	92	95	99
$γ = \arg \max_{γ} E F_C (γ)$	99	100	100

© 2017 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 (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Tian, C.; Jiang, Z.; Chen, W.; Wang, D. Adaptive Microwave Staring Correlated Imaging for Targets Appearing in Discrete Clusters. Sensors 2017, 17, 2409. https://doi.org/10.3390/s17102409

AMA Style

Tian C, Jiang Z, Chen W, Wang D. Adaptive Microwave Staring Correlated Imaging for Targets Appearing in Discrete Clusters. Sensors. 2017; 17(10):2409. https://doi.org/10.3390/s17102409

Chicago/Turabian Style

Tian, Chao, Zheng Jiang, Weidong Chen, and Dongjin Wang. 2017. "Adaptive Microwave Staring Correlated Imaging for Targets Appearing in Discrete Clusters" Sensors 17, no. 10: 2409. https://doi.org/10.3390/s17102409

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

Article Menu

Adaptive Microwave Staring Correlated Imaging for Targets Appearing in Discrete Clusters

Abstract

1. Introduction

2. MSCI Model

3. Adaptive Imaging Method for Targets Appearing in Discrete Clusters

3.1. Location of the Target Strips

3.2. Location of the Regions with Targets in a Target Strip

3.3. Target Reconstruction

3.3.1. Grid Refinement

3.3.2. Target Reconstruction Method

4. Simulation Results and Discussion

4.1. Location of the Target Strips

4.2. Target Reconstruction

4.3. Performance Simulations

5. Conclusions

Acknowledgments

Author Contributions

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI