dive4elements/river: artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/flood_duration/FloodDurationCalculator.java annotate

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/flood_duration/FloodDurationCalculator.java @ 9202:b4402594213b

More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2)

author	mschaefer
date	Mon, 02 Jul 2018 07:33:53 +0200
parents	1614cb14308f
children	3dae6b78e1da

rev	line source
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	1 /** Copyright (C) 2017 by Bundesanstalt für Gewässerkunde
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	2 * Software engineering by
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	3 * Björnsen Beratende Ingenieure GmbH
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	4 * Dr. Schumacher Ingenieurbüro für Wasser und Umwelt
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	5 *
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	6 * This file is Free Software under the GNU AGPL (>=v3)
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	7 * and comes with ABSOLUTELY NO WARRANTY! Check out the
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	8 * documentation coming with Dive4Elements River for details.
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	9 */
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	10 package org.dive4elements.river.artifacts.sinfo.flood_duration;
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	11
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	12 import java.util.ArrayList;
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	13 import java.util.Collection;
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	14 import java.util.HashMap;
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	15 import java.util.List;
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	16 import java.util.Map;
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	17 import java.util.Set;
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	18
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	19 import org.apache.commons.lang.math.DoubleRange;
9150 23945061daec gigantic refactoring: exporter, result, results gernotbelger parents: 9145 diff changeset	20 import org.dive4elements.artifacts.CallContext;
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	21 import org.dive4elements.river.artifacts.WINFOArtifact;
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	22 import org.dive4elements.river.artifacts.access.ComputationRangeAccess;
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	23 import org.dive4elements.river.artifacts.common.GeneralResultType;
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	24 import org.dive4elements.river.artifacts.common.ResultRow;
9150 23945061daec gigantic refactoring: exporter, result, results gernotbelger parents: 9145 diff changeset	25 import org.dive4elements.river.artifacts.model.Calculation;
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	26 import org.dive4elements.river.artifacts.model.Calculation.Problem;
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	27 import org.dive4elements.river.artifacts.model.CalculationResult;
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	28 import org.dive4elements.river.artifacts.model.WQKms;
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	29 import org.dive4elements.river.artifacts.sinfo.common.GaugeDurationValuesFinder;
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	30 import org.dive4elements.river.artifacts.sinfo.common.GaugeMainValueFinder;
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	31 import org.dive4elements.river.artifacts.sinfo.common.RiverInfoProvider;
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	32 import org.dive4elements.river.artifacts.sinfo.common.SInfoResultType;
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	33 import org.dive4elements.river.artifacts.sinfo.common.WQBaseTableFinder;
1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	34 import org.dive4elements.river.artifacts.sinfo.flood_duration.RiversideRadioChoice.RiversideChoiceKey;
1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	35 import org.dive4elements.river.model.Attribute.AttributeKey;
1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	36 import org.dive4elements.river.model.Gauge;
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	37 import org.dive4elements.river.model.MainValueType.MainValueTypeKey;
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	38 import org.dive4elements.river.model.sinfo.InfrastructureValue;
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	39
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	40 import gnu.trove.TDoubleArrayList;
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	41
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	42 /**
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	43 * Calculation of the result rows of the flood duration of the infrastructures in a river km range
1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	44 * and selected main value durations
1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	45 *
1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	46 * @author Matthias Schäfer
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	47 */
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	48 final class FloodDurationCalculator {
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	49
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	50 private final Collection<ResultRow> rows = new ArrayList<>();
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	51
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	52 private final RiverInfoProvider riverInfoProvider;
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	53
9150 23945061daec gigantic refactoring: exporter, result, results gernotbelger parents: 9145 diff changeset	54 private final CallContext context;
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	55
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	56
9150 23945061daec gigantic refactoring: exporter, result, results gernotbelger parents: 9145 diff changeset	57 public FloodDurationCalculator(final CallContext context, final RiverInfoProvider riverInfoProvider) {
23945061daec gigantic refactoring: exporter, result, results gernotbelger parents: 9145 diff changeset	58 this.context = context;
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	59 this.riverInfoProvider = riverInfoProvider;
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	60 }
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	61
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	62 /**
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	63 * Calculate the infrastructures flood duration result rows
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	64 */
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	65 public FloodDurationCalculationResult execute(final Calculation problems, final String label, final DoubleRange calcRange,
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	66 final RiversideChoiceKey riverside, final WINFOArtifact winfo) {
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	67
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	68 // Find all gauges of the calc range, and create the duration finders
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	69 final Map<Gauge, GaugeDurationValuesFinder> durFinders = new HashMap<>();
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	70 Gauge firstGauge = null;
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	71 for (final Gauge gauge : this.riverInfoProvider.getRiver().determineGauges(calcRange.getMinimumDouble(), calcRange.getMaximumDouble())) {
1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	72 durFinders.put(gauge, GaugeDurationValuesFinder.loadValues(gauge, problems));
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	73 if (firstGauge == null)
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	74 firstGauge = gauge;
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	75 }
9150 23945061daec gigantic refactoring: exporter, result, results gernotbelger parents: 9145 diff changeset	76
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	77 // Find all infrastructures within the calc range
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	78 final AttributeKey bankKey = riverside.getAttributeKey();
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	79 final List<InfrastructureValue> infras = InfrastructureValue.getValues(this.riverInfoProvider.getRiver(), calcRange.getMinimumDouble(),
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	80 calcRange.getMaximumDouble(), bankKey);
9150 23945061daec gigantic refactoring: exporter, result, results gernotbelger parents: 9145 diff changeset	81
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	82 // Merge all stations (range/step, borders of gauge ranges, infrastructures)
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	83 final Map<Double, InfrastructureValue> allStations = new HashMap<>();
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	84 final Map<Double, InfrastructureValue> secondBank = new HashMap<>(); // any second infrastructure in case of both-banks-option
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	85 addRangeStations(allStations, winfo);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	86 addGaugeLimits(allStations, durFinders.keySet(), calcRange.getMinimumDouble(), calcRange.getMaximumDouble());
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	87 addInfrastructures(allStations, secondBank, infras);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	88 final double[] stationsSorted = sortStations(allStations.keySet());
9150 23945061daec gigantic refactoring: exporter, result, results gernotbelger parents: 9145 diff changeset	89
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	90 // Calculate W and Q for all stations and the selected discharge states
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	91 final WQKms[] wqkmsArray = calculateWaterlevels(winfo, stationsSorted, problems);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	92
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	93 // Determine discharge state labels of the main values
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	94 final String[] mainValueLabels = findMainValueLabels(wqkmsArray, winfo.getQs(), firstGauge, problems);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	95
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	96 // Create a finder for Q in the {river}.wst km-w-q table
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	97 final WQBaseTableFinder wqFinder = WQBaseTableFinder.loadValues(this.riverInfoProvider.getRiver(), calcRange.getMinimumDouble(),
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	98 calcRange.getMaximumDouble(), problems);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	99
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	100 // Calculate the durations and create the result rows
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	101 for (int i = 0; i <= stationsSorted.length - 1; i++) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	102 final Gauge gauge = this.riverInfoProvider.getGauge(stationsSorted[i], true);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	103 final ResultRow row = createRow(stationsSorted[i], gauge, wqkmsArray, durFinders.get(gauge), i);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	104 if (allStations.containsKey(stationsSorted[i]) && (allStations.get(stationsSorted[i]) != null))
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	105 calculateInfrastructure(row, gauge, allStations.get(stationsSorted[i]), wqFinder, durFinders);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	106 this.rows.add(row);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	107 if (secondBank.containsKey(stationsSorted[i])) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	108 final ResultRow row2 = ResultRow.create(row);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	109 calculateInfrastructure(row2, gauge, secondBank.get(stationsSorted[i]), wqFinder, durFinders);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	110 this.rows.add(row2);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	111 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	112 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	113
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	114 return new FloodDurationCalculationResult(label, mainValueLabels, this.rows);
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	115 }
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	116
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	117 /**
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	118 * Adds to a stations map all stations corresponding to the active range and step
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	119 */
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	120 private void addRangeStations(final Map<Double, InfrastructureValue> allStations, final WINFOArtifact winfo) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	121 for (final double station : new ComputationRangeAccess(winfo).getKms())
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	122 allStations.put(Double.valueOf(station), null);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	123 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	124
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	125 /**
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	126 * Adds to a stations map all range limits of the gauges within the calc range
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	127 */
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	128 private void addGaugeLimits(final Map<Double, InfrastructureValue> allStations, final Set<Gauge> gauges, final double fromKm, final double toKm) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	129 for (final Gauge gauge : gauges) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	130 final Double kmA = Double.valueOf(gauge.getRange().getA().doubleValue());
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	131 final Double kmB = Double.valueOf(gauge.getRange().getB().doubleValue());
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	132 if (kmA > fromKm - 0.0001)
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	133 allStations.put(kmA, null);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	134 if (kmB < toKm + 0.0001)
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	135 allStations.put(kmB, null);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	136 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	137 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	138
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	139 /**
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	140 * Adds to a stations map all (first) infrastructures of a station, and the second, if any, to another map
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	141 */
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	142 private void addInfrastructures(final Map<Double, InfrastructureValue> allStations, final Map<Double, InfrastructureValue> secondBank,
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	143 final List<InfrastructureValue> infrastructures) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	144 for (final InfrastructureValue infrastructure : infrastructures) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	145 final Double station = infrastructure.getStation();
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	146 if (!allStations.containsKey(station) \|\| !(allStations.get(station) instanceof InfrastructureValue))
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	147 allStations.put(station, infrastructure);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	148 else
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	149 secondBank.put(station, infrastructure);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	150 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	151 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	152
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	153 /**
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	154 * Returns a double array with a sorted stations set
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	155 */
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	156 private double[] sortStations(final Set<Double> stations) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	157 final TDoubleArrayList sorted = new TDoubleArrayList();
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	158 for (final Double station : stations)
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	159 sorted.add(station.doubleValue());
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	160 sorted.sort();
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	161 return sorted.toNativeArray();
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	162 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	163
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	164 /**
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	165 * Calculates an array of w-q-longitudinal sections for all artifact W/Q options
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	166 */
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	167 private WQKms[] calculateWaterlevels(final WINFOArtifact winfo, final double[] stations, final Calculation problems) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	168 // REMARK aus TkhCalculation - move to WinfoArtifactWrapper?
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	169 // TODO das ist ziemlich langsam - durch den WQBaseTableFinder ersetzen? (vorher W-Optionen in Q umrechnen)
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	170 // (So funktioniert computeWaterlevelData wohl:
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	171 // Es sucht die Spalte(n) zum Bezugspegel-Q in der W-Q-Tabelle ({river}.wst in Wst etc.)
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	172 // und interpoliert für diese horizontale Tabellenposition jeweils die vertikale Tabellenposition der station;
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	173 // das ergibt das W einer station für einen Abflusszustand;
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	174 // bei Vorgabe eines Pegel-W wird vorher anhand der W-Q-Tabelle des Pegels ({gauge}.at in DischargeTable) das Q
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	175 // interpoliert;
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	176 // bei Vorgabe eines W auf freier Strecke wird wohl vorher noch die .wst-Interpolation eingesetzt.
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	177 final CalculationResult waterlevelData = winfo.computeWaterlevelData(stations);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	178
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	179 /* copy all problems */
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	180 final Calculation winfoProblems = waterlevelData.getReport();
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	181 final List<Problem> problems2 = winfoProblems.getProblems();
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	182 if (problems2 != null) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	183 for (final Problem problem : problems2) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	184 problems.addProblem(problem);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	185 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	186 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	187 return (WQKms[]) waterlevelData.getData();
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	188 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	189
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	190 /**
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	191 * Determines the discharge state labels for the selected Q or W values
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	192 */
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	193 private String[] findMainValueLabels(final WQKms[] wqkmsArray, final double[] qs, final Gauge gauge, final Calculation problems) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	194 final String[] mainValueLabels = new String[wqkmsArray.length];
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	195 if (wqkmsArray.length >= 1) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	196 // Labels like Q=123 or W=123
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	197 for (int i = 0; i <= wqkmsArray.length - 1; i++)
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	198 mainValueLabels[i] = wqkmsArray[i].getName();
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	199 // Replace labels for named main Q values
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	200 final GaugeMainValueFinder zoneFinder = GaugeMainValueFinder.loadValues(MainValueTypeKey.Q, gauge, problems);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	201 if ((zoneFinder != null) && (qs != null)) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	202 for (int i = 0; i <= qs.length - 1; i++)
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	203 mainValueLabels[i] = zoneFinder.findExactZoneName(qs[i], mainValueLabels[i]);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	204 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	205 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	206 return mainValueLabels;
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	207 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	208
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	209 /**
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	210 * Create a result row for a station and its gauge, and add w-q-values as selected
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	211 */
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	212 private ResultRow createRow(final Double station, final Gauge gauge, final WQKms[] wqkmsArray,
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	213 final GaugeDurationValuesFinder durationFinder, final int kmIndex) {
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	214
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	215 final ResultRow row = ResultRow.create();
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	216 row.putValue(GeneralResultType.station, station);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	217 row.putValue(SInfoResultType.infrastructuretype, null); // is replaced later for an infrastructure
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	218 row.putValue(SInfoResultType.floodDuration, Double.NaN); // is replaced later for an infrastructure
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	219 row.putValue(SInfoResultType.gaugeLabel, gauge.getName());
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	220 final String location = this.riverInfoProvider.getLocation(station);
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	221 row.putValue(SInfoResultType.location, location);
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	222
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	223 if (wqkmsArray.length >= 1) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	224 assert (wqkmsArray[0].getKm(kmIndex) == station.doubleValue());
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	225 row.putValue(SInfoResultType.waterlevel1, wqkmsArray[0].getW(kmIndex));
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	226 row.putValue(SInfoResultType.discharge1, wqkmsArray[0].getQ(kmIndex));
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	227 row.putValue(SInfoResultType.mainValue1Duration, underflowDaysToOverflowDays(durationFinder.getDuration(wqkmsArray[0].getQ(kmIndex))));
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	228 if (wqkmsArray.length >= 2) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	229 assert (wqkmsArray[1].getKm(kmIndex) == station.doubleValue());
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	230 row.putValue(SInfoResultType.waterlevel2, wqkmsArray[1].getW(kmIndex));
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	231 row.putValue(SInfoResultType.discharge2, wqkmsArray[1].getQ(kmIndex));
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	232 row.putValue(SInfoResultType.mainValue2Duration, underflowDaysToOverflowDays(durationFinder.getDuration(wqkmsArray[1].getQ(kmIndex))));
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	233 if (wqkmsArray.length >= 3) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	234 assert (wqkmsArray[2].getKm(kmIndex) == station.doubleValue());
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	235 row.putValue(SInfoResultType.waterlevel3, wqkmsArray[2].getW(kmIndex));
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	236 row.putValue(SInfoResultType.discharge3, wqkmsArray[2].getQ(kmIndex));
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	237 row.putValue(SInfoResultType.mainValue3Duration, underflowDaysToOverflowDays(durationFinder.getDuration(wqkmsArray[2].getQ(kmIndex))));
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	238 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	239 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	240 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	241 return row;
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	242 }
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	243
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	244 /**
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	245 * Calculate the result row fields for one infrastructure
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	246 */
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	247 private void calculateInfrastructure(final ResultRow row, final Gauge gauge, final InfrastructureValue infrastructure,
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	248 final WQBaseTableFinder wqFinder, final Map<Gauge, GaugeDurationValuesFinder> durFinders) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	249
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	250 final double q = wqFinder.getDischarge(infrastructure.getStation(), infrastructure.getHeight());
1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	251 final double qOut = Double.isInfinite(q) ? Double.NaN : q;
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	252 final double dur = underflowDaysToOverflowDays(durFinders.get(gauge).getDuration(q));
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	253 row.putValue(SInfoResultType.riverside, infrastructure.getAttributeKey().getName()); // TODO i18n
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	254 row.putValue(SInfoResultType.floodDuration, dur);
9176 1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	255 row.putValue(SInfoResultType.floodDischarge, qOut);
1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	256 row.putValue(SInfoResultType.infrastructureHeight, infrastructure.getHeight());
1614cb14308f Work on calculations for S-Info flood duration workflow mschaefer parents: 9170 diff changeset	257 row.putValue(SInfoResultType.infrastructuretype, infrastructure.getInfrastructure().getType().getName());
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	258 }
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	259
9202 b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	260 /**
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	261 * Translates underflow duration into overflow duration
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	262 */
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	263 private double underflowDaysToOverflowDays(final double underflowDays) {
b4402594213b More work on calculations and output for S-Info flood duration workflow (chart types 1 and 2) mschaefer parents: 9176 diff changeset	264 return 365 - underflowDays;
9145 e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	265 }
e6b63b2b41b9 sinfo.flood_duration pdf, csv, ui gernotbelger parents: diff changeset	266 }

Mercurial > dive4elements > river

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/flood_duration/FloodDurationCalculator.java @ 9202:b4402594213b