Mercurial > dive4elements > river
annotate artifacts/src/main/java/org/dive4elements/river/artifacts/uinfo/salix/SalixLineCalculation.java @ 9421:53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
author | mschaefer |
---|---|
date | Fri, 17 Aug 2018 11:13:31 +0200 |
parents | 6e7094368e97 |
children | 853f2dafc16e |
rev | line source |
---|---|
8996 | 1 /** Copyright (C) 2017 by Bundesanstalt für Gewässerkunde |
2 * Software engineering by | |
3 * Björnsen Beratende Ingenieure GmbH | |
4 * Dr. Schumacher Ingenieurbüro für Wasser und Umwelt | |
5 * | |
6 * This file is Free Software under the GNU AGPL (>=v3) | |
7 * and comes with ABSOLUTELY NO WARRANTY! Check out the | |
8 * documentation coming with Dive4Elements River for details. | |
9 */ | |
10 package org.dive4elements.river.artifacts.uinfo.salix; | |
11 | |
9309 | 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 | 14 import java.util.List; |
15 import java.util.NavigableMap; | |
16 import java.util.TreeMap; | |
8996 | 17 |
9321
a978b601a034
Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents:
9318
diff
changeset
|
18 import org.apache.commons.lang.math.DoubleRange; |
8996 | 19 import org.dive4elements.artifacts.CallContext; |
9309 | 20 import org.dive4elements.river.artifacts.access.RangeAccess; |
9397 | 21 import org.dive4elements.river.artifacts.common.AbstractResultType; |
9361
2aec052d4088
Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents:
9321
diff
changeset
|
22 import org.dive4elements.river.artifacts.common.GeneralResultType; |
8996 | 23 import org.dive4elements.river.artifacts.model.Calculation; |
24 import org.dive4elements.river.artifacts.model.CalculationResult; | |
25 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
|
26 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
|
27 import org.dive4elements.river.artifacts.sinfo.tkhstate.BedHeightsFinder; |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
28 import org.dive4elements.river.artifacts.sinfo.tkhstate.BedHeightsUtils; |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
29 import org.dive4elements.river.artifacts.sinfo.tkhstate.DefaultBedHeights; |
8996 | 30 import org.dive4elements.river.artifacts.sinfo.util.CalculationUtils; |
31 import org.dive4elements.river.artifacts.sinfo.util.RiverInfo; | |
32 import org.dive4elements.river.artifacts.uinfo.UINFOArtifact; | |
9309 | 33 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
|
34 import org.dive4elements.river.model.BedHeight; |
8996 | 35 import org.dive4elements.river.model.River; |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
36 import org.dive4elements.river.utils.Formatter; |
8996 | 37 |
38 /** | |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
39 * 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
|
40 * |
8996 | 41 * @author Domenico Nardi Tironi |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
42 * @author Matthias Schäfer |
8996 | 43 * |
44 */ | |
45 final class SalixLineCalculation { | |
46 | |
47 private final CallContext context; | |
48 | |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
49 private Calculation problems; |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
50 |
8996 | 51 public SalixLineCalculation(final CallContext context) { |
52 this.context = context; | |
53 } | |
54 | |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
55 /** |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
56 * 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
|
57 */ |
8996 | 58 public CalculationResult calculate(final UINFOArtifact uinfo) { |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
59 this.problems = new Calculation(); |
8996 | 60 |
61 final String calcModeLabel = Resources.getMsg(this.context.getMeta(), uinfo.getCalculationMode().name()); | |
62 final String user = CalculationUtils.findArtifactUser(this.context, uinfo); | |
63 | |
9070
611a523fc42f
VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents:
9069
diff
changeset
|
64 final SalixLineAccess accessSalix = new SalixLineAccess(uinfo); |
611a523fc42f
VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents:
9069
diff
changeset
|
65 |
611a523fc42f
VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents:
9069
diff
changeset
|
66 final River river = accessSalix.getRiver(); |
8996 | 67 final RiverInfo riverInfo = new RiverInfo(river); |
68 | |
9321
a978b601a034
Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents:
9318
diff
changeset
|
69 final DoubleRange range = accessSalix.getRange(); |
9068
a5cf8d7bff3c
access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents:
8996
diff
changeset
|
70 |
9321
a978b601a034
Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents:
9318
diff
changeset
|
71 final ScenarioType scenarioType = accessSalix.getScenario(); |
9070
611a523fc42f
VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents:
9069
diff
changeset
|
72 |
9397 | 73 final RiverInfoProvider riverInfoProvider = RiverInfoProvider.forRange(this.context, river, range, true); |
8996 | 74 |
9321
a978b601a034
Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents:
9318
diff
changeset
|
75 final SalixLineCalculationResults results = new SalixLineCalculationResults(calcModeLabel, user, riverInfo, range); |
a978b601a034
Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents:
9318
diff
changeset
|
76 |
a978b601a034
Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents:
9318
diff
changeset
|
77 final SalixLineCalculator calculator = new SalixLineCalculator(riverInfoProvider); |
9309 | 78 final NavigableMap<Double, List<Double>> rangeScenarios = buildRangeScenarios(accessSalix); |
79 | |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
80 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
|
81 buildAdditionalString(accessSalix), results); |
9243 | 82 |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
83 return new CalculationResult(results, this.problems); |
8996 | 84 } |
9243 | 85 |
9309 | 86 /** |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
87 * Builds a map of delta-Ws by from-km for the selected scenario |
9309 | 88 */ |
89 private NavigableMap<Double, List<Double>> buildRangeScenarios(final SalixLineAccess access) { | |
90 final NavigableMap<Double, List<Double>> rangeScenarios = new TreeMap<>(); | |
91 if (access.getScenario() == ScenarioType.REGIONAL) | |
92 fillRangeScenarios(rangeScenarios, access, access.getFromPart().doubleValue(), access.getToPart().doubleValue(), | |
93 access.getRegionalScenarioIntegers()); | |
94 else if (access.getScenario() == ScenarioType.SUPRAREGIONAL) | |
95 fillRangeScenarios(rangeScenarios, access.getSupraRegionalString()); | |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
96 else if (access.getScenario() == ScenarioType.HISTORICAL) |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
97 fillRangeScenarios(rangeScenarios, access, access.getFromPart().doubleValue(), access.getToPart().doubleValue(), access.getBedHeightId()); |
9309 | 98 else |
99 fillRangeScenarios(rangeScenarios, access); | |
9321
a978b601a034
Salix: Fixed ArrrayoutOfBoundsException; minor cleanup
gernotbelger
parents:
9318
diff
changeset
|
100 |
9309 | 101 return rangeScenarios; |
102 } | |
103 | |
104 /** | |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
105 * Fills a map of delta-Ws with only one 0-delta for the whole calc range (no scenario) |
9309 | 106 */ |
107 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final RangeAccess calcRange) { | |
108 final List<Double> nulls = new ArrayList<>(); | |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
109 nulls.add(null); |
9309 | 110 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls); |
111 } | |
112 | |
113 /** | |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
114 * Fills a map of delta-Ws by km-range from the regional scenario input data |
9309 | 115 */ |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
116 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
|
117 final double partFrom, final double partTo, final int[] deltaWs) { |
9309 | 118 final List<Double> nulls = new ArrayList<>(); |
119 final List<Double> dwsm = new ArrayList<>(); | |
120 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
|
121 nulls.add(null); |
9309 | 122 dwsm.add(deltaWs[i] / 100.0); |
123 } | |
124 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls); | |
125 rangeScenarios.put(Double.valueOf(partFrom - 0.0001), dwsm); | |
126 rangeScenarios.put(Double.valueOf(partTo + 0.0001), nulls); | |
127 } | |
128 | |
129 /** | |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
130 * Fills a map of delta-Ws by km-range from the supraregional scenario input data |
9309 | 131 * (the zones input by the user cover the calc range completely) |
132 */ | |
133 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final String zones) { | |
134 final List<SalixZone> parts = SalixZone.parse(zones); | |
135 for (final SalixZone part : parts) { | |
136 final List<Double> dwsm = new ArrayList<>(); | |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
137 if (part.getDwsplValue() == 0) |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
138 dwsm.add(null); |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
139 else |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
140 dwsm.add(part.getDwsplValue() / 100.0); |
9309 | 141 rangeScenarios.put(Double.valueOf(part.getFromKm().doubleValue() - 0.0001), dwsm); |
142 } | |
143 } | |
144 | |
9316
72b3270e1568
U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents:
9309
diff
changeset
|
145 /** |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
146 * 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
|
147 */ |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
148 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
|
149 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
|
150 |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
151 // Find relevant default bed-heights |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
152 final River river = calcRange.getRiver(); |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
153 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
|
154 if (defaultBedHeights.isEmpty()) |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
155 return; |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
156 final DoubleRange scenarioRange = new DoubleRange(partFrom, partTo); |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
157 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
|
158 final Collection<BedHeightsFinder> currentFinders = new ArrayList<>(allFinders); |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
159 |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
160 // Add historical bed-heights |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
161 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
|
162 allFinders.add(historicalFinder); |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
163 final Collection<Double> stations = BedHeightsUtils.extractStationCollection(allFinders, true); |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
164 final List<Double> nulls = new ArrayList<>(); |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
165 nulls.add(null); |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
166 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
|
167 for (final Double station : stations) { |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
168 rangeScenarios.put(station, new ArrayList<Double>()); |
9421
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
169 final double delta = bedHeightDifference(station, currentFinders, historicalFinder); |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
170 if (Double.isNaN(delta)) { |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
171 rangeScenarios.get(station).add(null); |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
172 if (!this.problems.hasProblems()) { |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
173 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
|
174 this.problems.addProblem(station, msg); |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
175 } |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
176 } |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
177 else |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
178 rangeScenarios.get(station).add(Double.valueOf(delta)); |
9394
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 rangeScenarios.put(Double.valueOf(partTo + 0.0001), nulls); |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
181 } |
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 * 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
|
185 */ |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
186 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
|
187 double currentMSH = Double.NaN; |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
188 for (final BedHeightsFinder bhf : currentFinders) { |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
189 currentMSH = bhf.getMeanBedHeight(station); |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
190 if (!Double.isNaN(currentMSH)) |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
191 break; |
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
192 } |
9421
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
193 if (Double.isNaN(currentMSH)) |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
194 return Double.NaN; |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
195 final double historicalMSH = historicalFinder.getMeanBedHeight(station); |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
196 if (Double.isNaN(historicalMSH)) |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
197 return Double.NaN; |
53e26734e4d2
Fixed NaN handling in salix historical scenarion calculation, added problem report
mschaefer
parents:
9397
diff
changeset
|
198 return Formatter.roundFlowDepth(historicalMSH).subtract(Formatter.roundFlowDepth(currentMSH)).doubleValue(); |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
199 } |
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 * 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
|
203 */ |
72b3270e1568
U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents:
9309
diff
changeset
|
204 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
|
205 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
|
206 if (access.getScenario() == ScenarioType.REGIONAL) { |
72b3270e1568
U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents:
9309
diff
changeset
|
207 final int[] deltaws = access.getRegionalScenarioIntegers(); |
72b3270e1568
U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents:
9309
diff
changeset
|
208 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
|
209 if (deltaws[i] != 0) |
72b3270e1568
U-Info salix line with regional and supraregional scenario calculation and chart display
mschaefer
parents:
9309
diff
changeset
|
210 labels.add(Integer.toString(deltaws[i]) + " cm"); |
9391 | 211 } 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
|
212 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
|
213 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
|
214 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
|
215 return labels.toArray(new String[labels.size()]); |
9243 | 216 } |
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 /** |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
219 * 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
|
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 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
|
222 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
|
223 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
|
224 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
|
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 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
|
227 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
|
228 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
|
229 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
|
230 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
|
231 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
|
232 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
|
233 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
|
234 sep = ", "; |
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 } |
2aec052d4088
Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents:
9321
diff
changeset
|
237 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
|
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 return ""; |
2aec052d4088
Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents:
9321
diff
changeset
|
240 } |
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 /** |
9394
439699ff9b2d
Added U-Info iota (prev. salix) calculation for historical scenario
mschaefer
parents:
9391
diff
changeset
|
243 * 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
|
244 */ |
2aec052d4088
Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents:
9321
diff
changeset
|
245 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
|
246 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
|
247 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
|
248 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
|
249 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
|
250 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
|
251 sep = ", "; |
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 if (access.getScenario() == ScenarioType.HISTORICAL) { |
9391 | 256 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
|
257 } |
2aec052d4088
Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents:
9321
diff
changeset
|
258 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
|
259 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
|
260 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
|
261 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
|
262 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
|
263 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
|
264 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
|
265 sep = ", "; |
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 } |
2aec052d4088
Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents:
9321
diff
changeset
|
268 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
|
269 } |
2aec052d4088
Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents:
9321
diff
changeset
|
270 return ""; |
2aec052d4088
Refactoring on SalixLineCalculationResult etc., calculation using MQ etc., scenario (regional+supra) csv columns and meta data added
mschaefer
parents:
9321
diff
changeset
|
271 } |
9397 | 272 |
273 /** | |
274 * Fetches a iota or waterlevel height of a station from a salix calculation result | |
275 */ | |
276 public double fetchStationHeight(final double station, final AbstractResultType resultType, final SalixLineCalculationResult result) { | |
277 | |
278 final Calculation problems = new Calculation(); | |
279 | |
280 final SalixLineCalculator calculator = new SalixLineCalculator(null); | |
281 return calculator.fetchStationHeight(problems, station, resultType, result); | |
282 } | |
8996 | 283 } |