gemma: pkg/mesh/polygon.go annotate

annotate pkg/mesh/polygon.go @ 4827:f4abfd0ee8ad remove-octree-debris

Renamed octree package to mesh as there is no octree any more.

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Tue, 05 Nov 2019 14:30:22 +0100
parents	pkg/octree/polygon.go@ea570f43d7a9
children	866eae1bd888

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"
2498 843f39b9327e Handle hole more precisely. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2495 diff changeset	20 "log"
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	21 "math"
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	22
2492 10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	23 "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	24
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
2492 10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	47 func (ls lineSegment) Rect(interface{}) ([]float64, []float64) {
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
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	67 return min[:], max[:]
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 {
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	74 index := rtree.New(nil)
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 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	86 index.Insert(ls)
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	87 }
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 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	91 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	92
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	93 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	94 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	95 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	96 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	97 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	98
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	99 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	100 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	101 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	102 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	103
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	104 // 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	105 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	106 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	107
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	108 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	109
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	110 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	111
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	112 // 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	113 // - 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	114 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	115 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	116 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	117 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	118 } else {
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	119 // - 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	120 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	121 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	122 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	123 } else {
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	124 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	125 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	126 }
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 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	129 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	130 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	131 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	132 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	133 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	134 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	135 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	136 }
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 // 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	140 // - 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	141 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	142 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	143 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	144 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	145 } else {
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	146 // - 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	147 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	148 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	149 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	150 } else {
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	151 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	152 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	153 }
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 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	156 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	157 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	158 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	159 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	160 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	161 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	162 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	163 }
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 // 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	167 // 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	168 // 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	169 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	170 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	171
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	172 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	173
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	174 p0 := ls[:2]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	175 p1 := ls[2:4]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	176 p2 := o[:2]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	177 p3 := o[2:4]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	178
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	179 s10x := p1[0] - p0[0]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	180 s10y := p1[1] - p0[1]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	181 s32x := p3[0] - p2[0]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	182 s32y := p3[1] - p2[1]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	183
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	184 den := s10xs32y - s32xs10y
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	185
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	186 if den == 0 {
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	187 return false
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	188 }
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 denPos := den > 0
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	191
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	192 s02x := p0[0] - p2[0]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	193 s02y := p0[1] - p2[1]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	194
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	195 sNum := s10xs02y - s10ys02x
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	196 if sNum < 0 == denPos {
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	197 return false
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	198 }
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 tNum := s32xs02y - s32ys02x
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	201 if tNum < 0 == denPos {
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	202 return false
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	203 }
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 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	206 return false
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	207 }
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 // t := tNum / den
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	210 // 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	211 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	212 }
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
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214 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	215
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	216 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	217 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	218 }
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 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	221 var intersects bool
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	222 index.Search(box, func(item rtree.Item) 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	223 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	224 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	225 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	226 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	227 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	228 })
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	229 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	230 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	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 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	233
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	234 // 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	235 // 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	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 // 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	238 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	239 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	240 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	241 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	242 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	243 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	244 }
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	245
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	246 // 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	247 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	248 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	249 }
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	250 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	251 }
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	252
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	253 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	254 box := t.BBox()
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	255 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	256 var intersects bool
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 index.Search(box, func(item rtree.Item) bool {
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	258 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	259 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	260 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	261 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	262 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	263 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	264 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	265 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	266 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	267 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	268 }
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 }
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 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	271 })
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 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	273 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	274 }
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 // 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	277 // 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	278 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	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 // 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	281 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	282 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	283 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	284 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	285 }
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 }
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	287 }
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
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 // 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	290 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	291 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	292 }
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	293 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	294 }
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	295
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	296 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	297 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	298 }
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 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) 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	301 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	302 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	303 }
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
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 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	306 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	307 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	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
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 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	311
5c80a33edd44 Fixed handling of none-closed polygons in containment test. Sascha L. Teichmann <teichmann@intevation.de> parents: 4729 diff changeset	312 n := s.length()
5c80a33edd44 Fixed handling of none-closed polygons in containment test. Sascha L. Teichmann <teichmann@intevation.de> parents: 4729 diff changeset	313 if n < 3 {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	314 return false
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	315 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	316
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	317 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	318 eX, eY := s.point(n - 1)
2498 843f39b9327e Handle hole more precisely. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2495 diff changeset	319
4764 5c80a33edd44 Fixed handling of none-closed polygons in containment test. Sascha L. Teichmann <teichmann@intevation.de> parents: 4729 diff changeset	320 const eps = 0.0000001
5c80a33edd44 Fixed handling of none-closed polygons in containment test. Sascha L. Teichmann <teichmann@intevation.de> parents: 4729 diff changeset	321
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	322 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	323 // It's not closed!
5c80a33edd44 Fixed handling of none-closed polygons in containment test. Sascha L. Teichmann <teichmann@intevation.de> parents: 4729 diff changeset	324 return false
5c80a33edd44 Fixed handling of none-closed polygons in containment test. Sascha L. Teichmann <teichmann@intevation.de> parents: 4729 diff changeset	325 }
5c80a33edd44 Fixed handling of none-closed polygons in containment test. Sascha L. Teichmann <teichmann@intevation.de> parents: 4729 diff changeset	326
5c80a33edd44 Fixed handling of none-closed polygons in containment test. Sascha L. Teichmann <teichmann@intevation.de> parents: 4729 diff changeset	327 var inside bool
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	328
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	329 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	330 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	331 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	332 inside = !inside
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	333 }
4764 5c80a33edd44 Fixed handling of none-closed polygons in containment test. Sascha L. Teichmann <teichmann@intevation.de> parents: 4729 diff changeset	334 sX, sY = eX, eY
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	335 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	336
4764 5c80a33edd44 Fixed handling of none-closed polygons in containment test. Sascha L. Teichmann <teichmann@intevation.de> parents: 4729 diff changeset	337 return inside
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	338 }
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 // 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	341 // 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	342 // 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	343 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	344
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	345 // 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	346 // 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	347 if sY > eY {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	348 // 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	349 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	350 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	351
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	352 // Move the point's y coordinate
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	353 // 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	354 // 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	355 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	356 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	357 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	358
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	359 // 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	360 if pY < sY \|\| pY > eY {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	361 return false
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
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 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	365 if pX > sX {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	366 return false
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 pX < eX {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	369 return true
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	370 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	371 } 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	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 false
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	374 }
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	375 if pX < sX {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	376 return true
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	377 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	378 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	379
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	380 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	381 diagSlope := (eY - sY) / (eX - sX)
2490 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 return raySlope >= diagSlope
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	384 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	385
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	386 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	387 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	388 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	389 }
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 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	391 }
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
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 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	394 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	395 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	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 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	398 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	399 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	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 }
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
3650 01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	403 func (p *Polygon) AsWKB() []byte {
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	404
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	405 size := 1 + 4 + 4
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	406 for _, r := range p.rings {
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	407 size += 4 + len(r)*8
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
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	410 buf := bytes.NewBuffer(make([]byte, 0, size))
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	411
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	412 binary.Write(buf, binary.LittleEndian, wkb.NDR)
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	413 binary.Write(buf, binary.LittleEndian, wkb.Polygon)
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	414 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	415
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	416 for _, r := range p.rings {
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	417 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	418 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	419 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	420 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	421 }
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	422 }
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	423
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	424 return buf.Bytes()
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
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	427 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	428
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	429 r := bytes.NewReader(data)
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	430
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	431 endian, err := r.ReadByte()
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 var order binary.ByteOrder
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 switch {
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	436 case err != nil:
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	437 return err
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	438 case endian == wkb.NDR:
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	439 order = binary.LittleEndian
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	440 case endian == wkb.XDR:
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	441 order = binary.BigEndian
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	442 default:
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	443 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	444 }
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	445
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	446 var geomType uint32
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	447 err = binary.Read(r, order, &geomType)
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 switch {
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	450 case err != nil:
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	451 return err
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	452 case geomType != wkb.Polygon:
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	453 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	454 }
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	455
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	456 var numRings uint32
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	457 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	458 return err
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
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	461 rngs := make([]ring, numRings)
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	462
2498 843f39b9327e Handle hole more precisely. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2495 diff changeset	463 log.Printf("info: Number of rings: %d\n", len(rngs))
843f39b9327e Handle hole more precisely. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2495 diff changeset	464
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	465 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	466 var numVertices uint32
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	467 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	468 return err
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	469 }
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	470
2498 843f39b9327e Handle hole more precisely. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2495 diff changeset	471 log.Printf("info: Number of vertices in ring %d: %d\n", rng, numVertices)
843f39b9327e Handle hole more precisely. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2495 diff changeset	472
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	473 numVertices *= 2
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	474 vertices := make([]float64, numVertices)
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	475
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	476 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	477 var lat, lon uint64
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	478 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	479 return err
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	480 }
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	481 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	482 return err
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	483 }
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	484 vertices[v] = math.Float64frombits(lat)
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	485 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	486 }
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 rngs[rng] = vertices
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	489 }
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	490
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	491 p.rings = rngs
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	492
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	493 return nil
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	494 }
3650 01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	495
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	496 func (p *Polygon) FromLineStringZ(ls LineStringZ) {
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	497 r := make([]float64, 2*len(ls))
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	498 var pos int
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	499 for i := range ls {
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	500 r[pos+0] = ls[i].X
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	501 r[pos+1] = ls[i].Y
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	502 pos += 2
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	503 }
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	504 p.rings = []ring{r}
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	505 }

Mercurial > gemma

annotate pkg/mesh/polygon.go @ 4827:f4abfd0ee8ad remove-octree-debris