Mercurial > dive4elements > river

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/uinfo/salix/SalixLineCalculation.java @ 9394:439699ff9b2d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added U-Info iota (prev. salix) calculation for historical scenario
author mschaefer
date Fri, 10 Aug 2018 17:31:46 +0200
parents 2da486c7c05f
children 6e7094368e97
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;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
24 import org.dive4elements.river.artifacts.resources.Resources;
9295
385b52ccde23 Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents: 9249
diff changeset
25 import org.dive4elements.river.artifacts.sinfo.common.RiverInfoProvider;
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;
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
35 import org.dive4elements.river.utils.Formatter;
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
36
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
37 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
38 * 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
39 *
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
40 * @author Domenico Nardi Tironi
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
41 * @author Matthias Schäfer
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
42 *
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
43 */
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
44 final class SalixLineCalculation {
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
45
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
46 private final CallContext context;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
47
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
48 private Calculation problems;
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
49
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
50 public SalixLineCalculation(final CallContext context) {
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
51 this.context = context;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
52 }
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
53
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
54 /**
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
55 * 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
56 */
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
57 public CalculationResult calculate(final UINFOArtifact uinfo) {
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
58 this.problems = new Calculation();
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
59
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
60 final String calcModeLabel = Resources.getMsg(this.context.getMeta(), uinfo.getCalculationMode().name());
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
61 final String user = CalculationUtils.findArtifactUser(this.context, uinfo);
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
62
9070
611a523fc42f VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents: 9069
diff changeset
63 final SalixLineAccess accessSalix = new SalixLineAccess(uinfo);
611a523fc42f VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents: 9069
diff changeset
64
611a523fc42f VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents: 9069
diff changeset
65 final River river = accessSalix.getRiver();
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
66 final RiverInfo riverInfo = new RiverInfo(river);
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
67
9321
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
68 final DoubleRange range = accessSalix.getRange();
9068
a5cf8d7bff3c access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents: 8996
diff changeset
69
9321
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
70 final ScenarioType scenarioType = accessSalix.getScenario();
9070
611a523fc42f VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents: 9069
diff changeset
71
9321
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
72 final RiverInfoProvider riverInfoProvider = RiverInfoProvider.forRange(this.context, river, range);
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
73
9321
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
74 final SalixLineCalculationResults results = new SalixLineCalculationResults(calcModeLabel, user, riverInfo, range);
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
75
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
76 final SalixLineCalculator calculator = new SalixLineCalculator(riverInfoProvider);
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
77 final NavigableMap<Double, List<Double>> rangeScenarios = buildRangeScenarios(accessSalix);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
78
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
79 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
80 buildAdditionalString(accessSalix), results);
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
81
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
82 return new CalculationResult(results, this.problems);
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
83 }
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
84
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
85 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
86 * 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
87 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
88 private NavigableMap<Double, List<Double>> buildRangeScenarios(final SalixLineAccess access) {
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
89 final NavigableMap<Double, List<Double>> rangeScenarios = new TreeMap<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
90 if (access.getScenario() == ScenarioType.REGIONAL)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
91 fillRangeScenarios(rangeScenarios, access, access.getFromPart().doubleValue(), access.getToPart().doubleValue(),
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
92 access.getRegionalScenarioIntegers());
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
93 else if (access.getScenario() == ScenarioType.SUPRAREGIONAL)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
94 fillRangeScenarios(rangeScenarios, access.getSupraRegionalString());
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
95 else if (access.getScenario() == ScenarioType.HISTORICAL)
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
96 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
97 else
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
98 fillRangeScenarios(rangeScenarios, access);
9321
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
99
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
100 return rangeScenarios;
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
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
103 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
104 * 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
105 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
106 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
107 final List<Double> nulls = new ArrayList<>();
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
108 nulls.add(null);
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
109 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls);
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
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
112 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
113 * 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
114 */
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
115 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final RangeAccess calcRange,
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
116 final double partFrom, final double partTo, final int[] deltaWs) {
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
117 final List<Double> nulls = new ArrayList<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
118 final List<Double> dwsm = new ArrayList<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
119 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
120 nulls.add(null);
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
121 dwsm.add(deltaWs[i] / 100.0);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
122 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
123 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
124 rangeScenarios.put(Double.valueOf(partFrom - 0.0001), dwsm);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
125 rangeScenarios.put(Double.valueOf(partTo + 0.0001), nulls);
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
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
128 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
129 * 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
130 * (the zones input by the user cover the calc range completely)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
131 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
132 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
133 final List<SalixZone> parts = SalixZone.parse(zones);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
134 for (final SalixZone part : parts) {
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
135 final List<Double> dwsm = new ArrayList<>();
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
136 if (part.getDwsplValue() == 0)
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
137 dwsm.add(null);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
138 else
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
139 dwsm.add(part.getDwsplValue() / 100.0);
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
140 rangeScenarios.put(Double.valueOf(part.getFromKm().doubleValue() - 0.0001), dwsm);
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 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
143
9316
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
144 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
145 * 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
146 */
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
147 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final RangeAccess calcRange,
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
148 final double partFrom, final double partTo, final int historicalBedHeightId) {
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
149
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
150 // Find relevant default bed-heights
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
151 final River river = calcRange.getRiver();
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
152 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
153 if (defaultBedHeights.isEmpty())
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
154 return;
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
155 final DoubleRange scenarioRange = new DoubleRange(partFrom, partTo);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
156 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
157 final Collection<BedHeightsFinder> currentFinders = new ArrayList<>(allFinders);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
158
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
159 // Add historical bed-heights
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
160 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
161 allFinders.add(historicalFinder);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
162 final Collection<Double> stations = BedHeightsUtils.extractStationCollection(allFinders, true);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
163 final List<Double> nulls = new ArrayList<>();
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
164 nulls.add(null);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
165 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
166 for (final Double station : stations) {
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
167 rangeScenarios.put(station, new ArrayList<Double>());
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
168 rangeScenarios.get(station).add(Double.valueOf(bedHeightDifference(station.doubleValue(), currentFinders, historicalFinder)));
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
169 }
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
170 rangeScenarios.put(Double.valueOf(partTo + 0.0001), nulls);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
171 }
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
172
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
173 /**
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
174 * 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
175 */
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
176 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
177 double currentMSH = Double.NaN;
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
178 for (final BedHeightsFinder bhf : currentFinders) {
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
179 currentMSH = bhf.getMeanBedHeight(station);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
180 if (!Double.isNaN(currentMSH))
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
181 break;
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 return Formatter.roundFlowDepth(historicalFinder.getMeanBedHeight(station)).subtract(Formatter.roundFlowDepth(currentMSH)).doubleValue();
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
184 }
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
185
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
186 /**
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
187 * 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
188 */
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
189 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
190 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
191 if (access.getScenario() == ScenarioType.REGIONAL) {
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
192 final int[] deltaws = access.getRegionalScenarioIntegers();
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
193 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
194 if (deltaws[i] != 0)
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
195 labels.add(Integer.toString(deltaws[i]) + " cm");
9391
2da486c7c05f salix historical sounding workflow changed
gernotbelger
parents: 9361
diff changeset
196 } 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
197 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
198 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
199 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
200 return labels.toArray(new String[labels.size()]);
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
201 }
9361
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
202
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
203 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
204 * 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
205 */
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
206 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
207 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
208 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
209 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
210 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
211 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
212 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
213 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
214 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
215 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
216 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
217 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
218 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
219 sep = ", ";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
220 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
221 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
222 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
223 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
224 return "";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
225 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
226
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
227 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
228 * 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
229 */
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
230 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
231 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
232 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
233 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
234 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
235 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
236 sep = ", ";
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 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
239 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
240 if (access.getScenario() == ScenarioType.HISTORICAL) {
9391
2da486c7c05f salix historical sounding workflow changed
gernotbelger
parents: 9361
diff changeset
241 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
242 }
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.SUPRAREGIONAL) {
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 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
247 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
248 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
249 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
250 sep = ", ";
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 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
253 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
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 return "";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
256 }
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
257 }
http://dive4elements.wald.intevation.org