Abstract
We present in this paper a recursive-in-order least-squares (LS) algorithm to compute efficiently the parameters of a 2-D Gaussian Markov random field (GMRF) model. The algorithm is based on the fact that the least-squares estimation of the parameters of a 2-D noncausal GMRF model is consistent and the coefficient matrix in the normal equation has near-to-block-Toeplitz structure. Hence, it has a Levinson-like form for the updating of model parameters by introducing auxiliary variables. Moreover, this paper proposes the concept ofrecursive path for 2-D recursive-in-order algorithms, and points out that there exists a tradeoff between fast computation of the parameters and accurate choice of model support; a compromise recursive path is then suggested where the orders change alternately in two directions. The computational complexity of the developed algorithm is analyzed, and the results show that the algorithm is more efficient when either the image size or the model support is larger. It is found that the total number of multiplications (mps) involved in the new algorithm is only about 14% of that in the conventional LS method when the image size is 512 × 512 and the neighbor set of the model is a 17 × 17 window. Computer simulation results using the recursive-in-order algorithm developed in this paper and the conventional LS method are given to verify the correctness of the new algorithm.
Similar content being viewed by others
References
C. Therien, T. Quatieri, and D. Dudgeon, Statistical model-based algorithms for image analysis,Proc. IEEE, Vol. 74, pp. 532–551, April, 1986.
A. Rosenfeld (editor),Image Models, Academic Press, New York, 1981.
R. Chellappa and A. A. Sawchuk (editors),Digital Image Processing and Analysis: Vols. 1 and 2, IEEE Computer Society Press, 1985.
R. Chellappa, Two-dimensional discrete Gaussian Markov random field models for image processing,Progress in Pattern Recognition 2, L. N. Kanal and A. Rosenfeld, eds., North Holland Publ. Co., New York, NY, 1985, pp. 79–112.
R. Chellappa and R. L. Kashyap, Digital image restoration using spatial interaction models,IEEE Trans. Acoustics Speech Signal Processing, Vol. ASSP-30, No. 3, pp. 461–472, 1982.
S. Geman and G. Geman, Stochastic relaxation: Gibbs distribution, and the Bayesian restoration of images,IEEE Trans. Pattern Analysis and Machine Intelligence, Vol. PAMI-6, No. 6, pp. 721–741, 1984.
G. R. Cross and A. K. Jain, Markov random field texture models,IEEE Trans. Pattern Analysis and Machine Intelligence, Vol. PAMI-5, pp. 25–39, 1983.
M. Hassner and H. Elliott, Image segmentation using simple Markov random field models,Computer Graphics and Image Processing, Vol. 20, pp. 357–370, April, 1980.
H. Derin, H. Elliot, R. Cristi, and D. Geman, Bayes smoothing algorithms for segmentation of binary images modeled by Markov random fields,IEEE Trans. Pattern Analysis and Machine Intelligence, Vol. PAMI-6, No. 6, pp. 707–720, 1984.
M. Hassner and J. Sklansky, The use of Markov random field as model of textures,Computer Graphics and Image Processing, Vol. 12, pp. 357–370, April, 1980.
P. Bouthemy and P. Lalaande, Motion detection in a sequence using Gibbs distributions,Proc. of ICASSP, pp. 1651–1654, 1989.
B. Chalmond, Image restoration using an estimated Markov model,Signal Processing, Vol. 15, pp. 115–129, 1988.
R. L. Kashyap and R. Chellappa, Estimation and choice of neighbors in spatial-interaction models of images,IEEE Trans. Information Theory, Vol. IT-29, No. 1, pp. 60–72, January, 1983.
P. Y. Zhao and Z. Y. He, Efficient order recursive LS algorithms for 2D SAR model fitting,Proc. of ICASSP, pp. 1200–1203,1988.
N. Kalouptsidis, G. Carayannis, D. Manolakis, and E. Koukoutsis, Efficient recursive in order least squares FIR filtering and prediction,IEEE Trans. Acoustics Speech Signal Processing, Vol. ASSP-33, No. 4, pp. 1175–1186, October, 1985.
G. A. Glentis and N. Kalouptsidis, Efficient order recursive algorithms for multichannel least squares filtering,IEEE Trans. Signal Processing, Vol. 40, No. 6, pp. 1354–1374, June, 1992.
C. R. Zou,Two-dimensional random field models — recursive least squares parameter estimation algorithms and applications to image processing, Doctoral dissertation of Southeast University, Nanjing, China, 1990.
Author information
Authors and Affiliations
Additional information
This work was supported by the NSERC under Grants A-4070 and A-7739, and by the FCAR, Grant H-70.
Rights and permissions
About this article
Cite this article
Zou, C.R., Plotkin, E.I., Swamy, M.N.S. et al. Recursive-in-order least-squares parameter estimation algorithm for 2-D noncausal Gaussian Markov random field model. Circuits Systems and Signal Process 14, 87–110 (1995). https://doi.org/10.1007/BF01183750
Received:
Accepted:
Issue Date:
DOI: https://doi.org/10.1007/BF01183750