dive4elements/river: flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/sq/Fitting.java annotate

annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/sq/Fitting.java @ 3552:1df6984628c3

S/Q: Extented the result data model of the S/Q calculation to store the curve coefficients for each iteration step of the outlier elimination. flys-artifacts/trunk@5146 c6561f87-3c4e-4783-a992-168aeb5c3f6f

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Fri, 27 Jul 2012 12:36:09 +0000
parents	49fe2ed03c12
children	8d0f06b76e09

rev	line source
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	1 package de.intevation.flys.artifacts.model.sq;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	2
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	3 import de.intevation.flys.artifacts.math.fitting.Function;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	4
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	5 import java.util.ArrayList;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	6 import java.util.Iterator;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	7 import java.util.List;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	8
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	9 import org.apache.commons.math.MathException;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	10
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11 import org.apache.commons.math.optimization.fitting.CurveFitter;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	13 import org.apache.commons.math.optimization.general.LevenbergMarquardtOptimizer;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	14
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15 import org.apache.log4j.Logger;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17 public class Fitting
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	18 implements Outlier.Callback
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	19 {
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	20 private static Logger log = Logger.getLogger(Fitting.class);
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	21
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	22 public interface Callback {
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	23
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	24 void afterIteration(
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	25 double [] parameters,
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	26 SQ [] measurements,
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	27 SQ [] outliers,
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	28 double standardDeviation,
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	29 double chiSqr);
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	30 } // interfacte
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	31
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32 protected Function function;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	34 protected double [] coeffs;
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	35
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	36 protected de.intevation.flys.artifacts.math.Function instance;
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	37
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	38 protected List<SQ> remainings;
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	39 protected List<SQ> outliers;
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	40
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	41 protected double standardDeviation;
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	42 protected double stdDevFactor;
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	43 protected double chiSqr;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	44
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	45 protected Callback callback;
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	46
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	47 public Fitting() {
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	48 remainings = new ArrayList<SQ>();
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	49 outliers = new ArrayList<SQ>();
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	50 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	51
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 public Fitting(Function function, double stdDevFactor) {
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	53 this();
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54 this.function = function;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55 this.stdDevFactor = stdDevFactor;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	58 public Function getFunction() {
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	59 return function;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	60 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	61
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	62 public void setFunction(Function function) {
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	63 this.function = function;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	64 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	65
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	66 public double getStdDevFactor() {
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	67 return stdDevFactor;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70 public void setStdDevFactor(double stdDevFactor) {
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	71 this.stdDevFactor = stdDevFactor;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	72 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	73
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	74 @Override
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	75 public void initialize(Iterator<SQ> good) throws MathException {
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	76
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	77 LevenbergMarquardtOptimizer lmo =
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	78 new LevenbergMarquardtOptimizer();
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	79
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	80 CurveFitter cf = new CurveFitter(lmo);
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	81 while (good.hasNext()) {
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	82 SQ sq = good.next();
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	83 cf.addObservedPoint(sq.getQ(), sq.getS());
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	84 }
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	85
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	86 coeffs = cf.fit(
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	87 function, function.getInitialGuess());
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	88
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	89 instance = function.instantiate(coeffs);
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	90
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	91 chiSqr = lmo.getChiSquare();
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	92
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	93 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	94
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	95 @Override
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	96 public double eval(SQ sq) {
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	97 double s = instance.value(sq.q);
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	98 return sq.s - s;
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	99 }
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	100
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	101 @Override
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	102 public void outlier(SQ sq) {
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	103 outliers.add(sq);
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	104 }
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	105
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	106 @Override
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	107 public void remaining(SQ sq) {
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	108 remainings.add(sq);
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	109 }
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	110
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	111 @Override
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	112 public void standardDeviation(double standardDeviation) {
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	113 this.standardDeviation = standardDeviation;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	114 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	115
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	116 @Override
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	117 public void iterationFinished() {
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	118 if (log.isDebugEnabled()) {
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	119 log.debug("iterationFinished ----");
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	120 log.debug(" num remainings: " + remainings.size());
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	121 log.debug(" num outliers: " + outliers.size());
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	122 log.debug(" standardDeviation: " + standardDeviation);
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	123 log.debug(" Chi^2: " + chiSqr);
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	124 log.debug("---- iterationFinished");
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	125 }
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	126 callback.afterIteration(
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	127 coeffs,
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	128 remainings.toArray(new SQ[remainings.size()]),
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	129 outliers.toArray(new SQ[outliers.size()]),
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	130 standardDeviation,
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	131 chiSqr);
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	132 remainings.clear();
1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	133 outliers.clear();
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	134 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	135
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	136 protected static final List<SQ> onlyValid(List<SQ> sqs) {
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	137
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	138 List<SQ> good = new ArrayList<SQ>(sqs.size());
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	139
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140 for (SQ sq: sqs) {
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	141 if (sq.isValid()) {
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142 good.add(sq);
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	143 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	144 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	145
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	146 return good;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	147 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	149 public boolean fit(List<SQ> sqs, Callback callback) {
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	150
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	151 sqs = onlyValid(sqs);
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	152
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	153 if (sqs.size() < 2) {
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	154 log.warn("Too less points for fitting.");
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	155 return false;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	156 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	157
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	158 this.callback = callback;
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	160 try {
3552 1df6984628c3 S/Q: Extented the result data model of the S/Q calculation to Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3190 diff changeset	161 Outlier.detectOutliers(this, sqs, stdDevFactor);
3188 1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	162 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	163 catch (MathException me) {
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	164 log.warn(me);
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	165 return false;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	166 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	167
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	168 return true;
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	170 }
1e46ced2bb57 SQ: Added fitting shell for SQ curves. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	171 // 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/sq/Fitting.java @ 3552:1df6984628c3