Mercurial > gemma
annotate pkg/mesh/polygon.go @ 5520:05db984d3db1
Improve performance of bottleneck area calculation
Avoid buffer calculations by replacing them with simple distance comparisons
and calculate the boundary of the result geometry only once per iteration.
In some edge cases with very large numbers of iterations, this reduced
the runtime of a bottleneck import by a factor of more than twenty.
author | Tom Gottfried <tom@intevation.de> |
---|---|
date | Thu, 21 Oct 2021 19:50:39 +0200 |
parents | 5f47eeea988d |
children | 1222b777f51f |
rev | line source |
---|---|
2488
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
1 // This is Free Software under GNU Affero General Public License v >= 3.0 |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
2 // without warranty, see README.md and license for details. |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
3 // |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
4 // SPDX-License-Identifier: AGPL-3.0-or-later |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
5 // License-Filename: LICENSES/AGPL-3.0.txt |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
6 // |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
7 // Copyright (C) 2018 by via donau |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
8 // – Österreichische Wasserstraßen-Gesellschaft mbH |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
9 // Software engineering by Intevation GmbH |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
10 // |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
11 // Author(s): |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
12 // * Sascha L. Teichmann <sascha.teichmann@intevation.de> |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
13 |
4827
f4abfd0ee8ad
Renamed octree package to mesh as there is no octree any more.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
4769
diff
changeset
|
14 package mesh |
2488
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
15 |
2489
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
16 import ( |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
17 "bytes" |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
18 "encoding/binary" |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
19 "fmt" |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
20 "math" |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
21 |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
22 "github.com/tidwall/rtree" |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
23 |
5490
5f47eeea988d
Use own logging package.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
5320
diff
changeset
|
24 "gemma.intevation.de/gemma/pkg/log" |
2489
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
25 "gemma.intevation.de/gemma/pkg/wkb" |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
26 ) |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
27 |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
28 type ( |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
29 ring []float64 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
30 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
31 Polygon struct { |
2574
2833ff156cb2
Morphological differences: Moved loading of clipping polygon into octree package.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2498
diff
changeset
|
32 rings []ring |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
33 indices []*rtree.RTree |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
34 } |
2488
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
35 |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
36 IntersectionType byte |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
37 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
38 lineSegment []float64 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
39 ) |
2488
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
40 |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
41 const ( |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
42 IntersectionInside IntersectionType = iota |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
43 IntersectionOutSide |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
44 IntersectionOverlaps |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
45 ) |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
46 |
5320
866eae1bd888
Upgraded the dependencies of the Gemma server.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
4827
diff
changeset
|
47 func (ls lineSegment) Rect() ([2]float64, [2]float64) { |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
48 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
49 var min, max [2]float64 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
50 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
51 if ls[0] < ls[2] { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
52 min[0] = ls[0] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
53 max[0] = ls[2] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
54 } else { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
55 min[0] = ls[2] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
56 max[0] = ls[0] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
57 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
58 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
59 if ls[1] < ls[3] { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
60 min[1] = ls[1] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
61 max[1] = ls[3] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
62 } else { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
63 min[1] = ls[3] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
64 max[1] = ls[1] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
65 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
66 |
5320
866eae1bd888
Upgraded the dependencies of the Gemma server.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
4827
diff
changeset
|
67 return min, max |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
68 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
69 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
70 func (p *Polygon) Indexify() { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
71 indices := make([]*rtree.RTree, len(p.rings)) |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
72 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
73 for i := range indices { |
5320
866eae1bd888
Upgraded the dependencies of the Gemma server.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
4827
diff
changeset
|
74 index := new(rtree.RTree) |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
75 indices[i] = index |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
76 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
77 rng := p.rings[i] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
78 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
79 for i := 0; i < len(rng); i += 2 { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
80 var ls lineSegment |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
81 if i+4 <= len(rng) { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
82 ls = lineSegment(rng[i : i+4]) |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
83 } else { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
84 ls = []float64{rng[i], rng[i+1], rng[0], rng[1]} |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
85 } |
5320
866eae1bd888
Upgraded the dependencies of the Gemma server.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
4827
diff
changeset
|
86 min, max := ls.Rect() |
866eae1bd888
Upgraded the dependencies of the Gemma server.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
4827
diff
changeset
|
87 index.Insert(min, max, ls) |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
88 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
89 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
90 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
91 p.indices = indices |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
92 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
93 |
2494
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
94 func (ls lineSegment) intersects(a Box2D) bool { |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
95 p1x := ls[0] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
96 p1y := ls[1] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
97 p2x := ls[2] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
98 p2y := ls[3] |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
99 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
100 left := a.X1 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
101 right := a.X2 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
102 top := a.Y1 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
103 bottom := a.Y2 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
104 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
105 // The direction of the ray |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
106 dx := p2x - p1x |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
107 dy := p2y - p1y |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
108 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
109 min, max := 0.0, 1.0 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
110 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
111 var t0, t1 float64 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
112 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
113 // Left and right sides. |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
114 // - If the line is parallel to the y axis. |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
115 if dx == 0 { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
116 if p1x < left || p1x > right { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
117 return false |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
118 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
119 } else { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
120 // - Make sure t0 holds the smaller value by checking the direction of the line. |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
121 if dx > 0 { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
122 t0 = (left - p1x) / dx |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
123 t1 = (right - p1x) / dx |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
124 } else { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
125 t1 = (left - p1x) / dx |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
126 t0 = (right - p1x) / dx |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
127 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
128 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
129 if t0 > min { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
130 min = t0 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
131 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
132 if t1 < max { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
133 max = t1 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
134 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
135 if min > max || max < 0 { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
136 return false |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
137 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
138 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
139 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
140 // The top and bottom side. |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
141 // - If the line is parallel to the x axis. |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
142 if dy == 0 { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
143 if p1y < top || p1y > bottom { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
144 return false |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
145 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
146 } else { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
147 // - Make sure t0 holds the smaller value by checking the direction of the line. |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
148 if dy > 0 { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
149 t0 = (top - p1y) / dy |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
150 t1 = (bottom - p1y) / dy |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
151 } else { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
152 t1 = (top - p1y) / dy |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
153 t0 = (bottom - p1y) / dy |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
154 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
155 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
156 if t0 > min { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
157 min = t0 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
158 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
159 if t1 < max { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
160 max = t1 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
161 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
162 if min > max || max < 0 { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
163 return false |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
164 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
165 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
166 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
167 // The point of intersection |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
168 // ix = p1x + dx*min |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
169 // iy = p1y + dy*min |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
170 return true |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
171 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
172 |
2494
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
173 func (ls lineSegment) intersectsLineSegment(o lineSegment) bool { |
2495
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
174 |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
175 p0 := ls[:2] |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
176 p1 := ls[2:4] |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
177 p2 := o[:2] |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
178 p3 := o[2:4] |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
179 |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
180 s10x := p1[0] - p0[0] |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
181 s10y := p1[1] - p0[1] |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
182 s32x := p3[0] - p2[0] |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
183 s32y := p3[1] - p2[1] |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
184 |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
185 den := s10x*s32y - s32x*s10y |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
186 |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
187 if den == 0 { |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
188 return false |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
189 } |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
190 |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
191 denPos := den > 0 |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
192 |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
193 s02x := p0[0] - p2[0] |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
194 s02y := p0[1] - p2[1] |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
195 |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
196 sNum := s10x*s02y - s10y*s02x |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
197 if sNum < 0 == denPos { |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
198 return false |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
199 } |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
200 |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
201 tNum := s32x*s02y - s32y*s02x |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
202 if tNum < 0 == denPos { |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
203 return false |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
204 } |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
205 |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
206 if sNum > den == denPos || tNum > den == denPos { |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
207 return false |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
208 } |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
209 |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
210 // t := tNum / den |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
211 // intersection at( p0[0] + (t * s10x), p0[1] + (t * s10y) ) |
98bc023750cf
Impelemented line segment / line segment intersection test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2494
diff
changeset
|
212 return true |
2494
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
213 } |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
214 |
2488
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
215 func (p *Polygon) IntersectionBox2D(box Box2D) IntersectionType { |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
216 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
217 if len(p.rings) == 0 { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
218 return IntersectionOutSide |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
219 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
220 |
5320
866eae1bd888
Upgraded the dependencies of the Gemma server.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
4827
diff
changeset
|
221 min, max := box.Rect() |
866eae1bd888
Upgraded the dependencies of the Gemma server.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
4827
diff
changeset
|
222 |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
223 for _, index := range p.indices { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
224 var intersects bool |
5320
866eae1bd888
Upgraded the dependencies of the Gemma server.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
4827
diff
changeset
|
225 index.Search(min, max, func(_, _ [2]float64, item interface{}) bool { |
2494
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
226 if item.(lineSegment).intersects(box) { |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
227 intersects = true |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
228 return false |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
229 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
230 return true |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
231 }) |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
232 if intersects { |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
233 return IntersectionOverlaps |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
234 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
235 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
236 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
237 // No intersection -> check inside or outside |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
238 // if an abritrary point is inside or not. |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
239 |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
240 // Check holes first: inside a hole means outside. |
2498
843f39b9327e
Handle hole more precisely.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2495
diff
changeset
|
241 if len(p.rings) > 1 { |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
242 for _, hole := range p.rings[1:] { |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
243 if contains(hole, box.X1, box.Y1) { |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
244 return IntersectionOutSide |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
245 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
246 } |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
247 } |
2494
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
248 |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
249 // Check shell |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
250 if contains(p.rings[0], box.X1, box.Y1) { |
2492
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
251 return IntersectionInside |
10681749371d
Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2490
diff
changeset
|
252 } |
2488
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
253 return IntersectionOutSide |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
254 } |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
255 |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
256 func (p *Polygon) IntersectionWithTriangle(t *Triangle) IntersectionType { |
2494
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
257 box := t.BBox() |
5320
866eae1bd888
Upgraded the dependencies of the Gemma server.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
4827
diff
changeset
|
258 min, max := box.Rect() |
2494
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
259 for _, index := range p.indices { |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
260 var intersects bool |
5320
866eae1bd888
Upgraded the dependencies of the Gemma server.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
4827
diff
changeset
|
261 index.Search(min, max, func(_, _ [2]float64, item interface{}) bool { |
2494
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
262 ls := item.(lineSegment) |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
263 other := make(lineSegment, 4) |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
264 for i := range t { |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
265 other[0] = t[i].X |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
266 other[1] = t[i].Y |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
267 other[2] = t[(i+1)%len(t)].X |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
268 other[3] = t[(i+1)%len(t)].Y |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
269 if ls.intersectsLineSegment(other) { |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
270 intersects = true |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
271 return false |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
272 } |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
273 } |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
274 return true |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
275 }) |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
276 if intersects { |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
277 return IntersectionOverlaps |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
278 } |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
279 } |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
280 // No intersection -> check inside or outside |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
281 // if an abritrary point is inside or not. |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
282 pX, pY := t[0].X, t[0].Y |
2494
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
283 |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
284 // Check holes first: inside a hole means outside. |
2498
843f39b9327e
Handle hole more precisely.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2495
diff
changeset
|
285 if len(p.rings) > 1 { |
2494
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
286 for _, hole := range p.rings[1:] { |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
287 if contains(hole, pX, pY) { |
2494
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
288 return IntersectionOutSide |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
289 } |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
290 } |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
291 } |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
292 |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
293 // Check shell |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
294 if contains(p.rings[0], pX, pY) { |
2494
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
295 return IntersectionInside |
a727e0426240
More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2492
diff
changeset
|
296 } |
2488
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
297 return IntersectionOutSide |
cb55d7eaaa36
Started with the idea to clip an octree by an bounding polygon.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
298 } |
2489
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
299 |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
300 func (rng ring) length() int { |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
301 return len(rng) / 2 |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
302 } |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
303 |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
304 func (rng ring) point(i int) (float64, float64) { |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
305 i *= 2 |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
306 return rng[i], rng[i+1] |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
307 } |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
308 |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
309 type segments interface { |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
310 length() int |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
311 point(int) (float64, float64) |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
312 } |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
313 |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
314 func contains(s segments, pX, pY float64) bool { |
4764
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
315 |
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
316 n := s.length() |
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
317 if n < 3 { |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
318 return false |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
319 } |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
320 |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
321 sX, sY := s.point(0) |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
322 eX, eY := s.point(n - 1) |
2498
843f39b9327e
Handle hole more precisely.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2495
diff
changeset
|
323 |
4764
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
324 const eps = 0.0000001 |
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
325 |
4765
db6c2955ee31
Fixed logical error in detecting if a ring is closed or not.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
4764
diff
changeset
|
326 if math.Abs(sX-eX) > eps || math.Abs(sY-eY) > eps { |
4764
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
327 // It's not closed! |
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
328 return false |
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
329 } |
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
330 |
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
331 var inside bool |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
332 |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
333 for i := 1; i < n; i++ { |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
334 eX, eY := s.point(i) |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
335 if intersectsWithRaycast(pX, pY, sX, sY, eX, eY) { |
4764
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
336 inside = !inside |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
337 } |
4764
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
338 sX, sY = eX, eY |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
339 } |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
340 |
4764
5c80a33edd44
Fixed handling of none-closed polygons in containment test.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4729
diff
changeset
|
341 return inside |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
342 } |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
343 |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
344 // Using the raycast algorithm, this returns whether or not the passed in point |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
345 // Intersects with the edge drawn by the passed in start and end points. |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
346 // Original implementation: http://rosettacode.org/wiki/Ray-casting_algorithm#Go |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
347 func intersectsWithRaycast(pX, pY, sX, sY, eX, eY float64) bool { |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
348 |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
349 // Always ensure that the the first point |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
350 // has a y coordinate that is less than the second point |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
351 if sY > eY { |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
352 // Switch the points if otherwise. |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
353 sX, sY, eX, eY = eX, eY, sX, sY |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
354 } |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
355 |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
356 // Move the point's y coordinate |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
357 // outside of the bounds of the testing region |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
358 // so we can start drawing a ray |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
359 for pY == sY || pY == eY { |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
360 pY = math.Nextafter(pY, math.Inf(1)) |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
361 } |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
362 |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
363 // If we are outside of the polygon, indicate so. |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
364 if pY < sY || pY > eY { |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
365 return false |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
366 } |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
367 |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
368 if sX > eX { |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
369 if pX > sX { |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
370 return false |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
371 } |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
372 if pX < eX { |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
373 return true |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
374 } |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
375 } else { |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
376 if pX > eX { |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
377 return false |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
378 } |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
379 if pX < sX { |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
380 return true |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
381 } |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
382 } |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
383 |
4729
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
384 raySlope := (pY - sY) / (pX - sX) |
1137c5a18242
Virtualized point in polygon test with an interface to be usable for contours, too.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
4706
diff
changeset
|
385 diagSlope := (eY - sY) / (eX - sX) |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
386 |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
387 return raySlope >= diagSlope |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
388 } |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
389 |
3658
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
390 func (p *Polygon) NumVertices(ring int) int { |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
391 if ring < 0 || ring >= len(p.rings) { |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
392 return 0 |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
393 } |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
394 return len(p.rings[ring]) / 2 |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
395 } |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
396 |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
397 func (p *Polygon) Vertices(ring int, fn func(float64, float64)) { |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
398 if ring < 0 || ring >= len(p.rings) { |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
399 return |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
400 } |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
401 rng := p.rings[ring] |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
402 for i := 0; i < len(rng); i += 2 { |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
403 fn(rng[i+0], rng[i+1]) |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
404 } |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
405 } |
1c3df921361d
Handle th case that a boundary polygon is uploaded along side with the single beam scan.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3650
diff
changeset
|
406 |
3650
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
407 func (p *Polygon) AsWKB() []byte { |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
408 |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
409 size := 1 + 4 + 4 |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
410 for _, r := range p.rings { |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
411 size += 4 + len(r)*8 |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
412 } |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
413 |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
414 buf := bytes.NewBuffer(make([]byte, 0, size)) |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
415 |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
416 binary.Write(buf, binary.LittleEndian, wkb.NDR) |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
417 binary.Write(buf, binary.LittleEndian, wkb.Polygon) |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
418 binary.Write(buf, binary.LittleEndian, uint32(len(p.rings))) |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
419 |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
420 for _, r := range p.rings { |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
421 binary.Write(buf, binary.LittleEndian, uint32(len(r)/2)) |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
422 for i := 0; i < len(r); i += 2 { |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
423 binary.Write(buf, binary.LittleEndian, math.Float64bits(r[i+0])) |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
424 binary.Write(buf, binary.LittleEndian, math.Float64bits(r[i+1])) |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
425 } |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
426 } |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
427 |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
428 return buf.Bytes() |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
429 } |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
430 |
2489
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
431 func (p *Polygon) FromWKB(data []byte) error { |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
432 |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
433 r := bytes.NewReader(data) |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
434 |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
435 endian, err := r.ReadByte() |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
436 |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
437 var order binary.ByteOrder |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
438 |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
439 switch { |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
440 case err != nil: |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
441 return err |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
442 case endian == wkb.NDR: |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
443 order = binary.LittleEndian |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
444 case endian == wkb.XDR: |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
445 order = binary.BigEndian |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
446 default: |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
447 return fmt.Errorf("unknown byte order %x", endian) |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
448 } |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
449 |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
450 var geomType uint32 |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
451 err = binary.Read(r, order, &geomType) |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
452 |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
453 switch { |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
454 case err != nil: |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
455 return err |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
456 case geomType != wkb.Polygon: |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
457 return fmt.Errorf("unknown geometry type %x", geomType) |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
458 } |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
459 |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
460 var numRings uint32 |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
461 if err = binary.Read(r, order, &numRings); err != nil { |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
462 return err |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
463 } |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
464 |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
465 rngs := make([]ring, numRings) |
2489
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
466 |
5490
5f47eeea988d
Use own logging package.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
5320
diff
changeset
|
467 log.Infof("number of rings: %d\n", len(rngs)) |
2498
843f39b9327e
Handle hole more precisely.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2495
diff
changeset
|
468 |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
469 for rng := uint32(0); rng < numRings; rng++ { |
2489
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
470 var numVertices uint32 |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
471 if err = binary.Read(r, order, &numVertices); err != nil { |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
472 return err |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
473 } |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
474 |
5490
5f47eeea988d
Use own logging package.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
5320
diff
changeset
|
475 log.Infof("number of vertices in ring %d: %d\n", rng, numVertices) |
2498
843f39b9327e
Handle hole more precisely.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2495
diff
changeset
|
476 |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
477 numVertices *= 2 |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
478 vertices := make([]float64, numVertices) |
2489
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
479 |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
480 for v := uint32(0); v < numVertices; v += 2 { |
2489
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
481 var lat, lon uint64 |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
482 if err = binary.Read(r, order, &lat); err != nil { |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
483 return err |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
484 } |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
485 if err = binary.Read(r, order, &lon); err != nil { |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
486 return err |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
487 } |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
488 vertices[v] = math.Float64frombits(lat) |
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
489 vertices[v+1] = math.Float64frombits(lon) |
2489
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
490 } |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
491 |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
492 rngs[rng] = vertices |
2489
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
493 } |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
494 |
2490
c9164ff98871
Started with point in polygon check.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2489
diff
changeset
|
495 p.rings = rngs |
2489
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
496 |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
497 return nil |
4f292ff74d9e
Added a WKB reader for clipping polygon.
Sascha L. Teichmann <teichmann@intevation.de>
parents:
2488
diff
changeset
|
498 } |
3650
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
499 |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
500 func (p *Polygon) FromLineStringZ(ls LineStringZ) { |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
501 r := make([]float64, 2*len(ls)) |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
502 var pos int |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
503 for i := range ls { |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
504 r[pos+0] = ls[i].X |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
505 r[pos+1] = ls[i].Y |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
506 pos += 2 |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
507 } |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
508 p.rings = []ring{r} |
01ce3ba9b0d0
Fixed generating of random points.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
2587
diff
changeset
|
509 } |