Mercurial > dive4elements > river
annotate artifacts/src/main/java/org/dive4elements/river/artifacts/uinfo/salix/SalixLineCalculation.java @ 9309:9a9f076d5716
Work on U-Info salix line calculation
author | mschaefer |
---|---|
date | Wed, 25 Jul 2018 19:29:25 +0200 |
parents | 385b52ccde23 |
children | 72b3270e1568 |
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; |
9243 | 13 import java.util.Iterator; |
9309 | 14 import java.util.List; |
9243 | 15 import java.util.Map; |
9309 | 16 import java.util.NavigableMap; |
17 import java.util.TreeMap; | |
8996 | 18 |
19 import org.apache.commons.lang.math.DoubleRange; | |
20 import org.dive4elements.artifacts.CallContext; | |
9309 | 21 import org.dive4elements.river.artifacts.access.RangeAccess; |
8996 | 22 import org.dive4elements.river.artifacts.model.Calculation; |
23 import org.dive4elements.river.artifacts.model.CalculationResult; | |
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; |
8996 | 26 import org.dive4elements.river.artifacts.sinfo.util.CalculationUtils; |
27 import org.dive4elements.river.artifacts.sinfo.util.RiverInfo; | |
28 import org.dive4elements.river.artifacts.uinfo.UINFOArtifact; | |
9309 | 29 import org.dive4elements.river.artifacts.uinfo.salix.SalixLineAccess.ScenarioType; |
8996 | 30 import org.dive4elements.river.model.River; |
31 | |
32 /** | |
33 * @author Domenico Nardi Tironi | |
34 * | |
35 */ | |
36 final class SalixLineCalculation { | |
37 | |
38 private final CallContext context; | |
39 | |
40 public SalixLineCalculation(final CallContext context) { | |
41 this.context = context; | |
42 } | |
43 | |
44 public CalculationResult calculate(final UINFOArtifact uinfo) { | |
45 final Calculation problems = new Calculation(); | |
46 | |
47 final String calcModeLabel = Resources.getMsg(this.context.getMeta(), uinfo.getCalculationMode().name()); | |
48 final String user = CalculationUtils.findArtifactUser(this.context, uinfo); | |
49 | |
9070
611a523fc42f
VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents:
9069
diff
changeset
|
50 final SalixLineAccess accessSalix = new SalixLineAccess(uinfo); |
611a523fc42f
VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents:
9069
diff
changeset
|
51 |
611a523fc42f
VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents:
9069
diff
changeset
|
52 final River river = accessSalix.getRiver(); |
8996 | 53 final RiverInfo riverInfo = new RiverInfo(river); |
54 | |
9309 | 55 // final DoubleRange range = accessSalix.getRange(); |
56 // final boolean useScenario = accessSalix.getUseScenario(); | |
57 // final String selectedScenario = accessSalix.getScenario(); | |
58 // final Double fromPart = accessSalix.getFromPart(); | |
59 // final Double toPart = accessSalix.getToPart(); | |
60 // | |
61 // final boolean useSCenario = accessSalix.getUseScenario(); | |
62 // final String scenario = accessSalix.getScenario(); | |
9068
a5cf8d7bff3c
access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents:
8996
diff
changeset
|
63 // calculation_mode |
a5cf8d7bff3c
access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents:
8996
diff
changeset
|
64 // ld_from , ld_to |
a5cf8d7bff3c
access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents:
8996
diff
changeset
|
65 // use_scenario (boolean) |
a5cf8d7bff3c
access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents:
8996
diff
changeset
|
66 // ld_from_part; ld_to_part |
a5cf8d7bff3c
access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents:
8996
diff
changeset
|
67 // scenario_selection (mögliche Werte:"scenarioType.option1" "scenarioType.option2" "scenarioType.option3" |
9309 | 68 final RiverInfoProvider riverInfoProvider = RiverInfoProvider.forRange(this.context, river, accessSalix.getRange()); |
9068
a5cf8d7bff3c
access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents:
8996
diff
changeset
|
69 |
9070
611a523fc42f
VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents:
9069
diff
changeset
|
70 // FIXME: real calculation |
611a523fc42f
VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents:
9069
diff
changeset
|
71 |
9309 | 72 final SalixLineCalculationResults results = new SalixLineCalculationResults(calcModeLabel, user, riverInfo, accessSalix.getRange()); |
8996 | 73 |
9295
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
74 final SalixLineCalculator calculator = new SalixLineCalculator(this.context, riverInfoProvider); |
9309 | 75 final NavigableMap<Double, List<Double>> rangeScenarios = buildRangeScenarios(accessSalix); |
76 | |
77 calculator.execute(problems, uinfo, rangeScenarios, accessSalix.getScenario(), results); | |
9243 | 78 |
9295
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
79 // final Collection<ResultRow> rows = new ArrayList<>(); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
80 // SalixLineCalculationNoScenarioResult result = null; // verzweigung etwas ungünstig. möglicherweise auch die |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
81 // Abstraktion. ist erstmal nur ne idee |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
82 // final ResultRow row1 = ResultRow.create(). // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
83 // putValue(GeneralResultType.station, 100).// |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
84 // putValue(UInfoResultType.salixline, 28).// |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
85 // putValue(UInfoResultType.salix_delta_mw, 2); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
86 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
87 // if (!useScenario) { |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
88 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
89 // rows.add(row1); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
90 // result = new SalixLineCalculationNoScenarioResult("Ergebnis 1", null, rows); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
91 // } else { |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
92 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
93 // if (scenario.equals("scenarioType.option1")) { // REGIONAL |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
94 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
95 // final int[] scenarios = accessSalix.getRegionalScenarioIntegers(); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
96 // final List<SalixScenario> list = new ArrayList<>(); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
97 // for (final int scen_val : scenarios) { |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
98 // list.add(new SalixScenario(scen_val, 666.));// TODO: replace 666 by real calculated value |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
99 // } |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
100 // row1.putValue(UInfoResultType.customMultiRowColSalixRegionalValue_Dwspl, list);// |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
101 // rows.add(row1); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
102 // result = new SalixLineCalculationRegionalResult("Ergebnis 1 regional test", null, rows, scenarios); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
103 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
104 // } else if (scenario.equals("scenarioType.option2")) { // SUPRA-REGIONAL |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
105 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
106 // final String supraRegional = accessSalix.getSupraRegionalString(); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
107 // final List<SalixZone> list = SalixZone.parse(supraRegional); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
108 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
109 // final Map<DoubleRange, SalixScenario> rangeScenarioMap = new HashMap<>(); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
110 // // make double range |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
111 // for (int i = 0; i < list.size(); i++) { |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
112 // final SalixZone zone = list.get(i); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
113 // final double upper = i < list.size() - 1 ? (zone.getUpperFromTo() - 0.0001) : zone.getUpperFromTo() + 0.0001;// |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
114 // "halboffenes Intervall |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
115 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
116 // final DoubleRange zonerange = new DoubleRange((double) zone.getLowerFromTo(), upper); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
117 // final double salixValue = 666.;// TODO: calculate the salix value |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
118 // final SalixScenario salixscen = new SalixScenario(zone.getDwsplValue(), salixValue); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
119 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
120 // rangeScenarioMap.put(zonerange, salixscen); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
121 // } |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
122 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
123 // // make calculation |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
124 // double currentKm = range.getMinimumDouble(); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
125 // final double step = 0.1; // TODO: get from global setting? |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
126 // while (currentKm < range.getMaximumDouble()) { |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
127 // final ResultRow rowSupraRegional = ResultRow.create(). // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
128 // putValue(GeneralResultType.station, currentKm).// |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
129 // putValue(UInfoResultType.salixline, 28).// |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
130 // putValue(UInfoResultType.salix_delta_mw, 2); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
131 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
132 // final SalixScenario scenarioCurrentKm = findScenarioByKm(currentKm, rangeScenarioMap); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
133 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
134 // if (scenarioCurrentKm != null) { // should not happen, scenarioCurrentKm == null -> BUG |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
135 // rowSupraRegional.putValue(UInfoResultType.salix_line_scenario, scenarioCurrentKm.getSalixValue()); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
136 // rowSupraRegional.putValue(UInfoResultType.salix_line_scenario_dwspl, scenarioCurrentKm.getDwspl()); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
137 // } |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
138 // rows.add(rowSupraRegional); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
139 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
140 // currentKm = currentKm + step; |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
141 // } |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
142 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
143 // result = new SalixLineCalculationSupraRegionalResult("Ergebnis 1 supra regional test", null, rows); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
144 // |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
145 // } else if (scenario.equals("scenarioType.option3")) { // HISTORICAL |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
146 // row1.putValue(UInfoResultType.salixlinehist, 66).// |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
147 // putValue(UInfoResultType.salix_line_scenario_dwspl, 88); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
148 // rows.add(row1); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
149 // result = new SalixLineCalculationHistoricalResult("Ergebnis 1 historical test", null, rows); |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
150 // } |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
151 // } |
385b52ccde23
Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents:
9249
diff
changeset
|
152 // results.addResult(result, problems); |
8996 | 153 return new CalculationResult(results, problems); |
154 } | |
9243 | 155 |
9309 | 156 /** |
157 * Build a map of delta-Ws by km-range for the selected scenario | |
158 */ | |
159 private NavigableMap<Double, List<Double>> buildRangeScenarios(final SalixLineAccess access) { | |
160 final NavigableMap<Double, List<Double>> rangeScenarios = new TreeMap<>(); | |
161 if (access.getScenario() == ScenarioType.REGIONAL) | |
162 fillRangeScenarios(rangeScenarios, access, access.getFromPart().doubleValue(), access.getToPart().doubleValue(), | |
163 access.getRegionalScenarioIntegers()); | |
164 else if (access.getScenario() == ScenarioType.SUPRAREGIONAL) | |
165 fillRangeScenarios(rangeScenarios, access.getSupraRegionalString()); | |
166 // TODO else if (access.getScenario().equals(ScenarioType.HISTORICAL.getKey())) | |
167 // historisches Szenario aus MSH etc. | |
168 else | |
169 fillRangeScenarios(rangeScenarios, access); | |
170 return rangeScenarios; | |
171 } | |
172 | |
173 /** | |
174 * Fill a map of delta-Ws with only one 0-delta for the whole calc range (no scenario) | |
175 */ | |
176 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final RangeAccess calcRange) { | |
177 final List<Double> nulls = new ArrayList<>(); | |
178 nulls.add(0.0); | |
179 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls); | |
180 } | |
181 | |
182 /** | |
183 * Fill a map of delta-Ws by km-range from the regional scenario input data | |
184 */ | |
185 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final RangeAccess calcRange, final double partFrom, | |
186 final double partTo, final int[] deltaWs) { | |
187 final List<Double> nulls = new ArrayList<>(); | |
188 final List<Double> dwsm = new ArrayList<>(); | |
189 for (int i = 0; i <= deltaWs.length - 1; i++) { | |
190 nulls.add(0.0); | |
191 dwsm.add(deltaWs[i] / 100.0); | |
192 } | |
193 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls); | |
194 rangeScenarios.put(Double.valueOf(partFrom - 0.0001), dwsm); | |
195 rangeScenarios.put(Double.valueOf(partTo + 0.0001), nulls); | |
196 } | |
197 | |
198 /** | |
199 * Fill a map of delta-Ws by km-range from the supraregional scenario input data | |
200 * (the zones input by the user cover the calc range completely) | |
201 */ | |
202 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final String zones) { | |
203 final List<SalixZone> parts = SalixZone.parse(zones); | |
204 for (final SalixZone part : parts) { | |
205 final List<Double> dwsm = new ArrayList<>(); | |
206 dwsm.add(part.getDwsplValue() / 100.0); | |
207 rangeScenarios.put(Double.valueOf(part.getFromKm().doubleValue() - 0.0001), dwsm); | |
208 } | |
209 } | |
210 | |
9243 | 211 private SalixScenario findScenarioByKm(final double km, final Map<DoubleRange, SalixScenario> rangeScenarioMap) { |
212 final Iterator<DoubleRange> rangeIterator = rangeScenarioMap.keySet().iterator(); | |
213 while (rangeIterator.hasNext()) { | |
214 final DoubleRange range = rangeIterator.next(); | |
215 if (range.containsDouble(km + 0.0001)) { | |
216 return rangeScenarioMap.get(range); | |
217 } | |
218 } | |
219 return null; | |
220 } | |
8996 | 221 } |