Querying approximate shortest paths in anisotropic regions Siu-Wing - - PowerPoint PPT Presentation

querying approximate shortest paths in anisotropic regions
SMART_READER_LITE
LIVE PREVIEW

Querying approximate shortest paths in anisotropic regions Siu-Wing - - PowerPoint PPT Presentation

Querying approximate shortest paths in anisotropic regions Siu-Wing Cheng (HKUST) Hyeon-Suk Na (Soongsil University) Antoine Vigneron (INRA Jouy-en-Josas) Yajun Wang (HKUST) The weighted region problem cost = weight length t w=1 w=2 s


slide-1
SLIDE 1

Querying approximate shortest paths in anisotropic regions

Siu-Wing Cheng (HKUST) Hyeon-Suk Na (Soongsil University) Antoine Vigneron (INRA Jouy-en-Josas) Yajun Wang (HKUST)

slide-2
SLIDE 2

t s cost = weight × length w=1 w=2 Shortest path

The weighted region problem

slide-3
SLIDE 3

t s cost = weight × length w=1 w=2 Shortest path

The weighted region problem

slide-4
SLIDE 4
  • Triangulation with n vertices
  • Each face has a weight
  • Find an approximate shortest path with cost at most 1 + ǫ

times the minimum s t

The weighted region problem

slide-5
SLIDE 5

No (strongly polynomial) FPTAS is known

slide-6
SLIDE 6

No (strongly polynomial) FPTAS is known Known results: polynomial in n, 1/ǫ, and other parameters.

slide-7
SLIDE 7

No (strongly polynomial) FPTAS is known Known results: polynomial in n, 1/ǫ, and other parameters.

  • ρ: maximum weight/ minimum weight
  • θ: minimum angle in the triangulation
slide-8
SLIDE 8

No (strongly polynomial) FPTAS is known Known results: polynomial in n, 1/ǫ, and other parameters.

  • ρ: maximum weight/ minimum weight
  • θ: minimum angle in the triangulation

Known results

  • Mitchell and Papadimitriou (1987): O(n8 log(nNρ/ǫ))
  • Aleksandrov, Maheshwari, and Sacks:

O(Cn/√ǫ × polylog(n, 1/ǫ)), C depends on ρ, θ and other parameters.

  • Sun and Reif: O(C′n/ǫ × polylog(n, 1/ǫ)),

C′ depends on θ, other parameters, but not ρ.

  • Cheng et al. (2007): O(n3ρ/ǫ × polylog(n, ρ, 1/ǫ))
slide-9
SLIDE 9

Query problem

s

slide-10
SLIDE 10

Query problem

s Build a data structure

slide-11
SLIDE 11

Query problem

s Build a data structure q Given a query point q, find a (1 + ε)-approximation of the cost of the shortest path from s to q.

slide-12
SLIDE 12

Query problem

s Build a data structure q Given a query point q, find a (1 + ε)-approximation of the cost of the shortest path from s to q.

slide-13
SLIDE 13

Query problem

s Build a data structure q Given a query point q, find a (1 + ε)-approximation of the cost of the shortest path from s to q. Output an approximate shortest path in time linear in its size

slide-14
SLIDE 14

Previous results

Aleksandrov et al.: Space and preprocessing time depend on several parameters including n, ε, ρ and θ.

slide-15
SLIDE 15

Previous results

Aleksandrov et al.: Space and preprocessing time depend on several parameters including n, ε, ρ and θ. Applies to polyhedral surfaces

slide-16
SLIDE 16

Previous results

Aleksandrov et al.: Space and preprocessing time depend on several parameters including n, ε, ρ and θ. Applies to polyhedral surfaces

Our results

Query time: O(log(ρn/ε)) Space: O(ρ2n4/ε2 log(ρn/ε)) Preprocessing time: O(ρ2n4/ε2 log2(ρn/ε))

slide-17
SLIDE 17

Previous results

Aleksandrov et al.: Space and preprocessing time depend on several parameters including n, ε, ρ and θ. Applies to polyhedral surfaces

Our results

Query time: O(log(ρn/ε)) Space: O(ρ2n4/ε2 log(ρn/ε)) Preprocessing time: O(ρ2n4/ε2 log2(ρn/ε)) Applies to anisotropic distance functions

slide-18
SLIDE 18

Convex distance function

Unit ball B: set of points at distance 1 from O O B

slide-19
SLIDE 19

Convex distance function

Unit ball B: set of points at distance 1 from O O x y x + λB dB(x, y) = λ B

slide-20
SLIDE 20

Convex distance function

Unit ball B: set of points at distance 1 from O O x y x + λB dB(x, y) = λ B Shortest path is a straight line segment

slide-21
SLIDE 21

Model

1 1/ρ B One convex distance function per face

slide-22
SLIDE 22

Model

1 1/ρ B One convex distance function per face Speed is in interval [1/ρ, 1]

slide-23
SLIDE 23

Model

1 1/ρ B One convex distance function per face Speed is in interval [1/ρ, 1] Cost of a length ℓ path is in [ℓ, ρℓ].

slide-24
SLIDE 24

s t

Exemple of a non-polygonal shortest path

slide-25
SLIDE 25

s t

Exemple of a non-polygonal shortest path

slide-26
SLIDE 26

s t

Exemple of a non-polygonal shortest path

There exists a rectifiable (=finite length) shortest path

slide-27
SLIDE 27

s t

Exemple of a non-polygonal shortest path

There exists a rectifiable (=finite length) shortest path There exists a (1 + ε)-approximate shortest path with O(ρn2/ε) edges.

slide-28
SLIDE 28

s t

Static (non query) algorithm

slide-29
SLIDE 29

s t Circle with radius ρ d(s, t) and center s

Static (non query) algorithm

slide-30
SLIDE 30

s t Circle with radius ρ d(s, t) and center s Steiner points with uniform spacing δ = 1/poly(ρ, n, 1/ε)

Static (non query) algorithm

slide-31
SLIDE 31

s t Circle with radius ρ d(s, t) and center s Steiner points with uniform spacing δ = 1/poly(ρ, n, 1/ε)

Static (non query) algorithm

slide-32
SLIDE 32

Observation

s q δ

slide-33
SLIDE 33

Observation

s q shortest path δ

slide-34
SLIDE 34

Observation

s q shortest path approximate shortest path δ

slide-35
SLIDE 35

Observation

s q shortest path approximate shortest path δ

slide-36
SLIDE 36

Observation

s q shortest path approximate shortest path δ additive error at most δ

slide-37
SLIDE 37

Observation

(1 + ε) approximation when q is in an annulus s q t radius ratio=poly(ρ, n, 1/ε)

slide-38
SLIDE 38

Answering queries within the annulus

s Additively weighted Voronoi diagram

slide-39
SLIDE 39

Answering queries within the annulus

s Weight= cost of approximate shortest path Additively weighted Voronoi diagram

slide-40
SLIDE 40

Answering queries within the annulus

s Weight= cost of approximate shortest path Additively weighted Voronoi diagram q

slide-41
SLIDE 41

Answering queries within the annulus

s Weight= cost of approximate shortest path Additively weighted Voronoi diagram q

slide-42
SLIDE 42

Answering queries within the annulus

s Weight= cost of approximate shortest path Additively weighted Voronoi diagram q

slide-43
SLIDE 43

Answering queries within the annulus

s Weight= cost of approximate shortest path Additively weighted Voronoi diagram q

slide-44
SLIDE 44

Overview

For each vertex vi, we construct a data structure for the annulus with radius ri = d(s, vi)/2ρ and Ri =poly(ρ, n, 1/ε)d(s, vi)

slide-45
SLIDE 45

Overview

For each vertex vi, we construct a data structure for the annulus with radius ri = d(s, vi)/2ρ and Ri =poly(ρ, n, 1/ε)d(s, vi) s vi ri Ri q

slide-46
SLIDE 46

Overview

For each vertex vi, we construct a data structure for the annulus with radius ri = d(s, vi)/2ρ and Ri =poly(ρ, n, 1/ε)d(s, vi) s vi ri Ri q We handle the neighborhood of s through scaling

slide-47
SLIDE 47

Overview

For each vertex vi, we construct a data structure for the annulus with radius ri = d(s, vi)/2ρ and Ri =poly(ρ, n, 1/ε)d(s, vi) s vi ri Ri q We handle the neighborhood of s through scaling We handle the empty space (if any) through scaling and perturbation

slide-48
SLIDE 48

Close range queries s

v1 q

slide-49
SLIDE 49

Close range queries s

v1 q q′

slide-50
SLIDE 50

Close range queries s

v1 q q′

slide-51
SLIDE 51

Close range queries s

v1 q q′

slide-52
SLIDE 52

Between two annuli s

q

slide-53
SLIDE 53

Between two annuli s

q

slide-54
SLIDE 54

Between two annuli s

q q′ Perturbation ϕ maps black points → blue points and q → q′

slide-55
SLIDE 55

Between two annuli s

q q′ Find approx shortest path in data structure for the blue subdivision

slide-56
SLIDE 56

Between two annuli s

q q′ Find approx shortest path in data structure for the blue subdivision Clip the path

slide-57
SLIDE 57

Between two annuli s

q q′ Find approx shortest path in data structure for the blue subdivision Clip the path Apply ϕ−1

slide-58
SLIDE 58

Between two annuli s

q q′ Find approx shortest path in data structure for the blue subdivision Clip the path Apply ϕ−1 Go straight to s

slide-59
SLIDE 59

Conclusion

We can also handle obstacles

slide-60
SLIDE 60

Conclusion

We can also handle obstacles 2 point queries?

slide-61
SLIDE 61

Conclusion

We can also handle obstacles 2 point queries?