gemma: pkg/mesh/meshserialize

annotate pkg/mesh/meshserialize_v2.go @ 5694:3bc15e38c7e8 sr-v2

Typo fix

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Mon, 12 Feb 2024 14:43:31 +0100
parents	d920f0fa2f04
children	ef80748ae4f3

rev	line source
5680 a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	1 // This is Free Software under GNU Affero General Public License v >= 3.0
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	2 // without warranty, see README.md and license for details.
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	3 //
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	4 // SPDX-License-Identifier: AGPL-3.0-or-later
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	5 // License-Filename: LICENSES/AGPL-3.0.txt
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	6 //
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	7 // Copyright (C) 2024 by via donau
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	8 // – Österreichische Wasserstraßen-Gesellschaft mbH
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	9 // Software engineering by Intevation GmbH
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	10 //
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11 // Author(s):
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12 // * Sascha L. Teichmann <sascha.teichmann@intevation.de>
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	13
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	14 package mesh
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15
5684 536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	16 import (
5686 c33a5354328d Add stubs for de/serialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5685 diff changeset	17 "bufio"
5684 536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	18 "cmp"
5687 8ff842858434 More stubs for serializing v2. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5686 diff changeset	19 "encoding/binary"
5686 c33a5354328d Add stubs for de/serialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5685 diff changeset	20 "errors"
c33a5354328d Add stubs for de/serialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5685 diff changeset	21 "io"
5684 536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	22 "math"
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	23 "slices"
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	24 "sort"
5688 6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	25
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	26 "gemma.intevation.de/gemma/pkg/common"
5684 536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	27 )
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	28
5688 6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	29 const quantScale = 1_000
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	30
5690 da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	31 func quant(x float64) int64 { return int64(math.Round(x * quantScale)) }
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	32 func unquant(x int64) float64 { return float64(x) / quantScale }
5688 6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	33
5686 c33a5354328d Add stubs for de/serialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5685 diff changeset	34 func (s *STRTree) serializeV2(w io.Writer) error {
5692 d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	35 return serializer(w, s,
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	36 serializeHeader(),
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	37 serializeVersion(2),
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	38 (*STRTree).serializeTinV2,
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	39 (*STRTree).serializeEntries,
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	40 (*STRTree).serializeIndexV2,
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	41 (*STRTree).serializeBBoxesV2,
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	42 )
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	43 }
5690 da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	44
5692 d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	45 func (s *STRTree) serializeTinV2(w io.Writer) error {
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	46 return s.tin.serializeV2(w)
5687 8ff842858434 More stubs for serializing v2. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5686 diff changeset	47 }
8ff842858434 More stubs for serializing v2. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5686 diff changeset	48
8ff842858434 More stubs for serializing v2. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5686 diff changeset	49 func (t *Tin) serializeV2(w io.Writer) error {
5692 d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	50 return serializer(w, t,
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	51 (*Tin).serializeEPSG,
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	52 (*Tin).serializeExtent,
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	53 (*Tin).serializeVerticesV2,
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	54 (*Tin).serializeTrianglesV2,
d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	55 )
5688 6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	56 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	57
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	58 func zigZag(buf []byte, w io.Writer, x int64) error {
5689 159291854063 Simplify zug zig store call. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5688 diff changeset	59 _, err := w.Write(buf[:binary.PutVarint(buf, x)])
159291854063 Simplify zug zig store call. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5688 diff changeset	60 return err
5688 6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	61 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	62
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	63 func (t *Tin) serializeVerticesV2(w io.Writer) error {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	64 vertices := t.Vertices
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	65 if err := binary.Write(w, binary.LittleEndian, uint32(len(vertices))); err != nil {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	66 return err
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	67 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	68 m := t.Min
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	69 buf := make([]byte, binary.MaxVarintLen64)
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	70 for dim := range 3 {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	71 delta := common.Delta()
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	72 for i := range vertices {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	73 var v float64
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	74 switch dim {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	75 case 0:
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	76 v = vertices[i].X - m.X
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	77 case 1:
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	78 v = vertices[i].Y - m.Y
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	79 case 2:
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	80 v = vertices[i].Z - m.Z
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	81 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	82 q := quant(v)
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	83 d := delta(q)
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	84 if err := zigZag(buf, w, d); err != nil {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	85 return err
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	86 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	87 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	88 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	89 return nil
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	90 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	91
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	92 func (t *Tin) serializeTrianglesV2(w io.Writer) error {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	93 triangles := t.Triangles
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	94 if err := binary.Write(w, binary.LittleEndian, uint32(len(triangles))); err != nil {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	95 return err
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	96 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	97 buf := make([]byte, binary.MaxVarintLen64)
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	98 delta := common.Delta()
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	99 for _, tri := range triangles {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	100 for _, idx := range tri {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	101 d := delta(int64(idx))
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	102 if err := zigZag(buf, w, d); err != nil {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	103 return err
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	104 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	105 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	106 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	107 return nil
5687 8ff842858434 More stubs for serializing v2. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5686 diff changeset	108 }
8ff842858434 More stubs for serializing v2. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5686 diff changeset	109
8ff842858434 More stubs for serializing v2. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5686 diff changeset	110 func (s *STRTree) serializeIndexV2(w io.Writer) error {
5688 6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	111 index := s.index
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	112 if err := binary.Write(w, binary.LittleEndian, uint32(len(index))); err != nil {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	113 return err
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	114 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	115 buf := make([]byte, binary.MaxVarintLen64)
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	116 delta := common.Delta()
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	117 for _, idx := range index {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	118 d := delta(int64(idx))
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	119 if err := zigZag(buf, w, d); err != nil {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	120 return err
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	121 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	122 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	123 return nil
5687 8ff842858434 More stubs for serializing v2. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5686 diff changeset	124 }
8ff842858434 More stubs for serializing v2. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5686 diff changeset	125
8ff842858434 More stubs for serializing v2. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5686 diff changeset	126 func (s *STRTree) serializeBBoxesV2(w io.Writer) error {
5688 6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	127 boxes := s.bboxes
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	128 if err := binary.Write(w, binary.LittleEndian, uint32(len(boxes))); err != nil {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	129 return err
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	130 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	131 var (
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	132 m = s.Min()
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	133 x1 = func(b Box2D) int64 { return quant(b.X1 - m.X) }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	134 y1 = func(b Box2D) int64 { return quant(b.Y1 - m.Y) }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	135 width = func(b Box2D) int64 { return quant(b.X2 - b.X1) }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	136 height = func(b Box2D) int64 { return quant(b.Y2 - b.Y1) }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	137 buf = make([]byte, binary.MaxVarintLen64)
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	138 )
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	139 for _, proj := range []func(Box2D) int64{x1, y1, width, height} {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	140 delta := common.Delta()
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	141 for _, b := range boxes {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	142 v := delta(proj(b))
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	143 if err := zigZag(buf, w, v); err != nil {
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	144 return err
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	145 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	146 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	147 }
6281c18b109f Finsh serializing v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5687 diff changeset	148 return nil
5686 c33a5354328d Add stubs for de/serialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5685 diff changeset	149 }
c33a5354328d Add stubs for de/serialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5685 diff changeset	150
c33a5354328d Add stubs for de/serialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5685 diff changeset	151 func (s STRTree) deserializeV2(r bufio.Reader) error {
5690 da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	152 // The header is already skipped here.
5692 d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	153 return serializer(r, s,
5691 9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	154 (*STRTree).deserializeTinV2,
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	155 (*STRTree).deserializeEntries,
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	156 (*STRTree).deserializeIndexV2,
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	157 (*STRTree).deserializeBBoxesV2,
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	158 )
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	159 }
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	160
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	161 func (s STRTree) deserializeTinV2(r bufio.Reader) error {
5690 da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	162 s.tin = new(Tin)
5691 9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	163 return s.tin.deserializeV2(r)
5690 da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	164 }
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	165
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	166 func (s STRTree) deserializeBBoxesV2(r bufio.Reader) error {
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	167 // TODO: implement me!
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	168 return errors.New("not implemented, yet")
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	169 }
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	170
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	171 func (s STRTree) deserializeIndexV2(r bufio.Reader) error {
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	172 // TODO: implement me!
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	173 return errors.New("not implemented, yet")
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	174 }
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	175
da9720b4aa42 More stubs and refactoring to deserialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5689 diff changeset	176 func (t Tin) deserializeV2(r bufio.Reader) error {
5692 d920f0fa2f04 User generic serializer for serialization, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5691 diff changeset	177 return serializer(r, t,
5691 9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	178 (*Tin).deserializeEPSG,
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	179 (*Tin).deserializeExtent,
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	180 (*Tin).deserializeVerticesV2,
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	181 (*Tin).deserializeTrianglesV2,
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	182 )
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	183 }
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	184
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	185 func (t Tin) deserializeVerticesV2(r bufio.Reader) error {
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	186 // TODO: implement me!
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	187 return errors.New("not implemented, yet")
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	188 }
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	189
9d2e74225104 Some refactoring to clean up deserialization. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5690 diff changeset	190 func (t Tin) deserializeTrianglesV2(r bufio.Reader) error {
5686 c33a5354328d Add stubs for de/serialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5685 diff changeset	191 // TODO: implement me!
c33a5354328d Add stubs for de/serialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5685 diff changeset	192 return errors.New("not implemented, yet")
c33a5354328d Add stubs for de/serialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5685 diff changeset	193 }
c33a5354328d Add stubs for de/serialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5685 diff changeset	194
c33a5354328d Add stubs for de/serialize v2 meshes. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5685 diff changeset	195 // optimizeForSerializationV2 prepares the mesh to be stored in V2.
5680 a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	196 func (s *STRTree) optimizeForSerializationV2() {
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	197 s.removeUnusedTriangles()
5684 536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	198 s.sortVertices()
5682 33499bd1b829 Sort indices in spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5680 diff changeset	199 s.sortIndices()
5680 a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	200 }
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	201
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	202 // removeUnusedTriangles removes all triangles from the
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	203 // TIN that are not registered in the spatial index.
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	204 func (s *STRTree) removeUnusedTriangles() {
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	205 used := make([]bool, len(s.tin.Triangles))
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	206 unused := len(used)
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	207 s.allTriangleIndices(func(indices []int32) {
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	208 for _, idx := range indices {
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	209 used[idx] = true
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	210 unused-- // They only occur once in the spatial index.
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	211 }
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	212 })
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	213 if unused <= 0 {
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214 return
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	215 }
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	216 newTris := make([][]int32, 0, len(s.tin.Triangles)-unused)
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	217 remap := map[int32]int32{}
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	218 for idx, tri := range s.tin.Triangles {
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	219 if used[idx] {
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	220 remap[int32(idx)] = int32(len(newTris))
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	221 newTris = append(newTris, tri)
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	222 }
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	223 }
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	224 s.tin.Triangles = newTris
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	225 // Update the spatial index as the gaps are now closed.
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	226 s.allTriangleIndices(func(indices []int32) {
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	227 for i, idx := range indices {
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	228 indices[i] = remap[idx]
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	229 }
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	230 })
a87900c0fd11 Remove triangles that are not registered in the spatial index. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	231 }
5684 536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	232
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	233 // sortVertices so that the deltas between neighbors are minimized.
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	234 func (s *STRTree) sortVertices() {
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	235 vertices := s.tin.Vertices
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	236 indices := make([]int32, len(vertices))
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	237 for i := range len(indices) {
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	238 indices[i] = int32(i)
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	239 }
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	240
5685 0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	241 // Take aspect ratio of bbox into account.
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	242 width := s.tin.Max.X - s.tin.Min.X
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	243 height := s.tin.Max.Y - s.tin.Min.Y
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	244 if width == 0 \|\| height == 0 {
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	245 return
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	246 }
5684 536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	247
5685 0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	248 // v = len(vertices)
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	249 // r = max(w, h)/min(w, h)
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	250 // n * t = v <=> n = v/t
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	251 // n / t = r <=> n = r*t
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	252 // v/t = r*t
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	253 // v = r*t^2
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	254 // v/r = t^2
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	255 // t = sqrt(v/r)
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	256 r := float64(max(width, height)) / float64(min(width, height))
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	257 t := int(math.Ceil(math.Sqrt(float64(len(vertices)) / r)))
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	258
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	259 var d1, d2 func(int32, int32) int
5684 536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	260
5685 0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	261 if width > height {
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	262 // Sort in X first and slices alternating up and down
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	263 slices.SortFunc(indices, func(a, b int32) int {
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	264 return cmp.Compare(vertices[a].X, vertices[b].X)
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	265 })
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	266 d1 = func(a, b int32) int { return cmp.Compare(vertices[a].Y, vertices[b].Y) }
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	267 d2 = func(a, b int32) int { return cmp.Compare(vertices[b].Y, vertices[a].Y) }
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	268 } else {
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	269 // Sort in Y first and slices alternating left and right
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	270 slices.SortFunc(indices, func(a, b int32) int {
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	271 return cmp.Compare(vertices[a].Y, vertices[b].Y)
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	272 })
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	273 d1 = func(a, b int32) int { return cmp.Compare(vertices[a].X, vertices[b].X) }
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	274 d2 = func(a, b int32) int { return cmp.Compare(vertices[b].X, vertices[a].X) }
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	275 }
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	276
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	277 for p := indices; len(p) > 0; d1, d2 = d2, d1 {
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	278 l := min(len(p), t)
0ee8ace01b60 When sorting vertices by x/y take aspectio ratio of bbox into acount. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5684 diff changeset	279 slices.SortStableFunc(p[:l], d2)
5684 536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	280 p = p[l:]
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	281 }
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	282
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	283 // We really need to to sort the vertices actually, too.
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	284 sort.Slice(vertices, func(i, j int) bool { return indices[i] < indices[j] })
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	285
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	286 reindices := make([]int32, len(indices))
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	287 for i, idx := range indices {
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	288 reindices[idx] = int32(i)
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	289 }
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	290
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	291 // Update the indices in the triangles.
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	292 for _, tri := range s.tin.Triangles {
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	293 for i, idx := range tri {
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	294 tri[i] = reindices[idx]
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	295 }
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	296 }
536e842d9bfa Reorder vertices in tins to minimize delta distances. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 5682 diff changeset	297 }

Mercurial > gemma

annotate pkg/mesh/meshserialize_v2.go @ 5694:3bc15e38c7e8 sr-v2