Sparse signal processing on estimation grid with constant information distance applied in radar

Radar obtains its parameters on a grid whose design supports resolution of underlying radar processing. Existing radar exploits a regular grid although the resolution changes with stronger echoes at shorter ranges. We compute the radar resolution from the intrinsic geometrical structure of data models that is characterized in terms of the Fisher information metric. Based on the information-based approach, we design an estimation grid whose cells have a constant Fisher information distance. In addition, we explore how this information-based grid can suit radar processing in practice and propose information-based processing on such an irregular estimation grid by applying the sparse signal processing from compressive sensing. Accordingly, the grid was adjusted to the sensing incoherence needed in sparse signal processing by setting a lower bound for the cell size. Our approach enables an adaptive estimation grid that can be adjusted with respect to the available resolution, the desired sensing incoherence, available computational power, and required operational priorities. The information-based design and processing are illustrated in a one-dimensional case of range estimation.


Introduction
A radar grid is designed to support the resolution of underlying radar processing.Radar resolution is the minimum distance between two objects that radar still can resolve.Besides the effective sensing bandwidth, the signal-to-noise ratio (SNR) is also crucial in the ability to resolve close objects (e.g., [1]).In existing radar, the cell size is constant, i.e., the estimation grid is regular because of the real-time computations done by fast Fourier transform (FFT).However, the resolution is not constant because it increases with stronger echoes at closer ranges as well known from the radar equation (e.g., [2]).
In this paper, we continue exploring a practical combination of information geometry (IG) and compressive sensing (CS) in a basic radar example of range estimation (started in [3,4]).IG is applied to design a resolution *Correspondence: radmila.pribic@nl.thalesgroup.com 2 Sensors Advanced Developments Delft, Thales Nederland B.V., Delft 2628 XH, The Netherlands Full list of author information is available at the end of the article cell in a radar system as related not only to the sensing bandwidth but rather to information distances on the statistical manifold of the Fisher information matrix of a data model.In a range-only case, this implies that closer ranges would have smaller cells.Our resolution analysis differs from previous work (e.g., [5,6]) as we include the range-dependent amplitude and use raw radar measurements throughout the whole analysis.In addition, we also explore how such an irregular grid suits radar processing and propose applying sparse signal processing (SSP) that is nowadays a major part of CS.
SSP designates a sparse model-based refinement of existing signal processing.In radar, matched filtering remains vital within SSP but followed by the 1 -norm iterative refinements promoting the sparsity [7].SSP is a major part in the back end of a sensor with CS, while its front end enables compressive acquisition of measurements.CS is optimized to the information in received measurements through a specific view based on the two assumptions: sparsity of processing results and the sensing incoherence (e.g., [8,9]).
IG raises a new approach to stochastic signal processing by treating the stochastic inferences as structures in differential geometry (e.g., [5,10]).The intrinsic geometrical structure of measurement models is conveniently characterized in terms of the Fisher information metric.Accordingly, resolution of sensors can be based on information distances on such statistical manifolds.
Both IG and CS have a potential to improve radar processing because of the emphasis on the information density in received data instead of the much larger sensing bandwidth.Accordingly, conventional processing can be improved if the demands of data acquisition and signal processing are optimized to the information content in radar measurements.
In our work, we have been exploring the innovative IG-CS combination that offers the new flexibility in radar design.In [3], we had focussed only on design of the IG-based grid.In [4], we proceeded by focusing on a combination of the IG-based grid design with SSP.The emphasis was on the effects of the IG-based estimation grid to the sensing incoherence (typical for SSP or CS).In this journal paper, we provide a more detailed analysis, express the Fisher information metric in more generic terms of the waveform specific constant, and indicate the grid adaptivity leading to a coherence-adjusted irregular grid.Moreover, we give numerical results to illustrate the theoretical analysis in a practical radar example.
In Section 2, typical radar processing is recalled first followed by the information-based approaches of IG and CS.In Section 3, the analysis is illustrated with simulations results in a range-only case.In Section 4, conclusions are drawn, and future work is indicated.

Typical radar processing
Raw radar measurements gathered in a vector y are described by nonrandom signals s in receiver (thermal) noise z in a model: y = s + z.If s is described as in CS, by a linear model: s = Ax, it writes as by a sensing matrix A, a sparse radar profile x, signals Ax , and complex Gaussian noise z with zero mean and equal variances γ , p(z|γ ) ∝ exp −|z| 2 /γ .The sensing matrix A contains a radar signal model with a desired profile x that is known from the physics.For example, in range, columns of A are time-delayed replicas of the transmitted waveform that is usually a linear chirp (linear frequency modulation (LFM)) because of the optimal processing gain.Typical radar processing starts with matched filtering (MF) because of the optimal SNR of a single target in radar echoes (e.g., [2]).The MF profile x MF , x MF = A H y (H denoting hermitian) is computed using FFT and also decoupled for angle(s), range, and Doppler.Accordingly, both observation and estimation grids are regular.The time resolution of x MF is proportional to the signal bandwidth B and √ SNR, while a grid cell is usually 1/B large (as given by the Nyquist sampling of complex signals).

CS and sparse signal processing
In CS, x from Equation 1 is assumed sparse (or compressible).In a Bayesian framework, it can be formalized by a multivariate Laplace prior p(x|λ), p(x|λ) ∝ exp (−λ ||x|| 1 ).The maximum a posteriori (MAP) estimator of x, written as gives the usual SSP from CS with the 1 norm ||x|| 1 promoting the sparsity and the 2 norm y − Ax for minimizing the noise, together with a threshold h that balances between the two tasks (e.g., [8,11]).An undetermined system can be solved by SSP from Equation 2, i.e., M measurements in y can be enough for N outputs in x, M < N, because of the sparsity, i.e., only K nonzeros in x, K < M, and incoherence of A, i.e., a low inner product between its different columns.There are three measures for the incoherence of a matrix: the mutual coherence, the restricted isometry property (RIP), and null space property (NSP); the last two being more difficult to compute (e.g.[8]).
The mutual coherence c(A) of a matrix A is given by the largest inner product of columns of A, denoted as a n , with n = 1, . . ., N and normalized, i.e.,|a i | = 1, as: In radar, the straight physical nature of a sensing matrix A suits CS and the incoherence well.The mutual coherence c(A) can be viewed as the largest off-diagonal value of A H A, i.e., the highest sidelobe (e.g., [12]).The sparsity of x is bounded by c(A), e.g., |x| 0 < 1 + 1/c(A), where |x| 0 is the (pseudo) norm 0 , i.e., a number of nonzeros in x.The number of measurements M needed to solve x is also related to c(A) as The mutual coherence c(A) is also connected to the RIP of a sensing matrix A in a linear model.The matrix , where C N K indicates a set of all complex vectors of size N and sparsity K (e.g., [8]).The RIP(K , ε) is satisfied for any K < 1 + 1/c(A).This relates the two measures of the sensing quality by a matrix A. http://asp.eurasipjournals.com/content/2014/1/78Practical CS in radar prefers stochastic SSP when treating noise, prior knowledge on signals or their data acquisition, and when providing results ( [13]).In a way, stochastic SSP is giving a fresh boost to radar initiated in the 1950s ( [14]).Accordingly, we use the algorithm complex fast Laplace (CFL, [15] and [13] for complex signals).In CFL, the prior p(x|λ) is built from a complex Gaussian prior for x and a hyperprior for the variance of x.CFL refines actually the MF profile x MF in a number of iterations by selecting significant elements based on increase in the assessed posterior in each element of x.For the optimal processing gain, measurements y are to be gathered over the whole observation time and space, supporting the SSP of a profile x over all its parameters: range, doppler and angle(s).The matrix A as well as the SSP size becomes huge but well arranged, what makes SSP even achievable in real time (as in [13], CFL implemented in range-Doppler SSP on GPU).
When SSP is applied on an irregular grid, the effects on the incoherence of the sensing matrix A should be first considered.When the cell size decreases, the columns of sensing matrix A become more coherent because they come closer to each other.This limits a minimum grid cell size, as the incoherence is a requirement for SSP.Thus, IG gives us possible grid cells while SSP poses the lower bound to the cell size.
In addition, the computational demands of SSP applied on an irregular grid should also be considered.The number of cells can become equal or even smaller as optimized to the possible resolution, desired incoherence, and available processing power.Finally, operational needs will also benefit from the grid adaptivity.For example, closer ranges can have a higher priority in the radar operation.
How SSP on irregular grid can suit radar processing, and make it even more efficient, is illustrated in a range-only case in a case study in Section 3. Note that with the same number of cells, closer targets are resolved.

Information geometry
IG is the study of manifolds in the parameter space of probability distributions, using the tools of differential geometry [5].These spaces are generally non-Euclidean, which basically implies that the scalar product of two vectors x and y, is redefined by a metric tensor g ij (e.g.[16]), i.e., x, y Here, metric g ij may depend on the location θ on the manifold, so g ij → g ij ( θ).A consequence of the metric tensor is that the actual length of curve is different from the length in Euclidean space.For an arbitrary, differentiable curve ψ(t) = (ψ 1 (t) . . .ψ n (t)), that is on the manifold equipped with metric g, denoted as ( , g), the length is given by When ψ(t) is (locally) the shortest paths between two points, ψ(t) is called a geodesic of ( , g): this is the extension of the notion of a straight line to non-Euclidean spaces.
A clear example of a geodesic in a non-Euclidean space is the shortest path on the spherical surface.There are no straight lines here: the shortest path between two points on the spherical surface is the shorter great circle arc (see Figure 1).
In IG, the metric tensor is defined by the Fisher information matrix.The Fisher information matrix is defined for some likelihood function p(y| θ) and parameter vector θ as the expectation value Figure 1 A geodesic on a spherical surface, a shorter great circle arc.http://asp.eurasipjournals.com/content/2014/1/78 In this paper, we will study the basic case for range only.For this one-dimensional case, the Fisher information metric reduces to Distances on this one-dimensional statistical manifold can simply be calculated from ds 2 = G(τ )dτ 2 (e.g., [5]).We will derive an expression for the Fisher information metric G(τ ) for arbitrary waveforms in the next part.

Information distance cell design
We derive the Fisher information metric G(τ ) to design and construct grid cells with a constant information distance as they support radar resolution.In the derivation of G(τ ), we closely follow the approach from [6].How we construct the grid cells differs from the previous works (e.g., [5] and [6]) as we directly incorporate the delaydependent amplitude and use raw radar measurements in the data model.This is incorporated by letting the signal amplitude a be a function of delay, i.e., a → a(τ ).The signal amplitude a (that is independent of τ ) is our reference case.The received signal is modeled in continuous form as where s(t − τ ) is a replica of the transmitted signal delayed by τ and z(t) is the zero-mean, complex Gaussian noise with variance γ .The measurement y is complex Gaussian, defined by The probability density function p(y|τ ) is given by where τ is the parameter of interest.The Fisher information metric is given by Equation 9, or In order to avoid certain convergence problems [17], we divide the likelihood function to obtain a likelihood ratio (y|τ ), which is again a likelihood function: Because is independent of τ , we know that so we can continue to find the Fisher information metric based on (y|τ ).We then take the natural logarithm of the likelihood function (y|τ ), Differentiating twice with respect to τ results in Taking the expectation value and multiplying by −1 gives the Fisher information metric Note that for the reference case, this result reduces to In previous work, the result G 0 was derived and the substitution a → a(τ ) was made to incorporate the delay into the Fisher information.The difference from current work is that the delay dependence in the Fisher information metric G(τ ) was derived by incorporating the fundamental radar equation in the model of raw data.

Delay-dependent amplitude and phase-modulated waveforms
From the fundamental radar equation (e.g., [18]), we know that the received power scales with τ −4 .This means that the amplitude of the received signal can be written as with b some complex parameter that may include gain, RCS, etc.The delay-dependent Fisher information metric G(τ ) from Equation 19 can then be written as In order to apply the previous equation to a test case, we define the unit energy replica of the transmitted waveform with square envelope and phase modulation (PM) ϕ(t) as where T p is the pulse duration and ϕ(t) is a real, differentiable function on the interval [ 0, T p ] that represents the phase of the signal.Note that using this construction, the frequency and phase modulation is controlled by ϕ(t) and the amplitude by a(τ ).The first derivative is calculated as so G 0 and G(τ ) may be written as and where the constant C is specific for a waveform and its PM, defined as Information distances on both manifolds are calculated intrinsically different.For the reference case, it is a Euclidean metric, i.e., while the distance on the new manifold, with delaydependent amplitude control is calculated as With these results, the grid cells will be constructed in the next part.

Choice of the cell size
Using the expression in Equation 29, we have to construct resolution cells.For this, we need to determine the information distance of one grid cell δ.We start with the calculation of the total information distance of the listening interval, i.e. the time interval between T p and pulse repetition time (PRT), d(T p , PRT) and the choice of an appropriate number of grid cells N cells .Combining these, the size of the information grid cell is determined through Note that in the current analysis, the relative size of the grid cell with respect to the total PRT δ/d(T p , PRT) is independent of parameters b and γ because they are held constant.
An alternative method is to (experimentally) determine at what distance the conventional resolution (from the reference grid) is 'correct'.The conventional resolution cell at this position m (and thus given on the interval [τ m , τ m +1]) determines the size of the information resolution cell Another possibility is to use the method from [5].This paper uses information distance in a Monte Carlo method for determining whether there are one or two targets present in a small area.The critical distance for the model radar in distinguishing the two targets with some empirical probability of error is defined as the resolution cell.
When the value for δ is determined, the new resolution grid is constructed by starting at τ 1 = T p and iterating along the PRT Note that when Equation 30 is used, this results in an integer number of grid cells on the listening interval.
Processing on the grid will be done with SSP.For this, the coherence of the sensing matrix c(A) must be low.This is not the case when the grid cells are too close together, which leads to a lower bound in grid cell size.This coherence-adjusted grid is designed for the practical case in Section 3.

LFM pulse
We apply the above results to the test case of a linear chirp (LFM pulse) of bandwidth B. The phase is given by Now, constant C for this waveform is calculated from Equation 27 as This gives expressions for both metrics, and The Fisher information metrics G 0 and G(τ ) are shown in Figure 2a.
Distances on the statistical manifold of the delaydependent case are calculated as By combining the last result with Equations 30 and 32, the cells were constructed and their size is shown in Figure 2b.

Numerical results
The information-based design and processing with the adaptive grid will be demonstrated using simulated data in the basic case of range only in pulse radar.We consider the LFM pulse from Equation 33 and set parameter b to be |b| = 1, the pulse length T p = 25, the bandwidth B equal to the normalized sampling frequency f s , B = f s = 1, PRT = 108, and variance γ = 1.Note that no compressive acquisition is involved yet but regular Nyquist sampling in the observation grid (as we focus here on the estimation grid).
We first determine the total information distance during the listening time on the PRT: d(T p , PRT).Then, we determine a desirable number of cells N cells : for this particular case, we set N cells = PRT − T p = 83, which is the same number of cells as for the reference case.From this, in Equations 30 and 32, we calculate the constant information distance δ and construct the grid cells, which are found in Figure 3.Note that our analysis holds for all the radar types providing G(τ ) from Equation 26with the waveform specific constant C in and the parameter b that includes the gain, RCS, etc.The sensing matrices A and corresponding coherences c(A) are illustrated in Figure 4. SSP on the grids is shown in Figure 5.
The reference grid (Figure 3) causes the coherence to be low and its maximum is constant for all ranges (Figure 4a).The information-based grid (Figure 3) leads to higher coherence for close ranges because of the smaller cells (Figure 4b).An optimal smallest grid size is to be designed as a compromise between the information-based distance (Equation 29), the incoherence of the sensing matrix (Equation 3), operational needs, and the available processing power.
The coherence of the sensing matrix c(A) is taken into consideration after the grid cells are determined with the constant information distance δ.For very small grid cells, the coherence is very high.An attempt in lowering the coherence is made by setting the minimum cell size to half the cell size of the reference case.For the test case, this means that for all cells with τ < 37 the cell size is fixed, as is seen in Figure 3 and Figure 4c.
In a range-only case, we show an example of the radar operation where closer ranges may have a higher priority.In our test case, there are three pairs of targets at different ranges from radar, where each pair consists of two targets that are equally distant from each other at every of  Note that in (c), the sensing matrix contains less elements because a lower bound on cell size has been applied in order to limit the coherence, hence the lower number of column combinations (i, j).
the three different ranges (Figure 5).The distance between the two targets in a pair is equal to one unit cell from the reference grid.On the reference grid (Figure 5a), the targets cannot be resolved at any range.SSP on both irregular grids (Figure 5b,c) enables resolving the two targets at closer ranges but on the coherence-adjusted grid with less cells.The coherence-adjusted grid (Figure 5c) is an example how the grid can be adapted by limiting its cell size to the half as it is small enough for a given scenario, and moreover, it lowers coherence (Equation 3) as well as the computational demands.We applied an ad hoc cell size that illustrates the design compromise well.

Conclusions
An adaptive estimation grid is created as a result of a practical combination of IG and CS in a basic radar case of range estimation.IG is applied to design a resolution cell in a radar system as related to the Fisher information distances on the statistical manifold of a data model with raw measurements and the range-dependent amplitude.In addition, we also explore how such an irregular estimation grid suits radar processing in practice, propose applying SSP from CS, and adjust the grid accordingly for this type of processing.The IG-based grid and its processing by SSP enable the adaptive grid adjustments that ensure better radar performance and more efficient radar processing because the processing power and resolution are optimized over ranges.Thus, such an adaptive grid compromises the IG-based grid with the CS-based sensing incoherence and naturally, also with the available processing power and operational demands.In existing radar processing, there is no flexibility in the FFT-based grid.
Future work on the combination of IG and CS will include extension of the optimal grid design to the range-Doppler-angle scenarios, as well as waveform design and other effects of the observation grid caused by compressive acquisition (typical for CS).

Figure 2
Figure 2 Fisher information metrics and grid cell sze.(a) The Fisher information metrics of the amplitude-dependent (constant information distance) model and the reference case.(b) Grid cell size as a function of distance for both models.

Figure 3
Figure 3 Cell sizes for the reference, information-based, and coherence-adjusted grids.

Figure 4
Figure 4Sensing matrices and corresponding coherences.Sensing matrix (left) and coherence plot (right) for the reference grid (a), the constant information distance grid (b), and the coherence-adjusted grid (c).Note that in (c), the sensing matrix contains less elements because a lower bound on cell size has been applied in order to limit the coherence, hence the lower number of column combinations (i, j).

Figure 5
Figure 5 Results of SSP on the reference grid (a), information-based grid (b), and coherence-adjusted grid (c).The closer targets are only resolved on the irregular grids.The coherence-adjusted grid gives the best results.