Mercurial > dive4elements > river
view artifacts/src/main/java/org/dive4elements/river/artifacts/model/WQ.java @ 8098:09725b65955a
Add new and simplyfied SedimentLoadFacet
The SedimentLoadFacet is intended to work with the Measurement stations.
It uses the same mechanismn to access the Mesurement station values
as the calculation does.
SedimentLoadLS values need a different facet that will come soon.
author | Andre Heinecke <andre.heinecke@intevation.de> |
---|---|
date | Fri, 15 Aug 2014 18:27:19 +0200 |
parents | d3a4b0d5bcd8 |
children | 839032ac1523 |
line wrap: on
line source
/* Copyright (C) 2011, 2012, 2013 by Bundesanstalt für Gewässerkunde * Software engineering by Intevation GmbH * * This file is Free Software under the GNU AGPL (>=v3) * and comes with ABSOLUTELY NO WARRANTY! Check out the * documentation coming with Dive4Elements River for details. */ package org.dive4elements.river.artifacts.model; import org.dive4elements.river.utils.DoubleUtil; import gnu.trove.TDoubleArrayList; import java.math.BigDecimal; import java.util.regex.Matcher; import java.util.regex.Pattern; import org.apache.log4j.Logger; public class WQ extends W { public static final Pattern NUMBERS_PATTERN = Pattern.compile("\\D*(\\d++.\\d*)\\D*"); private static Logger log = Logger.getLogger(WQ.class); protected TDoubleArrayList qs; public WQ() { this(""); } public WQ(String name) { super(name); qs = new TDoubleArrayList(); } public WQ(int capacity) { this(capacity, ""); } public WQ(int capacity, String name) { super(capacity, name); qs = new TDoubleArrayList(capacity); } public WQ(double [] qs, double [] ws) { this(qs, ws, ""); } public WQ(double [] qs, double [] ws, String name) { super(name); this.ws = new TDoubleArrayList(ws); this.qs = new TDoubleArrayList(qs); } public WQ(TDoubleArrayList qs, TDoubleArrayList ws, String name) { super(name); this.ws = ws; this.qs = qs; } public Double getRawValue() { if (name == null || name.length() == 0) { // this should never happen return null; } Matcher m = NUMBERS_PATTERN.matcher(name); if (m.matches()) { String raw = m.group(1); try { return Double.valueOf(raw); } catch (NumberFormatException nfe) { // do nothing } } return null; } public void add(double w, double q) { ws.add(w); qs.add(q); } public double getQ(int idx) { return qs.getQuick(idx); } @Override public double [] get(int idx) { return get(idx, new double [2]); } @Override public double [] get(int idx, double [] dst) { dst[0] = ws.getQuick(idx); dst[1] = qs.getQuick(idx); return dst; } public double [] getQs() { return qs.toNativeArray(); } @Override public void removeNaNs() { DoubleUtil.removeNaNs(new TDoubleArrayList [] { ws, qs }); } /** Returns either a modified copy or the same Object with fixed W values. * If a conversion takes place converted is set to true */ public static WQ getFixedWQforExportAtGauge(WQ wq, BigDecimal datum) { if (wq.getReferenceSystem() == W.CENTIMETER_AT_GAUGE) { // Do nothing return wq; } // If we convert we work on a copy to avoid side effects. WQ ret = new WQ(wq.size(), wq.getName()); ret.setReferenceSystem(W.CENTIMETER_AT_GAUGE); // When we convert and have a datum we have a calculated // result at a gauge so we must subtract the datum. double subtractDatum = datum == null ? 0 : datum.doubleValue(); double [] data = new double[8]; for (int i = 0, WQ = wq.size(); i < WQ; i++) { wq.get(i, data); ret.add((data[0] - subtractDatum)*100d, data[1]); } log.debug("Converted W values to centimeter and substracted: " + subtractDatum); return ret; } } // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :