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

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/tkhcalculation/FlowVelocityModelKmValueFinder.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/FlowVelocityModelKmValueFinder.java@89f3c5462a16
children	45f1ad66560e

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: 8898 diff changeset	11 package org.dive4elements.river.artifacts.sinfo.tkhcalculation;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	12
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	13 import java.util.ArrayList;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	14 import java.util.List;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	15
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	16 import org.apache.commons.lang.math.DoubleRange;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	17 import org.apache.log4j.Logger;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	18 import org.dive4elements.river.artifacts.math.Linear;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	19 import org.dive4elements.river.artifacts.math.Utils;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	20 import org.dive4elements.river.artifacts.sinfo.flowdepth.FlowVelocityKmModelValues;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	21 import org.dive4elements.river.backend.SessionHolder;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	22 import org.dive4elements.river.model.River;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	23 import org.hibernate.SQLQuery;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	24 import org.hibernate.Session;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	25 import org.hibernate.type.StandardBasicTypes;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	26
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	27 import gnu.trove.TDoubleArrayList;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	28
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	29 /**
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	30 * Searchable sorted km array with parallel FlowVelocityKmModelValues array and linear interpolation for km and the
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	31 * model values between the array elements.<br />
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	32 * {@link loadValues} loads all the model values for a given km range of a river.<br />
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	33 * {@link findKmQValues} then searches a km in the values table or the nearest including km interval, resp.
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	34 * The v and tau values for a given discharge are either found directly or also interpolated linearly.<br />
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	35 *
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	36 * (Created based on a copy of FlowVelocityMeasurementFactory.)
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	37 *
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	38 * @author Matthias Schäfer
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	39 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	40 final class FlowVelocityModelKmValueFinder {
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	41 /*** FIELDS ***/
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	42
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	43 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	44 * Private log to use here.
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	45 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	46 private static Logger log = Logger.getLogger(FlowVelocityModelKmValueFinder.class);
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	47
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	48 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	49 * Query for a range of stations of a river with all their q, main-v and tau values.<br />
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	50 * (Might be several 10000 rows if many stations and large q range)
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	51 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	52 private static final String SQL_SELECT_ALL = //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	53 "SELECT fvmv.station AS station, fvmv.q AS q, fvmv.main_channel AS vmain, fvmv.shear_stress AS tau"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	54 + " FROM (discharge_zone dz INNER JOIN flow_velocity_model fvm ON dz.id = fvm.discharge_zone_id)"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	55 + " INNER JOIN flow_velocity_model_values fvmv ON fvm.id = fvmv.flow_velocity_model_id"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	56 + " WHERE (dz.river_id = :river_id) AND (fvmv.station BETWEEN :kmfrom - 0.0001 AND :kmto + 0.0001)"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	57 /* + " WHERE (dz.river_id = :river_id) AND (fvmv.q BETWEEN :qmin AND :qmax)" */
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	58 + " ORDER BY fvmv.station ASC, fvmv.q ASC";
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	59
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	60 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	61 * Query for a river's max km below a limit with all its q, main-v and tau values.
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	62 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	63 private static final String SQL_SELECT_KMLOWER = //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	64 "SELECT fvmv.station AS station, fvmv.q AS q, fvmv.main_channel AS vmain, fvmv.shear_stress AS tau" + " FROM flow_velocity_model_values fvmv"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	65 + " INNER JOIN (SELECT MAX(fvmvi.station) AS kmmax"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	66 + " FROM(discharge_zone dz INNER JOIN flow_velocity_model fvm ON dz.id = fvm.discharge_zone_id)"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	67 + " INNER JOIN flow_velocity_model_values fvmvi ON fvm.id = fvmvi.flow_velocity_model_id"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	68 + " WHERE (dz.river_id = :river_id) AND (fvmvi.station < :kmfrom - 0.0001)) finf ON fvmv.station = finf.kmmax"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	69 + " ORDER BY fvmv.q ASC";
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	70
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	71 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	72 * Query for a river's min km above a limit with all its q, main-v and tau values.
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	73 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	74 private static final String SQL_SELECT_KMUPPER = //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	75 "SELECT fvmv.station AS station, fvmv.q AS q, fvmv.main_channel AS vmain, fvmv.shear_stress AS tau" + " FROM flow_velocity_model_values fvmv"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	76 + " INNER JOIN (SELECT MIN(fvmvi.station) AS kmmin"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	77 + " FROM(discharge_zone dz INNER JOIN flow_velocity_model fvm ON dz.id = fvm.discharge_zone_id)"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	78 + " INNER JOIN flow_velocity_model_values fvmvi ON fvm.id = fvmvi.flow_velocity_model_id"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	79 + " WHERE (dz.river_id = :river_id) AND (fvmvi.station > :kmto + 0.0001)) fsup ON fvmv.station = fsup.kmmin"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	80 + " ORDER BY fvmv.q ASC";
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	81
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	82 // /**
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	83 // * Query to select all km-q-v-tau of a river that are the q maxima below a q limit
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	84 // */
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	85 // private static final String SQL_SELECT_QLOWER =
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	86 // "SELECT fvmv.station AS station, fvmv.q AS q, fvmv.main_channel AS vmain, fvmv.shear_stress AS tau"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	87 // + " FROM flow_velocity_model_values fvmv"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	88 // + " INNER JOIN (SELECT fv2.station, MAX(fv2.q) AS q"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	89 // + " FROM (discharge_zone dz INNER JOIN flow_velocity_model fvm ON dz.id = fvm.discharge_zone_id)"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	90 // + " INNER JOIN flow_velocity_model_values fv2 ON fvm.id = fv2.flow_velocity_model_id"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	91 // + " WHERE (dz.river_id = :river_id) AND (fv2.q < :qlim) GROUP BY fv2.station) qx"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	92 // + " ON (fvmv.station=qx.station) AND (fvmv.q=qx.q)"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	93 // + " ORDER BY fvmv.station ASC";
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	94 //
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	95 // /**
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	96 // * Query to select all km-q-v-tau of a river that are the q minima above a q limit
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	97 // */
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	98 // private static final String SQL_SELECT_QUPPER =
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	99 // "SELECT fvmv.station AS station, fvmv.q AS q, fvmv.main_channel AS vmain, fvmv.shear_stress AS tau"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	100 // + " FROM flow_velocity_model_values fvmv"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	101 // + " INNER JOIN (SELECT fv2.station, MIN(fv2.q) AS q"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	102 // + " FROM (discharge_zone dz INNER JOIN flow_velocity_model fvm ON dz.id = fvm.discharge_zone_id)"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	103 // + " INNER JOIN flow_velocity_model_values fv2 ON fvm.id = fv2.flow_velocity_model_id"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	104 // + " WHERE (dz.river_id = :river_id) AND (fv2.q > :qlim) GROUP BY fv2.station) qx"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	105 // + " ON (fvmv.station=qx.station) AND (fvmv.q=qx.q)"
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	106 // + " ORDER BY fvmv.station ASC";
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	107
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	108 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	109 * Kms of the loaded river range
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	110 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	111 private final TDoubleArrayList kms = new TDoubleArrayList();
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	112
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	113 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	114 * For each km in kms a list of q-v-tau-tupels
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	115 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	116 private final List<FlowVelocityKmModelValues> values = new ArrayList<>();
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	117
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	118 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	119 * Searched km of the last findKmValue
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	120 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	121 private double findKm;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	122
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	123 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	124 * kms and values index of the interval start found by the last findKmValue
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	125 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	126 private int leftIndexFound = -1;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	127
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	128 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	129 * kms and values index of the interval end found by the last findKmValue
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	130 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	131 private int rightIndexFound = -1;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	132
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	133 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	134 * Q of the last findKmQValues
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	135 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	136 private double findQ;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	137
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	138 /*** METHODS ***/
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	139
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	140 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	141 * Discharge of the last {@link findKmQValue}
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	142 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	143 public double getFindQ() {
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	144 return this.findQ;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	145 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	146
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	147 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	148 * Velocity of the last {@link findKmQValues}
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	149 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	150 public double getFindVmainFound() {
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	151 if (this.leftIndexFound < 0)
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	152 return Double.NaN;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	153 else if (this.leftIndexFound == this.rightIndexFound)
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	154 return getLeftValues().getVmainFound();
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	155 else
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	156 return Linear.linear(this.findKm, getLeftValues().getKm(), getRightValues().getKm(), getLeftValues().getVmainFound(),
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	157 getRightValues().getVmainFound());
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	158 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	159
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	160 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	161 * Shear stress tau of the last {@link findKmQValues}
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	162 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	163 public double getFindTauFound() {
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	164 if (this.leftIndexFound < 0)
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	165 return Double.NaN;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	166 else if (this.leftIndexFound == this.rightIndexFound)
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	167 return getLeftValues().getTauFound();
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	168 else
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	169 return Linear.linear(this.findKm, getLeftValues().getKm(), getRightValues().getKm(), getLeftValues().getTauFound(), getRightValues().getTauFound());
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	170 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	171
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	172 /**
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	173 * Whether the discharge has been interpolated in the last {@link findKmQValues}
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	174 */
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	175 public boolean getFindIsQInterpolated() {
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	176 return (getLeftValues() != null) && (getLeftValues().getIsInterpolated() \|\| getRightValues().getIsInterpolated());
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	177 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	178
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: 8898 diff changeset	180 * Static constructor: queries a range of a river's kms with all their q-v-tau values.
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	181 *
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	182 * @return Whether the load has been successful the new instance, <code>null</code> otherwise.
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	183 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	184 public static FlowVelocityModelKmValueFinder loadValues(final River river, final DoubleRange kmRange, final DoubleRange qRange) {
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	185 // DB session
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	186 log.debug(String.format("loadValues km %.3f - %.3f / q %.1f - %.1f", kmRange.getMinimumDouble(), kmRange.getMaximumDouble(), qRange.getMinimumDouble(),
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	187 qRange.getMaximumDouble()));
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	188
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	189 final FlowVelocityModelKmValueFinder instance = new FlowVelocityModelKmValueFinder();
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	190
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	191 final TDoubleArrayList kms = instance.kms;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	192 final List<FlowVelocityKmModelValues> values = instance.values;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	193
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	194 final boolean isDemoValuesCorrection = river.getName().equalsIgnoreCase("beispielfluss");
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	195 final Session session = SessionHolder.HOLDER.get();
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	196
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	197 // Select km infimum
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	198 SQLQuery sqlQuery = session.createSQLQuery(SQL_SELECT_KMLOWER).addScalar("station", StandardBasicTypes.DOUBLE).addScalar("q", StandardBasicTypes.DOUBLE)
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	199 .addScalar("vmain", StandardBasicTypes.DOUBLE).addScalar("tau", StandardBasicTypes.DOUBLE);
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	200 sqlQuery.setParameter("river_id", river.getId());
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	201 sqlQuery.setParameter("kmfrom", kmRange.getMinimumDouble());
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	202 instance.addKms(sqlQuery.list(), isDemoValuesCorrection);
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	203
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	204 // Select km range
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	205 sqlQuery = session.createSQLQuery(SQL_SELECT_ALL).addScalar("station", StandardBasicTypes.DOUBLE).addScalar("q", StandardBasicTypes.DOUBLE)
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	206 .addScalar("vmain", StandardBasicTypes.DOUBLE).addScalar("tau", StandardBasicTypes.DOUBLE);
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	207 sqlQuery.setParameter("river_id", river.getId());
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	208 sqlQuery.setParameter("kmfrom", kmRange.getMinimumDouble());
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	209 sqlQuery.setParameter("kmto", kmRange.getMaximumDouble());
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	210 // sqlQuery.setParameter("qmin", qRange.getMinimumDouble());
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	211 // sqlQuery.setParameter("qmax", qRange.getMaximumDouble());
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	212
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	213 int kmcount = kms.size();
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	214 final int rowcount = instance.addKms(sqlQuery.list(), isDemoValuesCorrection);
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	215 kmcount = kms.size() - kmcount;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	216
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	217 // Select km supremum
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	218 sqlQuery = session.createSQLQuery(SQL_SELECT_KMUPPER).addScalar("station", StandardBasicTypes.DOUBLE).addScalar("q", StandardBasicTypes.DOUBLE)
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	219 .addScalar("vmain", StandardBasicTypes.DOUBLE).addScalar("tau", StandardBasicTypes.DOUBLE);
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	220 sqlQuery.setParameter("river_id", river.getId());
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	221 sqlQuery.setParameter("kmto", kmRange.getMaximumDouble());
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	222 final int supcnt = instance.addKms(sqlQuery.list(), isDemoValuesCorrection);
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	223
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	224 // Add copy of last km for search of max km value
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	225 if ((supcnt == 0) && (values.size() >= 1)) {
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	226 kms.add(kms.getQuick(kms.size()) + 0.0001);
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	227 values.add(new FlowVelocityKmModelValues(kms.getQuick(kms.size() - 1), values.get(values.size() - 1)));
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	228 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	229
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	230 // log.debug
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	231 if (values.size() - 1 >= 0) {
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	232 log.debug(String.format("loadValues %d: km %.3f - %d values", 0, values.get(0).getKm(), values.get(0).size()));
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	233
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	234 if (values.size() - 1 >= 1) {
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	235 log.debug(String.format("loadValues %d: km %.3f - %d values", 1, values.get(1).getKm(), values.get(1).size()));
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	236
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	237 if (values.size() - 1 >= 2)
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	238 log.debug("loadValues ...");
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	239
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	240 if (values.size() - 2 >= 3)
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	241 log.debug(String.format("loadValues %d: km %.3f - %d values", values.size() - 2, values.get(values.size() - 2).getKm(),
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	242 values.get(values.size() - 2).size()));
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	243
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	244 if (values.size() - 1 >= 3)
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	245 log.debug(String.format("loadValues %d: km %.3f - %d values", values.size() - 1, values.get(values.size() - 1).getKm(),
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	246 values.get(values.size() - 1).size()));
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	247 }
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	248 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	249
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	250 log.debug(String.format("loadValues %d kms, %d values loaded", kmcount, rowcount));
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	251
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	252 if (kms.size() == 0)
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	253 return null;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	254
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	255 return instance;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	256 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	257
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	258 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	259 * Adds the km-q-v-tau values of a query result row to the last km of the list, or a new one resp.
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	260 *
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	261 * @return Number of rows
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	262 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	263 private int addKms(final List<Object[]> rows, final boolean isDemoValuesCorrection) {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	264 for (final Object[] row : rows) {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	265 if ((this.kms.size() == 0) \|\| !Utils.epsilonEquals(this.kms.get(this.kms.size() - 1), (double) row[0], 0.0001)) {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	266 this.kms.add((double) row[0]);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	267 this.values.add(new FlowVelocityKmModelValues(this.kms.get(this.kms.size() - 1)));
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	268 }
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	269 if (isDemoValuesCorrection)
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	270 // "Verfremdung" der v-Werte etwas korrigieren (Originalwerte wurden mit Zufallswert zwischen 10 und 20 multipliziert)
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	271 this.values.get(this.values.size() - 1).addValues((double) row[1], ((double) row[2]) / 10, (double) row[3]);
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	272 else
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	273 this.values.get(this.values.size() - 1).addValues((double) row[1], (double) row[2], (double) row[3]);
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	274 }
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	275 return rows.size();
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	276 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	277
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	278 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	279 * Searches a km and finds or interpolates the velocity and shear stress values for a discharge<br />
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	280 * The values may be got via {@link getVmainFound} etc.
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	281 *
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	282 * @return Whether values have been found
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	283 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	284 public boolean findKmQValues(final double km, final double q) {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	285 this.findQ = q;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	286 if (!searchKm(km))
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	287 return false;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	288 if (this.leftIndexFound == this.rightIndexFound) {
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	289 // Exact km match
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	290 final double qfound = getLeftValues().findQ(q);
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	291 log.debug(String.format("findKmQValues km %.3f q %.0f = %.0f (%d)", km, q, qfound, this.leftIndexFound));
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	292 return !Double.isNaN(qfound);
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	293 } else {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	294 final double[] qfound = { getLeftValues().findQ(q), getRightValues().findQ(q) };
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	295 log.debug(String.format("findKmQValues km %.3f q %.0f = %.0f (%d, %.3f) - %.0f (%d, %.3f)", km, q, qfound[0], this.leftIndexFound,
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	296 getLeftValues().getKm(), qfound[1], this.rightIndexFound, getRightValues().getKm()));
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	297 return !Double.isNaN(qfound[0]) && !Double.isNaN(qfound[1]);
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	298 }
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	299 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	300
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	301 /**
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	302 * Searches a km
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	303 *
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	304 * @return Whether the km was within the supported range
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	305 */
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	306 private boolean searchKm(final double km) {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	307 this.findKm = km;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	308 this.leftIndexFound = -1;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	309 this.rightIndexFound = -1;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	310 if ((this.kms == null) \|\| (this.kms.size() == 0))
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	311 return false;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	312 int i = this.kms.binarySearch(km);
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	313 if (i >= 0) {
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	314 // Exact km match
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	315 this.leftIndexFound = i;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	316 this.rightIndexFound = i;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	317 return true;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	318 } else {
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	319 // Out of range or within km interval
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	320 if (i < 0)
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	321 i = -i - 1;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	322 if ((i <= 0) \|\| (i >= this.kms.size()))
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	323 return false;
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	324 this.leftIndexFound = i - 1;
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	325 this.rightIndexFound = i;
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	326 return true;
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	327 }
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	328 }
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	329
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	330 private FlowVelocityKmModelValues getLeftValues() {
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	331 return this.values.get(this.leftIndexFound);
8898 89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	332 }
89f3c5462a16 Implemented S-INFO Flowdepth TKH calculation mschaefer parents: diff changeset	333
8915 d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	334 private FlowVelocityKmModelValues getRightValues() {
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	335 return this.values.get(this.rightIndexFound);
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	336 }
d9dbf0b74bc2 Refaktoring of flow depth calculation, extracting tkh part. First implementation of tkh calculation. gernotbelger parents: 8898 diff changeset	337 }

Mercurial > dive4elements > river

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