gemma: cmd/isoareas/main.go comparison

Increased resolution.

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Fri, 27 Sep 2019 18:20:13 +0200
parents	a3f1d92b8597
children	befb94e3a953

comparison

equal deleted inserted replaced

-:a3f1d92b8597
+:6247f5a42226
 	return nil
 }
 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 {
 		if cuts[i].arc >= origLen {
 			// already has a connected arc assigned.
 			cuts[i].arc = nidx
 			for _, j := range modified {
 				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 {
 	cuts := make([][]indexedArc, len(heights))
 	ratio := (max.X - min.X) / (max.Y - min.Y)
 	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)
 	py := linear(min.Y, float64(height), max.Y, 0)

Mercurial > gemma