如何使用python在2D中计算多边形IoU

2021-07-21 10:07:07 浏览数 (1)

本文W3Cschool和大家介绍一下如何使用python在2D中计算多边形IoU。假设多边形不是自相交的，即圆圆周围的点的顺序是单调的，那么我相信有一个相对简单的方法来确定 IoU 值，而不需要一个一般的形状包。

  1. 假设每个多边形的点是顺时针在圆圈中排列的。如果我们发现签名区域为负值，我们可以通过增加角度 w.r.t x 轴或倒车点来确保这一点。
  2. 将两个多边形的点合并到单个列表中，跟踪每个点属于哪个多边形。我们还需要能够确定每个点在原始多边形中的上一点和下一点。L
  3. 通过增加角度来排序w.r.t x轴。L
  4. 如果输入多边形相交，则从一个多边形到另一个多边形的过渡次数将大于两个。L
  5. 遍历。如果连续点属于不同的多边形，则第一点与其下一点和第二点之间的线的交叉点及其前一点将属于两个多边形之间的交点。L
  6. 将步骤 4 中确定的每个点添加到新的多边形中。将按顺序遇到积分。II
  7. 每个多边形的面积之和将等于其联盟加上交叉口的区域，因为这将计算两次。
  8. 因此，将的价值由两个多边形的面积之和减去面积之和来表示。IoUII

唯一需要的几何形状是使用Shoelace 公式计算简单多边形的面积，并确定步骤 5 所需的两条线段之间的交汇点。

这里有一些Java代码（Ideone）来说明 - 你也许可以使它在Python更紧凑。

double[][] coords = {{-0.708, 0.707, 0.309, -0.951, 0.587, -0.809},
                       {1, 0, 0, 1, -1, 0, 0, -1, 0.708, -0.708}};

double areaSum = 0;
List<CPoint> pts = new ArrayList<>();
for(int p=0; p<coords.length; p++)
{
    List<CPoint> poly = new ArrayList<>();
    for(int j=0; j<coords[p].length; j+=2)
    {
        poly.add(new CPoint(p, coords[p][j], coords[p][j+1]));
    }
    
    double area = area(poly);
    if(area < 0)
    {
        area = -area;
        Collections.reverse(poly);
    }
    areaSum += area;
    
    pts.addAll(poly);

    int n = poly.size();
    for(int i=0, j=n-1; i<n; j=i++)
    {
        poly.get(i).prev = poly.get(j);
        poly.get(j).next = poly.get(i);             
    }
}       
        
pts.sort((a, b) -> Double.compare(a.theta, b.theta));
        
List<Point2D> intersections = new ArrayList<>();
int n = pts.size();
for(int i=0, j=n-1; i<n; j=i++)
{
    if(pts.get(j).id != pts.get(i).id)
    {
        intersections.add(intersect(pts.get(j), pts.get(j).next, pts.get(i).prev, pts.get(i)));
    }
}

double areaInt = area(intersections);
double iou = areaInt/(areaSum - areaInt);
System.out.println(iou);

输出：

0.12403616470027268

和支持代码：

static class CPoint extends Point2D.Double
{
    int id;
    double theta;
    CPoint prev, next;
    
    public CPoint(int id, double x, double y)
    {
        super(x, y);
        this.id = id;
        theta = Math.atan2(y, x);
        if(theta < 0) theta = 2*Math.PI + theta;
    }
}   

static double area(List<? extends Point2D> poly)
{
    double area = 0;
    for(int i=0, j=poly.size()-1; i<poly.size(); j=i++)
        area += (poly.get(j).getX() * poly.get(i).getY()) - (poly.get(i).getX() * poly.get(j).getY());
    return Math.abs(area)/2;
}

// https://rosettacode.org/wiki/Find_the_intersection_of_two_lines#Java
static Point2D intersect(Point2D p1, Point2D p2, Point2D p3, Point2D p4)
{
  double a1 = p2.getY() - p1.getY();
  double b1 = p1.getX() - p2.getX();
  double c1 = a1 * p1.getX() + b1 * p1.getY();

  double a2 = p4.getY() - p3.getY();
  double b2 = p3.getX() - p4.getX();
  double c2 = a2 * p3.getX() + b2 * p3.getY();
Python

0 人点赞