Anda perlu beralih di beberapa level. ( Pembaruan : Saya telah mengedit untuk menghapus semua loop "untuk", kecuali untuk satu pemahaman daftar )
# imports used throughout this example
from shapely.geometry import Point
from shapely.ops import cascaded_union
from itertools import combinations
# Here are your input shapes (circles A, B, C)
A = Point(3, 6).buffer(4)
B = Point(6, 2).buffer(4)
C = Point(1, 2).buffer(4)
# list the shapes so they are iterable
shapes = [A, B, C]
Pertama, Anda membutuhkan penyatuan semua persimpangan (menggunakan penyatuan bertingkat ), menggunakan pasangan kombinasi dari setiap bentuk. Kemudian Anda menghapus (melalui difference
) persimpangan dari penyatuan semua bentuk.
# All intersections
inter = cascaded_union([pair[0].intersection(pair[1]) for pair in combinations(shapes, 2)])
# Remove from union of all shapes
nonoverlap = cascaded_union(shapes).difference(inter)
Berikut ini nonoverlap
tampilannya (melalui JTS Test Builder):