dive4elements/river: flys-artifacts/src/main/java/de/intevation/flys/exports/ATWriter.java annotate

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Thu, 23 Jun 2011 15:50:44 +0000
parents
children	39d191f011dc

rev	line source
729 078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	1 package de.intevation.flys.exports;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	2
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	3 import java.io.IOException;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	4 import java.io.Writer;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	5 import java.io.PrintWriter;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	6
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	7 import java.util.Locale;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	8
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	9 import de.intevation.flys.artifacts.model.WQ;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	10
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11 import org.apache.commons.math.analysis.interpolation.SplineInterpolator;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	13 import org.apache.commons.math.analysis.polynomials.PolynomialSplineFunction;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	14
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15 import org.apache.commons.math.ArgumentOutsideDomainException;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17 import org.apache.log4j.Logger;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	18
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	19 public class ATWriter
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	20 {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	21 private static Logger logger = Logger.getLogger(ATWriter.class);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	22
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	23 public static final int COLUMNS = 10;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	24
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	25 public static final String EMPTY = " ";
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	26
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	27 protected double minW;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28 protected double maxW;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29 protected double minQ;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30 protected double maxQ;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	31
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32 protected PolynomialSplineFunction qFunc;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	34 public ATWriter() {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	35 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	36
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37 public ATWriter(WQ wq) throws IllegalArgumentException {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	39 int maxIndex = maxIncreasingWIndex(wq);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	40
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	41 if (maxIndex < 1) { // Only first w can be written out.
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	42 minW = maxW = wq.getW(0);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	43 minQ = maxQ = wq.getQ(0);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	44 return;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	45 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	46
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	47 double [] ws = new double[maxIndex+1];
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	48 double [] qs = new double[ws.length];
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	49
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	50 for (int i = 0; i < ws.length; ++i) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	51 ws[i] = wq.getW(i);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 qs[i] = wq.getQ(i);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	53 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55 SplineInterpolator interpolator = new SplineInterpolator();
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57 qFunc = interpolator.interpolate(ws, qs);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	58
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	59 minW = wq.getW(0);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	60 maxW = wq.getW(maxIndex);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	61 minQ = wq.getQ(0);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	62 maxQ = wq.getQ(maxIndex);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	63 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	64
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	65 public static int maxIncreasingWIndex(WQ wq) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	66
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	67 int N = wq.size();
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69 if (N < 2) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70 return N-1;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	71 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	72
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	73 double last = wq.getW(0);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	74
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	75 for (int i = 1; i < N; ++i) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	76 double current = wq.getW(i);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	77 if (current <= last) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	78 return i-1;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	79 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	80 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	81
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	82 return N-1;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	83 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	84
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	85 public double getQ(double w) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	86
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	87 try {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	88 return qFunc.value(w);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	89 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	90 catch (ArgumentOutsideDomainException aode) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	91 // should not happen
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	92 logger.warn("spline interpolation failed", aode);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	93 return w <= minW ? minQ : maxQ;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	94 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	95 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	96
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	97 protected static void printQ(PrintWriter out, double q) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	98 String format;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	99 if (q < 1) format = " % 8.3f";
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	100 else if (q < 10) format = " % 8.2f";
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	101 else if (q < 100) format = " % 8.1f";
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	102 else {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	103 format = " % 8.0f";
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	104 if (q > 1000) q = Math.rint(q/10)*10;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	105 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	106 out.printf(Locale.US, format, q);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	107 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	108
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	109 public void write(Writer writer) throws IOException {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	110
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	111 PrintWriter out = new PrintWriter(writer);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	112
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	113 double range = COLUMNS;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	114
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	115 double rest = Math.abs(minW % range);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	116
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	117 int startW = (int)(minW - rest);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	118 if (rest >= 1d) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	119 out.printf(Locale.US, "%8d", startW);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	120 int columns = (int)Math.floor(rest);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	121 for (int i = columns; i < COLUMNS; ++i) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	122 out.print(EMPTY);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	123 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	124 for (int i = 0, N = COLUMNS-columns; i < N; ++i) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	125 printQ(out, getQ(startW + (COLUMNS-i-1)));
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	126 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	127 startW += COLUMNS;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	128 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	129
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	130 int col = 0;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	131 for (int w = startW; w <= maxW; ++w) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	132 if (col == 0) out.printf(Locale.US, "%8d", w);
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	133
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	134 printQ(out, getQ(w));
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	135
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	136 if (++col >= COLUMNS) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	137 out.println();
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	138 col = 0;
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	139 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	141
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142 if (col > 0) {
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	143 out.println();
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	144 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	145
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	146 out.flush();
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	147 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148 }
078934923acb New AT writer for results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	149 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :

729

078934923acb New AT writer for results.