使用 Shapely 和 Geopy 计算多边形最远坐标及距离（海里）-Python教程-PHP中文网

使用 Shapely 和 Geopy 计算多边形最远坐标及距离（海里）

霞舞

发布： 2025-11-02 13:08:29

原创

462人浏览过

使用 shapely 和 geopy 计算多边形最远坐标及距离（海里）

本文档旨在指导如何使用 Python 的 Shapely 和 Geopy 库，计算给定多边形中最远的两个坐标点，并计算它们之间的距离，结果以海里为单位。我们将介绍必要的库、实现步骤以及示例代码，帮助读者理解并解决类似问题。

1. 简介

在地理空间数据处理中，经常需要计算多边形内的最远两点距离。这在诸如确定覆盖范围、优化路径规划等场景中非常有用。本教程将展示如何利用 Python 的 Shapely 库处理多边形几何对象，并结合 Geopy 库计算地理坐标之间的距离，最终以海里为单位输出结果。

2. 准备工作

在开始之前，请确保已安装以下 Python 库：

Shapely: 用于处理几何对象。
Geopy: 用于地理编码、逆地理编码和计算地理坐标之间的距离。

可以使用 pip 命令安装这些库：

pip install shapely geopy

登录后复制

3. 实现步骤

以下是计算多边形最远坐标及距离的详细步骤：

3.1 导入必要的库

首先，导入 shapely.geometry 用于创建多边形对象，并从 geopy.distance 导入 distance 函数，用于计算坐标之间的距离。

import shapely.geometry as sg
from geopy.distance import distance

登录后复制

3.2 创建多边形对象

千图设计室AI海报

千图网旗下的智能海报在线设计平台

查看详情

使用 shapely.geometry.Polygon 类创建一个多边形对象。多边形的顶点坐标以列表形式传递。

# 创建多边形对象
t1 = sg.Polygon([(-74.418225663382, 39.36239030236737), # Atlantic City, NJ
                 (-74.27880733397238, 39.71055595453288),
                 (-74.75681303480502, 40.219387193292164),
                 (-75.4705021020208, 40.60356289498688),
                 (-76.88460230031765, 40.264996135212186), # Harrisburg, PA
                 (-74.418225663382, 39.36239030236737)])

登录后复制

3.3 寻找最远坐标对

遍历多边形的所有顶点对，计算它们之间的距离。使用 geopy.distance.distance 函数计算距离，并通过 .nm 属性获取海里数。

# 初始化变量以存储最远坐标和距离
furthest_pair = None
max_distance = 0

# 遍历所有坐标对
for i, p1 in enumerate(t1.exterior.coords):
    for j, p2 in enumerate(t1.exterior.coords[i+1:]):
        dist = distance(p1, p2).nm  # 计算距离，单位为海里
        if dist > max_distance:
            max_distance = dist
            furthest_pair = (p1, p2)

登录后复制

3.4 输出结果

打印找到的最远坐标对以及它们之间的距离（海里）。

# 打印结果
print("最远坐标:", furthest_pair)
print("它们之间的距离:", max_distance, "海里")

登录后复制

4. 完整代码示例

import shapely.geometry as sg
from geopy.distance import distance

# 创建多边形对象
t1 = sg.Polygon([(-74.418225663382, 39.36239030236737), # Atlantic City, NJ
                 (-74.27880733397238, 39.71055595453288),
                 (-74.75681303480502, 40.219387193292164),
                 (-75.4705021020208, 40.60356289498688),
                 (-76.88460230031765, 40.264996135212186), # Harrisburg, PA
                 (-74.418225663382, 39.36239030236737)])

# 初始化变量以存储最远坐标和距离
furthest_pair = None
max_distance = 0

# 遍历所有坐标对
for i, p1 in enumerate(t1.exterior.coords):
    for j, p2 in enumerate(t1.exterior.coords[i+1:]):
        dist = distance(p1, p2).nm  # 计算距离，单位为海里
        if dist > max_distance:
            max_distance = dist
            furthest_pair = (p1, p2)

# 打印结果
print("最远坐标:", furthest_pair)
print("它们之间的距离:", max_distance, "海里")

登录后复制