On two-dimensional Markov spectral estimation

We give systematic methods for two-dimensional (2-D) spectral estimation from raw data using Gaussian Markov random field (MRF) models. The MRF models considered in this paper characterize the statistical dependence of the observation y(s) on its neighbors in all directions. Due to the modeling assumption, the 2-D spectrum is an explicit function of the parameters of the MRF model. Thus, the spectral estimation problem reduces to that of estimating the appropriate structure and the parameters of the MRF model. A tractable algorithm for maximum likelihood (ML) estimation of parameters is obtained by using doubly periodic boundary conditions for the given data. Asymptotically consistent decision rules are used to choose the appropriate neighbor set of the MRF model. The MRF spectral estimate obtained using the doubly periodic boundary conditions has the property that its inverse discrete Fourier transforms are in perfect agreement with the periodic sample correlation values of the given observations, over a symmetric window identical to the structure of the MRF model. Using this result, and the fact that the 2-D maximum entropy spectral (MES) estimate has a structure similar to the MRF power spectrum, we show that the MRF spectral estimate developed in this paper is a good approximation to the MES estimate. This result provides a reasonable basis (i.e., maximum likelihood) on which to justify the use of the MES estimate with correlation estimates obtained from a finite piece of data samples.

[1]  J. Capon.  High-resolution frequency-wavenumber spectrum analysis , 1969 .

[2]  John W. Woods,et al.  Two-dimensional discrete Markovian fields , 1972, IEEE Trans. Inf. Theory.

[3]  P. Moran.  A Gaussian Markovian process on a square lattice , 1973, Journal of Applied Probability.

[4]  H. Akaike.  A new look at the statistical model identification , 1974 .

[5]  J. Besag,et al.  On the estimation and testing of spatial interaction in Gaussian lattice processes , 1975 .

[6]  John W. Woods,et al.  Two-dimensional Markov spectral estimation , 1976, IEEE Trans. Inf. Theory.

[7]  Azriel Rosenfeld,et al.  A Comparative Study of Texture Measures for Terrain Classification , 1975, IEEE Transactions on Systems, Man, and Cybernetics.

[8]  R. Kashyap.  Inconsistency of the AIC rule for estimating the order of autoregressive models , 1980 .

[9]  M. Ekstrom,et al.  Multidimensional spectral factorization and unilateral autoregressive models , 1980 .

[10]  Ramalingam Chellappa,et al.  Decision rules for choice of neighbors in random field models of images , 1981 .

[11]  Dag Tjostheim,et al.  Autoregressive Modeling And Spectral Analysis of Array Data in the Plane , 1981, IEEE Transactions on Geoscience and Remote Sensing.

[12]  A.K. Jain,et al.  Advances in mathematical models for image processing , 1981, Proceedings of the IEEE.

[13]  T. Ulrych,et al.  HIGH RESOLUTION 2-DIMENSIONAL POWER SPECTRAL ESTIMATION , 1981 .

[14]  Jae Lim,et al.  A new algorithm for two-dimensional maximum entropy power spectrum estimation , 1981 .

[15]  J. Cadzow,et al.  Two-dimensional spectral estimation , 1981 .

[16]  R. L. Kashyap,et al.  Analysis and Synthesis of Image Patterns by Spatial Interaction Models , 1981 .

[17]  J.H. McClellan,et al.  Multidimensional spectral estimation , 1982, Proceedings of the IEEE.

[18]  Rama Chellappa,et al.  Estimation and choice of neighbors in spatial-interaction models of images , 1983, IEEE Trans. Inf. Theory.