gemma: pkg/octree/tin.go annotate

annotate pkg/octree/tin.go @ 971:f9fb6c399f3f

Moved generating of tins to octree package.

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Thu, 18 Oct 2018 11:52:13 +0200
parents	a4fe07a21ba7
children	a244b18cb916

rev	line source
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	1 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	2
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	3 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	4 "bytes"
971 f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	5 "context"
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	6 "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	7 "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	8 "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	9 "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	10 "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	11 "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	12 "math"
971 f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	13 "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	14 )
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 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	17 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	18 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	19 )
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
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 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	22
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 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	24 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	25 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	26 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	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 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	29 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	30 }
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
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 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	33 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	34
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 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	36
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 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	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	39 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	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 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	42 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	43 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	44 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	45 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	46 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	47 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	48 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	49 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	50 }
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 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	53 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	54
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 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	56 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	57 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	58 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	59 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	60 }
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
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 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	63 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	64 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	65 }
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
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 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	68
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 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	70
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 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	72
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 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	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 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	76 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	77 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	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 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	80 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	81 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	82 }
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 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	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 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	87 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	88
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 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	90
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 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	92 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	93 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	94 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	95 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	96 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	97 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	98 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	99 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	100 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	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
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 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	104 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	105 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	106 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	107 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	108 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	109 }
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
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 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	112 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	113 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	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 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	116
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 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	118 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	119 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	120 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	121 }
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 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	124 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	125 }
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
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 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	128 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	129 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	130 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	131 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	132 }
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 }
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 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	135 // 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	136 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	137 }
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 // 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	139 // 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	140 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	141 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	142 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	143 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	144 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	145 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	146 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	147 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	148 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	149 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	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 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	152 }
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 }
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 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	155 }
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
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 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	158 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	159
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 *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	161 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	162 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	163 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	164 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	165 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	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
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 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	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
971 f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	171 const (
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	172 tinSQLPrefix = `WITH trans AS (
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	173 SELECT
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	174 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	175 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	176 FROM waterway.sounding_results
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	177 `
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	178 tinSQLSuffix = `
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	179 ),
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	180 triangles AS (
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	181 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	182 SELECT (ST_Dump(
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	183 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	184 FROM trans) t
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	185 )
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	186 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	187 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	188
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	189 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	190 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	191 )
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 func GenerateTin(
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	194 conn *sql.Conn,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	195 ctx context.Context,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	196 bottleneck string,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	197 date time.Time,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	198 epsg uint32,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	199 ) (*Tin, error) {
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	200 var tin Tin
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	201 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	202 switch {
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	203 case err == sql.ErrNoRows:
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	204 return nil, nil
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	205 case err != nil:
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	206 return nil, err
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	207 }
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	208 tin.EPSG = epsg
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	209 return &tin, nil
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	210 }
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	211
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	212 func GenerateTinByID(
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	213 conn *sql.Conn,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	214 ctx context.Context,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	215 id int64,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	216 epsg uint32,
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	217 ) (*Tin, error) {
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	218 var tin Tin
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	219 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	220 switch {
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	221 case err == sql.ErrNoRows:
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	222 return nil, nil
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	223 case err != nil:
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	224 return nil, err
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	225 }
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	226 tin.EPSG = epsg
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	227 return &tin, nil
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	228 }
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	229
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	230 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	231 if raw == nil {
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	232 return nil
f9fb6c399f3f Moved generating of tins to octree package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 968 diff changeset	233 }
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	234 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	235 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	236 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	237 }
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	238 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	239 }
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	240
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	241 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	242
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 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	244 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	245 }
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	246
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 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	248 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	249 }
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 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	251 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	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 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	255 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	256 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	257 }
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 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	260 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	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 }
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 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	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 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	267 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	268 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	269 }
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 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	272 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	273 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	274 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	275 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	276 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	277 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	278 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	279 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	280 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	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 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	284 }
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 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	286 }
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 }
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 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	289 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	290
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 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	292 }

Mercurial > gemma

annotate pkg/octree/tin.go @ 971:f9fb6c399f3f