网上可以搜索到相关的论文和代码,但是实际测试发现代码不完整,在多边形比较复杂或者特殊的时候所得的结果不正确,发现原来是代码缺少了一条边导致的,只需要加上最后一条边即可。
以下是修改之后的完整代码:
bool IsPolyClockwise(std::vector<point2d>& vPts)
{
//沿着多边形的边求曲线积分,若积分为正,则是沿着边界曲线正方向(逆时针),反之为顺时针
double d = 0;
const size_t nSize = vPts.size();
for (int i = 0; i < nSize - 1; ++i)
{
d += -0.5 * (vPts[i + 1].y + vPts[i].y)*(vPts[i + 1].x - vPts[i].x);
}
//这条边不能忘记
d += -0.5 * (vPts[0].y + vPts[nSize - 1].y)*(vPts[0].x - vPts[nSize - 1].x);
//小于零为顺时针,大于零为逆时针
return d < 0.0;;
}
感谢原作者的分享和参考。