Mercurial > dive4elements > river
diff artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/BedloadDiameterResult.java @ 8581:073ea4bcea58
(issue1755) Interpolate BedQuality Results
This adds an interpolation function to each various bedQuality result class.
Imho this is ok as the interpolation function can be seen as part of the result.
The interpolation function is initalized on first use and can be accessed through
get.*Interpol functions.
author | Andre Heinecke <andre.heinecke@intevation.de> |
---|---|
date | Mon, 16 Mar 2015 15:36:38 +0100 |
parents | af13ceeba52a |
children |
line wrap: on
line diff
--- a/artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/BedloadDiameterResult.java Mon Mar 16 14:37:30 2015 +0100 +++ b/artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/BedloadDiameterResult.java Mon Mar 16 15:36:38 2015 +0100 @@ -9,13 +9,22 @@ package org.dive4elements.river.artifacts.model.minfo; import org.dive4elements.river.artifacts.model.DateRange; +import org.dive4elements.river.utils.DoubleUtil; + import gnu.trove.TDoubleArrayList; +import org.apache.commons.math.analysis.polynomials.PolynomialSplineFunction; + +import org.apache.commons.math.ArgumentOutsideDomainException; public class BedloadDiameterResult extends BedQualityDiameterResult { protected TDoubleArrayList diameter; + protected PolynomialSplineFunction interpol; + + /** Set to true if this result can't be interpolated.*/ + protected boolean nonInterpolResult; public BedloadDiameterResult( String type, @@ -25,6 +34,8 @@ ) { super (type, km); this.diameter = diameter; + interpol = null; + nonInterpolResult = false; } public double getDiameter(int ndx) { @@ -41,6 +52,25 @@ return Double.NaN; } + public double getDiameterInterpol(double km) { + if (nonInterpolResult) { + return Double.NaN; + } + if (interpol == null) { + interpol = DoubleUtil.getLinearInterpolator(kms, diameter); + if (interpol == null) { + nonInterpolResult = true; + return Double.NaN; + } + } + try { + return interpol.value(km); + } catch (ArgumentOutsideDomainException e) { + /* This is expected for many results. */ + return Double.NaN; + } + } + public double[][] getDiameterData() { return new double[][] { kms.toNativeArray(),