On convergence of the Horn and Schunck optical-flow estimation method.

The purpose of this study is to prove convergence results for the Horn and Schunck optical-flow estimation method. Horn and Schunck stated optical-flow estimation as the minimization of a functional. When discretized, the corresponding Euler-Lagrange equations form a linear system of equations We write explicitly this system and order the equations in such a way that its matrix is symmetric positive definite. This property implies the convergence Gauss-Seidel iterative resolution method, but does not afford a conclusion on the convergence of the Jacobi method. However, we prove directly that this method also converges. We also show that the matrix of the linear system is block tridiagonal. The blockwise iterations corresponding to this block tridiagonal structure converge for both the Jacobi and the Gauss-Seidel methods, and the Gauss-Seidel method is faster than the (sequential) Jacobi method.

Constraining active contour evolution via lie groups of transformation.

We present a novel approach to constraining the evolution of active contours used in image analysis. The proposed approach constrains the final curve obtained at convergence of curve evolution to be related to the initial curve from which evolution begins through an element of a desired Lie group of plane transformations. Constraining curve evolution in such a way is important in numerous tracking applications where the contour being tracked in a certain frame is known to be related to the contour in the previous frame through a geometric transformation such as translation, rotation, or affine transformation, for example. It is also of importance in segmentation applications where the region to be segmented is known up to a geometric transformation. Our approach is based on suitably modifying the Euler-Lagrange descent equations by using the correspondence between Lie groups of plane actions and their Lie algebras of infinitesimal generators, and thereby ensures that curve evolution takes place on an orbit of the chosen transformation group while remaining a descent equation of the original functional. The main advantage of our approach is that it does not necessitate any knowledge of nor any modification to the original curve functional and is extremely straightforward to implement. Our approach therefore stands in sharp contrast to other approaches where the curve functional is modified by the addition of geometric penalty terms. We illustrate our algorithm on numerous real and synthetic examples.

Multiple motion segmentation with level sets.

Segmentation of motion in an image sequence is one of the most challenging problems in image processing, while at the same time one that finds numerous applications. To date, a wealth of approaches to motion segmentation have been proposed. Many of them suffer from the local nature of the models used. Global models, such as those based on Markov random fields, perform, in general, better. In this paper, we propose a new approach to motion segmentation that is based on a global model. The novelty of the approach is twofold. First, inspired by recent work of other researchers we formulate the problem as that of region competition, but we solve it using the level set methodology. The key features of a level set representation, as compared to active contours, often used in this context, are its ability to handle variations in the topology of the segmentation and its numerical stability. The second novelty of the paper is the formulation in which, unlike in many other motion segmentation algorithms, we do not use intensity boundaries as an accessory; the segmentation is purely based on motion. This permits accurate estimation of motion boundaries of an object even when its intensity boundaries are hardly visible. Since occasionally intensity boundaries may prove beneficial, we extend the formulation to account for the coincidence of motion and intensity boundaries. In addition, we generalize the approach to multiple motions. We discuss possible discretizations of the evolution (PDE) equations and we give details of an initialization scheme so that the results could be duplicated. We show numerous experimental results for various formulations on natural images with either synthetic or natural motion.