Mercurial > dive4elements > river

comparison artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/flowdepth/SoilKindKmValueFinder.java @ 8898:89f3c5462a16

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Implemented S-INFO Flowdepth TKH calculation
author mschaefer
date Thu, 22 Feb 2018 12:03:31 +0100
parents
children
comparison
equal deleted inserted replaced
8891:f431aec10d2c 8898:89f3c5462a16
1 /* Copyright (C) 2017 by Bundesanstalt für Gewässerkunde
2 * Software engineering by
3 * Björnsen Beratende Ingenieure GmbH
4 * Dr. Schumacher Ingenieurbüro für Wasser und Umwelt
5 *
6 * This file is Free Software under the GNU AGPL (>=v3)
7 * and comes with ABSOLUTELY NO WARRANTY! Check out the
8 * documentation coming with Dive4Elements River for details.
9 */
10
11 package org.dive4elements.river.artifacts.sinfo.flowdepth;
12
13 import java.util.ArrayList;
14 import java.util.List;
15
16 import org.apache.commons.lang.math.DoubleRange;
17 import org.apache.commons.math.ArgumentOutsideDomainException;
18 import org.apache.log4j.Logger;
19 import org.dive4elements.river.model.River;
20
21 import gnu.trove.TDoubleArrayList;
22
23 /**
24 * @author matthias
25 *
26 */
27 public class SoilKindKmValueFinder
28 {
29 /**
30 * Private log to use here.
31 */
32 private static Logger log = Logger.getLogger(SoilKindKmValueFinder.class);
33
34 private TDoubleArrayList kms;
35
36 private List<SoilKind> values;
37
38 /***** METHODS *****/
39
40 /**
41 * Searches a km with its soil kind
42 */
43 public SoilKind findSoilKind(double km) throws ArgumentOutsideDomainException {
44 if ((kms == null) || (kms.size() == 0))
45 throw new ArgumentOutsideDomainException(km, Double.NaN, Double.NaN);
46 int i = kms.binarySearch(km);
47 if (i >= 0) {
48 // Exact km match
49 return values.get(i);
50 }
51 else {
52 // Out of range or within km interval
53 if (i < 0)
54 i = -i - 1;
55 if ((i <= 0) || (i >= kms.size()))
56 throw new ArgumentOutsideDomainException(km, Double.NaN, Double.NaN);
57 if (km <= ((kms.get(i-1) + kms.get(i)) / 2))
58 return values.get(i-1);
59 else
60 return values.get(i);
61 }
62 }
63
64 /**
65 * Loads the range of the river's kms with their soil kind.
66 * @return Whether the load has been successful
67 */
68 public boolean loadValues(River river, DoubleRange kmRange) {
69 kms = new TDoubleArrayList();
70 values = new ArrayList<SoilKind>();
71 //FIXME Echte Daten aus der Datenbank abfragen
72 addKmKind(0, SoilKind.starr);
73 addKmKind(15.7, SoilKind.mobil);
74 addKmKind(15.8, SoilKind.mobil);
75 addKmKind(15.9, SoilKind.starr);
76 addKmKind(108.7, SoilKind.mobil);
77 addKmKind(108.8, SoilKind.mobil);
78 addKmKind(108.9, SoilKind.starr);
79 addKmKind(119.1, SoilKind.mobil);
80 addKmKind(119.4, SoilKind.mobil);
81 addKmKind(119.5, SoilKind.starr);
82 addKmKind(128.3, SoilKind.mobil);
83 addKmKind(128.9, SoilKind.mobil);
84 addKmKind(129, SoilKind.starr);
85 addKmKind(133.1, SoilKind.mobil);
86 addKmKind(135.9, SoilKind.mobil);
87 addKmKind(136, SoilKind.starr);
88 addKmKind(136.5, SoilKind.mobil);
89 addKmKind(139.9, SoilKind.mobil);
90 addKmKind(140, SoilKind.starr);
91 addKmKind(140.5, SoilKind.mobil);
92 addKmKind(165, SoilKind.mobil);
93 addKmKind(165.1, SoilKind.starr);
94 addKmKind(165.9, SoilKind.mobil);
95 addKmKind(180.8, SoilKind.mobil);
96 addKmKind(180.9, SoilKind.starr);
97 addKmKind(182, SoilKind.mobil);
98 addKmKind(221.3, SoilKind.mobil);
99 return true;
100 }
101
102 private void addKmKind(double km, SoilKind kind) {
103 kms.add(km);
104 values.add(kind);
105 }
106 }
http://dive4elements.wald.intevation.org