欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

计算几何——多边形面积

程序员文章站 2022-04-01 17:55:48
...

规定顶点逆时针方向的多边形面积为正

计算几何——多边形面积

取原点0(0,0)

S(a,b,c,d)=S(o,a,b)+S(o,b,c)+S(o,c,d)+S(o,d,a)

图中o在外部,obc,ocd, oda全部是正的,加起来构成了5变形obcda的面积,但oab是负面积,因此最后刚好等于abcd的面积。

而如果o点在多边形内部,则四个三角形全是正面积,等式正好成立。

代码

#include<cmath>


class point{
 	public:
	double x;
	double y;
	point(double x_=0,double y_=0):x(x_),y(y_){} 
	friend const point operator+(const point& p1,const point& p2){
		return point(p1.x+p2.x,p1.y+p2.y);
	};
	friend const point operator-(const point& p1,const point& p2){
		return point(p1.x-p2.x,p1.y-p2.y);
	};
	friend const point operator*(const point& p,const double& m){
		return point(p.x*m,p.y*m);
	};
	friend const point operator*(const double& m,const point& p){
		return point(p.x*m,p.y*m);
	};
	friend const point operator/(const point& p,const double& m){
		return point(p.x/m,p.y/m);
	};
	friend ostream& operator <<(ostream& out,point& a){
		printf("(%lf,%lf)",a.x,a.y);
		return out;
	};
};
typedef point vect2;//重命名,向量也是用坐标表示 

class line{
	public:
	point start;
	point end; 
	line(point s=point(0,0),point e=point(0,0)):start(s),end(e){}
};
double cross(point O,point A,point B){//叉乘 
	double oa_x=A.x-O.x;
	double oa_y=A.y-O.y;
	double ob_x=B.x-O.x;
	double ob_y=B.y-O.y;
	return oa_x*ob_y-oa_y*ob_x;
}
double polyArea(point *ex,int num){//求多边形面积 
	if(num<3)return 0.0;
	double s=0;
	point O;O.x=0;O.y=0;//取原点为
	for(int i=0;i<num;i++)
		s+=cross(O,ex[i],ex[(i+1)%num]);
	return fabs(s/2.0);//取绝对值 
}