面积法判断多边形顺逆时针

面积法判断多边形顺逆时针

// 输入：按照顺序输入多边形的轮廓顶点
// 输出：0代表数据有问题，-1代表逆时针，1代表顺时针
int estimateContourType(std::vector<cv::Point2f> &contour) {
  int n = contour.size();
  float area = 0.0;
  if (n < 3) {
  // 多边形有问题
    return 0;
  }
  
  // 先将轮廓变成凸包，默认逆时针（通过索引方式）
  std::vector<int> index;
  cv::convexHull(contour, index);
  // 保证原来的顺序
  std::sort(index.begin(), index.end());
  n = index.size();
  for (int i = 0; i < index.size(); i++) {
    auto p0 = contour[index[i]];
    auto p1 = contour[index[(i + 1) % n]];
    area += (0.5 * (p1.y + p0.y) * (p1.x - p0.x));
  }

  cout << "contour type area=" << area << endl;
  if (area < 0.0f) {
    cout << "逆时针" << endl;
    return -1;
  }
  cout << "顺时针" << endl;
  return 1;
}