多边形最远坐标对（最大线性范围）查找及海里距离计算

需求：针对任意给定的多边形，找到距离最远的坐标对（Maximum Linear Extent），并计算这两个点之间的海里距离。

原代码问题分析 #

原代码存在以下关键错误：

• 导入geopy.distance as distance后，循环中用同名变量覆盖了该模块，导致无法调用地理距离计算方法。
• t1.exterior.coords返回的是元组类型的经纬度对，元组没有distance()方法，直接调用会抛出AttributeError。
• 原代码计算的是平面距离，而非地理坐标的大地距离，结果并非海里单位。

修正后的代码 #

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

# 创建多边形对象
t1 = sg.Polygon([
    (-74.418225663382, 39.36239030236737),  # 大西洋城，NJ
    (-74.27880733397238, 39.71055595453288),
    (-74.75681303480502, 40.219387193292164),
    (-75.4705021020208, 40.60356289498688),
    (-76.88460230031765, 40.264996135212186),  # 哈里斯堡，PA
    (-74.418225663382, 39.36239030236737)
])

# 初始化最远点对和最大距离
furthest_pair = None
max_nm_distance = 0.0

# 遍历所有顶点对（去掉最后一个重复的闭合点）
coords = list(t1.exterior.coords)[:-1]
for i, (lon1, lat1) in enumerate(coords):
    for (lon2, lat2) in coords[i+1:]:
        # 计算两点间的大地距离，直接获取海里单位结果
        distance = geo_dist.geodesic((lat1, lon1), (lat2, lon2)).nautical
        if distance > max_nm_distance:
            max_nm_distance = distance
            furthest_pair = ((lon1, lat1), (lon2, lat2))

# 输出结果
print("最远坐标对:", furthest_pair)
print("海里距离:", round(max_nm_distance, 1))
print("预期值: 126.0 NM")

关键改进说明 #

• 避免模块名冲突：将地理距离模块命名为geo_dist，防止循环变量覆盖导入的模块。
• 正确计算地理距离：使用geopy.distance.geodesic计算大地距离，通过.nautical属性直接获取海里单位结果。
• 优化顶点列表：移除多边形闭合时重复的最后一个顶点，减少无效计算。

大顶点量场景优化（可选） #

如果多边形顶点极多，遍历所有点对效率较低，可以通过最小旋转矩形快速定位最远点对：

# 优化方法：通过最小旋转矩形获取最远点对
min_rect = t1.minimum_rotated_rectangle
rect_coords = list(min_rect.exterior.coords)[:-1]

# 计算矩形四条边的海里距离，找到最长边
edge_distances = []
for i in range(4):
    lon1, lat1 = rect_coords[i]
    lon2, lat2 = rect_coords[(i+1)%4]
    dist = geo_dist.geodesic((lat1, lon1), (lat2, lon2)).nautical
    edge_distances.append((dist, (rect_coords[i], rect_coords[(i+1)%4])))

# 提取最长边的端点和距离
max_edge = max(edge_distances, key=lambda x: x[0])
print("优化方法-最远坐标对:", max_edge[1])
print("优化方法-海里距离:", round(max_edge[0], 1))

内容的提问来源于stack exchange，提问作者Tarak Pandya