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

新闻中心 NEWS CENTER

您当前位置：首页 > 新闻中心 > 网络学院

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

2025-11-02

浏览次数：次

返回列表

使用 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 创建多边形对象

Musho

AI网页设计Figma插件

76 查看详情 Musho

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

# 创建多边形对象
t1 = sg.Polygon([(-74.418225663382, 39.36239030236737), # Atlantic City, NJ
                 (-74.27880733397238, 39.71055595453288),
                 (-74.75681303480502, 40.219387193292164),
                 (-75.4705021020258, 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.4705021020258, 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, "海里")