gemma: pkg/mesh/vertex.go annotate

annotate pkg/mesh/vertex.go @ 5710:37c8feeecb4d

Merged branch sr-v2 into default.

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Tue, 20 Feb 2024 21:28:56 +0100
parents	148abae1fcd0
children	61eff98a40ae

rev	line source
1017 a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	1 // This is Free Software under GNU Affero General Public License v >= 3.0
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	2 // without warranty, see README.md and license for details.
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	3 //
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	4 // SPDX-License-Identifier: AGPL-3.0-or-later
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	5 // License-Filename: LICENSES/AGPL-3.0.txt
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	6 //
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	7 // Copyright (C) 2018 by via donau
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	8 // – Österreichische Wasserstraßen-Gesellschaft mbH
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	9 // Software engineering by Intevation GmbH
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	10 //
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	11 // Author(s):
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	12 // * Sascha L. Teichmann <sascha.teichmann@intevation.de>
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 988 diff changeset	13
4827 f4abfd0ee8ad Renamed octree package to mesh as there is no octree any more. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4768 diff changeset	14 package mesh
674 120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	16 import (
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	17 "bytes"
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	18 "encoding/binary"
2529 45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	19 "fmt"
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	20 "io"
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	21 "math"
768 3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	22 "sort"
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	23
5490 5f47eeea988d Use own logging package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5428 diff changeset	24 "gemma.intevation.de/gemma/pkg/log"
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	25 "gemma.intevation.de/gemma/pkg/wkb"
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	26 )
674 120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	27
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	28 type (
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	29 // Point is a 2D point.
2572 7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2529 diff changeset	30 Point struct {
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2529 diff changeset	31 X float64
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2529 diff changeset	32 Y float64
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2529 diff changeset	33 }
7686c7c23506 Morphological differences: Moved some code into octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2529 diff changeset	34
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	35 // Vertex is a 3D vertex.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	36 Vertex struct {
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	37 X float64
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	38 Y float64
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	39 Z float64
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	40 }
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	41
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	42 // Triangle is a triangle consisting of three vertices.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	43 Triangle [3]Vertex
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	44
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	45 // Line is a line defined by first vertex on that line
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	46 // and the second being the direction.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	47 Line [2]Vertex
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	48
3953 4233570de212 Building octrees: Don't cause havok by building ain overly large tree if the points are not really 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3770 diff changeset	49 // Box is a 3D box.
4233570de212 Building octrees: Don't cause havok by building ain overly large tree if the points are not really 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3770 diff changeset	50 Box [2]Vertex
4233570de212 Building octrees: Don't cause havok by building ain overly large tree if the points are not really 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3770 diff changeset	51
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	52 // MultiPointZ is a set of vertices.
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	53 MultiPointZ []Vertex
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	54
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	55 // LineStringZ is a line string formed of vertices.
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	56 LineStringZ []Vertex
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	57
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	58 // MultiLineStringZ is a set of line strings.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	59 MultiLineStringZ []LineStringZ
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	60
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	61 // Box2D is 2D area from (X1, Y1) to (X2, Y2).
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	62 Box2D struct {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	63 X1 float64
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	64 Y1 float64
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	65 X2 float64
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	66 Y2 float64
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	67 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	68
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	69 // Plane2D is a 2D plane (a line in 2D space).
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	70 Plane2D struct {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	71 A float64
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	72 B float64
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	73 C float64
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	74 }
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	75
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	76 // Plane3D is a 3D plane.
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	77 Plane3D struct {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	78 A float64
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	79 B float64
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	80 C float64
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	81 D float64
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	82 }
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	83 )
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	84
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	85 // Plane3D returns the plane in which
b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	86 // the three points of the triangles are in.
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	87 func (t *Triangle) Plane3D() Plane3D {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	88
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	89 v0 := t[1].Sub(t[0])
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	90 v1 := t[2].Sub(t[0])
5703 d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	91 n := v0.Cross(v1)
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	92
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	93 // If the length of normal is near zero assume we have
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	94 // a plane constant in z with an average z value of
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	95 // the three vertices.
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	96 // This should protect us from the effects of collinear
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	97 // geometries.
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	98 l := n.Length()
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	99 if l < 1e-7 {
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	100 sum := t[0].Z + t[1].Z + t[2].Z
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	101 return Plane3D{
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	102 A: 0,
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	103 B: 0,
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	104 C: 1,
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	105 D: -sum / 3,
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	106 }
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	107 }
d2ccf6bb6940 Make plane eval for z-values of triangles numerial more robust. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5601 diff changeset	108 n = n.Scale(1 / l)
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	109
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	110 // xnx+ yny+ z*nz + d = 0
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	111 // d = - (xnx+ yny + z*nz)
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	112 d := -t[0].Dot(n)
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	113 return Plane3D{
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	114 A: n.X,
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	115 B: n.Y,
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	116 C: n.Z,
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	117 D: d,
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	118 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	119 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	120
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	121 // BBox calculates the 2D (in X/Y plane) bounding box
b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	122 // of the triangle.
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	123 func (t *Triangle) BBox() Box2D {
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	124 minX := math.Min(math.Min(t[0].X, t[1].X), t[2].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	125 maxX := math.Max(math.Max(t[0].X, t[1].X), t[2].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	126 minY := math.Min(math.Min(t[0].Y, t[1].Y), t[2].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	127 maxY := math.Max(math.Max(t[0].Y, t[1].Y), t[2].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	128 return Box2D{
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	129 X1: minX, Y1: minY,
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	130 X2: maxX, Y2: maxY,
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	131 }
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	132 }
a727e0426240 More on triangle/clipping polygon intersection. TODO: line segment / line segment intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2492 diff changeset	133
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	134 // Inside test if b is completely inside a.
4726 c91e759007da Started with stacking shells and holes of generated polygons correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4554 diff changeset	135 func (a Box2D) Inside(b Box2D) bool {
c91e759007da Started with stacking shells and holes of generated polygons correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4554 diff changeset	136 return a.X1 >= b.X1 && a.X2 <= b.X2 &&
c91e759007da Started with stacking shells and holes of generated polygons correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4554 diff changeset	137 a.Y1 >= b.Y1 && a.Y2 <= b.Y2
c91e759007da Started with stacking shells and holes of generated polygons correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4554 diff changeset	138 }
c91e759007da Started with stacking shells and holes of generated polygons correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4554 diff changeset	139
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	140 // Size calculates the area of the box.
4768 a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4726 diff changeset	141 func (a Box2D) Size() (float64, float64) {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4726 diff changeset	142 return a.X2 - a.X1, a.Y2 - a.Y1
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4726 diff changeset	143 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4726 diff changeset	144
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	145 // Empty returns true if the box has no area.
4768 a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4726 diff changeset	146 func (a Box2D) Empty() bool {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4726 diff changeset	147 const eps = 0.0000001
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4726 diff changeset	148 return math.Abs(a.X2-a.X1) < eps &&
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4726 diff changeset	149 math.Abs(a.Y2-a.Y1) < eps
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4726 diff changeset	150 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4726 diff changeset	151
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	152 // Z calculates the Z value for a given X/Y value.
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	153 func (p Plane3D) Z(x, y float64) float64 {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	154 // p.Ax + p.By + p.C*z + p.D = 0
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	155 return -(p.Ax + p.By + p.D) / p.C
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	156 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	157
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	158 // Eval evalutes the plane equation for a given vertex.
3770 71164b817d6e Calculate distance from planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3655 diff changeset	159 func (p Plane3D) Eval(v Vertex) float64 {
71164b817d6e Calculate distance from planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3655 diff changeset	160 return p.Av.X + p.Bv.Y + p.C*v.Z + p.D
71164b817d6e Calculate distance from planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3655 diff changeset	161 }
71164b817d6e Calculate distance from planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3655 diff changeset	162
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	163 // Normalize constructs a new vertex with unit length for a given vertex.
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	164 func (v Vertex) Normalize() Vertex {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	165 s := 1 / v.Length()
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	166 return Vertex{
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	167 X: s * v.X,
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	168 Y: s * v.Y,
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	169 Z: s * v.Z,
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	170 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	171 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	172
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	173 // Dot returns the 3D dot product of the two given vertices.
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	174 func (v Vertex) Dot(w Vertex) float64 {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	175 return v.Xw.X + v.Yw.Y + v.Z*w.Z
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	176 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	177
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	178 // Length return the length of the vertex.
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	179 func (v Vertex) Length() float64 {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	180 return math.Sqrt(v.Dot(v))
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	181 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	182
5419 202715173935 Do clipping correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5418 diff changeset	183 // Box2D constructs a Box2D of this vertex.
202715173935 Do clipping correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5418 diff changeset	184 func (v Vertex) Box2D() Box2D {
202715173935 Do clipping correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5418 diff changeset	185 return Box2D{
202715173935 Do clipping correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5418 diff changeset	186 X1: v.X,
202715173935 Do clipping correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5418 diff changeset	187 Y1: v.Y,
202715173935 Do clipping correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5418 diff changeset	188 X2: v.X,
202715173935 Do clipping correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5418 diff changeset	189 Y2: v.Y,
202715173935 Do clipping correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5418 diff changeset	190 }
202715173935 Do clipping correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5418 diff changeset	191 }
202715173935 Do clipping correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5418 diff changeset	192
2483 620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	193 func area(a, b, c Vertex) float64 {
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	194 return (b.Y-a.Y)(c.X-b.X) - (b.X-a.X)(c.Y-b.Y)
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	195 }
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	196
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	197 func inCircle(a, b, c, p Vertex) bool {
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	198 dx := a.X - p.X
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	199 dy := a.Y - p.Y
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	200 ex := b.X - p.X
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	201 ey := b.Y - p.Y
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	202 fx := c.X - p.X
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	203 fy := c.Y - p.Y
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	204
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	205 ap := dxdx + dydy
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	206 bp := exex + eyey
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	207 cp := fxfx + fyfy
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	208
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	209 return dx(eycp-bpfy)-dy(excp-bpfx)+ap(exfy-ey*fx) < 0
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	210 }
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	211
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	212 func circumradius(a, b, c Vertex) float64 {
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	213 dx := b.X - a.X
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	214 dy := b.Y - a.Y
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	215 ex := c.X - a.X
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	216 ey := c.Y - a.Y
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	217
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	218 bl := dxdx + dydy
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	219 cl := exex + eyey
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	220 d := dxey - dyex
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	221
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	222 x := (eybl - dycl) * 0.5 / d
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	223 y := (dxcl - exbl) * 0.5 / d
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	224
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	225 r := xx + yy
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	226
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	227 if bl == 0 \|\| cl == 0 \|\| d == 0 \|\| r == 0 {
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	228 return infinity
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	229 }
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	230
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	231 return r
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	232 }
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	233
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	234 func circumcenter(a, b, c Vertex) Vertex {
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	235 dx := b.X - a.X
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	236 dy := b.Y - a.Y
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	237 ex := c.X - a.X
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	238 ey := c.Y - a.Y
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	239
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	240 bl := dxdx + dydy
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	241 cl := exex + eyey
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	242 d := dxey - dyex
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	243
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	244 x := a.X + (eybl-dycl)*0.5/d
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	245 y := a.Y + (dxcl-exbl)*0.5/d
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	246
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	247 return Vertex{X: x, Y: y}
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	248 }
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	249
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	250 func polygonArea(points []Vertex) float64 {
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	251 var result float64
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	252 for i, p := range points {
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	253 q := points[(i+1)%len(points)]
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	254 result += (p.X - q.X) * (p.Y + q.Y)
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	255 }
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	256 return result / 2
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	257 }
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	258
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	259 // Distance2D returns the distance of the two vertices in the X/Y plane.
2483 620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	260 func (v Vertex) Distance2D(w Vertex) float64 {
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	261 return math.Hypot(v.X-w.X, v.Y-w.Y)
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	262 }
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	263
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	264 // Distance returns the distance of the two vertices.
3650 01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3646 diff changeset	265 func (v Vertex) Distance(w Vertex) float64 {
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3646 diff changeset	266 v = v.Sub(w)
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3646 diff changeset	267 return math.Sqrt(v.Dot(v))
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3646 diff changeset	268 }
01ce3ba9b0d0 Fixed generating of random points. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3646 diff changeset	269
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	270 // Minimize adjust this vertex v to hold the minimum
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	271 // values component-wise of v and w.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	272 func (v *Vertex) Minimize(w Vertex) {
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	273 if w.X < v.X {
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	274 v.X = w.X
674 120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	275 }
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	276 if w.Y < v.Y {
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	277 v.Y = w.Y
674 120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	278 }
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	279 if w.Z < v.Z {
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	280 v.Z = w.Z
674 120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	281 }
120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	282 }
120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	283
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	284 // Maximize adjust this vertex v to hold the maximum
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	285 // values component-wise of v and w.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	286 func (v *Vertex) Maximize(w Vertex) {
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	287 if w.X > v.X {
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	288 v.X = w.X
674 120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	289 }
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	290 if w.Y > v.Y {
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	291 v.Y = w.Y
674 120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	292 }
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	293 if w.Z > v.Z {
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	294 v.Z = w.Z
674 120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	295 }
120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	296 }
120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	297
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	298 // SquaredDistance2D returns the squared distances of the
b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	299 // two given vertices in the X/Y plane.
2483 620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	300 func (v Vertex) SquaredDistance2D(w Vertex) float64 {
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	301 dx := v.X - w.X
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	302 dy := v.Y - w.Y
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	303 return dxdx + dydy
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	304 }
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	305
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	306 // Sub2D returns (v - w) component-wise.
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	307 func (v Vertex) Sub2D(w Vertex) Vertex {
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	308 return Vertex{
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	309 X: v.X - w.X,
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	310 Y: v.Y - w.Y,
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	311 }
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	312 }
620038ade708 Incorporated fogleman's fast Delaunay triangulation adjuted to our vertex model. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2472 diff changeset	313
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	314 // Sub returns (v - w) component-wise.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	315 func (v Vertex) Sub(w Vertex) Vertex {
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	316 return Vertex{
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	317 v.X - w.X,
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	318 v.Y - w.Y,
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	319 v.Z - w.Z,
674 120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	320 }
120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	321 }
120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	322
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	323 // Add returns (v + w) component-wise.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	324 func (v Vertex) Add(w Vertex) Vertex {
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	325 return Vertex{
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	326 v.X + w.X,
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	327 v.Y + w.Y,
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	328 v.Z + w.Z,
687 be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	329 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	330 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	331
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	332 // Scale returns s*v component-wise.
730 4c05bdbf8e4b Renamed scale in Vertex to uppercase to make it public. Sascha L. Teichmann <teichmann@intevation.de> parents: 729 diff changeset	333 func (v Vertex) Scale(s float64) Vertex {
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	334 return Vertex{
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	335 s * v.X,
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	336 s * v.Y,
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	337 s * v.Z,
687 be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	338 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	339 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	340
3953 4233570de212 Building octrees: Don't cause havok by building ain overly large tree if the points are not really 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3770 diff changeset	341 // Interpolate returns a function that return s*b[1] + b[0]
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	342 // component-wise.
3953 4233570de212 Building octrees: Don't cause havok by building ain overly large tree if the points are not really 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3770 diff changeset	343 func (b Box) Interpolate() func(Vertex) Vertex {
4233570de212 Building octrees: Don't cause havok by building ain overly large tree if the points are not really 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3770 diff changeset	344 v1, v2 := b[0], b[1]
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	345 v2 = v2.Sub(v1)
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	346 return func(s Vertex) Vertex {
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	347 return Vertex{
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	348 v2.X*s.X + v1.X,
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	349 v2.Y*s.Y + v1.Y,
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	350 v2.Z*s.Z + v1.Z,
674 120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	351 }
120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	352 }
120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	353 }
120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	354
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	355 // Less returns if one of v component is less than the
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	356 // corresponing component in w.
1329 ea2143adc6d3 Named method recievers consistently to make golint happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	357 func (v Vertex) Less(w Vertex) bool {
ea2143adc6d3 Named method recievers consistently to make golint happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	358 return v.X < w.X \|\| v.Y < w.Y \|\| v.Z < w.Z
674 120a82bd9953 octree: Added loading frame for contour generation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	359 }
687 be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	360
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	361 // NewLine return a line of point/direction.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	362 func NewLine(p1, p2 Vertex) Line {
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	363 return Line{
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	364 p2.Sub(p1),
687 be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	365 p1,
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	366 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	367 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	368
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	369 // Eval returns the vertex for t*l[0] + l[1].
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	370 func (l Line) Eval(t float64) Vertex {
730 4c05bdbf8e4b Renamed scale in Vertex to uppercase to make it public. Sascha L. Teichmann <teichmann@intevation.de> parents: 729 diff changeset	371 return l[0].Scale(t).Add(l[1])
687 be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	372 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	373
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	374 // IntersectHorizontal returns the intersection point
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	375 // for a given z value.
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	376 func (l Line) IntersectHorizontal(h float64) Vertex {
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	377 t := (h - l[1].Z) / l[0].Z
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	378 return l.Eval(t)
687 be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	379 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	380
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	381 func side(z, h float64) int {
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	382 switch {
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	383 case z < h:
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	384 return -1
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	385 case z > h:
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	386 return +1
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	387 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	388 return 0
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	389 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	390
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	391 // Dot2 calculates the 2D dot product of the vertices.
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	392 func (v Vertex) Dot2(w Vertex) float64 {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	393 return v.Xw.X + v.Yw.Y
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	394 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	395
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	396 // Contains returns true if the given point is inside the triangle.
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	397 func (t *Triangle) Contains(x, y float64) bool {
4554 23236657b032 Use Sub2D in triangles 2D contains. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4550 diff changeset	398 v0 := t[2].Sub2D(t[0])
23236657b032 Use Sub2D in triangles 2D contains. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4550 diff changeset	399 v1 := t[1].Sub2D(t[0])
23236657b032 Use Sub2D in triangles 2D contains. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4550 diff changeset	400 v2 := Vertex{X: x, Y: y}.Sub2D(t[0])
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	401
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	402 dot00 := v0.Dot2(v0)
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	403 dot01 := v0.Dot2(v1)
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	404 dot02 := v0.Dot2(v2)
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	405 dot11 := v1.Dot2(v1)
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	406 dot12 := v1.Dot2(v2)
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	407
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	408 // Compute barycentric coordinates
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	409 invDenom := 1 / (dot00dot11 - dot01dot01)
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	410 u := (dot11dot02 - dot01dot12) * invDenom
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	411 v := (dot00dot12 - dot01dot02) * invDenom
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	412
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	413 // Check if point is in triangle
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	414 return u >= 0 && v >= 0 && u+v < 1
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	415 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	416
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	417 // IntersectHorizontal calculates the line string that
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	418 // results when cutting a triangle a a certain height.
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	419 // Can be empty (on intersection),
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	420 // one vertex (only touching an vertex) or
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	421 // two vertices (real intersection).
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	422 func (t *Triangle) IntersectHorizontal(h float64) LineStringZ {
687 be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	423 sides := [3]int{
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	424 side(t[0].Z, h),
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	425 side(t[1].Z, h),
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	426 side(t[2].Z, h),
687 be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	427 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	428
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	429 var points LineStringZ
689 614135d69823 octree: prepare storing the lines to file in contouring tool. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 687 diff changeset	430
687 be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	431 for i := 0; i < 3; i++ {
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	432 j := (i + 1) % 3
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	433 si := sides[i]
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	434 sj := sides[j]
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	435
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	436 switch {
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	437 case si == 0 && sj == 0:
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	438 // both on plane
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	439 points = append(points, t[i], t[j])
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	440 case si == 0 && sj != 0:
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	441 // first on plane
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	442 points = append(points, t[i])
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	443 case si != 0 && sj == 0:
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	444 // second on plane
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	445 points = append(points, t[j])
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	446 case si == sj:
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	447 // both on same side
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	448 default:
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	449 // real intersection
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	450 v := NewLine(t[i], t[j]).IntersectHorizontal(h)
687 be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	451 points = append(points, v)
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	452 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	453 }
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	454
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	455 return points
be90ab542aa7 octree: contouring: Do the math to calculate the intersection points of the triangles and the planes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 674 diff changeset	456 }
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	457
768 3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	458 func linearScale(x1, y1, x2, y2 float64) func(Vertex) float64 {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	459 dx := x2 - x1
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	460 dy := y2 - y1
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	461
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	462 switch {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	463 case dx != 0:
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	464 return func(v Vertex) float64 {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	465 return (v.X - x1) / dx
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	466 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	467 case dy != 0:
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	468 return func(v Vertex) float64 {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	469 return (v.Y - y1) / dy
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	470 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	471 default:
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	472 return func(Vertex) float64 {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	473 return 0
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	474 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	475 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	476 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	477
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	478 // BBox calcultes the 2D bounding box in the X/Y plane
b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	479 // of the given line string.
3639 2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	480 func (ls LineStringZ) BBox() Box2D {
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	481
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	482 min := Vertex{math.MaxFloat64, math.MaxFloat64, math.MaxFloat64}
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	483 max := Vertex{-math.MaxFloat64, -math.MaxFloat64, -math.MaxFloat64}
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	484
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	485 for _, v := range ls {
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	486 min.Minimize(v)
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	487 max.Maximize(v)
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	488 }
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	489
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	490 return Box2D{
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	491 X1: min.X,
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	492 Y1: min.Y,
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	493 X2: max.X,
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	494 Y2: max.Y,
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	495 }
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	496 }
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	497
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	498 // Area calculated the area of the line string.
3639 2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	499 func (ls LineStringZ) Area() float64 {
3646 810b28f59b8b Generate random points for second mesh. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3639 diff changeset	500 return polygonArea(ls)
3639 2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	501 }
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	502
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	503 // Reverse reverses the the vertices of this line string in place.
3639 2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	504 func (ls LineStringZ) Reverse() {
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	505 for i, j := 0, len(ls)-1; i < j; i, j = i+1, j-1 {
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	506 ls[i], ls[j] = ls[j], ls[i]
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	507 }
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	508 }
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	509
768 3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	510 func (ls LineStringZ) order(position func(Vertex) float64) {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	511 type posVertex struct {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	512 pos float64
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	513 v Vertex
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	514 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	515 positions := make([]posVertex, len(ls))
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	516 for i, v := range ls {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	517 positions[i] = posVertex{position(v), v}
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	518 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	519 sort.Slice(positions, func(i, j int) bool {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	520 return positions[i].pos < positions[j].pos
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	521 })
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	522 for i := range positions {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	523 ls[i] = positions[i].v
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	524 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	525 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	526
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	527 // EpsEquals returns true if v and w are equal component-wise
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	528 // with the values within a epsilon range.
768 3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	529 func (v Vertex) EpsEquals(w Vertex) bool {
769 ba2007b746ef Fixed nasty typo. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 768 diff changeset	530 const eps = 1e-5
768 3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	531 return math.Abs(v.X-w.X) < eps &&
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	532 math.Abs(v.Y-w.Y) < eps && math.Abs(v.Z-w.Z) < eps
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	533 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	534
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	535 // EpsEquals2D returns true if v and w are equal component-wise
4550 aa2d0006e742 Write iso lines between classes to SVG for debugging. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4542 diff changeset	536 // in the X/Y plane with the values within a epsilon range.
aa2d0006e742 Write iso lines between classes to SVG for debugging. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4542 diff changeset	537 func (v Vertex) EpsEquals2D(w Vertex) bool {
aa2d0006e742 Write iso lines between classes to SVG for debugging. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4542 diff changeset	538 const eps = 1e-5
aa2d0006e742 Write iso lines between classes to SVG for debugging. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4542 diff changeset	539 return math.Abs(v.X-w.X) < eps &&
aa2d0006e742 Write iso lines between classes to SVG for debugging. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4542 diff changeset	540 math.Abs(v.Y-w.Y) < eps
aa2d0006e742 Write iso lines between classes to SVG for debugging. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4542 diff changeset	541 }
aa2d0006e742 Write iso lines between classes to SVG for debugging. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4542 diff changeset	542
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	543 // JoinOnLine joins the the elements of a given multi line string
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	544 // under the assumption that the segments are all on the line segment
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	545 // from (x1, y1) to (x2, y2).
763 d05bc3e34338 More infrastructure for the octree driven cross section endpoint. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 730 diff changeset	546 func (mls MultiLineStringZ) JoinOnLine(x1, y1, x2, y2 float64) MultiLineStringZ {
768 3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	547
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	548 position := linearScale(x1, y1, x2, y2)
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	549
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	550 type posLineString struct {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	551 pos float64
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	552 line LineStringZ
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	553 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	554
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	555 positions := make([]posLineString, 0, len(mls))
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	556
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	557 for _, ls := range mls {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	558 if len(ls) == 0 {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	559 continue
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	560 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	561 // order the atoms first
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	562 ls.order(position)
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	563 positions = append(positions, posLineString{position(ls[0]), ls})
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	564 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	565
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	566 sort.Slice(positions, func(i, j int) bool {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	567 return positions[i].pos < positions[j].pos
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	568 })
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	569
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	570 out := make(MultiLineStringZ, 0, len(positions))
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	571
796 0cc97135717c More fixes to vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 795 diff changeset	572 var ignored int
0cc97135717c More fixes to vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 795 diff changeset	573
768 3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	574 for i := range positions {
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	575 curr := positions[i].line
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	576 if l := len(out); l > 0 {
796 0cc97135717c More fixes to vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 795 diff changeset	577 last := out[l-1]
0cc97135717c More fixes to vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 795 diff changeset	578
0cc97135717c More fixes to vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 795 diff changeset	579 if last[len(last)-1].EpsEquals(curr[0]) {
768 3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	580 out[l-1] = append(last[:len(last)-1], curr...)
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	581 continue
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	582 }
796 0cc97135717c More fixes to vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 795 diff changeset	583 if position(last[len(last)-1]) > position(curr[0]) {
0cc97135717c More fixes to vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 795 diff changeset	584 ignored++
0cc97135717c More fixes to vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 795 diff changeset	585 continue
0cc97135717c More fixes to vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 795 diff changeset	586 }
768 3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	587 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	588 out = append(out, curr)
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	589 }
3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	590
5490 5f47eeea988d Use own logging package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5428 diff changeset	591 log.Infof("ignored parts: %d\n", ignored)
796 0cc97135717c More fixes to vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 795 diff changeset	592
768 3219e7e34953 Join neighbored linestring segments in octree cross sections and sort them along the input line. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 763 diff changeset	593 return out
763 d05bc3e34338 More infrastructure for the octree driven cross section endpoint. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 730 diff changeset	594 }
d05bc3e34338 More infrastructure for the octree driven cross section endpoint. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 730 diff changeset	595
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	596 // Write writes a Vertex as three 64 bit values in little endian order
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	597 // to the given writer.
729 b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	598 func (v *Vertex) Write(w io.Writer) error {
b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	599 if err := binary.Write(
b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	600 w, binary.LittleEndian, math.Float64bits(v.X)); err != nil {
b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	601 return err
b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	602 }
b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	603 if err := binary.Write(
b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	604 w, binary.LittleEndian, math.Float64bits(v.Y)); err != nil {
b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	605 return err
b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	606 }
b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	607 return binary.Write(
b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	608 w, binary.LittleEndian, math.Float64bits(v.Z))
b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	609 }
b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	610
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	611 // Read fills this vertex with three 64 bit values stored as
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	612 // little endian from the given reader.
729 b0bd242ff821 Removed vertex duplicate. Sascha L. Teichmann <teichmann@intevation.de> parents: 727 diff changeset	613 func (v *Vertex) Read(r io.Reader) error {
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	614 var buf [8]byte
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	615 b := buf[:]
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	616 if _, err := io.ReadFull(r, b); err != nil {
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	617 return nil
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	618 }
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	619 v.X = math.Float64frombits(binary.LittleEndian.Uint64(b))
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	620 if _, err := io.ReadFull(r, b); err != nil {
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	621 return nil
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	622 }
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	623 v.Y = math.Float64frombits(binary.LittleEndian.Uint64(b))
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	624 if _, err := io.ReadFull(r, b); err != nil {
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	625 return nil
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	626 }
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	627 v.Z = math.Float64frombits(binary.LittleEndian.Uint64(b))
726 5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	628 return nil
5af9ab39e715 Renamed a few types to uppercase names to prepare the move to the octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 718 diff changeset	629 }
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	630
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	631 // AsWKB returns the WKB representation of the given multi line string.
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	632 func (mls MultiLineStringZ) AsWKB() []byte {
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	633
924 c8146132059e Precalculate WKB size for MultiLineStringZ to avoid reallocations. Sascha L. Teichmann <teichmann@intevation.de> parents: 824 diff changeset	634 // pre-calculate size to avoid reallocations.
c8146132059e Precalculate WKB size for MultiLineStringZ to avoid reallocations. Sascha L. Teichmann <teichmann@intevation.de> parents: 824 diff changeset	635 size := 1 + 4 + 4
c8146132059e Precalculate WKB size for MultiLineStringZ to avoid reallocations. Sascha L. Teichmann <teichmann@intevation.de> parents: 824 diff changeset	636 for _, ml := range mls {
c8146132059e Precalculate WKB size for MultiLineStringZ to avoid reallocations. Sascha L. Teichmann <teichmann@intevation.de> parents: 824 diff changeset	637 size += 1 + 4 + 4 + len(ml)38
c8146132059e Precalculate WKB size for MultiLineStringZ to avoid reallocations. Sascha L. Teichmann <teichmann@intevation.de> parents: 824 diff changeset	638 }
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	639
924 c8146132059e Precalculate WKB size for MultiLineStringZ to avoid reallocations. Sascha L. Teichmann <teichmann@intevation.de> parents: 824 diff changeset	640 buf := bytes.NewBuffer(make([]byte, 0, size))
c8146132059e Precalculate WKB size for MultiLineStringZ to avoid reallocations. Sascha L. Teichmann <teichmann@intevation.de> parents: 824 diff changeset	641
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	642 binary.Write(buf, binary.LittleEndian, wkb.NDR)
63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	643 binary.Write(buf, binary.LittleEndian, wkb.MultiLineStringZ)
924 c8146132059e Precalculate WKB size for MultiLineStringZ to avoid reallocations. Sascha L. Teichmann <teichmann@intevation.de> parents: 824 diff changeset	644 binary.Write(buf, binary.LittleEndian, uint32(len(mls)))
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	645
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	646 for _, ml := range mls {
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	647 binary.Write(buf, binary.LittleEndian, wkb.NDR)
63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	648 binary.Write(buf, binary.LittleEndian, wkb.LineStringZ)
924 c8146132059e Precalculate WKB size for MultiLineStringZ to avoid reallocations. Sascha L. Teichmann <teichmann@intevation.de> parents: 824 diff changeset	649 binary.Write(buf, binary.LittleEndian, uint32(len(ml)))
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	650 for _, p := range ml {
924 c8146132059e Precalculate WKB size for MultiLineStringZ to avoid reallocations. Sascha L. Teichmann <teichmann@intevation.de> parents: 824 diff changeset	651 binary.Write(buf, binary.LittleEndian, math.Float64bits(p.X))
c8146132059e Precalculate WKB size for MultiLineStringZ to avoid reallocations. Sascha L. Teichmann <teichmann@intevation.de> parents: 824 diff changeset	652 binary.Write(buf, binary.LittleEndian, math.Float64bits(p.Y))
c8146132059e Precalculate WKB size for MultiLineStringZ to avoid reallocations. Sascha L. Teichmann <teichmann@intevation.de> parents: 824 diff changeset	653 binary.Write(buf, binary.LittleEndian, math.Float64bits(p.Z))
727 41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	654 }
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	655 }
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	656
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	657 return buf.Bytes()
41c8dc61f38f Moved octree loading stuff to octree package. Sascha L. Teichmann <teichmann@intevation.de> parents: 726 diff changeset	658 }
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	659
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	660 // AsWKB2D returns the WKB representation of the given multi line string
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	661 // leaving the z component out.
925 15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	662 func (mls MultiLineStringZ) AsWKB2D() []byte {
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	663
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	664 // pre-calculate size to avoid reallocations.
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	665 size := 1 + 4 + 4
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	666 for _, ml := range mls {
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	667 size += 1 + 4 + 4 + len(ml)28
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	668 }
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	669
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	670 buf := bytes.NewBuffer(make([]byte, 0, size))
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	671
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	672 binary.Write(buf, binary.LittleEndian, wkb.NDR)
63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	673 binary.Write(buf, binary.LittleEndian, wkb.MultiLineString)
925 15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	674 binary.Write(buf, binary.LittleEndian, uint32(len(mls)))
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	675
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	676 for _, ml := range mls {
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	677 binary.Write(buf, binary.LittleEndian, wkb.NDR)
63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	678 binary.Write(buf, binary.LittleEndian, wkb.LineString)
925 15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	679 binary.Write(buf, binary.LittleEndian, uint32(len(ml)))
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	680 for _, p := range ml {
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	681 binary.Write(buf, binary.LittleEndian, math.Float64bits(p.X))
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	682 binary.Write(buf, binary.LittleEndian, math.Float64bits(p.Y))
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	683 }
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	684 }
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	685
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	686 return buf.Bytes()
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	687 }
15bf101e1522 Send 2D X/Y vertices to the database directly instead of sending 3D data and dropping the Z value afterwards. Sascha L. Teichmann <teichmann@intevation.de> parents: 924 diff changeset	688
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	689 // CCW returns true if this line string is oriented counter clockwise.
4542 56f4e8cbfab7 Sew together triangles that totally inside a single class. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3953 diff changeset	690 func (ls LineStringZ) CCW() bool {
56f4e8cbfab7 Sew together triangles that totally inside a single class. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3953 diff changeset	691 var sum float64
56f4e8cbfab7 Sew together triangles that totally inside a single class. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3953 diff changeset	692 for i, v1 := range ls {
56f4e8cbfab7 Sew together triangles that totally inside a single class. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3953 diff changeset	693 v2 := ls[(i+1)%len(ls)]
56f4e8cbfab7 Sew together triangles that totally inside a single class. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3953 diff changeset	694 sum += (v2.X - v1.X) * (v2.Y + v1.Y)
56f4e8cbfab7 Sew together triangles that totally inside a single class. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3953 diff changeset	695 }
56f4e8cbfab7 Sew together triangles that totally inside a single class. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3953 diff changeset	696 return sum > 0
56f4e8cbfab7 Sew together triangles that totally inside a single class. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3953 diff changeset	697 }
56f4e8cbfab7 Sew together triangles that totally inside a single class. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3953 diff changeset	698
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	699 // Join joins two lines leaving the first of the second out.
928 0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	700 func (ls LineStringZ) Join(other LineStringZ) LineStringZ {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	701 nline := make(LineStringZ, len(ls)+len(other)-1)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	702 copy(nline, ls)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	703 copy(nline[len(ls):], other[1:])
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	704 return nline
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	705 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	706
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	707 // Merge merges line segments of a given multi line string
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	708 // by finding common start and end vertices.
928 0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	709 func (mls MultiLineStringZ) Merge() MultiLineStringZ {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	710
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	711 var out MultiLineStringZ
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	712
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	713 min := Vertex{math.MaxFloat64, math.MaxFloat64, math.MaxFloat64}
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	714
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	715 for _, line := range mls {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	716 for _, v := range line {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	717 min.Minimize(v)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	718 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	719 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	720
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	721 type point struct {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	722 x int64
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	723 y int64
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	724 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	725
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	726 const precision = 1e7
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	727
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	728 quant := func(v Vertex) point {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	729 return point{
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	730 x: int64(math.Round((v.X - min.X) * precision)),
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	731 y: int64(math.Round((v.Y - min.Y) * precision)),
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	732 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	733 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	734
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	735 heads := make(map[point]*[]LineStringZ)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	736
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	737 for _, line := range mls {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	738 if len(line) < 2 {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	739 out = append(out, line)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	740 continue
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	741 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	742 head := quant(line[0])
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	743 tail := quant(line[len(line)-1])
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	744 if head == tail { // its already a ring
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	745 out = append(out, line)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	746 continue
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	747 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	748
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	749 if hs := heads[tail]; hs != nil {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	750 l := len(*hs)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	751 last := (*hs)[l-1]
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	752 if l == 1 {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	753 delete(heads, tail)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	754 } else {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	755 (*hs)[l-1] = nil
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	756 hs = (hs)[:l-1]
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	757 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	758 line = line.Join(last)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	759
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	760 if head == quant(line[len(line)-1]) { // its a ring
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	761 out = append(out, line)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	762 continue
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	763 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	764 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	765
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	766 if hs := heads[head]; hs != nil {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	767 hs = append(hs, line)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	768 } else {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	769 heads[head] = &[]LineStringZ{line}
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	770 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	771 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	772
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	773 again:
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	774 for head, lines := range heads {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	775 for i, line := range *lines {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	776 tail := quant(line[len(line)-1])
940 be7b83638ec8 Simplified line merging in generating contour lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 928 diff changeset	777 for hs := heads[tail]; hs != nil && len(*hs) > 0; hs = heads[tail] {
928 0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	778 l := len(*hs)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	779 last := (*hs)[l-1]
940 be7b83638ec8 Simplified line merging in generating contour lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 928 diff changeset	780 (*hs)[l-1] = nil
be7b83638ec8 Simplified line merging in generating contour lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 928 diff changeset	781 hs = (hs)[:l-1]
928 0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	782 line = line.Join(last)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	783
940 be7b83638ec8 Simplified line merging in generating contour lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 928 diff changeset	784 if tail = quant(line[len(line)-1]); head == tail { // its a ring
928 0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	785 out = append(out, line)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	786 // remove from current lines
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	787 copy((lines)[i:], (lines)[i+1:])
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	788 (lines)[len(lines)-1] = nil
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	789 lines = (lines)[:len(*lines)-1]
940 be7b83638ec8 Simplified line merging in generating contour lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 928 diff changeset	790 goto again
928 0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	791 }
940 be7b83638ec8 Simplified line merging in generating contour lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 928 diff changeset	792 // overwrite in current lines
be7b83638ec8 Simplified line merging in generating contour lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 928 diff changeset	793 (*lines)[i] = line
928 0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	794 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	795 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	796 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	797
988 7dfd3db94e6d In preparation of persisting import jobs logging is done through an interface. Sascha L. Teichmann <teichmann@intevation.de> parents: 968 diff changeset	798 // rings := len(out)
928 0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	799
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	800 // The rest are open line strings.
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	801 for _, lines := range heads {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	802 for _, line := range *lines {
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	803 out = append(out, line)
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	804 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	805 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	806
5490 5f47eeea988d Use own logging package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5428 diff changeset	807 // log.Debugf("segments before/after merge: %d/%d (%d rings)\n",
988 7dfd3db94e6d In preparation of persisting import jobs logging is done through an interface. Sascha L. Teichmann <teichmann@intevation.de> parents: 968 diff changeset	808 // len(mls), len(out), rings)
928 0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	809
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	810 return out
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	811 }
0a0013fbda4a octree2contour: Merge multi linestrings in user land before sending them to database. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 925 diff changeset	812
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	813 // Rect returns the bounding box of this box as separated coordinates.
5320 866eae1bd888 Upgraded the dependencies of the Gemma server. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4854 diff changeset	814 func (a Box2D) Rect() ([2]float64, [2]float64) {
866eae1bd888 Upgraded the dependencies of the Gemma server. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4854 diff changeset	815 return [2]float64{a.X1, a.Y1}, [2]float64{a.X2, a.Y2}
2492 10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2483 diff changeset	816 }
10681749371d Implemented the BBox/clipping polygon test. TODO: triangle/clipping polygon test. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2483 diff changeset	817
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	818 // Intersects checks if two Box2Ds intersect.
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	819 func (a Box2D) Intersects(b Box2D) bool {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	820 return !(a.X2 < a.X1 \|\| a.X2 < b.X1 \|\|
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	821 a.Y2 < a.Y1 \|\| a.Y2 < b.Y1)
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	822 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	823
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	824 // Contains returns true if the given point is inside the box.
2466 a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	825 func (a Box2D) Contains(x, y float64) bool {
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	826 return a.X1 <= x && x <= a.X2 &&
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	827 a.Y1 <= y && y <= a.Y2
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	828 }
a1e751c08c56 Calculate difference on single core. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1692 diff changeset	829
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	830 // Xi returns the i-th x component.
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	831 func (a Box2D) Xi(i int) float64 {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	832 if i == 0 {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	833 return a.X1
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	834 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	835 return a.X2
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	836 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	837
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	838 // Yi returns the i-th y component.
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	839 func (a Box2D) Yi(i int) float64 {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	840 if i == 0 {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	841 return a.Y1
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	842 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	843 return a.Y2
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	844 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	845
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	846 // Union calculates the united bounding box of the two given boxes.
2512 2768f74d54ab Added an STRTree for the triangulation. Should be better suited for filtering out the clipped triangles. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	847 func (a Box2D) Union(b Box2D) Box2D {
2768f74d54ab Added an STRTree for the triangulation. Should be better suited for filtering out the clipped triangles. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	848 return Box2D{
2768f74d54ab Added an STRTree for the triangulation. Should be better suited for filtering out the clipped triangles. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	849 X1: math.Min(a.X1, b.X1),
2768f74d54ab Added an STRTree for the triangulation. Should be better suited for filtering out the clipped triangles. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	850 Y1: math.Min(a.Y1, b.Y1),
2768f74d54ab Added an STRTree for the triangulation. Should be better suited for filtering out the clipped triangles. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	851 X2: math.Max(a.X2, b.X2),
2768f74d54ab Added an STRTree for the triangulation. Should be better suited for filtering out the clipped triangles. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	852 Y2: math.Max(a.Y2, b.Y2),
2768f74d54ab Added an STRTree for the triangulation. Should be better suited for filtering out the clipped triangles. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	853 }
2768f74d54ab Added an STRTree for the triangulation. Should be better suited for filtering out the clipped triangles. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	854 }
2768f74d54ab Added an STRTree for the triangulation. Should be better suited for filtering out the clipped triangles. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2494 diff changeset	855
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	856 // Area returns the area of the box.
3639 2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	857 func (a Box2D) Area() float64 {
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	858 return (a.X2 - a.X1) * (a.Y2 - a.Y1)
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	859 }
2a4216c81e7b Extract the removed triangles from first triangulation, too. Useful to build a artifical DEM for second pass. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2572 diff changeset	860
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	861 // NewPlane2D creates a new Plane2D from two given points.
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	862 func NewPlane2D(x1, y1, x2, y2 float64) Plane2D {
778 9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	863 b := x2 - x1
9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	864 a := -(y2 - y1)
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	865
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	866 l := math.Sqrt(aa + bb)
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	867 a /= l
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	868 b /= l
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	869
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	870 // ax1 + by1 + c = 0
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	871 c := -(ax1 + by1)
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	872 return Plane2D{a, b, c}
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	873 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	874
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	875 // Eval determines the distance of a given point
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	876 // from the plane. The sign of the result indicates
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	877 // the sideness.
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	878 func (p Plane2D) Eval(x, y float64) float64 {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	879 return p.Ax + p.By + p.C
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	880 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	881
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	882 const epsPlane = 1e-5
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	883
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	884 func sides(s int, x float64) int {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	885 if math.Signbit(x) {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	886 return s \| 2
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	887 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	888 return s \| 1
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	889 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	890
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	891 // IntersectsPlane checks if a Box2D intersects with
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	892 // a given Plane2D.
774 c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	893 func (a Box2D) IntersectsPlane(p Plane2D) bool {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	894 var s int
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	895 for i := 0; i < 2; i++ {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	896 x := a.Xi(i)
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	897 for j := 0; j < 2; j++ {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	898 y := a.Yi(j)
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	899 v := p.Eval(x, y)
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	900 if math.Abs(v) < epsPlane {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	901 //log.Printf("on line")
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	902 return true
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	903 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	904 if s = sides(s, v); s == 3 {
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	905 //log.Printf("... on both sides (djt)")
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	906 return true
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	907 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	908 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	909 }
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	910 //log.Printf("side: %d\n", s)
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	911 return false
c55771b7c502 Moved Box2D and Plane2D into vertex.go and made to API public. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 769 diff changeset	912 }
778 9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	913
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	914 // Cross calculates the cross product of two vertices.
1329 ea2143adc6d3 Named method recievers consistently to make golint happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	915 func (v Vertex) Cross(w Vertex) Vertex {
778 9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	916 return Vertex{
1329 ea2143adc6d3 Named method recievers consistently to make golint happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	917 v.Yw.Z - v.Zw.Y,
ea2143adc6d3 Named method recievers consistently to make golint happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	918 v.Zw.X - v.Xw.Z,
ea2143adc6d3 Named method recievers consistently to make golint happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	919 v.Xw.Y - v.Yw.X,
778 9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	920 }
9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	921 }
9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	922
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	923 // Intersection calcultes the 2D intersection point of
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	924 // two Plane2Ds. If they do not intersect the returned
f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	925 // bool flags is set to false.
1329 ea2143adc6d3 Named method recievers consistently to make golint happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	926 func (p Plane2D) Intersection(o Plane2D) (float64, float64, bool) {
778 9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	927
1329 ea2143adc6d3 Named method recievers consistently to make golint happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	928 u1 := Vertex{p.A, p.B, p.C}
ea2143adc6d3 Named method recievers consistently to make golint happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	929 u2 := Vertex{o.A, o.B, o.C}
778 9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	930
1329 ea2143adc6d3 Named method recievers consistently to make golint happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	931 plane := u1.Cross(u2)
778 9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	932
1329 ea2143adc6d3 Named method recievers consistently to make golint happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	933 if plane.Z == 0 {
778 9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	934 return 0, 0, false
9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	935 }
9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	936
1329 ea2143adc6d3 Named method recievers consistently to make golint happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1017 diff changeset	937 return plane.X / plane.Z, plane.Y / plane.Z, true
778 9be20bd0f131 Fixed a bug with 2d planes (lines) found while working on line intersections. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 774 diff changeset	938 }
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	939
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	940 // VerticalLine is a 2D line segment.
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	941 type VerticalLine struct {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	942 x1 float64
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	943 y1 float64
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	944 x2 float64
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	945 y2 float64
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	946
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	947 line Plane2D
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	948 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	949
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	950 // NewVerticalLine creates a new 2D line segment.
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	951 func NewVerticalLine(x1, y1, x2, y2 float64) *VerticalLine {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	952 return &VerticalLine{
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	953 x1: x1,
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	954 y1: y1,
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	955 x2: x2,
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	956 y2: y2,
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	957 line: NewPlane2D(x1, y1, x2, y2),
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	958 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	959 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	960
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	961 func onPlane(x float64) bool { return math.Abs(x) < epsPlane }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	962
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	963 func relative(v1, v2 Vertex) func(x, y float64) float64 {
795 bd5f38eb6153 Vertical Triangle intersection: Fixed function to determine where intersection point Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 792 diff changeset	964 ls := linearScale(v1.X, v1.Y, v2.X, v2.Y)
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	965 return func(x, y float64) float64 {
795 bd5f38eb6153 Vertical Triangle intersection: Fixed function to determine where intersection point Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 792 diff changeset	966 return ls(Vertex{x, y, 0})
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	967 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	968 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	969
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	970 func interpolate(a, b float64) func(float64) float64 {
9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	971 return func(x float64) float64 {
9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	972 return (b-a)*x + a
9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	973 }
9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	974 }
9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	975
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	976 func linear(v1, v2 Vertex) func(t float64) Vertex {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	977 return func(t float64) Vertex {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	978 return Vertex{
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	979 (v2.X-v1.X)*t + v1.X,
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	980 (v2.Y-v1.Y)*t + v1.Y,
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	981 (v2.Z-v1.Z)*t + v1.Z,
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	982 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	983 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	984 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	985
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	986 func inRange(a float64) bool { return 0 <= a && a <= 1 }
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	987
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	988 // Intersection intersects a line segment with a triangle.
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	989 func (vl VerticalLine) Intersection(t Triangle) LineStringZ {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	990
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	991 var out LineStringZ
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	992
815 01019d4c8359 Cross section: Made triangle intersection less spammy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 813 diff changeset	993 //defer func() { log.Printf("length out: %d\n", len(out)) }()
796 0cc97135717c More fixes to vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 795 diff changeset	994
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	995 edges:
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	996 for i := 0; i < 3 && len(out) < 2; i++ {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	997 j := (i + 1) % 3
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	998 edge := NewPlane2D(t[i].X, t[i].Y, t[j].X, t[j].Y)
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	999
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1000 s1 := edge.Eval(vl.x1, vl.y1)
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1001 s2 := edge.Eval(vl.x2, vl.y2)
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1002
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1003 o1 := onPlane(s1)
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1004 o2 := onPlane(s2)
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1005
815 01019d4c8359 Cross section: Made triangle intersection less spammy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 813 diff changeset	1006 // log.Printf("s1, s2: %t %t (%f %f)\n", o1, o2, s1, s2)
795 bd5f38eb6153 Vertical Triangle intersection: Fixed function to determine where intersection point Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 792 diff changeset	1007
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1008 switch {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1009 case o1 && o2:
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1010 pos := relative(t[i], t[j])
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1011 t1 := pos(vl.x1, vl.y1)
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1012 t2 := pos(vl.x2, vl.y2)
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1013
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1014 r1 := inRange(t1)
9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1015 r2 := inRange(t2)
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1016
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1017 switch {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1018 case r1 && r2:
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1019 lin := linear(t[i], t[j])
9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1020 out = append(out, lin(t1), lin(t2))
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1021 return out
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1022
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1023 case !r1 && !r2: // the hole edge
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1024 out = append(out, t[i], t[j])
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1025 return out
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1026 case !r1:
789 d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1027 if t1 < 0 {
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1028 // below first -> clip by first
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1029 lin := linear(t[i], t[j])
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1030 out = append(out, t[i], lin(t2))
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1031 } else {
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1032 // above second -> clip by second
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1033 lin := linear(t[i], t[j])
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1034 out = append(out, lin(t2), t[j])
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1035 }
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1036 return out
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1037 case !r2:
789 d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1038 if t2 < 0 {
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1039 // below first -> clip by first
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1040 lin := linear(t[i], t[j])
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1041 out = append(out, t[i], lin(t1))
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1042 } else {
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1043 // above second -> clip by second
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1044 lin := linear(t[i], t[j])
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1045 out = append(out, lin(t1), t[j])
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1046 }
d43e61044ad8 Fill in last remaining corner cases in vertical triangle interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 788 diff changeset	1047 return out
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1048 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1049
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1050 case o1:
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1051 t1 := relative(t[i], t[j])(vl.x1, vl.y1)
9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1052 if !inRange(t1) {
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1053 continue edges
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1054 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1055 out = append(out, linear(t[i], t[j])(t1))
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1056
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1057 case o2:
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1058 t2 := relative(t[i], t[j])(vl.x2, vl.y2)
9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1059 if !inRange(t2) {
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1060 continue edges
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1061 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1062 out = append(out, linear(t[i], t[j])(t2))
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1063
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1064 default:
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1065 x, y, intersects := vl.line.Intersection(edge)
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1066 if !intersects {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1067 continue edges
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1068 }
795 bd5f38eb6153 Vertical Triangle intersection: Fixed function to determine where intersection point Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 792 diff changeset	1069
815 01019d4c8359 Cross section: Made triangle intersection less spammy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 813 diff changeset	1070 // log.Println("Intersection -----------------------------")
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1071 t1 := relative(t[i], t[j])(x, y)
815 01019d4c8359 Cross section: Made triangle intersection less spammy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 813 diff changeset	1072 // log.Printf("relative pos: %f\n", t1)
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1073 if !inRange(t1) {
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1074 continue edges
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1075 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1076
815 01019d4c8359 Cross section: Made triangle intersection less spammy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 813 diff changeset	1077 // log.Println("valid ***************")
795 bd5f38eb6153 Vertical Triangle intersection: Fixed function to determine where intersection point Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 792 diff changeset	1078
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1079 z := interpolate(t[j].Z, t[i].Z)(t1)
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1080 n := Vertex{x, y, z}
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1081
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1082 if math.Signbit(s1) != math.Signbit(s2) {
815 01019d4c8359 Cross section: Made triangle intersection less spammy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 813 diff changeset	1083 // log.Println("\ton different sides")
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1084 // different sides -> vertex on edge.
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1085 out = append(out, n)
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1086 } else { // same side -> inside. Find peer
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1087 if len(out) > 0 { // we already have our peer
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1088 out = append(out, n)
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1089 return out
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1090 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1091
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1092 for k := 0; k < 3; k++ {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1093 if k == i {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1094 continue
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1095 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1096 l := (k + 1) % 3
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1097 other := NewPlane2D(t[k].X, t[k].Y, t[l].X, t[l].Y)
813 1a808929c2c5 Cross sections: Fixed problem that certain inputs don't deliver data. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 796 diff changeset	1098 xo, yo, intersects := vl.line.Intersection(other)
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1099 if !intersects {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1100 continue
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1101 }
813 1a808929c2c5 Cross sections: Fixed problem that certain inputs don't deliver data. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 796 diff changeset	1102 t2 := relative(t[k], t[l])(xo, yo)
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1103 if !inRange(t2) {
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1104 continue
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1105 }
813 1a808929c2c5 Cross sections: Fixed problem that certain inputs don't deliver data. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 796 diff changeset	1106 zo := interpolate(t[k].Z, t[l].Z)(t2)
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1107
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1108 m := Vertex{xo, yo, zo}
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1109
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1110 var xn, yn, xf, yf float64
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1111
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1112 // Which is nearer to current edge?
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1113 if math.Abs(s1) < math.Abs(s2) {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1114 xn, yn = vl.x1, vl.y1
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1115 xf, yf = vl.x2, vl.y2
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1116 } else {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1117 xn, yn = vl.x2, vl.y2
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1118 xf, yf = vl.x1, vl.y1
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1119 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1120
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1121 if onPlane(other.Eval(xn, yn)) {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1122 // triangle intersect in only point
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1123 // treat as no intersection
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1124 break edges
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1125 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1126
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1127 pos := relative(n, m)
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1128
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1129 tzn := pos(xn, yn)
9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1130 tzf := pos(xf, yf)
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1131
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1132 if !inRange(tzn) {
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1133 // if near is out of range far is, too.
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1134 return out
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1135 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1136
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1137 lin := interpolate(n.Z, m.Z)
9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1138
9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1139 if inRange(tzf) {
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1140 m.X = xf
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1141 m.Y = yf
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1142 m.Z = lin(tzf)
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1143 } // else m is clipping
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1144
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1145 n.X = xn
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1146 n.Y = yn
788 9f3a4a60dc04 Cleaned up interpolation mess in vertical triangle intersection. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 787 diff changeset	1147 n.Z = lin(tzn)
787 3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1148
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1149 out = append(out, n, m)
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1150 return out
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1151 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1152 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1153 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1154 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1155
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1156 // supress single point touches.
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1157 if len(out) == 1 {
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1158 out = out[:0]
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1159 }
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1160
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1161 return out
3d927e06b92c Triangle intersection. WIP. Currently interpolation is messed up. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 778 diff changeset	1162 }
959 6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1163
5706 148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1164 // QuantizeXY quantize the X/Y values to scale and
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1165 // removes duplicates in place. The z value of
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1166 // duplicates will be averaged pairwise.
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1167 func (mpz MultiPointZ) QuantizeXY(scale float64) MultiPointZ {
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1168 type qpoint struct {
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1169 x int64
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1170 y int64
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1171 }
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1172 m := make(map[qpoint]float64, len(mpz))
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1173 for _, p := range mpz {
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1174 k := qpoint{
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1175 x: int64(math.Round(p.X * scale)),
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1176 y: int64(math.Round(p.Y * scale)),
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1177 }
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1178 if z, ok := m[k]; ok {
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1179 m[k] = (z + p.Z) * 0.5
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1180 } else {
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1181 m[k] = p.Z
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1182 }
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1183 }
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1184 i := 0
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1185 invScale := 1 / scale
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1186 for k, z := range m {
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1187 mpz[i] = Vertex{
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1188 X: float64(k.x) * invScale,
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1189 Y: float64(k.y) * invScale,
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1190 Z: z,
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1191 }
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1192 i++
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1193 }
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1194 return mpz[:i]
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1195 }
148abae1fcd0 Quantize xyz in points in X/Y to 1/QuantScale. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5703 diff changeset	1196
5428 b8d5f1cd15fb Simplified classification. Needs testing. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5419 diff changeset	1197 // Filter returns a copy removed the vertices which
b8d5f1cd15fb Simplified classification. Needs testing. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5419 diff changeset	1198 // don't pass the filter test.
b8d5f1cd15fb Simplified classification. Needs testing. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5419 diff changeset	1199 func (mpz MultiPointZ) Filter(filter func(Vertex) bool) MultiPointZ {
b8d5f1cd15fb Simplified classification. Needs testing. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5419 diff changeset	1200 n := make(MultiPointZ, 0, len(mpz))
b8d5f1cd15fb Simplified classification. Needs testing. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5419 diff changeset	1201 for _, v := range mpz {
b8d5f1cd15fb Simplified classification. Needs testing. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5419 diff changeset	1202 if filter(v) {
b8d5f1cd15fb Simplified classification. Needs testing. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5419 diff changeset	1203 n = append(n, v)
5418 e89ff1894bb4 Don't clip the points against the border polygon. The implemented way was wrong. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5416 diff changeset	1204 }
e89ff1894bb4 Don't clip the points against the border polygon. The implemented way was wrong. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5416 diff changeset	1205 }
5428 b8d5f1cd15fb Simplified classification. Needs testing. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5419 diff changeset	1206 return n
5418 e89ff1894bb4 Don't clip the points against the border polygon. The implemented way was wrong. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5416 diff changeset	1207 }
e89ff1894bb4 Don't clip the points against the border polygon. The implemented way was wrong. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5416 diff changeset	1208
5601 1222b777f51f Made golint finally happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5490 diff changeset	1209 // MinMax returns the extend of the point set.
5428 b8d5f1cd15fb Simplified classification. Needs testing. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5419 diff changeset	1210 func (mpz MultiPointZ) MinMax() (Vertex, Vertex) {
5413 99f32f0dc70c Added tools to classify vertices by their heights with class breaks. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5320 diff changeset	1211 min := Vertex{math.MaxFloat64, math.MaxFloat64, math.MaxFloat64}
99f32f0dc70c Added tools to classify vertices by their heights with class breaks. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5320 diff changeset	1212 max := Vertex{-math.MaxFloat64, -math.MaxFloat64, -math.MaxFloat64}
5428 b8d5f1cd15fb Simplified classification. Needs testing. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5419 diff changeset	1213 for _, v := range mpz {
5413 99f32f0dc70c Added tools to classify vertices by their heights with class breaks. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5320 diff changeset	1214 min.Minimize(v)
99f32f0dc70c Added tools to classify vertices by their heights with class breaks. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5320 diff changeset	1215 max.Maximize(v)
99f32f0dc70c Added tools to classify vertices by their heights with class breaks. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5320 diff changeset	1216 }
99f32f0dc70c Added tools to classify vertices by their heights with class breaks. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5320 diff changeset	1217 return min, max
99f32f0dc70c Added tools to classify vertices by their heights with class breaks. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5320 diff changeset	1218 }
99f32f0dc70c Added tools to classify vertices by their heights with class breaks. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5320 diff changeset	1219
1692 f4dcbe8941a1 Octree: Resolved the remaing golint issues with this package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1691 diff changeset	1220 // AsWKB returns a WKB representation of the given point cloud.
959 6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1221 func (mpz MultiPointZ) AsWKB() []byte {
6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1222 size := 1 + 4 + 4 + len(mpz)(1+4+38)
6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1223
6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1224 buf := bytes.NewBuffer(make([]byte, 0, size))
6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1225
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	1226 binary.Write(buf, binary.LittleEndian, wkb.NDR)
63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2466 diff changeset	1227 binary.Write(buf, binary.LittleEndian, wkb.MultiPointZ)
959 6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1228 binary.Write(buf, binary.LittleEndian, uint32(len(mpz)))
6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1229
2472 db0e4ab57977 Turn new point cloud into WKB form of MultiPointZ. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2471 diff changeset	1230 perPoint := bytes.NewBuffer(make([]byte, 0, 1+4))
db0e4ab57977 Turn new point cloud into WKB form of MultiPointZ. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2471 diff changeset	1231 binary.Write(perPoint, binary.LittleEndian, wkb.NDR)
db0e4ab57977 Turn new point cloud into WKB form of MultiPointZ. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2471 diff changeset	1232 binary.Write(perPoint, binary.LittleEndian, wkb.PointZ)
db0e4ab57977 Turn new point cloud into WKB form of MultiPointZ. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2471 diff changeset	1233 hdr := perPoint.Bytes()
db0e4ab57977 Turn new point cloud into WKB form of MultiPointZ. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2471 diff changeset	1234
959 6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1235 for _, p := range mpz {
2472 db0e4ab57977 Turn new point cloud into WKB form of MultiPointZ. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2471 diff changeset	1236 buf.Write(hdr)
959 6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1237 binary.Write(buf, binary.LittleEndian, math.Float64bits(p.X))
6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1238 binary.Write(buf, binary.LittleEndian, math.Float64bits(p.Y))
6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1239 binary.Write(buf, binary.LittleEndian, math.Float64bits(p.Z))
6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1240 }
6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1241
6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1242 return buf.Bytes()
6ab012d0f0c2 Started with writing an importer job for sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 940 diff changeset	1243 }
2529 45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1244
4854 b9599e5e8004 Added missing doc strings to the geometry models in the mesh package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4833 diff changeset	1245 // FromWKB de-serializes this multi point z geometry from a WKB representation.
2529 45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1246 func (mpz *MultiPointZ) FromWKB(data []byte) error {
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1247
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1248 r := bytes.NewReader(data)
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1249
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1250 endian, err := r.ReadByte()
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1251
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1252 var order binary.ByteOrder
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1253
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1254 switch {
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1255 case err != nil:
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1256 return err
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1257 case endian == wkb.NDR:
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1258 order = binary.LittleEndian
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1259 case endian == wkb.XDR:
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1260 order = binary.BigEndian
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1261 default:
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1262 return fmt.Errorf("unknown byte order %x", endian)
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1263 }
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1264
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1265 var geomType uint32
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1266 err = binary.Read(r, order, &geomType)
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1267
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1268 switch {
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1269 case err != nil:
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1270 return err
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1271 case geomType != wkb.MultiPointZ:
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1272 return fmt.Errorf("unknown geometry type %x", geomType)
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1273 }
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1274
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1275 var numPoints uint32
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1276 if err := binary.Read(r, order, &numPoints); err != nil {
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1277 return err
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1278 }
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1279
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1280 points := make(MultiPointZ, numPoints)
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1281
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1282 for i := range points {
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1283 endian, err = r.ReadByte()
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1284
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1285 switch {
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1286 case err != nil:
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1287 return err
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1288 case endian == wkb.NDR:
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1289 order = binary.LittleEndian
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1290 case endian == wkb.XDR:
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1291 order = binary.BigEndian
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1292 default:
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1293 return fmt.Errorf("unknown byte order %x", endian)
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1294 }
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1295
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1296 err = binary.Read(r, order, &geomType)
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1297
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1298 switch {
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1299 case err != nil:
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1300 return err
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1301 case geomType != wkb.PointZ:
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1302 return fmt.Errorf("unknown geometry type %x", geomType)
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1303 }
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1304
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1305 var x, y, z uint64
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1306 if err = binary.Read(r, order, &x); err != nil {
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1307 return err
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1308 }
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1309 if err = binary.Read(r, order, &y); err != nil {
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1310 return err
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1311 }
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1312 if err = binary.Read(r, order, &z); err != nil {
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1313 return err
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1314 }
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1315 points[i] = Vertex{
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1316 X: math.Float64frombits(x),
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1317 Y: math.Float64frombits(y),
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1318 Z: math.Float64frombits(z),
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1319 }
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1320 }
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1321
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1322 *mpz = points
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1323
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1324 return nil
45d51a49f191 SR import: Use own triangulation and clipping when importing sounding results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2512 diff changeset	1325 }

Mercurial > gemma

annotate pkg/mesh/vertex.go @ 5710:37c8feeecb4d