本篇文章小编给大家分享一下Python计算任意多边形间的重叠面积代码示例,文章代码介绍的很详细,小编觉得挺不错的,现在分享给大家供大家参考,有需要的小伙伴们可以来看看。
1. shapely工具箱
判断任意两个多边形的面积是否有交集有函数A.intersects(B)实现,若A和B有交集,返回为TRUE,若没有交集,返回为FALSE;
计算任意两个多边形的面积有函数A.intersection(B).area,输出直接为A和B的交集的面积。
那么问题在于,如何将多边形顶点的坐标,换为组成多边形的所有内部点的坐标。
python有个工具箱shapely。用于解决多边形有关问题。其中有个多边形填充函数Polygon.
shapely的安装方法为:
pip install Shapely
但是不知道为什么,用该方法安装时,一直报错,后来直接从https://www.*l*fd.*uci.edu/~gohlke/pythonlibs/#shapely (直接ctrl+F搜索shapely快速转到下载位置)下载了该文件,用地址名安装。
安装成功后,即可用下面示意的程序计算面积。
2. 程序
import numpy as np
import time
from shapely.geometry import Polygon # 多边形
import scipy.io as io
def Cal_area_2poly(data1,data2):
"""
任意两个图形的相交面积的计算
:param data1: 当前物体
:param data2: 待比较的物体
:return: 当前物体与待比较的物体的面积交集
"""
poly1 = Polygon(data1).convex_hull # Polygon:多边形对象
poly2 = Polygon(data2).convex_hull
if not poly1.intersects(poly2):
inter_area = 0 # 如果两四边形不相交
else:
inter_area = poly1.intersection(poly2).area # 相交面积
return inter_area
data1 = [] # 带比较的第一个物体的顶点坐标
data2 = [] #待比较的第二个物体的顶点坐标
area = Cal_area_2poly(data1,data2)