dive4elements/river: artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/tkhcalculation/BedQualityD50KmValueFinder.java annotate

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/tkhcalculation/BedQualityD50KmValueFinder.java @ 8915:d9dbf0b74bc2

Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation.

author	gernotbelger
date	Wed, 28 Feb 2018 17:27:15 +0100
parents	artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/flowdepth/BedQualityD50KmValueFinder.java@37ff7f435912
children	29442c03c6e3

rev	line source
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	1 /* Copyright (C) 2017 by Bundesanstalt für Gewässerkunde
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	2 * Software engineering by
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	3 * Björnsen Beratende Ingenieure GmbH
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	4 * Dr. Schumacher Ingenieurbüro für Wasser und Umwelt
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	5 *
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	6 * This file is Free Software under the GNU AGPL (>=v3)
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	7 * and comes with ABSOLUTELY NO WARRANTY! Check out the
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	8 * documentation coming with Dive4Elements River for details.
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	9 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	10
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	11 package org.dive4elements.river.artifacts.sinfo.tkhcalculation;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	12
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	13 import java.util.Calendar;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	14 import java.util.Date;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	15 import java.util.List;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	16
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	17 import org.apache.commons.lang.math.DoubleRange;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	18 import org.apache.commons.math.ArgumentOutsideDomainException;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	19 import org.apache.commons.math.analysis.interpolation.LinearInterpolator;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	20 import org.apache.commons.math.analysis.polynomials.PolynomialSplineFunction;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	21 import org.apache.log4j.Logger;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	22 import org.dive4elements.river.backend.SedDBSessionHolder;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	23 import org.dive4elements.river.model.River;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	24 import org.hibernate.SQLQuery;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	25 import org.hibernate.Session;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	26 import org.hibernate.type.StandardBasicTypes;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	27
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	28 /**
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	29 * Searchable sorted km array with parallel bed measurements value array and linear interpolation for km and d50 between
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	30 * the array elements.<br />
8901 0a900d605d52 S-INFO Flowdepth work on TKH calculation mschaefer parents: 8898 diff changeset	31 * <br />
0a900d605d52 S-INFO Flowdepth work on TKH calculation mschaefer parents: 8898 diff changeset	32 * See comment of SQL command on how the values are filtered and aggregated.
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	33 *
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	34 * @author Matthias Schäfer
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	35 *
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	36 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	37 final class BedQualityD50KmValueFinder {
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	38
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	39 /*** INNER CLASSES ***/
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	40
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	41 /**
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	42 * A bed measurements aggregate with its d50 characteristic grain diameter
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	43 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	44 private static class D50Measurement {
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	45 private double km;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	46
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	47 public double getKm() {
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	48 return this.km;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	49 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	50
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	51 // private Date mindate;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	52
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	53 // public Date getMinDate() {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	54 // return this.mindate;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	55 // }
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	56
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	57 // private Date maxdate;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	58
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	59 // public Date getMaxDate() {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	60 // return this.maxdate;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	61 // }
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	62
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	63 private int cnt;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	64
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	65 public int getCnt() {
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	66 return this.cnt;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	67 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	68
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	69 // private double mindepth;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	70
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	71 // public double getMinDepth() {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	72 // return this.mindepth;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	73 // }
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	74
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	75 // private double maxdepth;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	76
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	77 // public double getMaxDepth() {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	78 // return this.maxdepth;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	79 // }
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	80
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	81 private double d50;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	82
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	83 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	84 * D50 in m
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	85 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	86 public double getD50() {
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	87 return this.d50;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	88 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	89
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	90 // /**
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	91 // * Parameter constructor
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	92 // */
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	93 // public D50Measurement(final double km, final Date mindate, final Date maxdate, final int cnt, final double mindepth,
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	94 // final double maxdepth,
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	95 // final double d50mm) {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	96 // this.km = km;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	97 // this.mindate = mindate;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	98 // this.maxdate = maxdate;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	99 // this.cnt = cnt;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	100 // this.mindepth = mindepth;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	101 // this.maxdepth = maxdepth;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	102 // this.d50 = d50mm / 1000;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	103 // }
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	104
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	105 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	106 * Query result row constructor
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	107 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	108 public D50Measurement(final Object[] tuple, final String[] aliases) {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	109 this.km = 0;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	110 // this.mindate = null;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	111 // this.maxdate = null;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	112 this.cnt = 0;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	113 // this.mindepth = Double.NaN;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	114 // this.maxdepth = Double.NaN;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	115 this.d50 = Double.NaN;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	116 for (int i = 0; i < tuple.length; ++i) {
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	117 if (tuple[i] == null)
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	118 continue;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	119 switch (aliases[i]) {
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	120 case "km":
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	121 this.km = ((Number) tuple[i]).doubleValue();
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	122 break;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	123 // case "mindate":
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	124 // this.mindate = (Date) tuple[i];
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	125 // break;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	126 // case "maxdate":
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	127 // this.maxdate = (Date) tuple[i];
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	128 // break;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	129 case "cnt":
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	130 this.cnt = ((Number) tuple[i]).intValue();
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	131 break;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	132 // case "mindepth":
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	133 // this.mindepth = ((Number) tuple[i]).doubleValue();
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	134 // break;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	135 // case "maxdepth":
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	136 // this.maxdepth = ((Number) tuple[i]).doubleValue();
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	137 // break;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	138 case "d50":
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	139 this.d50 = ((Number) tuple[i]).doubleValue() / 1000; // mm to m
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	140 break;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	141 default:
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	142 break;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	143 }
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	144 }
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	145 }
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	146 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	147
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	148 /*** FIELDS ***/
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	149
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	150 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	151 * Private log to use here.
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	152 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	153 private static Logger log = Logger.getLogger(BedQualityD50KmValueFinder.class);
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	154
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	155 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	156 * Query that aggregates by km for a km range and a time period all sub layer bed measurements with their d50<br />
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	157 * <br />
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	158 * A km may have bed measurements for multiple dates, multiple distances from the river bank, and multiple depth layers.
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	159 * The query filters by km range, time period and layer (sub layer: below bed to max. 50 cm depth).
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	160 * Those measurements are then grouped by km, and the D50 aggregated as average value.
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	161 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	162 private static final String SQL_BED_D50_SUBLAYER_MEASUREMENT =
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	163 "SELECT t.km, MIN(t.datum) AS mindate, MAX(t.datum) AS maxdate, COUNT(*) AS cnt," //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	164 + " MIN(p.tiefevon) AS mindepth, MAX(p.tiefebis) AS maxdepth, AVG(a.d50) AS d50" //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	165 + " FROM sohltest t INNER JOIN station s ON t.stationid = s.stationid" //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	166 + " INNER JOIN gewaesser g ON s.gewaesserid = g.gewaesserid" //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	167 + " INNER JOIN sohlprobe p ON t.sohltestid = p.sohltestid" //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	168 + " INNER JOIN siebanalyse a ON p.sohlprobeid = a.sohlprobeid" //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	169 + " WHERE (g.name = :name) AND (s.km BETWEEN :fromkm - 0.0001 AND :tokm + 0.0001)" //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	170 + " AND (p.tiefevon > 0.0) AND (p.tiefebis <= 0.5)" //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	171 + " AND (t.datum BETWEEN :fromdate AND :todate)" //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	172 + " GROUP BY t.km" //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	173 + " ORDER BY t.km"; //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	174
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	175 private static final String[] SQL_BED_D50_SELECT_ALIAS = { "km", "mindate", "maxdate", "cnt", "mindepth", "maxdepth", "d50" };
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	176
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	177 /**
8911 37ff7f435912 SINFO Flowdepth: more error checks, d50 interpolation, avoid negative tkh mschaefer parents: 8901 diff changeset	178 * Real linear interpolator for kms and d50 values
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	179 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	180 private final PolynomialSplineFunction interpolator;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	181
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	182 /*** METHODS ***/
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	183
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	184 private BedQualityD50KmValueFinder(final double[] kms, final double[] values) {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	185 this.interpolator = new LinearInterpolator().interpolate(kms, values);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	186 }
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	187
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	188 /**
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	189 * Sohlbeschaffenheit (D50 Korndurchmesser aus Seddb)
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	190 * Abhängig von Peiljahr
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	191 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	192 public static BedQualityD50KmValueFinder loadBedMeasurements(final River river, final DoubleRange kmRange, final int soundingYear, final int validYears) {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	193
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	194 /* construct valid measurement time range */
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	195 final Calendar cal = Calendar.getInstance();
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	196 cal.clear();
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	197
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	198 cal.set(soundingYear - validYears, 0, 1);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	199 final Date startTime = cal.getTime();
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	200
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	201 cal.set(soundingYear + validYears, 11, 31);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	202 final Date endTime = cal.getTime();
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	203
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	204 log.debug(String.format("loadValues km %.3f - %.3f %tF - %tF", kmRange.getMinimumDouble(), kmRange.getMaximumDouble(), startTime, endTime));
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	205 final Session session = SedDBSessionHolder.HOLDER.get();
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	206 final SQLQuery sqlQuery = session.createSQLQuery(SQL_BED_D50_SUBLAYER_MEASUREMENT).addScalar("km", StandardBasicTypes.DOUBLE)
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	207 .addScalar("mindate", StandardBasicTypes.DATE).addScalar("maxdate", StandardBasicTypes.DATE).addScalar("cnt", StandardBasicTypes.INTEGER)
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	208 .addScalar("mindepth", StandardBasicTypes.DOUBLE).addScalar("maxdepth", StandardBasicTypes.DOUBLE).addScalar("d50", StandardBasicTypes.DOUBLE);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	209 final String seddbRiver = river.nameForSeddb();
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	210 sqlQuery.setString("name", seddbRiver);
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	211 sqlQuery.setDouble("fromkm", kmRange.getMinimumDouble());
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	212 sqlQuery.setDouble("tokm", kmRange.getMaximumDouble());
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	213 sqlQuery.setDate("fromdate", startTime);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	214 sqlQuery.setDate("todate", endTime);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	215
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	216 final List<Object[]> rows = sqlQuery.list();
8911 37ff7f435912 SINFO Flowdepth: more error checks, d50 interpolation, avoid negative tkh mschaefer parents: 8901 diff changeset	217 final double[] kms = new double[rows.size()];
37ff7f435912 SINFO Flowdepth: more error checks, d50 interpolation, avoid negative tkh mschaefer parents: 8901 diff changeset	218 final double[] values = new double[rows.size()];
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	219 D50Measurement measurement;
8911 37ff7f435912 SINFO Flowdepth: more error checks, d50 interpolation, avoid negative tkh mschaefer parents: 8901 diff changeset	220 int i = -1;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	221 for (final Object[] row : rows) {
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	222 measurement = new D50Measurement(row, SQL_BED_D50_SELECT_ALIAS);
8911 37ff7f435912 SINFO Flowdepth: more error checks, d50 interpolation, avoid negative tkh mschaefer parents: 8901 diff changeset	223 i++;
37ff7f435912 SINFO Flowdepth: more error checks, d50 interpolation, avoid negative tkh mschaefer parents: 8901 diff changeset	224 kms[i] = measurement.getKm();
37ff7f435912 SINFO Flowdepth: more error checks, d50 interpolation, avoid negative tkh mschaefer parents: 8901 diff changeset	225 values[i] = measurement.getD50();
37ff7f435912 SINFO Flowdepth: more error checks, d50 interpolation, avoid negative tkh mschaefer parents: 8901 diff changeset	226 log.debug(String.format("loadValues km %.3f d50(mm) %.1f count %d", kms[i], values[i], measurement.getCnt()));
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	227 }
8911 37ff7f435912 SINFO Flowdepth: more error checks, d50 interpolation, avoid negative tkh mschaefer parents: 8901 diff changeset	228 try {
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	229 return new BedQualityD50KmValueFinder(kms, values);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	230 }
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	231 catch (final Exception e) {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	232 e.printStackTrace();
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	233 return null;
8911 37ff7f435912 SINFO Flowdepth: more error checks, d50 interpolation, avoid negative tkh mschaefer parents: 8901 diff changeset	234 }
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	235 }
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	236
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	237 /**
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	238 * Returns the d50 value interpolated according to a km
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	239 *
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	240 * @return d50 (mm) of the km, or NaN
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	241 */
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	242 public double findD50(final double km) throws ArgumentOutsideDomainException {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	243 return this.interpolator.value(km);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	244 /*
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	245 * ohne interpolation:
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	246 * if ((kms == null) \|\| (kms.size() == 0))
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	247 * return Double.NaN;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	248 * int i = kms.binarySearch(km);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	249 * if (i >= 0)
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	250 * return values.get(i);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	251 * i = -i - 1;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	252 * if ((i - 1 >= 0) && Utils.epsilonEquals(km, kms.get(i - 1), 0.0001))
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	253 * return values.get(i - 1);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	254 * else if ((i >= 0) && (i <= kms.size() - 1) && Utils.epsilonEquals(km, kms.get(i), 0.0001))
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	255 * return values.get(i);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	256 * else
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	257 * return Double.NaN;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	258 */
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	259 }
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8911 diff changeset	260 }

Mercurial > dive4elements > river

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/tkhcalculation/BedQualityD50KmValueFinder.java @ 8915:d9dbf0b74bc2