Mercurial > dive4elements > river

annotate artifacts/src/main/java/org/dive4elements/river/exports/fixings/FixATWriter.java @ 8202:e4606eae8ea5

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
sed src/**/*.java 's/logger/log/g'
author Sascha L. Teichmann <teichmann@intevation.de>
date Fri, 05 Sep 2014 12:58:17 +0200
parents f4b3861712a4
children 47199406994a
rev   line source
5863
4897a58c8746 River artifacts: Added new copyright headers.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 5838
diff changeset
1 /* Copyright (C) 2011, 2012, 2013 by Bundesanstalt für Gewässerkunde
4897a58c8746 River artifacts: Added new copyright headers.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 5838
diff changeset
2 * Software engineering by Intevation GmbH
4897a58c8746 River artifacts: Added new copyright headers.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 5838
diff changeset
3 *
5994
af13ceeba52a Removed trailing whitespace.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 5863
diff changeset
4 * This file is Free Software under the GNU AGPL (>=v3)
5863
4897a58c8746 River artifacts: Added new copyright headers.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 5838
diff changeset
5 * and comes with ABSOLUTELY NO WARRANTY! Check out the
5994
af13ceeba52a Removed trailing whitespace.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 5863
diff changeset
6 * documentation coming with Dive4Elements River for details.
5863
4897a58c8746 River artifacts: Added new copyright headers.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 5838
diff changeset
7 */
4897a58c8746 River artifacts: Added new copyright headers.
Sascha L. Teichmann <teichmann@intevation.de>
parents: 5838
diff changeset
8
5831
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 3905
diff changeset
9 package org.dive4elements.river.exports.fixings;
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
10
5831
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 3905
diff changeset
11 import org.dive4elements.artifacts.CallMeta;
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
12
5831
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 3905
diff changeset
13 import org.dive4elements.river.artifacts.math.fitting.Function;
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
14
5831
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 3905
diff changeset
15 import org.dive4elements.river.artifacts.model.Parameters;
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
16
5831
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 3905
diff changeset
17 import org.dive4elements.river.artifacts.resources.Resources;
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 3905
diff changeset
18
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 3905
diff changeset
19 import org.dive4elements.river.exports.ATWriter;
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
20
6325
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
21 import org.dive4elements.river.model.Gauge;
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
22 import org.dive4elements.river.model.River;
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
23
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
24 import java.io.IOException;
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
25 import java.io.PrintWriter;
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
26 import java.io.Writer;
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
27
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
28 import java.util.Locale;
7687
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value.
Andre Heinecke <aheinecke@intevation.de>
parents: 6637
diff changeset
29 import java.util.Arrays;
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
30
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
31 import org.apache.log4j.Logger;
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
32
3905
8bbb9e173297 Cosmetics, docs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 3643
diff changeset
33 /** Export Fixation Analysis Results to AT. */
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
34 public class FixATWriter
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
35 {
8202
e4606eae8ea5 sed src/**/*.java 's/logger/log/g'
Sascha L. Teichmann <teichmann@intevation.de>
parents: 7687
diff changeset
36 /** Private log. */
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
37 private static Logger log = Logger.getLogger(FixATWriter.class);
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
38
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
39 public static final String I18N_HEADER_KEY =
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
40 "fix.export.at.header";
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
41
6325
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
42 public static final String I18N_GAUGE_HEADER_KEY =
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
43 "fix.export.at.gauge.header";
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
44
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
45 public static final String I18N_HEADER_DEFAULT =
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
46 "Exported fixings discharge curve for {0} {0}-km: {1}";
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
47
6325
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
48 public static final String I18N_GAUGE_HEADER_DEFAULT =
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
49 "Exported fixings discharge curve for {0}, gauge: {1} datum[{3}] = {2}";
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
50
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
51 public static final String [] Q_MAX_COLUMN = new String [] { "max_q" };
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
52
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
53 private static final int MAX_ITERATIONS = 10000;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
54 private static final double EPSILON = 1e-8;
3643
171db4d5d3cb FixA: AT export: prevent some numerical problems with steep functions around zero.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3642
diff changeset
55 private static final double MIN_Q = 1e-4;
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
56
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
57 protected Function function;
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
58 protected Parameters parameters;
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
59
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
60 public FixATWriter() {
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
61 }
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
62
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
63 public FixATWriter(Function function, Parameters parameters) {
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
64 this.function = function;
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
65 this.parameters = parameters;
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
66 }
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
67
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
68 public void write(
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
69 Writer writer,
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
70 CallMeta meta,
6325
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
71 River river,
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
72 double km
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
73 )
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
74 throws IOException {
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
75 PrintWriter out = new PrintWriter(writer);
6325
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
76
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
77 int subtractPNP = 0;
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
78 // Special case handling for at's at gauges
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
79 Gauge gauge = river.determineGaugeByPosition(km);
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
80 if (Math.abs(km - gauge.getStation().doubleValue()) < 1e-4) {
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
81 printGaugeHeader(out, meta, river, gauge);
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
82 subtractPNP = (int)Math.round(gauge.getDatum().doubleValue() * 100);
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
83 } else {
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
84 printHeader(out, meta, river.getName(), km);
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
85 }
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
86
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
87 double [] coeffs = parameters.interpolate(
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
88 "km", km, function.getParameterNames());
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
89
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
90 double [] qMax = parameters.interpolate(
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
91 "km", km, Q_MAX_COLUMN);
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
92
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
93 if (coeffs == null || qMax == null) {
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
94 log.debug("No data found at km " + km + ".");
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
95 return;
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
96 }
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
97
5831
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 3905
diff changeset
98 org.dive4elements.river.artifacts.math.Function funcInst =
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
99 function.instantiate(coeffs);
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
100
3642
83c0735092a3 FixA: Increase Q max about 5% if storing as AT file.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3641
diff changeset
101 // Increase Q max about 5%.
83c0735092a3 FixA: Increase Q max about 5% if storing as AT file.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3641
diff changeset
102 qMax[0] += Math.abs(qMax[0])*0.05;
83c0735092a3 FixA: Increase Q max about 5% if storing as AT file.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3641
diff changeset
103
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
104 double wMax = funcInst.value(qMax[0]);
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
105
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
106 if (Double.isNaN(wMax) || wMax < 0d) {
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
107 log.debug("function '" + function.getName() +
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
108 "' eval failed at " + wMax);
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
109 return;
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
110 }
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
111
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
112 Function inverse = function.getInverse();
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
113
5831
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 3905
diff changeset
114 org.dive4elements.river.artifacts.math.Function invInst =
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
115 inverse.instantiate(coeffs);
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
116
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
117 double wMin = minW(invInst, wMax, qMax[0]);
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
118
6637
65a5deaee769 Subtract PNP before calculating wcm from wMinCM (issue1155)
Andre Heinecke <aheinecke@intevation.de>
parents: 6325
diff changeset
119 double wMinCM = wMin * 100d - subtractPNP;
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
120 double wMaxCM = wMax * 100d;
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
121
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
122 if ((wMinCM - (int)wMinCM) > 0d) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
123 wMinCM = (int)wMinCM + 1d;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
124 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
125
7687
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value.
Andre Heinecke <aheinecke@intevation.de>
parents: 6637
diff changeset
126 int wRow = ((int)wMinCM / 10) * 10;
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value.
Andre Heinecke <aheinecke@intevation.de>
parents: 6637
diff changeset
127
6637
65a5deaee769 Subtract PNP before calculating wcm from wMinCM (issue1155)
Andre Heinecke <aheinecke@intevation.de>
parents: 6325
diff changeset
128 double w = (wMinCM + subtractPNP) / 100.0;
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
129
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
130 int wcm = ((int)wMinCM) % 10;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
131
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
132 if (log.isDebugEnabled()) {
7687
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value.
Andre Heinecke <aheinecke@intevation.de>
parents: 6637
diff changeset
133 log.debug("km: " + km);
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
134 log.debug("wMinCM: " + wMinCM);
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
135 log.debug("wMaxCM: " + wMaxCM);
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
136 log.debug("wcm: " + wcm);
6637
65a5deaee769 Subtract PNP before calculating wcm from wMinCM (issue1155)
Andre Heinecke <aheinecke@intevation.de>
parents: 6325
diff changeset
137 log.debug("subtractPNP: " + subtractPNP);
7687
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value.
Andre Heinecke <aheinecke@intevation.de>
parents: 6637
diff changeset
138 log.debug("coeffs: " + Arrays.toString(coeffs));
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value.
Andre Heinecke <aheinecke@intevation.de>
parents: 6637
diff changeset
139 log.debug("function description: " + inverse.getDescription());
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
140 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
141
6637
65a5deaee769 Subtract PNP before calculating wcm from wMinCM (issue1155)
Andre Heinecke <aheinecke@intevation.de>
parents: 6325
diff changeset
142 out.printf(Locale.US, "%8d", wRow);
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
143
6209
1642e08d14d2 Fix empty column logic in FixAT export
Andre Heinecke <aheinecke@intevation.de>
parents: 6208
diff changeset
144 for (int i = 0; i < wcm; i++) {
1642e08d14d2 Fix empty column logic in FixAT export
Andre Heinecke <aheinecke@intevation.de>
parents: 6208
diff changeset
145 out.print(ATWriter.EMPTY);
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
146 }
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
147
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
148 for (;;) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
149 while (wcm++ < 10) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
150 if (w > wMax) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
151 break;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
152 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
153 double q = invInst.value(w);
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
154 if (Double.isNaN(w)) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
155 out.print(ATWriter.EMPTY);
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
156 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
157 else {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
158 ATWriter.printQ(out, q);
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
159 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
160 w += 0.01d;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
161 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
162 out.println();
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
163 if (w > wMax) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
164 break;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
165 }
6637
65a5deaee769 Subtract PNP before calculating wcm from wMinCM (issue1155)
Andre Heinecke <aheinecke@intevation.de>
parents: 6325
diff changeset
166 out.printf(Locale.US, "%8d", (wRow += 10));
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
167 wcm = 0;
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
168 }
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
169
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
170 out.flush();
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
171 }
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
172
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
173 protected void printHeader(
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
174 PrintWriter out,
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
175 CallMeta meta,
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
176 String river,
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
177 double km
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
178 ) {
6208
715e010d40de Add * in the first line of AT Exports.
Andre Heinecke <aheinecke@intevation.de>
parents: 5994
diff changeset
179 out.println("*" + Resources.format(
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
180 meta,
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
181 I18N_HEADER_KEY,
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
182 I18N_HEADER_DEFAULT,
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
183 river, km));
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
184 }
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
185
6325
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
186 protected void printGaugeHeader(
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
187 PrintWriter out,
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
188 CallMeta meta,
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
189 River river,
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
190 Gauge gauge
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
191 ) {
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
192 out.println("*" + Resources.format(
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
193 meta,
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
194 I18N_GAUGE_HEADER_KEY,
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
195 I18N_GAUGE_HEADER_DEFAULT,
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
196 new Object[] { river.getName(), gauge.getName(),
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
197 gauge.getDatum(), river.getWstUnit().getName() }));
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
198 }
ee701de9c875 Add special case handling for Fixing AT export at a gauge
Andre Heinecke <aheinecke@intevation.de>
parents: 6209
diff changeset
199
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
200 private static double minW(
5831
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 3905
diff changeset
201 org.dive4elements.river.artifacts.math.Function function,
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
202 double maxW,
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
203 double maxQ
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
204 ) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
205 double stepWidth = 10d;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
206
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
207 double lastW = maxW;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
208 double lastQ = maxQ;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
209
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
210 for (int i = 0; i < MAX_ITERATIONS; ++i) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
211 double w = lastW - stepWidth;
7687
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value.
Andre Heinecke <aheinecke@intevation.de>
parents: 6637
diff changeset
212
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value.
Andre Heinecke <aheinecke@intevation.de>
parents: 6637
diff changeset
213 if (w <= 0) {
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value.
Andre Heinecke <aheinecke@intevation.de>
parents: 6637
diff changeset
214 return 0;
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value.
Andre Heinecke <aheinecke@intevation.de>
parents: 6637
diff changeset
215 }
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value.
Andre Heinecke <aheinecke@intevation.de>
parents: 6637
diff changeset
216
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
217 double q = function.value(w);
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
218
3643
171db4d5d3cb FixA: AT export: prevent some numerical problems with steep functions around zero.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3642
diff changeset
219 if (Double.isNaN(q) || q > lastQ || q < MIN_Q) {
3641
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
220 if (stepWidth < EPSILON) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
221 break;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
222 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
223 stepWidth *= 0.5d;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
224 continue;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
225 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
226
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
227 lastW = w;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
228 lastQ = q;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
229 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
230
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
231 return lastW;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents: 3220
diff changeset
232 }
3220
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
233 }
0c8a6145098b FixA: Added AT writer for fitted functions.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff changeset
234 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :
http://dive4elements.wald.intevation.org