dive4elements/river: flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/fixings/FixCalculation.java annotate

annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/fixings/FixCalculation.java @ 2783:cbb513a8f548

FixA: km indexed datastructure for delta w/t results. flys-artifacts/trunk@4522 c6561f87-3c4e-4783-a992-168aeb5c3f6f

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Fri, 25 May 2012 13:49:27 +0000
parents	c1f2e792704a
children	f9e62e76ee1f

rev	line source
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	1 package de.intevation.flys.artifacts.model.fixings;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	2
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	3 import de.intevation.flys.artifacts.FixationArtifactAccess;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	4
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	5 import de.intevation.flys.artifacts.math.fitting.Function;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	6 import de.intevation.flys.artifacts.math.fitting.FunctionFactory;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	7
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	8 import de.intevation.flys.artifacts.model.Calculation;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	9 import de.intevation.flys.artifacts.model.CalculationResult;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	10 import de.intevation.flys.artifacts.model.FixingsColumn;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11 import de.intevation.flys.artifacts.model.FixingsColumnFactory;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12 import de.intevation.flys.artifacts.model.FixingsOverview;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	13 import de.intevation.flys.artifacts.model.FixingsOverviewFactory;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	14 import de.intevation.flys.artifacts.model.Parameters;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	16 import de.intevation.flys.artifacts.model.FixingsOverview.Range;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17 import de.intevation.flys.artifacts.model.FixingsOverview.Fixing;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	18 import de.intevation.flys.artifacts.model.FixingsOverview.IdFilter;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	19 import de.intevation.flys.artifacts.model.FixingsOverview.Fixing.Filter;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	20 import de.intevation.flys.artifacts.model.FixingsOverview.AndFilter;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	21 import de.intevation.flys.artifacts.model.FixingsOverview.DateRangeFilter;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	22 import de.intevation.flys.artifacts.model.FixingsOverview.SectorRangeFilter;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	23
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	24 import de.intevation.flys.utils.DoubleUtil;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	25
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	26 import java.util.ArrayList;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	27 import java.util.List;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	28 import java.util.Date;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30 import org.apache.commons.math.MathException;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	31
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32 import org.apache.commons.math.optimization.fitting.CurveFitter;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	34 import org.apache.commons.math.optimization.general.LevenbergMarquardtOptimizer;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	35
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	36 import org.apache.log4j.Logger;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38 public class FixCalculation
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	39 extends Calculation
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	40 {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	41 private static Logger logger = Logger.getLogger(FixCalculation.class);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	42
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	43 protected String river;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	44 protected double from;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	45 protected double to;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	46 protected double step;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	47 protected boolean preprocessing;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	48 protected String function;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	49 protected int [] events;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	50 protected long [][] analysisPeriods;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	51 protected int qSectorStart;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	52 protected int qSectorEnd;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	53
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54 public FixCalculation() {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57 public FixCalculation(FixationArtifactAccess access) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	58
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	59 String river = access.getRiver();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	60 Double from = access.getFrom();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	61 Double to = access.getTo();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	62 Double step = access.getStep();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	63 String function = access.getFunction();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	64 int [] events = access.getEvents();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	65 long [][] analysisPeriods = access.getAnalysisPeriods();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	66 Integer qSectorStart = access.getQSectorStart();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	67 Integer qSectorEnd = access.getQSectorEnd();
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69 if (river == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	71 addProblem("fix.missing.river");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	72 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	73
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	74 if (from == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	75 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	76 addProblem("fix.missing.from");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	77 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	78
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	79 if (to == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	80 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	81 addProblem("fix.missing.to");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	82 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	83
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	84 if (step == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	85 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	86 addProblem("fix.missing.step");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	87 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	88
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	89 if (function == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	90 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	91 addProblem("fix.missing.function");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	92 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	93
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	94 if (events == null \|\| events.length < 1) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	95 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	96 addProblem("fix.missing.events");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	97 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	98
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	99 if (analysisPeriods == null \|\| analysisPeriods.length < 1) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	100 // TODO: i18n
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	101 addProblem("fix.missing.analysis.periods");
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	102 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	103
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	104 if (qSectorStart == null) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	105 // TODO: i18n
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	106 addProblem("fix.missing.qstart.sector");
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	107 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	108
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	109 if (qSectorEnd == null) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	110 // TODO: i18n
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	111 addProblem("fix.missing.qend.sector");
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	112 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	113
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	114 if (!hasProblems()) {
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	115 this.river = river;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	116 this.from = from;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	117 this.to = to;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	118 this.step = step;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	119 this.function = function;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	120 this.events = events;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	121 this.analysisPeriods = analysisPeriods;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	122 this.qSectorStart = qSectorStart;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	123 this.qSectorEnd = qSectorEnd;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	124 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	125 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	126
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	127 public CalculationResult calculate() {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	128
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	129 FixingsOverview overview =
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	130 FixingsOverviewFactory.getOverview(river);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	131
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	132 if (overview == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	133 addProblem("fix.no.overview.available");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	134 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	135
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	136 Function func = FunctionFactory.getInstance()
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	137 .getFunction(function);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	138
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	139 if (func == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	141 addProblem("fix.invalid.function.name");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	143
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	144 if (hasProblems()) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	145 return new CalculationResult(this);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	146 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	147
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148 FixingsColumnFactory fcf = FixingsColumnFactory.getInstance();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	149
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	150 List<FixingsColumn> dataColumns =
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	151 new ArrayList<FixingsColumn>(events.length);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	152
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	153 for (int eventId: events) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	154 IdFilter idFilter = new IdFilter(eventId);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	155
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	156 List<Fixing.Column> columns = overview.filter(null, idFilter);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	157 if (columns.isEmpty()) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	158 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159 addProblem("fix.missing.column", eventId);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	160 continue;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	161 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	162 FixingsColumn dataColumn = fcf.getColumnData(columns.get(0));
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	163 if (dataColumn == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	164 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	165 addProblem("fix.cannot.load.data", eventId);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	166 continue;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	167 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	168 dataColumns.add(dataColumn);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	170
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	171 if (dataColumns.size() < 2) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	172 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	173 addProblem("fix.too.less.data.columns");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174 return new CalculationResult(this);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	175 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	177 double [] kms = DoubleUtil.explode(from, to, step);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	178
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	179 double [] ws = new double[dataColumns.size()];
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	180 double [] qs = new double[ws.length];
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	181
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182 String [] parameterNames = func.getParameterNames();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	183
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	184 Parameters results =
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	185 new Parameters(createColumnNames(parameterNames));
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	186
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	187 boolean invalid = false;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	188
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	189 for (int i = 0; i < kms.length; ++i) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	190 double km = kms[i];
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	191
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	192 for (int j = 0; j < ws.length; ++j) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	193 FixingsColumn column = dataColumns.get(j);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	194 qs[j] = column.getQ(km);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	195 boolean interpolated = column.getW(km, ws, j);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	196 // TODO: mark as interpolated.
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	197 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	198
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	199 // TODO: Do preprocessing here!
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	200 double [] parameters = fit(func, km, ws, qs);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	201 if (parameters == null) { // Problems are reported already.
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	202 continue;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	203 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	204
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	205 int row = results.newRow();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	206
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	207 results.set(row, "km", km);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	208 for (int j = 0; j < parameters.length; ++j) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	209 if (Double.isNaN(parameters[j])) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	210 invalid = true;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	211 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	212 else {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	213 results.set(row, parameterNames[j], parameters[j]);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	215 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	216 // TODO: Calculate statistics, too!
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	217 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	218
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	219 if (invalid) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	220 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	221 addProblem("fix.invalid.values");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	222 results.removeNaNs();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	223 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	224
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	225 // Calculate Delta W/t
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	226 DeltaWTsKM deltaWTsKM = calculateDeltaWTs(
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	227 func,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	228 overview,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	229 results);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	230
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	231 return new CalculationResult(deltaWTsKM, this);
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	232 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	233
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	234 public DeltaWTsKM calculateDeltaWTs(
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	235 Function function,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	236 FixingsOverview overview,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	237 Parameters results
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	238 ) {
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	239 DeltaWTsKM deltaWTsKM = new DeltaWTsKM(results.size());
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	240
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	241 Column [][] analysisColumns = getAnalysisColumns(overview);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	242
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	243 int [] parameterIndices =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	244 results.columnIndices(function.getParameterNames());
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	245
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	246 double [] parameterValues =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	247 new double[parameterIndices.length];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	248
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	249 double [] ow = new double[1];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	250
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	251 for (int i = 0, N = results.size(); i < N; ++i) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	252 double km = results.get(i, "km");
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	253 results.get(i, parameterIndices, parameterValues);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	254
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	255 DeltaWTsKM.KM dwtkm = new DeltaWTsKM.KM(km);
cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	256 deltaWTsKM.add(dwtkm);
cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	257
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	258 // This is the paraterized function for a given km.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	259 de.intevation.flys.artifacts.math.Function instance =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	260 function.instantiate(parameterValues);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	261
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	262 // Evaluate all columns for all analysis periods.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	263 for (int j = 0; j < analysisColumns.length; ++j) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	264 Column [] periodColumns = analysisColumns[j];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	265
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	266 for (int k = 0; k < periodColumns.length; ++k) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	267 Column pc = periodColumns[k];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	268
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	269 // Q from real data.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	270 double q = pc.data.getQ(km);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	271 if (Double.isNaN(q)) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	272 continue;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	273 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	274
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	275 // Calculate W from function.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	276 double nw = instance.value(q);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	277 if (Double.isNaN(nw)) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	278 continue;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	279 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	280
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	281 // W from real data.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	282 pc.data.getW(km, ow);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	283
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	284 if (Double.isNaN(ow[0])) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	285 continue;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	286 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	287
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	288 double deltaW = ow[0] - nw;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	289
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	290 DeltaWT deltaWT = new DeltaWT(
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	291 deltaW,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	292 pc.meta.getStartTime(),
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	293 pc.meta.getDescription());
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	294
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	295 dwtkm.add(deltaWT);
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	296 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	297 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	298 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	299
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	300 return deltaWTsKM;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	301 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	302
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	303 /** Helper class to bundle the meta information of a column
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	304 * and the real data.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	305 */
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	306 protected static class Column {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	307
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	308 protected Fixing.Column meta;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	309 protected FixingsColumn data;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	310
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	311 public Column() {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	312 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	313
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	314 public Column(Fixing.Column meta, FixingsColumn data) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	315 this.meta = meta;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	316 this.data = data;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	317 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	318 } // class Column
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	319
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	320 /** Fetch meta and data columns for analysis periods. */
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	321 protected Column [][] getAnalysisColumns(FixingsOverview overview) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	322
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	323 Column columns [][] = new Column[analysisPeriods.length][];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	324
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	325 Range range = new Range(from, to);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	326 SectorRangeFilter sectorRangeFilter =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	327 new SectorRangeFilter(qSectorStart, qSectorEnd);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	328
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	329 FixingsColumnFactory fcf = FixingsColumnFactory.getInstance();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	330
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	331 for (int i = 0; i < columns.length; ++i) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	332
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	333 // Construct filter for period.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	334 long [] period = analysisPeriods[i];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	335
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	336 AndFilter filter = new AndFilter();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	337
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	338 DateRangeFilter dateRangeFilter =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	339 new DateRangeFilter(
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	340 new Date(period[0]),
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	341 new Date(period[1]));
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	342
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	343 filter.add(dateRangeFilter);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	344 filter.add(sectorRangeFilter);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	345
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	346 List<Fixing.Column> metaCols = overview.filter(range, filter);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	347
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	348 ArrayList<Column> cols = new ArrayList<Column>(metaCols.size());
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	349
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	350 // Only use columns which have data.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	351 for (Fixing.Column meta: metaCols) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	352 FixingsColumn data = fcf.getColumnData(meta);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	353 if (data != null) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	354 cols.add(new Column(meta, data));
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	355 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	356 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	357 columns[i] = cols.toArray(new Column[cols.size()]);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	358 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	359
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	360 return columns;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	361 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	362
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	363 protected static String [] createColumnNames(String [] parameters) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	364 String [] result = new String[parameters.length + 1];
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	365 result[0] = "km";
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	366 // TODO: Add statistic columns, too.
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	367 System.arraycopy(parameters, 0, result, 1, parameters.length);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	368 return result;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	369 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	370
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	371 protected double [] fit(
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	372 Function function,
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	373 double km,
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	374 double [] ws,
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	375 double [] qs
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	376 ) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	377 LevenbergMarquardtOptimizer lmo = new LevenbergMarquardtOptimizer();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	378 CurveFitter cf = new CurveFitter(lmo);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	379
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	380 boolean missingWs = false;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	381 boolean missingQs = false;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	382
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	383 for (int i = 0; i < ws.length; ++i) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	384 boolean ignore = false;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	385 if (Double.isNaN(ws[i])) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	386 ignore = true;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	387 if (!missingWs) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	388 missingWs = true;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	389 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	390 addProblem(km, "fix.missing.w");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	391 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	392 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	393 if (Double.isNaN(qs[i])) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	394 ignore = true;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	395 if (!missingQs) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	396 missingQs = true;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	397 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	398 addProblem(km, "fix.missing.q");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	399 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	400 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	401 if (!ignore) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	402 cf.addObservedPoint(ws[i], qs[i]);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	403 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	404 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	405
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	406 try {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	407 return cf.fit(function, function.getInitialGuess());
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	408 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	409 catch (MathException me) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	410 addProblem(km, "fix.fitting.failed");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	411 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	412
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	413 return null;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	414 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	415 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	416 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :

Mercurial > dive4elements > river

annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/fixings/FixCalculation.java @ 2783:cbb513a8f548