Research Article  Open Access
Xiaoli Zhou, Hongqiang Wang, Yongqiang Cheng, Yuliang Qin, "OffGrid Radar Coincidence Imaging Based on Variational Sparse Bayesian Learning", Mathematical Problems in Engineering, vol. 2016, Article ID 1782178, 12 pages, 2016. https://doi.org/10.1155/2016/1782178
OffGrid Radar Coincidence Imaging Based on Variational Sparse Bayesian Learning
Abstract
Radar coincidence imaging (RCI) is a highresolution staring imaging technique motivated by classical optical coincidence imaging. In RCI, sparse reconstruction methods are commonly used to achieve better imaging result, while the performance guarantee is based on the general assumption that the scatterers are located at the prediscretized gridcell centers. However, the widely existing offgrid problem degrades the RCI performance considerably. In this paper, an algorithm based on variational sparse Bayesian learning (VSBL) is developed to solve the offgrid RCI. Applying Taylor expansion, the unknown true dictionary is approximated accurately to a linear model. Then target reconstruction is reformulated as a joint sparse recovery problem that recovers three groups of sparse coefficients over three known dictionaries with the constraint of the common support shared by the groups. VSBL is then applied to solve the problem by assigning appropriate priors to the three groups of coefficients. Results of numerical experiments demonstrate that the algorithm can achieve outstanding reconstruction performance and yield superior performance both in suppressing noise and in adapting to offgrid error.
1. Introduction
Motivated by classical coincidence imaging, radar coincidence imaging (RCI) is a novel highresolution imaging technique which has been realized in optical systems [1–3]. The RCI can realize staring imaging without the requirement of the target relative motion and operate under the nonideal observing geometry of forwardlooking/staring and shorten the imaging time to even a single pulse width, with significant potentials for resolution enhancement, interference, and jamming suppression. The basic principle of RCI is to excite timespace independent signals in the imaging area. Thus the spatial variety of wavefront is increased, and the scatterers within a beam then reflect independent signals associated with their respective positions, so the superresolution within a beam emerges.
To achieve better imaging performance, sparse reconstruction is commonly used, and the continuous space needs to be discretized to a fine grid and targetscattering centers are assumed to be exactly located at the center of these prediscretized grids [3]. Based on the model, the independent detecting signals at different gridcell centers can be formed as the atoms of sparse representation dictionary. Meanwhile, the scatterers of target are often distributed sparsely in most radar imaging applications; thus sparse reconstruction and compressive sensing (CS) [4, 5] are suitable for RCI by exploiting the sparsity of target. In sparse reconstruction theory, signal reconstruction depends on presetting an appropriate sparsifying dictionary which is supported on the prediscretized grids and defines the signal sparsity. However, as the scatterers distributed in a continuous scene are generally located off the gridcell centers, the offgrid problem emerges no matter how fine the grid is [6], and the performance of RCI would degrade significantly.
The effect of general dictionary mismatch, which is the direct consequence of offgrid effect, is analyzed in [7–10]. This mismatch causes the performance of conventional sparse reconstruction methods to degrade considerably [3, 6, 7, 11–14]. An intuitive way to sidestep offgrid effect is to work directly on the continuous parameter space, that is, atomic norm minimization approach [15], continuous basis pursuit (CBP) [16]. Considering the offgrid problem, several algorithms have been proposed to alleviate the effect. One simple approach is to use refinement strategy and decrease the grid size [17]. Nevertheless, a finer grid may enhance the coherence between the columns of dictionary and increase the computational complexity and instability of reconstruction [14]. Modeling the offgrid problem as a multiplicative perturbation, the sparse total least squares (STLS) [18] and joint correlationparameterization (JCP) [3] algorithms are proposed. To explore the structure of dictionary mismatch, the supportconstrained orthogonal matching pursuit (SCOMP) [13] and joint sparse signal recovery methods [19] are proposed based on the firstorder Taylor expansion to utilize the support constraint. Lately, from the sparse Bayesian learning (SBL) perspective, several approaches are proposed, such as offgrid sparse Bayesian inference (OGSBI) [20], sparse adaptive calibration recovery via iterative maximum a posteriori (SACRiMAP) [6], and block SBL (BSBL) [21]. The merit of SBL is its flexibility in modeling sparse signals that can not only promote the sparsity but also exploit the possible structure of the signal to be recovered [20].
Although the performance of OGSBI and SACRiMAP is outstanding, the prior distribution of offgrid error (OGE) is not fully utilized; the point estimate rather than distribution of OGE is obtained. Adopting the ideas of variational expectationmaximization (EM) and variational Bayesian inference (VBI) described in [7, 22, 23], the offgrid RCI is investigated in the framework of variational sparse Bayesian learning (VSBL). As for SBL, inference in SBL model is not tractable in closed form; thus approximations are needed, such as maximizing the marginal likelihood (MML), EMSBL, and VBI [24]. Compared to other approximations, VBI has several advantages when applied to SBL. First, the distributions rather than point estimates of the unobserved variables can be obtained. Second, VSBL allows obtaining analytical approximations to the posterior distributions of interest even when their exact expressions are intractable. Finally, the VSBL methodology allows using different prior distributions and is a deterministic approximate inference framework that can be applied to many models.
In this paper, based on the firstorder Taylor expansion, the offgrid RCI model in rangeazimuth space is reformed to be sparsely approximated using atoms from three different dictionaries, and meanwhile the three groups of approximation coefficients share the same support. By assigning appropriate priors to the approximation coefficients, such a model can be conveniently manipulated to recover sparse coefficients under the VSBL framework. Compared with the aforementioned algorithm discussed above, the proposed algorithm imposes a group structure on the coefficient vector and explores the groupsparse structure using VSBL. Numerical experiments show that the algorithm realizes the target reconstruction robustly and achieves both highresolution and outstanding imaging quality in the presence of offgrid scatterers and is also simple to implement.
The rest of the paper is organized as follows. In Section 2, the offgrid RCI model in the rangeazimuth space is presented. Section 3 presents the offgrid variational sparse Bayesian learning (OGVSBL) algorithm in detail. In Section 4, the performance of the presented algorithm is verified by numerical examples and compared with some existing sparse reconstruction methods. Finally, some conclusions and future work are discussed in Section 5.
Notations used in this paper are as follows. We use boldface lowercase letters for vectors and boldface uppercase letters for matrices. , , , and denote the transpose, conjugate transpose, inverse, and pseudoinverse of a vector or matrix, respectively. and are the Hadamard product and vectorization operation, separately. and denote the norm and norm of a vector. is the determinant of a matrix; is a matrix with the elements of a vector on the main diagonal. Finally, denotes the expectation of a variable.
2. Problem Formulation
2.1. Signal Model
As a novel imaging technique, RCI has shown its potentials in highresolution, staring, and instantaneous imaging [1]. In RCI, the target is illuminated by timespace independent signals, and then the echo components reflected by scatterers at different positions are mutually independent, which could result in the superresolution within a beam. The RCI can be realized by a multipleinput multipleoutput (MIMO) radar system to transmit timeindependent and grouporthogonal waveforms (e.g., stochastic waveforms) [1]. While compared with conventional MIMO radar which focuses on multiple paths or multiple observation angles, RCI needs the interference of transmitted waveforms to make the wavefront show spatial fluctuations, thus the spatial variety of detecting signals increases. In addition, the components of each path are separated utilizing the waveform orthogonality in conventional MIMO radar, while the components are not separated in the whole RCI procedure.
In this paper, a RCI system with transmitters and receivers is considered; each transmitter emits an independent stochastic waveform. Without loss of generality, the target is assumed to be composed of several ideal point scattering centers without the spatial property for an enough high carrier frequency, which is widely used in the imaging radar system. Furthermore, the scatterers are widely separated to provide some performance guarantee.
The target scene is considered to be a 2D space using polar coordinate as illustrated in Figure 1. The sparsebased RCI methods discretize the continuous target scene and generate a number of gridcells. Thus the scene is discretized in rangeazimuth space with azimuth cells, range cells, and associated cell size , . Hence, the number of gridcells is . Assume that the scatterers are initially located at the gridcell centers; thus the th scatterer is located at the rangeazimuth pair . As the scatterers possess nonzero complex scattering coefficient which is proportional to the radar cross section (RCS), the associated means that there is no scatterer at the th gridcell center. Denote by the index set of scatterer locations.
Note that the echo is a linear combination of all the scatterers reflected waveforms from all the transmitters. As it is unknown a priori which is nonzero, we sum over all possible locations, although contributions to the echo only result from the scatterers. Then, the echo at the th receiver can be expressed aswhere is the signal emitted by the th transmitter, is an independent complex Gaussian random process with zeromean, and variance denotes the background noise at the th receiver. is the relative propagation delay corresponding to the th transmitter and th receiver with respect to the th scatterer, which is defined aswhere and are the position vectors of the th transmitter and th receiver, respectively. is the speed of wave propagation.
For the sake of simplicity, the RCI formula needs a detecting signal [1], which can be simply structured as
Then, the echo can be expressed as the superposition of the detecting signals
Sample the echo and obtain samples; then we stack and , via the vectorization operation , into dimensional column vectors and , respectively. And
The column vector is the atom corresponding to the th gridcell. Then we arrange into an dimensional matrix , which is often referred to as a dictionary matrix and defines the basic elements of sparse representation.
Denote by an unknown vector to be recovered. Then the measurement model reduces to (6), which is a familiar linear model used in most applications of sparse modeling:
can be deduced from the prior knowledge of the measurement model. Although the scattering coefficient vector can be derived via the correlation between the echo and dictionary matrix directly, the sparse reconstruction methods are preferred to improve the resolution.
2.2. OffGrid RCI Model
The imaging equation presented in (6) is derived under the grid matching condition. However, in most cases of radar imaging, the scatterers are distributed in a continuous scene. Regardless of how finely the imaging scene is discretized, the scatterers may not lie at the gridcell centers, and then the offgrid problem yields. Considering the offgrid problem, the actual location of the th scatterer can be expressed aswhere is the OGE. denotes the normalized OGE which varies from to .
The offgrid effect can be compensated by estimating the perturbation . In other words, the offgrid RCI is a joint estimation problem to recover associated with the nonzero . As the estimation of is a nonlinear optimization problem which is analytically intractable, linearization through Taylor expansion is commonly used. Thus the true unknown dictionary is approximated, and then an accurate linear model for sparse reconstruction is established to reduce the OGE in the conventional fixeddictionary setting.
Since the OGE is restricted in one gridcell which is significantly small, then we make the following approximation by the firstorder Taylor expansion:
Denote , , , , and , . The model in (6) could be concisely rewritten as
Alternatively, (9) can be expressed aswhere and . Clearly, and are also sparse vectors and share the same support with . Hence, (10) means that the imaging equation can be sparsely approximated using three known dictionaries . Now, the problem of reconstructing the target image becomes recovering , , and simultaneously based on (10). We approach this joint sparse recovery problem under the VSBL framework.
3. Target Reconstruction: Principle and Algorithm
The key to reconstruct the target image successfully is that all scatterers fall on the prediscretized gridcell centers. However, the generally existing offgrid problem jeopardizes the reconstruction performance seriously. In this section, VSBL is modified and then applied to offgrid RCI by reconstructing , , and sequentially. This modified algorithm alleviates the offgrid effect and realizes the grid matching which means that the scatterers falling off the gridcell centers are captured by the closest gridcell centers.
3.1. TwoLayer Hierarchical Prior Model
We begin with the specification of the graphical representation of model priors for model (10), as depicted in Figure 2.
To recover , , and in (10), we model them as hidden variables and assign them appropriate priors. Typically, a twolevel hierarchical GaussianGamma prior is imposed on sparse coefficients to induce sparsity and also allows for a convenient conjugateexponential analysis. Thus the following hierarchical prior is assigned:where denotes the complex Gaussian distribution with mean vector and covariance matrix . Consider , , is the prior variance of , and is the Gamma distribution with parameters and . The prior parameters are also called sparsity parameters, being inversely proportional to the width of the pdf. Naturally, a large value of will drive the corresponding weight to zero, thus inducing a sparse solution. As , , , and the prior variance of is ; then the prior variances of and can be set to be to accommodate the worst case, where the true scatterer is at the middle of the two consecutive grids. Besides, this setting satisfies the constraint of common support among , , and . Consequently, the priors of and conditioned on can be expressed as and , respectively.
In addition, the measurement noise is assumed to be independent and complex Gaussian with zeromean and variance , to allow conjugateexponential analysis:
3.2. OffGrid Variational Sparse Bayesian Learning
Based on the model presented above, (10) becomes a model containing only hidden variables and no parameters. Then we employ the SBL to realize the target reconstruction. However, the derivation of the posterior is generally intractable in the SBL framework. Thus VBI is adopted in this case to obtain the variational solution to SBL.
VBI is a methodology that can deal with complicated Bayesian models [23], aiming at finding maximum likelihood estimates of unknown parameters involved in a Bayesian model with hidden (i.e., unobserved) variables. VBI first uses variational methods to find approximate posterior distributions of hidden variables and then updates the estimates of unknown parameters based on the posterior densities of hidden variables. Specifically, the variational solution is obtained by maximizing a variational lower bound on a logevidence [25, 26], which is equivalent to minimize the KullbackLeibler divergence (KLD) between the approximating pdf and the posterior pdf . Moreover, a structured mean field approximation can be applied to compute the approximate posterior :where automatically follows from the independence assumption (pdf factorization) among , , , and posterior approximations. Meanwhile, VBI obtains the posterior densities via [24]where is the set of hidden variables, denotes the posterior density for , is the joint pdf defined as , denotes the expectation with respect to , and is a normalized constant.
By using (14), the variational approximating factors of can be obtained as , , , , and . The parameters of the approximating factors, that is, variational parameters, can be computed as follows:where , , , , , represents the th main diagonal element of a matrix, and is the th element of a vector. Using the property of Gamma distribution, the expected values can be computed as
Based on the above results, the posterior distributions of can be obtained by iteratively performing (15)–(20) until convergence. During that process, the parameters , , , and are typically fixed at values close to zero to assume noninformative priors for and . Alternatively, the update for and can be changed as follows:
The updates in (21) are observed to lead to much faster convergence although they do not benefit from the guarantee of local optimum of object function. It is remarkable that each can be interpreted as a measure of how “welldetermined” its corresponding parameters , , and are by the data. When is large, where , , and are highly constrained by the prior, and it follows that ; otherwise, , , and fit the data, .
Due to the convexity of the variational lower bound in the approximating factors , , , , and , we can update these factors in any order while keeping the other factors fixed. However, the update order is important since different update orders might lead to different local optima of the variational lower bound.
In addition, we terminate the algorithm if or the maximum number of iterations is reached, where is a userdefined tolerance and the superscript refers to the iteration.
3.3. Fast Implementation
The above algorithm suffers from two main drawbacks: high computational complexity and low convergence rate. As each iteration consists of matrix inversion and matrixvector multiplication in (15)–(17), the costs of the two operations are and , respectively. Hence, the procedure is computationally expensive to implement when is large. This exists generally since a dense grid is required in principle to reduce the error incurred by the linear approximation. In order to overcome the drawbacks, a fast implementing approach is developed.
In conventional SBL, many of the prior precisions of coefficients typically take on quite large values upon convergence. This implies that the corresponding coefficients are quite small and in turn the contributions of the corresponding atoms are negligible; thus these atoms can be pruned reasonably to realize sparsity [27]. Then we maintain this feature in the algorithm described above.
From (15)–(17), it can be shown that when is large enough, , , , , , and will be close to 0, which means that the contribution of the th grid to signal synthesis can be neglected, and thus this grid could be removed from the current grid set. Based on this observation, grid pruning can be conducted through the threshold of . Specifically, denoting a large threshold , the pruning of the current grid set can be achieved viawhere denotes the grid set after pruning. It should be mentioned that, due to the grid pruning, the posterior mean vectors (, , and ), covariance matrix (, , and ), and the dictionaries (, , and ) should be also pruned accordingly.
Because of the grid pruning, the computational burden is reduced and the convergence speed is improved, the sparsity of the reconstruction results is enhanced, and meanwhile the common support shared by the scattering coefficients and OGE is guaranteed. Thus the reconstruction performance would be better, whereas the threshold needs to be selected carefully; otherwise the grid would be pruned mistakenly.
3.4. Algorithm Summarization
Based on the above discussions, the procedure of the proposed algorithm, that is, offgrid variational sparse Bayesian learning (OGVSBL), is given in Algorithm 1.
Next, we provide the following comments to provide further insight into the proposed algorithm.
(1) The Advantages. Similar to other sparsebased or CSbased algorithms, the OGVSBL algorithm performs well for highresolution RCI under the scenarios of low signaltonoise ratio (SNR) with only a few measurements. While compared with many sparsebased algorithms, the proposed algorithm is a datadriven learning process which needs none of the parameters tuned manually. Furthermore, the OGVSBL algorithm utilizes the hierarchical modeling procedure to encode signal sparsity which could obtain the approximate posterior distribution and is regarded as a full Bayesian method. The statistical information is used to enhance the estimation performance and avoid converging to a shallow local minimum during the learning procedure, due to the utilization of higher order statistical information.
(2) Comparison with Other Algorithms. As discussed in Section 1, various methods could solve the offgrid problem, such as bandexcluded locally optimized orthogonal matching pursuit (BLOOMP) [17], JCP [3], STLS [18], and block sparse Bayesian learning (BSBL) [21]. As a greedy algorithm, BLOOMP modifies the search strategy of OMP and improves the performance of OMP in the presence of offgrid scatterers. JCP models the offgrid imaging equation as , where denotes the modeling error caused by offgrid effect. Then JCP uses the correlation method to estimate the modeling error and then iteratively modifies the results of the parameterized method. STLS seeks to solve the nonconvex optimization problem:where is a regularization parameter and denotes the dictionary perturbation caused by offgrid effect. From the Bayesian perspective, the STLS yield an MAP solution if is Gaussian which cannot be satisfied in the offgrid RCI. These algorithms are designed for RCI with modeling error (including OGE), but they do not consider the structure of the dictionary perturbation from statistic viewpoint. It is noteworthy that and share the same support with ; thus the unknown variables can regroup into nonoverlapping subblocks and the th subblock is . Then is block sparse and can be solved in BSBL framework. Both BSBL and OGVSBL formulate the offgrid RCI using Bayesian hierarchical prior modeling, but their implementation processes are different.
4. Discussion and Numerical Analysis
In this section, the performance of the proposed offgrid RCI algorithm is analyzed by conducting a wide range of numerical experiments. The imaging performance is indicated by two criteria, that is, the reconstructed sparsity and relative imaging error (RIE) which is expressed as .
Consider an Xband RCI system with carrier frequency of 10 GHz. Both the transmitters and receivers are configured as a uniform linear array (ULA) with the numbers of elements and interelement spacing m, respectively. The transmitters transmit frequencyhopping (FH) waveforms, described in [28, 29], with the minimum FH interval MHz and maximal hopping frequency MHz. The waveforms are comprised by hopping frequencies with the FH duration ns. With the above parameters, the grid size for conventional imaging approaches can be computed, while we set a smaller grid size for the superresolving ability of sparse reconstruction approaches. A rangeazimuth target space, covering 4 m × 0.04 rad, is discretized to 40 × 40 gridcells. Further, there supposed to be six widely separated point scatterers which are randomly located within the gridcells.
In the proposed OGVSBL algorithm, we initialize , , , and . For the termination condition, we set and the maximum number of iterations to . The threshold for grid pruning is .
4.1. Performance of OGVSBL
To illuminate the validity of the proposed OGVSBL algorithm, we conduct a numerical simulation where VSBL algorithm is implemented as a comparison. Figure 3 shows the imaging results. The white circles in Figure 3(a) represent the actual locations of scatterers and the red points represent the gridcells closest to the scatterers. In this example, the external noise is fixed at dB. Figures 3(b) and 3(c) plot the reconstructed images with respect to ongrid and offgrid targets, respectively, where VSBL is utilized. When the scatterers fall on the gridcell centers, VSBL can reconstruct the target perfectly, while for offgrid target, the image is defocused as shown in Figure 3(c), where many spurious scatterers exist and the signal energy spills over the imaging plane. Compared with Figure 3(c), the image shown in Figure 3(d) is improved significantly. Six scatterers are reconstructed much better without any obvious spurious scatterer, mainly as the prior knowledge of sparsity and the OGE are fully utilized. Thus OGVSBL can be used to image the offgrid target.
(a)
(b)
(c)
(d)
Figure 4 depicts the RIEs for both the VSBL and OGVSBL decrease during the imaging process. Both algorithms are converged after about 20 iterations, indicating the validity of the algorithms. As the prior knowledge of both the sparsity and the distributions of scattering coefficient and the OGE are fully utilized, the OGVSBL algorithm improves the imaging performance by about 8 dB from the RIE perspective, compared with VSBL.
As OGVSBL is a sparse reconstructiontype algorithm whose performance depends on the sparse prior of target, we test the performance under the sparsity degree of target. In the experiment, OGVSBL is implemented for different sparsity degrees varying from 5 to 35, and VSBL is also simulated for comparison. For each sparsity degree, the scattering centers are randomly located within the gridcells and the RIE and reconstructed sparsity are recorded over 100 independent trials.
Figure 5 shows the reconstructed results. It can be observed from Figure 5(a) that the performance lies in the sparsity property of a target for both algorithms and the performance degrades as the sparsity degree increases. Besides, OGVSBL reconstructs the target precisely without spurious scatterer when the sparsity degree is not higher than 15 under the simulation parameters. For targets with higher sparsity degree, it is observed that OGVSBL generates slightly more scatterers than the true sparsity degree but still obtains much sparser image than VSBL. As seen from Figure 5(b), the RIE reconstructed by OGVSBL increases with the sparsity degree, which means the image quality tends to be worse. Then we can conclude that performance of OGVSBL is affected by the sparsity degree of the target, since it is based on the prior knowledge of sparsity.
(a)
(b)
4.2. Performance Comparison with Other Algorithms
In this subsection, the performance of OGVSBL is compared with that of other algorithms which take into account offgrid effect. These algorithms for comparison include BLOOMP, JCP, STLS, and BSBL algorithms. All these algorithms realize the grid matching which means that the offgrid scatterers are captured by the closest gridcell centers.
Figure 6 shows the original image and target images reconstructed by the five algorithms. For BLOOMP, JCP, and STLS, the reconstructed images are blurry and the signal energy spills over the target space as depicted in Figures 6(b)–6(d). Apart from the six strong scatterers, many spurious scatterers exist around the actual scatterers, which will lead to failure in picking up the scatterers without other prior knowledge (i.e., the number and scattering coefficients of scatterers). Comparably, the image reconstructed by BSBL is much sparser as presented in Figure 6(e), but a few spurious scatterers still exist. In contrast with Figures 6(b)–6(e), the image shown in Figure 6(f) is improved significantly. Six scatterers are reconstructed much better without any obvious spurious scatterer, and their support is recovered precisely, mainly as the prior knowledge of sparsity and the OGE are fully exploited. Thus, for offgrid RCI, OGVSBL shows the best performance among the algorithms tested above.
(a)
(b)
(c)
(d)
(e)
(f)
Next, we vary the SNR to evaluate the performance of OGVSBL by numerical experiments. To provide the lower bound for the listed algorithms, we define the Oracle solution for grid matching, which is the least squares solution within the true support set [30]; that is,where is the dictionary associated with the support set and is the complement of the support set .
The results are shown in Figure 7 and Table 1. Monte Carlo experiments are conducted to test the robustness of the listed algorithms to different SNRs varying from −20 to 40 dB. For each SNR, we conduct 100 independent trials and record the RIE and reconstructed sparsity. Figure 7 shows the RIEs under different SNRs. It follows that as the SNR increases, the RIEs decrease quickly to an asymptotic region where the RIEs keep nearly unchanged, implying that the performance is improved and the algorithms are sensitive to noise in low SNR regimes, while, among the five algorithms, OGVSBL shows the best performance. Besides, the RIE obtained by OGVSBL is much closer to that obtained by Oracle solution and approaches the lower bound defined by Oracle estimator especially in high SNR regimes. The reconstructed sparsity is shown in Table 1. Algorithms like BLOOMP, JCP, and STLS tend to generate nonsparse and defocused images in the presence of offgrid scatterers, since they do not exploit the offgrid problem accurately. Therefore, the signal energy spills over the target space and the spurious scatterers emerge, especially in low SNR regimes. As presented in Table 1, BSBL and OGVSBL can achieve much sparser image due to the utilization of the hierarchical Bayesian modeling; thus the image quality is improved significantly.

To make further comparisons, the performance of OGVSBL under different grid sizes is investigated, comparing with the other algorithms. In this experiment, we assume that the grid number and actual scatterer distribution are invariant; the discretized grid size is changed by grid refinement. The refinement factor , defined as the proportion of the refinement grid size and initial grid size, is used as a metric of grid size. We assign dB and conduct 50 independent trials to achieve the averaged RIEs.
Table 2 presents the simulation result where the RIE is shown to be sensitive to the refinement factor. It can be shown from Table 2 that the refinement factor induces the highest reconstruction error; all algorithms except JCP have the worst performance deterioration. As the refinement factor increases, the gridcell becomes larger which yields a high coherence between the columns of the dictionary. Then the RIP condition for reliable sparse recovery is guaranteed, and the RIEs would be lower. When the refinement factor becomes smaller, the RIEs also decrease because the modeling error caused by offgrid problem decreases accordingly, which is indicated by the Oracle solution. Furthermore, a smaller grid size would decrease the approximation error due to high order terms of Taylor expansion, which also render a smaller RIE consequently. Besides, the RIE obtained by OGVSBL is much lower than that obtained by other algorithms, which shows the performance advantage of OGVSBL. Consequently, the grid size should be chosen carefully when the imaging scene is discretized.

5. Conclusion
In this paper, we investigate the offgrid RCI based on VSBL framework. For sparsebased RCI, the offgrid effect would lead to the imaging performance deterioration significantly. To improve the imaging performance, we approximate the offgrid RCI model linearly using the Taylor expansion, and then an accurate model for sparse approximation is established. Based on this model, target reconstruction is recast as a joint sparse recovery problem that recovers three groups of coefficients over three known dictionaries subject to the constraint that the three groups share the same support. Such a problem can be readily manipulated from a Bayesian perspective; thus we developed an algorithm based on VSBL framework to solve the problem. Numerical experiments illustrate that this algorithm can exhibit outstanding performance, which shows the potential for the algorithm to be applied in practical radar system.
Although grid pruning is applied, the proposed algorithm still suffers from high computational complexity especially when the dictionary has large numbers of atoms in the case of 3D RCI or dense discretized grids. Then developing fast versions of OGVSBL are of great interest, and the fast VSBL scheme in [15] could be a candidate. Meanwhile, OGVSBL is sensitive to the RCI model parameters; then the effects of model parameters would be investigated. Besides, the algorithm achieves grid matching, while the exact locations of scatterers are ignored. Thus OGE estimation should be executed to reconstruct the target exactly, which would make the imaging results more outstanding. These open issues are beyond the scope of this paper and will be investigated deeply in the near future.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (nos. 61302149 and 61302142) and Research Fund for the Doctoral Program of Higher Education of China (20124307110013).
References
 D. Li, X. Li, Y. Cheng, Y. Qin, and H. Wang, “Radar coincidence imaging: an instantaneous imaging technique with stochastic signals,” IEEE Transactions on Geoscience and Remote Sensing, vol. 52, no. 4, pp. 2261–2277, 2014. View at: Publisher Site  Google Scholar
 D. Li, X. Li, Y. Cheng, Y. Qin, and H. Wang, “Radar coincidence imaging in the presence of targetmotioninduced error,” Journal of Electronic Imaging, vol. 23, no. 2, Article ID 023014, 2014. View at: Publisher Site  Google Scholar
 D. Li, X. Li, Y. Cheng, Y. Qin, and H. Wang, “Radar coincidence imaging under grid mismatch,” ISRN Signal Processing, vol. 2014, Article ID 987803, 8 pages, 2014. View at: Publisher Site  Google Scholar
 D. L. Donoho, “Compressed sensing,” IEEE Transactions on Information Theory, vol. 52, no. 4, pp. 1289–1306, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 S. Qaisar, R. M. Bilal, W. Iqbal, M. Naureen, and S. Lee, “Compressive sensing: from theory to applications, a survey,” Journal of Communications and Networks, vol. 15, no. 5, pp. 443–456, 2013. View at: Publisher Site  Google Scholar
 X. He, C. Liuc, B. Liu, and D. Wang, “Sparse frequency diverse MIMO radar imaging for offgrid target based on adaptive iterative MAP,” Remote Sensing, vol. 5, no. 2, pp. 631–647, 2013. View at: Publisher Site  Google Scholar
 L. Hu, J.X. Zhou, Z.G. Shi, and Q. Fu, “An EMbased approach for compressed sensing using dynamic dictionaries,” Journal of Electronics & Information Technology, vol. 34, no. 11, pp. 2554–2560, 2012. View at: Publisher Site  Google Scholar
 Y. Chi, L. L. Scharf, A. Pezeshki, and A. R. Calderbank, “Sensitivity to basis mismatch in compressed sensing,” IEEE Transactions on Signal Processing, vol. 59, no. 5, pp. 2182–2195, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 X. Han, H. Zhang, and G. Li, “Fast algorithms for sparse recovery with perturbed dictionary,” http://arxiv.org/abs/1111.6237. View at: Google Scholar
 I. Kyriakides, R. Pribic, H. Sar, and N. At, “GRID matching in Monte Carlo Bayesian compressive sensing,” in Proceedings of the 16th International Conference of Information Fusion (FUSION '13), pp. 2103–2109, Istanbul, Turkey, July 2013. View at: Google Scholar
 M. A. C. Tuncer and A. C. Gurbuz, “Analysis of unknown velocity and target off the grid problems in compressive sensing based subsurface imaging,” in Proceedings of the 36th IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '11), pp. 2880–2883, Prague, Czech Republic, May 2011. View at: Publisher Site  Google Scholar
 T. L. Hansen, M. A. Badiu, B. H. Fleury, and B. D. Rao, “A sparse Bayesian learning algorithm with dictionary parameter estimation,” in Proceedings of the IEEE 8th Sensor Array and Multichannel Signal Processing Workshop (SAM '14), pp. 385–388, A Coruña, Spain, June 2014. View at: Publisher Site  Google Scholar
 A. Fannjiang and H.C. Tseng, “Compressive radar with offgrid targets: a perturbation approach,” Inverse Problems, vol. 29, no. 5, Article ID 054008, 2013. View at: Publisher Site  Google Scholar
 A. C. Gurbuz, O. Teke, and O. Arikan, “Sparse groundpenetrating radar imaging method for offthegrid target problem,” Journal of Electronic Imaging, vol. 22, no. 2, Article ID 12317SS, 2013. View at: Publisher Site  Google Scholar
 G. Tang, B. N. Bhaskar, P. Shah, and B. Recht, “Compressed sensing off the grid,” IEEE Transactions on Information Theory, vol. 59, no. 11, pp. 7465–7490, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 C. Ekanadham, D. Tranchina, and E. P. Simoncelli, “Recovery of sparse translationinvariant signals with continuous basis pursuit,” IEEE Transactions on Signal Processing, vol. 59, no. 10, pp. 4735–4744, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 A. Fannjiang and W. Liao, “Coherencepatternguided compressive sensing with unresolved grids,” SIAM Journal on Imaging Sciences, vol. 5, no. 1, pp. 179–202, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 H. Zhu, G. B. Giannakis, and G. Leus, “Weighted and structured sparse total leastsquares for perturbed compressive sampling,” in Proceedings of the 36th IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '11), pp. 3792–3795, Prague, Czech Republic, May 2011. View at: Publisher Site  Google Scholar
 Z. Tan, P. Yang, and A. Nehorai, “Joint sparse recovery method for compressed sensing with structured dictionary mismatches,” IEEE Transactions on Signal Processing, vol. 62, no. 19, pp. 4997–5008, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Yang, L. Xie, and C. Zhang, “Offgrid direction of arrival estimation using sparse Bayesian inference,” IEEE Transactions on Signal Processing, vol. 61, no. 1, pp. 38–43, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 X. Zhou, H. Wang, Y. Cheng, Y. Qin, and X. Xu, “Offgrid radar coincidence imaging based on block sparse Bayesian learning,” in Proceedings of the IEEE Workshop on Signal Processing Systems (SiPS '15), pp. 1–5, Hangzhou, China, October 2015. View at: Google Scholar
 L. Hu, Z. Shi, J. Zhou, and Q. Fu, “Compressive highrangeresolution radar imaging using dynamic dictionaries,” IET Radar, Sonar & Navigation, vol. 7, no. 5, pp. 497–507, 2013. View at: Publisher Site  Google Scholar
 L. Hu, J. Zhou, Z. Shi, and Q. Fu, “A fast and accurate reconstruction algorithm for compressed sensing of complex sinusoids,” IEEE Transactions on Signal Processing, vol. 61, no. 22, pp. 5744–5754, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 D. G. Tzikas, A. C. Likas, and N. P. Galatsanos, “The variational approximation for Bayesian inference: life after the EM algorithm,” IEEE Signal Processing Magazine, vol. 25, no. 6, pp. 131–146, 2008. View at: Publisher Site  Google Scholar
 C. M. Bishop and M. E. Tipping, “Variational relevance vector machines,” in Proceedings of the 16th Conference on Uncertainty in Artificial Intelligence (UAI '00), pp. 46–53, San Francisco, Calif, USA, 2000. View at: Google Scholar
 T. Buchgraber and D. Shutin, “Distributed variational sparse Bayesian learning for sensor networks,” in Proceedings of the 22nd IEEE International Workshop on Machine Learning for Signal Processing (MLSP '12), pp. 1–6, Santander, Spain, September 2012. View at: Publisher Site  Google Scholar
 M. E. Tipping, “Sparse Bayesian learning and the relevance vector machine,” Journal of Machine Learning Research, vol. 1, no. 3, pp. 211–244, 2001. View at: Publisher Site  Google Scholar  MathSciNet
 C.Y. Chen and P. P. Vaidyanathan, “MIMO radar ambiguity properties and optimization using frequencyhopping waveforms,” IEEE Transactions on Signal Processing, vol. 56, no. 12, pp. 5926–5936, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 V. U. Reddy, S. Badrinath, and A. Srinivas, “Lowcomplexity design of frequencyhopping codes for MIMO radar for arbitrary doppler,” EURASIP Journal on Advances in Signal Processing, vol. 2010, Article ID 319065, 14 pages, 2010. View at: Publisher Site  Google Scholar
 S. Kwon, J. Wang, and B. Shim, “Multipath matching pursuit,” IEEE Transactions on Information Theory, vol. 60, no. 5, pp. 2986–3001, 2014. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2016 Xiaoli Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.