python 笔记：shapely（形状篇）

主要是点（point）、线（linestring）、面（surface）

1 基本方法和属性

2 Points

2.1 创建点

使用 Shapely 创建点时，可以直接提供 x 和 y 的坐标值，或者提供一个包含 x 和 y 值的元组

from shapely import Point

# 使用位置坐标值创建点
point = Point(0.0, 0.0)

# 使用点元组参数创建点
q = Point((0.0, 0.0))

2.2 属性和方法 

3 LineStrings

• 接受一个有序的2个或更多个 (x, y[, z]) 点元组序列
• 构造出的 LineString 对象代表点之间的一个或多个连接的线性样条。有序序列中允许重复的点，但可能会导致性能损失，应该避免。LineString 可以自我交叉

3.1 方法和属性

4 LinearRing

• 受一个有序的 (x, y[, z]) 点元组序列
• 序列可以通过在第一个和最后一个索引处传递相同的值来显式关闭。否则，序列将通过将第一个元组复制到最后一个索引来隐式关闭
• LinearRing 有零面积和非零长度
• 与 LineString 一样，有序序列中允许重复的点，但可能会导致性能损失，应该避免。LinearRing 不能自我交叉，也不能在单个点上接触自身

coords
area	始终为0
length	长度
bounds

5 polygon

• Polygon 类的构造函数接受两个位置参数。
  • 第一个参数是一个有序的 (x, y[, z]) 点元组序列，这与 LinearRing 的情况完全相同。
  • 第二个参数是一个可选的，无序的环状序列序列，指定了特征的内部边界或“洞”（同时在第一个LinearRing里面的部分挖掉，不在第一个

from shapely import Polygon

poly=Polygon([(0,2),(2,2),(2,0)],
             [[(0.5, 1.5), (1.5, 1.5), (1.5, 0.5), (0.5, 0.5)]])
poly

6 box

方形多边形

from shapely.geometry import box

box(0,0,1,2)

7 形状的集合

7.1 GeometryCollection

• 当两个几何对象的交集包含多种类型的几何元素时，例如 LineString 和 Point，Shapely 会使用 GeometryCollection 类型来表示这些结果

from shapely.geometry import LineString

a = LineString([(0, 0), (1, 1), (1,2), (2,2)])
b = LineString([(0, 0), (1, 1), (2,1), (2,2)])
x = a.intersection(b)
print(x)
#GEOMETRYCOLLECTION (LINESTRING (0 0, 1 1), POINT (2 2))

7.1.1 geoms

 通过 geoms 属性访问 GeometryCollection 中的成员

list(x.geoms)
#[<LINESTRING (0 0, 1 1)>, <POINT (2 2)>]

7.2 集合中的点：MultiPoint 类

MultiPoint 类用于表示点的集合。

构造函数接受一个点的序列，其中每个点可以是 (x, y[, z]) 形式的元组。

from shapely import MultiPoint
points = MultiPoint([(0.0, 0.0), (1.0, 1.0)])
print(points.area)  
# 输出 0.0
print(points.length)  
# 输出 0.0

7.3  集合中的线：MultiLineString 类

MultiLineString 类表示线的集合。构造函数接受一系列类似线的序列或对象。

一个 MultiLineString 对象的面积为零，长度为非零

from shapely import MultiLineString
coords = [((0, 0), (1, 1)), ((-1, 0), (1, 0))]
lines = MultiLineString(coords)
lines

7.4 多边形中的点 MultiPolygon

from shapely import Polygon
from shapely import MultiPolygon
polygon1 = Polygon([(0, 0), (1, 1), (1, 0),(0,1)])
polygon2 = Polygon([(0,0),(0,2),(1,1)])

multiPolygon=MultiPolygon([polygon1,polygon2])
multiPolygon

area
length
bounds
geoms