Abstract-Robust analysis of vector fields has been established as an important tool for deriving insights from the complex systems these fields model. Traditional analysis and visualization techniques rely primarily on computing streamlines through numerical integration. The inherent numerical errors of such approaches are usually ignored, leading to inconsistencies that cause unreliable visualizations and can ultimately prevent in-depth analysis. We propose a new representation for vector fields on surfaces that replaces numerical integration through triangles with maps from the triangle boundaries to themselves. This representation, called edge maps, permits a concise description of flow behaviors and is equivalent to computing all possible streamlines at a user defined error threshold. Independent of this error streamlines computed using edge maps are guaranteed to be consistent up to floating point precision, enabling the stable extraction of features such as the topological skeleton. Furthermore, our representation explicitly stores spatial and temporal errors which we use to produce more informative visualizations. This work describes the construction of edge maps, the error quantification, and a refinement procedure to adhere to a user defined error bound. Finally, we introduce new visualizations using the additional information provided by edge maps to indicate the uncertainty involved in computing streamlines and topological structures.Index Terms-Vector Fields, Error Quantification, Edge Maps.
!
MOTIVATIONSV ECTOR fields are a common form of simulation data appearing in a wide variety of applications ranging from computational fluid dynamics (CFD) and weather prediction to engineering design. Visualizing and analyzing the flow behavior of these fields can help provide critical insights into simulated physical processes. However, achieving a consistent and rigorous interpretation of vector fields is difficult, in part because traditional numerical techniques for integration do not preserve the expected invariants of vector fields.To better understand this challenge inherent in traditional numerical techniques, we reconsider the most common way to store vector fields. Both a discretization of the domain of the field (often in the form of a triangulated mesh) as well as a set of sample vectors (defined at the vertices of the mesh) are required. The vector field on the interior of a triangle is approximated by interpolating vector values from the samples at the triangle's corners. Subsequently, computing properties that require integrating these vector values presents a significant computational challenge. For example, consider computing the flow paths (streamlines) of massless particles that travel using the instantaneous velocity defined by the field. Naive integration techniques may violate the property that every two of these paths are expected to be pairwise disjoint