gemma: pkg/mesh/polygon.go annotate

annotate pkg/mesh/polygon.go @ 5591:0011f50cf216 surveysperbottleneckid

Removed no longer used alternative api for surveys/ endpoint. As bottlenecks in the summary for SR imports are now identified by their id and no longer by the (not guarantied to be unique!) name, there is no longer the need to request survey data by the name+date tuple (which isn't reliable anyway). So the workaround was now reversed.

author	Sascha Wilde <wilde@sha-bang.de>
date	Wed, 06 Apr 2022 13:30:29 +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 := s10xs32y - s32xs10y
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 := s10xs02y - s10ys02x
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 := s32xs02y - s32ys02x
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 }

Mercurial > gemma

annotate pkg/mesh/polygon.go @ 5591:0011f50cf216 surveysperbottleneckid