Maximum-likelihood parameter estimation of the harmonic, evanescent, and purely indeterministic components of discrete homogeneous random fields

This paper presents a maximum-likelihood solution to the general problem of fitting a parametric model to observations from a single realization of a two-dimensional (2-D) homogeneous random field with mixed spectral distribution. On the basis of a 2-D Wold-like decomposition, the field is represented as a sum of mutually orthogonal components of three types: purely indeterministic, harmonic, and evanescent. The suggested algorithm involves a two-stage procedure. In the first stage, we obtain a suboptimal initial estimate for the parameters of the spectral support of the evanescent and harmonic components. In the second stage, we refine these initial estimates by iterative maximization of the conditional likelihood of the observed data, which is expressed as a function of only the parameters of the spectral supports of the evanescent and harmonic components. The solution for the unknown spectral supports of the harmonic and evanescent components reduces the problem of solving for the other unknown parameters of the field to a linear least squares. The Cramer-Rao lower bound on the accuracy of jointly estimating the parameters of the different components is derived, and it is shown that the bounds on the purely indeterministic and deterministic components are decoupled. Numerical evaluation of the bounds provides some insight into the effects of various parameters on the achievable estimation accuracy. The performance of the maximum-likelihood algorithm is illustrated by Monte Carlo simulations and is compared with the Cramer-Rao bound