Scalar Representation of 2D Steady Vector Fields

Holger Theisel,Michael Motejat,Janos Zimmermann,Christian Rössl

2024-01-12

Abstract:We introduce a representation of a 2D steady vector field ${\mathbf v}$ by two scalar fields $a$, $b$, such that the isolines of $a$ correspond to stream lines of ${\mathbf v}$, and $b$ increases with constant speed under integration of ${\mathbf v}$. This way, we get a direct encoding of stream lines, i.e., a numerical integration of ${\mathbf v}$ can be replaced by a local isoline extraction of $a$. To guarantee a solution in every case, gradient-preserving cuts are introduced such that the scalar fields are allowed to be discontinuous in the values but continuous in the gradient. Along with a piecewise linear discretization and a proper placement of the cuts, the fields $a$ and $b$ can be computed. We show several evaluations on non-trivial vector fields.

Graphics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: how to represent a two - dimensional steady vector field \(v\) by two scalar fields \(a\) and \(b\), so that these scalar fields can directly encode streamlines and can solve integration and connectivity problems in local look - up, thus avoiding the error accumulation brought by numerical integration.

Specifically, the author hopes to find two smooth and differentiable scalar fields \(a(x)\) and \(b(x)\) that satisfy the following conditions:
1. \(\nabla a\neq0\) (Ensure that the gradient of \(a\) is not zero)
2. \(v^{T}\nabla a = 0\) (Ensure that the contour lines of \(a\) correspond to the streamlines of \(v\))
3. \(v^{T}\nabla b=1\) (Ensure that \(b\) increases at a unit speed along the streamlines)

These conditions ensure that the streamline integral can be simplified to the problem of the intersection of the contour lines of \(a\) and \(b\), that is:
\[a(\phi(x,\tau)) = a(x),\quad b(\phi(x,\tau)) = b(x)+\tau\]

However, in general, the scalar field \(a\) that satisfies the above conditions does not exist. To overcome this problem, the author introduced gradient - preserving cuts and special treatment of critical points, thus allowing to find scalar fields that meet the requirements under certain specific conditions.

In addition, the author also proposed a strictly local image - based flow visualization technique without any numerical integration. This method not only improves the computational efficiency but also reduces the error accumulation caused by numerical integration.

In summary, the main goal of this paper is to improve the visualization of two - dimensional steady vector fields by introducing a new scalar field representation method, with special attention to solving integration and connectivity problems.