【图文并茂】使用Shapely计算点面关系&距离

已于 2023-11-15 22:18:41 修改
阅读量1.1k 收藏 4
点赞数 1
文章标签： python
于 2023-10-06 16:40:50 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ymzhu385/article/details/133613409
版权

文章目录


作者：小猪快跑

基础数学&计算数学，从事优化领域5年+，主要研究方向：MIP求解器、整数规划、随机规划、智能优化算法

本文档介绍如何使用 Shapely Python 包 计算几何点面关系&距离。

如有错误，欢迎指正。如有更好的算法，也欢迎交流！！！——@小猪快跑

相关文献

基础

Point、LineString、Polygon的通用属性

object.area   返回对象的面积
object.bounds 返回包含对象的x、y的最大值最小值的元组
object.length 返回对象的长度
object.geom_type 返回图形的类型 字符串
object.distance(other) 返回两个对象的最小距离
object.hausdorff_distance(other) 返回两个对象之间的哈多夫距离
object.representative_point() 返回一个保证在对象上的点

创建Point、LineString、Polygon：

from shapely import Point, Polygon, GeometryCollection,LineString

point = Point(3, 3)
polygon = Polygon([(0, 5), (1, 1), (3, 0)])
circ = Point(4, 0).buffer(2)
line = LineString([(0, 0), (2, 2)])
polygon.intersection(circ)

GeometryCollection([point, polygon, circ,line])

在这里插入图片描述

Point对象

point = Point(0, 1)  # 创建点对象，可以是三维的
print(point.area) # 0.0
print(point.length)  # 0.0
print(point.bounds)  # (0.0, 1.0, 0.0, 1.0)
print(point.geom_type)  # Point
print(point.coords)  # <shapely.coords.CoordinateSequence object at 0x00000204DDC2B370>
print(point.coords[:])  # [(0.0, 1.0)]
print(list(point.coords))  # [(0.0, 1.0)]
print(point.x, point.y) # 0.0 1.0

LineString对象

# 可以自己相交
line = LineString([(0, 0), (1, 1)])
print(line.area)  # 0.0
print(line.length)  # 1.4142135623730951
print(line.bounds)  # (0.0, 0.0, 1.0, 1.0)
print(line.geom_type)  # LineString
print(line.coords)  # <shapely.coords.CoordinateSequence object at 0x00000204DEA0C7F0>
print(line.coords[:], type(line.coords[:]))  # [(0.0, 0.0), (1.0, 1.0)] <class 'list'>

Polygon对象

# 可以传入两个参数，第一个是有序的外边缘轮廓，第二个是无序的内边缘轮廓，多边形内所有的环形最多只能有一个焦点，否则都是无效的
polygon = Polygon([(0, 0), (-1, 1), (2, 1), (3, -2)])
print(polygon.area)  # 5.0
print(polygon.length)  # 11.182042498005464
print(polygon.bounds)  # (-1.0, -2.0, 3.0, 1.0)
x_min, y_min, x_max, y_max = polygon.bounds
print(x_min, y_min, x_max, y_max)  # -1.0 -2.0 3.0 1.0
print(polygon.geom_type)  # Polygon
print(polygon.exterior.coords)  # <shapely.coords.CoordinateSequence object at 0x00000208A0C5D510>
print(polygon.exterior.coords[:])  # [(0.0, 0.0), (-1.0, 1.0), (2.0, 1.0), (3.0, -2.0), (0.0, 0.0)]
print(polygon.interiors)  # <shapely.geometry.polygon.InteriorRingSequence object at 0x00000208A0C5D510>
print(polygon.interiors[:])  # []

一元判定

has_z 判断是否有z坐标

p = Point(0, 0)
print(p.has_z)  # False
p = Point(0, 0, 0)
print(p.has_z)  # True

is_ccw 判断是否是逆时针的

lr = LinearRing([(1, 0), (1, 1), (0, 0)])
print(lr.is_ccw)  # True
lr = LinearRing([(1, 0), (0, 0), (1, 1)])
print(lr.is_ccw)  # False

is_ring 闭合返回TRUE

line = LineString([(0, 0), (1, 1), (1, 0)])
print(line.is_ring)  # False
line = LineString([(0, 0), (1, 1), (1, 0), (0, 0)])
print(line.is_ring)  # True
lr = LinearRing([(0, 0), (1, 1), (1, 0)])
print(lr.is_ring)  # True
lr = LinearRing([(0, 0), (1, 1), (1, 0), (0, 0)])
print(lr.is_ring)  # True

在这里插入图片描述

is_empty 返回TRUE如果内部和边界都是空

p = Point()
print(p.is_empty)  # True
p = Point(0, 0)
print(p.is_empty)  # False

is_simple 自己不相交，返回TRUE

line = LineString([(0, 0), (1, 1), (1, -1), (0, 1)])
print(line.is_simple)  # False
line = LineString([(0, 0), (0, 1), (1, 1), (1, 0)])
print(line.is_simple)  # True

在这里插入图片描述

二元判定

object.__eq__(other)  # 类型和坐标点相同返回TRUE,和坐标点顺序有关
object.equals(other)  # boundary, interior, and exterior完全相同返回TRUE，和坐标点顺序无关
object.almost_equals(other[, decimal=6])  # 近似相等返回TRUE decimal是小数点的位数完全相同,和坐标点的顺序有关
object.contains(other)  # other里没有点在object的exterior，且other的interior里至少有一个点在object的interior,和坐标点的顺序无关
object.crosses(other)  # object的interior与other的interior相交，且不包含他,该判断只适用于图形与线之间的判定
object.disjoint(other)  # object的interior和boundary 和other的interior和boundary都不想交 返回True，完全不相交才返回True，有一个点都不行
object.intersects(other)  # object的interior或者boundary和other的interior或者boundary相交 返回TRUE
object.overlaps(other)  # object的interior或者boundary和other的interior或者boundary相交，且不包含他, 返回TRUE
object.touches(other)  # other和object至少有一个共同的点，且他们的interior不相交，也就是说在边界上

contains

a = Point(1, 1)
b = Point(3, 3)
print(a.contains(b))
a = Point(2, 2)
b = Point(2, 2)
print(a.contains(b))
a = Point(2, 1)
b = LineString([(1, 2), (3, 2)])
print(a.contains(b))
a = Point(2, 2)
b = LineString([(1, 2), (3, 2)])
print(a.contains(b))
a = Point(3, 1)
b = Point(2, 2).buffer(1)
print(a.contains(b))
a = Point(2, 2)
b = Point(2, 2).buffer(1)
print(a.contains(b))
a = Point(2, 1)
b = Point(2, 2).buffer(1)
print(a.contains(b))
a = LineString([(1, 1), (3, 1)])
b = LineString([(1, 3), (3, 3)])
print(a.contains(b))
a = LineString([(1, 2), (3, 2)])
b = LineString([(2, 1), (2, 3)])
print(a.contains(b))
a = LineString([(1, 0.5), (3, 0.5)])
b = Point(2, 2).buffer(1)
print(a.contains(b))
a = LineString([(1, 1), (3, 1)])
b = Point(2, 2).buffer(1)
print(a.contains(b))
a = LineString([(1, 1), (2, 2)])
b = Point(2, 2).buffer(1)
print(a.contains(b))
a = Point(1, 1).buffer(1)
b = Point(3, 3).buffer(1)
print(a.contains(b))
a = Polygon([(2, 2), (2, 3), (3, 3), (3, 2)])
b = Polygon([(1, 1), (1, 2), (2, 2), (2, 1)])
print(a.contains(b))
a = Point(1, 1).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.contains(b))
a = Point(2, 2).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.contains(b))

在这里插入图片描述
在这里插入图片描述

crosses

a = Point(1, 1)
b = Point(3, 3)
print(a.crosses(b))
a = Point(2, 2)
b = Point(2, 2)
print(a.crosses(b))
a = Point(2, 1)
b = LineString([(1, 2), (3, 2)])
print(a.crosses(b))
a = Point(2, 2)
b = LineString([(1, 2), (3, 2)])
print(a.crosses(b))
a = Point(3, 1)
b = Point(2, 2).buffer(1)
print(a.crosses(b))
a = Point(2, 2)
b = Point(2, 2).buffer(1)
print(a.crosses(b))
a = Point(2, 1)
b = Point(2, 2).buffer(1)
print(a.crosses(b))
a = LineString([(1, 1), (3, 1)])
b = LineString([(1, 3), (3, 3)])
print(a.crosses(b))
a = LineString([(1, 2), (3, 2)])
b = LineString([(2, 1), (2, 3)])
print(a.crosses(b))
a = LineString([(1, 0.5), (3, 0.5)])
b = Point(2, 2).buffer(1)
print(a.crosses(b))
a = LineString([(1, 1), (3, 1)])
b = Point(2, 2).buffer(1)
print(a.crosses(b))
a = LineString([(1, 1), (2, 2)])
b = Point(2, 2).buffer(1)
print(a.crosses(b))
a = Point(1, 1).buffer(1)
b = Point(3, 3).buffer(1)
print(a.crosses(b))
a = Polygon([(2, 2), (2, 3), (3, 3), (3, 2)])
b = Polygon([(1, 1), (1, 2), (2, 2), (2, 1)])
print(a.crosses(b))
a = Point(1, 1).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.crosses(b))
a = Point(2, 2).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.crosses(b))

在这里插入图片描述

disjoint

a = Point(1, 1)
b = Point(3, 3)
print(a.disjoint(b))
a = Point(2, 2)
b = Point(2, 2)
print(a.disjoint(b))
a = Point(2, 1)
b = LineString([(1, 2), (3, 2)])
print(a.disjoint(b))
a = Point(2, 2)
b = LineString([(1, 2), (3, 2)])
print(a.disjoint(b))
a = Point(3, 1)
b = Point(2, 2).buffer(1)
print(a.disjoint(b))
a = Point(2, 2)
b = Point(2, 2).buffer(1)
print(a.disjoint(b))
a = Point(2, 1)
b = Point(2, 2).buffer(1)
print(a.disjoint(b))
a = LineString([(1, 1), (3, 1)])
b = LineString([(1, 3), (3, 3)])
print(a.disjoint(b))
a = LineString([(1, 2), (3, 2)])
b = LineString([(2, 1), (2, 3)])
print(a.disjoint(b))
a = LineString([(1, 0.5), (3, 0.5)])
b = Point(2, 2).buffer(1)
print(a.disjoint(b))
a = LineString([(1, 1), (3, 1)])
b = Point(2, 2).buffer(1)
print(a.disjoint(b))
a = LineString([(1, 1), (2, 2)])
b = Point(2, 2).buffer(1)
print(a.disjoint(b))
a = Point(1, 1).buffer(1)
b = Point(3, 3).buffer(1)
print(a.disjoint(b))
a = Polygon([(2, 2), (2, 3), (3, 3), (3, 2)])
b = Polygon([(1, 1), (1, 2), (2, 2), (2, 1)])
print(a.disjoint(b))
a = Point(1, 1).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.disjoint(b))
a = Point(2, 2).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.disjoint(b))

在这里插入图片描述

intersects

a = Point(1, 1)
b = Point(3, 3)
print(a.intersects(b))
a = Point(2, 2)
b = Point(2, 2)
print(a.intersects(b))
a = Point(2, 1)
b = LineString([(1, 2), (3, 2)])
print(a.intersects(b))
a = Point(2, 2)
b = LineString([(1, 2), (3, 2)])
print(a.intersects(b))
a = Point(3, 1)
b = Point(2, 2).buffer(1)
print(a.intersects(b))
a = Point(2, 2)
b = Point(2, 2).buffer(1)
print(a.intersects(b))
a = Point(2, 1)
b = Point(2, 2).buffer(1)
print(a.intersects(b))
a = LineString([(1, 1), (3, 1)])
b = LineString([(1, 3), (3, 3)])
print(a.intersects(b))
a = LineString([(1, 2), (3, 2)])
b = LineString([(2, 1), (2, 3)])
print(a.intersects(b))
a = LineString([(1, 0.5), (3, 0.5)])
b = Point(2, 2).buffer(1)
print(a.intersects(b))
a = LineString([(1, 1), (3, 1)])
b = Point(2, 2).buffer(1)
print(a.intersects(b))
a = LineString([(1, 1), (2, 2)])
b = Point(2, 2).buffer(1)
print(a.intersects(b))
a = Point(1, 1).buffer(1)
b = Point(3, 3).buffer(1)
print(a.intersects(b))
a = Polygon([(2, 2), (2, 3), (3, 3), (3, 2)])
b = Polygon([(1, 1), (1, 2), (2, 2), (2, 1)])
print(a.intersects(b))
a = Point(1, 1).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.intersects(b))
a = Point(2, 2).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.intersects(b))

在这里插入图片描述

overlaps

a = Point(1, 1)
b = Point(3, 3)
print(a.overlaps(b))
a = Point(2, 2)
b = Point(2, 2)
print(a.overlaps(b))
a = Point(2, 1)
b = LineString([(1, 2), (3, 2)])
print(a.overlaps(b))
a = Point(2, 2)
b = LineString([(1, 2), (3, 2)])
print(a.overlaps(b))
a = Point(3, 1)
b = Point(2, 2).buffer(1)
print(a.overlaps(b))
a = Point(2, 2)
b = Point(2, 2).buffer(1)
print(a.overlaps(b))
a = Point(2, 1)
b = Point(2, 2).buffer(1)
print(a.overlaps(b))
a = LineString([(1, 1), (3, 1)])
b = LineString([(1, 3), (3, 3)])
print(a.overlaps(b))
a = LineString([(1, 2), (3, 2)])
b = LineString([(2, 1), (2, 3)])
print(a.overlaps(b))
a = LineString([(1, 0.5), (3, 0.5)])
b = Point(2, 2).buffer(1)
print(a.overlaps(b))
a = LineString([(1, 1), (3, 1)])
b = Point(2, 2).buffer(1)
print(a.overlaps(b))
a = LineString([(1, 1), (2, 2)])
b = Point(2, 2).buffer(1)
print(a.overlaps(b))
a = Point(1, 1).buffer(1)
b = Point(3, 3).buffer(1)
print(a.overlaps(b))
a = Polygon([(2, 2), (2, 3), (3, 3), (3, 2)])
b = Polygon([(1, 1), (1, 2), (2, 2), (2, 1)])
print(a.overlaps(b))
a = Point(1, 1).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.overlaps(b))
a = Point(2, 2).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.overlaps(b))

在这里插入图片描述

touches

a = Point(1, 1)
b = Point(3, 3)
print(a.touches(b))
a = Point(2, 2)
b = Point(2, 2)
print(a.touches(b))
a = Point(2, 1)
b = LineString([(1, 2), (3, 2)])
print(a.touches(b))
a = Point(2, 2)
b = LineString([(1, 2), (3, 2)])
print(a.touches(b))
a = Point(3, 1)
b = Point(2, 2).buffer(1)
print(a.touches(b))
a = Point(2, 2)
b = Point(2, 2).buffer(1)
print(a.touches(b))
a = Point(2, 1)
b = Point(2, 2).buffer(1)
print(a.touches(b))
a = LineString([(1, 1), (3, 1)])
b = LineString([(1, 3), (3, 3)])
print(a.touches(b))
a = LineString([(1, 2), (3, 2)])
b = LineString([(2, 1), (2, 3)])
print(a.touches(b))
a = LineString([(1, 0.5), (3, 0.5)])
b = Point(2, 2).buffer(1)
print(a.touches(b))
a = LineString([(1, 1), (3, 1)])
b = Point(2, 2).buffer(1)
print(a.touches(b))
a = LineString([(1, 1), (2, 2)])
b = Point(2, 2).buffer(1)
print(a.touches(b))
a = Point(1, 1).buffer(1)
b = Point(3, 3).buffer(1)
print(a.touches(b))
a = Polygon([(2, 2), (2, 3), (3, 3), (3, 2)])
b = Polygon([(1, 1), (1, 2), (2, 2), (2, 1)])
print(a.touches(b))
a = Point(1, 1).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.touches(b))
a = Point(2, 2).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.touches(b))

在这里插入图片描述

计算距离

distance

from shapely.geometry import Point, LineString, Polygon

a = Point(1, 1)
b = Point(3, 3)
print(a.distance(b))
a = Point(2, 2)
b = Point(2, 2)
print(a.distance(b))
a = Point(2, 1)
b = LineString([(1, 2), (3, 2)])
print(a.distance(b))
a = Point(2, 2)
b = LineString([(1, 2), (3, 2)])
print(a.distance(b))
a = Point(3, 1)
b = Point(2, 2).buffer(1)
print(a.distance(b))
a = Point(2, 2)
b = Point(2, 2).buffer(1)
print(a.distance(b))
a = Point(2, 1)
b = Point(2, 2).buffer(1)
print(a.distance(b))
a = LineString([(1, 1), (3, 1)])
b = LineString([(1, 3), (3, 2)])
print(a.distance(b))
a = LineString([(1, 2), (3, 2)])
b = LineString([(2, 1), (2, 3)])
print(a.distance(b))
a = LineString([(1, 0.5), (3, 0.5)])
b = Point(2, 2).buffer(1)
print(a.distance(b))
a = LineString([(1, 1), (3, 1)])
b = Point(2, 2).buffer(1)
print(a.distance(b))
a = LineString([(1, 1), (2, 2)])
b = Point(2, 2).buffer(1)
print(a.distance(b))
a = Point(1, 1).buffer(1)
b = Point(3, 3).buffer(1)
print(a.distance(b))
a = Polygon([(2, 2), (2, 3), (3, 3), (3, 2)])
b = Polygon([(1, 1), (1, 2), (2, 2), (2, 1)])
print(a.distance(b))
a = Point(1, 1).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.distance(b))
a = Point(2, 2).buffer(1)
b = Point(2, 2).buffer(1.5)
print(a.distance(b))

在这里插入图片描述

点面关系图代码

# figures.py
from math import sqrt
from shapely import affinity

GM = (sqrt(5)-1.0)/2.0
W = 8.0
H = W*GM
SIZE = (W, H)

BLUE = '#6699cc'
GRAY = '#999999'
DARKGRAY = '#333333'
YELLOW = '#ffcc33'
GREEN = '#339933'
RED = '#ff3333'
BLACK = '#000000'


def add_origin(ax, geom, origin):
    x, y = xy = affinity.interpret_origin(geom, origin, 2)
    ax.plot(x, y, 'o', color=GRAY, zorder=1)
    ax.annotate(str(xy), xy=xy, ha='center',
                textcoords='offset points', xytext=(0, 8))


def set_limits(ax, x0, xN, y0, yN):
    ax.set_xlim(x0, xN)
    ax.set_xticks(range(x0, xN+1))
    ax.set_ylim(y0, yN)
    ax.set_yticks(range(y0, yN+1))
    ax.set_aspect("equal")

# main.py
import matplotlib.pyplot as plt
from shapely.geometry import Point, LineString, Polygon
from shapely.plotting import plot_polygon, plot_points, plot_line
from figures import BLUE, GRAY, set_limits

fig = plt.figure(1, figsize=(9, 9), dpi=300)
fig.subplots_adjust(wspace=0.5, hspace=0.5)  # 调整边距和子图的间距

ax = fig.add_subplot(4, 4, 1)
a = Point(1, 1)
b = Point(3, 3)
plot_points(a, ax=ax, color=GRAY)
plot_points(b, ax=ax, color=GRAY)
plot_points(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 2)
a = Point(2, 2)
b = Point(2, 2)
plot_points(a, ax=ax, color=GRAY)
plot_points(b, ax=ax, color=GRAY)
plot_points(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 3)
a = Point(2, 1)
b = LineString([(1, 2), (3, 2)])
plot_points(a, ax=ax, color=GRAY)
plot_line(b, ax=ax, color=GRAY)
plot_points(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 4)
a = Point(2, 2)
b = LineString([(1, 2), (3, 2)])
plot_points(a, ax=ax, color=GRAY)
plot_line(b, ax=ax, color=GRAY)
plot_points(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 5)
a = Point(3, 1)
b = Point(2, 2).buffer(1)
plot_points(a, ax=ax, color=GRAY)
plot_polygon(b, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_points(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 6)
a = Point(2, 2)
b = Point(2, 2).buffer(1)
plot_points(a, ax=ax, color=GRAY)
plot_polygon(b, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_points(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 7)
a = Point(2, 1)
b = Point(2, 2).buffer(1)
plot_points(a, ax=ax, color=GRAY)
plot_polygon(b, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_points(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 8)
a = LineString([(1, 1), (3, 1)])
b = LineString([(1, 3), (3, 3)])
plot_line(a, ax=ax, color=GRAY)
plot_line(b, ax=ax, color=GRAY)
plot_points(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 9)
a = LineString([(1, 2), (3, 2)])
b = LineString([(2, 1), (2, 3)])
plot_line(a, ax=ax, color=GRAY)
plot_line(b, ax=ax, color=GRAY)
plot_points(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 10)
a = LineString([(1, 0.5), (3, 0.5)])
b = Point(2, 2).buffer(1)
plot_line(a, ax=ax, color=GRAY)
plot_polygon(b, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_points(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 11)
a = LineString([(1, 1), (3, 1)])
b = Point(2, 2).buffer(1)
plot_line(a, ax=ax, color=GRAY)
plot_polygon(b, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_line(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 12)
a = LineString([(1, 1), (2, 2)])
b = Point(2, 2).buffer(1)
plot_line(a, ax=ax, color=GRAY)
plot_polygon(b, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_line(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 13)
a = Point(1, 1).buffer(1)
b = Point(3, 3).buffer(1)
plot_polygon(a, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_polygon(b, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_points(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 14)
a = Polygon([(2, 2), (2, 3), (3, 3), (3, 2)])
b = Polygon([(1, 1), (1, 2), (2, 2), (2, 1)])
plot_polygon(a, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_polygon(b, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_points(a.intersection(b), ax=ax, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 15)
a = Point(1, 1).buffer(1)
b = Point(2, 2).buffer(1.5)
plot_polygon(a, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_polygon(b, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_polygon(a.intersection(b), ax=ax, add_points=False, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

ax = fig.add_subplot(4, 4, 16)
a = Point(2, 2).buffer(1)
b = Point(2, 2).buffer(1.5)
plot_polygon(a, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_polygon(b, ax=ax, add_points=False, color=GRAY, alpha=0.2)
plot_polygon(a.intersection(b), ax=ax, add_points=False, color=BLUE)
ax.set_title(f'a.touches(b):{a.touches(b)}')
set_limits(ax, 0, 4, 0, 4)

plt.savefig('output.png')
plt.show()
小猪快跑爱摄影
关注
  • 1
    点赞
  • 踩
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
shapely包的使用
qq_27009517的博客
04-27 1227
1. Point,LineStrings 共有的变量和方法 Point LineStrings 引入包 from shapely.geometry import Point from shapely.geometry import LineString 共有的变量和方法 object.area 　　Returns the area (float) of the object. object...
shapely_dll.rar
10-25
windows下pip安装shapely库时，需要将geos.dll与geos_c.dll放置到C:\Windows\System32下，即可安装成功。
Python模块之Shapely
热门推荐
Python打杂工程师
06-01 3万+
Shapely通过Python的ctypes模块，对平面特征进行集合理论分析和操作。进行地理信息系统研究、数据处理、空间分析的常用模块
python 笔记：shapely（形状篇）
最新发布
qq_40206371的博客
02-12 2909
主要是点（point）、线（linestring）、面（surface）
距离(distance)算法小结
xiao1_1bing的博客
12-04 2万+
18种和“距离(distance)”、“相似度(similarity)”相关的量的小结 在计算机人工智能领域，距离(distance)、相似度(similarity)是经常出现的基本概念，它们在自然语言处理、计算机视觉等子领域有重要的应用，而这些概念又大多源于数学领域的度量(metric)、测度(measure)等概念。  这里拮取其中18种做下小结备忘，也借机熟悉markdown的数学公式语法...
shapely库的用法，高效处理点、线、面的几何关系和相关延申（GeoPandas 库）python
qq_44442727的博客
08-08 4849
它提供了一套丰富的函数和方法，用于读取、写入、操作和分析地理空间数据，包括点、线、多边形等几何对象。以上是一些常用的多边形处理方法，你可以根据具体需求使用 Shapely 提供的函数和方法来处理你的多边形数据。以上仅是 GeoPandas 库的一些基础用法示例，该库还提供了更多强大的功能，如地理空间分析、数据合并、投影转换等。: Shapely 是一个优秀的 Python 库，用于进行几何对象的创建、操作和分析。以上仅是 Shapely 库的一些基础用法示例，该库还提供了丰富的几何操作和分析功能。
python shapely.geometry.polygon任意两个四边形的IOU计算实例
09-17
这可以通过计算所有点的平均经度和纬度得到中心点，然后找出与中心点距离最远的点来实现。这部分未在提供的代码中详细展开，但可以通过类似的方法实现，例如计算所有点到中心点的欧氏距离，并找出最大值。 总结来说...
python Shapely使用指南详解
09-17
主要介绍了python Shapely使用指南详解,需要的朋友可以参考下
Python求两点之间的直线距离(2种实现方法)
09-19
直线长度类则通过两个点的坐标差计算距离，利用`math.sqrt()`函数求平方根，从而得到距离。以下是如何使用这两个类的示例： ```python import math class Point: def __init__(self, x=0, y=0): self.x = x ...
shapely 1.7.0
09-20
这些对象支持丰富的操作，如面积计算、距离测量、几何对象的相交、合并、分割等，为地理信息系统的开发提供了强大支持。 Shapely库的另一个关键特性是对WKT（Well-Known Text）和WKB（Well-Known Binary）两种空间...
Python Shapely指导手册：Python Shapely Documentation Release 1.6.4
12-20
Shapely是一个Python包，使用GEOS库中的函数进行矢量数据构建与分析，Shapely指导手册包含该包中函数功能的介绍与使用，适合相关分析人员参考使用。
shapely and geos break在distance方法
weixin_34337381的博客
11-22 97
问题 from shapely.geometry import Point print Point(0,0).distance(Point(1,1)) 报错 python: GeometryComponentFilter.cpp:34: virtual void geos::geom::GeometryComponentFilter::filter_ro(const geos:...
python寻找距离最近的点_python – 在Shapely中查找多边形上最近点的坐标
weixin_39542340的博客
12-08 3926
说我有以下多边形和点：>>> poly = Polygon([(0, 0), (2, 8), (14, 10), (6, 1)])>>> point = Point(12, 4)我可以计算点到多边形的距离……>>> dist = point.distance(poly)>>> print(dist)2.49136439561…...
python shapely库 的简单使用（点，线，面的构建以及基本函数）
weixin_50880191的博客
08-31 4267
一、点 1.点的构造 from shapely.geometry import Point point = Point(0.0, 0.0) Point((0.0, 0.0)) 2. 点的区域、长度 Point.area Point.length 点的长度和区域面积都为0 3.点的边界框 Point.bounds 返回(minx, miny, maxx, maxy) (0.0, 0.0, 0.0, 0.0) 4.点的...
Shapely解决error
sunjingsujing的博客
01-07 702
https://blog.csdn.net/u010205128/article/details/81003845
python shapely
XINGBAIDE的博客
07-27 706
https://shapely.readthedocs.io/en/latest/manual.html
python基础编程：python Shapely使用指南详解
python基础编程开发教程
03-22 9642
这篇文章主要介绍了python Shapely使用指南详解,需要的朋友可以参考下 Shapely是一个Python库，用于操作和分析笛卡尔坐标系中的几何对象。 引入包 from shapely.geometry import Point from shapely.geometry import LineString 共有的变量和方法 object.area 　　Returns the ar...
Python地理数据处理库shapely支持函数总结
日进一步
12-21 1万+
Shapely是一个Python库，用于操作和分析笛卡尔坐标系中的几何对象。 本文通过部分示例介绍了空间处理库Shape的部分概念与操作函数。 官方文档：https://shapely.readthedocs.io/en/latest/manual.html#introduction 1、函数列表 函数类型 算子名称 Shapely包提供...
shapely 计算两个点坐标之间的距离
07-17
在Shapely中，你可以使用`distance`方法来计算两个点之间的距离。以下是一个示例代码： ```python from shapely.geometry import Point # 创建两个点对象 point1 = Point(0, 0) point2 = Point(3, 4) # 计算两个点之间的距离 distance = point1.distance(point2) # 打印距离 print(distance) ``` 在这个示例中，我们首先导入`Point`类，然后创建了两个点对象`point1`和`point2`，其中`(0, 0)`和`(3, 4)`分别是点的坐标。 接下来，我们使用`distance`方法计算`point1`和`point2`之间的距离，并将结果保存在`distance`变量中。 最后，我们使用`print`语句打印出计算得到的距离。 请确保你已经安装了Shapely库，并按照你的实际需求替换点的坐标。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值