通过坐标点位，计算多边形面积

多边形面积

本文使用三角形面积累计法计算多边形面积，就是将多边形按照一个顶点，分割成多个三角形，计算三角形的面积，累加，得到多边形的面积。

当然，这个算法也有一些缺点，当这个多边形比较奇怪的时候，不能将多边形按照一个顶点分割成多个三角形面积之和，比如：

本文先按照常规的多边形进行处理，特殊情况之后再分析。

栗子

由于基本的测量点位坐标都是dat文本格式的，本文采用读取txt/dat格式的文本点位数据，更符合测量的操作。

data的数据如下所示：

这里简单以五个点组成的五边形为例，保存到path中。

def points(path):
    list_x = []
    list_y = []
    with open(path,'r',encoding='utf-8') as fp:
        for each in fp.readlines():
            pat = re.compile(',')
            x,y = float(pat.split(each.strip())[0]),float(pat.split((each.strip()))[1])
            list_x.append(x)
            list_y.append(y)
    return list_x,list_y

读取的数据分为两种，X、Y，为列表格式储存。

#三点组成的三角形的面积
def cal_triangle_area(a_x,a_y,b_x,b_y,c_x,c_y):
    triangle_area = 0.5*((b_x-a_x)*(c_y-a_y)-(b_y-a_y)*(c_x-a_x))
    return abs(triangle_area)

def poly(list_x,list_y):
    area = 0
    if len(list_x) < 3:
        return 0
    for i in range(len(list_x)-2):
        area += cal_triangle_area(list_x[0],list_y[0],list_x[i+1],list_y[i+1],list_x[i+2],list_y[i+2])
    return abs(area)

先写了一个三角形的面积计算函数，然后将多边形分成N-2个三角形，当N<3时，不能组成多边形，返回0，反之，计算按(X[0],Y[0])为顶点的分割的多个三角形的面积之和。

完整代码如下：

import re
import matplotlib.pyplot as plt

def points(path):
    list_x = []
    list_y = []
    with open(path,'r',encoding='utf-8') as fp:
        for each in fp.readlines():
            pat = re.compile(',')
            x,y = float(pat.split(each.strip())[0]),float(pat.split((each.strip()))[1])
            list_x.append(x)
            list_y.append(y)
    return list_x,list_y

#三点组成的三角形的面积
def cal_triangle_area(a_x,a_y,b_x,b_y,c_x,c_y):
    triangle_area = 0.5*((b_x-a_x)*(c_y-a_y)-(b_y-a_y)*(c_x-a_x))
    return abs(triangle_area)

def poly(list_x,list_y):
    area = 0
    if len(list_x) < 3:
        return 0
    for i in range(len(list_x)-2):
        area += cal_triangle_area(list_x[0],list_y[0],list_x[i+1],list_y[i+1],list_x[i+2],list_y[i+2])
    return abs(area)


if __name__ == '__main__':
    print('the area polyon is : {} '.format(poly(*points(path='./points.txt'))))
    plt.scatter(*points(path='./points.txt'),color='red')
    plt.show()

这里除了分割多边形需要优化之外，还有就是组成多边形的点位，这些点位在文件储存中，必须按照顺序储存，不然计算出的多边形可能不是我们想要的多边形的面积。可以试着改变一下点位顺序试试看。

这里还没想好怎么优化这个bug，对于分割的多边形的三角形，可以按照上面一个图例，寻找一些思路，就是判断分割的三角形内部是否有多边形的点位，如果三角形内部有顶点，则该部分的面积为分割的三角形的面积，减去内部的一个三角形的面积，这个方法可以自己动手试试看。

或者也可以先计算多边形的重心或中心等，或者内部的随机一个点位，将多边形以重心或中心或随机点位分割成一个个的三角形，计算按这个顶点分割的三角形的面积之和。