python计算图形面积的方法_Python计算任意多边形面积算法

本文介绍了一种计算多边形面积的方法,通过连接多边形各顶点与一个虚拟点P,计算三角形面积的和。代码实现中,先判断三角形是顺时针还是逆时针,然后计算并累加面积。最后,给出一个具体的四边形面积计算示例。
摘要由CSDN通过智能技术生成

多边形面积求解的方法很多,其中比较多见的就是增加一个点P,然后分别连接多边形各个定点与P点,然后计算每个三角形的符号面积(面积有正负之分),求和就可以计算出面积。

0818b9ca8b590ca3270a3433284dd417.png

鉴于上面的方法,我们也可以用另外一种更简单的方法,将P点直接放在V0点,那么我们要求的面积就是V0V1V2,VoV2V3,...,V0Vn-2Vn-3的面积和。此处注意,如果我们假设逆时针的三角形(如上图V0V2V3)为正,那么顺时针的面积(如上图V0V1V2)就为负。

代码如下:

#------------------------------------------------------------------------------- # Name:        多边形面积计算 # Purpose: # # Author:      Administrator # # Created:     20-02-2013 # Copyright:   (c) Administrator 2013 # Licence:     #------------------------------------------------------------------------------- import math class Point():     def __init__(self,x,y):         self.x = x         self.y = y def GetAreaOfPolyGon(points):     '''计算多边形面积值        points:多边形的点集,每个点为Point类型        返回:多边形面积'''     area = 0     if(len(points)<3):         raise Exception("至少需要3个点才有面积")     p1 = points[0]     for i in range(1,len(points)-1):         p2 = points[1]         p3 = points[2]         #计算向量         vecp1p2 = Point(p2.x - p1.x,p2.y - p1.y)         vecp2p3 = Point(p3.x - p2.x,p3.y - p2.y)         #判断顺时针还是逆时针,顺时针面积为正,逆时针面积为负         vecMult = vecp1p2.x*vecp2p3.y - vecp1p2.y*vecp2p3.x         sign = 0         if(vecMult>0):             sign = 1         elif(vecMult < 0):             sign = -1         triArea = GetAreaOfTriangle(p1,p2,p3)*sign         area+=triArea     return abs(area) def GetAreaOfTriangle(p1,p2,p3):     '''计算三角形面积'''     area = 0     p1p2 = GetLineLength(p1,p2)     p2p3 = GetLineLength(p2,p3)     p3p1 = GetLineLength(p3,p1)     s = (p1p2 + p2p3 + p3p1)/2     area = s*(s-p1p2)*(s-p2p3)*(s-p3p1)     area = math.sqrt(area)     return area def GetLineLength(p1,p2):     '''计算边长'''     length = math.pow((p1.x-p2.x),2) + math.pow((p1.y-p2.y),2)     length = math.sqrt(length)     return length def main():     p1 = Point(1,1)     p2 = Point(2,1)     p3 = Point(2,2)     p4 = Point(1,2)     points = [p1,p2,p3,p4]     area = GetAreaOfPolyGon(points)     print(math.ceil(area))     assert math.ceil(area)==1 if __name__ == '__main__':     main()

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值