Mercurial > dive4elements > river
view flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/sq/SQOutlierCurveFacet.java @ 4357:11cf29eb002c
FacetTypes: New type: exreme_wq_curve_base (for the Stuetz-Q-Poits).
author | Felix Wolfsteller <felix.wolfsteller@intevation.de> |
---|---|
date | Fri, 02 Nov 2012 11:58:11 +0100 |
parents | 1df6984628c3 |
children |
line wrap: on
line source
package de.intevation.flys.artifacts.model.sq; import de.intevation.artifactdatabase.state.Facet; import de.intevation.artifacts.Artifact; import de.intevation.artifacts.CallContext; import de.intevation.flys.artifacts.FLYSArtifact; import de.intevation.flys.artifacts.math.fitting.Function; import de.intevation.flys.artifacts.math.fitting.FunctionFactory; import de.intevation.flys.artifacts.model.CalculationResult; import de.intevation.flys.artifacts.model.DataFacet; import de.intevation.flys.artifacts.model.FacetTypes; import de.intevation.flys.artifacts.model.Parameters; import de.intevation.flys.artifacts.states.DefaultState.ComputeType; import org.apache.log4j.Logger; public class SQOutlierCurveFacet extends DataFacet implements FacetTypes { private static final Logger log = Logger.getLogger(SQOutlierCurveFacet.class); public static final String FUNCTION = "sq-pow"; public static final int BITMASK_ITERATION = (1 << 16) - 1; private int fractionIdx; public SQOutlierCurveFacet() { } public SQOutlierCurveFacet( int idx, int fractionIdx, String name, String description, String hash, String stateId ) { super(idx, name, description, ComputeType.ADVANCE, hash, stateId); this.fractionIdx = fractionIdx; } @Override public Object getData(Artifact artifact, CallContext context) { log.debug("SQOutlierCurveFacet.getData"); if (!(artifact instanceof FLYSArtifact)) { return null; } FLYSArtifact flys = (FLYSArtifact) artifact; CalculationResult res = (CalculationResult) flys.compute( context, ComputeType.ADVANCE, false); int idx = this.index >> 16; int iter = this.index & BITMASK_ITERATION; SQResult[] results = (SQResult[]) res.getData(); SQFractionResult result = results[idx].getFraction(fractionIdx); Parameters params = result.getParameters(iter); if (params == null) { return null; } Function func = FunctionFactory.getInstance().getFunction(FUNCTION); String [] paramNames = func.getParameterNames(); double [] coeffs = params.get(0, paramNames); if (log.isDebugEnabled()) { for (int i = 0; i < paramNames.length; i++) { log.debug("retrieved parameter " + paramNames[i] + " = " + coeffs[i]); } } de.intevation.flys.artifacts.math.Function mf = func.instantiate(coeffs); double [] extent = result.getQExtent(); return new SQFunction(mf, extent[0], extent[1]); } @Override public Facet deepCopy() { SQOutlierCurveFacet copy = new SQOutlierCurveFacet(); copy.set(this); return copy; } } // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :