gemma: pkg/octree/tin.go annotate

annotate pkg/octree/tin.go @ 1130:42617bba8709

Go from 50cm to 10cm stepwidth in calculating the contorlines.

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Wed, 07 Nov 2018 17:54:31 +0100
parents	a244b18cb916
children	d753ce6cf588

rev	line source
1017 a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 971 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: 971 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: 971 diff changeset	3 //
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 971 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: 971 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: 971 diff changeset	6 //
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 971 diff changeset	7 // Copyright (C) 2018 by via donau
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 971 diff changeset	8 // – Österreichische Wasserstraßen-Gesellschaft mbH
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 971 diff changeset	9 // Software engineering by Intevation GmbH
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 971 diff changeset	10 //
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 971 diff changeset	11 // Author(s):
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 971 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: 971 diff changeset	13
968 a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	14 package octree
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16 import (
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17 "bytes"
971 f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	18 "context"
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	19 "database/sql"
968 a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	20 "encoding/binary"
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	21 "errors"
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	22 "fmt"
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	23 "io"
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	24 "log"
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	25 "math"
971 f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	26 "time"
968 a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	27 )
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29 var (
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30 errNoByteSlice = errors.New("Not a byte slice")
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	31 errTooLessPoints = errors.New("Too less points")
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32 )
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	34 const wgs84 = 4326
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	35
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	36 type Tin struct {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37 EPSG uint32
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38 Vertices []Vertex
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	39 Triangles [][]int32
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	40
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	41 Min Vertex
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	42 Max Vertex
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	43 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	44
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	45 func (t *Tin) FromWKB(data []byte) error {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	46 log.Printf("data length %d\n", len(data))
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	47
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	48 r := bytes.NewReader(data)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	49
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	50 endian, err := r.ReadByte()
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	51
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 var order binary.ByteOrder
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	53
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54 switch {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55 case err != nil:
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57 case endian == wkbNDR:
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	58 order = binary.LittleEndian
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	59 case endian == wkbXDR:
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	60 order = binary.BigEndian
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	61 default:
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	62 return fmt.Errorf("unknown byte order %x", endian)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	63 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	64
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	65 var geomType uint32
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	66 err = binary.Read(r, order, &geomType)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	67
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68 switch {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69 case err != nil:
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	71 case geomType != wkbTinZ:
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	72 return fmt.Errorf("unknown geometry type %x", geomType)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	73 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	74
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	75 var num uint32
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	76 if err = binary.Read(r, order, &num); err != nil {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	77 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	78 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	79
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	80 vertices := make([]Vertex, 0, 100000)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	81
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	82 var v Vertex
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	83
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	84 v2i := make(map[Vertex]int32, 100000)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	85
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	86 var indexPool []int32
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	87
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	88 allocIndices := func() []int32 {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	89 if len(indexPool) == 0 {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	90 indexPool = make([]int32, 381024)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	91 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	92 ids := indexPool[:3]
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	93 indexPool = indexPool[3:]
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	94 return ids
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	95 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	96
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	97 var triangles [][]int32
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	98
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	99 min := Vertex{math.MaxFloat64, math.MaxFloat64, math.MaxFloat64}
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	100 max := Vertex{-math.MaxFloat64, -math.MaxFloat64, -math.MaxFloat64}
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	101
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	102 for i := uint32(0); i < num; i++ {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	103
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	104 endian, err = r.ReadByte()
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	105 switch {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	106 case err != nil:
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	107 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	108 case endian == wkbNDR:
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	109 order = binary.LittleEndian
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	110 case endian == wkbXDR:
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	111 order = binary.BigEndian
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	112 default:
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	113 return fmt.Errorf("unknown byte order %x", endian)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	114 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	115
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	116 err = binary.Read(r, order, &geomType)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	117 switch {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	118 case err != nil:
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	119 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	120 case geomType != wkbTriangleZ:
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	121 return fmt.Errorf("unknown geometry type %d", geomType)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	122 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	123
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	124 var rings uint32
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	125 if err = binary.Read(r, order, &rings); err != nil {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	126 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	127 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	128 triangle := allocIndices()
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	129
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	130 for ring := uint32(0); ring < rings; ring++ {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	131 var npoints uint32
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	132 if err = binary.Read(r, order, &npoints); err != nil {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	133 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	134 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	135
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	136 if npoints < 3 {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	137 return errTooLessPoints
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	138 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	139
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140 for p := uint32(0); p < npoints; p++ {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	141 var x, y, z uint64
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142 for _, addr := range []*uint64{&x, &y, &z} {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	143 if err = binary.Read(r, order, addr); err != nil {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	144 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	145 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	146 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	147 if p >= 3 \|\| ring >= 1 {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148 // Don't store the forth point.
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	149 continue
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	150 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	151 // Do this conversion later to spare reflect calls
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	152 // and allocs in binary.Read.
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	153 v.X = math.Float64frombits(x)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	154 v.Y = math.Float64frombits(y)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	155 v.Z = math.Float64frombits(z)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	156 idx, found := v2i[v]
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	157 if !found {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	158 idx = int32(len(vertices))
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159 v2i[v] = idx
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	160 vertices = append(vertices, v)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	161 min.Minimize(v)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	162 max.Maximize(v)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	163 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	164 triangle[p] = idx
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	165 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	166 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	167 triangles = append(triangles, triangle)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	168 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	170 log.Printf("bbox: [[%f, %f], [%f, %f]]\n",
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	171 min.X, min.Y, max.X, max.Y)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	172
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	173 *t = Tin{
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174 EPSG: wgs84,
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	175 Vertices: vertices,
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176 Triangles: triangles,
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	177 Min: min,
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	178 Max: max,
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	179 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	180
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	181 return nil
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	183
971 f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	184 const (
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	185 tinSQLPrefix = `WITH trans AS (
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	186 SELECT
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	187 ST_Buffer(ST_Transform(area::geometry, $1::int), 0.001) AS area,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	188 ST_Transform(point_cloud::geometry, $1::int) AS point_cloud
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	189 FROM waterway.sounding_results
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	190 `
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	191 tinSQLSuffix = `
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	192 ),
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	193 triangles AS (
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	194 SELECT t.geom AS geom, ST_MakePolygon(ST_ExteriorRing(t.geom)) AS poly FROM (
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	195 SELECT (ST_Dump(
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	196 ST_DelaunayTriangles(point_cloud, 0, 2))).geom
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	197 FROM trans) t
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	198 )
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	199 SELECT ST_AsBinary(ST_Collect(triangles.geom)) FROM triangles, trans
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	200 WHERE ST_Covers(trans.area, triangles.poly)`
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	201
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	202 loadTinSQL = tinSQLPrefix + `WHERE bottleneck_id = $2 AND date_info = $3` + tinSQLSuffix
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	203 loadTinByIDSQL = tinSQLPrefix + `WHERE id = $2` + tinSQLSuffix
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	204 )
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	205
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	206 func GenerateTin(
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	207 conn *sql.Conn,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	208 ctx context.Context,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	209 bottleneck string,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	210 date time.Time,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	211 epsg uint32,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	212 ) (*Tin, error) {
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	213 var tin Tin
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	214 err := conn.QueryRowContext(ctx, loadTinSQL, epsg, bottleneck, date).Scan(&tin)
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	215 switch {
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	216 case err == sql.ErrNoRows:
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	217 return nil, nil
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	218 case err != nil:
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	219 return nil, err
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	220 }
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	221 tin.EPSG = epsg
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	222 return &tin, nil
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	223 }
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	224
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	225 func GenerateTinByID(
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	226 conn *sql.Conn,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	227 ctx context.Context,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	228 id int64,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	229 epsg uint32,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	230 ) (*Tin, error) {
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	231 var tin Tin
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	232 err := conn.QueryRowContext(ctx, loadTinByIDSQL, epsg, id).Scan(&tin)
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	233 switch {
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	234 case err == sql.ErrNoRows:
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	235 return nil, nil
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	236 case err != nil:
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	237 return nil, err
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	238 }
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	239 tin.EPSG = epsg
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	240 return &tin, nil
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	241 }
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	242
968 a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	243 func (t *Tin) Scan(raw interface{}) error {
971 f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	244 if raw == nil {
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	245 return nil
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	246 }
968 a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	247 data, ok := raw.([]byte)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	248 if !ok {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	249 return errNoByteSlice
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	250 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	251 return t.FromWKB(data)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	252 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	253
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	254 func (t *Tin) Serialize(w io.Writer) error {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	255
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	256 if err := binary.Write(w, binary.LittleEndian, t.EPSG); err != nil {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	257 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	258 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	259
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	260 if err := t.Min.Write(w); err != nil {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	261 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	262 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	263 if err := t.Max.Write(w); err != nil {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	264 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	265 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	266
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	267 if err := binary.Write(
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	268 w, binary.LittleEndian, uint32(len(t.Vertices))); err != nil {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	269 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	270 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	271
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	272 for _, v := range t.Vertices {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	273 if err := v.Write(w); err != nil {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	274 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	275 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	276 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	277 log.Printf("vertices %d (%d)\n", len(t.Vertices), len(t.Vertices)38)
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	278
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	279 if err := binary.Write(
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	280 w, binary.LittleEndian, uint32(len(t.Triangles))); err != nil {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	281 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	282 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	283
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	284 var buf [binary.MaxVarintLen32]byte
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	285 var written int
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	286 var last int32
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	287 for _, triangle := range t.Triangles {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	288 for _, idx := range triangle {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	289 value := idx - last
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	290 n := binary.PutVarint(buf[:], int64(value))
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	291 for p := buf[:n]; len(p) > 0; p = p[n:] {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	292 var err error
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	293 if n, err = w.Write(p); err != nil {
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	294 return err
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	295 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	296 written += n
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	297 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	298 last = idx
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	299 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	300 }
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	301 log.Printf("compressed tin indices in bytes: %d (%d)\n",
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	302 written, 34len(t.Triangles))
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	303
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	304 return nil
a4fe07a21ba7 Moved octree builder into octree package to be reusable by the sounding result import job. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	305 }

Mercurial > gemma

annotate pkg/octree/tin.go @ 1130:42617bba8709