A Proof of Convergence of the Horn-Schunck Optical Flow Algorithm in Arbitrary Dimension

The Horn and Schunck (HS) method, which amounts to the Jacobi iterative scheme in the interior of the image, was one of the first optical flow algorithms. In this article, we prove the convergence of the HS method, whenever the problem is well-posed. Our result is shown in the framework of a generalization of the HS method in dimension n ≥ 1, with a broad definition of the discrete Laplacian. In this context, the condition for the convergence is that the intensity gradients are not all contained in a same hyperplane. Two other articles ([17] and [13]) claimed to solve this problem in the case n = 2, but it appears that both of these proofs are erroneous. Moreover, we explain why some standard results on the convergence of the Jacobi method do not apply for the HS problem, unless n = 1. It is also shown that the convergence of the HS scheme implies the convergence of the Gauss-Seidel and SOR schemes for the HS problem.

[1]  L. Schwartz.  Cours d'analyse , 1963 .

[2]  R. Varga,et al.  Block diagonally dominant matrices and generalizations of the Gerschgorin circle theorem , 1962 .

[3]  Du-Ming Tsai,et al.  Low-contrast surface inspection of mura defects in liquid crystal displays using optical flow-based motion analysis , 2010, Machine Vision and Applications.

[4]  Richard P. Wildes,et al.  Recovering Estimates of Fluid Flow from Image Sequence Data , 2000, Comput. Vis. Image Underst..

[5]  Tianshu Liu,et al.  Fluid flow and optical flow , 2008, Journal of Fluid Mechanics.

[6]  Christoph Schnörr,et al.  Determining optical flow for irregular domains by minimizing quadratic functionals of a certain class , 1991, International Journal of Computer Vision.

[7]  Amar Mitiche,et al.  On convergence of the Horn and Schunck optical-flow estimation method , 2004, IEEE Transactions on Image Processing.

[8]  Anne E. Trefethen,et al.  On the Usage of GPUs for Efficient Motion Estimation in Medical Image Sequences , 2011, Int. J. Biomed. Imaging.

[9]  David J. Evans,et al.  Parallel S.O.R. iterative methods , 1984, Parallel Comput..

[10]  Joachim Weickert,et al.  Variational Optic Flow Computation with a Spatio-Temporal Smoothness Constraint , 2001, Journal of Mathematical Imaging and Vision.

[11]  Berthold K. P. Horn,et al.  Determining Optical Flow , 1981, Other Conferences.

[12]  On the refinements of the Jensen-Steffensen inequality , 2011 .

[13]  Roberto Bagnara,et al.  A Unified Proof for the Convergence of Jacobi and Gauss-Seidel Methods , 1995, SIAM Rev..

[14]  Geoffrey G. Zhang,et al.  Intrathoracic tumour motion estimation from CT imaging using the 3D optical flow method. , 2004, Physics in medicine and biology.

[15]  Elsa D. Angelini,et al.  Coronary Occlusion Detection with 4D Optical Flow Based Strain Estimation on 4D Ultrasound , 2009, FIMH.

[16]  Anand Kumar.  A discretization of the n–dimensional Laplacian for a dimension–independent stability limit , 2001, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[17]  Antonin Chambolle,et al.  A First-Order Primal-Dual Algorithm for Convex Problems with Applications to Imaging , 2011, Journal of Mathematical Imaging and Vision.

[18]  Louis A. Hageman,et al.  Iterative Solution of Large Linear Systems. , 1971 .

[19]  J. Jensen.  Sur les fonctions convexes et les inégalités entre les valeurs moyennes , 1906 .

[20]  K. R. James.  Convergence of Matrix Iterations Subject to Diagonal Dominance , 1973 .

[21]  Atsushi Imiya,et al.  A Convergence Proof for the Horn-Schunck Optical-Flow Computation Scheme Using Neighborhood Decomposition , 2008, IWCIA.

[22]  Gerhard Wellein,et al.  Efficient multicore-aware parallelization strategies for iterative stencil computations , 2010, J. Comput. Sci..

[23]  Otmar Scherzer,et al.  Optical Flow on Evolving Surfaces with an Application to the Analysis of 4D Microscopy Data , 2013, SSVM.

[24]  Richard Szeliski,et al.  A Database and Evaluation Methodology for Optical Flow , 2007, 2007 IEEE 11th International Conference on Computer Vision.