dive4elements/river: artifacts/src/main/java/org/dive4elements/river/exports/ATWriter.java comparison

comparison artifacts/src/main/java/org/dive4elements/river/exports/ATWriter.java @ 6304:152e1ba8e5bd

Rework ATWriter/Exporter - Make the write function static as both data and gauge information is necessary for writing. - Use the rivers WST unit for the gauge datum - Handle the case that we are at a gauge - Differentiate between calculated values - Always work with w in CM and thus avoid double number format problems

author	Andre Heinecke <aheinecke@intevation.de>
date	Wed, 12 Jun 2013 16:07:28 +0200
parents	02ea1c3fba8c
children	3ac657d60e60

comparison

equal deleted inserted replaced

-:5fda7065b74a
+:152e1ba8e5bd
 "export.discharge.curve.at.header";
 public static final String I18N_AT_GAUGE_HEADER =
 "export.discharge.curve.at.gauge.header";
+public static final String I18N_AT_CALC_GAUGE_HEADER =
+"export.discharge.curve.at.gauge.calc.header";
 public static final String EMPTY = "         ";
-protected double minW;
+public static double getQ(int w, UnivariateRealFunction qFunc) {
-protected double maxW;
-protected double minQ;
-protected double maxQ;
-protected UnivariateRealFunction qFunc;
-public ATWriter() {
-}
-public ATWriter(WQ wq) throws IllegalArgumentException {
-int [] bounds = wq.longestIncreasingWRangeIndices();
-if (logger.isDebugEnabled()) {
-logger.debug("exporting w between indices " +
-bounds[0] + " and " + bounds[1] + " (" +
-wq.getW(bounds[0]) + ", " + wq.getW(bounds[1]));
-}
-if (bounds[1]-bounds[0] < 1) { // Only first w can be written out.
-minW = maxW = wq.getW(bounds[0]);
-minQ = maxQ = wq.getQ(bounds[0]);
-// constant function
-qFunc = new PolynomialFunction(new double [] { minQ });
-return;
-}
-double [] ws = new double[bounds[1]-bounds[0]];
-double [] qs = new double[ws.length];
-for (int i = 0; i < ws.length; ++i) {
-int idx = bounds[0]+i;
-ws[i] = wq.getW(idx);
-qs[i] = wq.getQ(idx);
-}
-qFunc = ws.length < 3
-? new LinearInterpolator().interpolate(ws, qs)
-: new SplineInterpolator().interpolate(ws, qs);
-minW = wq.getW(bounds[0]);
-maxW = wq.getW(bounds[1]);
-minQ = wq.getQ(bounds[0]);
-maxQ = wq.getQ(bounds[1]);
-}
-public double getQ(double w) {
 try {
-return qFunc.value(w);
+double val = qFunc.value(w);
+return val;
 }
 catch (FunctionEvaluationException aode) {
 // should not happen
-logger.warn("spline interpolation failed", aode);
+logger.error("spline interpolation failed", aode);
-return w <= minW ? minQ : maxQ;
+return Double.NaN;
 }
 }
 public static void printQ(PrintWriter out, double q) {
 String format;
 if (q > 1000d) q = Math.rint(q/10d)*10d;
 }
 out.printf(Locale.US, format, q);
 }
+protected static void printCalculatedGaugeHeader(
+PrintWriter out,
+CallMeta    callMeta,
+String      river,
+double      km,
+String      gName,
+BigDecimal  datum,
+Date        date,
+String      unit
+) {
+out.print("*" + Resources.getMsg(
+callMeta,
+I18N_AT_CALC_GAUGE_HEADER,
+I18N_AT_CALC_GAUGE_HEADER,
+new Object[] { river, gName, datum, unit } ));
+out.print("\r\n");
+}
 protected static void printGaugeHeader(
 PrintWriter out,
 CallMeta    callMeta,
 String      river,
 double      km,
 String      gName,
 BigDecimal  datum,
-Date        date
+Date        date,
+String      unit
 ) {
-DateFormat f = DateFormat.getDateInstance();
+DateFormat f = DateFormat.getDateInstance(DateFormat.MEDIUM,
+Resources.getLocale(callMeta));
 out.print("*" + Resources.getMsg(
 callMeta,
 I18N_AT_GAUGE_HEADER,
 I18N_AT_GAUGE_HEADER,
-new Object[] { river, gName, f.format(date), datum } ));
+new Object[] { river, gName, f.format(date), datum, unit} ));
 out.print("\r\n");
 }
 protected static void printHeader(
 PrintWriter out,
 I18N_AT_HEADER,
 new Object[] { river, km } ));
 out.print("\r\n");
 }
-public void write(
+public static void write(
+WQ values,
 Writer writer,
 CallMeta meta,
 String river,
 double km,
 String gName,
 BigDecimal datum,
 Date date,
-boolean wOutUnitIsMeter)
+String unit)
 throws IOException
 {
+int minW;
+int maxW;
+double minQ;
+double maxQ;
+UnivariateRealFunction qFunc;
+WQ wq = WQ.getFixedWQforExportAtGauge(values, datum);
+// If we converted to centimeter we know that the WQ table is
+// calculated because of the assumption that all calculations
+// are in Meter and only the discharge tables data is in meter.
+boolean isCalculation = wq.getReferenceSystem() != values.getReferenceSystem();
+int [] bounds = wq.longestIncreasingWRangeIndices();
+if (logger.isDebugEnabled()) {
+logger.debug("exporting " + (isCalculation ? "calculated " : "") +
+"w between indices " + bounds[0] + " and " + bounds[1] + " (" +
+(int)Math.ceil(wq.getW(bounds[0])) + ", " +
+(int)Math.floor(wq.getW(bounds[1]))+ ")");
+}
+if (bounds[1]-bounds[0] < 1) { // Only first w can be written out.
+minW = maxW = (int)Math.round(wq.getW(bounds[0]));
+minQ = maxQ = wq.getQ(bounds[0]);
+// constant function
+qFunc = new PolynomialFunction(new double [] { minQ });
+return;
+}
+double [] ws = new double[bounds[1]-bounds[0]];
+double [] qs = new double[ws.length];
+for (int i = 0; i < ws.length; ++i) {
+int idx = bounds[0]+i;
+ws[i] = wq.getW(idx);
+qs[i] = wq.getQ(idx);
+}
+qFunc = ws.length < 3
+? new LinearInterpolator().interpolate(ws, qs)
+: new SplineInterpolator().interpolate(ws, qs);
+minW = (int)Math.ceil(wq.getW(bounds[0]));
+maxW = (int)Math.floor(wq.getW(bounds[1]));
+minQ = wq.getQ(bounds[0]);
+maxQ = wq.getQ(bounds[1]);
 PrintWriter out = new PrintWriter(writer);
 // A header is required, because the desktop version of FLYS will skip
 // the first row.
 if (gName != null) {
-printGaugeHeader(out, meta, river, km, gName, datum, date);
+printGaugeHeader(out, meta, river, km, gName, datum, date, unit);
+if (isCalculation) {
+printCalculatedGaugeHeader(out, meta, river, km, gName, datum, date, unit);
+}
 }
 else {
 printHeader(out, meta, river, km);
 }
-double rest = (minW * 100.0) % 10.0;
+int rest = minW % 10;
-double startW = Math.rint((minW - rest*0.01)*10.0)*0.1;
+int startW = minW - rest;
 if (logger.isDebugEnabled()) {
 logger.debug("startW: " + startW);
 logger.debug("rest: " + rest);
+logger.debug("maxW: " + maxW);
 }
 int col = 0;
-for (double w = startW; w <= maxW; w += 0.01) {
+for (int w = startW; w < maxW; w++) {
 if (col == 0) {
-if (wOutUnitIsMeter) {
+out.printf(Locale.US, "%8d", w);
-// For some crazy reason W's at a gauge should be in the
-// unit but W's exported elsewhere should be in Meter
-out.printf(Locale.US, "%5.2f", w);
-} else {
-out.printf(Locale.US, "%8d", (int)Math.round(w * 100d));
-}
 }
 if (w < minW) {
 out.print(EMPTY);
-}
+} else {
-else {
+double actQ = getQ(w, qFunc);
-printQ(out, getQ(w));
+if (Double.isNaN(actQ)) {
+// Can't happen™
+break;
+} else {
+printQ(out, actQ);
+}
 }
 if (++col >= COLUMNS) {
 out.print("\r\n");
 col = 0;

Mercurial > dive4elements > river

comparison artifacts/src/main/java/org/dive4elements/river/exports/ATWriter.java @ 6304:152e1ba8e5bd