dive4elements/gnv-client: gnv-artifacts/src/main/java/de/intevation/gnv/math/Interpolation2D.java annotate

annotate gnv-artifacts/src/main/java/de/intevation/gnv/math/Interpolation2D.java @ 436:6642ab6c583c

Added vectorizer rings callback which generates polygon datasets suitable for polygon plots. gnv-artifacts/trunk@484 c6561f87-3c4e-4783-a992-168aeb5c3f6f

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Sat, 26 Dec 2009 15:32:08 +0000
parents	0eed5749fd63
children	f42ed4f10b79

rev	line source
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	1 package de.intevation.gnv.math;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	2
365 f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	3 import java.util.ArrayList;
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	4 import java.util.List;
365 f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	5 import java.util.HashMap;
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	6 import java.util.Collections;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	7
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	8 import com.vividsolutions.jts.geom.Coordinate;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	9 import com.vividsolutions.jts.geom.Envelope;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	10
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11 import com.vividsolutions.jts.index.quadtree.Quadtree;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12
365 f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	13 import org.apache.log4j.Logger;
f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	14
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15 /**
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16 * @author Sascha L. Teichmann
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17 */
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	18 public final class Interpolation2D
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	19 {
365 f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	20 private static Logger log = Logger.getLogger(Interpolation2D.class);
f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	21
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	22 public interface Consumer {
416 04a242c67fe6 Added support of gap detection in horizontalcrossprofile charts. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 366 diff changeset	23 void interpolated(Coordinate point, boolean success);
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	24 } // interface Consumer
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	25
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	26 private Interpolation2D() {
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	27 }
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28
431 422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	29 public static final double [] calculateBuffer(
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	30 List <? extends Point2d> points
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	31 ) {
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	32 HashMap<Integer, ArrayList<Point2d>> iMap =
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	33 new HashMap<Integer, ArrayList<Point2d>>();
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	34
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	35 HashMap<Integer, ArrayList<Point2d>> jMap =
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	36 new HashMap<Integer, ArrayList<Point2d>>();
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	37
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	38 for (int k = points.size()-1; k >= 0; --k) {
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	39 Point2d p = points.get(k);
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	40
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	41 ArrayList<Point2d> jList = jMap.get(p.j);
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	42 ArrayList<Point2d> iList = jMap.get(p.i);
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	43
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	44 if (jList == null) {
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	45 iMap.put(p.j, jList = new ArrayList<Point2d>());
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	46 }
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	47 jList.add(p);
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	48
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	49 if (iList == null) {
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	50 jMap.put(p.i, iList = new ArrayList<Point2d>());
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	51 }
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	52 iList.add(p);
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	53 }
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	54
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	55 double dxMax = -Double.MAX_VALUE;
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	56 double dyMax = -Double.MAX_VALUE;
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	57
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	58 for (ArrayList<Point2d> v: jMap.values()) {
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	59 Collections.sort(v, Point2d.Y_COMPARATOR);
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	60 for (int i = 1, L = v.size(); i < L; ++i) {
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	61 double dy = Math.abs(v.get(i).x - v.get(i-1).x);
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	62 if (dy > dyMax) {
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	63 dyMax = dy;
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	64 }
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	65 }
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	66 }
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	67
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	68 dyMax += 1e-5d;
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	69
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	70 for (ArrayList<Point2d> v: iMap.values()) {
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	71 Collections.sort(v, Point2d.X_COMPARATOR);
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	72 for (int i = 1, L = v.size(); i < L; ++i) {
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	73 double dx = Math.abs(v.get(i).x - v.get(i-1).x);
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	74 if (dx > dxMax) {
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	75 dxMax = dx;
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	76 }
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	77 }
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	78 }
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	79
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	80 dxMax += 1e-5d;
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	81
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	82 return new double [] { dxMax, dyMax };
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	83 }
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	84
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	85 public static void interpolate(
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	86 List<? extends Coordinate> path,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	87 List<? extends Point2d> points,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	88 double from,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	89 double to,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	90 int steps,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	91 Metrics metrics,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	92 Consumer consumer
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	93 ) {
434 0eed5749fd63 Implemented the raster interpolation for the 'Profilschnitt'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 431 diff changeset	94 boolean debug = log.isDebugEnabled();
0eed5749fd63 Implemented the raster interpolation for the 'Profilschnitt'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 431 diff changeset	95
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	96 int N = path.size();
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	97 int M = points.size();
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	98
434 0eed5749fd63 Implemented the raster interpolation for the 'Profilschnitt'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 431 diff changeset	99 if (debug) {
416 04a242c67fe6 Added support of gap detection in horizontalcrossprofile charts. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 366 diff changeset	100 log.debug("Size of path: " + N);
04a242c67fe6 Added support of gap detection in horizontalcrossprofile charts. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 366 diff changeset	101 log.debug("Size of points: " + M);
04a242c67fe6 Added support of gap detection in horizontalcrossprofile charts. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 366 diff changeset	102 }
365 f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	103
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	104 if (M < 1 \|\| N < 2) { // nothing to do
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	105 return;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	106 }
365 f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	107
431 422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	108 double [] buffer = calculateBuffer(points);
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	109 double dxMax = buffer[0];
422275fc9927 Refactored the XYColumn and Point2d code a bit to be more reusable in 3D. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 416 diff changeset	110 double dyMax = buffer[1];
365 f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	111
434 0eed5749fd63 Implemented the raster interpolation for the 'Profilschnitt'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 431 diff changeset	112 if (debug) {
416 04a242c67fe6 Added support of gap detection in horizontalcrossprofile charts. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 366 diff changeset	113 log.debug("buffer size: " + dxMax + " / " + dyMax);
04a242c67fe6 Added support of gap detection in horizontalcrossprofile charts. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 366 diff changeset	114 }
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	115
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	116 // put into spatial index to speed up finding neighbors.
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	117 Quadtree spatialIndex = new Quadtree();
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	118
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	119 for (int i = 0; i < M; ++i) {
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	120 Point2d p = points.get(i);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	121 spatialIndex.insert(p.envelope(), p);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	122 }
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	123
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	124 LinearToMap linearToMap = new LinearToMap(
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	125 path, from, to, metrics);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	126
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	127 double dP = (to - from)/steps;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	128
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	129 Coordinate center = new Coordinate();
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	130
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	131 Envelope queryBuffer = new Envelope();
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	132
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	133 Point2d [] neighbors = new Point2d[4];
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	134
365 f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	135 int missedInterpolations = 0;
f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	136 int interpolations = 0;
f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	137
f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	138 for (double p = from; p <= to; p += dP) {
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	139 if (!linearToMap.locate(p, center)) {
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140 continue;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	141 }
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142 queryBuffer.init(
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	143 center.x - dxMax, center.x + dxMax,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	144 center.y - dyMax, center.y + dyMax);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	145
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	146 List potential = spatialIndex.query(queryBuffer);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	147
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148 L1Comparator invL1 = new L1Comparator(center);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	149 Collections.sort(potential, invL1);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	150
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	151 neighbors[0] = neighbors[1] =
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	152 neighbors[2] = neighbors[3] = null;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	153
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	154 /* bit code of neighbors
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	155 0---1
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	156 \| x \|
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	157 2---3
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	158 */
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	160 int mask = 0;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	161 // reversed order is essential here!
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	162 for (int i = potential.size()-1; i >= 0; --i) {
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	163 Point2d n = (Point2d)potential.get(i);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	164 int code = n.x > center.x ? 1 : 0;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	165 if (n.y > center.y) code \|= 2;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	166 neighbors[code] = n;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	167 mask \|= 1 << code;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	168 }
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	170 int numNeighbors = Integer.bitCount(mask);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	171
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	172 // only interpolate if we have all four neighbors
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	173 // and we do not have any gaps.
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174 if (numNeighbors == 4
366 086e3af38b96 Fixed index error in i-gap detection in 2D interpolation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 365 diff changeset	175 && !neighbors[0].hasIGap(neighbors[1])
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176 && !neighbors[1].hasJGap(neighbors[3])
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	177 && !neighbors[3].hasIGap(neighbors[2])
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	178 && !neighbors[2].hasJGap(neighbors[0])
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	179 ) {
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	180 double z1 = interpolate(
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	181 neighbors[0].x, neighbors[0].z,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182 neighbors[1].x, neighbors[1].z,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	183 center.x);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	184 double z2 = interpolate(
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	185 neighbors[2].x, neighbors[2].z,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	186 neighbors[3].x, neighbors[3].z,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	187 center.x);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	188 double y1 = interpolate(
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	189 neighbors[0].x, neighbors[0].y,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	190 neighbors[1].x, neighbors[1].y,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	191 center.x);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	192 double y2 = interpolate(
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	193 neighbors[2].x, neighbors[2].y,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	194 neighbors[3].x, neighbors[3].y,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	195 center.x);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	196 center.z = interpolate(
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	197 y1, z1,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	198 y2, z2,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	199 center.y);
416 04a242c67fe6 Added support of gap detection in horizontalcrossprofile charts. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 366 diff changeset	200 consumer.interpolated(center, true);
365 f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	201 ++interpolations;
f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	202 }
f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	203 else {
416 04a242c67fe6 Added support of gap detection in horizontalcrossprofile charts. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 366 diff changeset	204 consumer.interpolated(center, false);
365 f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	205 ++missedInterpolations;
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	206 }
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	207 }
365 f66088a43ecc Added horizontal crossprofile charts to chart pallet. Fixed some bugs before interpolation. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 361 diff changeset	208
434 0eed5749fd63 Implemented the raster interpolation for the 'Profilschnitt'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 431 diff changeset	209 if (debug) {
416 04a242c67fe6 Added support of gap detection in horizontalcrossprofile charts. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 366 diff changeset	210 log.debug("interpolations: " +
04a242c67fe6 Added support of gap detection in horizontalcrossprofile charts. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 366 diff changeset	211 interpolations + " / " + missedInterpolations);
04a242c67fe6 Added support of gap detection in horizontalcrossprofile charts. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 366 diff changeset	212 }
361 aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	213 }
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	215 public static final double interpolate(
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	216 double x1, double y1,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	217 double x2, double y2,
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	218 double x
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	219 ) {
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	220 if (x2 == x1) {
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	221 return (y1 + y2)*0.5d;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	222 }
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	223 double m = (y2-y1)/(x2-x1);
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	224 double b = y1 - m*x1;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	225 return m*x + b;
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	226 }
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	227 }
aec85d00d82c Added code to do 2D interpolations along a digitied track on the map. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	228 // 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/math/Interpolation2D.java @ 436:6642ab6c583c