Mercurial > dive4elements > river

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/uinfo/salix/SalixLineCalculation.java @ 9597:5395c6d4ca50

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Softwaretests...20181219 7.3: no interpolation of missing bed heights for Uinfo/Salix historical scenario and B&U/Bzws
author mschaefer
date Tue, 05 Feb 2019 15:47:58 +0100
parents b9c87bbff6a4
children
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;
9529
429b62373633 Reference gauge for salix line always the river's first gauge, all others 'not-in-range'
mschaefer
parents: 9504
diff changeset
34 import org.dive4elements.river.model.Gauge;
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
35 import org.dive4elements.river.model.River;
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
9529
429b62373633 Reference gauge for salix line always the river's first gauge, all others 'not-in-range'
mschaefer
parents: 9504
diff changeset
72 final Gauge firstUpstreamGauge = river.firstUpstreamGauge();
429b62373633 Reference gauge for salix line always the river's first gauge, all others 'not-in-range'
mschaefer
parents: 9504
diff changeset
73 final DoubleRange maxRange = new DoubleRange(firstUpstreamGauge.getRange().getA(), river.getKmUp() ? -99999.999 : 99999.999);
429b62373633 Reference gauge for salix line always the river's first gauge, all others 'not-in-range'
mschaefer
parents: 9504
diff changeset
74 final RiverInfoProvider riverInfoProvider1 = RiverInfoProvider.forRange(this.context, river, maxRange);
429b62373633 Reference gauge for salix line always the river's first gauge, all others 'not-in-range'
mschaefer
parents: 9504
diff changeset
75 final RiverInfoProvider riverInfoProvider = riverInfoProvider1.forReferenceRange(maxRange, false);
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
76
9321
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
77 final SalixLineCalculationResults results = new SalixLineCalculationResults(calcModeLabel, user, riverInfo, range);
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
78
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
79 final SalixLineCalculator calculator = new SalixLineCalculator(riverInfoProvider);
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
80 final NavigableMap<Double, List<Double>> rangeScenarios = buildRangeScenarios(accessSalix);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
81
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
82 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
83 buildAdditionalString(accessSalix), results);
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
84
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
85 return new CalculationResult(results, this.problems);
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
86 }
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
87
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
88 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
89 * 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
90 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
91 private NavigableMap<Double, List<Double>> buildRangeScenarios(final SalixLineAccess access) {
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
92 final NavigableMap<Double, List<Double>> rangeScenarios = new TreeMap<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
93 if (access.getScenario() == ScenarioType.REGIONAL)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
94 fillRangeScenarios(rangeScenarios, access, access.getFromPart().doubleValue(), access.getToPart().doubleValue(),
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
95 access.getRegionalScenarioIntegers());
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
96 else if (access.getScenario() == ScenarioType.SUPRAREGIONAL)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
97 fillRangeScenarios(rangeScenarios, access.getSupraRegionalString());
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
98 else if (access.getScenario() == ScenarioType.HISTORICAL)
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
99 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
100 else
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
101 fillRangeScenarios(rangeScenarios, access);
9321
a978b601a034 Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents: 9318
diff changeset
102
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
103 return rangeScenarios;
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
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
106 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
107 * 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
108 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
109 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
110 final List<Double> nulls = new ArrayList<>();
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
111 nulls.add(null);
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
112 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
113 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
114
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
115 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
116 * 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
117 */
9499
853f2dafc16e VegetationZones in CrossSectionsDiagram
gernotbelger
parents: 9421
diff changeset
118 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final RangeAccess calcRange, final double partFrom,
853f2dafc16e VegetationZones in CrossSectionsDiagram
gernotbelger
parents: 9421
diff changeset
119 final double partTo, final int[] deltaWs) {
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
120 final List<Double> nulls = new ArrayList<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
121 final List<Double> dwsm = new ArrayList<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
122 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
123 nulls.add(null);
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
124 dwsm.add(deltaWs[i] / 100.0);
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 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
127 rangeScenarios.put(Double.valueOf(partFrom - 0.0001), dwsm);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
128 rangeScenarios.put(Double.valueOf(partTo + 0.0001), nulls);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
129 }
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 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
132 * 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
133 * (the zones input by the user cover the calc range completely)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
134 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
135 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
136 final List<SalixZone> parts = SalixZone.parse(zones);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
137 for (final SalixZone part : parts) {
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
138 final List<Double> dwsm = new ArrayList<>();
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
139 if (part.getDwsplValue() == 0)
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
140 dwsm.add(null);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
141 else
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
142 dwsm.add(part.getDwsplValue() / 100.0);
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
143 rangeScenarios.put(Double.valueOf(part.getFromKm().doubleValue() - 0.0001), dwsm);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
144 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
145 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
146
9316
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
147 /**
9573
b9c87bbff6a4 mean bed height -> mean bed LEVEL
gernotbelger
parents: 9529
diff changeset
148 * Fetches historical and reference bed levels and fills a map of delta-MSHs for all fetched stations in the calc range
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
149 */
9499
853f2dafc16e VegetationZones in CrossSectionsDiagram
gernotbelger
parents: 9421
diff changeset
150 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final RangeAccess calcRange, final double partFrom,
853f2dafc16e VegetationZones in CrossSectionsDiagram
gernotbelger
parents: 9421
diff changeset
151 final double partTo, final int historicalBedHeightId) {
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
152
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
153 // Find relevant default bed-heights
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
154 final River river = calcRange.getRiver();
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
155 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
156 if (defaultBedHeights.isEmpty())
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
157 return;
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
158 final DoubleRange scenarioRange = new DoubleRange(partFrom, partTo);
9597
5395c6d4ca50 Softwaretests...20181219 7.3: no interpolation of missing bed heights for Uinfo/Salix historical scenario and B&U/Bzws
mschaefer
parents: 9573
diff changeset
159 final Collection<BedHeightsFinder> allFinders = BedHeightsFinder.createScenarioBedHeights(this.problems, scenarioRange, defaultBedHeights);
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
160 final Collection<BedHeightsFinder> currentFinders = new ArrayList<>(allFinders);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
161
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
162 // Add historical bed-heights
9597
5395c6d4ca50 Softwaretests...20181219 7.3: no interpolation of missing bed heights for Uinfo/Salix historical scenario and B&U/Bzws
mschaefer
parents: 9573
diff changeset
163 final BedHeightsFinder historicalFinder = BedHeightsFinder.forId(this.problems, historicalBedHeightId, scenarioRange, false);
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
164 allFinders.add(historicalFinder);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
165 final Collection<Double> stations = BedHeightsUtils.extractStationCollection(allFinders, true);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
166 final List<Double> nulls = new ArrayList<>();
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
167 nulls.add(null);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
168 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
169 for (final Double station : stations) {
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
170 rangeScenarios.put(station, new ArrayList<Double>());
9421
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
171 final double delta = bedHeightDifference(station, currentFinders, historicalFinder);
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
172 if (Double.isNaN(delta)) {
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
173 rangeScenarios.get(station).add(null);
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
174 if (!this.problems.hasProblems()) {
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
175 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
176 this.problems.addProblem(station, msg);
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
177 }
9499
853f2dafc16e VegetationZones in CrossSectionsDiagram
gernotbelger
parents: 9421
diff changeset
178 } else
9421
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
179 rangeScenarios.get(station).add(Double.valueOf(delta));
9394
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 rangeScenarios.put(Double.valueOf(partTo + 0.0001), nulls);
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
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
184 /**
9573
b9c87bbff6a4 mean bed height -> mean bed LEVEL
gernotbelger
parents: 9529
diff changeset
185 * Gets the difference of a historical bed level against a current one for a station
9394
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 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
188 double currentMSH = Double.NaN;
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
189 for (final BedHeightsFinder bhf : currentFinders) {
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
190 currentMSH = bhf.getMeanBedHeight(station);
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
191 if (!Double.isNaN(currentMSH))
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
192 break;
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
193 }
9421
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
194 if (Double.isNaN(currentMSH))
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
195 return Double.NaN;
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
196 final double historicalMSH = historicalFinder.getMeanBedHeight(station);
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
197 if (Double.isNaN(historicalMSH))
53e26734e4d2 Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents: 9397
diff changeset
198 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
199 return (historicalMSH - currentMSH);
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
200 }
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
201
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
202 /**
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
203 * 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
204 */
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
205 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
206 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
207 if (access.getScenario() == ScenarioType.REGIONAL) {
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
208 final int[] deltaws = access.getRegionalScenarioIntegers();
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
209 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
210 if (deltaws[i] != 0)
72b3270e1568 U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents: 9309
diff changeset
211 labels.add(Integer.toString(deltaws[i]) + " cm");
9391
2da486c7c05f salix historical sounding workflow changed
gernotbelger
parents: 9361
diff changeset
212 } 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
213 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
214 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
215 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
216 return labels.toArray(new String[labels.size()]);
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
217 }
9361
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
218
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
219 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
220 * 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
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 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
223 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
224 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
225 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
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 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
228 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
229 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
230 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
231 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
232 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
233 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
234 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
235 sep = ", ";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
236 }
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 ranges;
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 return "";
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
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
243 /**
9394
439699ff9b2d Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents: 9391
diff changeset
244 * 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
245 */
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
246 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
247 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
248 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
249 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
250 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
251 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
252 sep = ", ";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
253 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
254 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
255 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
256 if (access.getScenario() == ScenarioType.HISTORICAL) {
9391
2da486c7c05f salix historical sounding workflow changed
gernotbelger
parents: 9361
diff changeset
257 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
258 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
259 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
260 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
261 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
262 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
263 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
264 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
265 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
266 sep = ", ";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
267 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
268 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
269 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
270 }
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
271 return "";
2aec052d4088 Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents: 9321
diff changeset
272 }
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
273 }
http://dive4elements.wald.intevation.org