轉(zhuǎn)載自
對(duì)于凸多邊形，很容易計(jì)算，如下圖，以多邊形的某一點(diǎn)為頂點(diǎn)，將其劃分成幾個(gè)三角形，計(jì)算這些三角形的面積，然后加起來即可。已知三角形頂點(diǎn)坐標(biāo)，三角形面積可以利用向量的叉乘來計(jì)算。

對(duì)于凹多邊形，如果還是按照上述方法劃分成三角形，如下圖，多邊形的面積 = S_ABC + S_ACD + S_ADE, 這個(gè)面積明顯超過多邊形的面積。



該算法還可以優(yōu)化一下，對(duì)上面的式子合并一下同類項(xiàng)

S = S_OAB + S_OBC + S_OCD + S_ODE + S_OEA =
0.5*(A_y*(E_x-B_x) + B_y*(A_x-C_x) + C_y*(B_x-D_x) + D_y*(C_x-E_x) + E_y*(D_x-A_x))

struct Point2d {double x;double y;Point2d(double xx, double yy): x(xx), y(yy){} };//計(jì)算任意多邊形的面積，頂點(diǎn)按照順時(shí)針或者逆時(shí)針方向排列 double ComputePolygonArea(const vector<Point2d> &points) {int point_num = points.size();if(point_num < 3)return 0.0;double s = 0;for(int i = 0; i < point_num; ++i)s += points[i].x * points[(i+1)%point_num].y - points[i].y * points[(i+1)%point_num].x;return fabs(s/2.0); } struct Point2d {double x;double y;Point2d(double xx, double yy): x(xx), y(yy){} };//計(jì)算任意多邊形的面積，頂點(diǎn)按照順時(shí)針或者逆時(shí)針方向排列 double ComputePolygonArea(const vector<Point2d> &points) {int point_num = points.size();if(point_num < 3)return 0.0;double s = points[0].y * (points[point_num-1].x - points[1].x);for(int i = 1; i < point_num; ++i)s += points[i].y * (points[i-1].x - points[(i+1)%point_num].x);return fabs(s/2.0); }

總結(jié)

以上是生活随笔為你收集整理的计算任意多边形的面积的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。