一、多边形边界顺序

多边形可以由一个点集 { v 1 , v 2 , . . . , v n } \{v_1,v_2,...,v_n\} {v1,v2,...,vn} 表示,构成多边形的点集确定,多边形边界的顺序也就确定了

有关多边形边界的顺序的示意图,如下图所示:

【计算几何】判断多边形边界顺逆时针 & C++代码实现-LMLPHP

有时候关于数据格式有规定,要求多边形边界必须为顺时针或者逆时针。

因此,我们需要对多边形边界的顺序进行判断,如果不符合要求,则需要将构成多边形的点集进行反转


二、数学原理

2.1 Green公式

Green公式揭示了平面区域的二重积分和封闭曲线上的线积分的关系:沿着多边形的边求曲线积分,若积分为正,则是沿着边界曲线正方向(逆时针),反之为顺时针,且所得绝对值为多边形面积

判断多边形是逆时针还是顺时针,就更简单了,直接计算多边形面积(不取绝对值):

  • 如果面积为正值,是逆时针
  • 如果面积为负值,是顺时针

计算的时候要注意,多边形的边界是个环,首尾闭合,最后一点和起点,坐标相同,别少个节点。

2.2 鞋带公式

鞋带公式也叫高斯面积公式,是一种数学算法,可求确定区域的一个简单多边形的面积。

所谓“简单多边形”,可以是凹、或凸多边形,但原则上边与边之间不能有交叉。

鞋带公式如下:

S = 0.5 ⋅ ∣ ∑ i = 1 n ( x i ⋅ y i + 1 − y i ⋅ x i + 1 ) ∣ S = 0.5 · |\sum^n_{i=1}{(x_i·y_{i+1}-y_i·x_{i+1})}| S=0.5i=1n(xiyi+1yixi+1)

因为判断环的方向只需要知道多边形面积的正负,所以最终使用的“简化版鞋带公式”为:

S ′ = ∑ i = 1 n ( x i ⋅ y i + 1 − y i ⋅ x i + 1 ) S' = \sum^n_{i=1}{(x_i·y_{i+1}-y_i·x_{i+1})} S=i=1n(xiyi+1yixi+1)


三、代码实现

// 点 对象
struct Vertex{
    Vertex(double px=0.0, double py=0.0) {
        x = px;
        y = py;
    }
    double x;
    double y;
};

// 多边形对象
class Poly
{
public:
	vector<Vertex> vertex; // 顶点集
};

// 判断多边形顺逆时针(true:顺时针,false:逆时针)
bool judgePolyDirection(Poly poly){
	double s = 0;
	for (int i = 0; i < poly.vertex.size(); i++){
		int j = i + 1 < poly.vertex.size() ? i +1:0;
		s += (poly.vertex[i].x * poly.vertex[j].y - poly.vertex[i].y * poly.vertex[j].x);
	}
	return s < 0;
}
05-30 23:20