Mercurial > dive4elements > river

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/uinfo/salix/SalixLineCalculation.java @ 9504:76c0665888a3

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
No rounding during calculation (Meilenstein-2 2.4.2 and 2.9), delta-w-cm as double (for historical scenario)
author mschaefer
date Fri, 28 Sep 2018 10:13:09 +0200
parents 853f2dafc16e
children 429b62373633
rev   line source
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
1 /** Copyright (C) 2017 by Bundesanstalt für Gewässerkunde
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
2 * Software engineering by
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
3 * Björnsen Beratende Ingenieure GmbH
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
4 * Dr. Schumacher Ingenieurbüro für Wasser und Umwelt
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
5 *
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
6 * This file is Free Software under the GNU AGPL (>=v3)
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
7 * and comes with ABSOLUTELY NO WARRANTY! Check out the
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
8 * documentation coming with Dive4Elements River for details.
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
9 */
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
10 package org.dive4elements.river.artifacts.uinfo.salix;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
11
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
12 import java.util.ArrayList;
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
13 import java.util.Collection;
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
14 import java.util.List;
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
15 import java.util.NavigableMap;
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
16 import java.util.TreeMap;
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
17
9321
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
18 import org.apache.commons.lang.math.DoubleRange;
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
19 import org.dive4elements.artifacts.CallContext;
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
20 import org.dive4elements.river.artifacts.access.RangeAccess;
9361
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
21 import org.dive4elements.river.artifacts.common.GeneralResultType;
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
22 import org.dive4elements.river.artifacts.model.Calculation;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
23 import org.dive4elements.river.artifacts.model.CalculationResult;
9499
853f2dafc16e VegetationZones in CrossSectionsDiagram
gernotbelger
parents: 9421
diff changeset
24 import org.dive4elements.river.artifacts.model.river.RiverInfoProvider;
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
25 import org.dive4elements.river.artifacts.resources.Resources;
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
26 import org.dive4elements.river.artifacts.sinfo.tkhstate.BedHeightsFinder;
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
27 import org.dive4elements.river.artifacts.sinfo.tkhstate.BedHeightsUtils;
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
28 import org.dive4elements.river.artifacts.sinfo.tkhstate.DefaultBedHeights;
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
29 import org.dive4elements.river.artifacts.sinfo.util.CalculationUtils;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
30 import org.dive4elements.river.artifacts.sinfo.util.RiverInfo;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
31 import org.dive4elements.river.artifacts.uinfo.UINFOArtifact;
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
32 import org.dive4elements.river.artifacts.uinfo.salix.SalixLineAccess.ScenarioType;
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
33 import org.dive4elements.river.model.BedHeight;
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
34 import org.dive4elements.river.model.River;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
35
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
36 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
37 * Calculation of a iota (former salix) longitudinal section, optionally with a delta scenario
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
38 *
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
39 * @author Domenico Nardi Tironi
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
40 * @author Matthias Schäfer
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
41 *
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
42 */
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
43 final class SalixLineCalculation {
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
44
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
45 private final CallContext context;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
46
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
47 private Calculation problems;
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
48
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
49 public SalixLineCalculation(final CallContext context) {
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
50 this.context = context;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
51 }
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
52
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
53 /**
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
54 * Calculates the iota longitudinal section and delta scenario of a uinfo artifact
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
55 */
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
56 public CalculationResult calculate(final UINFOArtifact uinfo) {
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
57 this.problems = new Calculation();
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
58
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
59 final String calcModeLabel = Resources.getMsg(this.context.getMeta(), uinfo.getCalculationMode().name());
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
60 final String user = CalculationUtils.findArtifactUser(this.context, uinfo);
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
61
9070
611a523fc42f VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents: 9069
diff changeset
62 final SalixLineAccess accessSalix = new SalixLineAccess(uinfo);
611a523fc42f VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents: 9069
diff changeset
63
611a523fc42f VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents: 9069
diff changeset
64 final River river = accessSalix.getRiver();
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
65 final RiverInfo riverInfo = new RiverInfo(river);
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
66
9321
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
67 final DoubleRange range = accessSalix.getRange();
9068
a5cf8d7bff3c access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents: 8996
diff changeset
68
9321
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
69 final ScenarioType scenarioType = accessSalix.getScenario();
9070
611a523fc42f VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents: 9069
diff changeset
70
9397
6e7094368e97 Added gauge name column for uinfo iota
mschaefer
parents: 9394
diff changeset
71 final RiverInfoProvider riverInfoProvider = RiverInfoProvider.forRange(this.context, river, range, true);
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
72
9321
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
73 final SalixLineCalculationResults results = new SalixLineCalculationResults(calcModeLabel, user, riverInfo, range);
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
74
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
75 final SalixLineCalculator calculator = new SalixLineCalculator(riverInfoProvider);
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
76 final NavigableMap<Double, List<Double>> rangeScenarios = buildRangeScenarios(accessSalix);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
77
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
78 calculator.execute(this.problems, uinfo, rangeScenarios, scenarioType, buildScenarioLabels(accessSalix), buildPartialRangeString(accessSalix),
9361
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
79 buildAdditionalString(accessSalix), results);
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
80
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
81 return new CalculationResult(results, this.problems);
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
82 }
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
83
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
84 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
85 * Builds a map of delta-Ws by from-km for the selected scenario
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
86 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
87 private NavigableMap<Double, List<Double>> buildRangeScenarios(final SalixLineAccess access) {
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
88 final NavigableMap<Double, List<Double>> rangeScenarios = new TreeMap<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
89 if (access.getScenario() == ScenarioType.REGIONAL)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
90 fillRangeScenarios(rangeScenarios, access, access.getFromPart().doubleValue(), access.getToPart().doubleValue(),
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
91 access.getRegionalScenarioIntegers());
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
92 else if (access.getScenario() == ScenarioType.SUPRAREGIONAL)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
93 fillRangeScenarios(rangeScenarios, access.getSupraRegionalString());
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
94 else if (access.getScenario() == ScenarioType.HISTORICAL)
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
95 fillRangeScenarios(rangeScenarios, access, access.getFromPart().doubleValue(), access.getToPart().doubleValue(), access.getBedHeightId());
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
96 else
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
97 fillRangeScenarios(rangeScenarios, access);
9321
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
98
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
99 return rangeScenarios;
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
100 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
101
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
102 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
103 * Fills a map of delta-Ws with only one 0-delta for the whole calc range (no scenario)
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
104 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
105 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final RangeAccess calcRange) {
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
106 final List<Double> nulls = new ArrayList<>();
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
107 nulls.add(null);
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
108 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
109 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
110
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
111 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
112 * Fills a map of delta-Ws by km-range from the regional scenario input data
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
113 */
9499
853f2dafc16e VegetationZones in CrossSectionsDiagram
gernotbelger
parents: 9421
diff changeset
114 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final RangeAccess calcRange, final double partFrom,
853f2dafc16e VegetationZones in CrossSectionsDiagram
gernotbelger
parents: 9421
diff changeset
115 final double partTo, final int[] deltaWs) {
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
116 final List<Double> nulls = new ArrayList<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
117 final List<Double> dwsm = new ArrayList<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
118 for (int i = 0; i <= deltaWs.length - 1; i++) {
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
119 nulls.add(null);
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
120 dwsm.add(deltaWs[i] / 100.0);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
121 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
122 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
123 rangeScenarios.put(Double.valueOf(partFrom - 0.0001), dwsm);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
124 rangeScenarios.put(Double.valueOf(partTo + 0.0001), nulls);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
125 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
126
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
127 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
128 * Fills a map of delta-Ws by km-range from the supraregional scenario input data
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
129 * (the zones input by the user cover the calc range completely)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
130 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
131 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final String zones) {
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
132 final List<SalixZone> parts = SalixZone.parse(zones);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
133 for (final SalixZone part : parts) {
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
134 final List<Double> dwsm = new ArrayList<>();
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
135 if (part.getDwsplValue() == 0)
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
136 dwsm.add(null);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
137 else
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
138 dwsm.add(part.getDwsplValue() / 100.0);
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
139 rangeScenarios.put(Double.valueOf(part.getFromKm().doubleValue() - 0.0001), dwsm);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
140 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
141 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
142
9316
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
143 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
144 * Fetches historical and reference bed heights and fills a map of delta-MSHs for all fetched stations in the calc range
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
145 */
9499
853f2dafc16e VegetationZones in CrossSectionsDiagram
gernotbelger
parents: 9421
diff changeset
146 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final RangeAccess calcRange, final double partFrom,
853f2dafc16e VegetationZones in CrossSectionsDiagram
gernotbelger
parents: 9421
diff changeset
147 final double partTo, final int historicalBedHeightId) {
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
148
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
149 // Find relevant default bed-heights
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
150 final River river = calcRange.getRiver();
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
151 final List<BedHeight> defaultBedHeights = new DefaultBedHeights(river).getBedHeights(this.problems);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
152 if (defaultBedHeights.isEmpty())
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
153 return;
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
154 final DoubleRange scenarioRange = new DoubleRange(partFrom, partTo);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
155 final Collection<BedHeightsFinder> allFinders = BedHeightsFinder.createTkhBedHeights(this.problems, scenarioRange, defaultBedHeights);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
156 final Collection<BedHeightsFinder> currentFinders = new ArrayList<>(allFinders);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
157
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
158 // Add historical bed-heights
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
159 final BedHeightsFinder historicalFinder = BedHeightsFinder.forId(this.problems, historicalBedHeightId, scenarioRange);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
160 allFinders.add(historicalFinder);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
161 final Collection<Double> stations = BedHeightsUtils.extractStationCollection(allFinders, true);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
162 final List<Double> nulls = new ArrayList<>();
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
163 nulls.add(null);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
164 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
165 for (final Double station : stations) {
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
166 rangeScenarios.put(station, new ArrayList<Double>());
9421
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
167 final double delta = bedHeightDifference(station, currentFinders, historicalFinder);
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
168 if (Double.isNaN(delta)) {
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
169 rangeScenarios.get(station).add(null);
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
170 if (!this.problems.hasProblems()) {
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
171 final String msg = Resources.getMsg(this.context.getMeta(), "uinfo_salix_calc.warning.missing_bedheights");
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
172 this.problems.addProblem(station, msg);
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
173 }
9499
853f2dafc16e VegetationZones in CrossSectionsDiagram
gernotbelger
parents: 9421
diff changeset
174 } else
9421
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
175 rangeScenarios.get(station).add(Double.valueOf(delta));
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
176 }
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
177 rangeScenarios.put(Double.valueOf(partTo + 0.0001), nulls);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
178 }
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
179
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
180 /**
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
181 * Gets the difference of a historical bed height against a current one for a station
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
182 */
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
183 private double bedHeightDifference(final double station, final Collection<BedHeightsFinder> currentFinders, final BedHeightsFinder historicalFinder) {
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
184 double currentMSH = Double.NaN;
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
185 for (final BedHeightsFinder bhf : currentFinders) {
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
186 currentMSH = bhf.getMeanBedHeight(station);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
187 if (!Double.isNaN(currentMSH))
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
188 break;
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
189 }
9421
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
190 if (Double.isNaN(currentMSH))
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
191 return Double.NaN;
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
192 final double historicalMSH = historicalFinder.getMeanBedHeight(station);
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
193 if (Double.isNaN(historicalMSH))
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
194 return Double.NaN;
9504
76c0665888a3 No rounding during calculation (Meilenstein-2 2.4.2 and 2.9), delta-w-cm as double (for historical scenario)
mschaefer
parents: 9499
diff changeset
195 return (historicalMSH - currentMSH);
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
196 }
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
197
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
198 /**
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
199 * Builds the list of delta-w labels for the scenario type
9316
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
200 */
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
201 private String[] buildScenarioLabels(final SalixLineAccess access) {
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
202 final List<String> labels = new ArrayList<>();
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
203 if (access.getScenario() == ScenarioType.REGIONAL) {
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
204 final int[] deltaws = access.getRegionalScenarioIntegers();
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
205 for (int i = 0; i <= deltaws.length - 1; i++)
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
206 if (deltaws[i] != 0)
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
207 labels.add(Integer.toString(deltaws[i]) + " cm");
9391
2da486c7c05f salix historical sounding workflow changed
gernotbelger
parents: 9361
diff changeset
208 } else if (access.getScenario() == ScenarioType.SUPRAREGIONAL)
9361
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
209 labels.add(Resources.getMsg(this.context.getMeta(), "uinfo_salix_scenario_supraregional"));
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
210 else if (access.getScenario() == ScenarioType.HISTORICAL)
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
211 labels.add(Resources.getMsg(this.context.getMeta(), "uinfo_salix_scenario_historical"));
9316
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
212 return labels.toArray(new String[labels.size()]);
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
213 }
9361
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
214
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
215 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
216 * Builds the km range string for the scenario type
9361
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
217 */
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
218 private String buildPartialRangeString(final SalixLineAccess access) {
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
219 if ((access.getScenario() == ScenarioType.REGIONAL) || (access.getScenario() == ScenarioType.HISTORICAL)) {
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
220 return String.format("%s - %s", GeneralResultType.station.exportValue(this.context, access.getFromPart().doubleValue()),
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
221 GeneralResultType.station.exportValue(this.context, access.getToPart().doubleValue()));
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
222 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
223 if (access.getScenario() == ScenarioType.SUPRAREGIONAL) {
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
224 String ranges = "";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
225 String sep = "";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
226 final List<SalixZone> parts = SalixZone.parse(access.getSupraRegionalString());
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
227 for (final SalixZone part : parts) {
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
228 if (part.getDwsplValue() != 0) {
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
229 ranges = ranges + sep + String.format("%s - %s", GeneralResultType.station.exportValue(this.context, part.getFromKm().doubleValue()),
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
230 GeneralResultType.station.exportValue(this.context, part.getToKm().doubleValue()));
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
231 sep = ", ";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
232 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
233 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
234 return ranges;
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
235 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
236 return "";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
237 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
238
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
239 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
240 * Builds the delta w or time string for the scenario type
9361
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
241 */
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
242 private String buildAdditionalString(final SalixLineAccess access) {
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
243 if (access.getScenario() == ScenarioType.REGIONAL) {
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
244 String deltas = "";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
245 String sep = "";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
246 for (final int d : access.getRegionalScenarioIntegers()) {
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
247 deltas = deltas + sep + Integer.toString(d);
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
248 sep = ", ";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
249 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
250 return deltas;
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
251 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
252 if (access.getScenario() == ScenarioType.HISTORICAL) {
9391
2da486c7c05f salix historical sounding workflow changed
gernotbelger
parents: 9361
diff changeset
253 return String.valueOf(access.getYearEpoch());
9361
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
254 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
255 if (access.getScenario() == ScenarioType.SUPRAREGIONAL) {
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
256 String deltas = "";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
257 String sep = "";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
258 final List<SalixZone> parts = SalixZone.parse(access.getSupraRegionalString());
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
259 for (final SalixZone part : parts) {
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
260 if (part.getDwsplValue() != 0) {
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
261 deltas = deltas + sep + Integer.toString(part.getDwsplValue());
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
262 sep = ", ";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
263 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
264 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
265 return deltas;
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
266 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
267 return "";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
268 }
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
269 }
http://dive4elements.wald.intevation.org