gemma: pkg/octree/polygon.go annotate

annotate pkg/octree/polygon.go @ 4729:1137c5a18242 stack-polygons

Virtualized point in polygon test with an interface to be usable for contours, too.

author	Sascha L. Teichmann <teichmann@intevation.de>
date	Thu, 17 Oct 2019 23:33:20 +0200
parents	a38d846d9fd5
children	5c80a33edd44

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
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	14 package octree
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"
2574 2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	18 "context"
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	19 "database/sql"
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	20 "encoding/binary"
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	21 "fmt"
2498 843f39b9327e Handle hole more precisely. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2495 diff changeset	22 "log"
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	23 "math"
2574 2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	24 "time"
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	25
2492 10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	26 "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	27
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	28 "gemma.intevation.de/gemma/pkg/wkb"
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	29 )
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	30
2492 10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	31 type (
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	32 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	33
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	34 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	35 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	36 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	37 }
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38
2492 10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	39 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	40
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	41 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	42 )
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	43
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	44 const (
2574 2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	45 clippingPolygonSQL = `
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	46 WITH joined AS (
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	47 SELECT
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	48 sr.area AS area,
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	49 sr.date_info AS date_info
2587 857bb070b9f1 Morphological differences: Bottleneck_id is an integer represented as a string. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2574 diff changeset	50 FROM waterway.sounding_results sr
857bb070b9f1 Morphological differences: Bottleneck_id is an integer represented as a string. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2574 diff changeset	51 WHERE sr.bottleneck_id = $1
2574 2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	52 )
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	53 SELECT ST_AsBinary(
4706 a38d846d9fd5 morphology difference: buffer common boundary polygon by 10cm to fix numerical problems leading to unexpected holes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3658 diff changeset	54 ST_Buffer(ST_intersection(
2574 2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	55 (SELECT ST_Transform(area::geometry, $2::int) FROM joined WHERE date_info = $3::date),
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	56 (SELECT ST_Transform(area::geometry, $2::int) FROM joined WHERE date_info = $4::date)
4706 a38d846d9fd5 morphology difference: buffer common boundary polygon by 10cm to fix numerical problems leading to unexpected holes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3658 diff changeset	57 ), 0.1)
2574 2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	58 ) AS area
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	59 `
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	60 )
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	61
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	62 const (
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	63 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	64 IntersectionOutSide
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	65 IntersectionOverlaps
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	66 )
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	67
2574 2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	68 func LoadClippingPolygon(
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	69 ctx context.Context,
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	70 conn *sql.Conn,
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	71 epsg uint32,
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	72 bottleneck string,
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	73 first, second time.Time,
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	74 ) (*Polygon, error) {
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	75
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	76 var clip []byte
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	77
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	78 if err := conn.QueryRowContext(
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	79 ctx, clippingPolygonSQL,
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	80 bottleneck,
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	81 epsg,
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	82 first, second,
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	83 ).Scan(&clip); err != nil {
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	84 return nil, err
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	85 }
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	86
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	87 var polygon Polygon
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	88 if err := polygon.FromWKB(clip); err != nil {
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	89 return nil, err
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	90 }
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	91 polygon.Indexify()
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	92 return &polygon, nil
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	93 }
2833ff156cb2 Morphological differences: Moved loading of clipping polygon into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2498 diff changeset	94
2492 10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	95 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	96
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	97 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	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 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	100 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	101 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	102 } else {
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	103 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	104 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	105 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	106
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	107 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	108 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	109 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	110 } else {
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	111 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	112 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	113 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	114
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	115 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	116 }
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 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	119 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	120
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	121 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	122 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	123 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	124
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	125 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	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 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	128 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	129 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	130 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	131 } else {
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	132 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	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 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	135 }
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 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	139 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	140
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	141 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	142 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	143 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	144 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	145 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	146
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	147 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	148 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	149 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	150 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	151
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	152 // 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	153 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	154 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	155
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	156 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	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 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	159
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	160 // 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	161 // - 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	162 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	163 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	164 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	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 } else {
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	167 // - 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	168 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	169 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	170 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	171 } else {
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	172 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	173 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	174 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	175
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	176 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	177 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	178 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	179 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	180 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	181 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	182 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	183 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	184 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	185 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	186
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	187 // 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	188 // - 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	189 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	190 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	191 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	192 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	193 } else {
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	194 // - 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	195 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	196 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	197 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	198 } else {
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	199 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	200 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	201 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	202
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	203 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	204 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	205 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	206 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	207 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	208 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	209 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	210 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	211 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	212 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	213
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	214 // 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	215 // 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	216 // 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	217 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	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
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	220 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	221
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	222 p0 := ls[:2]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	223 p1 := ls[2:4]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	224 p2 := o[:2]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	225 p3 := o[2:4]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	226
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	227 s10x := p1[0] - p0[0]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	228 s10y := p1[1] - p0[1]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	229 s32x := p3[0] - p2[0]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	230 s32y := p3[1] - p2[1]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	231
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	232 den := s10xs32y - s32xs10y
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	233
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	234 if den == 0 {
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	235 return false
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	236 }
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	237
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	238 denPos := den > 0
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	239
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	240 s02x := p0[0] - p2[0]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	241 s02y := p0[1] - p2[1]
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	242
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	243 sNum := s10xs02y - s10ys02x
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	244 if sNum < 0 == denPos {
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	245 return false
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	246 }
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	247
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	248 tNum := s32xs02y - s32ys02x
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	249 if tNum < 0 == denPos {
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	250 return false
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	251 }
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	252
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	253 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	254 return false
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	255 }
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	256
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	257 // t := tNum / den
98bc023750cf Impelemented line segment / line segment intersection test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	258 // 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	259 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	260 }
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
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	262 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	263
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	264 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	265 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	266 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	267
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	268 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	269 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	270 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	271 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	272 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	273 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	274 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	275 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	276 })
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	277 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	278 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	279 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	280 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	281
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	282 // 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	283 // 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	284
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	285 // 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	286 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	287 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	288 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	289 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	290 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	291 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2490 diff changeset	292 }
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	293
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	294 // 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	295 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	296 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	297 }
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	298 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	299 }
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	300
cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	301 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	302 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	303 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	304 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	305 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	306 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	307 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	308 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	309 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	310 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	311 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	312 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	313 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	314 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	315 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	316 }
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	317 }
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	318 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	319 })
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	320 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	321 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	322 }
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	323 }
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	324 // 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	325 // 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	326 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	327
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	328 // 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	329 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	330 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	331 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	332 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	333 }
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	334 }
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	335 }
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	336
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	337 // 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	338 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	339 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	340 }
2488 cb55d7eaaa36 Started with the idea to clip an octree by an bounding polygon. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	341 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	342 }
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	343
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	344 func (rng ring) closed() bool {
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	345 return (len(rng) / 2) >= 3
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	346 }
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
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	348 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	349 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	350 }
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	351
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	352 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	353 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	354 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	355 }
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
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	357 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	358 closed() bool
1137c5a18242 Virtualized point in polygon test with an interface to be usable for contours, too. Sascha L. Teichmann <teichmann@intevation.de> parents: 4706 diff changeset	359 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	360 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	361 }
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	362
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	363 func contains(s segments, pX, pY float64) bool {
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 !s.closed() {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	365 return false
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	366 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	367
4729 1137c5a18242 Virtualized point in polygon test with an interface to be usable for contours, too. Sascha L. Teichmann <teichmann@intevation.de> parents: 4706 diff changeset	368 n := s.length()
1137c5a18242 Virtualized point in polygon test with an interface to be usable for contours, too. Sascha L. Teichmann <teichmann@intevation.de> parents: 4706 diff changeset	369
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	370 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	371 eX, eY := s.point(n - 1)
2498 843f39b9327e Handle hole more precisely. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2495 diff changeset	372
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	373 contains := intersectsWithRaycast(pX, pY, sX, sY, eX, eY)
2490 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 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	376 sX, sY := s.point(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	377 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	378 if intersectsWithRaycast(pX, pY, sX, sY, eX, eY) {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	379 contains = !contains
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	380 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	381 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	382
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	383 return contains
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
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	386 // 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	387 // 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	388 // 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	389 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	390
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	391 // 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	392 // 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	393 if sY > eY {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	394 // 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	395 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	396 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	397
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	398 // Move the point's y coordinate
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	399 // 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	400 // 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	401 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	402 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	403 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	404
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	405 // 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	406 if pY < sY \|\| pY > eY {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	407 return false
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	408 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	409
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	410 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	411 if pX > sX {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	412 return false
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	413 }
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	414 if pX < eX {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	415 return true
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	416 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	417 } 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	418 if pX > eX {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	419 return false
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	420 }
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	421 if pX < sX {
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	422 return true
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	423 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	424 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	425
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	426 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	427 diagSlope := (eY - sY) / (eX - sX)
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	428
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	429 return raySlope >= diagSlope
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	430 }
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	431
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	432 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	433 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	434 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	435 }
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	436 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	437 }
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	438
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	439 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	440 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	441 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	442 }
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	443 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	444 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	445 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	446 }
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	447 }
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	448
3650 01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	449 func (p *Polygon) AsWKB() []byte {
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	450
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	451 size := 1 + 4 + 4
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	452 for _, r := range p.rings {
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	453 size += 4 + len(r)*8
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	454 }
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	455
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	456 buf := bytes.NewBuffer(make([]byte, 0, size))
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	457
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	458 binary.Write(buf, binary.LittleEndian, wkb.NDR)
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	459 binary.Write(buf, binary.LittleEndian, wkb.Polygon)
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	460 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	461
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	462 for _, r := range p.rings {
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	463 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	464 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	465 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	466 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	467 }
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	468 }
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	469
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	470 return buf.Bytes()
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	471 }
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	472
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	473 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	474
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	475 r := bytes.NewReader(data)
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	476
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	477 endian, err := r.ReadByte()
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	478
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	479 var order binary.ByteOrder
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 switch {
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	482 case err != nil:
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	483 return err
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	484 case endian == wkb.NDR:
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	485 order = binary.LittleEndian
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	486 case endian == wkb.XDR:
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	487 order = binary.BigEndian
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	488 default:
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	489 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	490 }
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	491
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	492 var geomType uint32
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	493 err = binary.Read(r, order, &geomType)
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	494
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	495 switch {
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	496 case err != nil:
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	497 return err
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	498 case geomType != wkb.Polygon:
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	499 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	500 }
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	501
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	502 var numRings uint32
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	503 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	504 return err
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	505 }
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	506
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	507 rngs := make([]ring, numRings)
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	508
2498 843f39b9327e Handle hole more precisely. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2495 diff changeset	509 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	510
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	511 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	512 var numVertices uint32
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	513 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	514 return err
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	515 }
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	516
2498 843f39b9327e Handle hole more precisely. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2495 diff changeset	517 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	518
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	519 numVertices *= 2
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	520 vertices := make([]float64, numVertices)
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	521
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	522 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	523 var lat, lon uint64
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	524 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	525 return err
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	526 }
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	527 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	528 return err
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	529 }
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	530 vertices[v] = math.Float64frombits(lat)
c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	531 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	532 }
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	533
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	534 rngs[rng] = vertices
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	535 }
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	536
2490 c9164ff98871 Started with point in polygon check. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2489 diff changeset	537 p.rings = rngs
2489 4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	538
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	539 return nil
4f292ff74d9e Added a WKB reader for clipping polygon. Sascha L. Teichmann <teichmann@intevation.de> parents: 2488 diff changeset	540 }
3650 01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	541
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	542 func (p *Polygon) FromLineStringZ(ls LineStringZ) {
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	543 r := make([]float64, 2*len(ls))
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	544 var pos int
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	545 for i := range ls {
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	546 r[pos+0] = ls[i].X
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	547 r[pos+1] = ls[i].Y
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	548 pos += 2
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	549 }
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	550 p.rings = []ring{r}
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2587 diff changeset	551 }

Mercurial > gemma

annotate pkg/octree/polygon.go @ 4729:1137c5a18242 stack-polygons