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

最新推荐文章于 2025-01-03 17:32:22 发布
最新推荐文章于 2025-01-03 17:32:22 发布
阅读量3.6k 收藏 14
点赞数
CC 4.0 BY-SA版权
分类专栏： numpy OCR 文章标签： polygon面积
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/fanzonghao/article/details/103260817
本文介绍了使用Shapely库计算不规则多边形面积的方法，包括逆时针坐标系下的多边形调整为顺时针，以及通过分解梯形计算面积。提供了Python代码示例，适用于处理地理空间数据和计算机图形学中的面积计算问题。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

一．利用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

 

ArcObject java代码Polygon逆时针转顺时针
挣钱花388的博客
05-25 251
ArcObject java代码Polygon逆时针转顺时针
Polygon-Area
05-17
多边形区域 我提出了一个简单的程序来应对个人挑战，以找到任何2D多边形内的区域。 这是通过找到多边形的最高3个点，计算这三个点内的面积，然后删除多边形的最高点并重复进行的。 注意：此计算包括原始多边形内任何多边形的面积。 IE，如果此程序以其内部具有三角形的正方形作为输入参数运行，则正方形的面积和三角形的面积都将包括在计算中。 我认为这是一个错误，并希望将来进行修复。
判断一个多边形是顺时针还是逆时针的方法
weixin_34161083的博客
05-18 1416
1、关于如何判定多边形是顺时针还是逆时针对于凸多边形而言，只需对某一个点计算叉积 = ((xi - xi-1),(yi - yi-1)) x ((xi+1 - xi),(yi+1 - yi)) = (xi - xi-1) * (yi+1 - yi) - (yi - yi-1) * (xi+1 - xi) 如果上式的值为正，逆时针；为负则是顺时针。 而对于一般的简单多边形，则需对于多边形的每一个点计...
计算多边形（polygon）面积的算法原理和python实现
Enjoy Coding
11-12 8241
计算多边形（polygon）的面积其实很简单，假设是n边形，分割为n-2个三角形，积分别计算每个三角形面积后累加得多边形面积。 如下图所示： 那么如何计算三角形的面积呢？ △ABC的面积是“向量AB”和“向量AC”两个向量叉乘的绝对值的一半。其正负表示三角形顶点是在右手系还是左手系。 代码如下： import math import re VERTICES_FILE = "vertices.txt" class Point: def __init__(self, x, y):.
Area
当你有能力站在巅峰时
04-23 446
D - Area Crawling in process...Crawling failedTime Limit:1000MS    Memory Limit:10000KB    64bit IO Format:%I64d & %I64u SubmitStatus Description You are going to compute the area of a
polygon_area:计算具有有序顶点列表的 2D 多边形的面积。 顶点是实值的。-matlab开发
05-29
A =polygon_area(M) 计算由包含在矩阵 M 中的顶点（顺时针或逆时针）有序列表的多边形包围的区域，其中 x 坐标在第一列中给出，y 坐标在第二列中给出柱子。
如何判断一个多边形的环是逆时针还是顺时针
吴道简的博客
08-18 2011
OGC关于数据格式有规定，多边形的外环应为逆时针，内环应为顺时针。GIS相关的库包，在对数据进行空间操作的时候，也基于这个规定，将数据标准化为外环逆时针，内环顺时针之后，再进行空间计算和分析。 这种标准的数据结构，会将计算量减少一半。不过我现在看的大部分面数据，都不是严格执行这个规定的。数据不规范的话，数据量越大，处理时间越长，效率越低。所以我一般拿到数据的时候，都会把数据进行规范化标准化。 重点：数据处理中，数据标准化是前置步骤。 如何判断多边形是逆时针还是顺时针，比较简..
C#实现多边形方向与顶点凹凸性判断方法
在`Polygon`类中，可以实现一个方法`IsClockwise`来判断多边形的方向，通过计算所有相邻顶点向量的叉积并求和，如果和大于零，则多边形是顺时针方向的，否则为逆时针方向。对于凹凸性的判断，可以实现一个方法`...
计算几何学 | 点的内包 | Polygon-Point Containment | C/C++实现
qq_40998706的博客
02-17 1164
问题描述 对于多边形g与点p，p位于g内时输出2，p位于g的边上时输出1，其余情况输出0。 多边形g由其顶点的序列p1,p2,...,pnp_1,p_2,...,p_np1​,p2​,...,pn​表示，相邻两点pip_ipi​、pi+1p_{i+1}pi+1​（1 ≤ i ≤ n-1）相连构成g的边。另外，点pnp_npn​和p1p_1p1​相连的线也是g的边。 请注意，g不一定是凸多边形。 输...
计算两个多边形的重叠面积（C++）
qq_43561784的博客
04-19 1591
【代码】计算两个多边形的重叠面积（C++）
多边形Polygon
跟随技术的脚步-linolzhang的专栏
02-21 3608
多边形 大家都很熟悉，我们把他定义为 多条边首尾连接的封闭图形，关于多边形的算法有很多，比如 是否为凸多边形、直线与多边形求交、填充算法、三角化等，具体的算法原理作者就不再细讲，直接给出作者以前写的具体代码，大家作为参考：         封装的Math工具类（MathTool.h），这里也用到了之前文章里的（Tuple2.h和Tuple3.h） /* Math工具类 - 包含基本数学算法
多边形（Polygon）
最新发布
欢迎来到《技术探索》，这是一个专注于游戏开发技术的博客。在这里，我们将深入探讨游戏引擎、图形渲染、人工智能、物理模拟等领域的最新技术和最佳实践。无论您是初学者还是经验丰富的开发者，我们都希望为您提供有价值的见解和实用的技巧。
01-03 1794
多边形是由有限个线段（边）连接而成的闭合图形。每个线段的端点称为多边形的顶点。边：连接两个顶点的线段。顶点：多边形的角点。闭合性：多边形的边连接形成一个封闭的区域。计算多边形的面积。计算多边形的周长。判断多边形是凸多边形还是凹多边形。判断一个点是否在多边形内部。使用可视化多边形。这些功能可以为后续的图形处理、游戏开发或几何计算提供基础。你可以根据需要进一步扩展这个类，例如添加更多的几何操作、支持不同类型的多边形等。通过以上实现，我们可以有效地判断一个点是否在多边形内。
CGL_3_A——Area（多边形面积）
creator平
08-16 319
Area For a given polygong, computes the area of the polygon. gis represented by a sequence of pointsp1,p2,...,pnwhere line segments connectingpiandpi+1(1 ≤i≤n−1) are sides ofg. The l...
python计算多边形面积
Double_V的博客
02-18 6600
本文提供一个简单的方法计算多边形面积，参考维基百科 实现代码： def polygon_area(polygon): """ compute polygon area polygon: list with shape [n, 2], n is the number of polygon points """ area = 0 q = polygon[-1] for p in polygon: area += p[0] * q[1] -
js判断多边形的坐标点是顺时针还是逆时针的两种方法
qq_34447899的博客
06-28 3364
js判断多边形的坐标点是顺时针还是逆时针的两种方法关键算法通过极值点即其相邻点的走向算出多边形走向通过计算一边的角度和来推算多边形走向 由于arcgis对顺时针生成的多边形面积计算为正，逆时针生成的为负，所以产生了需要判断多边形坐标是顺时针还是逆时针的需求。 关键算法 首先我们可以通过三点构成的两个矢量的叉乘结果 k 来判断这三点的走向是顺时针还是逆时针的。 如果k<0，则为顺时针， 如果k...
计算多边形的面积
lh3325251325的专栏
11-02 981
  已知凸多边的各个顶点的坐标可以计算其面积，详细如下：   代码如下：// Polygon_area.cpp : 定义控制台应用程序的入口点。//#include "stdafx.h"#include#include#includeusing namespace std;struct point{ doubl
多边形顶点的顺逆时针判断
AlieZ的博客
07-20 4722
方法一：格林公式
Java、正多边形的面积
m0_62659797的博客
01-21 945
正多边形是一个具有n条边的多边形，它每条边的长度都相等，而且所有角的度数也相等（即多边形既等边又等角）。计算正多边形面积的公式是： 这里，s是边长。编写一个程序，提示用户输入边的个数以及正多边形的边长，然后显示它的面积。 package pack2; import java.util.Scanner; public class AreaOfPolygon { public static void main(String[] args) { try(Scanner inp...
判断多边形顺时针还是逆时针--根据多边形面积正负
u012138730的专栏
04-10 8174
最近在看recast&amp;detour源码的时候有遇到许多数学上的算法问题，特此记录，以便以后查看。问题：给出一个点序列，判断按这个点序列连接组成的多边形是逆时针还是顺时针。计算面积--鞋带公式：其中：    公式1【源码中所用的公式】           公式2 边的和【见参考中 overflow 链接】去掉上述计算S中的绝对值，让其有正负，判断逆顺：    点序为顺， 面积为负值。    ...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值