R中SDM空间点与栅格掩膜不匹配问题求助

• 检查栅格掩膜的范围与分辨率
  WorldClim的30s分辨率栅格默认是全球范围，若观测点集中在局部区域，全局绘图会导致点被压缩到视觉角落。可裁剪栅格到观测点的边界范围后再验证：

  # 提取观测点的外接矩形范围
  extent_pt <- extent(pt)
  # 裁剪掩膜到该范围
  mask_SA_cropped <- crop(mask_SA, extent_pt)
  # 重新绘图验证匹配情况
  plot(mask_SA_cropped)
  plot(pt, pch = 16, col = "red", add = T)
  

• 修正掩膜的无数据值判断逻辑
  WorldClim栅格的标准无数据值是-9999而非-500，用> -500创建掩膜会包含大量无效区域，导致视觉上的不匹配。改用栅格自带的NA值判断：

  # 基于原栅格的NA值创建有效区域掩膜
  mask_SA <- !is.na(clim_SA$AnnMeanTemp)
  

• 排查观测数据的数值异常
  确认经纬度未颠倒后，需检查观测点是否存在超出合理范围的异常值（如经度不在-180_{180、纬度不在-90}90），这类异常点会导致绘图偏移：

  # 查看经纬度的数值范围
  summary(obs_data[, c("decimalLongitude", "decimalLatitude")])
  # 过滤异常值
  obs_data_clean <- obs_data[obs_data$decimalLongitude >= -180 & obs_data$decimalLongitude <= 180 &
                            obs_data$decimalLatitude >= -90 & obs_data$decimalLatitude <= 90, ]
  # 重新创建空间点
  pt <- SpatialPoints(obs_data_clean[, c("decimalLongitude", "decimalLatitude")], 
                      proj4string = CRS("+proj=longlat +datum=WGS84"))
  

• 更新依赖包并尝试sf空间对象
  R 4.3.2版本较新，旧版raster/sp包可能存在兼容性问题。建议更新raster包，同时尝试用更现代的sf包处理空间数据：

  # 安装并加载sf包
  install.packages("sf")
  library(sf)
  # 将观测数据转换为sf空间点对象（EPSG:4326对应WGS84）
  pt_sf <- st_as_sf(obs_data, coords = c("decimalLongitude", "decimalLatitude"), crs = 4326)
  # 绘图验证
  plot(mask_SA)
  plot(st_geometry(pt_sf), pch = 16, col = "red", add = T)
  

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