Mercurial > dive4elements > river

comparison artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/tkhcalculation/BedQualityD50KmValueFinder.java @ 9081:a0feac6ea1ba

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed d50 query, added log.trace
author mschaefer
date Fri, 25 May 2018 09:10:19 +0200
parents 5294114b1df4
children b5dbcbaed42b
comparison
equal deleted inserted replaced
9080:89c1bb914f4e 9081:a0feac6ea1ba
55 * 55 *
56 * If PostgreSQL would support a median aggregate function like Oracle does, the aggregation could be placed into this 56 * If PostgreSQL would support a median aggregate function like Oracle does, the aggregation could be placed into this
57 * query. 57 * query.
58 */ 58 */
59 private static final String SQL_BED_D50_SUBLAYER_MEASUREMENT = // 59 private static final String SQL_BED_D50_SUBLAYER_MEASUREMENT = //
60 "SELECT t.km, t.datum, p.tiefevon, p.tiefebis, a.d50" // 60 "SELECT s.km, t.datum, p.tiefevon, p.tiefebis, a.d50" //
61 + " FROM sohltest t INNER JOIN station s ON t.stationid = s.stationid" // 61 + " FROM sohltest t INNER JOIN station s ON t.stationid = s.stationid" //
62 + " INNER JOIN gewaesser g ON s.gewaesserid = g.gewaesserid" // 62 + " INNER JOIN gewaesser g ON s.gewaesserid = g.gewaesserid" //
63 + " INNER JOIN sohlprobe p ON t.sohltestid = p.sohltestid" // 63 + " INNER JOIN sohlprobe p ON t.sohltestid = p.sohltestid" //
64 + " INNER JOIN siebanalyse a ON p.sohlprobeid = a.sohlprobeid" // 64 + " INNER JOIN siebanalyse a ON p.sohlprobeid = a.sohlprobeid" //
65 + " WHERE (g.name = :name) AND (s.km BETWEEN :fromkm - 0.0001 AND :tokm + 0.0001)" // 65 + " WHERE (g.name = :name) AND (s.km BETWEEN :fromkm - 0.0001 AND :tokm + 0.0001)" //
119 final TDoubleArrayList kms = new TDoubleArrayList(); 119 final TDoubleArrayList kms = new TDoubleArrayList();
120 final TDoubleArrayList values = new TDoubleArrayList(); 120 final TDoubleArrayList values = new TDoubleArrayList();
121 final TDoubleArrayList kmd50s = new TDoubleArrayList(); 121 final TDoubleArrayList kmd50s = new TDoubleArrayList();
122 122
123 // Median aggregate d50 values for each km 123 // Median aggregate d50 values for each km
124 for (int i = 0; i <= rows.size() - 1; i++) { 124 if (rows != null) {
125 kmd50s.add((double) rows.get(i)[4]); 125 for (int i = 0; i <= rows.size() - 1; i++) {
126 if (((i == rows.size() - 1) || !Utils.epsilonEquals((double) rows.get(i)[0], (double) rows.get(i + 1)[0], 0.0001))) { 126 log.trace("loadValues rows(" + i + ") " + rows.get(i)[0] + " " + rows.get(i)[1] + " " + rows.get(i)[2] + " " + rows.get(i)[3] + " "
127 final int k = kmd50s.size() / 2; 127 + rows.get(i)[4]);
128 values.add(((k + k < kmd50s.size()) ? kmd50s.get(k) : (kmd50s.get(k - 1) + kmd50s.get(k)) / 2) / 1000); 128 kmd50s.add((double) rows.get(i)[4]);
129 kms.add((double) rows.get(i)[0]); 129 if (((i == rows.size() - 1) || !Utils.epsilonEquals((double) rows.get(i)[0], (double) rows.get(i + 1)[0], 0.0001))) {
130 log.debug(String.format("loadValues km %.3f d50(mm) %.1f count %d", kms.get(kms.size() - 1), values.get(values.size() - 1) * 1000, 130 final int k = kmd50s.size() / 2;
131 kmd50s.size())); 131 values.add(((k + k < kmd50s.size()) ? kmd50s.get(k) : (kmd50s.get(k - 1) + kmd50s.get(k)) / 2) / 1000);
132 kmd50s.clear(); 132 kms.add((double) rows.get(i)[0]);
133 log.debug(String.format("loadValues km %.3f d50(mm) %.1f count %d", kms.get(kms.size() - 1), values.get(values.size() - 1) * 1000,
134 kmd50s.size()));
135 kmd50s.clear();
136 }
133 } 137 }
134 } 138 }
135 if (kms.size() >= 1) 139 if (kms.size() >= 1)
136 log.debug(String.format("loadValues: %d kms found from %.3f to %.3f", kms.size(), kms.get(0), kms.get(kms.size() - 1))); 140 log.debug(String.format("loadValues: %d kms found from %.3f to %.3f", kms.size(), kms.get(0), kms.get(kms.size() - 1)));
137 else 141 else
http://dive4elements.wald.intevation.org