# Radar Coincidence Imaging under Grid Mismatch

Authors: Dongze Li; Xiang Li; Yongqiang Cheng; Yuliang Qin; Hongqiang Wang
Journal: ISRN Signal Processing
 (2014)

---

## Abstract

Radar coincidence imaging is an instantaneous imaging technique which does not depend on the relative motion between targets and radars. High-resolution, fine-quality images can be obtained using a single pulse either for stationary targets or for complexly maneuvering ones. There are two image-reconstruction algorithms used for radar coincidence imaging, that is, the correlation method and the parameterized method. In comparison with the former, the parameterized method can achieve much higher resolution but is seriously sensitive to grid mismatch. In the presence of grid mismatch, neither of the two algorithms can obtain recognizable high-resolution images. The above problem largely limits the applicability of radar coincidence imaging in actual imaging scenes where grid mismatch generally exists. This paper proposes a joint correlation-parameterization algorithm, which uses the correlation method to estimate the grid-mismatch error and then iteratively modifies the results of the parameterized method. The proposed algorithm can achieve high resolution with fine imagery quality under the grid mismatch. Examples are provided to illustrate the improvement of the proposed method.


---

## Body

## 1. Introduction

Radar coincidence imaging is developed as the extension of classical coincidence imaging in microwave radar systems [1]. Classical coincidence imaging, which is realized in thermal optical systems, is a method to nonlocally image an object with high resolution by producing sharply fluctuating signals in the imaging plane [1, 2]. Similarly, the essential principle of radar coincidence imaging is to produce time-space independent signals in the detecting area. Then radar signals at different positions have mutually independent waveforms. Illuminated by such signals, target-scattering centers within a radar beam will also reflect echoes of independent waveforms associated with their respective positions. Therefore, the echo component of each scattering center can be extracted from the receiving signal and then correlated to their respective positions. Consequently, scattering centers within a beam can be resolved and the spatial distribution of the target scattering centers can be obtained. Obviously, it is quite different from the radar imaging techniques based on the range-Doppler (RD) principle, where scattering centers are resolved based on the analysis of time-delay and Doppler frequency [3].Because of the different imaging principles, the radar coincidence imaging technique has two advantages over most of the RD imaging methods. Firstly, radar coincidence imaging does not depend on the aspect-angle integration or the Doppler gradient to achieve high azimuth resolution. Thus, it does not require relative motion between radars and targets and can obtain images of the targets which remain stationary with respect to radars. Furthermore, radar coincidence imaging can achieve high-resolution using a single pulse. The extremely short imaging time, which is shorter than a pulse width, considerably decreases the impact of the noncooperative motion to imagery quality. Therefore, radar coincidence imaging can obtain high-resolution, fine-quality images either for stationary targets or for the ones in complex maneuvers.In radar coincidence imaging, the target area needs to be discretized to a grid and target-scattering centers are assumed to be located at the grid points. If scattering centers are located off the grid points, then the grid mismatch yields. There are two main image-reconstruction methods used in radar coincidence imaging. One is the correlation method with lower resolution. The other is the parameterized method which can produce much higher resolution but is too sensitive under grid mismatch to give recognizable target images. Furthermore, current algorithms that are applicable to solve the sensitivity to basis mismatch are unfortunately ineffective for the grid mismatch in radar coincidence imaging [4, 5]. This limitation seriously affects the applicability of the imaging method in actual scenarios where the grid mismatch generally exists. This paper proposes a joint correlation-parameterization method for image reconstruction which can achieve both high resolution and good imagery quality under grid mismatch.The paper is organized as follows. Section2 is devoted to the basic principles of radar coincidence imaging. Section 3 analyzes the impact of grid mismatch and proposes the joint correlation-parameterization method for image reconstruction. Section 4 concludes the work.


## 2. Basic Principles

Coherent signals, which are widely employed by most of the imaging radars, generally produce detecting signals that show significant spatial correlation, as shown in Figure1(b). However, the essence of radar coincidence imaging is to produce time-space independent radar signals in the detecting area. In other words, the detecting signals at different instants or in different positions are independent of each other. As shown in Figure 1(a), at an arbitrary instant, the time-space independent detecting signals fluctuate sharply in the imaging plane and its instantaneous spatial distribution presents remarkable variety.(a) The spatial distribution of time-space detecting signals. (b) The spatial distribution produced by coherent transmitted signals.
(a)
(b)Generally, the target location can be firstly estimated based on the detection and localization techniques [6, 7], which is defined as the center of the target area. Then a local coordinate system is established in the target area center, as shown in Figure 2. The target area is discretized to a grid consisting of L small rectangles of uniform size and shape. Each small rectangle is the grid cell and is approximated by its own center. Thus, target-scattering centers are actually assumed to be initially located at the grid-cell centers. The discrete target area is expressed as I={r1,r2,…,rL}, where rl is the position vector of the lth grid cell center. There are N transmitters and a receiver, whose position vectors are Rn and Rr, respectively. The transmitted signal of the nth transmitter is denoted as Stn(t) and the receiving signal is Sr(t). The detecting signals distributed in the target area are labeled as SI(r,t), where r is the position vector of an arbitrary grid cell within I.Figure 2
Geometry of the target area grid.The ideal transmitted signals for radar coincidence imaging are supposed to be group-orthogonal and time-independent as denoted in (1), which can produce the time-space independent SI(r,t) via single transmitting [1]. Consider
(1)RT(n1,n2;t1,t2)=∫Stn1(t-t1)Stn2(t-t2)dt=δ(n1-n2,t1-t2),
where Stn(t)=rect(t/Tp)·stn(t), stn(t) is the envelope function and Tp is the pulse width. If the antenna number is more than 2, for example, N=5, the detecting signal produced by such transmitted signals will show approximate time-space independence [1], expressed as follows:
(2)RI(r,r′;τ,τ′)=∫SI(r,t-τ)SI*(r′,t-τ′)dt~Nδ(r-r′,τ′-τ).
Then, we express the receiving signal as the superposition of the detecting signals. Consider
(3)Sr(t)=∑l=1LσlSI(rl,t-|rl-Rr|c),
where σl is the scattering coefficient of the scattering center located at the lth grid cell and for the grid cell without target-scattering center σl=0. For the sake of simplicity, the coincidence imaging formulism needs a reference signal [1], which herein can be simply structured using SI(r,t) as follows:
(4)S(r,t)=SI(r,t-|r-Rr|c).
Consequently, (3) becomes
(5)Sr(t)=∑l=1LσlS(rl,t).
Then, the scattering coefficient of an arbitrary grid cell rx can be explicitly obtained via the correlation between the receiving signal and S(rx,t),
(6)∫Sr(t)S*(rx,t)dt=∫(∑l=1Lσl·SI(rl,t-|rl-Rr|c))·SI*(rx,t-|rx-Rr|c)dt=∑l=1Lσl·RI(rl,rx;|rl-Rr|c,|rx-Rr|c)~∑l=1Lσl·Nδ(rl-rx)=N·σx.
That is,
(7)σx~1N∫Sr(t)S*(rx,t)dt.Note that the detecting signalSI(r,t) or S(r,t) can be computed based on the known transmitted waveforms. Thus, the scattering coefficient of each grid cell can be derived via the correlation between the receiving signal and the reference signal. Therefore, the target scene can be recovered via conducting the correlation in (7) for the whole target area grid. This method to reconstruct target images is defined as the correlation method.The excellent point-to-point relationship in (7) means a high resolution, but it requires that the transmitted signals have the perfect time-independence as presented in (1). Unfortunately, microwave transmitted signals are far from this requirement and their time-independent degree is inadequate to achieve the high resolution represented in (7) via the correlation method [1]. Therefore, to improve the resolution of radar coincidence imaging, the parameterized method is employed which is less constrained by the signal time-independence. This method uses the relationship between the receiving signal and the reference signal, as shown in (5), to structure a radar coincidence imaging equation, given as follows:
(8)Sr=S·σ,[Sr(t1)Sr(t2)⋮Sr(tK)]=[S(r1,t1)S(r2,t1)⋯S(rL,t1)S(r1,t2)S(r2,t2)⋯S(rL,t2)⋮⋮⋯⋮S(r1,tK)S(r2,tK)⋯S(rL,tK)]·[σ1σ2⋮σL],
where K=L, S is the reference signal matrix,Sr is the vector of the receiving signal, and σ is the unknown vector of the scattering coefficient. The columns and the rows of S basically represent the detecting signals in different positions and at different instants, respectively. Hence, the incoherent property of S is basically determined by the time-space independence of the detecting signals. Therefore, the independent characteristic of the detecting signal will ensure a full-rank S. Consequently, the target scattering-coefficient vector can be uniquely recovered as σ=S-1·Sr.


## 3. Impact of the Grid Mismatch and Image Reconstruction

The imaging equation presented in (8) is derived under the grid-match condition. Generally, regardless of how finely the target area is gridded, target-scattering centers may not lie in the grid-cell centers. Assume that there exists the grid mismatch. The scattering-center position vector actually is rl+Δrl instead of the assumed rl for thelth grid cell. Thus, the actual receiving signal is Sr(t)=∑l=1LσlS(rl+Δrl,t). Consequently, the imaging equation should be Sr=(S+εs)·σ, where εs is caused by the difference between S(rl,t) and S(rl+Δrl,t). As εs and σ are both unknown variables, the above imaging equation can be simplified as
(9)Sr=S·σ+εm,
where εm=εs·σ and denotes the modeling error caused by grid mismatch. Then (9) represents the well-known linear signal model. Thus, the actual scattering-coefficient vector should be derived as σ=S-1·(Sr-εm). The previous solution derived from (8) is relabeled as σ^. That is, Sr=S·σ^, where εm is neglected because of the grid-match assumption. Certainly, the imaging quality of the parameterized method will be decreased due to the modeling error.To illustrate the grid-mismatch impact, we give an example to show how the position bias affects the modeling error and how the modeling error affects the imaging quality. In the quantitative manner, the imaging quality discussed here is firstly indicated by the relative imaging error, expressed as∥σ^-σ∥/∥σ∥. Then the modeling error is normalized as the signal-error-ratio, denoted as ζ=20log⁡10∥Sr∥/∥εm∥. The position bias is defined as ∑l=1L|Δrl|/lc, where lc is the grid-cell perimeter.The example employs anN-transmitter 1-receiver array, consisting of 5 antennas. The arrangement of the antennas and the target is shown in Figure 3. The target area is 2 m × 2 m and is discretized to 40 × 40 grid cells, each of which is 5 cm × 5 cm. The group-orthogonal and time-independent transmitted signals herein are generated by modulating the sinusoid signals on amplitude via mutually independent Gaussian noises. The carrier frequency is 9.5 GHz, the pulse width is 50 μs, the bandwidth is 1 GHz, and the sampling frequency is 2 GHz. The results are shown in Figure 4.Figure 3
The arrangement of antennas and the target model.Figure 4
The grid-mismatch impact to imaging quality of the parameterized method. (a) The imaging error versusζ. (b) The imaging result when ζ=26 dB. (c) The imaging result when ζ=46 dB. (d) The signal-error-ratio versus the position bias.In this example, the target image in Figure4(b) is badly blurred beyond reorganization when ζ=26dB. The target image in Figure 4(c) is recognizable when the imaging error is 0.9, which requires that ζ should not be lower than 46 dB. It demonstrates that the parameterized method has a low tolerance to the modeling error. A recognizable image with a small imaging error requires a small modeling error or a high signal-error-ratio. Then the question is “will the position bias generally produce a high-level signal-error-ratio to guarantee low imaging error?” Unfortunately, as shown in Figure 3(d), even a very small position bias of 0.001 in this example can cause a large modeling error, which only provides 20 dB single-error-ratio. Moreover, the position bias less than 0.02 will generate a negative ζ. It shows that the modeling error is seriously sensitive to the position bias. The position bias of 0.04 can cause −4 dB signal-error-ratio in radar coincidence imaging equation, with no mention of a much larger bias in practice. In this case, modeling error almost overwhelms the receiving signal in the imaging equation. The experiment indicates that the recovery, simply depending on the imaging equation, might impossibly provide a good estimation under such a poor signal-error-ratio condition.Consider another example to compare the correlation method with the parameterized method under grid mismatch. In the example, target images will be reconstructed via the two methods with and without the grid mismatch. Herein, the grid mismatch leads toζ=-4dB. Imaging results are shown in Figure 5.The imaging results recovered by different methods. (a) The correlation-method result without grid mismatch. (b) The parameterized-method result without grid mismatch. (c) The correlation-method result whenζ=-4 dB. (d) The parameterized-method result when ζ=-4 dB.
(a)
(b)
(c)
(d)As shown in Figure5, the parameterized method provides higher resolution than the correlation method in the case of no grid mismatch. On the other side, when ζ=-4dB, the imaging quality of the correlation method is almost unchanged but the parameterized-method result gets degraded and almost gives no information on the target shape. It is obvious that the correlation method is less sensitive to the modeling error than the parameterized method but has much lower resolution.In conclusion, neither of the two image-reconstruction methods can achieve both high resolution and good imaging quality in the presence of grid mismatch. According to their respective limitations, there are two ways to solve the image reconstruction under grid mismatch. One way is to improve the resolution of the correlation method. The other one is to decrease the grid-mismatch impact in the parameterized method. As stated previously, however, microwave signal does not have adequate time-independence in nature, resulting in a limited resolution. Thus, the latter might be a possible way to solve the problem of grid mismatch.A direct idea to decrease grid-mismatch impact in the parameterized method is to estimate the modeling error. Note that the modeling error is relevant to the scattering center number, scattering intensity, and the position-bias level of all scattering centers. If there is no prior information of the target shape, the modeling error will be quite difficult to estimate. Then for the targets that are unknown in advance, the obtainable shape information could be sought from the imaging results of the two reconstruction methods. Thus, a possible source to provide the desired knowledge would be the correlation-method result, considering the seriously distorted images of the parameterized method caused by grid mismatch.The knowledge that can be used by the imaging equation should be quantitative and precise. However, the result of the correlation method, as observed in Figure5, presents a blurry target shape. The image cannot precisely indicate the target scattering-center location or provide the exact value of the scattering coefficient. Nonetheless, it indeed provides valuable information. That is, the correlation-method result can denote part of the grid cells which do not contain target-scattering centers. In Figures 5(a) and 5(c), only a minority of image patches, which constitute the target contour, have large values. By contrast, a majority of other image patches have far too small values, which implies the corresponding grid cells have no scattering centers. Hence, these corresponding elements in the scattering-coefficient vector are estimated to be zero. Such quantitative knowledge will be helpful for the estimation of the modeling error. Detailed implementation of image reconstruction using the idea will be described in the following paragraphs.Firstly, still under the grid-match assumption, we estimate the target scattering-coefficient vector via the two reconstruction methods:σ^c is derived via the correlation method and σ^ is obtained via the parameterized method. According to Sr=S·σ^ and (9), we have the following relationship:
(10)S·σ^-S·σ=εm⟶σ^-σ=S-1εm,
where S-1 is written as [α1,α2,…,αL]T, αl being the lth row-vector of S-1. The next step is to select the grid cells without scattering centers according to σ^c. In comparison with the average value of σ^c, part of the elements in σ^c, which have far too small values, are determined to be 0. The index of these zero-elements is denoted as Λ={l1,l2,…,lM}; that is, σlm=0,lm∈Λ. According to Λ, we define SΛ-1=[αl1,αl2,…,αlM]T, σ^Λ=[σ^l1,σ^l2,…,σ^lM]T, and σΛ=[σl1,σl2,…,σlM]T. Based on (10), we have
(11)σ^Λ-σΛ=SΛ-1·εm.
As σΛ=0, finally σ^Λ=SΛ-1·εp. Therefore, the modeling noise is estimated as
(12)ε^m=SΛ†σ^Λ,
where SΛ† is the pseudoinverse matrix of SΛ-1. By substituting ε^m to (9), the new target scattering-coefficient vector obtained is
(13)σ^(1)=S-1(Sr-ε^m).If the result remains unsatisfactory, we can repeat the computation of (12) and (13) using the new σ^(1) and ε^m. That is, update the receiving signal as Sr(1)=Sr-ε^m. Then derive σ^Λ(1) according to Λ and update the modeling error as ε^m(1)=SΛ†σ^Λ(1). Finally, the estimation is updated as σ^(2)=S-1(Sr(1)-ε^m(1)). The update can be performed iteratively until the result is satisfactory. This image-reconstruction method using both the correlation method and the parameterized method is defined as the joint correlation-parameterization method. The detailed steps are given in Table 1.Table 1
The flow of the joint correlation-parameterization method.
Step 1
Set the iteration step maximumTmax⁡ and the stop limitation η, k=0.
Step 2
Discretize the target area and compute the reference signal.
Step 3
Obtain the scattering-coefficient vectorσ^c using the correlation method.
Step 4
Set a thresholdσT to determine Λ={l1,l2,…,lM}, where σ^c(lm)<σT.
Step 5
Compose the coincidence imaging equation. Obtain the initial estimationσ^(k) using the parameterized method. Set  Sr(k)=Sr.
Step 6
Compute the inverse matrix of the reference matrixS, and then according to Λ obtain SΛ-1=[αl1,αl2,…,αlM]T.
Step 7
Deriveσ^Λ(k)=[σ^l1(k),σ^l2(k),…,σ^lM(k)]T.
Step 8
Compute the modeling noiseε^(k)m=SΛ†σ^(k)Λ, k=k+1.
Step 9
Update the scattering-coefficient vector asσ^(k)=S-1(Sr(k-1)-ε^m(k-1)).
Step 10
Stop if∥σ^(k)-σ^(k-1)∥≤η or k≥Tmax⁡. Otherwise Sr(k)=Sr(k-1)-ε^m(k), and go to Step 7.To examine the correlation-parameterization method, target images are reconstructed according to the iteration process of Table1, where the grid mismatch is the same as that in Figure 5(d). The imaging results are given in Figure 6.The imaging results of the joint correlation-parameter method. (a) The imaging result of the first iteration. (b) The imaging result of the second iteration. (c) The imaging result of the third iteration.
(a)
(b)
(c)As shown in Figure6, a high-quality target image can be obtained through three iterations. The result of the first iteration has presented the target contour, though the contour is light and the imaging error remains remarkable. After the second iteration, the imaging quality is largely improved. By contrast, the imaging quality of Figure 6(b) is much better than the result in Figure 5(d). The third iteration maintains the good imagery quality. Therefore, two iterations actually give a high-resolution and fine-quality target image. It indicates that the joint correlation-parameterization method satisfies the requirement for both the resolution and the robustness in the presence of the modeling noise caused by grid mismatch.


## 4. Conclusions

Radar coincidence imaging can achieve excellent high-resolution target images on the condition of grid match, but the imagery quality gets degraded beyond recognition in the presence of the modeling error caused by grid mismatch. Therefore, the paper proposes the joint correlation-parameterization method for image reconstruction. The proposed algorithm iteratively modifies the parameterized-method result with the estimated modeling error, which is obtained based on the correlation-method result. Consequently, the grid-mismatch impact on the imaging quality is considerably reduced. The example shows that the joint correlation-parameterization method can achieve high resolution and maintain the robustness under grid mismatch.



---

*Source: 987803-2014-04-22.xml*