使用sf包绘制墨西哥经纬度点时位置偏移的解决方法

核心问题排查方向 #

导致点位错位的常见原因（尤其是对比geopandas正常的情况）：

• 初始sf对象创建时错误设置了CRS（直接用投影CRS解析经纬度坐标）
• 经纬度列顺序搞反（sf默认要求经度在前，纬度在后，geopandas对顺序的容错性更高）
• 底图与点位的CRS不统一

修正后的完整代码示例 #

# 安装依赖包（首次运行）
install.packages(c("sf", "ggplot2", "rnaturalearth"))

# 加载包
library(sf)
library(ggplot2)
library(rnaturalearth)

# 1. 读取CSV数据（假设列名为latitude和longitude，根据你的实际列名修改）
mex_points <- read.csv("mexico_points.csv")

# 2. 创建sf对象：先以WGS84（EPSG:4326）初始化（关键！经纬度属于地理坐标）
# 注意coords参数顺序：先经度（longitude），后纬度（latitude）
mex_sf <- st_as_sf(
  mex_points,
  coords = c("longitude", "latitude"),  # 必须匹配你的CSV列名
  crs = 4326  # 经纬度的标准CRS
)

# 3. 转换到你确定的墨西哥区域投影CRS
# 替换为你实际找到的EPSG码或proj4字符串（比如墨西哥常用的EPSG:6372）
target_crs <- 6372
mex_sf_proj <- st_transform(mex_sf, crs = target_crs)

# 4. 获取并转换墨西哥底图到同一CRS
mex_map <- ne_countries(country = "mexico", returnclass = "sf")
mex_map_proj <- st_transform(mex_map, crs = target_crs)

# 5. 绘图验证
ggplot() +
  geom_sf(data = mex_map_proj, fill = "#f0e68c", color = "black") +
  geom_sf(data = mex_sf_proj, color = "#ff0000", size = 2) +
  labs(title = "墨西哥点位分布图", x = "", y = "") +
  theme_minimal()

关键细节说明 #

1. CRS设置顺序：
   经纬度是地理坐标，必须先以WGS84（EPSG:4326）创建sf对象，再转换到目标投影CRS。如果直接用投影CRS初始化，sf会把经纬度数值当成投影坐标解析，导致完全错位。
2. 经纬度顺序：
   sf的coords参数严格要求经度在前，纬度在后。如果你的CSV是纬度列在前，一定要调整coords的顺序（比如c("lat", "lon")改成c("lon", "lat")）。
3. CRS统一性：
   底图和点位必须使用完全相同的CRS，否则会出现图层错位。

额外排查步骤 #

如果仍有问题：

• 用st_coordinates(mex_sf)检查点位坐标是否在墨西哥范围内（大致经度-118_{-86，纬度14}33）
• 用st_crs(mex_sf)和st_crs(mex_map_proj)确认CRS是否匹配
• 尝试交换CSV中经纬度列的顺序后重新创建sf对象

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