Mercurial > dive4elements > river

annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/extreme/ExtremeCalculation.java @ 4385:10e277c2fe0f flys-2.9.4

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Some fixes to make the calculation work.
author Sascha L. Teichmann <teichmann@intevation.de>
date Sun, 04 Nov 2012 23:23:07 +0100
parents d095b4267772
children 73200f5462fa
rev   line source
3775
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
1 package de.intevation.flys.artifacts.model.extreme;
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
2
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
3 import org.apache.commons.logging.Log;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
4 import org.apache.commons.logging.LogFactory;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
5
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
6 import org.apache.commons.math.MathException;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
7
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
8 import org.apache.commons.math.optimization.fitting.CurveFitter;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
9
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
10 import org.apache.commons.math.optimization.general.LevenbergMarquardtOptimizer;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
11
3775
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
12 import de.intevation.flys.artifacts.access.ExtremeAccess;
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
13
4355
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
14 import de.intevation.flys.artifacts.math.Linear;
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
15 //import de.intevation.flys.artifacts.math.Utils;
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
16
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
17 import de.intevation.flys.artifacts.math.fitting.Function;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
18 import de.intevation.flys.artifacts.math.fitting.FunctionFactory;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
19
3775
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
20 import de.intevation.flys.artifacts.model.Calculation;
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
21 import de.intevation.flys.artifacts.model.CalculationResult;
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
22 import de.intevation.flys.artifacts.model.RangeWithValues;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
23 import de.intevation.flys.artifacts.model.RiverFactory;
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
24 import de.intevation.flys.artifacts.model.WQKms;
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
25 import de.intevation.flys.artifacts.model.WstValueTable;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
26 import de.intevation.flys.artifacts.model.WstValueTableFactory;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
27
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
28 import de.intevation.flys.model.River;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
29
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
30 import de.intevation.flys.utils.DoubleUtil;
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
31 import de.intevation.flys.utils.KMIndex;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
32
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
33 import gnu.trove.TDoubleArrayList;
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
34
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
35 import java.util.List;
3775
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
36
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
37 import java.awt.geom.Line2D;
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
38
4040
363445873737 Documentation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 3785
diff changeset
39 /** Calculate extrapolated W. */
3775
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
40 public class ExtremeCalculation
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
41 extends Calculation
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
42 {
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
43 private static final Log log =
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
44 LogFactory.getLog(ExtremeCalculation.class);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
45
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
46 protected String river;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
47 protected String function;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
48 protected double from;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
49 protected double to;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
50 protected double step;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
51 protected double percent;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
52 protected List<RangeWithValues> ranges;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
53
3775
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
54 public ExtremeCalculation() {
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
55 }
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
56
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
57 public ExtremeCalculation(ExtremeAccess access) {
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
58 String river = access.getRiver();
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
59 String function = access.getFunction();
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
60 Double from = access.getFrom();
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
61 Double to = access.getTo();
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
62 Double step = access.getStep();
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
63 Double percent = access.getPercent();
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
64 List<RangeWithValues> ranges = access.getRanges();
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
65
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
66 if (river == null) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
67 // TODO: i18n
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
68 addProblem("extreme.no.river");
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
69 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
70
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
71 if (function == null) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
72 // TODO: i18n
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
73 addProblem("extreme.no.function");
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
74 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
75
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
76 if (from == null) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
77 // TODO: i18n
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
78 addProblem("extreme.no.from");
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
79 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
80
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
81 if (to == null) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
82 // TODO: i18n
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
83 addProblem("extreme.no.to");
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
84 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
85
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
86 if (step == null) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
87 // TODO: i18n
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
88 addProblem("extreme.no.step");
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
89 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
90
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
91 if (percent == null) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
92 // TODO: i18n
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
93 addProblem("extreme.no.percent");
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
94 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
95
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
96 if (ranges == null) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
97 // TODO: i18n
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
98 addProblem("extreme.no.ranges");
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
99 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
100
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
101 if (!hasProblems()) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
102 this.river = river;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
103 this.function = function;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
104 this.from = Math.min(from, to);
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
105 this.to = Math.max(from, to);
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
106 this.step = Math.max(0.001d, Math.abs(step)/1000d);
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
107 this.percent = Math.max(0d, Math.min(100d, percent));
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
108 this.ranges = ranges;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
109 }
3775
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
110 }
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
111
4059
c63f0b4ac1b4 Documentation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 4040
diff changeset
112
c63f0b4ac1b4 Documentation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 4040
diff changeset
113 /** Calculate an extreme curve (extrapolate). */
3775
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
114 public CalculationResult calculate() {
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
115
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
116 WstValueTable wst = null;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
117
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
118 River river = RiverFactory.getRiver(this.river);
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
119 if (river == null) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
120 // TODO: i18n
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
121 addProblem("extreme.no.such.river", this.river);
3775
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
122 }
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
123 else {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
124 wst = WstValueTableFactory.getTable(river);
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
125 if (wst == null) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
126 // TODO: i18n
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
127 addProblem("extreme.no.wst.table");
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
128 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
129 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
130
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
131 Function function =
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
132 FunctionFactory.getInstance().getFunction(this.function);
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
133 if (function == null) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
134 // TODO: i18n
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
135 addProblem("extreme.no.such.function", this.function);
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
136 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
137
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
138 return hasProblems()
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
139 ? new CalculationResult(this)
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
140 : innerCalculate(wst, function);
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
141 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
142
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
143 protected String wqkmsName(int i) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
144 StringBuilder sb = new StringBuilder("W(");
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
145 boolean already = false;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
146 for (RangeWithValues r: ranges) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
147 double [] values = r.getValues();
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
148 if (i < values.length) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
149 if (already) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
150 sb.append(", ");
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
151 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
152 else {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
153 already = true;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
154 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
155 // TODO: i18n
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
156 sb.append(values[i]);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
157 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
158 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
159 return sb.append(')').toString();
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
160 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
161
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
162 protected WQKms [] allocWQKms() {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
163 int max = 0;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
164 for (RangeWithValues r: ranges) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
165 double [] values = r.getValues();
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
166 if (values.length > max) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
167 max = values.length;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
168 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
169 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
170 WQKms [] wqkms = new WQKms[max];
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
171 for (int i = 0; i < max; ++i) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
172 wqkms[i] = new WQKms(wqkmsName(i));
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
173 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
174 return wqkms;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
175 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
176
4059
c63f0b4ac1b4 Documentation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 4040
diff changeset
177
c63f0b4ac1b4 Documentation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 4040
diff changeset
178 /** Calculate an extreme curve (extrapolate). */
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
179 protected CalculationResult innerCalculate(
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
180 WstValueTable wst,
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
181 Function function
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
182 ) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
183 RangeWithValues range = null;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
184
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
185 double [] chiSqr = { 0d };
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
186
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
187 KMIndex<Curve> curves = new KMIndex<Curve>();
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
188 WQKms [] wqkms = allocWQKms();
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
189
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
190 boolean debug = log.isDebugEnabled();
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
191
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
192 from = DoubleUtil.round(from);
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
193 to = DoubleUtil.round(to);
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
194
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
195 for (double km = from; km <= to; km = DoubleUtil.round(km+step)) {
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
196
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
197 if (debug) {
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
198 log.debug("km: " + km);
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
199 }
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
200
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
201 if (range == null || !range.inside(km)) {
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
202 for (RangeWithValues r: ranges) {
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
203 if (r.inside(km)) {
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
204 range = r;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
205 break;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
206 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
207 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
208 // TODO: i18n
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
209 addProblem(km, "extreme.no.range.inner");
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
210 continue;
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
211 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
212
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
213 double [][] wqs = wst.interpolateTabulated(km);
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
214 if (wqs == null) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
215 // TODO: i18n
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
216 addProblem(km, "extreme.no.raw.data");
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
217 continue;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
218 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
219
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
220 // XXX: This should not be necessary for model data.
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
221 if (!DoubleUtil.isValid(wqs)) {
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
222 // TODO: i18n
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
223 addProblem(km, "extreme.invalid.data");
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
224 continue;
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
225 }
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
226
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
227 double [][] fitWQs = extractPointsToFit(wqs);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
228 if (fitWQs == null) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
229 // TODO: i18n
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
230 addProblem(km, "extreme.too.less.points");
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
231 continue;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
232 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
233
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
234 double [] coeffs = doFitting(function, fitWQs, chiSqr);
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
235 if (coeffs == null) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
236 // TODO: i18n
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
237 addProblem(km, "extreme.fitting.failed");
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
238 continue;
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
239 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
240
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
241 Curve curve = new Curve(
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
242 wqs[1], wqs[0],
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
243 function.getName(),
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
244 coeffs,
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
245 chiSqr[0]);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
246
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
247 curves.add(km, curve);
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
248
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
249 double [] values = range.getValues();
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
250
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
251 int V = Math.min(values.length, wqkms.length);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
252 for (int i = 0; i < V; ++i) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
253 double q = values[i];
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
254 double w = curve.value(q);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
255 if (Double.isNaN(w)) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
256 // TODO: i18n
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
257 addProblem(km, "extreme.evaluate.failed", values[i]);
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
258 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
259 else {
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
260 wqkms[i].add(w, q, km);
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
261 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
262 }
3785
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
263 }
a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from
Ingo Weinzierl <ingo.weinzierl@intevation.de>
parents: 3775
diff changeset
264
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
265 ExtremeResult result = new ExtremeResult(curves, wqkms);
3775
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
266 return new CalculationResult(result, this);
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
267 }
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
268
4355
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
269 protected double [] doFitting(
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
270 Function function,
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
271 double [][] wqs,
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
272 double [] chiSqr
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
273 ) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
274 LevenbergMarquardtOptimizer lmo = null;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
275
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
276 double [] coeffs = null;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
277
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
278 double [] ws = wqs[0];
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
279 double [] qs = wqs[1];
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
280
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
281 for (double tolerance = 1e-10; tolerance < 1e-3; tolerance *= 10d) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
282 lmo = new LevenbergMarquardtOptimizer();
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
283 lmo.setCostRelativeTolerance(tolerance);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
284 lmo.setOrthoTolerance(tolerance);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
285 lmo.setParRelativeTolerance(tolerance);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
286
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
287 CurveFitter cf = new CurveFitter(lmo);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
288
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
289 for (int i = 0; i < ws.length; ++i) {
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
290 cf.addObservedPoint(qs[i], ws[i]);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
291 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
292
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
293 try {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
294 coeffs = cf.fit(function, function.getInitialGuess());
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
295 break;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
296 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
297 catch (MathException me) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
298 if (log.isDebugEnabled()) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
299 log.debug("tolerance " + tolerance + " + failed.");
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
300 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
301 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
302 }
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
303 if (coeffs != null) {
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
304 chiSqr[0] = lmo.getChiSquare();
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
305 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
306 return coeffs;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
307 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
308
4355
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
309 protected double [][] extractPointsToFit(double [][] wqs) {
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
310
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
311 double [] ws = wqs[0];
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
312 double [] qs = wqs[1];
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
313
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
314 int N = Math.min(ws.length, qs.length);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
315
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
316 if (N < 2) {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
317 log.warn("Too less points for fitting");
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
318 return null;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
319 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
320
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
321 double q2 = qs[N-1];
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
322 double w2 = ws[N-1];
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
323 double q1 = qs[N-2];
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
324 double w1 = ws[N-2];
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
325
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
326 boolean ascending = w2 > w1;
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
327
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
328 TDoubleArrayList ows = new TDoubleArrayList();
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
329 TDoubleArrayList oqs = new TDoubleArrayList();
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
330
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
331 oqs.add(q2); oqs.add(q1);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
332 ows.add(w2); ows.add(w1);
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
333
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
334 int lastDir = -2;
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
335
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
336 for (int i = N-3; i >= 0; --i) {
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
337 double q = qs[i];
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
338 double w = ws[i];
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
339
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
340 if ((ascending && w > w1) || (!ascending && w < w1)) {
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
341 break;
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
342 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
343
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
344 int dir = Line2D.relativeCCW(q2, w2, q1, w1, q, w);
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
345 //int dir = Utils.relativeCCW(q2, w2, q1, w1, q, w);
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
346 if (lastDir == -2) {
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
347 lastDir = dir;
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
348 }
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
349 else if (lastDir != dir) {
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
350 break;
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
351 }
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
352
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
353 oqs.add(q);
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
354 ows.add(w);
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
355 w2 = w1;
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
356 q2 = q1;
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
357 w1 = w;
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
358 q1 = q;
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
359 }
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
360
4355
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
361 oqs.reverse();
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
362 ows.reverse();
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
363
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
364 boolean debug = log.isDebugEnabled();
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
365 if (debug) {
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
366 log.debug("from table: " + N);
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
367 log.debug("after trim: " + oqs.size());
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
368 }
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
369
4355
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
370 cutPercent(ows, oqs);
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
371
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
372 if (debug) {
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
373 log.debug("after percent cut: " + oqs.size());
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
374 }
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
375
4259
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
376 return new double [][] {
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
377 ows.toNativeArray(),
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
378 oqs.toNativeArray()
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
379 };
5cc9453456a7 First complete but untested version of the 'Auslagerung extremer Wasserspiegellagen' calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4059
diff changeset
380 }
4355
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
381
4385
10e277c2fe0f Some fixes to make the calculation work.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4355
diff changeset
382
4355
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
383 protected void cutPercent(TDoubleArrayList ws, TDoubleArrayList qs) {
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
384 int N = qs.size();
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
385 if (percent <= 0d || N == 0) {
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
386 return;
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
387 }
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
388
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
389 double minQ = qs.getQuick(0);
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
390 double maxQ = qs.getQuick(N-1);
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
391 double factor = Math.min(Math.max(0d, percent/100d), 1d);
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
392 double cutQ = Linear.weight(factor, minQ, maxQ);
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
393 int cutIndex = 0;
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
394 for (; cutIndex < N; ++cutIndex) {
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
395 double q = qs.getQuick(cutIndex);
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
396 if (minQ < maxQ) {
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
397 if (q > cutQ) {
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
398 break;
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
399 }
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
400 }
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
401 else {
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
402 if (q < cutQ) {
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
403 break;
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
404 }
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
405 }
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
406 }
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
407 ws.remove(0, cutIndex);
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
408 qs.remove(0, cutIndex);
d095b4267772 Added missing percentage cut off in extreme calculation.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 4294
diff changeset
409 }
3775
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
410 }
d52c4ca93ffb More infrastructure stuff for the "Auslagerung extremer Wasserspiegellagen".
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
411 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf-8 :
http://dive4elements.wald.intevation.org