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

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Mon, 21 Dec 2009 18:00:54 +0000
parents
children	6e8364e766fa

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 /**
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	4 * @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	5 */
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	6 public class Raster
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	7 {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	8 public static class Kernel {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	9
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	10 public static final class Coeff {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11 public int di;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12 public int dj;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	13 public double c;
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 Coeff() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	18 public Coeff(int di, int dj, double c) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	19 this.di = di;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	20 this.dj = dj;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	21 this.c = c;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	22 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	23
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	24 public String toString() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	25 return String.valueOf(c);
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
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28 } // class Coeff
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30 protected Coeff [] coeffs;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	31 protected int width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33 public Kernel() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	34 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	35
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	36 public Kernel(Coeff [] coeffs, int width) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37 this.coeffs = coeffs;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38 this.width = width;
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
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	41 public static double gauss(double x, double y, double sigma) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	42 double s2 = sigma*sigma;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	43 return 1.0d/(2.0dMath.PIs2)Math.exp(-(xx + yy)/(2.0ds2));
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	44 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	45
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	46 public static Kernel createGauss() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	47 return createGauss(1.0d, 5);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	48 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	49
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	50 public static Kernel createGauss(double sigma, int width) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	51 Coeff [] coeffs = new Coeff[width*width];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 double sum = 0d;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	53 for (int j = 0; j < width; ++j) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54 int y = -width/2 + j;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55 for (int i = 0; i < width; ++i) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56 int x = -width/2 + i;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57 double c = gauss(x, y, sigma);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	58 coeffs[j*width + i] = new Coeff(x, y, c);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	59 sum += c;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	60 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	61 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	62 sum = 1.0/sum;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	63 for (int i = 0; i < coeffs.length; ++i) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	64 coeffs[i].c *= sum;
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 return new Kernel(coeffs, width);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	67 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69 public double fold(Access access, int i, int j) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	71 double s = 0.0d;
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 double unused = 0d;
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 for (int k = 0; k < coeffs.length; ++k) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	76 Coeff coeff = coeffs[k];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	77 double v = access.get(i + coeff.di, j + coeff.dj);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	78 if (Double.isNaN(v)) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	79 unused += coeff.c;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	80 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	81 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	82 s += v * coeff.c;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	83 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	84 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	85
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	86 return s*(1.0d - unused);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	87 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	88
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	89 public String toString() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	90 StringBuilder sb = new StringBuilder();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	91
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	92 int height = coeffs.length/width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	93
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	94 for (int j = 0; j < height; ++j) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	95 if (j > 0) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	96 sb.append(System.getProperty("line.separator"));
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	97 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	98 for (int i = 0; i < width; ++i) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	99 if (i > 0) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	100 sb.append("\t");
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	101 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	102 sb.append(coeffs[j*width + i]);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	103 }
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
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	106 return sb.toString();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	107 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	108 } // class Kernel
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	109
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	110 public interface Access {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	111
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	112 double get(int i, int j);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	113
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	114 } // interface Access
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	115
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	116 public interface ValueToIndex {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	117
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	118 int toIndex(double value);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	119
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	120 } // interface ValueToIndex;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	121
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	122 public static class IsoClasses
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	123 implements ValueToIndex
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	124 {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	125 protected double m;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	126 protected double b;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	127 protected int numClasses;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	128
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	129 public IsoClasses(Raster raster, int numClasses) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	130
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	131 this.numClasses = numClasses;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	132
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	133 double min = Double.MAX_VALUE;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	134 double max = -Double.MAX_VALUE;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	135
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	136 double [] src = raster.raster;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	137
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	138 for (int i = 0; i < src.length; ++i) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	139 double x = src[i];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140 if (!Double.isNaN(x)) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	141 if (x < min) min = x;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142 if (x > max) max = x;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	143 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	144 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	145
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	146 /* f(min) = 0, f(max) = numClasses - 1
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	147
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148 I: 0 = m*min + b
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	149 II: numClasses - 1 = m*max + b
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	150
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	151 II - I:
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	152 numClasses - 1 = m*(max - min)
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	153
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	154 m = (numClasses - 1)/(max - min)
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	155 b = m*min
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	156 */
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	157
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	158 if (max == min) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159 m = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	160 b = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	161 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	162 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	163 m = (numClasses - 1)/(max - min);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	164 b = m*min;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	165 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	166 }
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 public int toIndex(double value) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169 return Double.isNaN(value)
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	170 ? -1
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	171 : Math.min(Math.max(0, (int)Math.round(m*value + b)), numClasses-1);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	172 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	173 } // class IsoClasses
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	175 protected double [] raster;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176 protected int 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 public Raster() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	179 this(new double[0], 0);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	180 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	181
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182 public Raster(double [] raster, int width) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	183 this.raster = raster;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	184 this.width = width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	185 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	186
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	187 public int getWidth() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	188 return width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	189 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	190
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	191 public int getHeight() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	192 return raster.length / width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	193 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	194
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	195 public int [] toIndexed(ValueToIndex valueToIndex) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	196 int [] dst = new int[raster.length];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	197 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	198 dst[i] = valueToIndex.toIndex(raster[i]);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	199 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	200 return dst;
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
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	203 public Raster create(Kernel kernel) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	204 double [] dst = new double[raster.length];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	205 Raster r = new Raster(dst, width);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	206 r.apply(kernel, continueBorder());
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	207 return r;
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
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	210 public void apply(Kernel kernel, Access access) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	211 int height = getHeight();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	212 for (int j = 0; j < height; ++j) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	213 int row = j*width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214 for (int i = 0; i < width; ++i) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	215 raster[row + i] = kernel.fold(access, i, j);
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 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	218 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	219
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	220 public Access continueBorder() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	221 return new Access() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	222 int height = getHeight()-1;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	223 public double get(int i, int j) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	224 if (i < 0) i = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	225 else if (i >= width) i = width-1;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	226 if (j < 0) j = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	227 else if (j > height) j = height;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	228 return raster[j*width + i];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	229 }
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 }
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 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8:

424

21fbd254db71 Added support for converting 2D rasters into polygons.