dive4elements/river: artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/BedQualityCalculation.java annotate

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/BedQualityCalculation.java @ 9801:1d7a72a50183 3.2.x tip

Assume Compose V2, consistently

author	Tom Gottfried <tom@intevation.de>
date	Thu, 23 Nov 2023 10:14:13 +0100
parents	0a5239a1e46e
children

rev	line source
5863 4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	1 /* Copyright (C) 2011, 2012, 2013 by Bundesanstalt für Gewässerkunde
4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	2 * Software engineering by Intevation GmbH
4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	3 *
5994 af13ceeba52a Removed trailing whitespace. Sascha L. Teichmann <teichmann@intevation.de> parents: 5863 diff changeset	4 * This file is Free Software under the GNU AGPL (>=v3)
5863 4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	5 * and comes with ABSOLUTELY NO WARRANTY! Check out the
5994 af13ceeba52a Removed trailing whitespace. Sascha L. Teichmann <teichmann@intevation.de> parents: 5863 diff changeset	6 * documentation coming with Dive4Elements River for details.
5863 4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	7 */
4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	8
5831 bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 4806 diff changeset	9 package org.dive4elements.river.artifacts.model.minfo;
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	10
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	11 import gnu.trove.TDoubleArrayList;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	12
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	13 import java.util.LinkedList;
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	14 import java.util.List;
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	15 import java.util.Map;
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	16
9726 0a5239a1e46e Upgrade to Log4j 2 Tom Gottfried <tom@intevation.de> parents: 8856 diff changeset	17 import org.apache.logging.log4j.Logger;
0a5239a1e46e Upgrade to Log4j 2 Tom Gottfried <tom@intevation.de> parents: 8856 diff changeset	18 import org.apache.logging.log4j.LogManager;
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	19
5831 bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 4806 diff changeset	20 import org.dive4elements.river.artifacts.access.BedQualityAccess;
bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 4806 diff changeset	21 import org.dive4elements.river.artifacts.model.Calculation;
bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 4806 diff changeset	22 import org.dive4elements.river.artifacts.model.CalculationResult;
bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 4806 diff changeset	23 import org.dive4elements.river.artifacts.model.DateRange;
bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 4806 diff changeset	24 import org.dive4elements.river.backend.SedDBSessionHolder;
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	25
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	26
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	27 public class BedQualityCalculation extends Calculation {
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	28
9726 0a5239a1e46e Upgrade to Log4j 2 Tom Gottfried <tom@intevation.de> parents: 8856 diff changeset	29 private static final Logger log = LogManager
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	30 .getLogger(BedQualityCalculation.class);
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	31
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	32 protected String river;
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	33 protected double from;
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	34 protected double to;
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	35 protected List<String> bedDiameter;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	36 protected List<String> bedloadDiameter;
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	37 protected List<DateRange> ranges;
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	38
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	39 public BedQualityCalculation() {
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	40 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	41
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	42 public CalculationResult calculate(BedQualityAccess access) {
8202 e4606eae8ea5 sed src/*/.java 's/logger/log/g' Sascha L. Teichmann <teichmann@intevation.de> parents: 7707 diff changeset	43 log.info("BedQualityCalculation.calculate");
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	44
7261 a56fe3bc6700 Refactoring: Let RiverAccess.getRiver return an River. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6257 diff changeset	45 String river = access.getRiverName();
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	46 Double from = access.getFrom();
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	47 Double to = access.getTo();
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	48 List<String> bedDiameter = access.getBedDiameter();
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	49 List<String> bedloadDiameter = access.getBedloadDiameter();
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	50 List<DateRange> ranges = access.getDateRanges();
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	51
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	52 if (river == null) {
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	53 // TODO: i18n
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	54 addProblem("minfo.missing.river");
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	55 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	56
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	57 if (from == null) {
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	58 // TODO: i18n
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	59 addProblem("minfo.missing.from");
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	60 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	61
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	62 if (to == null) {
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	63 // TODO: i18n
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	64 addProblem("minfo.missing.to");
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	65 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	66
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	67 if (ranges == null) {
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	68 // TODO: i18n
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	69 addProblem("minfo.missing.periods");
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	70 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	71
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	72 if (!hasProblems()) {
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	73 this.river = river;
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	74 this.from = from;
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	75 this.to = to;
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	76 this.ranges = ranges;
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	77 this.bedDiameter = bedDiameter;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	78 this.bedloadDiameter = bedloadDiameter;
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	79
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	80 SedDBSessionHolder.acquire();
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	81 try {
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	82 return internalCalculate();
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	83 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	84 finally {
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	85 SedDBSessionHolder.release();
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	86 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	87 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	88
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	89 return new CalculationResult();
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	90 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	91
8657 4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	92 /** Adds non empty values to a result and adds Problems for empty ones.*/
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	93 protected void addValuesToResult(BedQualityResult result,
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	94 BedQualityResultValue[] values) {
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	95 for (BedQualityResultValue value: values) {
8686 9cbbaedc6cba (issue1755) Remove redundant conditional checks. Tom Gottfried <tom@intevation.de> parents: 8668 diff changeset	96 if (!value.isInterpolateable()) {
8657 4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	97 if (value.isDiameterResult()) {
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	98 addProblem("bedquality.missing.diameter." +
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	99 value.getType(), value.getName().toUpperCase(),
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	100 result.getDateRange().getFrom(),
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	101 result.getDateRange().getTo());
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	102 } else {
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	103 addProblem("bedquality.missing." + value.getName() + "." +
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	104 value.getType(), result.getDateRange().getFrom(),
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	105 result.getDateRange().getTo());
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	106 }
8686 9cbbaedc6cba (issue1755) Remove redundant conditional checks. Tom Gottfried <tom@intevation.de> parents: 8668 diff changeset	107 if (!value.isEmpty() && !value.isNaN()) {
8668 fc5fcec91b73 (issue1755) Add Problem also if a result is not interpolateable Andre Heinecke <andre.heinecke@intevation.de> parents: 8667 diff changeset	108 // we want to keep single point results
fc5fcec91b73 (issue1755) Add Problem also if a result is not interpolateable Andre Heinecke <andre.heinecke@intevation.de> parents: 8667 diff changeset	109 result.add(value);
fc5fcec91b73 (issue1755) Add Problem also if a result is not interpolateable Andre Heinecke <andre.heinecke@intevation.de> parents: 8667 diff changeset	110 }
8657 4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	111 } else {
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	112 result.add(value);
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	113 }
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	114 }
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	115 }
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	116
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	117
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	118 protected CalculationResult internalCalculate() {
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	119
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	120 List<BedQualityResult> results = new LinkedList<BedQualityResult>();
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	121 // Calculate for all time periods.
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	122 for (DateRange dr : ranges) {
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	123 QualityMeasurements loadMeasurements =
3785 a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3784 diff changeset	124 QualityMeasurementFactory.getBedloadMeasurements(
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	125 river,
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	126 from,
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	127 to,
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	128 dr.getFrom(),
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	129 dr.getTo());
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	130 QualityMeasurements bedMeasurements =
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	131 QualityMeasurementFactory.getBedMeasurements(
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	132 river,
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	133 from,
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	134 to,
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	135 dr.getFrom(),
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	136 dr.getTo());
3784 9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	137 BedQualityResult result = new BedQualityResult();
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	138 result.setDateRange(dr);
8592 36faef4f8acb (issue1755) Add interpolated facets and update themes accordingly Andre Heinecke <andre.heinecke@intevation.de> parents: 8587 diff changeset	139 if (!bedDiameter.isEmpty()) {
36faef4f8acb (issue1755) Add interpolated facets and update themes accordingly Andre Heinecke <andre.heinecke@intevation.de> parents: 8587 diff changeset	140 log.debug("Bed diameter is not empty + " + bedDiameter);
8856 5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	141 addValuesToResult(
5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	142 result, calculateBedParameter(bedMeasurements));
3780 7fa38f8bcd8d Avoid NPE if the user selects onlyone sort of characteristic diameter. Raimund Renkert <raimund.renkert@intevation.de> parents: 3772 diff changeset	143 for (String bd : bedDiameter) {
8856 5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	144 addValuesToResult(
5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	145 result, calculateBed(bedMeasurements, bd));
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	146 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	147 }
8592 36faef4f8acb (issue1755) Add interpolated facets and update themes accordingly Andre Heinecke <andre.heinecke@intevation.de> parents: 8587 diff changeset	148 if (!bedloadDiameter.isEmpty()) {
3780 7fa38f8bcd8d Avoid NPE if the user selects onlyone sort of characteristic diameter. Raimund Renkert <raimund.renkert@intevation.de> parents: 3772 diff changeset	149 for (String bld : bedloadDiameter) {
8856 5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	150 addValuesToResult(
5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	151 result, calculateBedload(loadMeasurements, bld));
3780 7fa38f8bcd8d Avoid NPE if the user selects onlyone sort of characteristic diameter. Raimund Renkert <raimund.renkert@intevation.de> parents: 3772 diff changeset	152 }
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	153 }
3784 9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	154 results.add(result);
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	155 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	156
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	157 return new CalculationResult(
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	158 results.toArray(new BedQualityResult[results.size()]), this);
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	159 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	160
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	161 private BedQualityResultValue[] calculateBedParameter(
07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	162 QualityMeasurements qm
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	163 ) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	164 List<Double> kms = qm.getKms();
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	165 TDoubleArrayList location = new TDoubleArrayList();
3784 9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	166 QualityMeasurements capFiltered = filterCapMeasurements(qm);
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	167 QualityMeasurements subFiltered = filterSubMeasurements(qm);
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	168 TDoubleArrayList porosityCap = new TDoubleArrayList();
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	169 TDoubleArrayList porositySub = new TDoubleArrayList();
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	170 TDoubleArrayList densityCap = new TDoubleArrayList();
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	171 TDoubleArrayList densitySub = new TDoubleArrayList();
3785 a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3784 diff changeset	172
3784 9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	173 for(double km : kms) {
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	174 double[] pCap = calculatePorosity(capFiltered, km);
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	175 double[] pSub = calculatePorosity(subFiltered, km);
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	176 double[] dCap = calculateDensity(capFiltered, pCap);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	177 double[] dSub = calculateDensity(subFiltered, pSub);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	178
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	179 double pCapRes = 0d;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	180 double pSubRes = 0d;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	181 double dCapRes = 0d;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	182 double dSubRes = 0d;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	183 for (int i = 0; i < pCap.length; i++) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	184 pCapRes += pCap[i];
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	185 dCapRes += dCap[i];
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	186 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	187 for (int i = 0; i < pSub.length; i++) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	188 pSubRes += pSub[i];
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	189 dSubRes += dSub[i];
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	190 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	191 location.add(km);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	192 porosityCap.add((pCapRes / pCap.length) * 100 );
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	193 porositySub.add((pSubRes / pSub.length) * 100);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	194 densityCap.add((dCapRes / dCap.length) / 1000);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	195 densitySub.add((dSubRes / dSub.length) / 1000);
3785 a5f65e8983be Merged revisions 5501-5502,5504-5508,5511-5513,5516-5519 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3784 diff changeset	196
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	197 }
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	198 return new BedQualityResultValue[] {
07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	199 new BedQualityResultValue("porosity",
8856 5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	200 new double[][] {location.toNativeArray(),
5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	201 porositySub.toNativeArray()},
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	202 "sublayer"),
07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	203 new BedQualityResultValue("porosity",
8856 5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	204 new double[][] {location.toNativeArray(),
5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	205 porosityCap.toNativeArray()},
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	206 "toplayer"),
07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	207 new BedQualityResultValue("density",
8856 5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	208 new double[][] {location.toNativeArray(),
5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	209 densitySub.toNativeArray()},
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	210 "sublayer"),
07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	211 new BedQualityResultValue("density",
8856 5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	212 new double[][] {location.toNativeArray(),
5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	213 densityCap.toNativeArray()},
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	214 "toplayer")};
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	215 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	216
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	217 protected BedQualityResultValue[] calculateBed(
3784 9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	218 QualityMeasurements qm,
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	219 String diameter
3784 9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	220 ) {
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	221 List<Double> kms = qm.getKms();
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	222 TDoubleArrayList location = new TDoubleArrayList();
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	223 TDoubleArrayList avDiameterCap = new TDoubleArrayList();
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	224 TDoubleArrayList avDiameterSub = new TDoubleArrayList();
8563 59b4f6cf1cf4 Move static operation out of loop. Andre Heinecke <andre.heinecke@intevation.de> parents: 8202 diff changeset	225 QualityMeasurements capFiltered = filterCapMeasurements(qm);
59b4f6cf1cf4 Move static operation out of loop. Andre Heinecke <andre.heinecke@intevation.de> parents: 8202 diff changeset	226 QualityMeasurements subFiltered = filterSubMeasurements(qm);
59b4f6cf1cf4 Move static operation out of loop. Andre Heinecke <andre.heinecke@intevation.de> parents: 8202 diff changeset	227
3784 9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	228 for (double km : kms) {
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	229 //Filter cap and sub measurements.
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	230 List<QualityMeasurement> cm = capFiltered.getMeasurements(km);
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	231 List<QualityMeasurement> sm = subFiltered.getMeasurements(km);
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	232
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	233 double avCap = calculateAverage(cm, diameter);
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	234 double avSub = calculateAverage(sm, diameter);
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	235 location.add(km);
4642 ca18d8d06218 Multiplicate diameter with 1000 to have mm instead of m. Raimund Renkert <rrenkert@intevation.de> parents: 3974 diff changeset	236 avDiameterCap.add(avCap * 1000);// bring to mm.
ca18d8d06218 Multiplicate diameter with 1000 to have mm instead of m. Raimund Renkert <rrenkert@intevation.de> parents: 3974 diff changeset	237 avDiameterSub.add(avSub * 1000);
3784 9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	238 }
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	239 return new BedQualityResultValue[] {
07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	240 new BedQualityResultValue(diameter,
8856 5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	241 new double[][] {location.toNativeArray(),
5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	242 avDiameterSub.toNativeArray()},
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	243 "sublayer"),
07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	244 new BedQualityResultValue(diameter,
8856 5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	245 new double[][] {location.toNativeArray(),
5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	246 avDiameterCap.toNativeArray()},
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	247 "toplayer")};
3784 9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	248 }
9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	249
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	250 private double[] calculateDensity(
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	251 QualityMeasurements capFiltered,
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	252 double[] porosity
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	253 ) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	254 double[] density = new double[porosity.length];
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	255 for (int i = 0; i < porosity.length; i++) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	256 density[i] = (1 - porosity[i]) * 2650;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	257 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	258 return density;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	259 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	260
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	261 private double[] calculatePorosity(
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	262 QualityMeasurements capFiltered,
3784 9f9d5ada96ca Merged revisions 5498-5499 via svnmerge from Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 3780 diff changeset	263 double km
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	264 ) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	265 List<QualityMeasurement> list = capFiltered.getMeasurements(km);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	266 double[] results = new double[list.size()];
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	267 int i = 0;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	268 for (QualityMeasurement qm : list) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	269 double deviation = calculateDeviation(qm);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	270 double p = calculateP(qm);
7670 03dc5cc22e26 Update coefficients in porosity formula. Tom Gottfried <tom.gottfried@intevation.de> parents: 7261 diff changeset	271 double porosity = 0.362 - 0.088 * deviation + 0.219 * p;
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	272 results[i] = porosity;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	273 i++;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	274 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	275
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	276 return results;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	277 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	278
8657 4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	279 protected BedQualityResultValue[] calculateBedload(
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	280 QualityMeasurements qm,
8587 07c9ac22f611 (issue1755) Generalise BedQuality result handling Andre Heinecke <andre.heinecke@intevation.de> parents: 8563 diff changeset	281 String diameter
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	282 ) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	283 List<Double> kms = qm.getKms();
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	284 TDoubleArrayList location = new TDoubleArrayList();
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	285 TDoubleArrayList avDiameter = new TDoubleArrayList();
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	286 for (double km : kms) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	287 List<QualityMeasurement> measurements = qm.getMeasurements(km);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	288 double mid = calculateAverage(measurements, diameter);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	289 location.add(km);
6257 165f60f420cf Bring bedload average diameter to 'mm'. Raimund Renkert <rrenkert@intevation.de> parents: 6151 diff changeset	290 avDiameter.add(mid * 1000);
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	291 }
8657 4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	292 return new BedQualityResultValue[] {
4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	293 new BedQualityResultValue(diameter,
8856 5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	294 new double[][] {location.toNativeArray(),
5e38e2924c07 Fix code style. Tom Gottfried <tom@intevation.de> parents: 8686 diff changeset	295 avDiameter.toNativeArray()},
8657 4054c7fce43e (issue1755) Add problem / report handling to BedQualityCalculation Andre Heinecke <andre.heinecke@intevation.de> parents: 8592 diff changeset	296 "bedload")};
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	297 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	298
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	299 protected double calculateAverage(
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	300 List<QualityMeasurement> list,
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	301 String diameter
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	302 ) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	303 double av = 0;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	304 for (QualityMeasurement qm : list) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	305 av += qm.getDiameter(diameter);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	306 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	307 return av/list.size();
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	308 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	309
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	310 protected QualityMeasurements filterCapMeasurements(
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	311 QualityMeasurements qms
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	312 ) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	313 List<QualityMeasurement> result = new LinkedList<QualityMeasurement>();
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	314 for (QualityMeasurement qm : qms.getMeasurements()) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	315 if (qm.getDepth1() == 0d && qm.getDepth2() <= 0.3) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	316 result.add(qm);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	317 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	318 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	319 return new QualityMeasurements(result);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	320 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	321
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	322 protected QualityMeasurements filterSubMeasurements(
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	323 QualityMeasurements qms
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	324 ) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	325 List<QualityMeasurement> result = new LinkedList<QualityMeasurement>();
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	326 for (QualityMeasurement qm : qms.getMeasurements()) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	327 if (qm.getDepth1() > 0d && qm.getDepth2() <= 0.5) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	328 result.add(qm);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	329 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	330 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	331 return new QualityMeasurements(result);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	332 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	333
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	334 public double calculateDeviation(QualityMeasurement qm) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	335 Map<String, Double> dm = qm.getAllDiameter();
6151 63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	336 int size = dm.size();
63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	337
63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	338 double phiM = 0;
63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	339 double [] phis = new double[size];
63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	340 double [] ps = new double[size];
63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	341 double scale = -1d/Math.log(2d);
63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	342
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	343 int i = 0;
6151 63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	344 for (Map.Entry<String, Double> entry: dm.entrySet()) {
63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	345 double phi = scale*Math.log(entry.getValue());
63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	346 double p = calculateWeight(qm, entry.getKey());
63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	347 phiM += phi * p;
63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	348 ps[i] = p;
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	349 phis[i] = phi;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	350 i++;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	351 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	352
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	353 double sig = 0d;
6151 63d1e2a9b311 Bed quality calculation: Replaced another iteration over the key set + directly fetching the values withan iteration over the entry set. Sascha L. Teichmann <teichmann@intevation.de> parents: 6150 diff changeset	354 for (i = 0; i < size; i++) {
7707 cead6df5b7b3 Variance is defined as the sum of squared differences. Tom Gottfried <tom.gottfried@intevation.de> parents: 7670 diff changeset	355 sig += ps[i] * Math.pow((phis[i] - phiM), 2);
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	356 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	357 double deviation = Math.sqrt(sig);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	358 return deviation;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	359 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	360
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	361 protected double calculateP(QualityMeasurement qm) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	362 return calculateWeight(qm, "dmin");
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	363 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	364
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	365 public double calculateWeight(QualityMeasurement qm, String diameter) {
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	366 Map<String, Double> dm = qm.getAllDiameter();
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	367 double value = qm.getDiameter(diameter);
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	368
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	369 double sum = 0d;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	370 for (Double d : dm.values()) {
6150 1f057463765b Bed quality calculation: Fixed a calculation typo: s/sum =+ d/sum += d/ Sascha L. Teichmann <teichmann@intevation.de> parents: 5994 diff changeset	371 sum += d.doubleValue();
3769 728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	372 }
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	373 double weight = sum/100*value;
728ecd2afa20 Implemented bed quality calculation in minfo module. Raimund Renkert <raimund.renkert@intevation.de> parents: 3760 diff changeset	374 return weight;
3760 312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	375 }
312870fded7e Added result objects for minfo bed quality calculation. Raimund Renkert <raimund.renkert@intevation.de> parents: diff changeset	376 }
4806 750d72a16f4e BedQualityCalculation: Added vim-line. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 4642 diff changeset	377 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :

Mercurial > dive4elements > river

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/BedQualityCalculation.java @ 9801:1d7a72a50183 3.2.x tip