如何在R语言中准确计算国家间的国界线长度？

问题概述 #

在以国界线长度衡量国家文化距离的场景中，使用sf包处理naturalearthdata的世界shapefile时遇到两个问题：

1. 默认调用st_intersection()无法获取丹麦与德国的边界交集；
2. 设置model="closed"后得到的边界长度约111公里，与实际68公里偏差显著。

核心解决方案 #

1. 切换到平面几何处理（解决无交集问题）

当前环境启用了S2球面几何（sf_use_s2() is TRUE），球面几何对多边形边界的交集逻辑与平面几何存在差异，导致默认调用返回空结果。只需切换到GEOS平面几何即可解决：

suppressWarnings(library(sf))
sf_use_s2(FALSE)  # 关闭S2球面几何，启用GEOS平面几何

# 读取并预处理数据
world = read_sf("data/ne_10m_admin_0_countries/ne_10m_admin_0_countries.shp")
world = st_make_valid(world)

# 获取目标国家索引
germany_idx = which(world$ADMIN=="Germany")
denmark_idx = which(world$ADMIN=="Denmark")

# 计算边界交集
border = st_intersection(world$geometry[germany_idx], world$geometry[denmark_idx])

2. 转换到投影坐标系（解决长度偏差）

原数据使用WGS84地理坐标系（EPSG:4326），st_length()在此坐标系下计算的是球面距离，会产生系统偏差。需转换到适合欧洲区域的投影坐标系（如UTM 32N，EPSG:32632）后再计算长度：

# 转换到UTM 32N投影
world_proj = st_transform(world, crs = 32632)
germany_proj = world_proj$geometry[germany_idx]
denmark_proj = world_proj$geometry[denmark_idx]

# 重新计算交集并获取长度
border_proj = st_intersection(germany_proj, denmark_proj)
st_length(border_proj)  # 结果接近68000米，与实际值一致

潜在问题排查 #

• 海岸线悖论：此场景计算的是陆地国界线，而非海岸线，且naturalearth 10m精度数据已做简化，不会因细分程度导致大幅偏差；
• shapefile数据问题：naturalearth的10m行政边界数据可靠性较高，结合tmap绘图验证无异常，排除数据本身问题；
• 参数设置错误：model="closed"是球面几何下的参数，切换到平面几何后无需额外设置，该参数引入的水上边界计算是之前长度偏差的主要原因。

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