Mercurial > dive4elements > river
view flys-artifacts/src/main/java/de/intevation/flys/artifacts/WINFOArtifact.java @ 395:aae8f327425e
'berechnete Abflusskurve': generate ws from qs
flys-artifacts/trunk@1822 c6561f87-3c4e-4783-a992-168aeb5c3f6f
author | Sascha L. Teichmann <sascha.teichmann@intevation.de> |
---|---|
date | Wed, 04 May 2011 13:51:25 +0000 |
parents | befedd7629d5 |
children | eb22ffe4d74c |
line wrap: on
line source
package de.intevation.flys.artifacts; import java.util.List; import java.util.Vector; import org.w3c.dom.Document; import org.w3c.dom.Element; import org.w3c.dom.Node; import org.apache.log4j.Logger; import de.intevation.artifacts.ArtifactNamespaceContext; import de.intevation.artifacts.CallContext; import de.intevation.artifactdatabase.ProtocolUtils; import de.intevation.artifactdatabase.state.Output; import de.intevation.artifactdatabase.state.State; import de.intevation.artifactdatabase.state.StateEngine; import de.intevation.artifactdatabase.transition.TransitionEngine; import de.intevation.artifacts.common.utils.XMLUtils; import de.intevation.flys.model.Gauge; import de.intevation.flys.model.River; import de.intevation.flys.artifacts.states.DefaultState; import de.intevation.flys.artifacts.context.FLYSContext; import de.intevation.flys.artifacts.model.MainValuesFactory; import de.intevation.flys.artifacts.model.WQDay; import de.intevation.flys.artifacts.model.WQKms; import de.intevation.flys.artifacts.model.WstValueTable; /** * The default WINFO artifact. * * @author <a href="mailto:ingo.weinzierl@intevation.de">Ingo Weinzierl</a> */ public class WINFOArtifact extends FLYSArtifact { /** The logger for this class */ private static Logger logger = Logger.getLogger(WINFOArtifact.class); /** The name of the artifact.*/ public static final String ARTIFACT_NAME = "winfo"; /** XPath */ public static final String XPATH_STATIC_UI ="/art:result/art:ui/art:static"; /** * The default constructor. */ public WINFOArtifact() { } /** * This method returns a description of this artifact. * * @param data Some data. * @param context The CallContext. * * @return the description of this artifact. */ public Document describe(Document data, CallContext context) { logger.debug("Describe: the current state is: " + getCurrentStateId()); FLYSContext flysContext = null; if (context instanceof FLYSContext) { flysContext = (FLYSContext) context; } else { flysContext = (FLYSContext) context.globalContext(); } StateEngine stateEngine = (StateEngine) flysContext.get( FLYSContext.STATE_ENGINE_KEY); TransitionEngine transitionEngine = (TransitionEngine) flysContext.get( FLYSContext.TRANSITION_ENGINE_KEY); List<State> reachable = transitionEngine.getReachableStates( this, getCurrentState(context), stateEngine); Document description = XMLUtils.newDocument(); XMLUtils.ElementCreator creator = new XMLUtils.ElementCreator( description, ArtifactNamespaceContext.NAMESPACE_URI, ArtifactNamespaceContext.NAMESPACE_PREFIX); Element root = ProtocolUtils.createRootNode(creator); description.appendChild(root); State current = getCurrentState(context); ProtocolUtils.appendDescribeHeader(creator, root, identifier(), hash()); ProtocolUtils.appendState(creator, root, current); ProtocolUtils.appendReachableStates(creator, root, reachable); Element name = ProtocolUtils.createArtNode( creator, "name", new String[] { "value" }, new String[] { getName() }); Element ui = ProtocolUtils.createArtNode( creator, "ui", null, null); Element staticUI = ProtocolUtils.createArtNode( creator, "static", null, null); Element outs = ProtocolUtils.createArtNode( creator, "outputmodes", null, null); appendOutputModes(description, outs, context, identifier()); appendStaticUI(description, staticUI, context, identifier()); Element dynamic = current.describe( this, description, root, context, identifier()); if (dynamic != null) { ui.appendChild(dynamic); } ui.appendChild(staticUI); root.appendChild(name); root.appendChild(ui); root.appendChild(outs); return description; } /** * Returns the name of the concrete artifact. * * @return the name of the concrete artifact. */ public String getName() { return ARTIFACT_NAME; } protected void appendOutputModes( Document doc, Element outs, CallContext context, String uuid) { Vector<String> stateIds = getPreviousStateIds(); XMLUtils.ElementCreator creator = new XMLUtils.ElementCreator( doc, ArtifactNamespaceContext.NAMESPACE_URI, ArtifactNamespaceContext.NAMESPACE_PREFIX); FLYSContext flysContext = getFlysContext(context); StateEngine engine = (StateEngine) flysContext.get( FLYSContext.STATE_ENGINE_KEY); for (String stateId: stateIds) { logger.debug("Append output modes for state: " + stateId); State state = engine.getState(stateId); List<Output> list = state.getOutputs(); if (list == null || list.size() == 0) { continue; } ProtocolUtils.appendOutputModes(creator, outs, list); } try { DefaultState cur = (DefaultState) getCurrentState(context); if (cur.validate(this, context)) { List<Output> list = cur.getOutputs(); if (list != null && list.size() > 0) { logger.debug( "Append output modes for state: " + cur.getID()); ProtocolUtils.appendOutputModes(creator, outs, list); } } } catch (IllegalArgumentException iae) { // state is not valid, so we do not append its outputs. } } /** * This method appends the static data - that has already been inserted by * the user - to the static node of the DESCRIBE document. * * @param doc The document. * @param ui The root node. * @param context The CallContext. * @param uuid The identifier of the artifact. */ protected void appendStaticUI( Document doc, Node ui, CallContext context, String uuid) { Vector<String> stateIds = getPreviousStateIds(); FLYSContext flysContext = getFlysContext(context); StateEngine engine = (StateEngine) flysContext.get( FLYSContext.STATE_ENGINE_KEY); for (String stateId: stateIds) { logger.debug("Append static data for state: " + stateId); DefaultState state = (DefaultState) engine.getState(stateId); state = (DefaultState) fillState(state); ui.appendChild(state.describeStatic(doc, ui, context, uuid)); } } // // METHODS FOR RETRIEVING COMPUTED DATA FOR DIFFERENT CHART TYPES // /** * Returns the data that is computed by a waterlevel computation. * * @return an array of data triples that consist of W, Q and Kms. */ public WQKms[] getWaterlevelData() throws NullPointerException { logger.debug("WINFOArtifact.getWaterlevelData"); River river = getRiver(); if (river == null) { throw new NullPointerException("No river selected."); } double[] kms = getKms(); if (kms == null) { throw new NullPointerException("No Kms selected."); } double[] qs = getQs(); if (qs == null) { logger.debug("Determine Q values based on a set of W values."); double[] ws = getWs(); qs = getQsForWs(ws); } WstValueTable wst = WstValueTable.getTable(river); if (wst == null) { throw new NullPointerException("No Wst found for selected river."); } // TODO Introduce a caching mechanism here! return computeWaterlevelData(kms, qs, wst); } /** * Computes the data of a waterlevel computation based on the interpolation * in WstValueTable. * * @param kms The kilometer values. * @param qa The discharge values. * @param wst The WstValueTable used for the interpolation. * * @return an array of data triples that consist of W, Q and Kms. */ public static WQKms[] computeWaterlevelData( double[] kms, double[] qs, WstValueTable wst) { logger.info("WINFOArtifact.computeWaterlevelData"); WQKms[] wqkms = new WQKms[qs.length]; for (int i = 0; i < wqkms.length; i++) { wqkms[i] = new WQKms(kms.length); } double [] interpolatedW = new double[qs.length]; for (double km: kms) { wst.interpolateW(km, qs, interpolatedW); // TODO Modify the interpolation to return a better formed data // structure. for (int i = 0; i < interpolatedW.length; i++) { wqkms[i].add(interpolatedW[i], qs[i], km); } } return wqkms; } /** * Returns the data that is computed by a duration curve computation. * * @return the data computed by a duration curve computation. */ public WQDay getDurationCurveData() throws NullPointerException { logger.debug("WINFOArtifact.getDurationCurveData"); River r = getRiver(); if (r == null) { throw new NullPointerException("Cannot determine river."); } Gauge g = getGauge(); if (g == null) { throw new NullPointerException("Cannot determine gauge."); } double[] locations = getLocations(); if (locations == null) { throw new NullPointerException("Cannot determine location."); } WstValueTable wst = WstValueTable.getTable(r); if (wst == null) { throw new NullPointerException("No Wst found for selected river."); } // TODO Introduce a caching mechanism here! return computeDurationCurveData(g, wst, locations[0]); } /** * Computes the data used to create duration curves. * * @param gauge The selected gauge. * @param location The selected location. * * @return the computed data. */ public static WQDay computeDurationCurveData( Gauge gauge, WstValueTable wst, double location) { logger.info("WINFOArtifact.computeDurationCurveData"); Object[] obj = MainValuesFactory.getDurationCurveData(gauge); int[] days = (int[]) obj[0]; double[] qs = (double[]) obj[1]; double[] interpolatedW = new double[qs.length]; interpolatedW = wst.interpolateW(location, qs, interpolatedW); WQDay wqday = new WQDay(qs.length); for (int i = 0; i < days.length; i++) { wqday.add(days[i], interpolatedW[i], qs[i]); } return wqday; } /** * Returns the data that is computed by a discharge curve computation. * * @return the data computed by a discharge curve computation. */ public WQKms getComputedDischargeCurveData() throws NullPointerException { logger.debug("WINFOArtifact.getComputedDischargeCurveData"); River r = getRiver(); if (r == null) { throw new NullPointerException("Cannot determine river."); } double[] locations = getLocations(); if (locations == null) { throw new NullPointerException("Cannot determine location."); } WstValueTable wst = WstValueTable.getTable(r); if (wst == null) { throw new NullPointerException("No Wst found for selected river."); } // TODO Introduce a caching mechanism here! return computeDischargeCurveData(wst, locations[0]); } /** * Computes the data used to create computed discharge curves. * * @param wst The WstValueTable that is used for the interpolation. * @param location The location where the computation should be based on. * * @return an object that contains tuples of W/Q values at the specified * location. */ public static WQKms computeDischargeCurveData( WstValueTable wst, double location) { logger.info("WINFOArtifact.computeDischargeCurveData"); double[][] wqs = wst.interpolateWQ(location); if (wqs == null) { logger.error("Cannot compute discharge curve data."); return null; } double[] ws = wqs[0]; double[] qs = wqs[1]; WQKms wqkms = new WQKms(ws.length); for (int i = 0; i < ws.length; i++) { wqkms.add(ws[i], qs[i], location); } return wqkms; } } // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :