Mercurial > dive4elements > river
view flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/MainValuesWFacet.java @ 4798:39885bdfc6fc
Added calculation of the "Umhuellende" to calculation of "W fuer ungleichwertige Abfluesse".
This is done by figuring out the WST columns that imfold the data and then do
simple "gleichwertige" calculations from the start of the interval. This
is too much because only the Qs are needed for the "Umhuellende".
author | Sascha L. Teichmann <teichmann@intevation.de> |
---|---|
date | Sun, 13 Jan 2013 16:18:28 +0100 |
parents | 5b8919ef601d |
children |
line wrap: on
line source
package de.intevation.flys.artifacts.model; import java.util.ArrayList; import java.util.List; import org.apache.log4j.Logger; import de.intevation.artifacts.Artifact; import de.intevation.artifacts.CallContext; import de.intevation.artifacts.DataProvider; import de.intevation.artifactdatabase.state.DefaultFacet; import de.intevation.flys.artifacts.MainValuesArtifact; import de.intevation.flys.artifacts.math.Linear; import de.intevation.flys.jfree.FLYSAnnotation; import de.intevation.flys.jfree.StickyAxisAnnotation; /** * Facet to show Main W Values. */ public class MainValuesWFacet extends DefaultFacet implements FacetTypes { /** Own logger. */ private static Logger logger = Logger.getLogger(MainValuesWFacet.class); /** Do we want MainValues at Gauge (not interpolated)? */ protected boolean isAtGauge; /** Trivial Constructor. */ public MainValuesWFacet(String name, String description, boolean atGauge) { this.description = description; this.name = name; this.index = 0; this.isAtGauge = atGauge; } /** * Set the hit-point in W where a line drawn from the axis would hit the * curve in WQDay (if hit). * Employ linear interpolation. */ protected static void setHitPoint(WQDay wqday, StickyAxisAnnotation annotation) { int idx = 0; float w = annotation.getPos(); boolean wIncreases = wqday.getW(0) < wqday.getW(wqday.size()-1); if (wIncreases) { while (idx < wqday.size() && wqday.getW(idx) < w) { idx++; } } else { idx = wqday.size() -1; while (idx > 0 && wqday.getW(idx) > w) { idx--; } } double day = 0d; int mod = (wIncreases) ? -1 : +1; if (idx != 0 && idx < wqday.size()-1-mod) { day = Linear.linear(w, wqday.getW(idx +mod), wqday.getW(idx), wqday.getDay(idx+mod), wqday.getDay(idx)); annotation.setHitPoint((float)day); } else { logger.debug("StickyAnnotation does not hit wqday curve"); } } /** * Returns the data this facet provides. * * @param artifact the owner artifact. * @param context the CallContext (ignored). * * @return the data. */ @Override public Object getData(Artifact artifact, CallContext context) { MainValuesArtifact mvArtifact = (MainValuesArtifact) artifact; List<NamedDouble> ws = mvArtifact.getMainValuesW(isAtGauge); List<StickyAxisAnnotation> xy = new ArrayList<StickyAxisAnnotation>(); // Find whether a duration curve is on the blackboard. WQDay wqdays = null; List<DataProvider> providers = context. getDataProvider(DurationCurveFacet.BB_DURATIONCURVE); if (providers.size() < 1) { logger.warn("Could not find durationcurve data provider."); } else { wqdays = (WQDay) providers.get(0).provideData( DurationCurveFacet.BB_DURATIONCURVE, null, context); } for (NamedDouble w: ws) { StickyAxisAnnotation annotation = new StickyAxisAnnotation( w.getName(), (float) w.getValue(), StickyAxisAnnotation.SimpleAxis.Y_AXIS); xy.add(annotation); if (wqdays != null) { setHitPoint(wqdays, annotation); } } return new FLYSAnnotation(description, xy); } /** * Create a deep copy of this Facet. * @return a deep copy. */ @Override public MainValuesWFacet deepCopy() { MainValuesWFacet copy = new MainValuesWFacet(this.name, description, this.isAtGauge); copy.set(this); return copy; } } // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :