python画蝴蝶结_Shapely中的多边形相交错误：“Shapely.geos.TopologicalError:geos intersection操作生成了空几何体”...

我一直在试着调试这个问题，但无法这样做。我试图找到两个Polygon对象的交集。它在大多数情况下都可以工作，但对于以下情况，它会引发以下异常：P1 area: 13.125721955

P2 area: 1.0

Traceback (most recent call last):

File "geom2d.py", line 235, in

print p1.intersection(p2)

File "/usr/local/lib/python2.7/dist-packages/shapely/geometry/base.py", line 334, in intersection

return geom_factory(self.impl['intersection'](self, other))

File "/usr/local/lib/python2.7/dist-packages/shapely/topology.py", line 47, in __call__

"The operation '%s' produced a null geometry. Likely cause is invalidity of the geometry %s" % (self.fn.__name__, repr(this)))

shapely.geos.TopologicalError: The operation 'GEOSIntersection_r' produced a null geometry. Likely cause is invalidity of the geometry

代码如下。from shapely.geometry import Point,Polygon,MultiPolygon

poly1 = [(35.0041000000000011, -88.1954999999999956), (34.9917999999999978, -85.6068000000000069), (32.8404000000000025, -85.1756000000000029), (32.2593000000000032, -84.8927000000000049), (32.1535000000000011, -85.0341999999999985), (31.7946999999999989, -85.1358000000000033), (31.5199999999999996, -85.0438000000000045), (31.3384000000000000, -85.0836000000000041), (31.2092999999999989, -85.1069999999999993), (31.0023000000000017, -84.9943999999999988), (30.9953000000000003, -87.6008999999999958), (30.9422999999999995, -87.5926000000000045), (30.8538999999999994, -87.6256000000000057), (30.6744999999999983, -87.4072000000000031), (30.4404000000000003, -87.3687999999999931), (30.1463000000000001, -87.5240000000000009), (30.1545999999999985, -88.3863999999999947), (31.8938999999999986, -88.4742999999999995), (34.8937999999999988, -88.1020999999999930), (34.9478999999999971, -88.1721000000000004), (34.9106999999999985, -88.1461000000000041)]

poly2 = [(34.7998910000000024, -88.2202139999999986), (34.7998910000000024, -87.2202139999999986), (35.7998910000000024, -87.2202139999999986), (35.7998910000000024, -88.2202139999999986)]

p1 = Polygon(poly1)

p2 = Polygon(poly2)

print 'P1 area:',p1.area

print 'P2 area:',p2.area

print p1.intersection(p2)

因为它打印了两个多边形的区域，所以我假设多边形的形状是正确的。我还(不知何故)打印了第一个多边形，以确保它确实是一个简单的多边形。

有人能解释一下为什么我会有这个例外吗？

编辑：我打印的p1.is_有效，结果是错误的。有一些解释here。搜索字符串is_valid。上面说A valid Polygon may not possess any overlapping exterior or interior rings.

有人能解释一下这意味着什么吗？如果有可能的话？

顺便说一下，我还注意到，如果我从poly1中删除最后一个坐标，那么整个过程都会正常工作。也许整个坐标列表使多边形变得复杂。