5

我需要摆脱形状中的自相交。形状由点数组构成，因此该形状的所有线段都是线。（只有直线，没有曲线和圆弧）

以前，我尝试从这些点创建 Path2D，从中构造一个 Area，然后使用它的 PathIterator 我创建了几个 Path2D，它们不知何故是先前路径的子路径，因此自干涉消失了。但这不适用于某些路径 - 自相交仍然存在。

那么你能指出我可以找到做类似事情的好算法的地方吗？

编辑：我在任何地方都没有发现任何有用的东西，所以我编写了自己的算法。查看答案。

4

3 回答 3

1

如果Area不适合您，您可以尝试使用GLUtessellator将您分解Shape为一组三角形，或者（使用该GL_LINE_LOOP选项）只是边界边缘。

于 2011-01-31T12:48:00.813 回答
1

所以，由于我在网上找不到类似的东西，我编写了自己的算法。

它可能非常无效，但它对我来说足够快。

它是这样的：

/**
 * Takes a polygon, defined by a list of lines, and splits it into several
 * paths on points of intersection. If non-self-intersected path is passed in,
 * the same path is returned.
 * @param path
 * @return
 */
public static List<List<Line2D>> splitPath(List<Line2D> lines) {
    List<List<Line2D>> splitted = new ArrayList<List<Line2D>>();
    // find intersections.
    loop1:
    for (Line2D l1 : lines) {
        for (Line2D l2 : lines) {
            if (l1 == l2) continue;
            Point2D intr;
            if ((intr = linesIntersect(l1, l2)) != null) {
                // creating two splitted subpaths
                int i1 = lines.indexOf(l1);
                int i2 = lines.indexOf(l2);

                List<Line2D> subpath1 = new ArrayList<Line2D>();
                subpath1.addAll(lines.subList(0, i1));
                subpath1.add(new Line2D.Double(l1.getP1(), intr));
                subpath1.add(new Line2D.Double(intr, l2.getP2()));
                subpath1.addAll(lines.subList(i2 + 1, lines.size()));
                splitted.addAll(splitPath(subpath1));

                List<Line2D> subpath2 = new ArrayList<Line2D>();
                subpath2.add(new Line2D.Double(intr, l1.getP2()));
                subpath2.addAll(lines.subList(i1 + 1, i2));
                subpath2.add(new Line2D.Double(l2.getP1(), intr));
                splitted.addAll(splitPath(subpath2));
                break loop1;
            }
        }
    }
    if (splitted.size() > 0) {
        return splitted;
    } else {
        return Collections.singletonList(lines);
    }
}

/**
 * Returns point of intersection of this line segments.
 * @param l1
 * @param l2
 * @return
 */
public static Point2D linesIntersect(Line2D l1, Line2D l2) {
    if (l1.getP1().equals(l2.getP2()) || l1.getP2().equals(l2.getP1())) return null;
    Point2D inter = lineIntersection(l1, l2);
    if (inter == null) return null;
    double infS = HEADER.infS;
    double x = inter.getX();
    if (((l1.getX1() > l1.getX2()) ? (x + infS > l1.getX2() && x - infS < l1.getX1()) : (x - infS < l1.getX2() && x + infS > l1.getX1())) &&
           ((l2.getX1() > l2.getX2()) ? (x + infS > l2.getX2() && x - infS < l2.getX1()) : (x - infS < l2.getX2() && x + infS > l2.getX1()))) {
        return inter;
    } else {
        return null;
    }
}

/**
 * Returns point of lines intersection, or null if they are parallel.
 * @param l1
 * @param l2
 * @return
 */
public static Point2D lineIntersection(Line2D l1, Line2D l2) {
    double a1 = l1.getY2() - l1.getY1();
    double b1 = l1.getX1() - l1.getX2();
    double c1 = a1*l1.getX1() + b1*l1.getY1();
    double a2 = l2.getY2() - l2.getY1();
    double b2 = l2.getX1() - l2.getX2();
    double c2 = a2*l2.getX1() + b2*l2.getY1();
    double det = a1*b2 - a2*b1;
    if (det != 0) {
        double x = (b2*c1 - b1*c2)/det;
        double y = (a1*c2 - a2*c1)/det;
        return new Point2D.Double(x, y);
    } else {
        // lines are parallel
        return null;
    }
}
于 2011-01-31T15:02:50.183 回答
1

我为你的问题/答案添加了书签，以防我自己必须实现类似的东西，但后来我发现了GEOS项目，它有一个简单的方法来实现这一点。我从 Python/Django 调用 GEOS，但整个事情都是基于JTS（Java 拓扑套件），所以我将从那里开始并将以下 Python 视为伪代码。

基本上，如果一条线不是简单连接的（此处解释），UNION 操作会将一条线拆分为简单连接的部分，因此使用它的第一点对线进行 UNION 可以满足我们的需要：

line  = LineString(list_of_lines_x_y_coordinates)
# union with first point splits into MultiLineString containing segments
segments = line.union(line[0])
于 2011-02-28T09:34:55.413 回答