1 / 29

Representing Higher Order Vector Fields Singularities on Piecewise Linear Surfaces

Representing Higher Order Vector Fields Singularities on Piecewise Linear Surfaces. Wan Chiu Li Bruno Vallet Nicolas Ray Bruno Lévy IEEE Visualization 2006 Baltimore, USA alice.loria.fr. Outline. Introduction Discrete representation Singularities Encoding an existing vector field

delphina
Download Presentation

Representing Higher Order Vector Fields Singularities on Piecewise Linear Surfaces

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Representing Higher Order Vector Fields Singularities on Piecewise Linear Surfaces Wan Chiu Li Bruno Vallet Nicolas Ray Bruno Lévy IEEE Visualization 2006 Baltimore, USA alice.loria.fr

  2. Outline • Introduction • Discrete representation • Singularities • Encoding an existing vector field • LIC-based visualization • Conclusion

  3. Introduction • What is a vector field singularity ? • It is a 0 of the field • How can we characterize singularities ? • By their index = ∫ dθ -1 +1

  4. Introduction • What is a vector field singularity ? • It is a 0 of the field • How can we characterize singularities ? • By their index = ∫ dθ -2 +2

  5. Introduction • How can we visualize a singularity ? • Piecewise linear methods: index∊[-val/2+1, val/2+1] [Tricoche00]

  6. Introduction • How can we visualize a singularity ? • Piecewise linear methods: index∊[-val/2+1, val/2+1] • Higher order singularities: index∊ℤ

  7. Introduction • Basic idea: • 2D vectors are complex reiθ • Interpolate r and θ • Justification: • Singularity = (r =0, θ undefined) • Singularity index depend only on θ in neighborhood

  8. Discrete representation • On triangulated meshes • Dual vertex-edge encoding using polar coordinates: • Dual vertices: norm r∊ℝ, angle θ∊ℝ • Dual edges: period jump p∊ ℤ • 3 Step interpolation (0D, 1D, 2D) • Dual vertices=Facet centers (0D) • Dual edges (1D) • Subdivision simplex (2D)

  9. 0D r(v*) θ(v*) v* x(v*) • θ(v*) : measured from a reference vectorx(v*) • r(v*) : vector norm, basis independent

  10. 1D v* θ(v*) P v’* x(v*) e* θ(v’*) x(v’*) Linear interpolation: θ(P) = θ(v*) + ∆θ(e*)⍺(e*)t ⍺(e*): height ratio, ∆θ(e*): angular variation alonge*

  11. 1D v* θ(v*) v’* P x(v*) e* θ(v’*) x(v’*) Linear interpolation: θ(P) = θ(v*) + ∆θ(e*)(1-⍺(e*))t ⍺(e*): height ratio, ∆θ(e*): angular variation alonge*

  12. 1D v* Height ratio: ⍺(e*) =h/H 1-⍺(e*) =h’/H e* h H h’ v’*

  13. 1D Angular variation alonge* : Period Jump: ∆θ(e*)= ∫dθ = θ(B) -θ(A) + 2 πp(e*) e* B p(e*) = -1 e* A B p(e*) = 0 e* A B p(e*) = 1 e* A

  14. 1D

  15. 2D Subdivision simplex

  16. 2D Subdivision simplex

  17. 2D • A variant the side-vertex interpolation [Nielson79] • Linear along the side • Constant along a side-vertex path (=side value) side P’ P vertex

  18. Singularities • Singularities may occur only at vertices • Singularity indexdepends only on period jumps : • I(v)= ∫dθ = I0(v)+ ∑ p(e*) e*∊∂f* ∂f* f* v

  19. Singularities • Advantages: • Control over placement and index of singularities • Coherent with Poincare-Hopf index theorem • Index independent of the valence • Easy extension to fractional indices

  20. Extension to fractional indices • Fractional indices appear in N-symmetry vector fields : • Not vectors but equivalence class of vectors by ≡N • u≡Nv ⇔ ∃k | u=R(v, 2kπ/N) • Period jump and Indices are multiples of 1/N -1/2 +1/2

  21. Extension to fractional indices • Fractional indices appear in N-symmetry vector fields : • Not vectors but equivalence class of vectors by ≡N • u≡Nv ⇔ ∃k | u=R(v, 2kπ/N) • Period jump and Indices are multiples of 1/N -1/4 +1/4

  22. Encoding an Existing Vector Field • r(v*) : norm of the vector at facet center v* • θ(v*): choose one of the 3 edges and compute angle • 2πp(e*)=∆θ(e*) - ∠(x(v’*), x(v*)) - θ(v’*) + θ(v*) • Requires an interpolation or an analytic form: vxdvy – vydvx ∆θ= ∫dθ = ∫ ||v||2 e* e*

  23. Encoding an Existing Vector Field ∆θ= ∫dθ = ∫ (vxdvy – vydvx) / ||v||2 e* e* dθ v* v’* e*

  24. LIC-based Visualization • GPU accelerated • Works in image space • 3 passes [Laramee et. al. 03] [Van Wijk 03] Ensure geometric discontinuity (in depth buffer) Line integral convolution (in image space) Direction on the surface (in fragment shader)

  25. Results Index = -3 Index = 5

  26. Results

  27. Conclusion • Dual vertex-edge encoding • + 3 steps interpolation: • A very good candidate for visualizing non-linear vector fields on piecewise linear surfaces or 2d meshes. • Efficient and simple way to visualize arbitrary singularities • Easy generalization to fractional indices • Easier particularization to 2d fields

  28. Future work • Smooth non singular vertices • Topological operations • Trace streamlines • Extension to 3d vector fields

  29. alice.loria.fr Questions ?

More Related