dive4elements/gnv-client: gnv-artifacts/src/main/java/de/intevation/gnv/raster/Vectorizer.java annotate

annotate gnv-artifacts/src/main/java/de/intevation/gnv/raster/Vectorizer.java @ 458:92d6cf448598

Improved raster tile based height evaluation. gnv-artifacts/trunk@512 c6561f87-3c4e-4783-a992-168aeb5c3f6f

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Wed, 06 Jan 2010 10:09:51 +0000
parents	92b7ccbf6163
children	9a828e5a2390

rev	line source
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	1 package de.intevation.gnv.raster;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	2
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	3 import java.util.BitSet;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	4 import java.util.List;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	5 import java.util.ArrayList;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	6
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	7 import gnu.trove.TIntStack;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	8 import gnu.trove.TIntObjectHashMap;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	9
445 f42ed4f10b79 Fixed some bugs and create "Profilschnitt" polygons via configuration. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 437 diff changeset	10 import org.apache.log4j.Logger;
f42ed4f10b79 Fixed some bugs and create "Profilschnitt" polygons via configuration. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 437 diff changeset	11
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12 /**
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	13 * @author Sascha L. Teichmann (sascha.teichmann@intevation.de)
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	14 */
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15 public class Vectorizer
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16 {
445 f42ed4f10b79 Fixed some bugs and create "Profilschnitt" polygons via configuration. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 437 diff changeset	17 private static Logger log = Logger.getLogger(Vectorizer.class);
f42ed4f10b79 Fixed some bugs and create "Profilschnitt" polygons via configuration. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 437 diff changeset	18
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	19 public interface RingsHandler {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	20
436 6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	21 void handleRings(
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	22 List<Edge> rings,
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	23 int value,
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	24 int width,
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	25 int height);
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	26
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	27 } // interface RingsHandler
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29 public static final class Edge {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	31 public Edge prev;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32 public Edge next;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	34 public int a;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	35 public int b;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	36
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37 public Edge() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	39
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	40 public Edge(int a, int b) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	41 this.a = a;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	42 this.b = b;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	43 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	44
437 b624879d2902 Added vectorizer rings callback to generate iso lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 436 diff changeset	45 public Edge(Edge other) {
b624879d2902 Added vectorizer rings callback to generate iso lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 436 diff changeset	46 a = other.a;
b624879d2902 Added vectorizer rings callback to generate iso lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 436 diff changeset	47 b = other.b;
b624879d2902 Added vectorizer rings callback to generate iso lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 436 diff changeset	48 }
b624879d2902 Added vectorizer rings callback to generate iso lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 436 diff changeset	49
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	50 public void chain(Edge other, int found) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	51
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 if (found == a) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	53 other.next = this;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54 prev = other;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55 return;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	58 if (found == b) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	59 next = other;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	60 other.prev = this;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	61 return;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	62 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	63
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	64 throw new IllegalStateException("cannot chain");
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	65 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	66
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	67 public boolean isComplete() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68 Edge current = this;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69 do {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70 if (current.prev == null \|\| current.next == null) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	71 return false;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	72 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	73 current = current.next;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	74 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	75 while (current != this);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	76 return true;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	77 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	78
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	79 public int length() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	80 int length = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	81 Edge current = this;
437 b624879d2902 Added vectorizer rings callback to generate iso lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 436 diff changeset	82 do { ++length; }
b624879d2902 Added vectorizer rings callback to generate iso lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 436 diff changeset	83 while ((current = current.next) != null && current != this);
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	84 return length;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	85 }
436 6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	86
447 92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	87 public Edge head() {
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	88 Edge current = this;
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	89 while (current.prev != null) {
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	90 current = current.prev;
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	91 }
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	92 return current;
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	93 }
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	94
436 6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	95 public int hashCode() {
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	96 return (a << 16) \| b;
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	97 }
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	98
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	99 public boolean equals(Object other) {
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	100 Edge e = (Edge)other;
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	101 return a == e.a && b == e.b;
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	102 }
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	103 } // class Edge
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	104
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	105 protected static Edge simplify(Edge edge, int width) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	106
447 92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	107 Edge e1 = edge, start = edge;
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	108
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	109 int length = edge.length();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	110
437 b624879d2902 Added vectorizer rings callback to generate iso lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 436 diff changeset	111 if (length < 2) {
b624879d2902 Added vectorizer rings callback to generate iso lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 436 diff changeset	112 return e1;
b624879d2902 Added vectorizer rings callback to generate iso lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 436 diff changeset	113 }
b624879d2902 Added vectorizer rings callback to generate iso lines. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 436 diff changeset	114
447 92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	115 Edge e2 = edge.next;
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	116
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	117 int count = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	118
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	119 do {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	120 int e1x = e1.a % width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	121 int e1y = e1.a / width;
447 92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	122 int e2x = e1.b % width;
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	123 int e2y = e1.b / width;
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	124 int e3x = e2.b % width;
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	125 int e3y = e2.b / width;
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	126
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	127 if ((e1x == e2x && e2x == e3x && e1x == e3x)
447 92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	128 \|\| (e1y == e2y && e2y == e3y && e1y == e3y)) {
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	129 e1.b = e2.b;
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	130 Edge removed = e1.next;
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	131 e1.next = e2.next;
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	132 if (e1.next != null) {
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	133 e1.next.prev = e1;
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	134 }
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	135 e2 = e1.next;
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	136 count = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	137 --length;
447 92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	138 if (removed == start) {
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	139 start = e1;
92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	140 }
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	141 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	143 e1 = e1.next;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	144 e2 = e2.next;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	145 ++count;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	146 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	147 }
447 92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	148 while (length > 1 && e2 != null && count < length + 2);
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	149
447 92b7ccbf6163 Improved generation of iso lines in vertical cross section. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 445 diff changeset	150 return start;
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	151 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	152
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	153 protected int [] raster;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	154 protected int width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	155 protected TIntObjectHashMap openEdges;
436 6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	156 protected List<Edge> rings;
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	157 protected boolean simplify;
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	158
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159 public Vectorizer() {
436 6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	160 this(true);
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	161 }
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	162
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	163 public Vectorizer(boolean simplify) {
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	164 openEdges = new TIntObjectHashMap();
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	165 rings = new ArrayList<Edge>();
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	166 this.simplify = simplify;
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	167 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	168
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169 public Vectorizer(int [] raster, int width) {
445 f42ed4f10b79 Fixed some bugs and create "Profilschnitt" polygons via configuration. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 437 diff changeset	170 this(true, raster, width);
f42ed4f10b79 Fixed some bugs and create "Profilschnitt" polygons via configuration. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 437 diff changeset	171 }
f42ed4f10b79 Fixed some bugs and create "Profilschnitt" polygons via configuration. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 437 diff changeset	172
f42ed4f10b79 Fixed some bugs and create "Profilschnitt" polygons via configuration. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 437 diff changeset	173 public Vectorizer(boolean simplify, int [] raster, int width) {
f42ed4f10b79 Fixed some bugs and create "Profilschnitt" polygons via configuration. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 437 diff changeset	174 this(simplify);
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	175 this.raster = raster;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176 this.width = width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	177 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	178
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	179 public static final int tl(int i, int w) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	180 int x = i % w;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	181 int y = i / w;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182 return x + (w + 1)*y;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	183 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	184
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	185 public static final int tr(int i, int w) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	186 return tl(i, w) + 1;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	187 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	188
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	189 public static final int bl(int i, int w) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	190 return tl(i, w) + w + 1;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	191 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	192
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	193 public static final int br(int i, int w) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	194 return bl(i, w) + 1;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	195 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	196
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	197 protected void resetRegion() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	198 openEdges.clear();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	199 rings.clear();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	200 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	201
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	202 protected void emit(Edge edge) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	203
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	204 Edge otherA = (Edge)openEdges.remove(edge.a);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	205 if (otherA != null) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	206 otherA.chain(edge, edge.a);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	207 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	208
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	209 Edge otherB = (Edge)openEdges.remove(edge.b);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	210 if (otherB != null) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	211 otherB.chain(edge, edge.b);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	212 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	213
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214 if (edge.isComplete()) {
436 6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	215 rings.add(simplify ? simplify(edge, width + 1) : edge);
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	216 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	217 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	218 if (otherA == null) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	219 openEdges.put(edge.a, edge);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	220 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	221 if (otherB == null) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	222 openEdges.put(edge.b, edge);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	223 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	224 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	225 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	226
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	227 public int process(RingsHandler handler) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	228
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	229 BitSet visited = new BitSet(raster.length);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	230
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	231 TIntStack stack = new TIntStack();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	232
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	233 int regions = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	234
436 6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	235 int height = raster.length / width;
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	236
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	237 for (int i = 0; i < raster.length; ++i) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	238 if (visited.get(i)) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	239 continue;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	240 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	241
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	242 ++regions;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	243
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	244 int currentValue = raster[i];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	245 visited.set(i);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	246
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	247 int current = i;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	248 visited.set(current);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	249
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	250 for (;;) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	251 int tl = tl(current, width);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	252 int tr = tr(current, width);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	253 int bl = bl(current, width);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	254 int br = br(current, width);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	255
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	256 int t = current - width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	257
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	258 if (t < 0) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	259 emit(new Edge(tr, tl));
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	260 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	261 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	262 if (raster[t] != currentValue) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	263 emit(new Edge(tr, tl));
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	264 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	265 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	266 if (!visited.get(t)) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	267 visited.set(t);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	268 stack.push(t);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	269 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	270 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	271 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	272
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	273 int b = current + width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	274
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	275 if (b >= raster.length) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	276 emit(new Edge(bl, br));
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	277 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	278 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	279 if (raster[b] != currentValue) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	280 emit(new Edge(bl, br));
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	281 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	282 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	283 if (!visited.get(b)) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	284 visited.set(b);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	285 stack.push(b);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	286 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	287 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	288 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	289
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	290 int x = current % width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	291
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	292 if (x == 0) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	293 emit(new Edge(tl, bl));
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	294 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	295 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	296 int l = current - 1;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	297 if (raster[l] != currentValue) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	298 emit(new Edge(tl, bl));
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	299 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	300 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	301 if (!visited.get(l)) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	302 visited.set(l);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	303 stack.push(l);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	304 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	305 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	306 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	307
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	308 if (x == width - 1) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	309 emit(new Edge(br, tr));
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	310 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	311 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	312 int r = current + 1;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	313 if (raster[r] != currentValue) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	314 emit(new Edge(br, tr));
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	315 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	316 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	317 if (!visited.get(r)) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	318 visited.set(r);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	319 stack.push(r);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	320 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	321 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	322 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	323
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	324 if (stack.size() == 0) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	325 break;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	326 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	327
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	328 current = stack.pop();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	329 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	330
436 6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	331 handler.handleRings(
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	332 rings,
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	333 currentValue,
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	334 width + 1,
6642ab6c583c Added vectorizer rings callback which generates polygon datasets suitable Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	335 height + 1);
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	336
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	337 resetRegion();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	338 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	339
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	340 return regions;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	341 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	342 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	343 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8:

Mercurial > dive4elements > gnv-client

annotate gnv-artifacts/src/main/java/de/intevation/gnv/raster/Vectorizer.java @ 458:92d6cf448598