经纬度坐标系下terra与sf包面积计算准确性对比

核心差异原因 #

R中处理经纬度（EPSG:4326）面积计算的方式主要分为两类，差异源于是否考虑地球的球面特性：

• GEOS引擎（sf关闭S2/terra默认）：把经纬度当作平面坐标计算，属于平面几何近似，未考虑地球的球面曲率。这种方式计算速度快，但针对大范围、高纬度区域，误差会明显增大。
• S2几何库（sf默认启用）：基于球面几何模型直接计算，完全贴合地球的球面形状，能更精准反映真实地表面积，尤其适合跨大纬度、大跨度的区域计算。

结果对比与结论 #

从你的代码示例来看：

• 启用S2时，sf::st_area计算结果为166235 km²
• 关闭S2（GEOS）或使用terra::expanse，结果为165515 km²

两者的差异正是球面计算与平面近似的误差体现。S2的计算结果更准确，因为它修正了平面几何忽略地球曲率带来的偏差。如果你的业务场景对面积精度要求较高（比如大区域地理分析、自然资源核算等），建议使用sf启用S2的st_area；如果是小范围区域，两种方式差异极小，GEOS/terra的计算效率会更高。

验证代码 #

library(terra)
library(sf)

# 创建经纬度坐标系下的多边形
p <- vect('POLYGON ((2 2, 7 6, 4 9, 2 2))', crs = 'EPSG:4326')

# terra计算面积
expanse(p, unit = 'km')
# 输出：165515 km2

# sf启用S2计算面积
st_area(st_as_sf(p))/1e6
# 输出：166235 km2

# 关闭S2后，sf与terra结果一致
sf_use_s2(FALSE)
round((expanse(p, unit = 'km')) - (as.numeric(st_area(st_as_sf(p))/1e6)))
# 输出：0

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