使用R中mapboxapi生成宽纬度地图时纬度偏移、数据错位问题求助

我正尝试在R语言中通过mapboxapi调用Mapbox生成跨宽纬度范围的地图，当前经度定位基本准确，但直布罗陀海峡（约36°纬度）的位置存在偏差；添加EPSG:4326格式的数据后，错位情况更加严重。调用Mapbox的代码如下：

basemap <- layer_static_mapbox(location = c(-20,10,20,55),
                               style_url = "mapbox://styles/mapbox/streets-v12",
                               username = "User",
                               access_token = "Token",
                               attribution = T,
                               logo = T)

而使用ggplot的borders绘制的地图能精准匹配数据与地理要素：

ggplot(data, aes(lon, lat)) +
  borders(xlim = c(-10, 10), ylim = c(10, 55)) +
  geom_hline(yintercept = 36, lty = 2, color = "cornflowerblue") +
  geom_point()

请问如何实现地图与数据的精准匹配？是否需要重投影？

核心问题与解决方案 #

问题根源在于Mapbox静态地图默认使用Web墨卡托投影（EPSG:3857），而你的数据是WGS84（EPSG:4326），两者坐标系统不一致导致错位，高纬度区域的拉伸变形会更明显。必须通过统一投影解决：

1. 将EPSG:4326数据转换为EPSG:3857
   使用sf包的st_transform()函数把数据投影转换为Web墨卡托，再与Mapbox底图叠加：

   library(sf)
   library(mapboxapi)
   library(ggplot2)
   
   # 非sf格式数据先转换为sf对象
   data_sf <- st_as_sf(data, coords = c("lon", "lat"), crs = 4326)
   # 转换为Web墨卡托投影
   data_3857 <- st_transform(data_sf, crs = 3857)
   
   # 生成Mapbox底图（默认EPSG:3857）
   basemap <- layer_static_mapbox(
     location = c(-20,10,20,55),
     style_url = "mapbox://styles/mapbox/streets-v12",
     username = "User",
     access_token = "Token",
     attribution = TRUE,
     logo = TRUE
   )
   
   # 叠加转换后的数据
   basemap +
     geom_sf(data = data_3857, color = "red", size = 2) +
     geom_hline(yintercept = st_transform(st_sfc(st_point(c(0, 36)), crs = 4326), 3857)[[1]][2], 
                lty = 2, color = "cornflowerblue")
   

2. 强制Mapbox底图使用EPSG:4326
   若你的Mapbox权限支持，可在layer_static_mapbox中指定crs = 4326，使底图与数据使用同一投影，无需转换数据：

   basemap <- layer_static_mapbox(
     location = c(-20,10,20,55),
     style_url = "mapbox://styles/mapbox/streets-v12",
     username = "User",
     access_token = "Token",
     attribution = TRUE,
     logo = TRUE,
     crs = 4326  # 指定WGS84投影
   )
   
   # 直接叠加EPSG:4326的数据
   basemap +
     geom_point(data = data, aes(lon, lat), color = "red", size = 2) +
     geom_hline(yintercept = 36, lty = 2, color = "cornflowerblue")
   

3. 注意事项

   • 跨宽纬度范围时，Web墨卡托投影在高纬度会有明显拉伸，若追求无变形的地理呈现，优先选EPSG:4326，但部分Mapbox样式可能不兼容该投影。
   • 确保location参数的范围格式与所用CRS匹配：EPSG:4326是c(min_lon, min_lat, max_lon, max_lat)，EPSG:3857则是对应的米制坐标。

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