【C++】计算多边形的重心

可直接复制粘贴

#include <iostream>
#include <vector>
#include <math.h>

typedef struct point{
    double x = 0.;
    double y = 0.;
}Point;

/**
 *@brief 计算多边形的重心点
 *@param polygon 多边形的坐标列表
 *@return 重心坐标
*/
Point GetCenterOfGravity(std::vector<Point> polygon){
    if(polygon.size() <= 2)return Point();

    auto Area = [](Point p0,Point p1,Point p2){
        double area = 0.;
        area = p0.x * p1.y + p1.x * p2.y +
               p2.x * p0.y - p1.x * p0.y -
               p2.x * p1.y - p0.x * p2.y;
        return area / 2;
    };

    Point p0 = polygon[0];
    Point p1 = polygon[1];
    Point p2;
    double sumarea = 0, sumx = 0, sumy = 0;
    for (int i = 2; i < polygon.size(); i++)
    {
        p2 = polygon[i];
        double area = Area(p0, p1, p2);//求三角形的面积
        sumarea += area;
        sumx += (p0.x + p1.x + p2.x) * area; //求∑cx[i] * s[i]和∑cy[i] * s[i]
        sumy += (p0.y + p1.y + p2.y) * area;
        p1 = p2;//求总面积
    }

    Point barycenter;
    barycenter.x = sumx / sumarea / 3;
    barycenter.y = sumy / sumarea / 3;
    return barycenter;
}

int main()
{
    return 0;
}