gemma: pkg/models/cross.go annotate

annotate pkg/models/cross.go @ 1234:1a5564655f2a

refac: Sidebar reorganized In order to make context switches between administrative tasks which are map related and those which are system related, we now have a category "administration" and "systemadministration". The Riverbedmorphology does nothing than display the map, so it is renamed to that (map). In case the context of "systemadministration" is chosen, the "map" brings you just back to the map.

author	Thomas Junk <thomas.junk@intevation.de>
date	Tue, 20 Nov 2018 09:54:53 +0100
parents	a244b18cb916
children	19108851af13

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"
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	24 )
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	25
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	26 type (
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	27 GeoJSONFeature string
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	28 GeoJSONLineStringType string
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	29
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	30 GeoJSONDate struct{ time.Time }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	31
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	32 GeoJSONLineString struct {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	33 Type GeoJSONLineStringType `json:"type"`
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
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	36 GeoJSONCoordinate struct {
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	37 Lat float64
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	38 Lon float64
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	39 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	40
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	41 GeoJSONCoordinateZ struct {
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	42 Lat float64
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	43 Lon float64
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	44 Z float64
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	45 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	46
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	47 CrossSectionInputProperties struct {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	48 Bottleneck string `json:"bottleneck"`
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	49 Date GeoJSONDate `json:"date"`
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	50 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	51
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	52 GeoJSONLineCoordinates []GeoJSONCoordinate
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	53 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	54
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	55 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	56
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	57 CrossSectionInputGeometry struct {
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	58 Type GeoJSONLineStringType `json:"type"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	59 Coordinates GeoJSONLineCoordinates `json:"coordinates"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	60 }
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	61
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	62 CrossSectionInput struct {
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	63 Type GeoJSONFeature `json:"type"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	64 Geometry CrossSectionInputGeometry `json:"geometry"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	65 Properties CrossSectionInputProperties `json:"properties"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	66 }
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	67
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	68 CrossSectionOutputGeometry struct {
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	69 Type string `json:"type"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	70 Coordinates GeoJSONMultiLineCoordinatesZ `json:"coordinates"`
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	71 }
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	72
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	73 CrossSectionOutput struct {
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	74 Type string `json:"type"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	75 Geometry CrossSectionOutputGeometry `json:"geometry"`
7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	76 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	77 }
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	78 )
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	79
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	80 var (
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	81 errNoGeoJSONFeature = errors.New("Not a valid GeoJSON feature")
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	82 errNoGeoJSONLineStringType = errors.New("Not a valid GeoJSON linestring type")
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	83 errTooLessComponents = errors.New("Too less components in coordinate")
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	84 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	85 )
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	86
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	87 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	88 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	89 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	90 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	91 }
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 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	93 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	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 *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	96 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	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
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	99 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	100 var s string
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	101 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	102 return err
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	103 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	104 t, err := time.Parse("2006-01-02", s)
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	105 if err != nil {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	106 return err
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	107 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	108 *d = GeoJSONDate{t}
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	109 return nil
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
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	112 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	113 var pos []float64
633 6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	114 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	115 return err
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	116 }
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	117 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	118 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	119 }
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	120 *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	121 return nil
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	122 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	123
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	124 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	125 var s string
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	126 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	127 return err
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	128 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	129 if s != "Feature" {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	130 return errNoGeoJSONFeature
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 *t = GeoJSONFeature(s)
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	133 return nil
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
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	136 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	137 var s string
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	138 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	139 return err
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	140 }
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	141 if s != "LineString" {
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	142 return errNoGeoJSONLineStringType
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 *t = GeoJSONLineStringType(s)
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	145 return nil
6093016fac88 Cross sections: started with implementation of web service (WIP). Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	146 }
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	147
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	148 const (
808 e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	149 wkbXDR byte = 0
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	150 wkbNDR byte = 1
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	151 wkbLineString uint32 = 2
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	152 wkbLineStringZ uint32 = 1000 + 2
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	153 wkbMultiLineStringZ uint32 = 1000 + 5
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
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	156 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	157
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	158 size := 1 + 4 + 4 + len(lc)(1+4+28)
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	159
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	160 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	161
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	162 binary.Write(buf, binary.LittleEndian, wkbNDR)
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	163 binary.Write(buf, binary.LittleEndian, wkbLineString)
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, uint32(len(lc)))
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 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	167 c := &lc[i]
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	168 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	169 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	170 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	171
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	172 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	173 }
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 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	176 var buf bytes.Buffer
653 7aeacd7f150b Cross sections: Take and produce valid GeoJSON. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 652 diff changeset	177 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	178 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	179 }
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 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	182 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	183 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	184 return errNoByteSlice
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	185 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	186 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	187 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	188
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	189 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	190
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	191 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	192
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	193 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	194
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	195 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	196
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	197 switch {
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	198 case err != nil:
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	199 return err
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	200 case endian == wkbNDR:
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	201 order = binary.LittleEndian
643 756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	202 case endian == wkbXDR:
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	203 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	204 default:
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	205 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	206 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	207
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	208 var geomType uint32
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	209 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	210
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	211 switch {
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	212 case err != nil:
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	213 return err
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	214 case geomType != wkbLineStringZ:
645 794592cad75a Cross sections: Fixed SQL, endian, WKB and JSON encoding errors. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 643 diff changeset	215 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	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 var num uint32
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	219 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	220 return err
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	221 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	222
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	223 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	224
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	225 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	226 c := &coords[i]
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	227 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	228 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	229 return err
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 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	232 }
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	233
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	234 *lcz = coords
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	235
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	236 return nil
756f3fc62da6 Cross sections: Using the database in the web service. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 636 diff changeset	237 }
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	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 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	240 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	241 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	242 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	243 )
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	244 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	245 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	246 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	247 }
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	248
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	249 func (cz GeoJSONCoordinateZ) mid(other GeoJSONCoordinateZ) GeoJSONCoordinateZ {
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	250 return GeoJSONCoordinateZ{
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	251 Lon: (cz.Lon + other.Lon) * 0.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	252 Lat: (cz.Lat + other.Lat) * 0.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	253 Z: (cz.Z + other.Z) * 0.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	254 }
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	255 }
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	256
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	257 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	258
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 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	260
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	261 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	262
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	263 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	264 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	265
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	266 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	267 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	268 }
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	269
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	270 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	271 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	272 }
808 e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	273
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	274 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	275 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	276
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	277 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	278
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	279 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	280
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	281 switch {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	282 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	283 return err
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	284 case endian == wkbNDR:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	285 order = binary.LittleEndian
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	286 case endian == wkbXDR:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	287 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	288 default:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	289 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	290 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	291
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	292 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	293 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	294
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	295 switch {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	296 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	297 return err
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	298 case geomType != wkbMultiLineStringZ:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	299 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	300 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	301
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	302 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	303 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	304 return err
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	305 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	306
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	307 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	308
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	309 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	310 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	311
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	312 switch {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	313 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	314 return err
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	315 case endian == wkbNDR:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	316 order = binary.LittleEndian
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	317 case endian == wkbXDR:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	318 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	319 default:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	320 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	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 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	324 switch {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	325 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	326 return err
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	327 case geomType != wkbLineStringZ:
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	328 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	329 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	330
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	331 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	332 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	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
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	336 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	337 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	338 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	339 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	340 return err
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 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	343 return err
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 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	346 return err
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 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	349 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	350 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	351 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	352 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	353 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	354 }
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 *mls = lines
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	357
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	358 return nil
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	359 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	360
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	361 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	362 if src == nil {
3bb275235c89 Let GeoJSONMultiLineCoordinatesZ be nil when scanned from DB. Sascha L. Teichmann <teichmann@intevation.de> parents: 808 diff changeset	363 return nil
3bb275235c89 Let GeoJSONMultiLineCoordinatesZ be nil when scanned from DB. Sascha L. Teichmann <teichmann@intevation.de> parents: 808 diff changeset	364 }
808 e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	365 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	366 if !ok {
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	367 return errNoByteSlice
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	368 }
e10a8a682297 Cross sections: Clip result against area of the sounding result. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	369 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	370 }

Mercurial > gemma

annotate pkg/models/cross.go @ 1234:1a5564655f2a