Mercurial > dive4elements > river
annotate artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/collision/GaugeDischargeValuesFinder.java @ 9157:f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
author | mschaefer |
---|---|
date | Tue, 19 Jun 2018 14:19:32 +0200 |
parents | |
children |
rev | line source |
---|---|
9157
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
1 /** Copyright (C) 2017 by Bundesanstalt für Gewässerkunde |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
2 * Software engineering by |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
3 * Björnsen Beratende Ingenieure GmbH |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
4 * Dr. Schumacher Ingenieurbüro für Wasser und Umwelt |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
5 * |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
6 * This file is Free Software under the GNU AGPL (>=v3) |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
7 * and comes with ABSOLUTELY NO WARRANTY! Check out the |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
8 * documentation coming with Dive4Elements River for details. |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
9 */ |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
10 package org.dive4elements.river.artifacts.sinfo.collision; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
11 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
12 import org.apache.commons.lang.math.DoubleRange; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
13 import org.apache.commons.math.FunctionEvaluationException; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
14 import org.apache.commons.math.analysis.UnivariateRealFunction; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
15 import org.apache.commons.math.analysis.interpolation.LinearInterpolator; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
16 import org.dive4elements.river.artifacts.model.Calculation; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
17 import org.dive4elements.river.model.DischargeTable; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
18 import org.dive4elements.river.model.DischargeTableValue; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
19 import org.dive4elements.river.model.Gauge; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
20 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
21 import gnu.trove.TDoubleArrayList; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
22 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
23 /** |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
24 * Loading and search/interpolation of a gauge's discharge table |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
25 * |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
26 * @author Matthias Schäfer |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
27 * |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
28 */ |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
29 public final class GaugeDischargeValuesFinder { |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
30 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
31 /***** FIELDS *****/ |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
32 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
33 // private static Logger log = Logger.getLogger(GaugeDischargeValuesFinder.class); |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
34 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
35 private final Gauge gauge; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
36 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
37 private final Calculation problems; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
38 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
39 private final UnivariateRealFunction wInterpolator; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
40 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
41 private final DoubleRange wRange; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
42 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
43 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
44 /***** CONSTRUCTORS *****/ |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
45 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
46 private GaugeDischargeValuesFinder(final Gauge gauge, final Calculation problems, final DischargeTable dischargeTable) { |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
47 this.gauge = gauge; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
48 this.problems = problems; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
49 final TDoubleArrayList ws = new TDoubleArrayList(); |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
50 final TDoubleArrayList qs = new TDoubleArrayList(); |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
51 for (final DischargeTableValue v : DischargeTable.getValuesSortedByW(dischargeTable)) { |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
52 ws.add(v.getW().doubleValue()); |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
53 qs.add(v.getQ().doubleValue()); |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
54 } |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
55 if (!ws.isEmpty()) |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
56 this.wRange = new DoubleRange(ws.get(0), ws.get(ws.size() - 1)); |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
57 else |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
58 this.wRange = null; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
59 this.wInterpolator = new LinearInterpolator().interpolate(ws.toNativeArray(), qs.toNativeArray()); |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
60 } |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
61 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
62 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
63 /***** METHODS *****/ |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
64 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
65 /** |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
66 * Loads the the main discharge table of a gauge (GAUGE.at) |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
67 * |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
68 * @return The discharge table values finder of the gauge, or null |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
69 */ |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
70 public static GaugeDischargeValuesFinder loadValues(final Gauge gauge, final Calculation problems) { |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
71 final DischargeTable table = DischargeTable.getGaugeMainDischargeTable(gauge); |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
72 if ((table == null) || (table.getDischargeTableValues().size() == 0)) |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
73 return null; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
74 else { |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
75 return new GaugeDischargeValuesFinder(gauge, problems, table); |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
76 } |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
77 } |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
78 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
79 /** |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
80 * If this provider may return valid data at all. |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
81 */ |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
82 public boolean isValid() { |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
83 return (this.wInterpolator != null); |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
84 } |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
85 |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
86 /** |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
87 * Discharge for a W |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
88 * |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
89 * @return Q, or NegInf for w less than all, or PosInf for w greater then all, or NaN in case of exception |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
90 */ |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
91 public double getDischarge(final double w) { |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
92 try { |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
93 if (this.wInterpolator == null) |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
94 return Double.NaN; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
95 else if (w < this.wRange.getMinimumDouble()) |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
96 return Double.NEGATIVE_INFINITY; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
97 else if (w > this.wRange.getMaximumDouble()) |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
98 return Double.POSITIVE_INFINITY; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
99 else |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
100 return this.wInterpolator.value(w); |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
101 } |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
102 catch (@SuppressWarnings("unused") final FunctionEvaluationException e) { |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
103 // ignore exception because this can/will happen regularly |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
104 return Double.NaN; |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
105 } |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
106 } |
f9bb5d0a6ff3
Added the S-Info collision calculation and chart output
mschaefer
parents:
diff
changeset
|
107 } |