annotate artifacts/src/main/java/org/dive4elements/river/artifacts/uinfo/salix/SalixLineCalculation.java @ 9249:600e1ac42049

Historical Result
author gernotbelger
date Thu, 12 Jul 2018 14:15:36 +0200
parents c08d5cfa4981
children 385b52ccde23
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
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
12 import java.util.ArrayList;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
13 import java.util.Collection;
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
14 import java.util.HashMap;
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
15 import java.util.Iterator;
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
16 import java.util.List;
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
17 import java.util.Map;
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;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
21 import org.dive4elements.river.artifacts.common.GeneralResultType;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
22 import org.dive4elements.river.artifacts.common.ResultRow;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
23 import org.dive4elements.river.artifacts.model.Calculation;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
24 import org.dive4elements.river.artifacts.model.CalculationResult;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
25 import org.dive4elements.river.artifacts.resources.Resources;
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;
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
29 import org.dive4elements.river.artifacts.uinfo.commons.UInfoResultType;
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
9068
a5cf8d7bff3c access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents: 8996
diff changeset
55 final DoubleRange range = accessSalix.getRange();
a5cf8d7bff3c access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents: 8996
diff changeset
56 final boolean useScenario = accessSalix.getUseScenario();
a5cf8d7bff3c access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents: 8996
diff changeset
57 final String selectedScenario = accessSalix.getScenario();
9222
5a0e8d76582e uinfo.salix regional workflow
gernotbelger
parents: 9070
diff changeset
58 final Double fromPart = accessSalix.getFromPart();
5a0e8d76582e uinfo.salix regional workflow
gernotbelger
parents: 9070
diff changeset
59 final Double toPart = accessSalix.getToPart();
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
60
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
61 final boolean useSCenario = accessSalix.getUseScenario();
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
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"
a5cf8d7bff3c access uinfo.salx; tablePanel edit/non-edit merge
gernotbelger
parents: 8996
diff changeset
68
9070
611a523fc42f VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents: 9069
diff changeset
69 // FIXME: real calculation
611a523fc42f VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents: 9069
diff changeset
70
611a523fc42f VegetationZoneAccessHelper, VegetationTablePanels verbessert
gernotbelger
parents: 9069
diff changeset
71 final SalixLineCalculationResults results = new SalixLineCalculationResults(calcModeLabel, user, riverInfo, range);
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
72
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
73 final Collection<ResultRow> rows = new ArrayList<>();
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
74 SalixLineCalculationNoScenarioResult result = null; // verzweigung etwas ungünstig. möglicherweise auch die Abstraktion. ist erstmal nur ne idee
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
75 final ResultRow row1 = ResultRow.create(). //
9069
1ffd38826175 access uinfo.vegetationzones+inundation_duration
gernotbelger
parents: 9068
diff changeset
76 putValue(GeneralResultType.station, 100).//
1ffd38826175 access uinfo.vegetationzones+inundation_duration
gernotbelger
parents: 9068
diff changeset
77 putValue(UInfoResultType.salixline, 28).//
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
78 putValue(UInfoResultType.salix_delta_mw, 2);
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
79
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
80 if (!useScenario) {
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
81
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
82 rows.add(row1);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
83 result = new SalixLineCalculationNoScenarioResult("Ergebnis 1", null, rows);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
84 } else {
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
85
9246
c08d5cfa4981 some hibernate queries on bedheigts for salix
gernotbelger
parents: 9244
diff changeset
86 if (scenario.equals("scenarioType.option1")) { // REGIONAL
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
87
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
88 final int[] scenarios = accessSalix.getRegionalScenarioIntegers();
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
89 final List<SalixScenario> list = new ArrayList<>();
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
90 for (final int scen_val : scenarios) {
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
91 list.add(new SalixScenario(scen_val, 666.));// TODO: replace 666 by real calculated value
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
92 }
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
93 row1.putValue(UInfoResultType.customMultiRowColSalixRegionalValue_Dwspl, list);//
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
94 rows.add(row1);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
95 result = new SalixLineCalculationRegionalResult("Ergebnis 1 regional test", null, rows, scenarios);
9246
c08d5cfa4981 some hibernate queries on bedheigts for salix
gernotbelger
parents: 9244
diff changeset
96
c08d5cfa4981 some hibernate queries on bedheigts for salix
gernotbelger
parents: 9244
diff changeset
97 } else if (scenario.equals("scenarioType.option2")) { // SUPRA-REGIONAL
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
98
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
99 final String supraRegional = accessSalix.getSupraRegionalString();
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
100 final List<SalixZone> list = SalixZone.parse(supraRegional);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
101
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
102 final Map<DoubleRange, SalixScenario> rangeScenarioMap = new HashMap<>();
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
103 // make double range
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
104 for (int i = 0; i < list.size(); i++) {
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
105 final SalixZone zone = list.get(i);
9244
b9505aa78cdd commit error
gernotbelger
parents: 9243
diff changeset
106 final double upper = i < list.size() - 1 ? (zone.getUpperFromTo() - 0.0001) : zone.getUpperFromTo() + 0.0001;// "halboffenes Intervall
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
107
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
108 final DoubleRange zonerange = new DoubleRange((double) zone.getLowerFromTo(), upper);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
109 final double salixValue = 666.;// TODO: calculate the salix value
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
110 final SalixScenario salixscen = new SalixScenario(zone.getDwsplValue(), salixValue);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
111
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
112 rangeScenarioMap.put(zonerange, salixscen);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
113 }
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
114
9246
c08d5cfa4981 some hibernate queries on bedheigts for salix
gernotbelger
parents: 9244
diff changeset
115 // make calculation
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
116 double currentKm = range.getMinimumDouble();
9246
c08d5cfa4981 some hibernate queries on bedheigts for salix
gernotbelger
parents: 9244
diff changeset
117 final double step = 0.1; // TODO: get from global setting?
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
118 while (currentKm < range.getMaximumDouble()) {
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
119 final ResultRow rowSupraRegional = ResultRow.create(). //
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
120 putValue(GeneralResultType.station, currentKm).//
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
121 putValue(UInfoResultType.salixline, 28).//
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
122 putValue(UInfoResultType.salix_delta_mw, 2);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
123
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
124 final SalixScenario scenarioCurrentKm = findScenarioByKm(currentKm, rangeScenarioMap);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
125
9246
c08d5cfa4981 some hibernate queries on bedheigts for salix
gernotbelger
parents: 9244
diff changeset
126 if (scenarioCurrentKm != null) { // should not happen, scenarioCurrentKm == null -> BUG
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
127 rowSupraRegional.putValue(UInfoResultType.salix_line_scenario, scenarioCurrentKm.getSalixValue());
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
128 rowSupraRegional.putValue(UInfoResultType.salix_line_scenario_dwspl, scenarioCurrentKm.getDwspl());
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
129 }
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
130 rows.add(rowSupraRegional);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
131
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
132 currentKm = currentKm + step;
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
133 }
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
134
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
135 result = new SalixLineCalculationSupraRegionalResult("Ergebnis 1 supra regional test", null, rows);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
136
9246
c08d5cfa4981 some hibernate queries on bedheigts for salix
gernotbelger
parents: 9244
diff changeset
137 } else if (scenario.equals("scenarioType.option3")) { // HISTORICAL
9249
600e1ac42049 Historical Result
gernotbelger
parents: 9246
diff changeset
138 row1.putValue(UInfoResultType.salixlinehist, 66).//
600e1ac42049 Historical Result
gernotbelger
parents: 9246
diff changeset
139 putValue(UInfoResultType.salix_line_scenario_dwspl, 88);
600e1ac42049 Historical Result
gernotbelger
parents: 9246
diff changeset
140 rows.add(row1);
600e1ac42049 Historical Result
gernotbelger
parents: 9246
diff changeset
141 result = new SalixLineCalculationHistoricalResult("Ergebnis 1 historical test", null, rows);
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
142 }
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
143 }
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
144 results.addResult(result, problems);
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
145 return new CalculationResult(results, problems);
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
146 }
9243
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
147
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
148 private SalixScenario findScenarioByKm(final double km, final Map<DoubleRange, SalixScenario> rangeScenarioMap) {
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
149 final Iterator<DoubleRange> rangeIterator = rangeScenarioMap.keySet().iterator();
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
150 while (rangeIterator.hasNext()) {
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
151 final DoubleRange range = rangeIterator.next();
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
152 if (range.containsDouble(km + 0.0001)) {
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
153 return rangeScenarioMap.get(range);
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
154 }
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
155 }
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
156 return null;
da0bdbcd6f09 salix.supraRegional+Regional results
gernotbelger
parents: 9222
diff changeset
157 }
8996
fb9430250899 Work on uinfo
gernotbelger
parents:
diff changeset
158 }

http://dive4elements.wald.intevation.org