Mercurial > dive4elements > river
comparison artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/BedDiameterResult.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 |
comparison
equal
deleted
inserted
replaced
| 8580:d9f038b8e2ce | 8581:073ea4bcea58 |
|---|---|
| 6 * documentation coming with Dive4Elements River for details. | 6 * documentation coming with Dive4Elements River for details. |
| 7 */ | 7 */ |
| 8 | 8 |
| 9 package org.dive4elements.river.artifacts.model.minfo; | 9 package org.dive4elements.river.artifacts.model.minfo; |
| 10 | 10 |
| 11 import org.dive4elements.river.utils.DoubleUtil; | |
| 12 | |
| 11 import gnu.trove.TDoubleArrayList; | 13 import gnu.trove.TDoubleArrayList; |
| 12 | 14 |
| 15 import org.apache.commons.math.analysis.polynomials.PolynomialSplineFunction; | |
| 16 | |
| 17 import org.apache.commons.math.ArgumentOutsideDomainException; | |
| 13 | 18 |
| 14 public class BedDiameterResult | 19 public class BedDiameterResult |
| 15 extends BedQualityDiameterResult | 20 extends BedQualityDiameterResult |
| 16 { | 21 { |
| 17 protected TDoubleArrayList diameterCap; | 22 protected TDoubleArrayList diameterCap; |
| 18 protected TDoubleArrayList diameterSub; | 23 protected TDoubleArrayList diameterSub; |
| 24 | |
| 25 protected PolynomialSplineFunction interpolSub; | |
| 26 protected PolynomialSplineFunction interpolCap; | |
| 27 protected boolean nonInterpolSub; | |
| 28 protected boolean nonInterpolCap; | |
| 19 | 29 |
| 20 public BedDiameterResult ( | 30 public BedDiameterResult ( |
| 21 String type, | 31 String type, |
| 22 TDoubleArrayList diameterCap, | 32 TDoubleArrayList diameterCap, |
| 23 TDoubleArrayList diameterSub, | 33 TDoubleArrayList diameterSub, |
| 24 TDoubleArrayList km | 34 TDoubleArrayList km |
| 25 ) { | 35 ) { |
| 26 super(type, km); | 36 super(type, km); |
| 27 this.diameterCap = diameterCap; | 37 this.diameterCap = diameterCap; |
| 28 this.diameterSub = diameterSub; | 38 this.diameterSub = diameterSub; |
| 39 interpolSub = null; | |
| 40 nonInterpolSub = false; | |
| 41 interpolCap = null; | |
| 42 nonInterpolCap = false; | |
| 29 } | 43 } |
| 30 | 44 |
| 31 public double getDiameterCap(int ndx) { | 45 public double getDiameterCap(int ndx) { |
| 32 if (diameterCap != null) { | 46 if (diameterCap != null) { |
| 33 return this.diameterCap.get(ndx); | 47 return this.diameterCap.get(ndx); |
| 47 return diameterCap.get(kms.indexOf(km)); | 61 return diameterCap.get(kms.indexOf(km)); |
| 48 } | 62 } |
| 49 return Double.NaN; | 63 return Double.NaN; |
| 50 } | 64 } |
| 51 | 65 |
| 66 public double getDiameterCapInterpol(double km) { | |
| 67 if (nonInterpolCap) { | |
| 68 return Double.NaN; | |
| 69 } | |
| 70 if (interpolCap == null) { | |
| 71 interpolCap = DoubleUtil.getLinearInterpolator(kms, diameterCap); | |
| 72 if (interpolCap == null) { | |
| 73 nonInterpolCap = true; | |
| 74 return Double.NaN; | |
| 75 } | |
| 76 } | |
| 77 try { | |
| 78 return interpolCap.value(km); | |
| 79 } catch (ArgumentOutsideDomainException e) { | |
| 80 /* This is expected for many results. */ | |
| 81 return Double.NaN; | |
| 82 } | |
| 83 } | |
| 84 | |
| 52 public double getDiameterSub(double km) { | 85 public double getDiameterSub(double km) { |
| 53 if (kms.indexOf(km) >= 0) { | 86 if (kms.indexOf(km) >= 0) { |
| 54 return diameterSub.get(kms.indexOf(km)); | 87 return diameterSub.get(kms.indexOf(km)); |
| 55 } | 88 } |
| 56 return Double.NaN; | 89 return Double.NaN; |
| 90 } | |
| 91 | |
| 92 public double getDiameterSubInterpol(double km) { | |
| 93 if (nonInterpolSub) { | |
| 94 return Double.NaN; | |
| 95 } | |
| 96 if (interpolSub == null) { | |
| 97 interpolSub = DoubleUtil.getLinearInterpolator(kms, diameterSub); | |
| 98 if (interpolSub == null) { | |
| 99 nonInterpolSub = true; | |
| 100 return Double.NaN; | |
| 101 } | |
| 102 } | |
| 103 try { | |
| 104 return interpolSub.value(km); | |
| 105 } catch (ArgumentOutsideDomainException e) { | |
| 106 /* This is expected for many results. */ | |
| 107 return Double.NaN; | |
| 108 } | |
| 57 } | 109 } |
| 58 | 110 |
| 59 public double[][] getDiameterCapData() { | 111 public double[][] getDiameterCapData() { |
| 60 return new double[][] { | 112 return new double[][] { |
| 61 kms.toNativeArray(), | 113 kms.toNativeArray(), |