gemma: pkg/mesh/raster.go annotate

annotate pkg/mesh/raster.go @ 4889:ab6eb160cd29

Sounding result cache: Start background cache cleanup lazily only when the cache is really used.

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Sun, 09 Feb 2020 20:24:54 +0100
parents	f4abfd0ee8ad
children	34bc6041e61e

rev	line source
4768 a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	1 // This is Free Software under GNU Affero General Public License v >= 3.0
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	2 // without warranty, see README.md and license for details.
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	3 //
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	4 // SPDX-License-Identifier: AGPL-3.0-or-later
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	5 // License-Filename: LICENSES/AGPL-3.0.txt
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	6 //
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	7 // Copyright (C) 2019 by via donau
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	8 // – Österreichische Wasserstraßen-Gesellschaft mbH
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	9 // Software engineering by Intevation GmbH
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	10 //
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11 // Author(s):
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12 // * Sascha L. Teichmann <sascha.teichmann@intevation.de>
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	13
4827 f4abfd0ee8ad Renamed octree package to mesh as there is no octree any more. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4768 diff changeset	14 package mesh
4768 a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16 import (
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17 "log"
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	18 "math"
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	19 "runtime"
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	20 "sync"
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	21 "time"
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	22
4827 f4abfd0ee8ad Renamed octree package to mesh as there is no octree any more. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4768 diff changeset	23 "github.com/fogleman/contourmap"
f4abfd0ee8ad Renamed octree package to mesh as there is no octree any more. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 4768 diff changeset	24
4768 a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	25 "gemma.intevation.de/gemma/pkg/common"
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	26 "gemma.intevation.de/gemma/pkg/wkb"
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	27 )
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29 type Raster struct {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30 BBox Box2D
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	31 CellSize float64
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32 XCells int
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33 YCells int
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	34 Cells []float64
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	35 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	36
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37 const noData = -math.MaxFloat64
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	39 func NewRaster(bbox Box2D, cellSize float64) *Raster {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	40
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	41 width, height := bbox.Size()
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	42
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	43 log.Printf("info raster extent: %.2f / %.2f", width, height)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	44
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	45 xCells := int(math.Ceil(width / cellSize))
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	46 yCells := int(math.Ceil(height / cellSize))
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	47
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	48 log.Printf("info raster size: %d / %d\n", xCells, yCells)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	49
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	50 size := (xCells + 2) * (yCells + 2)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	51 cells := make([]float64, size)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 for i := range cells {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	53 cells[i] = noData
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55 return &Raster{
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56 BBox: bbox,
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57 CellSize: cellSize,
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	58 XCells: xCells,
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	59 YCells: yCells,
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	60 Cells: cells,
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	61 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	62 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	63
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	64 func (r *Raster) Rasterize(eval func(float64, float64) (float64, bool)) {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	65 var wg sync.WaitGroup
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	66
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	67 rows := make(chan int)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69 rasterRow := func() {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70 defer wg.Done()
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	71 quat := 0.25 * r.CellSize
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	72 for i := range rows {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	73 pos := (i+1)*(r.XCells+2) + 1
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	74 row := r.Cells[pos : pos+r.XCells]
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	75 py := r.BBox.Y1 + float64(i)*r.CellSize + r.CellSize/2
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	76 px := r.BBox.X1 + r.CellSize/2
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	77 y1 := py - quat
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	78 y2 := py + quat
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	79 for j := range row {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	80 var n int
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	81 var sum float64
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	82
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	83 if v, ok := eval(px-quat, y1); ok {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	84 sum = v
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	85 n = 1
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	86 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	87 if v, ok := eval(px-quat, y2); ok {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	88 sum += v
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	89 n++
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	90 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	91 if v, ok := eval(px+quat, y1); ok {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	92 sum += v
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	93 n++
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	94 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	95 if v, ok := eval(px+quat, y2); ok {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	96 sum += v
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	97 n++
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	98 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	99
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	100 if n > 0 {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	101 row[j] = sum / float64(n)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	102 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	103 px += r.CellSize
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	104 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	105 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	106 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	107
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	108 for n := runtime.NumCPU(); n >= 1; n-- {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	109 wg.Add(1)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	110 go rasterRow()
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	111 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	112
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	113 for i := 0; i < r.YCells; i++ {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	114 rows <- i
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	115 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	116 close(rows)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	117 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	118
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	119 func (r *Raster) Diff(eval func(float64, float64) (float64, bool)) {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	120 var wg sync.WaitGroup
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	121
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	122 rows := make(chan int)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	123
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	124 rasterRow := func() {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	125 defer wg.Done()
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	126 quat := 0.25 * r.CellSize
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	127 for i := range rows {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	128 pos := (i+1)*(r.XCells+2) + 1
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	129 row := r.Cells[pos : pos+r.XCells]
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	130 py := r.BBox.Y1 + float64(i)*r.CellSize + r.CellSize/2
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	131 px := r.BBox.X1 + r.CellSize/2
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	132 y1 := py - quat
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	133 y2 := py + quat
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	134 for j, old := range row {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	135 // only diff where values
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	136 if old == noData {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	137 px += r.CellSize
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	138 continue
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	139 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140 var n int
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	141 var sum float64
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	143 if v, ok := eval(px-quat, y1); ok {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	144 sum = v
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	145 n = 1
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	146 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	147 if v, ok := eval(px-quat, y2); ok {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148 sum += v
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	149 n++
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	150 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	151 if v, ok := eval(px+quat, y1); ok {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	152 sum += v
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	153 n++
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	154 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	155 if v, ok := eval(px+quat, y2); ok {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	156 sum += v
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	157 n++
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	158 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	160 if n > 0 {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	161 row[j] -= sum / float64(n)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	162 } else {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	163 row[j] = noData
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	164 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	165
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	166 px += r.CellSize
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	167 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	168 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	170
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	171 for n := runtime.NumCPU(); n >= 1; n-- {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	172 wg.Add(1)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	173 go rasterRow()
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	175
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176 for i := 0; i < r.YCells; i++ {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	177 rows <- i
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	178 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	179 close(rows)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	180 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	181
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182 func (r *Raster) ZExtent() (float64, float64, bool) {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	183 min, max := math.MaxFloat64, -math.MaxFloat64
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	184 for _, v := range r.Cells {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	185 if v == noData {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	186 continue
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	187 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	188 if v < min {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	189 min = v
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	190 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	191 if v > max {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	192 max = v
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	193 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	194 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	195 return min, max, min != math.MaxFloat64
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	196 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	197
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	198 func (r *Raster) Trace(heights []float64) []wkb.MultiPolygonGeom {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	199 start := time.Now()
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	200
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	201 tracer := contourmap.FromFloat64s(r.XCells+2, r.YCells+2, r.Cells)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	202
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	203 areas := make([]wkb.MultiPolygonGeom, len(heights))
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	204
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	205 reprojX := common.Linear(0.5, r.BBox.X1, 1.5, r.BBox.X1+r.CellSize)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	206 reprojY := common.Linear(0.5, r.BBox.Y1, 1.5, r.BBox.Y1+r.CellSize)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	207
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	208 var wg sync.WaitGroup
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	209
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	210 cnts := make(chan int)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	211
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	212 doContours := func() {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	213 defer wg.Done()
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214 for hIdx := range cnts {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	215 if c := tracer.Contours(heights[hIdx]); len(c) > 0 {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	216 areas[hIdx] = buildMultipolygon(c, reprojX, reprojY)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	217 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	218 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	219 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	220
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	221 for n := runtime.NumCPU(); n >= 1; n-- {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	222 wg.Add(1)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	223 go doContours()
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	224 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	225
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	226 for i := range heights {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	227 cnts <- i
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	228 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	229 close(cnts)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	230
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	231 wg.Wait()
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	232 log.Printf("info: Tracing areas took %v\n", time.Since(start))
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	233
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	234 return areas
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	235 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	236
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	237 type contour []contourmap.Point
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	238
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	239 type bboxNode struct {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	240 box Box2D
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	241 cnt contour
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	242 children []*bboxNode
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	243 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	244
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	245 func (cnt contour) contains(o contour) bool {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	246 return contains(cnt, o[0].X, o[0].Y) \|\|
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	247 contains(cnt, o[len(o)/2].X, o[len(o)/2].Y)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	248 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	249
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	250 func (cnt contour) length() int {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	251 return len(cnt)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	252 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	253
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	254 func (cnt contour) point(i int) (float64, float64) {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	255 return cnt[i].X, cnt[i].Y
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	256 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	257
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	258 func (cnt contour) bbox() Box2D {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	259 minX, minY := math.MaxFloat64, math.MaxFloat64
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	260 maxX, maxY := -math.MaxFloat64, -math.MaxFloat64
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	261
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	262 for _, p := range cnt {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	263 if p.X < minX {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	264 minX = p.X
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	265 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	266 if p.X > maxX {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	267 maxX = p.X
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	268 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	269 if p.Y < minY {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	270 minY = p.Y
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	271 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	272 if p.Y > maxY {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	273 maxY = p.Y
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	274 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	275 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	276 return Box2D{X1: minX, X2: maxX, Y1: minY, Y2: maxY}
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	277 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	278
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	279 func (bn *bboxNode) insert(cnt contour, box Box2D) {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	280 // check if children are inside new
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	281 var nr *bboxNode
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	282
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	283 for i, r := range bn.children {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	284 if r.box.Inside(box) && cnt.contains(r.cnt) {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	285 if nr == nil {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	286 nr = &bboxNode{box: box, cnt: cnt}
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	287 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	288 nr.children = append(nr.children, r)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	289 bn.children[i] = nil
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	290 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	291 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	292
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	293 // we have a new child
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	294 if nr != nil {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	295 // compact the list
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	296 for i := len(bn.children) - 1; i >= 0; i-- {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	297 if bn.children[i] == nil {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	298 if i < len(bn.children)-1 {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	299 copy(bn.children[i:], bn.children[i+1:])
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	300 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	301 bn.children[len(bn.children)-1] = nil
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	302 bn.children = bn.children[:len(bn.children)-1]
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	303 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	304 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	305 bn.children = append(bn.children, nr)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	306 return
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	307 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	308
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	309 // check if new is inside an old
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	310 for _, r := range bn.children {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	311 if box.Inside(r.box) && r.cnt.contains(cnt) {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	312 r.insert(cnt, box)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	313 return
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	314 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	315 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	316
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	317 // its a new child node.
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	318 nr = &bboxNode{box: box, cnt: cnt}
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	319 bn.children = append(bn.children, nr)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	320 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	321
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	322 func (bn *bboxNode) insertRoot(cnt contour) {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	323 bn.insert(cnt, cnt.bbox())
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	324 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	325
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	326 type bboxOutFunc func(contour, []contour)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	327
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	328 func (bn *bboxNode) generate(out bboxOutFunc) {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	329
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	330 var grands []*bboxNode
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	331
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	332 holes := make([]contour, len(bn.children))
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	333
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	334 for i, ch := range bn.children {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	335 holes[i] = ch.cnt
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	336 grands = append(grands, ch.children...)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	337 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	338 out(bn.cnt, holes)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	339
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	340 // the grand children are new polygons.
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	341 for _, grand := range grands {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	342 grand.generate(out)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	343 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	344 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	345
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	346 func (bn *bboxNode) generateRoot(out bboxOutFunc) {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	347 for _, r := range bn.children {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	348 r.generate(out)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	349 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	350 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	351
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	352 func buildMultipolygon(
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	353 cnts []contourmap.Contour,
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	354 reprojX, reprojY func(float64) float64,
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	355 ) wkb.MultiPolygonGeom {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	356
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	357 var forest bboxNode
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	358
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	359 for _, cnt := range cnts {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	360 forest.insertRoot(contour(cnt))
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	361 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	362
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	363 //log.Printf("cnts: %d roots: %d\n", len(cnts), len(bf.roots))
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	364
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	365 var mp wkb.MultiPolygonGeom
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	366
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	367 out := func(sh contour, hls []contour) {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	368
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	369 polygon := make(wkb.PolygonGeom, 1+len(hls))
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	370
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	371 // Handle shell
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	372 shell := make(wkb.LinearRingGeom, len(sh))
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	373 for j, pt := range sh {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	374 shell[j] = wkb.PointGeom{
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	375 X: reprojX(pt.X),
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	376 Y: reprojY(pt.Y),
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	377 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	378 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	379 if shell.CCW() {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	380 shell.Reverse()
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	381 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	382 polygon[0] = shell
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	383
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	384 // handle holes
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	385 for i, hl := range hls {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	386 hole := make(wkb.LinearRingGeom, len(hl))
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	387 for j, pt := range hl {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	388 hole[j] = wkb.PointGeom{
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	389 X: reprojX(pt.X),
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	390 Y: reprojY(pt.Y),
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	391 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	392 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	393 if !hole.CCW() {
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	394 hole.Reverse()
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	395 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	396 polygon[1+i] = hole
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	397 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	398
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	399 mp = append(mp, polygon)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	400 }
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	401
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	402 forest.generateRoot(out)
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	403
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	404 return mp
a2f16bbcc846 Morph differences: Directly raster A and subtract B as a raster. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	405 }

Mercurial > gemma

annotate pkg/mesh/raster.go @ 4889:ab6eb160cd29