Mercurial > gemma

changeset 4547:6247f5a42226 iso-areas

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Increased resolution.
author Sascha L. Teichmann <sascha.teichmann@intevation.de>
date Fri, 27 Sep 2019 18:20:13 +0200
parents a3f1d92b8597
children befb94e3a953
files cmd/isoareas/main.go
diffstat 1 files changed, 13 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/cmd/isoareas/main.go	Fri Sep 27 17:26:43 2019 +0200
+++ b/cmd/isoareas/main.go	Fri Sep 27 18:20:13 2019 +0200
@@ -134,6 +134,12 @@
 
 func connectArcs(tri *octree.Triangulation, cuts []indexedArc, arcs *[]octree.LineStringZ) {
 
+	unique := map[int32]struct{}{}
+	for i := range cuts {
+		unique[cuts[i].arc] = struct{}{}
+	}
+	before := len(unique)
+
 	origLen := int32(len(*arcs))
 
 	for i := range cuts {
@@ -189,8 +195,14 @@
 				cuts[j].arc = nidx
 			}
 		}
+	}
 
+	unique = map[int32]struct{}{}
+	for i := range cuts {
+		unique[cuts[i].arc] = struct{}{}
 	}
+	log.Printf("unique arcs: before: %d after %d\n",
+		before, len(unique))
 }
 
 func intersectTriangles(tri *octree.Triangulation, heights []float64) [][]octree.LineStringZ {
@@ -442,7 +454,7 @@
 
 	log.Printf("ratio: %.2f\n", ratio)
 
-	const width = 1000
+	const width = 50000
 	height := int(math.Ceil(width * ratio))
 
 	px := linear(min.X, 0, max.X, width)