/* * Copyright (c) 2016 Vivid Solutions. * * All rights reserved. This program and the accompanying materials * are made available under the terms of the Eclipse Public License v1.0 * and Eclipse Distribution License v. 1.0 which accompanies this distribution. * The Eclipse Public License is available at http://www.eclipse.org/legal/epl-v10.html * and the Eclipse Distribution License is available at * * http://www.eclipse.org/org/documents/edl-v10.php. */ package org.locationtech.jts.geomgraph.index; /** * @version 1.7 */ import java.util.ArrayList; import java.util.Collections; import java.util.Iterator; import java.util.List; import org.locationtech.jts.geom.Coordinate; import org.locationtech.jts.geomgraph.Edge; /** * Finds all intersections in one or two sets of edges, * using a simple x-axis sweepline algorithm. * While still O(n^2) in the worst case, this algorithm * drastically improves the average-case time. * * @version 1.7 */ public class SimpleSweepLineIntersector extends EdgeSetIntersector { List events = new ArrayList(); // statistics information int nOverlaps; public SimpleSweepLineIntersector() { } public void computeIntersections(List edges, SegmentIntersector si, boolean testAllSegments) { if (testAllSegments) add(edges, null); else add(edges); computeIntersections(si); } public void computeIntersections(List edges0, List edges1, SegmentIntersector si) { add(edges0, edges0); add(edges1, edges1); computeIntersections(si); } private void add(List edges) { for (Iterator i = edges.iterator(); i.hasNext(); ) { Edge edge = (Edge) i.next(); // edge is its own group add(edge, edge); } } private void add(List edges, Object edgeSet) { for (Iterator i = edges.iterator(); i.hasNext(); ) { Edge edge = (Edge) i.next(); add(edge, edgeSet); } } private void add(Edge edge, Object edgeSet) { Coordinate[] pts = edge.getCoordinates(); for (int i = 0; i < pts.length - 1; i++) { SweepLineSegment ss = new SweepLineSegment(edge, i); SweepLineEvent insertEvent = new SweepLineEvent(edgeSet, ss.getMinX(), null); events.add(insertEvent); events.add(new SweepLineEvent(ss.getMaxX(), insertEvent)); } } /** * Because DELETE events have a link to their corresponding INSERT event, * it is possible to compute exactly the range of events which must be * compared to a given INSERT event object. */ private void prepareEvents() { Collections.sort(events); // set DELETE event indexes for (int i = 0; i < events.size(); i++ ) { SweepLineEvent ev = (SweepLineEvent) events.get(i); if (ev.isDelete()) { ev.getInsertEvent().setDeleteEventIndex(i); } } } private void computeIntersections(SegmentIntersector si) { nOverlaps = 0; prepareEvents(); for (int i = 0; i < events.size(); i++ ) { SweepLineEvent ev = (SweepLineEvent) events.get(i); if (ev.isInsert()) { processOverlaps(i, ev.getDeleteEventIndex(), ev, si); } } } private void processOverlaps(int start, int end, SweepLineEvent ev0, SegmentIntersector si) { SweepLineSegment ss0 = (SweepLineSegment) ev0.getObject(); /** * Since we might need to test for self-intersections, * include current INSERT event object in list of event objects to test. * Last index can be skipped, because it must be a Delete event. */ for (int i = start; i < end; i++ ) { SweepLineEvent ev1 = (SweepLineEvent) events.get(i); if (ev1.isInsert()) { SweepLineSegment ss1 = (SweepLineSegment) ev1.getObject(); // don't compare edges in same group, if labels are present if (! ev0.isSameLabel(ev1)) { ss0.computeIntersections(ss1, si); nOverlaps++; } } } } }