Optical flow is one of the classical problems in computer vision, but it has recently also been adapted to applications from other fields, such as fluid mechanics and dynamical systems. If the goal is to analyze the dynamics of system whose evolution is governed by a flow field that is the gradient of a potential function -which describes many flows in fluid dynamics -it is natural to approach the optical flow problem by reconstructing the potential function, also called the stream function, rather than reconstructing the components of the flow directly. This alternate approach allows one to impose scientific priors, via regularization, directly on the flow itself rather than on its components independently. We demonstrate the stream function formulation of optical flow and its application to reconstructing an oceanic fluid flow driven by satellite measurements. It is also shown how these flow fields can be used to analyze mixing and mass transport in the fluid system being imaged.
Stream Function Optical FlowOptical flow is an image-based method for computing the apparent flow field governing a system evolution. Recently much work has been done using optical flow to analyze fluids, which requires a different physical model from the classical HornSchunk method [12]. Corpetti et al. [6,7] proposed an alternative based on the continuity equation, resulting in an optical flow energy of the formwhere Ω is the spatial domain, R(u, v) is an appropriate regularization scheme, and α > 0 is the regularization parameter. This formulation has become common for analyzing fluid flows, for example in connection with and in comparison to particle image velocimetry (PIV) [5,13], geophysical fluid flows [1,2], and atmospheric motion [4]. While it is possible to incorporate the full Navier-Stokes equations for fluid flows in order to more accurately model the evolution of a fluid system [8], the model (1) is robust and results in more straightforward computational schemes. Under the assumption that the flow being imaged is a "potential flow," i.e. that the flow field is the gradient of a potential, or stream function, it is natural to instead formulate the optical flow problem in terms of reconstructing the stream function directly. Within a Hamiltonian framework, the evolution equation governing the flow is thus I t = −div (I∇ H ψ), where ∇ H ψ = −ψ y , ψ x is the symplectic gradient. The resulting optical flow energy iswhere R(ψ) is a regularization of the stream function. There are two primary goals of regularizing a variational optimization problem. The first goal is to ensure that the optimization problem is well-posed, i.e. that a solution exists, is unique, and is stable with respect to perturbations in the input data. For optical flow problems, this need not happen in general. In particular, for nearly any physically appropriate regularization, there exist data for which the governing flow field is not unique. The second reason for regularization -and one of the main advantages to computing optical flow via (2) -is to im...