In [2]:
from shapely.geometry import LineString, MultiLineString, Point
coords = [((0, 0), (1, 1)), ((-1, 0), (1, 0))]
lines = MultiLineString(coords)
lines.boundary
Out[2]:
获取多元要素的每个部分，在Shapely 2.0以前使用 lines.boundary ，在版本2.0后，使用下面的方法：
In [3]:
len(lines.geoms)
Out[3]:
2
In [4]:
for geom in lines.geoms:
print(geom)
LINESTRING (0 0, 1 1) LINESTRING (-1 0, 1 0)
In [5]:
lines.boundary.boundary
Out[5]:
In [6]:
lines.boundary.boundary.is_empty
Out[6]:
True
object.centroid
返回几何对象的质心的代表。
In [7]:
LineString([(0, 0), (1, 1)]).centroid
Out[7]:
In [8]:
LineString([(0, 0), (1, 1)]).centroid.wkt
Out[8]:
'POINT (0.5 0.5)'
注意：
一个对象的质心可能是它点中的一个，但是并不保证所有对象的质心都是这样。
In [9]:
from shapely.geometry import LineString
line = LineString([(0, 0), (1, 1), (0, 2), (2, 2), (3, 1), (1, 0)])
dilated = line.buffer(0.5)
eroded = dilated.buffer(-0.3)
In [10]:
dilated
Out[10]:
In [11]:
eroded
Out[11]:
一个点的默认缓冲区是与99.8%的圆盘面积的多边形补丁。
In [12]:
from shapely.geometry import Point
p = Point(0, 0).buffer(10.0)
len(p.exterior.coords)
Out[12]:
65
In [13]:
p.area
Out[13]:
313.6548490545941
分辨率为1，缓冲区是一个矩形补丁。例如：
In [14]:
q = Point(0, 0).buffer(10.0, 1)
len(q.exterior.coords)
Out[14]:
5
In [15]:
q
Out[15]:
In [16]:
q.area
Out[16]:
200.0
通过设定距离为0，
buffer 可清除自接触或者自交多边形，例如经典的“领结”。
In [17]:
from shapely.geometry import Polygon
coords = [(0, 0), (0, 2), (1, 1), (2, 2), (2, 0), (1, 1), (0, 0)]
bowtie = Polygon(coords)
bowtie.is_valid
Out[17]:
False
In [18]:
bowtie
Out[18]:
In [19]:
clean = bowtie.buffer(0)
clean.is_valid
Out[19]:
True
In [20]:
clean
Out[20]:
In [21]:
len(clean.geoms)
Out[21]:
2
In [22]:
list(clean.geoms[0].exterior.coords)
Out[22]:
[(0.0, 0.0), (0.0, 2.0), (1.0, 1.0), (0.0, 0.0)]
In [23]:
clean.geoms[0].exterior
Out[23]:
In [24]:
list(clean.geoms[0].exterior.coords)
Out[24]:
[(0.0, 0.0), (0.0, 2.0), (1.0, 1.0), (0.0, 0.0)]
In [25]:
clean.geoms[0].exterior.coords
Out[25]:
<shapely.coords.CoordinateSequence at 0x7f7a742ff860>
In [26]:
list(clean.geoms[0].exterior.coords)
Out[26]:
[(0.0, 0.0), (0.0, 2.0), (1.0, 1.0), (0.0, 0.0)]
In [27]:
Point(0, 0).convex_hull
Out[27]:
In [28]:
from shapely.geometry import MultiPoint
MultiPoint([(0, 0), (1, 1)]).convex_hull
Out[28]:
In [29]:
from shapely.geometry import MultiPoint
MultiPoint([(0, 0), (1, 1), (1, -1)]).convex_hull
Out[29]:
In [30]:
from shapely.geometry import Point
Point(0, 0).envelope
Out[30]:
In [31]:
from shapely.geometry import MultiPoint
MultiPoint([(0, 0), (1, 1)]).envelope
Out[31]:
In [32]:
p = Point(0.0, 0.0)
x = p.buffer(1.0)
x.area
Out[32]:
3.136548490545939
In [33]:
len(x.exterior.coords)
Out[33]:
65
In [34]:
s = x.simplify(0.05, preserve_topology=False)
s.area
Out[34]:
3.061467458920719
In [35]:
len(s.exterior.coords)
Out[35]:
17
In [36]:
a = Point(1, 1).buffer(1.5)
b = Point(2, 1).buffer(1.5)
a.difference(b)
Out[36]:
In [37]:
a = Point(1, 1).buffer(1.5)
b = Point(2, 1).buffer(1.5)
a.intersection(b)
Out[37]:
object.symmetric_difference(other) 返回一个在这个对象不在另一个几何对象的点的代表，
并且其它对象的点不在这个几何对象中。
In [38]:
a = Point(1, 1).buffer(1.5)
b = Point(2, 1).buffer(1.5)
a.symmetric_difference(b)
Out[38]:
In [39]:
a = Point(1, 1).buffer(1.5)
b = Point(2, 1).buffer(1.5)
a.union(b)
Out[39]:
这些操作的语义随着几何对象类型的不同而变化。例如， 比较合并之后多边形的边界，与多边形边界的合并结果。
In [40]:
a.union(b).boundary
Out[40]:
In [41]:
a.boundary.union(b.boundary)
Out[41]:
注意：
union() 是一个找到许多对象的累积集合的昂贵的方式。
shapely.ops.cascaded_union() 是一个更有效的方法。