通过减去与另一个多边形的交集来创建新的形状多边形

Question

我有两个形状多边形实例（由lon,lat点组成），在不同的部分相交。我试图遍历，确定两个多边形之间是否有交集，然后创建一个新的多边形，排除该交集。从所附的图像来看，我基本上不希望红色的圆圈与黄色的轮廓重叠，我希望边缘正好是黄色轮廓开始的地方。

我已经尝试按照这里的说明，但它根本没有改变我的输出，而且我不想将它们合并到一个级联的联合中。我没有得到任何错误消息，但是当我将这些多边形添加到一个KML文件中（只是在python中进行原始文本操作，没有任何花哨的程序）时，它们仍然显示为圆形，没有任何修改。

# multipol1 and multipol2 are my shapely MultiPolygons
from shapely.ops import cascaded_union
from itertools import combinations
from shapely.geometry import Polygon,MultiPolygon

outmulti = []
for pol in multipoly1:
    for pol2 in multipoly2:
        if pol.intersects(pol2)==True:
            # If they intersect, create a new polygon that is
            # essentially pol minus the intersection
            intersection = pol.intersection(pol2)
            nonoverlap = pol.difference(intersection)
            outmulti.append(nonoverlap)

        else:
            # Otherwise, just keep the initial polygon as it is.
            outmulti.append(pol)

finalpol = MultiPolygon(outmulti)

Answer 1

我想您可以在这两个多边形之间使用symmetric_difference，通过与第二个多边形的差值组合来实现您想要做的事情（对称差值将为您带来两个多边形中不重叠的部分，在这些部分上通过差值删除了多边形2的部分）。我还没有测试过，但看起来像：

# multipol1 and multipol2 are my shapely MultiPolygons
from shapely.ops import cascaded_union
from itertools import combinations
from shapely.geometry import Polygon,MultiPolygon

outmulti = []
for pol in multipoly1:
    for pol2 in multipoly2:
        if pol.intersects(pol2)==True:
            # If they intersect, create a new polygon that is
            # essentially pol minus the intersection
            nonoverlap = (pol.symmetric_difference(pol2)).difference(pol2)
            outmulti.append(nonoverlap)

        else:
            # Otherwise, just keep the initial polygon as it is.
            outmulti.append(pol)

finalpol = MultiPolygon(outmulti)