＜基础数学＞ 多边形点集排序

多边形点集排序

一、问题描述

已知多边形点集P = {P1, P2, … , PN}，其排列顺序是杂乱，依次连接这N个点，无法形成确定的多边形，需要对点集P进行排序后，再绘制多边形。

二、数学原理

有两个二维向量$a(x_1,y_1)、b(x_2,y_2)$，两个向量的$a\times b=x_1\ast y_2-y_1\ast x_2$。$a\times b$可以判断旋转角度。如果$a\times b>0$时，那么$a$正旋转到$b$的角度为<180°，如果$a\times b<0$，那么$a$正旋转到$b$的角度为>180°，如果$a\times b=0$，那么a，b向量平行。

关于向量的叉乘右手定则判方向
$a\times b$的方向：四指由$a$开始，指向$b$，拇指的指向就是$a\times b$的方向，垂直于$a$和$b$所在的平面；
$b\times a$的方向：四指由$b$开始，指向$a$，拇指的指向就是$b\times a$的方向，垂直于$b$和$a$所在的平面；
$a\times b$的方向与$b\times a$的方向是相反的，且有：$a\times b=-b\times a$。

三、排序规则

点集排序过程中，关键在于如何定义点的大小关系。以按顺时针排序为例，算法步骤如下。

1. 计算点集的重心（质心）O，以重心作为旋转的中心点；
2. 根据向量叉积的定义，向量OA和OB的叉积大于0，即点B在点A逆时针方向。反之，点B在点A顺时针方向。

代码

//点集排序
//若点A大于点B,即点A在点B顺时针方向,返回true,否则返回false
bool PointCompare(const cv::Point &a, const cv::Point &b, const cv::Point &center)
{
    if (a.x >= 0 && b.x < 0)
        return true;
    if (a.x == 0 && b.x == 0)
        return a.y > b.y;
    //向量OA和向量OB的叉积
    int det = (a.x - center.x) * (b.y - center.y) - (b.x - center.x) * (a.y - center.y);
    if (det < 0)
        return true;
    if (det > 0)
        return false;
    //向量OA和向量OB共线，以距离判断大小
    int d1 = (a.x - center.x) * (a.x - center.x) + (a.y - center.y) * (a.y - center.y);
    int d2 = (b.x - center.x) * (b.x - center.y) + (b.y - center.y) * (b.y - center.y);
    return d1 > d2;
}
// 顺时针方向排序
void ClockwiseSortPoints(std::vector<cv::Point> &vPoints)
{
    //计算重心
    cv::Point center;
    int count_size = vPoints.size();
    double x = 0, y = 0;
    for (int i = 0; i < count_size; i++)
    {
        x += vPoints[i].x;
        y += vPoints[i].y;
    }
    center.x = (int)x / count_size;
    center.y = (int)y / count_size;
 
    //冒泡排序
    for (int i = 0; i < count_size - 1; i++)
    {
        for (int j = 0; j < count_size - i - 1; j++)
        {
            if (PointCompare(vPoints[j], vPoints[j + 1], center))
            {
                cv::Point tmp = vPoints[j];
                vPoints[j] = vPoints[j + 1];
                vPoints[j + 1] = tmp;
            }
        }
    }
 
    return;
}

参考文章

https://blog.csdn.net/beijingmake209/article/details/128563122
https://zhidao.baidu.com/question/528267710.html