Mercurial > dive4elements > gnv-client
view gnv-artifacts/src/main/java/de/intevation/gnv/raster/IsoPolygonSeriesProducer.java @ 453:537e663d6c0c
vertical cross section: Fixed problem with running to max depth on every column in interpolation.
gnv-artifacts/trunk@501 c6561f87-3c4e-4783-a992-168aeb5c3f6f
author | Sascha L. Teichmann <sascha.teichmann@intevation.de> |
---|---|
date | Mon, 04 Jan 2010 13:42:49 +0000 |
parents | 3cb2bea50456 |
children | f7038820df2e |
line wrap: on
line source
package de.intevation.gnv.raster; import java.util.List; import java.util.HashSet; import java.util.HashMap; import java.util.ArrayList; import java.util.Collection; import org.apache.log4j.Logger; import gnu.trove.TIntHashSet; import gnu.trove.TIntObjectHashMap; import gnu.trove.TDoubleArrayList; import gnu.trove.TObjectProcedure; import de.intevation.gnv.raster.Vectorizer.RingsHandler; import de.intevation.gnv.raster.Vectorizer.Edge; import de.intevation.gnv.math.IJKey; import de.intevation.gnv.jfreechart.PolygonSeries; import de.intevation.gnv.jfreechart.CompactXYItems; /** * @author Sascha L. Teichmann (sascha.teichmann@intevation.de) */ public class IsoPolygonSeriesProducer implements RingsHandler { private static Logger log = Logger.getLogger( IsoPolygonSeriesProducer.class); public static final Float LINE_WIDTH = Float.valueOf(0.1f); public interface AttributeGenerator { Object generateAttribute(int neighbor1, int neighbor2); } // interface AttributeGenerator protected HashMap<Edge, Integer> open; protected HashMap<IJKey, TIntObjectHashMap> commonOpen; protected HashMap<IJKey, ArrayList<Edge>> complete; protected int width; protected int height; protected double minX; protected double minY; protected double maxX; protected double maxY; public IsoPolygonSeriesProducer( double minX, double minY, double maxX, double maxY ) { this.minX = minX; this.minY = minY; this.maxX = maxX; this.maxY = maxY; open = new HashMap<Edge, Integer>(); commonOpen = new HashMap<IJKey, TIntObjectHashMap>(); complete = new HashMap<IJKey, ArrayList<Edge>>(); } public void handleRings( List<Edge> rings, int value, int width, int height ) { if (value == -1) { return; } this.width = width; this.height = height; Integer v = Integer.valueOf(value); for (Edge head: rings) { Edge current = head; do { Integer neighbor = open.remove(current); if (neighbor != null) { IJKey pair = new IJKey(value, neighbor.intValue()); pair.sort(); TIntObjectHashMap co = commonOpen.get(pair); if (co == null) { commonOpen.put(pair, co = new TIntObjectHashMap()); } Edge edge = new Edge(current); Edge otherA = (Edge)co.remove(edge.a); if (otherA != null) { otherA.chain(edge, edge.a); } Edge otherB = (Edge)co.remove(edge.b); if (otherB != null) { otherB.chain(edge, edge.b); } if (edge.isComplete()) { ArrayList list = complete.get(pair); if (list == null) { complete.put(pair, list = new ArrayList()); } list.add(Vectorizer.simplify(edge, width)); } else { if (otherA == null) { co.put(edge.a, edge); } if (otherB == null) { co.put(edge.b, edge); } } } else { Edge edge = new Edge(current.b, current.a); open.put(edge, v); } current = current.next; } while (current != head); } // for all rings } // handleRings public Collection<PolygonSeries> getSeries() { return getSeries(null); } public Collection<PolygonSeries> getSeries( AttributeGenerator attributeGenerator ) { ArrayList<PolygonSeries> series = new ArrayList<PolygonSeries>(); double b1 = minX; double m1 = width != 1 ? (maxX - minX)/(width-1) : 0d; double b2 = minY; double m2 = height != 1 ? (maxY - minY)/(height-1) : 0d; // join keys of complete and open polygons HashSet<IJKey> pairs = new HashSet<IJKey>(); for (IJKey key: complete.keySet()) { pairs.add(key); } for (IJKey key: commonOpen.keySet()) { pairs.add(key); } TDoubleArrayList vertices = new TDoubleArrayList(); for (IJKey key: pairs) { PolygonSeries ps = new PolygonSeries(); // process complete ArrayList<Edge> completeList = complete.get(key); if (completeList != null) { for (Edge head: completeList) { Edge current = head; do { vertices.add(m1*(current.a % width) + b1); vertices.add(m2*(current.a / width) + b2); } while ((current = current.next) != head); // add head again to close shape vertices.add(m1*(head.a % width) + b1); vertices.add(m2*(head.a / width) + b2); ps.addRing(new CompactXYItems(vertices.toNativeArray())); vertices.clear(); } } // process open TIntObjectHashMap map = commonOpen.get(key); if (map != null) { final ArrayList<Edge> headList = new ArrayList<Edge>(); map.forEachValue(new TObjectProcedure() { TIntHashSet headSet = new TIntHashSet(); public boolean execute(Object value) { Edge head = ((Edge)value).head(); if (headSet.add(head.a)) { headList.add(head); } return true; } }); for (Edge head: headList) { head = Vectorizer.simplify(head, width); Edge current = head, last = head; do { vertices.add(m1*(current.a % width) + b1); vertices.add(m2*(current.a / width) + b2); last = current; } while ((current = current.next) != null); // add b from tail vertices.add(m1*(last.b % width) + b1); vertices.add(m2*(last.b / width) + b2); ps.addRing(new CompactXYItems(vertices.toNativeArray())); vertices.clear(); } // for all in common open } // if map defined for key if (ps.getItemCount() > 0) { series.add(ps); if (attributeGenerator != null) { Object attribute = attributeGenerator .generateAttribute(key.i, key.j); if (attribute != null) { ps.setAttribute("label", attribute); } } ps.setAttribute("line.width", LINE_WIDTH); } } // for all pairs return series; } public void clear() { open.clear(); complete.clear(); commonOpen.clear(); } } // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8: