sgo R包转换EPSG27700至WGS84的结果精度疑问

差异根源 #

EPSG 27700（英国国家格网，基于OSGB36基准面）与EPSG 4326（WGS84，谷歌地图采用的基准面）属于不同大地基准体系，转换精度差异来自使用的转换模型：

• 你使用的在线工具采用了OSTN15（英国地形测量局2015转换模型），这是英国官方提供的精细网格转换方案，能修正局部基准面偏移，精度可达厘米级。
• sgo包默认的sgo_transform函数采用的是简化的七参数Helmert变换，这是一种全局近似转换方法，在英国部分区域会产生亚米级偏差，这就是两者结果细微差异的原因。

实现高精度批量转换的方案 #

若要在R中得到与在线工具一致的高精度结果，推荐使用支持OSTN15转换的sf包（需依赖PROJ 6.0及以上版本），代码示例如下：

library(sf)

# 构造EPSG 27700坐标系的点数据
locations <- st_sfc(
  st_point(c(448098.4741, 271753.1604)),
  st_point(c(448028.2853, 271665.0542)),
  crs = 27700
)

# 调用OSTN15模型转换至WGS84
locations_wgs84 <- st_transform(locations, crs = 4326, options = "+proj=ostn15")

# 查看转换结果
print(locations_wgs84)

补充说明 #

• 若坚持使用sgo包，可尝试通过transform参数自定义OSTN15转换参数，但需要额外导入官方转换网格文件，操作复杂度远高于sf包。
• OSTN15是英国区域坐标转换的官方标准，转换结果完全匹配官方工具精度，适合批量处理大规模TPOs数据。

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