Mercurial > dive4elements > gnv-client
view gnv-artifacts/src/main/java/de/intevation/gnv/state/profile/horizontal/HorizontalProfileMeshCrossOutputState.java @ 425:15b8e95fa8da
Horizontale Schnittprofile. Ignore invalid interpolation correctly.
gnv-artifacts/trunk@473 c6561f87-3c4e-4783-a992-168aeb5c3f6f
author | Sascha L. Teichmann <sascha.teichmann@intevation.de> |
---|---|
date | Mon, 21 Dec 2009 19:58:17 +0000 |
parents | 2402173a1490 |
children | bed9735adf84 |
line wrap: on
line source
/** * */ package de.intevation.gnv.state.profile.horizontal; import java.util.Arrays; import java.util.ArrayList; import java.util.Collection; import java.util.List; import java.util.Locale; import org.apache.log4j.Logger; import org.w3c.dom.Node; import com.vividsolutions.jts.io.ParseException; import de.intevation.artifactdatabase.Config; import de.intevation.gnv.artifacts.cache.CacheFactory; import de.intevation.gnv.geobackend.base.DefaultResult; import de.intevation.gnv.geobackend.base.DefaultResultDescriptor; import de.intevation.gnv.geobackend.base.Result; import de.intevation.gnv.geobackend.base.ResultDescriptor; import de.intevation.gnv.geobackend.base.query.QueryExecutor; import de.intevation.gnv.geobackend.base.query.QueryExecutorFactory; import de.intevation.gnv.geobackend.base.query.exception.QueryException; import de.intevation.gnv.math.Interpolation2D; import de.intevation.gnv.math.LinearMetrics; import de.intevation.gnv.math.Point2d; import de.intevation.gnv.chart.Chart; import de.intevation.gnv.chart.ChartLabels; import de.intevation.gnv.chart.HorizontalCrossProfileChart; import de.intevation.gnv.utils.DistanceCalculator; import de.intevation.gnv.utils.WKTUtils; import org.jfree.chart.ChartTheme; import com.vividsolutions.jts.geom.Coordinate; import com.vividsolutions.jts.geom.LineString; import com.vividsolutions.jts.io.WKTReader; /** * @author Tim Englich <tim.englich@intevation.de> * */ public class HorizontalProfileMeshCrossOutputState extends HorizontalProfileOutputState { /** * */ private static final long serialVersionUID = 2205958041745637263L; /** * the logger, used to log exceptions and additonaly information */ private static Logger log = Logger.getLogger(HorizontalProfileMeshCrossOutputState.class); private String ijkQueryID = null; /** * Constructor */ public HorizontalProfileMeshCrossOutputState() { super(); } /** * @see de.intevation.gnv.state.timeseries.TimeSeriesOutputState#setup(org.w3c.dom.Node) */ @Override public void setup(Node configuration) { super.setup(configuration); this.ijkQueryID = Config.getStringXPath(configuration,"queryID-ijk"); } @Override protected Chart getChart( ChartLabels chartLables, ChartTheme theme, Collection parameters, Collection measurements, Collection dates, Collection result, Locale locale, String uuid, boolean linesVisible, boolean shapesVisible ) { Chart chart = null; if (CACHE_CHART) { log.info("Try to get horizontalprofilemeshcross chart from cache."); chart = (Chart) getChartFromCache(uuid); } if (chart != null) return chart; log.info("Chart not in cache yet."); chart = new HorizontalCrossProfileChart( chartLables, theme, parameters, measurements, dates, result, null, locale, linesVisible, shapesVisible ); chart.generateChart(); if (CACHE_CHART) { log.info("Put chart into cache."); purifyChart(chart, uuid); } return chart; } @Override protected Collection<Result> getChartResult(String uuid) { log.debug("HorizontalProfileMeshCrossOutputState.getChartResult"); Collection<Result> result = null; if (CacheFactory.getInstance().isInitialized()) { String key = uuid + super.getID(); log.debug("Hash for Queryelements: " + key); net.sf.ehcache.Element value = CacheFactory.getInstance().getCache().get(key); if (value != null) { result = (Collection<Result>) (value.getObjectValue()); }else{ if (this.inputData.containsKey("mesh_linestring")){ try { LineString ls = (LineString)new WKTReader().read( inputData.get("mesh_linestring").getValue()); Coordinate[] coords = ls.getCoordinates(); String additionWhere = WKTUtils.worldCoordinatesToIndex( result, inputData, ijkQueryID ); String[] filterValues = generateFilterValuesFromInputData(); String[] addedFilterValues = new String[filterValues.length+1]; System.arraycopy( filterValues, 0, addedFilterValues, 0, filterValues.length); addedFilterValues[filterValues.length] = additionWhere; QueryExecutor queryExecutor = QueryExecutorFactory .getInstance() .getQueryExecutor(); result = process( Arrays.asList(coords), queryExecutor.executeQuery( queryID, addedFilterValues) ); } catch (ParseException e) { log.error(e,e); }catch (QueryException e) { log.error(e,e); } }else{ // TODO: definieren was passiert wenn kein linestring vorhanden ist. } if (CacheFactory.getInstance().isInitialized()) { CacheFactory.getInstance().getCache().put(new net.sf.ehcache.Element(key, result)); } } } return result; } @Override protected String createChartSubtitle(Locale locale, String uuid) { log.debug("create chart subtitle for horizontal crossprofile charts."); String subtitle = createTimePeriod(locale, uuid); return subtitle; } public static Collection<Result> process( List<Coordinate> path, Collection<Result> input ) { log.debug("------ number of points before processing: " + input.size()); ArrayList<Result> output = new ArrayList<Result>(); Result last = null; int [] diffColums = null; SectionHandler sectionHandler = null; for (Result result: input) { if (sectionHandler == null) { ResultDescriptor rd = result.getResultDescriptor(); diffColums = rd.getColumnIndices(DIFF_COLUMS); int columns = rd.getColumnCount(); DefaultResultDescriptor resultDescriptor = new DefaultResultDescriptor(); for (int j = 0; j < columns; ++j) { String columnName = rd.getColumnName(j); if (!blacklisted(columnName)) { resultDescriptor.addColumn( columnName, rd.getColumnClassName(j)); } } sectionHandler = new SectionHandler( path, output, resultDescriptor); sectionHandler.setPrototyp(result); } if (last != null && WKTUtils.different(last, result, diffColums)) { sectionHandler.finish(); sectionHandler.setPrototyp(result); } sectionHandler.handle(result); last = result; } if (sectionHandler != null) { sectionHandler.finish(); } log.debug("------ number of points after processing: " + output.size()); return output; } private static final boolean blacklisted(String column) { for (int i = 0; i < COLUMN_BLACKLIST.length; ++i) { if (COLUMN_BLACKLIST.equals(column)) { return true; } } return false; } private static final String [] DIFF_COLUMS = { "GROUP1", "GROUP2", "GROUP3" }; private static final String [] COLUMN_BLACKLIST = { "MEDIAN.MESHPOINT.JPOSITION", "MEDIAN.MESHPOINT.IPOSITION" }; public static final double EPSILON = 1e-5d; public static final int INTERPOLATION_STEPS = Integer.getInteger("interpolation.steps", 500).intValue(); public static final class SectionHandler implements Interpolation2D.Consumer { private ArrayList<Point2d> points; private List<Coordinate> path; private Collection<Result> output; private Result prototyp; private ResultDescriptor descriptor; private boolean lastWasSuccess; public SectionHandler() { } public SectionHandler( List<Coordinate> path, Collection<Result> output, ResultDescriptor descriptor ) { this.path = path; this.output = output; this.descriptor = descriptor; points = new ArrayList<Point2d>(); lastWasSuccess = true; } public void finish() { if (!points.isEmpty()) { double distance = WKTUtils.toKM( DistanceCalculator.calculateDistance(path)); if (distance > EPSILON) { Interpolation2D.interpolate( path, points, 0d, distance, INTERPOLATION_STEPS, LinearMetrics.INSTANCE, this); } points.clear(); } lastWasSuccess = true; } public void setPrototyp(Result prototyp) { this.prototyp = prototyp; } public void handle(Result result) { Coordinate coordinate = WKTUtils.toCoordinate(result.getString("SHAPE")); double value = result.getDouble("YORDINATE"); int iPos = result.getInteger("MEDIAN.MESHPOINT.JPOSITION"); int jPos = result.getInteger("MEDIAN.MESHPOINT.JPOSITION"); Point2d p = new Point2d( coordinate.x, coordinate.y, value, iPos, jPos); points.add(p); } public void interpolated(Coordinate coordinate, boolean success) { if (!success && !lastWasSuccess) { // only insert null if last was valid. // This prevents flooding the result set with nulls // if interpolating over a large gap. return; } DefaultResult result = new DefaultResult(descriptor); ResultDescriptor pd = prototyp.getResultDescriptor(); int pcolums = pd.getColumnCount(); for (int i = 0, j = 0; i < pcolums; ++i) { String colname = pd.getColumnName(i); if (blacklisted(colname)) { continue; } if (colname.equals("SHAPE")) { result.addColumnValue(j, WKTUtils.toWKT(coordinate)); } else if (colname.equals("YORDINATE")) { result.addColumnValue(j, success ? Double.valueOf(coordinate.z) : null); } else { result.addColumnValue(j, prototyp.getObject(i)); } ++j; } output.add(result); lastWasSuccess = success; } } }