Python shapely库的具体使用

from shapely.geometry import Point
point_1 = Point(1, 1)
# 传递的参数为坐标元组
point_2 = Point((2, 2))
print(point_1, point_2)
# 具体坐标
print('x =', point_1.x, 'y =', point_1.y)
# 长度面积
print('length =', point_1.length, 'area =', point_1.area)
# 边界框
print('bounds of point:', point_1.bounds)
# 坐标值
print('coords of point:', list(point_1.coords))

# 1.4142135623730951
point_1.distance(point_2)

from shapely.geometry import LineString, Point
line_1 = LineString([(0, 0), (0, 1), (1, 1)])
line_2 = LineString([Point(0, 0), Point(0, 1), Point(1, 1)])
print(line_1)
print(line_2)
# 长度面积
print('length =', line_1.length, 'area =', line_1.area)
# 边界框
print('bounds:', line_1.bounds)
# 坐标值-切片
print('coords:', line_1.coords[:])

from shapely.geometry import LinearRing, Point
ring_1 = LinearRing([(0, 0), (0, 1), (1, 1), (1, 0), (0, 0)])
# 起始点与终止点没有显式封闭
ring_2 = LinearRing([Point(0, 0), Point(0, 1), Point(1, 1), Point(1, 0)])
print(ring_1)
print(ring_2)
# 长度面积
print('length =', ring_2.length, 'area =', ring_2.area)
# 边界框
print('bounds:', ring_2.bounds)
# 坐标值
print('coords:', ring_2.coords[:])

from shapely.geometry import LinearRing, Point, Polygon
# 1. 坐标元组
polygon_1 = Polygon([(0, 0), (0, 1), (1, 1), (1, 0), (0, 0)])
# 2. Point
# 起始点与终止点没有显式封闭
polygon_2 = Polygon([Point(0, 0), Point(0, 1), Point(1, 1), Point(1, 0)])
# 3. LinearRing
polygon_3 = Polygon(LinearRing([(0, 0), (0, 1), (1, 1), (1, 0)]))
print("1.", polygon_1)
print("2.", polygon_2)
print("3.", polygon_3)
# 长度面积
print('length =', polygon_1.length, 'area =', polygon_1.area)
# 边界框
print('bounds:', polygon_1.bounds)
# 外环
print('exterior:', polygon_1.exterior)
# 外环的坐标值
print('coords of exterior:', polygon_1.exterior.coords[:])
# 所有内环
# 多边形没有"洞"时, 没有内环
print('interiors:', polygon_1.interiors[:])

from shapely.geometry import Polygon
ex = [(0, 0), (0, 3), (3, 3), (3, 0), (0, 0)]
hole = [(1, 1), (1, 2), (2, 2), (2, 1)]
polygon = Polygon(ex, holes=[hole])
# 长度面积
print('length =', polygon.length, 'area =', polygon.area)
# 边界框
print('bounds:', polygon.bounds)
# 外环
print('exterior:', polygon.exterior)
# 外环的坐标值
print('coords of exterior:', polygon.exterior.coords[:])
# 所有内环
print('interiors:', list(polygon.interiors))
# 第一个内环
print('first interior:', polygon.interiors[0])
# 第一个内环的坐标
print('coords of first interior:', polygon.interiors[0].coords[:])

from shapely.geometry import (
    LineString,
    MultiLineString,
    MultiPoint,
    MultiPolygon,
    Point,
    Polygon
)
# 坐标元组参数
mul_point_1 = MultiPoint([(1, 1), (2, 2)])
# 几何对象参数
mul_point_2 = MultiPoint([Point(1, 1), Point(2, 2)])
# 混合参数
mul_point_3 = MultiPoint([(1, 1), Point(2, 2)])
line_1 = LineString([(0, 0), (0, 1)])
line_2 = LineString([Point(1, 0), Point(1, 1)])
multi_line_1 = MultiLineString([[(0, 0), (0, 1)], [(1, 0), (1, 1)]])
multi_line_2 = MultiLineString([line_1, line_2])
polygon_1 = Polygon([(0, 0), (0, 1), (1, 1), (1, 0), (0, 0)])
polygon_2 = Polygon([Point(1, 1), Point(1, 2), Point(2, 2), Point(2, 1)])
# None 表示 holes=None
multi_polygon_1 = MultiPolygon([
    ([(0, 0), (0, 1), (1, 1), (1, 0), (0, 0)], None),
    ([(1, 1), (1, 2), (2, 2), (2, 1), (1, 1)], None)
])
multi_polygon_2 = MultiPolygon([polygon_1, polygon_2])

from shapely.geometry import GeometryCollection, LineString, Point, Polygon
line = LineString([Point(0, 1), Point(0, 2)])
polygon = Polygon([Point(1, 1), Point(1, 2), Point(2, 2), Point(2, 1)])
c = GeometryCollection([Point(0, 0), line, polygon])
print(c)
# 长度面积
print('length =', c.length, 'area =', c.area)
# 边界框
print('bounds:', c.bounds)

import shapely
from shapely import Point, wkb, wkt
p = Point(1, 1)
print("initial wkb:", wkb.dumps(p))
print("initial wkt:", wkt.dumps(p))

# 十六进制表示
print("wkb with hex:", wkb.dumps(p, hex=True))
# 移除多余的零
print("wkt trim zero:", wkt.dumps(p, trim=True))

print("to_wkb:", shapely.to_wkb(p))  # 等价于 p.wkb
print("to_wkt:", shapely.to_wkt(p))  # 等价于 p.wkt

print(p._repr_svg_())

import shapely
from shapely import wkb, wkt
wkb_i = '0101000000000000000000F03F000000000000F03F'
wkt_i = 'POINT (1 1)'
print(shapely.from_wkb(wkb_i))
print(shapely.from_wkt(wkt_i))
print(wkb.loads(wkb_i))
print(wkt.loads(wkt_i))