gemma: pkg/models/cross.go annotate

annotate pkg/models/cross.go @ 4847:4847ac70103a

Made staticcheck happy. - error strings should not end with punctuation or a newline (ST1005) - error strings should not be capitalized (ST1005)

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Mon, 18 Nov 2019 11:54:20 +0100
parents	a347468ea545
children	1222b777f51f

rev	line source
1017 a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	1 // This is Free Software under GNU Affero General Public License v >= 3.0
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	2 // without warranty, see README.md and license for details.
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	3 //
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	4 // SPDX-License-Identifier: AGPL-3.0-or-later
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	5 // License-Filename: LICENSES/AGPL-3.0.txt
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	6 //
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	7 // Copyright (C) 2018 by via donau
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	8 // – Österreichische Wasserstraßen-Gesellschaft mbH
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	9 // Software engineering by Intevation GmbH
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	10 //
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	11 // Author(s):
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	12 // * Sascha L. Teichmann <sascha.teichmann@intevation.de>
a244b18cb916 Added GNU Affero General Public License. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 895 diff changeset	13
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	14 package models
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	15
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	16 import (
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	17 "bytes"
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	18 "encoding/binary"
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	19 "encoding/json"
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	20 "errors"
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	21 "fmt"
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	22 "math"
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	23 "time"
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	24
3094 a347468ea545 Use common.DateFormat for parsing dates everywhere. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2471 diff changeset	25 "gemma.intevation.de/gemma/pkg/common"
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	26 "gemma.intevation.de/gemma/pkg/wkb"
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	27 )
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	28
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	29 type (
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	30 GeoJSONFeature string
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	31 GeoJSONLineStringType string
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	32
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	33 GeoJSONDate struct{ time.Time }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	34
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	35 GeoJSONLineString struct {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	36 Type GeoJSONLineStringType `json:"type"`
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	37 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	38
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	39 GeoJSONCoordinate struct {
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	40 Lat float64
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	41 Lon float64
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	42 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	43
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	44 GeoJSONCoordinateZ struct {
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	45 Lat float64
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	46 Lon float64
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	47 Z float64
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	48 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	49
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	50 CrossSectionInputProperties struct {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	51 Bottleneck string `json:"bottleneck"`
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	52 Date GeoJSONDate `json:"date"`
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	53 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	54
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	55 GeoJSONLineCoordinates []GeoJSONCoordinate
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	56 GeoJSONLineCoordinatesZ []GeoJSONCoordinateZ
636 7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	57
648 cf62cb84fa23 Cross sections: Added a naive O(N^2) algorithm to join the neighbored line segments. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	58 GeoJSONMultiLineCoordinatesZ []GeoJSONLineCoordinatesZ
cf62cb84fa23 Cross sections: Added a naive O(N^2) algorithm to join the neighbored line segments. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	59
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	60 CrossSectionInputGeometry struct {
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	61 Type GeoJSONLineStringType `json:"type"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	62 Coordinates GeoJSONLineCoordinates `json:"coordinates"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	63 }
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	64
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	65 CrossSectionInput struct {
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	66 Type GeoJSONFeature `json:"type"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	67 Geometry CrossSectionInputGeometry `json:"geometry"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	68 Properties CrossSectionInputProperties `json:"properties"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	69 }
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	70
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	71 CrossSectionOutputGeometry struct {
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	72 Type string `json:"type"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	73 Coordinates GeoJSONMultiLineCoordinatesZ `json:"coordinates"`
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	74 }
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	75
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	76 CrossSectionOutput struct {
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	77 Type string `json:"type"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	78 Geometry CrossSectionOutputGeometry `json:"geometry"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	79 Properties map[string]interface{} `json:"properties"`
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	80 }
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	81 )
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	82
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	83 var (
4847 4847ac70103a Made staticcheck happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3094 diff changeset	84 errNoGeoJSONFeature = errors.New("not a valid GeoJSON feature")
4847ac70103a Made staticcheck happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3094 diff changeset	85 errNoGeoJSONLineStringType = errors.New("not a valid GeoJSON linestring type")
4847ac70103a Made staticcheck happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3094 diff changeset	86 errTooLessComponents = errors.New("too less components in coordinate")
4847ac70103a Made staticcheck happy. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3094 diff changeset	87 errTooLessCoordinates = errors.New("too less coordinates")
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	88 )
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	89
636 7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	90 func (lc *GeoJSONLineCoordinates) UnmarshalJSON(data []byte) error {
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	91 var coords []GeoJSONCoordinate
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	92 if err := json.Unmarshal(data, &coords); err != nil {
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	93 return err
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	94 }
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	95 if len(coords) < 2 {
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	96 return errTooLessCoordinates
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	97 }
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	98 *lc = GeoJSONLineCoordinates(coords)
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	99 return nil
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	100 }
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	101
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	102 func (d *GeoJSONDate) UnmarshalJSON(data []byte) error {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	103 var s string
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	104 if err := json.Unmarshal(data, &s); err != nil {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	105 return err
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	106 }
3094 a347468ea545 Use common.DateFormat for parsing dates everywhere. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2471 diff changeset	107 t, err := time.Parse(common.DateFormat, s)
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	108 if err != nil {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	109 return err
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	110 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	111 *d = GeoJSONDate{t}
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	112 return nil
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	113 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	114
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	115 func (c *GeoJSONCoordinate) UnmarshalJSON(data []byte) error {
636 7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	116 var pos []float64
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	117 if err := json.Unmarshal(data, &pos); err != nil {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	118 return err
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	119 }
636 7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	120 if len(pos) < 2 {
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	121 return errTooLessComponents
7a867e48b3c2 Cross sections: Enforce minimal length of coordinates components and coordinates in input linestrings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	122 }
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	123 *c = GeoJSONCoordinate{Lat: pos[0], Lon: pos[1]}
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	124 return nil
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	125 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	126
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	127 func (t *GeoJSONFeature) UnmarshalJSON(data []byte) error {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	128 var s string
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	129 if err := json.Unmarshal(data, &s); err != nil {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	130 return err
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	131 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	132 if s != "Feature" {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	133 return errNoGeoJSONFeature
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	134 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	135 *t = GeoJSONFeature(s)
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	136 return nil
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	137 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	138
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	139 func (t *GeoJSONLineStringType) UnmarshalJSON(data []byte) error {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	140 var s string
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	141 if err := json.Unmarshal(data, &s); err != nil {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	142 return err
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	143 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	144 if s != "LineString" {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	145 return errNoGeoJSONLineStringType
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	146 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	147 *t = GeoJSONLineStringType(s)
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	148 return nil
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	149 }
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	150
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	151 func (lc GeoJSONLineCoordinates) AsWKB() []byte {
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	152
1683 509e8728e846 Fixed pre-calculated size of 2D WKB line strings. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1324 diff changeset	153 size := 1 + 4 + 4 + len(lc)(28)
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	154
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	155 buf := bytes.NewBuffer(make([]byte, 0, size))
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	156
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	157 binary.Write(buf, binary.LittleEndian, wkb.NDR)
63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	158 binary.Write(buf, binary.LittleEndian, wkb.LineString)
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	159 binary.Write(buf, binary.LittleEndian, uint32(len(lc)))
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	160
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	161 for i := range lc {
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	162 c := &lc[i]
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	163 binary.Write(buf, binary.LittleEndian, math.Float64bits(c.Lat))
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	164 binary.Write(buf, binary.LittleEndian, math.Float64bits(c.Lon))
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	165 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	166
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	167 return buf.Bytes()
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	168 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	169
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	170 func (cz GeoJSONCoordinateZ) MarshalJSON() ([]byte, error) {
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	171 var buf bytes.Buffer
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	172 fmt.Fprintf(&buf, "[%.8f,%.8f,%.8f]", cz.Lat, cz.Lon, cz.Z)
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	173 return buf.Bytes(), nil
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	174 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	175
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	176 func (lcz *GeoJSONLineCoordinatesZ) Scan(src interface{}) error {
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	177 data, ok := src.([]byte)
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	178 if !ok {
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	179 return errNoByteSlice
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	180 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	181 return lcz.FromWKB(data)
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	182 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	183
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	184 func (lcz *GeoJSONLineCoordinatesZ) FromWKB(data []byte) error {
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	185
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	186 r := bytes.NewReader(data)
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	187
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	188 endian, err := r.ReadByte()
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	189
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	190 var order binary.ByteOrder
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	191
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	192 switch {
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	193 case err != nil:
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	194 return err
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	195 case endian == wkb.NDR:
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	196 order = binary.LittleEndian
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	197 case endian == wkb.XDR:
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	198 order = binary.BigEndian
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	199 default:
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	200 return fmt.Errorf("unknown byte order %x", endian)
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	201 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	202
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	203 var geomType uint32
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	204 err = binary.Read(r, order, &geomType)
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	205
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	206 switch {
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	207 case err != nil:
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	208 return err
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	209 case geomType != wkb.LineStringZ:
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	210 return fmt.Errorf("unknown geometry type %x", geomType)
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	211 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	212
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	213 var num uint32
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	214 if err = binary.Read(r, order, &num); err != nil {
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	215 return err
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	216 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	217
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	218 coords := make(GeoJSONLineCoordinatesZ, num)
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	219
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	220 for i := range coords {
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	221 c := &coords[i]
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	222 for _, addr := range []*float64{&c.Lat, &c.Lon, &c.Z} {
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	223 if err = binary.Read(r, order, addr); err != nil {
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	224 return err
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	225 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	226 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	227 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	228
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	229 *lcz = coords
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	230
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	231 return nil
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	232 }
648 cf62cb84fa23 Cross sections: Added a naive O(N^2) algorithm to join the neighbored line segments. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	233
651 ce18231825a2 Cross sections: Small details on Equals Method of GeoJSONCoordinateZ (forget to commit). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 648 diff changeset	234 func (cz GeoJSONCoordinateZ) Equals(other GeoJSONCoordinateZ) bool {
648 cf62cb84fa23 Cross sections: Added a naive O(N^2) algorithm to join the neighbored line segments. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	235 const (
cf62cb84fa23 Cross sections: Added a naive O(N^2) algorithm to join the neighbored line segments. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	236 xyEps = 1e-7
cf62cb84fa23 Cross sections: Added a naive O(N^2) algorithm to join the neighbored line segments. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	237 zEps = 1e-5
cf62cb84fa23 Cross sections: Added a naive O(N^2) algorithm to join the neighbored line segments. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	238 )
651 ce18231825a2 Cross sections: Small details on Equals Method of GeoJSONCoordinateZ (forget to commit). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 648 diff changeset	239 return math.Abs(cz.Lat-other.Lat) < xyEps &&
648 cf62cb84fa23 Cross sections: Added a naive O(N^2) algorithm to join the neighbored line segments. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	240 math.Abs(cz.Lon-other.Lon) < xyEps &&
cf62cb84fa23 Cross sections: Added a naive O(N^2) algorithm to join the neighbored line segments. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	241 math.Abs(cz.Z-other.Z) < zEps
cf62cb84fa23 Cross sections: Added a naive O(N^2) algorithm to join the neighbored line segments. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	242 }
cf62cb84fa23 Cross sections: Added a naive O(N^2) algorithm to join the neighbored line segments. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	243
652 f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	244 func deg2rad(d float64) float64 { return d * math.Pi / 180.0 }
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	245
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	246 func (cz GeoJSONCoordinateZ) Distance(other GeoJSONCoordinateZ) float64 {
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	247
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	248 const EarthRadius = 6378137.0
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	249
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	250 dLat := deg2rad(other.Lat - cz.Lat)
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	251 dLng := math.Abs(deg2rad(other.Lon - cz.Lon))
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	252
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	253 if dLng > math.Pi {
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	254 dLng = 2*math.Pi - dLng
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	255 }
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	256
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	257 x := dLng * math.Cos(deg2rad((cz.Lat+other.Lat)/2.0))
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	258 return math.Sqrt(dLatdLat+xx) * EarthRadius
f5ecd1d72a6e Cross section: Replaced naive O(N^2) segment joining with a more clever one. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 651 diff changeset	259 }
808 e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	260
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	261 func (mls *GeoJSONMultiLineCoordinatesZ) FromWKB(data []byte) error {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	262 r := bytes.NewReader(data)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	263
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	264 var order binary.ByteOrder
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	265
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	266 endian, err := r.ReadByte()
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	267
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	268 switch {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	269 case err != nil:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	270 return err
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	271 case endian == wkb.NDR:
808 e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	272 order = binary.LittleEndian
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	273 case endian == wkb.XDR:
808 e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	274 order = binary.BigEndian
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	275 default:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	276 return fmt.Errorf("unknown byte order %x", endian)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	277 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	278
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	279 var geomType uint32
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	280 err = binary.Read(r, order, &geomType)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	281
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	282 switch {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	283 case err != nil:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	284 return err
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	285 case geomType != wkb.MultiLineStringZ:
808 e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	286 return fmt.Errorf("unknown geometry type %d", geomType)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	287 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	288
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	289 var numLines uint32
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	290 if err = binary.Read(r, order, &numLines); err != nil {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	291 return err
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	292 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	293
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	294 lines := make(GeoJSONMultiLineCoordinatesZ, numLines)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	295
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	296 for i := range lines {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	297 endian, err = r.ReadByte()
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	298
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	299 switch {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	300 case err != nil:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	301 return err
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	302 case endian == wkb.NDR:
808 e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	303 order = binary.LittleEndian
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	304 case endian == wkb.XDR:
808 e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	305 order = binary.BigEndian
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	306 default:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	307 return fmt.Errorf("unknown byte order %x", endian)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	308 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	309
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	310 err = binary.Read(r, order, &geomType)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	311 switch {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	312 case err != nil:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	313 return err
2471 63475c8e710e Moved WKB constants to own package. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1683 diff changeset	314 case geomType != wkb.LineStringZ:
808 e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	315 return fmt.Errorf("unknown geometry type %d", geomType)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	316 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	317
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	318 var numPoints uint32
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	319 if err = binary.Read(r, order, &numPoints); err != nil {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	320 return err
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	321 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	322
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	323 points := make(GeoJSONLineCoordinatesZ, numPoints)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	324 for j := range points {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	325 var lat, lon, z uint64
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	326 if err = binary.Read(r, order, &lat); err != nil {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	327 return err
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	328 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	329 if err = binary.Read(r, order, &lon); err != nil {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	330 return err
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	331 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	332 if err = binary.Read(r, order, &z); err != nil {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	333 return err
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	334 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	335 c := &points[j]
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	336 c.Lat = math.Float64frombits(lat)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	337 c.Lon = math.Float64frombits(lon)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	338 c.Z = math.Float64frombits(z)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	339 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	340 lines[i] = points
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	341 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	342
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	343 *mls = lines
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	344
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	345 return nil
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	346 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	347
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	348 func (mls *GeoJSONMultiLineCoordinatesZ) Scan(src interface{}) error {
895 3bb275235c89 Let GeoJSONMultiLineCoordinatesZ be nil when scanned from DB. Sascha L. Teichmann <teichmann@intevation.de> parents: 808 diff changeset	349 if src == nil {
3bb275235c89 Let GeoJSONMultiLineCoordinatesZ be nil when scanned from DB. Sascha L. Teichmann <teichmann@intevation.de> parents: 808 diff changeset	350 return nil
3bb275235c89 Let GeoJSONMultiLineCoordinatesZ be nil when scanned from DB. Sascha L. Teichmann <teichmann@intevation.de> parents: 808 diff changeset	351 }
808 e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	352 data, ok := src.([]byte)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	353 if !ok {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	354 return errNoByteSlice
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	355 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	356 return mls.FromWKB(data)
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	357 }

Mercurial > gemma

annotate pkg/models/cross.go @ 4847:4847ac70103a