计算polygon面积和判断顺逆时针方向的方法

一．利用shapely求polygon面积

　　import shapely
    from shapely.geometry import Polygon, MultiPoint  # 多边形
    # box1 = [2, 0, 4, 2, 2, 4, 0, 2, 0, 0]
    box1 = [2, 0, 4, 2, 2, 4, 0, 2, 2, 2]
    poly_box1 = Polygon(np.array(box1).reshape(-1,2))
    print(poly_box1)
    print(poly_box1.area)

二．逆时针调整为顺时针

１．四个点 

def polygon_area1():
    '''
    compute area of a polygon
    :param poly:
    :return:
    '''
    # # #顺时针
    # poly = np.array([[2, 0],
    #                  [4, 2],
    #                  [2, 4],
    #                  [0, 2]])
    # 逆时针
    poly = np.array([[2, 0],
                     [0, 2],
                     [2, 4],
                     [4, 2]])
    poly_h = poly.shape[0]
    edge = []
    for i in range(poly_h - 1):
        edge.append((poly[i][0] - poly[i+1][0])*(poly[i][1] + poly[i + 1][1]))
    edge.append((poly[poly_h-1][0] - poly[0][0])*(poly[poly_h-1][1]+poly[0][1]))
    # edge = [
    #     (poly[1][0] - poly[0][0]) * (poly[1][1] + poly[0][1]),
    #     (poly[2][0] - poly[1][0]) * (poly[2][1] + poly[1][1]),
    #     (poly[3][0] - poly[2][0]) * (poly[3][1] + poly[2][1]),
    #     (poly[0][0] - poly[3][0]) * (poly[0][1] + poly[3][1])
    # ]
    area = np.sum(edge) / 2.
    print('=====The first way======')
    print('area:', area)
    if area <0:
        index = [0] + [i for i in range(poly_h - 1, 0, -1)]
        fix_poly = poly[index, :]
        print('fix_poly:', fix_poly)
def polygon_area2():
    '''
    compute area of a polygon
    :param poly:
    :return:
    '''
    # # #顺时针
    # poly = np.array([[2, 0],
    #                  [4, 2],
    #                  [2, 4],
    #                  [0, 2]])
    # 逆时针
    poly = np.array([[2, 0],
                     [0, 2],
                     [2, 4],
                     [4, 2]])
    edge = []
    poly_h = poly.shape[0]
    for i in range(poly_h-1):
        edge.append(poly[i][0] * poly[i+1][1] - poly[i+1][0] * poly[i][1])
    edge.append(poly[poly_h-1][0] * poly[0][1] - poly[0][0] * poly[poly_h-1][1])
    print('=====The second way======')
    # print('edge:', edge)
    area = np.sum(edge)/2.
    print('area:', area)
    #如果是逆时针　调整为顺时针
    if area < 0:
        index = [0] + [i for i in range(poly_h-1, 0, -1)]
        # print(index)
        fix_poly = poly[index, :]
        print('fix_poly：', fix_poly)
if __name__ == '__main__':
    # show_image()
    # compute_IOU_()
    polygon_area1()
    polygon_area2()

２．多个点

def polygon_area1():
    '''
    compute area of a polygon
    :param poly:
    :return:
    '''
    # # #顺时针
    # poly = np.array([[2, 0],
    #                  [4, 2],
    #                  [2, 4],
    #                  [0, 2]])
    # 逆时针 凹多边形
    print('凹多边形')
    poly = np.array([[2, 0],
                     [2, 2],
                     [0, 2],
                     [2, 4],
                     [4, 2]])
    # 逆时针 凸多边形
    # print('凸多边形')
    # poly = np.array([[2, 0],
    #                  [0, 0],
    #                  [0, 2],
    #                  [2, 4],
    #                  [4, 2]])
    poly_h = poly.shape[0]
    edge = []
    for i in range(poly_h - 1):
        edge.append((poly[i][0] - poly[i+1][0])*(poly[i][1] + poly[i + 1][1]))
    edge.append((poly[poly_h-1][0] - poly[0][0])*(poly[poly_h-1][1]+poly[0][1]))
    # edge = [
    #     (poly[1][0] - poly[0][0]) * (poly[1][1] + poly[0][1]),
    #     (poly[2][0] - poly[1][0]) * (poly[2][1] + poly[1][1]),
    #     (poly[3][0] - poly[2][0]) * (poly[3][1] + poly[2][1]),
    #     (poly[0][0] - poly[3][0]) * (poly[0][1] + poly[3][1])
    # ]
    area = np.sum(edge) / 2.
    print('=====The first way======')
    print('area:', area)
    if area <0:
        print('此为逆时针,调整为顺时针')
        index = [0] + [i for i in range(poly_h - 1, 0, -1)]
        fix_poly = poly[index, :]
        print('fix_poly:', fix_poly)
def polygon_area2():
    '''
    compute area of a polygon
    :param poly:
    :return:
    '''
    # # #顺时针
    # poly = np.array([[2, 0],
    #                  [4, 2],
    #                  [2, 4],
    #                  [0, 2]])
    # 逆时针 凹多边形
    print('凹多边形')
    poly = np.array([[2, 0],
                     [2, 2],
                     [0, 2],
                     [2, 4],
                     [4, 2]])
    # 逆时针 凸多边形
    # print('凸多边形')
    # poly = np.array([[2, 0],
    #                  [0, 0],
    #                  [0, 2],
    #                  [2, 4],
    #                  [4, 2]])
    edge = []
    poly_h = poly.shape[0]
    for i in range(poly_h-1):
        edge.append(poly[i][0] * poly[i+1][1] - poly[i+1][0] * poly[i][1])
    edge.append(poly[poly_h-1][0] * poly[0][1] - poly[0][0] * poly[poly_h-1][1])
    print('=====The second way======')
    # print('edge:', edge)
    area = np.sum(edge)/2.
    print('area:', area)
    #如果是逆时针　调整为顺时针
    if area < 0:
        print('此为逆时针,调整为顺时针')
        index = [0] + [i for i in range(poly_h-1, 0, -1)]
        # print(index)
        fix_poly = poly[index, :]
        print('fix_poly：', fix_poly)
if __name__ == '__main__':
    # show_image()
    # compute_IOU_()
    polygon_area1()
    polygon_area2()

　　　　　　　　　　　

三.计算不规则多边形面积

分解成梯形，面积依次相加

import numpy as np
points = np.array([[0.72, 2.28],
                   [2.66, 4.71],
                   [5,	 3.5],
                   [3.63, 2.52],
                   [4,	1.6],
                   [1.9,	1]])
print(points.shape)
h, w = points.shape
areas = []
for i in range(0, points.shape[0]):
    area = (points[(i+1) % h][-1] + points[i][-1])*(points[(i+1) % h][0] - points[i][0])/2
    areas.append(area)
    print(points[i])
    print(points[(i+1) % h])
print('==areas:', areas)
print(sum(areas))

 

参考：https://www.cnblogs.com/TenosDoIt/p/4047211.html