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
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;
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
13 import java.util.Iterator;
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
14 import java.util.List;
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
15 import java.util.Map;
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
16 import java.util.NavigableMap;
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
17 import java.util.TreeMap;
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
18
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
19 import org.apache.commons.lang.math.DoubleRange;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
20 import org.dive4elements.artifacts.CallContext;
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
21 import org.dive4elements.river.artifacts.access.RangeAccess;
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
22 import org.dive4elements.river.artifacts.model.Calculation;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
23 import org.dive4elements.river.artifacts.model.CalculationResult;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
24 import org.dive4elements.river.artifacts.resources.Resources;
9295
385b52ccde23 Work on U-Info salix line calculation and chart (no scenario case)
mschaefer
parents: 9249
diff changeset
25 import org.dive4elements.river.artifacts.sinfo.common.RiverInfoProvider;
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
26 import org.dive4elements.river.artifacts.sinfo.util.CalculationUtils;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
27 import org.dive4elements.river.artifacts.sinfo.util.RiverInfo;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
28 import org.dive4elements.river.artifacts.uinfo.UINFOArtifact;
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
29 import org.dive4elements.river.artifacts.uinfo.salix.SalixLineAccess.ScenarioType;
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
30 import org.dive4elements.river.model.River;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
31
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
32 /**
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
33 * @author Domenico Nardi Tironi
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
34 *
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
35 */
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
36 final class SalixLineCalculation {
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
37
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
38 private final CallContext context;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
39
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
40 public SalixLineCalculation(final CallContext context) {
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
41 this.context = context;
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 public CalculationResult calculate(final UINFOArtifact uinfo) {
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
45 final Calculation problems = new Calculation();
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
46
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
47 final String calcModeLabel = Resources.getMsg(this.context.getMeta(), uinfo.getCalculationMode().name());
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
48 final String user = CalculationUtils.findArtifactUser(this.context, uinfo);
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
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
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
53 final RiverInfo riverInfo = new RiverInfo(river);
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
54
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
55 // final DoubleRange range = accessSalix.getRange();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
56 // final boolean useScenario = accessSalix.getUseScenario();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
57 // final String selectedScenario = accessSalix.getScenario();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
58 // final Double fromPart = accessSalix.getFromPart();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
59 // final Double toPart = accessSalix.getToPart();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
60 //
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
61 // final boolean useSCenario = accessSalix.getUseScenario();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
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
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
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
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
72 final SalixLineCalculationResults results = new SalixLineCalculationResults(calcModeLabel, user, riverInfo, accessSalix.getRange());
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
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
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
75 final NavigableMap<Double, List<Double>> rangeScenarios = buildRangeScenarios(accessSalix);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
76
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
77 calculator.execute(problems, uinfo, rangeScenarios, accessSalix.getScenario(), results);
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
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
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
153 return new CalculationResult(results, problems);
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
154 }
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
155
9309
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
156 /**
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
157 * Build a map of delta-Ws by km-range for the selected scenario
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
158 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
159 private NavigableMap<Double, List<Double>> buildRangeScenarios(final SalixLineAccess access) {
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
160 final NavigableMap<Double, List<Double>> rangeScenarios = new TreeMap<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
161 if (access.getScenario() == ScenarioType.REGIONAL)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
162 fillRangeScenarios(rangeScenarios, access, access.getFromPart().doubleValue(), access.getToPart().doubleValue(),
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
163 access.getRegionalScenarioIntegers());
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
164 else if (access.getScenario() == ScenarioType.SUPRAREGIONAL)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
165 fillRangeScenarios(rangeScenarios, access.getSupraRegionalString());
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
166 // TODO else if (access.getScenario().equals(ScenarioType.HISTORICAL.getKey()))
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
167 // historisches Szenario aus MSH etc.
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
168 else
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
169 fillRangeScenarios(rangeScenarios, access);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
170 return rangeScenarios;
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
171 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
172
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
173 /**
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
174 * Fill a map of delta-Ws with only one 0-delta for the whole calc range (no scenario)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
175 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
176 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
177 final List<Double> nulls = new ArrayList<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
178 nulls.add(0.0);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
179 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
180 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
181
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
182 /**
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
183 * Fill a map of delta-Ws by km-range from the regional scenario input data
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
184 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
185 private void fillRangeScenarios(final NavigableMap<Double, List<Double>> rangeScenarios, final RangeAccess calcRange, final double partFrom,
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
186 final double partTo, final int[] deltaWs) {
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
187 final List<Double> nulls = new ArrayList<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
188 final List<Double> dwsm = new ArrayList<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
189 for (int i = 0; i <= deltaWs.length - 1; i++) {
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
190 nulls.add(0.0);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
191 dwsm.add(deltaWs[i] / 100.0);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
192 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
193 rangeScenarios.put(Double.valueOf(calcRange.getLowerKm() - 0.0001), nulls);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
194 rangeScenarios.put(Double.valueOf(partFrom - 0.0001), dwsm);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
195 rangeScenarios.put(Double.valueOf(partTo + 0.0001), nulls);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
196 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
197
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
198 /**
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
199 * Fill a map of delta-Ws by km-range from the supraregional scenario input data
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
200 * (the zones input by the user cover the calc range completely)
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
201 */
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
202 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
203 final List<SalixZone> parts = SalixZone.parse(zones);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
204 for (final SalixZone part : parts) {
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
205 final List<Double> dwsm = new ArrayList<>();
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
206 dwsm.add(part.getDwsplValue() / 100.0);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
207 rangeScenarios.put(Double.valueOf(part.getFromKm().doubleValue() - 0.0001), dwsm);
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
208 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
209 }
9a9f076d5716 Work on U-Info salix line calculation
mschaefer
parents: 9295
diff changeset
210
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
211 private SalixScenario findScenarioByKm(final double km, final Map<DoubleRange, SalixScenario> rangeScenarioMap) {
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
212 final Iterator<DoubleRange> rangeIterator = rangeScenarioMap.keySet().iterator();
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
213 while (rangeIterator.hasNext()) {
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
214 final DoubleRange range = rangeIterator.next();
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
215 if (range.containsDouble(km + 0.0001)) {
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
216 return rangeScenarioMap.get(range);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
217 }
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
218 }
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
219 return null;
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
220 }
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
221 }

http://dive4elements.wald.intevation.org