A collaborative method for dense motion estimation of fluid flows
(back)

 

Description

Optical-flow methods dedicated to fluid motion have been recently introduced. These approaches have the main advantage to produce dense motion fields. However, as the underlying optimization are based on gradient descent techniques, the estimates are very sensitive to the initialization procedure. In particular, large displacement can only be recovered using a multiresolution strategy with the possible drawback of smoothing out small structure motions.
Fluid motion estimation techniques based on cross-correlation produce sparse motion fields. The estimation is in addition local and thus prone to erroneous spatial variability. However they are generally locally very robust to noise and may recover large displacement of some traceable features, that is to say features sufficiently contrasted and persisting over time on consecutive images.


To overcome the limitation induced by both optical-flow and correlation based methods, we propose to introduce an initialization scheme for a fluid dedicated dense motion estimation where the multiresolution scheme has been removed. The initial velocity field is estimated by minimization of a energy function composed of a 'div-curl' regularization and a data term. The latter constrains the dense estimated field to fit the sparse correlation vectors convoluted by bidimensional gaussian functions. The resulting vector field is then used to initialize the finest resolution of the fluid dedicated dense motion estimator introduced by [Corpetti et al., 2002].



Results

Comparison of multiresolution and collaborative schemes in the case of a image sequence depicting two-dimensional motion of a passive scalar. Image sequence has been provided by Cemagref, Rennes.
Figure: Above : vorticity provided by the DNS (left), vorticity estimation by the fluid flow dedicated multiresolution approach of [Corpetti et al., 2002] (right). Below : vorticity estimation after the first (left) and the second (right) level of the collaborative scheme.

In order to analyze visually the accuracy of the estimation method, the trajectories of uniformly spread points have been reconstructed using a fourth-order Runge-Kutta integration method. The following image sequence has been provided by the Laboratoire de Météorologie Dynamique (LMD) du CNRS.


Figure: Comparison between trajectories estimated without (a) and with the initialization procedure (b). These trajectories correspond to the highest layer
Image LayerClass
a) b)

References

  1. T. Corpetti, E. Mémin, and P. Pérez. Dense estimation of fluid flows. IEEE Trans. Pattern Anal. Machine Intell., 24(3):365– 380, 2002.
  2. P.Heas, E. Mémin, and N. Papadakis. Dense estimation of layer motions in the atmosphere, International Conference on Pattern Recognition, Hong-Kong, 2006.
  3. A. Szantai, P.Heas, E. Mémin, Comparison of atmospheric motion vectors and dense vector fields calculated from MSG images, 8th international winds workshop Beijing, China, 2006