求助：为何CSV点位无法叠加显示在SHP流域图层上？

以下是针对该问题的分步排查和解决方法：

• 检查坐标参考系（CRS）匹配性
  流域SHP与CSV转换的空间对象必须使用完全一致的CRS，否则点位会因坐标系统偏差无法显示在流域范围内。

  1. 先查看流域SHP的CRS：

     st_crs(watershed_shp)
     

  2. 将CSV转换为sf对象时，先指定原始经纬度的CRS（通常为WGS84，EPSG:4326），再转换为与流域一致的CRS：

     points_sf <- st_as_sf(csv_data, coords = c("long", "lat"), crs = 4326) %>%
       st_transform(st_crs(watershed_shp))
     

  确认两者CRS的EPSG代码和描述完全一致，避免细微差异导致的偏移。

• 验证经纬度数据有效性
  检查CSV的lat/long列是否为数值型，是否存在异常值：

  # 查看数据类型
  str(csv_data$lat)
  str(csv_data$long)
  # 检查缺失值
  sum(is.na(csv_data$lat))
  sum(is.na(csv_data$long))
  

  如果是字符型，先转换为数值：

  csv_data$lat <- as.numeric(csv_data$lat)
  csv_data$long <- as.numeric(csv_data$long)
  

  同时确认经纬度顺序未搞反（很多点位不显示是因为把lat和long的顺序写反了，导致点落在地图范围外）。

• 确认点位是否在流域范围内
  查看点位和流域的边界范围是否重叠：

  st_bbox(points_sf)
  st_bbox(watershed_shp)
  

  如果两者范围完全不重叠，可单独绘制点位确认是否存在：

  plot(st_geometry(points_sf))
  

  若单独绘制能看到点位，说明坐标或CRS转换有误，导致点不在流域视图范围内。

• 检查绘图代码逻辑
  确保绘图时先绘制流域图层，再叠加点位，且点位样式足够明显，同时必须添加add = TRUE参数实现叠加：

  # 先画流域
  plot(st_geometry(watershed_shp), col = "lightblue")
  # 叠加点位，用醒目样式
  plot(st_geometry(points_sf), col = "red", pch = 19, cex = 1.5, add = TRUE)
  

  若之前代码遗漏add = TRUE，会导致点位覆盖流域图层或单独绘制在新窗口。

• 排查版本兼容性问题
  检查sf包版本是否与之前正常运行时一致，新版本可能存在语法变化：

  packageVersion("sf")
  

  若版本差异较大，尝试改用兼容旧版本的写法，例如：

  points_sfc <- st_sfc(st_multipoint(cbind(csv_data$long, csv_data$lat)), crs = 4326)
  points_sf <- st_sf(data = csv_data, geometry = points_sfc) %>% st_transform(st_crs(watershed_shp))
  

内容的提问来源于stack exchange，提问作者Lauren.S