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 @ 812:41a7c49b5cb4

ISSUE 215 msg1170: Integrated additional Column into the ODV-Export to identify TimeSeries gnv-artifacts/trunk@896 c6561f87-3c4e-4783-a992-168aeb5c3f6f

author	Tim Englich <tim.englich@intevation.de>
date	Fri, 09 Apr 2010 13:00:38 +0000
parents	d766fe2d917a
children	05bf8534a35a

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 /**
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	4 * 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	5 * @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	6 */
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	7 public class Raster
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	8 {
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	9 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	10 * 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	11 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12 public static class Kernel {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	13
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	14 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	15 * The coefficients of the folding matrix.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	16 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17 public static final class Coeff {
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	18 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	19 * i delta index of this coefficient.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	20 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	21 public int di;
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 * i delta index of this coefficient.
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 int dj;
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 * The weight 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 double c;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	31 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	32 * Default constructor
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	33 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	34 public Coeff() {
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
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	37 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	38 * Constructor to set the delta position and the weight of
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	39 * the coefficient.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	40 * @param di The i delta index.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	41 * @param dj The j delta index.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	42 * @param c The weight of the index.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	43 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	44 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	45 this.di = di;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	46 this.dj = dj;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	47 this.c = c;
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
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	50 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	51 * Returns a string representation of this coefficient.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	52 * @return The string representation.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	53 */
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	54 @Override
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55 public String toString() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56 return String.valueOf(c);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	58
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	59 } // class Coeff
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	60
801 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 * The coefficients of the folding matrix.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	63 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	64 protected Coeff [] coeffs;
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	65 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	66 * The width of the kernel.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	67 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68 protected int width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	70 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	71 * Default constrcutor.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	72 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	73 public Kernel() {
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
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	76 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	77 * Constructor for a kernel with a given coefficient matrix
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	78 * with a given width.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	79 * @param coeffs The coefficients.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	80 * @param width The width of the kernel.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	81 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	82 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	83 this.coeffs = coeffs;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	84 this.width = width;
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
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	87 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	88 * Evaluates the 2D Gauss normal distribution a given x, y
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	89 * and a given sigma.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	90 * @param x The x location.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	91 * @param y The y location.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	92 * @param sigma The sigma of the distribution.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	93 * @return The value at point x, y.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	94 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	95 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	96 double s2 = sigma*sigma;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	97 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	98 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	99
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	100 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	101 * Creates a Gauss kernel with sigma = 1 and width 5.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	102 * @return The new kernel.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	103 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	104 public static Kernel createGauss() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	105 return createGauss(1.0d, 5);
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 a given sigma and width.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	110 * @param sigma The sigma value
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	111 * @param width The width of the kernel.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	112 * @return The new kernel.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	113 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	114 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	115 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	116 double sum = 0d;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	117 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	118 int y = -width/2 + j;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	119 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	120 int x = -width/2 + i;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	121 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	122 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	123 sum += c;
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 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	126 sum = 1.0/sum;
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 < coeffs.length; ++i) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	128 coeffs[i].c *= sum;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	129 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	130 return new Kernel(coeffs, width);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	131 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	132
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	133 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	134 * 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	135 * @param access The accessor to the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	136 * @param i The i raster index.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	137 * @param j The j raster index.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	138 * @return The folded value at index i, j.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	139 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140 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	141
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142 double s = 0.0d;
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 double unused = 0d;
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 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	147 Coeff coeff = coeffs[k];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148 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	149 if (Double.isNaN(v)) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	150 unused += coeff.c;
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 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	153 s += v * coeff.c;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	154 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	155 }
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 return s*(1.0d - unused);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	158 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	160 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	161 * Returns a string representation of this kernel.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	162 * @return The string representation.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	163 */
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	164 @Override
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	165 public String toString() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	166 StringBuilder sb = new StringBuilder();
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 int height = coeffs.length/width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	170 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	171 if (j > 0) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	172 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	173 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174 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	175 if (i > 0) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176 sb.append("\t");
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 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	179 }
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 return sb.toString();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	183 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	184 } // class Kernel
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	185
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	186 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	187 * Interface of allow special access patterns to the raster,
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	188 * especially at the border.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	189 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	190 public interface Access {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	191
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	192 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	193 * Returns the raster value at a given index point.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	194 * @param i The i index
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	195 * @param j The j index
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	196 * @return The value of the raster.
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 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	199
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	200 } // interface Access
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	201
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	202 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	203 * Interface to decouple the lookup which is needed
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	204 * to map double values to integers.
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 public interface ValueToIndex {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	207
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	208 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	209 * Returns the corresponding integer to a given value.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	210 * @param value The value to be queried.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	211 * @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	212 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	213 int toIndex(double value);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	215 } // interface ValueToIndex;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	216
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	217 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	218 * Brings raster to an integer representation with equal sized bins
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	219 * of integer vaules.
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 public static class IsoClasses
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	222 implements ValueToIndex
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	223 {
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	224 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	225 * Linear scaling factor.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	226 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	227 protected double m;
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	228 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	229 * Linear offset.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	230 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	231 protected double b;
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 * Number of classes.
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 int numClasses;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	236
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	237 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	238 * Constructor to build a lookup to transform
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	239 * an raster into an equals sized bin integer array.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	240 * @param raster The original raster
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	241 * @param numClasses The 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 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	244
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	245 this.numClasses = numClasses;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	246
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	247 double min = Double.MAX_VALUE;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	248 double max = -Double.MAX_VALUE;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	249
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	250 double [] src = raster.raster;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	251
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	252 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	253 double x = src[i];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	254 if (!Double.isNaN(x)) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	255 if (x < min) min = x;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	256 if (x > max) max = x;
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 }
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 /* 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	261
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	262 I: 0 = m*min + b
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	263 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	264
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	265 II - I:
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	266 numClasses - 1 = m*(max - min)
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 m = (numClasses - 1)/(max - min)
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	269 b = m*min
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	270 */
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	271
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	272 if (max == min) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	273 m = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	274 b = 0;
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 else {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	277 m = (numClasses - 1)/(max - min);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	278 b = m*min;
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 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	281
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	282 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	283 * Implements the double to integer lookup.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	284 * @param value The value to map.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	285 * @return The mapped value.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	286 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	287 public int toIndex(double value) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	288 return Double.isNaN(value)
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	289 ? -1
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	290 : 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	291 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	292 } // class IsoClasses
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	293
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	294 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	295 * Internal representation of the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	296 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	297 protected double [] raster;
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	298 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	299 * The width of the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	300 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	301 protected int width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	302
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	303 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	304 * Default constructor. Creates an 0x0 raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	305 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	306 public Raster() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	307 this(new double[0], 0);
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	308 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	309
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	310 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	311 * Creates a raster with given values and width.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	312 * @param raster
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	313 * @param width
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	314 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	315 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	316 this.raster = raster;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	317 this.width = width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	318 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	319
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	320 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	321 * The width of the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	322 * @return The width.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	323 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	324 public int getWidth() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	325 return 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 height of the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	330 * @return The height.
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 getHeight() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	333 return raster.length / 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 values of the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	338 * @return The values.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	339 */
495 6e8364e766fa Generate JTS geometries of "Horizontalschnitte" correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	340 public double [] getValues() {
6e8364e766fa Generate JTS geometries of "Horizontalschnitte" correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	341 return raster;
6e8364e766fa Generate JTS geometries of "Horizontalschnitte" correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 diff changeset	342 }
6e8364e766fa Generate JTS geometries of "Horizontalschnitte" correctly. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 424 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 * Creates an integer representation of this raster using a
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	346 * given value to index transformer.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	347 * @param valueToIndex The transformer to map double to integer values.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	348 * @return The new created integer representation.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	349 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	350 public int [] toIndexed(ValueToIndex valueToIndex) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	351 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	352 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	353 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	354 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	355 return dst;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	356 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	357
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	358 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	359 * 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	360 * Access to the raster is continuous at the border.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	361 * The original raster is not modified.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	362 * @param kernel The kernel to be applied.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	363 * @return The new raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	364 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	365 public Raster create(Kernel kernel) {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	366 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	367 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	368 r.apply(kernel, continueBorder());
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	369 return r;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	370 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	371
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	372 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	373 * 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	374 * Access to the raster is given.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	375 * @param kernel The kernel to be applied.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	376 * @param access The access method to the raster.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	377 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	378 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	379 int height = getHeight();
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	380 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	381 int row = j*width;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	382 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	383 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	384 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	385 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	386 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	387
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	388 /**
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	389 * Returns a continuous access to the raster. Values
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	390 * at the border are repeated til infinity if indices
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	391 * are accessed outside the defined area.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	392 * @return The accessor.
d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	393 */
424 21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	394 public Access continueBorder() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	395 return new Access() {
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	396 int height = getHeight()-1;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	397 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	398 if (i < 0) i = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	399 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	400 if (j < 0) j = 0;
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	401 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	402 return raster[j*width + i];
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	403 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	404 };
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	405 }
21fbd254db71 Added support for converting 2D rasters into polygons. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	406 }
801 d766fe2d917a More javadoc. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 798 diff changeset	407 // 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 @ 812:41a7c49b5cb4