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

annotate gnv-artifacts/src/main/java/de/intevation/gnv/raster/Raster.java @ 1129:ccfa07b88476

merged geo-backend

author	Thomas Arendsen Hein <thomas@intevation.de>
date	Fri, 28 Sep 2012 12:14:01 +0200
parents	f953c9a559d8
children

rev	line source
1115 f953c9a559d8 Added license file and license headers. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 836 diff changeset	1 /*
f953c9a559d8 Added license file and license headers. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 836 diff changeset	2 * Copyright (c) 2010 by Intevation GmbH
f953c9a559d8 Added license file and license headers. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 836 diff changeset	3 *
f953c9a559d8 Added license file and license headers. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 836 diff changeset	4 * This program is free software under the LGPL (>=v2.1)
f953c9a559d8 Added license file and license headers. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 836 diff changeset	5 * Read the file LGPL.txt coming with the software for details
f953c9a559d8 Added license file and license headers. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 836 diff changeset	6 * or visit http://www.gnu.org/licenses/ if it does not exist.
f953c9a559d8 Added license file and license headers. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 836 diff changeset	7 */
f953c9a559d8 Added license file and license headers. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 836 diff changeset	8
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	9 package de.intevation.gnv.raster;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	10
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11 /**
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	12 * A 2D double valued representation of a raster array.
780 c4156275c1e1 Bring @author javadoc tags in form '@author <a href="john.doe@example.com">John Doe</a>' Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 495 diff changeset	13 * @author <a href="mailto:sascha.teichmann@intevation.de">Sascha L. Teichmann</a>
424 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 Raster
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16 {
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	17 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	18 * A kernel can be applied to a raster to fold its values.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	19 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	20 public static class Kernel {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	21
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	22 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	23 * The coefficients of the folding matrix.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	24 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	25 public static final class Coeff {
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	26 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	27 * i delta index of this coefficient.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	28 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29 public int di;
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	30 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	31 * i delta index of this coefficient.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	32 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33 public int dj;
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	34 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	35 * The weight of this coefficient.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	36 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37 public double c;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	39 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	40 * Default constructor
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	41 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	42 public Coeff() {
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
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	45 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	46 * Constructor to set the delta position and the weight of
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	47 * the coefficient.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	48 * @param di The i delta index.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	49 * @param dj The j delta index.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	50 * @param c The weight of the index.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	51 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 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	53 this.di = di;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54 this.dj = dj;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55 this.c = c;
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
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	58 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	59 * Returns a string representation of this coefficient.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	60 * @return The string representation.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	61 */
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	62 @Override
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	63 public String toString() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	64 return String.valueOf(c);
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 } // class Coeff
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	69 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	70 * The coefficients of the folding matrix.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	71 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	72 protected Coeff [] coeffs;
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	73 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	74 * The width of the kernel.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	75 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	76 protected int width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	77
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	78 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	79 * Default constrcutor.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	80 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	81 public Kernel() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	82 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	83
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	84 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	85 * Constructor for a kernel with a given coefficient matrix
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	86 * with a given width.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	87 * @param coeffs The coefficients.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	88 * @param width The width of the kernel.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	89 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	90 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	91 this.coeffs = coeffs;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	92 this.width = 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
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	95 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	96 * Evaluates the 2D Gauss normal distribution a given x, y
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	97 * and a given sigma.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	98 * @param x The x location.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	99 * @param y The y location.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	100 * @param sigma The sigma of the distribution.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	101 * @return The value at point x, y.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	102 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	103 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	104 double s2 = sigma*sigma;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	105 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	106 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	107
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	108 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	109 * Creates a Gauss kernel with sigma = 1 and width 5.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	110 * @return The new kernel.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	111 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	112 public static Kernel createGauss() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	113 return createGauss(1.0d, 5);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	114 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	115
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	116 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	117 * Creates a Gauss kernel with a given sigma and width.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	118 * @param sigma The sigma value
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	119 * @param width The width of the kernel.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	120 * @return The new kernel.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	121 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	122 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	123 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	124 double sum = 0d;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	125 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	126 int y = -width/2 + j;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	127 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	128 int x = -width/2 + i;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	129 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	130 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	131 sum += c;
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 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	134 sum = 1.0/sum;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	135 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	136 coeffs[i].c *= sum;
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 return new Kernel(coeffs, width);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	139 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	141 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	142 * Evaluates this Kernel at the raster at raster index i, j.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	143 * @param access The accessor to the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	144 * @param i The i raster index.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	145 * @param j The j raster index.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	146 * @return The folded value at index i, j.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	147 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148 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	149
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	150 double s = 0.0d;
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 double unused = 0d;
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 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	155 Coeff coeff = coeffs[k];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	156 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	157 if (Double.isNaN(v)) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	158 unused += coeff.c;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	160 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	161 s += v * coeff.c;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	162 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	163 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	164
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	165 return s*(1.0d - unused);
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
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	168 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	169 * Returns a string representation of this kernel.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	170 * @return The string representation.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	171 */
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	172 @Override
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	173 public String toString() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174 StringBuilder sb = new StringBuilder();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	175
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176 int height = coeffs.length/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 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	179 if (j > 0) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	180 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	181 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182 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	183 if (i > 0) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	184 sb.append("\t");
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 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	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
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	190 return sb.toString();
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 } // class Kernel
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	193
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	194 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	195 * Interface of allow special access patterns to the raster,
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	196 * especially at the border.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	197 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	198 public interface Access {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	199
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	200 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	201 * Returns the raster value at a given index point.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	202 * @param i The i index
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	203 * @param j The j index
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	204 * @return The value of the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	205 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	206 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	207
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	208 } // interface Access
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	209
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	210 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	211 * Interface to decouple the lookup which is needed
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	212 * to map double values to integers.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	213 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214 public interface ValueToIndex {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	215
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	216 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	217 * Returns the corresponding integer to a given value.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	218 * @param value The value to be queried.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	219 * @return The integer value or -1 if no such value was found.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	220 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	221 int toIndex(double value);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	222
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	223 } // interface ValueToIndex;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	224
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	225 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	226 * Brings raster to an integer representation with equal sized bins
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	227 * of integer vaules.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	228 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	229 public static class IsoClasses
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	230 implements ValueToIndex
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	231 {
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	232 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	233 * Linear scaling factor.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	234 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	235 protected double m;
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	236 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	237 * Linear offset.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	238 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	239 protected double b;
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	240 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	241 * Number of classes.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	242 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	243 protected int numClasses;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	244
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	245 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	246 * Constructor to build a lookup to transform
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	247 * an raster into an equals sized bin integer array.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	248 * @param raster The original raster
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	249 * @param numClasses The number of classes.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	250 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	251 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	252
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	253 this.numClasses = numClasses;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	254
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	255 double min = Double.MAX_VALUE;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	256 double max = -Double.MAX_VALUE;
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 double [] src = raster.raster;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	259
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	260 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	261 double x = src[i];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	262 if (!Double.isNaN(x)) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	263 if (x < min) min = x;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	264 if (x > max) max = x;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	265 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	266 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	267
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	268 /* 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	269
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	270 I: 0 = m*min + b
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	271 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	272
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	273 II - I:
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	274 numClasses - 1 = m*(max - min)
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	275
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	276 m = (numClasses - 1)/(max - min)
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	277 b = m*min
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	278 */
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	279
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	280 if (max == min) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	281 m = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	282 b = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	283 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	284 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	285 m = (numClasses - 1)/(max - min);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	286 b = m*min;
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
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	290 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	291 * Implements the double to integer lookup.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	292 * @param value The value to map.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	293 * @return The mapped value.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	294 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	295 public int toIndex(double value) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	296 return Double.isNaN(value)
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	297 ? -1
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	298 : 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	299 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	300 } // class IsoClasses
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	301
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	302 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	303 * Internal representation of the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	304 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	305 protected double [] raster;
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	306 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	307 * The width of the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	308 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	309 protected int width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	310
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	311 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	312 * Default constructor. Creates an 0x0 raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	313 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	314 public Raster() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	315 this(new double[0], 0);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	316 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	317
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	318 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	319 * Creates a raster with given values and width.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	320 * @param raster
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	321 * @param width
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	322 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	323 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	324 this.raster = raster;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	325 this.width = width;
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
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	328 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	329 * The width of the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	330 * @return The width.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	331 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	332 public int getWidth() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	333 return width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	334 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	335
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	336 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	337 * The height of the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	338 * @return The height.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	339 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	340 public int getHeight() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	341 return raster.length / width;
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
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	344 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	345 * The values of the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	346 * @return The values.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	347 */
495 6e8364e766fa Generate JTS geometries of "Horizontalschnitte" correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	348 public double [] getValues() {
6e8364e766fa Generate JTS geometries of "Horizontalschnitte" correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	349 return raster;
6e8364e766fa Generate JTS geometries of "Horizontalschnitte" correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	350 }
6e8364e766fa Generate JTS geometries of "Horizontalschnitte" correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	351
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	352 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	353 * Creates an integer representation of this raster using a
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	354 * given value to index transformer.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	355 * @param valueToIndex The transformer to map double to integer values.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	356 * @return The new created integer representation.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	357 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	358 public int [] toIndexed(ValueToIndex valueToIndex) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	359 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	360 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	361 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	362 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	363 return dst;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	364 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	365
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	366 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	367 * Creates a new raster by applying a kernel to this raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	368 * Access to the raster is continuous at the border.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	369 * The original raster is not modified.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	370 * @param kernel The kernel to be applied.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	371 * @return The new raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	372 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	373 public Raster create(Kernel kernel) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	374 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	375 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	376 r.apply(kernel, continueBorder());
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	377 return r;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	378 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	379
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	380 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	381 * Creates a new raster by applying a kernel to this raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	382 * Access to the raster is given.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	383 * @param kernel The kernel to be applied.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	384 * @param access The access method to the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	385 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	386 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	387 int height = getHeight();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	388 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	389 int row = j*width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	390 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	391 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	392 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	393 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	394 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	395
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	396 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	397 * Returns a continuous access to the raster. Values
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	398 * at the border are repeated til infinity if indices
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	399 * are accessed outside the defined area.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	400 * @return The accessor.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	401 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	402 public Access continueBorder() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	403 return new Access() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	404 int height = getHeight()-1;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	405 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	406 if (i < 0) i = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	407 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	408 if (j < 0) j = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	409 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	410 return raster[j*width + i];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	411 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	412 };
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	413 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	414 }
836 05bf8534a35a Using unix line endings only. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 801 diff changeset	415 // 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/Raster.java @ 1129:ccfa07b88476