计算几何：关于求多边形的重心

计算几何：关于求多边形的重心

 

将一个多边形分解为以编号0的点为原点的若干个三角形

求出每个三角形的面积a， //即两向量的叉乘的模,a = ( x1 * y2 - x2 * y1 ) / 2;

以及每个三角形的重心的坐标(x,y)； //即三个点的坐标的平均数，x = (x1+x2+x3)/3; y = (y1+y2+y3)/3;

smux = Sigma( xi * ai );

sumy = Sigma( yi * ai );

suma = Sigma( ai );

那末多边形的中心坐标为( sumx / suma , sumy / suma )。

 

计算时，可以将求x和y时算平均数的“/3”提出来，最后再除。

code;

		double sumx = 0,sumy = 0,suma = 0;
		scanf("%d",&N); //N用来计数点的个数
		scanf("%lf%lf%lf%lf",x,y,x+1,y+1);
		for ( ; N-2 > 0; N -- )
		{
			scanf("%lf%lf",x+2,y+2);
			a = ( ( x[1] - x[0] ) * ( y[2] - y[0] ) - ( x[2] - x[0] ) * ( y[1] - y[0] ) ) / 2 ;
			suma += a;
			sumx += ( x[0] + x[1] + x[2] ) * a;
			sumy += ( y[0] + y[1] + y[2] ) * a;
			x[1] = x[2];
			y[1] = y[2];
		}
		printf("%.2lf %.2lf\n",sumx / (suma*3),sumy / (suma*3) );