Mercurial > dive4elements > river
annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/minfo/SedimentLoadFacet.java @ 5653:5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Allow gaps in graph (not sorting, living with NaNs).
author | Felix Wolfsteller <felix.wolfsteller@intevation.de> |
---|---|
date | Thu, 11 Apr 2013 10:52:56 +0200 |
parents | 4feda81c38bc |
children | 5bb179d4fd5f |
rev | line source |
---|---|
4372
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
1 package de.intevation.flys.artifacts.model.minfo; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
2 |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
3 import org.apache.log4j.Logger; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
4 |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
5 import java.util.ArrayList; |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
6 import java.util.Arrays; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
7 import java.util.Collections; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
8 import java.util.Comparator; |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
9 import java.util.List; |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
10 import java.util.Map; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
11 import java.util.TreeMap; |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
12 |
4372
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
13 import de.intevation.artifactdatabase.state.Facet; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
14 import de.intevation.artifacts.Artifact; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
15 import de.intevation.artifacts.CallContext; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
16 import de.intevation.flys.artifacts.FLYSArtifact; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
17 import de.intevation.flys.artifacts.model.CalculationResult; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
18 import de.intevation.flys.artifacts.model.DataFacet; |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
19 import de.intevation.flys.model.MeasurementStation; |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
20 import de.intevation.flys.artifacts.model.FacetTypes; |
4372
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
21 import de.intevation.flys.artifacts.states.DefaultState.ComputeType; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
22 |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
23 import de.intevation.flys.utils.FLYSUtils; |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
24 |
4372
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
25 |
5647
ddb2a4e982b8
Doc and cosmetics.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5645
diff
changeset
|
26 /** Facet to access various sediment loads. */ |
4372
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
27 public class SedimentLoadFacet |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
28 extends DataFacet |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
29 { |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
30 /** Very own logger. */ |
4372
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
31 private static Logger logger = Logger.getLogger(SedimentLoadFacet.class); |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
32 |
5647
ddb2a4e982b8
Doc and cosmetics.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5645
diff
changeset
|
33 /** Used as tolerance value when fetching measurement stations. */ |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
34 private static double EPSILON = 1e-5; |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
35 |
5647
ddb2a4e982b8
Doc and cosmetics.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5645
diff
changeset
|
36 |
4372
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
37 public SedimentLoadFacet() { |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
38 } |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
39 |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
40 public SedimentLoadFacet(int idx, String name, String description, |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
41 ComputeType type, String stateId, String hash) { |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
42 super(idx, name, description, type, hash, stateId); |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
43 } |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
44 |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
45 public Object getData(Artifact artifact, CallContext context) { |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
46 logger.debug("Get data for sediment load at index: " + index); |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
47 |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
48 FLYSArtifact flys = (FLYSArtifact) artifact; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
49 |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
50 CalculationResult res = (CalculationResult) flys.compute(context, hash, |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
51 stateId, type, false); |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
52 |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
53 Object[] data = |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
54 (SedimentLoadResult[]) res.getData(); // TODO CAST TO SPECIFIC CLASS |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
55 |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
56 List<MeasurementStation> allStations = FLYSUtils.getRiver(flys).getMeasurementStations(); |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
57 SedimentLoadResult result = data != null && data.length > index ? (SedimentLoadResult)data[index] : null; |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
58 if (result == null) { |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
59 return null; |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
60 } |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
61 |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
62 List<Double> sortedStarts = new ArrayList<Double>(); |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
63 // Filter stations according to type. |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
64 List<MeasurementStation> stations = new ArrayList<MeasurementStation>(); |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
65 for (MeasurementStation station: allStations) { |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
66 if (station.getRange() == null || station.getMeasurementType() == null) { |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
67 continue; |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
68 } |
5647
ddb2a4e982b8
Doc and cosmetics.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5645
diff
changeset
|
69 if (FacetTypes.IS.SEDIMENT_LOAD_NO_FLOAT(this.getName()) |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
70 && station.getMeasurementType().equals("Geschiebe")) { |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
71 stations.add(station); |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
72 sortedStarts.add(station.getStation()); |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
73 } |
5647
ddb2a4e982b8
Doc and cosmetics.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5645
diff
changeset
|
74 else if (!FacetTypes.IS.SEDIMENT_LOAD_NO_FLOAT(this.getName()) |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
75 && station.getMeasurementType().equals("Schwebstoff")) { |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
76 stations.add(station); |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
77 sortedStarts.add(station.getStation()); |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
78 } |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
79 } |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
80 Collections.sort(sortedStarts); |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
81 |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
82 // Access data according to type. |
5648
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
83 double[][] sd = getLoadData(result); |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
84 |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
85 // Sort by km. |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
86 TreeMap<Double, Double> sortData = new TreeMap<Double,Double>(); |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
87 |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
88 double[] km = sd[0]; |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
89 double[] load = sd[1]; |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
90 |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
91 for (int i = 0 ; i < km.length; i++) { |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
92 sortData.put(km[i], load[i]); |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
93 } |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
94 |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
95 double[][] values = new double[2][]; |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
96 values[0] = new double[km.length*3]; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
97 values[1] = new double[km.length*3]; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
98 |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
99 List<double[]> kmWithoutStation = new ArrayList<double[]>(); |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
100 |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
101 // Find station via its station (km). |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
102 // TODO use a binarySearch instead of linear absdiff approach |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
103 int i = 0; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
104 for (Map.Entry<Double, Double> entry: sortData.entrySet()) { |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
105 boolean matchFound = false; |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
106 // For now, ignore overlaps like (B> next A) |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
107 for (MeasurementStation station: stations) { |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
108 if (Math.abs(station.getStation() - entry.getKey()) < EPSILON) { |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
109 values[0][i*3] = station.getRange().getA().doubleValue() + EPSILON; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
110 values[1][i*3] = entry.getValue(); |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
111 values[0][i*3+1] = station.getRange().getB().doubleValue() - EPSILON; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
112 values[1][i*3+1] = entry.getValue(); |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
113 values[0][i*3+2] = station.getRange().getB().doubleValue(); |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
114 values[1][i*3+2] = entry.getValue(); |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
115 matchFound = true; |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
116 } |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
117 } |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
118 // Store points without match for later assessment. |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
119 if (!matchFound) { |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
120 kmWithoutStation.add(new double[] {entry.getKey(), entry.getValue(), i}); |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
121 } |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
122 i++; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
123 } |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
124 |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
125 // Find fitting measurement stations for values without match. |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
126 for (double misses[]: kmWithoutStation) { |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
127 int idc = Math.abs(Collections.binarySearch(sortedStarts, misses[0])); |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
128 double locationA = sortedStarts.get(idc-2); |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
129 double locationB = sortedStarts.get(idc-1); |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
130 values[0][(int)misses[2]*3] = locationA + EPSILON; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
131 values[1][(int)misses[2]*3] = misses[1]; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
132 values[0][(int)misses[2]*3+1] = locationB - EPSILON; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
133 values[1][(int)misses[2]*3+1] = misses[1]; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
134 values[0][(int)misses[2]*3+2] = locationB ; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
135 values[1][(int)misses[2]*3+2] = misses[1]; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
136 } |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
137 |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
138 |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
139 for (int x = 0; x < values[0].length-1; x++) { |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
140 // Correct measurement stationo overlap. |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
141 if (values[0][x] > values[0][x+1]) { |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
142 values[0][x+1] = values[0][x] + EPSILON; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
143 } |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
144 // Introduce gaps where no data in measurement station. |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
145 if (Math.abs(values[0][x+1] - values[0][x]) > 3*EPSILON |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
146 && values[1][x+1] != values[1][x]) { |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
147 values[0][x] = Double.NaN; |
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
148 values[1][x] = Double.NaN; |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
149 } |
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
150 } |
5653
5231e6b849ce
issue1077: Handle overlapping measurement station ranges, refactored Generator,
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5648
diff
changeset
|
151 |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
152 return values; |
4372
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
153 } |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
154 |
5645
696d710470f5
flys/issue1077: Show loads as step line, therefore transform data in
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4372
diff
changeset
|
155 |
5648
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
156 /** Get data according to type of facet. */ |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
157 private double[][] getLoadData(SedimentLoadResult result) { |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
158 if (getName().equals(FacetTypes.SEDIMENT_LOAD_SAND)) |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
159 return result.getSandData(); |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
160 else if (getName().equals(FacetTypes.SEDIMENT_LOAD_COARSE)) |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
161 return result.getCoarseData(); |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
162 else if (getName().equals(FacetTypes.SEDIMENT_LOAD_FINEMIDDLE)) |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
163 return result.getFineMiddleData(); |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
164 else if (getName().equals(FacetTypes.SEDIMENT_LOAD_SUSP_SAND)) |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
165 return result.getSuspSandData(); |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
166 else if (getName().equals(FacetTypes.SEDIMENT_LOAD_SUSP_SAND_BED)) |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
167 return result.getSuspSandBedData(); |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
168 else if (getName().equals(FacetTypes.SEDIMENT_LOAD_SUSP_SEDIMENT)) |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
169 return result.getSuspSedimentData(); |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
170 else if (getName().equals(FacetTypes.SEDIMENT_LOAD_TOTAL_LOAD)) |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
171 return result.getTotalLoadData(); |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
172 else if (getName().equals(FacetTypes.SEDIMENT_LOAD_TOTAL)) |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
173 return result.getTotalData(); |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
174 else { |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
175 logger.error("SedimentLoadFacet " + getName() + " cannot determine data type."); |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
176 return null; |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
177 } |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
178 } |
4feda81c38bc
SedimentLoadFacet: Minor refactor.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5647
diff
changeset
|
179 |
4372
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
180 /** Copy deeply. */ |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
181 @Override |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
182 public Facet deepCopy() { |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
183 SedimentLoadFacet copy = new SedimentLoadFacet(); |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
184 copy.set(this); |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
185 copy.type = type; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
186 copy.hash = hash; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
187 copy.stateId = stateId; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
188 return copy; |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
189 } |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
190 } |
19772b414d46
New facet and result set for sediment load.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
191 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 : |