|
8854
|
1 /* Copyright (C) 2017 by Bundesanstalt für Gewässerkunde |
|
8877
|
2 * Software engineering by |
|
|
3 * Björnsen Beratende Ingenieure GmbH |
|
8854
|
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.sinfo.flowdepth; |
|
|
11 |
|
|
12 import java.util.Collection; |
|
|
13 import java.util.List; |
|
|
14 |
|
8879
|
15 import org.dive4elements.artifacts.ArtifactDatabase; |
|
8854
|
16 import org.dive4elements.artifacts.CallContext; |
|
|
17 import org.dive4elements.river.artifacts.BedHeightsArtifact; |
|
|
18 import org.dive4elements.river.artifacts.model.Calculation; |
|
|
19 import org.dive4elements.river.artifacts.model.CalculationResult; |
|
|
20 import org.dive4elements.river.artifacts.model.LocationProvider; |
|
8877
|
21 import org.dive4elements.river.artifacts.model.QKms; |
|
8854
|
22 import org.dive4elements.river.artifacts.model.WKms; |
|
|
23 import org.dive4elements.river.artifacts.resources.Resources; |
|
|
24 import org.dive4elements.river.artifacts.sinfo.SINFOArtifact; |
|
|
25 import org.dive4elements.river.artifacts.sinfo.flowdepth.FlowDepthAccess.DifferencesPair; |
|
|
26 import org.dive4elements.river.artifacts.states.WDifferencesState; |
|
|
27 import org.dive4elements.river.model.BedHeight; |
|
|
28 import org.dive4elements.river.model.Gauge; |
|
|
29 import org.dive4elements.river.model.River; |
|
|
30 import org.dive4elements.river.utils.GaugeIndex; |
|
|
31 import org.dive4elements.river.utils.RiverUtils; |
|
|
32 |
|
|
33 class FlowDepthCalculation { |
|
|
34 |
|
|
35 private static final String CSV_NOT_IN_GAUGE_RANGE = "export.waterlevel.csv.not.in.gauge.range"; |
|
|
36 |
|
8877
|
37 private final CallContext context; |
|
8854
|
38 |
|
8877
|
39 public FlowDepthCalculation( final CallContext context ) { |
|
|
40 this.context = context; |
|
8854
|
41 } |
|
|
42 |
|
8877
|
43 public CalculationResult calculate(final SINFOArtifact sinfo) { |
|
8854
|
44 |
|
8879
|
45 /* |
|
|
46 * find the user of this artifact, sadly this is not part of the calling context, so instead we determine the |
|
|
47 * owner oft the artifact |
|
|
48 */ |
|
|
49 final ArtifactDatabase database = this.context.getDatabase(); |
|
|
50 final String user = database.findArtifactUser(sinfo.identifier()); |
|
8854
|
51 |
|
8877
|
52 /* access input data */ |
|
|
53 final FlowDepthAccess access = new FlowDepthAccess(sinfo); |
|
|
54 final River river = access.getRiver(); |
|
8854
|
55 |
|
8877
|
56 final Collection<DifferencesPair> diffPairs = access.getDifferencePairs(); |
|
|
57 |
|
|
58 final double from = access.getFrom(); |
|
|
59 final double to = access.getTo(); |
|
|
60 |
|
|
61 final boolean useTkh = access.isUseTransportBodies(); |
|
|
62 |
|
|
63 /* calculate results for each diff pair */ |
|
|
64 final Calculation problems = new Calculation(); |
|
|
65 |
|
|
66 final List<Gauge> gauges = river.determineGauges(from, to); |
|
|
67 final GaugeIndex gaugeIndex = new GaugeIndex(gauges); |
|
|
68 |
|
|
69 final String calcModeLabel = Resources.getMsg(this.context.getMeta(),sinfo.getCalculationMode().name() ); |
|
|
70 |
|
|
71 final FlowDepthCalculationResults results = new FlowDepthCalculationResults(calcModeLabel, user, river, from, to, useTkh); |
|
|
72 |
|
|
73 for (final DifferencesPair diffPair : diffPairs) { |
|
|
74 final FlowDepthCalculationResult result = calculateResult( river, from, to, diffPair, problems, gaugeIndex ); |
|
|
75 if( result != null ) |
|
|
76 results.addResult(result); |
|
|
77 } |
|
|
78 |
|
|
79 return new CalculationResult(results,problems); |
|
|
80 } |
|
|
81 |
|
|
82 private FlowDepthCalculationResult calculateResult(final River river, final double from, final double to, final DifferencesPair diffPair, final Calculation problems, final GaugeIndex gaugeIndex) { |
|
|
83 |
|
|
84 /* access real input data from database */ |
|
|
85 final String soundingId = diffPair.getSoundingId(); |
|
|
86 final String wstId = diffPair.getWstId(); |
|
|
87 |
|
|
88 final BedHeight bedHeight = loadBedHeight( soundingId, from, to ); |
|
|
89 if( bedHeight == null ) |
|
|
90 { |
|
|
91 final String message = Resources.format(this.context.getMeta(), "Failed to access sounding with id '{0}'", soundingId); |
|
|
92 problems.addProblem(message); |
|
|
93 return null; |
|
|
94 } |
|
|
95 |
|
|
96 final WKms wstKms = new WDifferencesState().getWKms(wstId, this.context, from, to); |
|
|
97 if( wstKms == null ) |
|
|
98 { |
|
|
99 final String message = Resources.format(this.context.getMeta(), "Failed to access waterlevel with id '{0}'", wstId); |
|
|
100 problems.addProblem(message); |
|
|
101 return null; |
|
|
102 } |
|
|
103 |
|
|
104 // FIXME: woher bekommen? |
|
|
105 final int wspYear = 0; |
|
8879
|
106 // FIXME: angeblich metadatum bestimmter wsps? |
|
|
107 // Könnte 'source' an 'wst_columns' sein |
|
|
108 final String wspSource = "FIXME"; |
|
|
109 // FIXME: Umsetzung IDENTISCH zu allen möglichen Arten wie ein WSPL berechnet wird.... |
|
|
110 final String wspGauge = "FIXME"; |
|
8877
|
111 |
|
|
112 final String wspLabel = wstKms.getName(); |
|
|
113 final String soundingLabel = bedHeight.getDescription(); |
|
|
114 final String label = String.format("%s - %s", wspLabel, soundingLabel); |
|
|
115 |
|
|
116 final BedHeightInfo sounding = BedHeightInfo.from(bedHeight); |
|
8879
|
117 final WstInfo wstInfo = new WstInfo(wspLabel, wspYear, wspSource, wspGauge); |
|
8877
|
118 |
|
|
119 final FlowDepthCalculationResult resultData = new FlowDepthCalculationResult(label, wstInfo, sounding); |
|
|
120 |
|
|
121 final String notinrange = Resources.getMsg(this.context.getMeta(), CSV_NOT_IN_GAUGE_RANGE, CSV_NOT_IN_GAUGE_RANGE); |
|
|
122 |
|
|
123 // TODO: prüfe diskretisierung wsp --> > 1000m --> Fehlermeldung |
|
|
124 |
|
|
125 // TODO: Berechnung der Transportkörperhöhen |
|
|
126 // - woher kommen die zusätzlichen eingangsdaten? sind das fixe daten pro gewässer? --> falls ja, warum nicht einmal berechnen und in db ablegen? |
|
|
127 |
|
|
128 // Benötigte Daten |
|
|
129 // - Abfluss / Station |
|
|
130 // - kein Abfluss --> Fehler |
|
|
131 if( !(wstKms instanceof QKms)) |
|
|
132 { |
|
|
133 final String message = Resources.format(this.context.getMeta(), "{0}: keine Abflussdaten vorhanden, Transportkörperhöhenberechnung nicht möglich", label); |
|
|
134 problems.addProblem(message); |
|
|
135 // TODO: keine Berechnung TKH |
|
|
136 } |
|
|
137 |
|
|
138 // - Sohlbeschaffenheit (D50 Korndurchmesser aus Seddb) |
|
|
139 // - Abhängig von Peiljahr |
|
|
140 // - kein D50 vorhanden --> Fehler |
|
|
141 // - Art der Gewässersohle (starr/mobil) |
|
|
142 |
|
|
143 |
|
|
144 final String bedHeightLabel = bedHeight.getDescription(); |
|
|
145 final String wstLabel = wstKms.getName(); |
|
|
146 |
|
|
147 // FIXME: basis der diskretisierung ist bedHeight, die wspl werden interpoliert |
|
|
148 for (int i = 0; i < wstKms.size(); i++) { |
|
|
149 |
|
|
150 final double km = wstKms.getKm(i); |
|
|
151 final double wst = wstKms.getW(i); |
|
|
152 // FIXME: interpolate from bedheights? |
|
|
153 final double meanBedHeight = 79.32; |
|
|
154 |
|
|
155 final double flowDepth = wst - meanBedHeight; |
|
|
156 |
|
|
157 final double discharge = wstKms instanceof QKms ? ((QKms) wstKms).getQ(i) : Double.NaN; |
|
|
158 |
|
|
159 // FIXME: calculate tkh |
|
|
160 final double tkh = 0; |
|
|
161 final double flowDepthTkh = flowDepth - tkh; |
|
|
162 |
|
|
163 // REMARK: access the location once only during calculation |
|
|
164 final String location = LocationProvider.getLocation(river.getName(), km); |
|
|
165 |
|
|
166 // REMARK: access the gauge once only during calculation |
|
|
167 // FIXME: copy specific handling from original wst |
|
|
168 final Gauge gauge = gaugeIndex.findGauge(km); |
|
|
169 final String gaugeLabel = gauge == null ? notinrange : gauge.getName(); |
|
|
170 |
|
|
171 resultData.addRow( km, flowDepth, flowDepthTkh, tkh, wst, discharge, wstLabel, gaugeLabel, meanBedHeight, bedHeightLabel, location ); |
|
|
172 } |
|
|
173 |
|
|
174 return resultData; |
|
|
175 } |
|
|
176 |
|
|
177 private BedHeight loadBedHeight(final String soundingId, final double from, final double to) { |
|
|
178 |
|
|
179 // FIXME: absolutely unbelievable.... |
|
|
180 // The way how bed-heights (and other data too) is accessed is different for nearly ever calculation-type throughout flys. |
|
|
181 // The knowledge on how to parse the datacage-ids is spread through the complete code-base... |
|
|
182 |
|
|
183 // We use here the way on how bed-heights are accessed by the BedDifferenceAccess/BedDifferenceCalculation, but this is plain random |
|
|
184 final String[] parts = soundingId.split(";"); |
|
|
185 |
|
|
186 final BedHeightsArtifact artifact = (BedHeightsArtifact) RiverUtils.getArtifact(parts[0], this.context); |
|
|
187 |
|
|
188 final Integer bedheightId = artifact.getDataAsInteger("height_id"); |
|
|
189 // FIXME: this only works with type 'single'; unclear on how to distinguish from epoch data (or whatever the other type means) |
|
|
190 // Luckily, the requirement is to only access 'single' data here. |
|
|
191 // final String bedheightType = artifact.getDataAsString("type"); |
|
|
192 |
|
|
193 // FIXME: BedDifferences uses this, but we also need the metadata of the BedHeight |
|
|
194 // FIXME: second absolutely awful thing: BedHeight is a hibernate binding class, accessing the database via hibernate stuff |
|
|
195 // BedHeightFactory uses its own (direct) way of accessing the data, with its own implemented data classes. |
|
|
196 //return BedHeightFactory.getHeight(bedheightType, bedheightId, from, to); |
|
|
197 |
|
|
198 return BedHeight.getBedHeightById(bedheightId); |
|
|
199 } |
|
8854
|
200 } |
