Python - 检查 Shapely Polygon 是否为矩形

Question

在Polygon上进行分割等操作后，我想验证它是一个矩形。
我试过simplify，然后计算coords的数量是否为5...

>>> from shapely.geometry import Polygon
>>> from shapely.ops import split
>>> 
>>> poly1 = Polygon([(0, 0), (0, 1), (0, 3), (2, 3), (2, 2), (2, 0), (0, 0)])
>>> 
>>> poly_check=poly1.simplify(0)
>>> if len(poly_check.exterior.coords)==5:
>>>     print 'Yes, it is a rectangle...'
>>> else:
>>>     print 'No, it is not a rectangle...'
>>> 
Yes, it is a rectangle...

但是如果起点在边的中间就不行了。

>>> #poly2 is actually a rectangle
>>> poly2 = Polygon([(0, 1), (0, 3), (2, 3), (2, 2), (2, 0), (0, 0), (0, 1)])
>>> 
>>> poly_check=poly2.simplify(0)
>>> if len(poly_check.exterior.coords)==5:
>>>     print 'Yes, it is a rectangle...'
>>> else:
>>>     print 'No, it is not a rectangle...'
>>> 
No, it is not a rectangle...

如何检查？

谢谢

Answer 1

如果多边形的面积与其最小有界矩形的面积相匹配，则多边形是一个矩形。这是一个已知的形状指标，称为矩形度。

if poly.area == poly.minimum_rotated_rectangle.area:
    return True

编辑：考虑到下面关于浮点错误的评论，您可以直接测量矩形并将 > .99 的所有内容视为矩形或进行近似比较（例如舍入）。

if (poly.area / poly.minimum_rotated_rectangle.area) > .99:
    return True

编辑 2：最好只使用 math.isclose 函数来确定两个变量的相等性。前面比较中的除法降低了比较的整体精度，我们在这里避免：

import math

if math.isclose(poly.minimum_rotated_rectangle.area, poly.area):
    return True