Analysis of Moving Object Imaging from Compressively Sensed SAR Data in the Presence of Dictionary Mismatch

We present compressed sensing (CS) synthetic aperture radar (SAR)moving target imaging in the presence of dictionarymismatch. Unlike existing work on CS SAR moving target imaging, we analyze the sensitivity of the imaging process to the mismatch and present an iterative scheme to cope with dictionary mismatch. We analyze and investigate the effects of mismatch in range and azimuth positions, as well as range velocity. The analysis reveals that the reconstruction error increases with the mismatch and range velocity mismatch is the major cause of error. Instead of using traditional Laplacian prior (LP), we use Gaussian-Bernoulli prior (GBP) for CS SAR imaging mismatch.The results show that the performance of GBP is much better than LP. We also provide the Cramer-Rao Bounds (CRB) that demonstrate theoretically the lowering of mean square error between actual and reconstructed result by using the GBP. We show that a combination of an upsampled dictionary and the GBP for reconstruction can deal with position mismatch effectively. We further present an iterative scheme to deal with the range velocity mismatch. Numerical and simulation examples demonstrate the accuracy of the analysis as well as the effectiveness of the proposed upsampling and iterative scheme.

Compressed sensing for synthetic aperture radar (SAR) is an active area of research for remote sensing.The use of CS based reconstruction can have an impact on the design of high resolution SAR systems as these systems encounter hardware design problems and require significant processing [12].CS has been applied for imaging of static objects in through-the-wall SAR imaging [13][14][15], tomographic SAR imaging [16][17][18], and SAR image formation with reduced data [19], where advantage is taken of the fact that the observed scenes are sparse.The static scenes may not always be sparse.The scenes containing a few strong intensity moving scatterers in a weak stationary background present an opportunity for CS application as they are inherently sparse.These moving targets suffer from position displacement and defocusing due to motion [20].The use of CS can help in reducing acquired data size as well as simultaneous motion parameter estimation imaging with reduced data.Sparsity can be further enhanced using clutter cancelation where the static parts of an observed scene are suppressed [21].
Compressed sensing for SAR moving object imaging has become an active area of research.References [22,23] apply CS for moving target parameter estimation by defining a dictionary based on the response of moving objects for different motion parameters.Both of these references use clutter cancelation to enhance sparsity.Reference [24] makes use of distributed CS applied to along-track interferometric SAR data for moving target imaging and shows that distributed CS can offer better performance with less samples compared

References
SAR moving target Prior Dictionary mismatch Remarks [6] N o G B P Y e s Superior performance of GBP over LP is shown by simulations [13,14] N o L P N o Apply CS for through-the-wall imaging [15] N o L P Y e s Apply CS for through-the-wall imaging.
Performance degradation due to pixel mismatch and wave propagation velocity shown by simulations [16][17][18] N o L P N o Apply CS for tomographic SAR imaging [19] N o L P N o Apply CS for focusing of static scenes [7] N o G B P Y e s Fast implementation of GBP reconstruction shows superior performance by simulations [21] N o L P N o Suggest clutter cancelation to enhance sparsity of a scene containing moving targets [22] Y e s L P N o Apply CS for motion estimation [23] Yes LP Yes Apply CS for motion estimation, performance degradation due to range velocity mismatch shown by simulations [24] Y e s L P N o Apply distributed CS for motion estimation [25] Yes LP Yes Apply CS for motion estimation; simulations show no performance degradation due to velocity mismatch [26,27] N o L P Y e s Performance degradation due to dictionary mismatch shown by simulations and theory to traditional CS.Reference [25] uses CS for moving target parameter estimation for mono-and multistatic SAR configurations and simulated data.These references show that CS can achieve imaging of moving objects as well as moving object parameter estimation when SAR data are sampled at a rate less than the traditional Nyquist sampling rate.
Compressed sensing reconstruction algorithms use a dictionary in which the reconstructed signal is assumed to be sparse.However, the dictionary in which the signal is actually sparse may be different and the resulting dictionary mismatch causes a performance degradation [26,27].In order to apply CS for practical applications, it is necessary to study the reconstruction performance degradation in the presence of dictionary mismatch.Reference [26] shows that dictionary mismatch can be seen equivalent to multiplicative noise.It also shows that reconstruction error increases linearly with mismatch.Reference [27] considers the effect of dictionary mismatch in CS reconstruction.It shows that, in case of using a Fourier dictionary, reconstruction performance degrades considerably when a mismatch exists.Due to this reason, it recommends examining the effects of mismatch on radar imaging.Reference [15] has shown performance degradation by means of imaging examples for static targets in the presence of mismatch in position and wave propagation velocity.The authors in [15] also state that they are extending the initial results presented in [28] for dealing with position mismatch in through-the-wall imaging.
According to the best of our knowledge, dictionary mismatch analysis has not been done theoretically for CS moving target SAR imaging in the presence of position and range velocity mismatch.A summary of the main features of the existing references is given in Table 1.It shows that, in the existing literature, the theoretical analysis of the effects of dictionary mismatch for moving target CS SAR imaging have not been carried out.Therefore, it remains an open problem.It further shows that a prior other than Laplacian prior (LP), for example, Gaussian-Bernoulli prior (GBP), for CS moving target imaging has not been used.Similarly, a theoretical analysis to show the advantage of the prior in dealing with dictionary mismatch is also missing.In [29], we have partially studied this problem and its effects for SAR and inverse SAR.We showed that dictionary generation using upsampled parameters is required to deal with errors arising due to mismatch in positions and range velocity.
The emphasis of this paper is to show the performance degradation in case of a target moving in the range direction.The dictionary mismatch arising due to discretization and dictionary size considerations causes performance degradation in terms of mean square error (MSE) between actual and reconstructed results, especially when there is a range velocity mismatch.We examine reasons for this degradation and also show theoretically and experimentally that using GBP for CS reconstruction compared to the traditionally used LP can compensate for some amount of mismatch.The motivation of using a different prior is to make use of extra information in improving reconstructed image quality as shown in [30].We propose to deal with CS SAR moving target imaging in the presence of dictionary mismatch due to positions and range velocity.The main contributions of this paper are as follows.
(1) We analyze dictionary mismatch and its effects theoretically, show MSE calculated from simulated SAR data for different types of mismatch in range and azimuth pixels as well as range velocity, and give parameter resolution limits for maintaining a reasonable level of reconstruction accuracy.We show that CS SAR moving target imaging is very sensitive to range velocity mismatch.
(2) We analyze the problem by means of Cramer-Rao Bounds (CRB) and show theoretically that reconstruction with Gaussian-Bernoulli prior (CSGBP) instead of traditional Laplacian prior (CSLP) can deal with some mismatch effectively.
(3) We present simulation results using CSGBP reconstruction and show that its use can lead to lower MSE, especially when the dictionary mismatch is small.This can be used to deal with position mismatch and reduce upsampling in positions that is required to counter mismatch effects.
(4) We also propose to reconstruct in the presence of range-velocity mismatch using an iterative scheme, where dictionaries with different range velocities are created efficiently.The contrast of the reconstructed result is maximized.
We would also like to point out that we deal specifically with the case of pulsed SAR.Any extension of dictionary mismatch effects and parameter resolution calculations to other types of SAR will need to take into account the difference in imaging mechanism; for example, in case of continuous wave SAR, it is known that range velocity creates a shift in the range direction, which is absent in pulsed SAR.Therefore, results for mismatch analysis and resolutions in range position and range velocity will need to take this additional shift into account.This paper is organized as follows.Section 2 presents the data model and formulation of moving target velocity estimation problem in case of CS SAR.Section 3 analyzes the effects of different kinds of dictionary mismatch, that is, range and azimuth positions and range velocity on CS SAR moving target imaging.Section 4 presents numerical and imaging examples to present the effects of dictionary mismatch in terms of MSE as well as the accuracy of the analysis and the effectiveness of the proposed method.Conclusions are given in Section 5.

System Model and Problem Formulation
In this paper,  denotes a scalar, x denotes a vector, and X denotes a matrix.We use X  and X  to denote conjugate transpose and transpose of X, respectively.The same notation is used for Greek characters; that is,  denotes a scalar,  denotes a vector, and Σ denotes a matrix.We use Σ  and Σ  to denote conjugate transpose and transpose of Σ, respectively.The function diag(x) represents a function that converts a vector x of size   × 1 into a diagonal matrix of size   ×   and det(X) represents the determinant of the matrix X.
Synthetic aperture radar consists of an antenna mounted on a moving platform [31].A pulsed SAR sends electromagnetic pulse () at a carrier frequency   and a chirp rate .The pulse length is denoted by   .This pulse is given as where and  =  1 , where  denotes measurement noise.Ψ 0 contains the response of each moving point in 1D form.This response for th moving point having th range velocity is given as [29] s   = [s   ( where The size of s   is   × 1. Ψ 0 is generated for initial velocity V  and final velocity V  for a total number of V  range velocities.The dictionary element corresponding to a velocity V 1   is as follows: The final dictionary Ψ 0 is stored in an   ×   matrix given as Due to undersampling, the problem of recovering  0 from s 0 becomes an underdetermined problem.We can solve this problem by including a-priori information for getting the solution; for example, select a solution such that the number of nonzero coefficients is the smallest.This can be expressed as follows: The number of nonzero coefficients is denoted by ‖ 0 ‖ 0 , known as  0 norm.However, this minimization problem is nonconvex, which means that finding a global solution is difficult or not guaranteed.In addition, it is computationally difficult to solve as it requires search over all possible combinations of the columns of ΦΨ 0 .To deal with these issues, we use  1 norm minimization.This minimization is a convex approximation of the  0 norm minimization if a property known as restricted isometric property (RIP) is satisfied.This property essentially means that the columns formed by the matrix ΦΨ 0 are sufficiently decorrelated with one another.The problem can be expressed as min      0    1 subject to s 0 = ΦΨ 0  0 .
In order to obtain a solution based on  1 norm minimization, we use Laplacian prior (LP) [32] as follows: If noise is Gaussian with variance  2  , the solution is obtained by where The solution can be written as that leads to Thus, by using LP, we include the  1 norm minimization in the solution.The parameter  gives weight to a priori sparse information.Equation ( 15) can be solved using different recovery methods, for example, linear programming and OMP.The reconstructed result σ0 is of size   × 1 and can be written as where each entry of σ0 shows the reconstructed reflectivity for each point in the scene for one velocity value; for example, σ1 0,1 represents the reflectivity for a point at position ( 1 ,  1 ) and having a velocity V 1  .The result can be rearranged into V  2D matrices, each having a size   ×   , to show the estimated reflectivities at different velocities for SAR.The matrices of size   ×   may also be summed to give a final focussed reconstructed result Σ0 , shown as follows: Dictionary mismatch can occur in the reconstruction process due to discretization of positions as well as range velocity; for example, instead of actual position of the scatterer (  ,   ) and velocity V   , the basis has elements corresponding to (  +Δ  ,   +Δ  ) and V   +ΔV   .Considering Ψ as the mismatched dictionary, (3) can be rewritten as and reconstruction using the mismatched dictionary Ψ causes the results to be decorrelated from σ0 shown as follows: Therefore, the effects of dictionary mismatch are related to the correlation between the mismatched and the original dictionary.In the next section, we examine the effects of this correlation on the reconstruction.Furthermore, we present solutions for the recovery of  0 that can be written as where Ψ 0 is the actual dictionary.We present solutions for calculation of σ0 using GBP that can reduce position mismatch effects, and propose an iterative scheme to recover σ0 in the presence of range velocity mismatch.

Effects of Position Mismatch.
We consider a chirp signal that is commonly used in imaging radars and show the effects of position mismatch on reconstruction.The reconstruction in the presence of mismatch depends upon the correlation between the original and the mismatched dictionaries, as given by (19).Therefore, any form of mismatch will cause erroneous results due to a correlation loss.This can be seen by taking the inner product of two chirp signals where  ∈ {0, Δ  , 2Δ  , . . ., ( Ñ−1)Δ  }.Δ  should be less than where  =  ÑΔ  .As the position mismatch increases, Δ  increases and, with the increase of Δ  , correlation decreases.Consequently, the amplitude of the reconstructed result will be reduced by a factor of |sinc(Δ  )|.The result of the correlation would be 0 when or Consequently, if the two chirp signals are displaced by Δ  with respect to each other, CS imaging will fail to reconstruct the correct position.An arbitrary element will be selected and cause a failure of the CS reconstruction model.Therefore, the smaller the distance between the dictionary elements, the smaller the mismatch and the better the reconstruction at the expense of larger dictionary size and higher number of computations.In effect, by making an oversampled dictionary, we can improve the reconstruction and this oversampling should be more than twice the sampling frequencies to reduce mismatch errors; that is, The implication of the above result will be studied in the next section.A moving target and the effects of mismatch in range and azimuth positions as well as range velocity on the moving target reconstruction are considered.[33]; that is: where Assuming that our dictionary is created with resolutions of   ,   , and 1 m/s in range position, azimuth position, and range velocity, respectively, the mismatch effects on a moving target can be divided into 3 categories as follows.
(i) A subpixel mismatch in range position represented as Δ  .This mismatch will lead to an equivalent shift of Δ  in the range position and an equivalent shift of −Δ  (V   sin   /) in the azimuth position.(ii) A subpixel mismatch in azimuth position represented as Δ  .This mismatch will lead to an equivalent shift of Δ  (V  sin   /) in the range position and an equivalent shift of −Δ  in the azimuth position.
(iii) A fraction of m/s mismatch in range velocity represented as ΔV   .This mismatch will lead to an equivalent shift of   (ΔV   /) sin   in the range position and an equivalent shift of −  (ΔV   /) sin   in the azimuth position.As an example, if a point in acquired raw data is at position (  ,   ) moving with a velocity V   + ΔV   , and the dictionary contains elements with velocity V   , the reconstructed estimate will be a point at position (  +   (ΔV   /) sin   ,   −   (ΔV   /) sin   ) instead of the true position of (  ,   ).As   is large, the effect on azimuth position will be more evident even when range velocity mismatch is small.The mismatch effects due to Δ  , Δ  , and ΔV   are summarized in Table 2.

Effects on
, where  1 and  1 are the pixel positions, can be summarized as follows.
(i) A mismatch of Δ 1 will cause a shift of Δ 1 in range position in the reconstructed result.The shift in azimuth position −Δ 1 (ΔV 1  /) sin  1 can be neglected as it is small due to the presence of  in the denominator.The result will be σ0 =  1  1 sinc ( leading to a loss of amplitude.
(ii) A mismatch of Δ 1 will cause a shift of Δ 1 in azimuth position in the reconstructed result.The shift in azimuth position Δ 1 (ΔV 1  /) sin  1 can be neglected due to the presence of  in the denominator.The result will be σ0 =  1  1 sinc ( leading to a loss of amplitude. (iii) A range velocity mismatch ΔV 1  causes a large shift in azimuth from the true position, given as The shift in the range position  1 (ΔV 1  /) sin  1 can be neglected.However, the shift in azimuth position cannot be neglected due to the presence of  1 in the numerator that is of the order of 10 3 or higher.It can be further divided into 2 parts as follows.
The reconstructed result will be σ0 =  1  1 sinc ( leading to a loss of amplitude and azimuth position shift. (iv) In order to avoid the loss in amplitude as well as azimuth mispositioning of the reconstructed result, the dictionary can be created with higher parameter resolution.The dictionary resolutions in range and azimuth positions and range velocity are such that any mismatch does not lead to a misselection of elements.This can be achieved if the dictionary resolutions are less than half the pixel sizes.This ensures that a correct pixel positions is selected.These criteria can be expressed as follows for range and azimuth positions: In case of velocity, the shift in azimuth position caused by range velocity mismatch should be less than half the pixel size; that is, or As this shift is large for a larger value of   , we choose the farthest slant-range distance  max to get a conservative estimate as follows: where the angle corresponding to  max is  max .This leads to The limit given by ( 32) is also applicable for compensating intrapixel displacements due to velocity mismatch.Please note that due to the presence of  max in the denominator, ΔV  is very small, which means that the dictionary needs to be created with very closely spaced velocity values.When there is a moving scene consisting of a number of points given as the reconstructed result in the presence of mismatch is as follows: The three sinc functions represent a loss in estimated amplitude due to the mismatch, whereas the second term in the delta function represents a pixel-level shift.
International Journal of Antennas and Propagation 7 3.4.Using CSGBP to Improve Performance in the Presence of Mismatch.In order to avoid errors due to dictionary mismatch, the dictionary needs to be created with upsampled positions and range velocity parameters.This high upsampling may not be feasible due to limited storage and computational complexity.We propose to reduce this high upsampling requirement by using a different prior as well as an iterative scheme.The chosen prior is GBP given as [6] where   0, is the th element of  0 moving with th velocity.The main motivation of using this prior is to utilize a priori information about sparsity and signal strength for image reconstruction. 0 can be assumed as -sparse that is represented by the probability of active elements  = /  in  0 .The prior assumes that the probability of active elements, that is, an entry of  0 being nonzero, is given by  and these active elements are represented by a Gaussian distribution with mean   and variance  2  .The probability of an inactive element is given by 1 − .The solution to recover  0 from s 0 for the prior can be obtained by rewriting (3) as follows: where The th entry of z 0 is 1 if the corresponding entry in  0 is 1.In this case,  0 can be recovered from s 0 in two steps as follows.
(1) The 1st step is the solution to the following problem: where and (s 0 | z 0 ) is given on the next page.For the sake of convenience, we define D = ΦΨ 0 and the covariance matrix is given as R =  2  D diag(z 0 )D  +  2  I.The solution can be further written as (2) The solution ẑ0 obtained from the 1st step is used to recover estimate of  0 by using least squares solution given as Furthermore, this model is suitable for man-made moving scatterers as they may be represented as consisting of a coherent mean part and variation of reflectivities can be represented by an incoherent part represented as variance; that is,  0 ∼ N(  ,  2  ).In addition, noise  can be assumed to be zero-mean Gaussian with variance  2  ; that is,  ∼ N(0,  2  ).This CSGBP model can be solved using the algorithms in [6] or [7].In [6], the raw data is correlated with each column of the matrix D, and the presence or the absence of an element is decided by hypothesis testing.This testing is based on the assumption that the signal is distributed according to the GBP and the noise has Gaussian distribution.In [7], an efficient method is proposed for finding a combination of active and inactive elements.

Analysis of CSGBP and CSLP Performance in the Presence of Dictionary Mismatch Using Cramer-Rao Bounds.
To show theoretically the advantage gained by using CSGBP reconstruction model given in (17) over CSLP model in (7), CRB of the vector σ0 estimated from data vector s 0 is calculated as the inverse of Fisher information matrix (FIM) J.We consider Φ to be identity matrix in (3) for the sake of convenience.The FIM bounds the estimation error in the following form: Φ is assumed to be an identity matrix for the sake of convenience.J is decomposed into two parts [34]; and J P represents prior information matrix whose individual elements are given as Making use of the explanation given in [35] and smooth approximation; that is, , the FIM is given as for the case where CSLP is used.When CSGBP is used, the FIM is As (50) contains more information compared to (49), J CSGBP in (50) will be larger and hence the estimation error will be lower that shows the improvement in performance.In case of a dictionary mismatch, using ( 18) and ( 19), ( 47) and ( 48) become and J P represents prior information matrix whose individual elements are given as Equation ( 49) becomes and ( 50) is rewritten as ) . ( When no dictionary mismatch is present, Ψ  Ψ 0 has a maximum value along the diagonal elements.In case of mismatch, the diagonal elements of Ψ  Ψ 0 decrease.Subsequently, J decreases leading to an increase in estimated error.It can be inferred that, due to the prior information in (54), the increase of estimated error in the presence of dictionary mismatch is less when CSGBP is used.This can be seen in Figure 1, where an identity matrix of size 512 × 512 pixels is used as Ψ 0 .Ψ is a mismatched basis that is decorrelated with Ψ 0 in varying proportions as follows: where  is the degree of correlation and the measure 1− can be seen equivalent to dictionary mismatch proportion.MSE is calculated using the expression It can be seen that using the model given in (39) lowers MSE that can help in countering effects of decorrelation arising due to dictionary mismatch.

Dealing with Range Velocity Mismatch Using Iterative CSGBP.
As outlined in the previous section, CSGBP can compensate for some mismatch, which can help in reducing upsampling requirements.However, it is still not possible to deal with range velocity mismatch using only CSGBP.In general, CS SAR moving target imaging is very sensitive to range velocity mismatch.To avoid any error due to range velocity mismatch, the dictionary should be created with a very high resolution in range velocity; for example, for typical SAR configurations, this resolution can be of the order of 0.01 m/s.Such a high upsampling requirement is not feasible due to limited memory requirements and very high computations.
In this section, we propose to compensate for velocity mismatch by creating a dictionary iteratively, with range velocities varying at each iteration.In order to reduce the computational time, we make use of the following approximation to (6): We can make use of this approximation to create s+1  (  ,   ) from s  (  ,   ) as follows: where ΔV  , = V  , − V +1 , .This allows us to create a dictionary with varying mismatch iteratively using already computed dictionaries.Using the approximation, we propose the following scheme to reconstruct SAR image in the presence of dictionary mismatch.
(1) Create a dictionary with range and azimuth positions at a subpixel resolution.This resolution is chosen so as to meet the upsampling requirements given by ( 31) and (32).We chose an upsampling factor of 4 in position, which means that the maximum mismatch that can occur is 1/8 of the pixel size.This process is carried out only once.
(2) Carry out CSGBP reconstruction using the dictionary created in Step 1. Due to the upsampling chosen in the range and azimuth directions, and, due to the fact that the range velocity mismatch does not affect the range position, the result contains correct range position as well as range velocity.There will be azimuth position displacements due to range velocity that will be compensated in the next steps.
(3) For each set of reconstructed   points belonging to the same range velocity V   , regenerate new dictionary elements at the selected range positions using (58) and a velocity increment of ⌊Δ   ,   ⌋. (4) Step 3 is repeated by incrementing the velocity in steps of ⌊Δ   ,   ⌋, until the reconstructed image is judged to be of the best quality for the   points.As a quality measure, contrast of the reconstructed vector is calculated as follows: where ⟨⋅⟩ is the averaging operator.
(5) Steps 3 and 4 are repeated for each velocity in the dictionary where moving points were detected in Step 2.

Numerical and Imaging Results
This section presents numerical and imaging results.We give examples with MSE calculated for different amounts of mismatch in range, azimuth, and range velocity for SAR data, followed by imaging examples.

Numerical Results.
The simulation parameters for SAR data are given in Table 3.A scene of size 50 m × 50 m or 12 × 70 pixels in range and azimuth directions is considered.Raw data corresponding to multiple points are simulated and 5% of range data are retained.Positions and amplitudes of these points are chosen randomly, whereas ground-range velocities are chosen randomly from a set of 7 velocities: {0, 3, 4, 5, −3, −4, −5} m/s.Performance in terms of dictionary mismatch is compared.For this purpose, data are generated using a dictionary Ψ 0 and CS reconstruction is carried out using a mismatched dictionary Ψ.The mismatch has a value of 0.01, followed by values from 0.1 to 0.7 with a step-size of 0.1.For range and azimuth pixels, the mismatch unit is pixel size, whereas, for range velocity, it is m/s.A series of simulation is carried out at a signal-to-clutter ratio (SCR) of 20 dB with randomly chosen positions and velocities of the moving targets.Reconstruction is carried out using CSLP and CSGBP and the resulting MSE between the original points and the reconstructed points are shown in Figure 2. MSE is calculated as follows: Four main parameters are used in CSGBP reconstruction: ,   ,   , and   , which are initially estimated by using a-priori information.The value of  is decided according to the ratio of supposed active scatterers to total number of scatterers present in the data, whereas the values of   ,   , and   are chosen based on SCR.They are then refined by trial and error to get the best results.In general, higher than required values of ,   , and   help in producing weak scatterers but lead to more side lobes, whereas a higher value of   suppresses weak scatterers.From Figure 2, the following observations can be made.
(i) In general, reasonable reconstruction is obtained when the effect of basis mismatch is less than 1/3 of a pixel size.
(ii) MSE is less in case of no range and azimuth pixel mismatch using CSGBP.Similarly, for a small mismatch in range and azimuth directions, the MSE level in case of CSGBP based reconstruction is less.Specifically, it can be remarked that although for the velocity mismatch, MSE increases when velocity mismatch reaches 0.1 m/s; however, in case of range and azimuth pixels mismatch, MSE is very small as long as pixel mismatch stays less than 0.3 of the pixel size.Thus, CSGBP can be used for better reconstruction and reduction of the dictionary size in practical scenarios, compared to CSLP based reconstruction, where the MSE is higher even in case of no dictionary mismatch.
(iii) MSE for range velocity is high using both methods.After the mismatch of 0.1 m/s, CSLP seems to give slightly lower MSE.The reason may be that CSGBP gives higher number of side lobes.Further simulations for the values of mismatch ranging from 0.01 to 0.1 in a step size of 0.01 m/s are shown in Figure 3.
It can be seen that MSE using CSGBP is still smaller than that using CSLP.The reason for not reporting any ill effects of mismatch in velocity in [25] may be that the amount of mismatch considered is small for the configuration that was studied.There are two types of moving targets that are considered in [25], a slow one and a fast one.The former target has a range velocity of 2.35 m/s, whereas the latter target has a range velocity of 28.15 m/s.The range velocity mismatch for the slow target is 0.85 m/s, whereas, for the fast target, it is 0.45 m/s.The amount of mismatch is small to have any effect on the reconstruction for the particular case.This can be seen from reconstruction results in Figure 5 of [25] that shows focussing assuming no motion.The slow object, despite having a mismatch of 2.35 m/s in the range direction, is still focussed at the same position.Our results show theoretically as well as experimentally that a mismatch in velocity can have a serious impact on reconstruction.
(iv) The error increases gradually for position mismatch but increases very rapidly for range velocity mismatch.The reason is that, in case of range velocity mismatch, a large shift arises in azimuth direction.This is due to the reason that   is of the order of 10  for example, for a velocity mismatch of 0.05 m/s,   = 2595 m, and   = 40 ∘ , there is a single-pixel shift between the original and reconstructed position.Thus, the reconstruction result will contain azimuth pixels shifted according to the mismatch, which leads to a sudden increase in MSE.As there is a total misalignment between actual and estimation positions, MSE rises and stays at a roughly constant maximum level.This is further demonstrated in Figure 4, where a reconstructed scene contains a single pixel shift with Amplitude Azimuth bins (a.u.) Results with mismatch Results without mismatch respect to the actual position.The reason is that there is a mismatch in range velocity of 0.05 m/s.A loss of amplitude and side lobes can be seen.
(v) The error in azimuth is more than that in range position in general, especially using CSLP.The reason is that there are more than one combination of   ,   , and   that lead to closely resembling values of   and   in (26); for example, for the parameters given in Table 3, we can see that   and   of 2594.5 m and

Results with mismatch
Results without mismatch  .2128m and V  = 5 m/s.Thus, it is possible that a dictionary mismatch will lead to selection of dictionary elements and subsequently, side lobes that are not in the immediate neighborhood.This is demonstrated in Figure 5, where there are 4 points at different azimuth positions having a velocity of 3 m/s.The mismatch is 0.1, 0.3, 0.5, and 0.7 of a pixel size.When CSLP is used to carry out reconstruction, only a single point is identified with a velocity of 3 m/s.This is shown in Figure 5(a).Two of the points are detected at shifted azimuth positions with a velocity of 4 m/s, as shown in Figure 5(b).The fourth point is not identified at all.In case of reconstruction using CSGBP, all of the four points are identified correctly as shown in Figure 5(c), albeit with higher side lobes.This also demonstrates the advantage offered by CSGBP by identifying correct positions and velocity even in the presence of pixel mismatch.Reconstruction is carried out in the presence of a mismatch of 1/8 of a pixel size in range and azimuth.Results using CSLP are shown in Figure 6(b), which shows that the point at (7,50) is not reconstructed correctly.CSGBP results shown in Figure 6(c) indicate that all the points are correctly reconstructed.This demonstrates the superior performance of CSGBP.Furthermore, results using CSGBP show side lobes in the vicinity of actual positions, whereas, in case of CSGBP, the side lobes appear at positions that are not in the vicinity of actual positions.A further example is shown with a mismatch of 0.4 m/s in range velocity.The original scene is shown in Figure 7(a), where there are closely spaced scatterers roughly in the middle of the scene.They have a velocity of 4.4 m/s, whereas the closest velocity in the dictionary is 4 m/s.Reconstruction using both CSLP and CSGBP shows shifted results due to the mismatch.Furthermore, results obtained using CSLP were obtained at a velocity of 5 m/s.Result obtained using iterative CSGBP is shown in Figure 7(d), where the points are located at their correct positions.The velocity in the dictionary is increased iteratively with a step size of 0.05 m/s, until the highest contrast is achieved.A plot of contrast with velocity is shown in Figure 8, where it can be seen that the contrast is the highest when the velocity in the dictionary matches the actual velocity.This shows that creating dictionary elements iteratively and using contrast to measure quality are effective methods for dealing with CS moving target imaging in the presence of range velocity mismatch.Another example is shown with a scene in Figure 9(a).The points are at positions of (2, 5), (5, 45), (6,34), (2, 70), (10,15), and (7,65).The point at (2, 5) has a velocity of −4.9 m/s, the point at (   A zi m ut h bi ns R a n g e b i n s 70 60 50 40 30 20   The result obtained by maximizing the contrast for the point at (6,34) is shown in Figure 9(d), where the points moving at 3.3 m/s are focussed at their true position.Some side lobes can be seen.Similarly, the point moving at −4.9 m/s is shown correctly focussed in Figure 9(e).Final result obtained using the calculated velocities is shown in Figure 9(f), where all the points are focussed at their true positions.Some side lobes can be observed.

Conclusion
In this paper, we studied compressed sensing (CS) synthetic aperture radar (SAR) moving target imaging in the presence of dictionary mismatch.We analyzed the sensitivity of the imaging process to range pixel, azimuth pixel, and range velocity mismatches.The mismatch analysis shows that the reconstruction error increases with mismatch and especially increases very rapidly in the presence of range velocity mismatch.Unlike existing references, we show that using a Gaussian-Bernoulli prior compared to the traditionally used Laplacian prior offers advantage in CS SAR imaging for dealing with small mismatch.This advantage is apparent in dealing with positions mismatch.We calculated Cramer-Rao bounds that demonstrate theoretically the lowering of mean square error between actual and reconstructed result by using the GBP.We show that creating an upsampled dictionary and using the GBP for reconstruction can deal with position mismatch.We also presented an iterative scheme to deal with the range velocity mismatch where dictionary elements are created efficiently.CS reconstruction is carried out at each iteration until the image contrast is maximized for each velocity.Numerical and imaging examples confirm the analysis and the effectiveness of the proposed upsampling and iterative scheme.

Figure 4 :
Figure 4: Effects of range velocity mismatch using CSGBP.A singlepixel shift and some side lobes can be seen.
Reconstruction in the presence of azimuth pixel mismatch using CSLP Azimuth bins (a.u.)(a)Effects of azimuth pixel mismatch using CSLP.Only one position is identified correctly with a velocity of 3 m/sReconstruction in the presence of azimuth pixel mismatch using CSLP Azimuth bins (a.u.)(b)Effects of azimuth pixel mismatch using CSLP.The scatterers are identified at a velocity of 4 m/s, instead of a velocity of 3 m/s.Scatterers are also shifted from their true positions Azimuth bins (a.u.)Reconstruction in the presence of azimuth pixel mismatch using CSGBP (c) Effects of azimuth pixel mismatch using CSGBP.The positions are identified correctly with the presence of higher side lobes

Figure 5 :
Figure 5: Effects of azimuth pixel mismatch using CSLP and CSGBP.

00Figure 6 :
Figure 6: Reconstruction in the presence of range and azimuth pixels mismatch using CSLP and CSGBP.

00 4 (
Reconstructed scene using CSLP (b) Reconstruction using CSLP.The points are shifted in azimuth direction Amplitude A zi m ut h bi ns R a n g e b i n s Reconstructed scene using CSGBP (c) Reconstruction using CSGBP.The points are shifted in azimuth direction Amplitude A zi m ut h bi ns R a n g e b i n s d) Reconstruction using iterative CSGBP.The points are at their original positions.A few side lobes can be seen.

Table 1 :
Comparison of existing references.
2 ,  3 , . . .,    .The signals are reflected from each scatterer in the observed scene.Let  0 be a sparse vector of size   × 1 that contains reflectivities for each point in the scene having different motion parameters.Ψ 0 is an   ×   matrix in which the signal is actually sparse and contains response of moving targets for every point in the scene with each considered motion parameter.Let Φ be a sampling matrix of size ×  , where  <   .

Table 2 :
Effect of mismatch on range and azimuth positions. (ΔV   /) sin   Azimuth Δ  (ΔV   /) sin   −Δ  Velocity   (ΔV   /) sin   −  (ΔV   /) sin   3.2.Effects of Range Position, Azimuth Position, and Range Velocity Mismatch on Reconstruction of a Moving Target.First, we consider the equivalent static position of a moving point.A moving point at an initial position of   ,   and having a velocity of V   can be equivalently seen as a static point with coordinates   and   and rotated with an angle Θ
Range time   : Central frequency : Slant-range positions : Azimuth positions  Speed of light : Azimuth time   : Incidence angle at range   , equal to cos −1 (ℎ/  ) b e ro fr a n g ep i x e l si nt h es c e n e   : Number of azimuth pixels in the scene   : Number of range pixels in raw data   : Number of azimuth pixels in raw data  V : Numberofrangevelocities   :  ×  V :   ×     :   ×   s 0 : Raw data from all the points in the scene arranged in 1D form s   : R a wd a t af o rth point moving with th velocity, arranged in 1D form s  (  ,   ): Raw data element for range time   and azimuth time      (  ): Radar-target distance for th point moving with th velocity Ψ :  : 4  / : Chirp rate : Sensor velocity ℎ: Antenna height   : Pulse length : Ground-range positions :