R语言循环加载SpatialPolygonsDataFrame RData文件异常求助

你遇到的问题其实是load()函数的特性导致的：它默认会把RData文件里的对象直接加载到当前工作环境中，而返回值只是这个对象的名称字符串，不是对象本身。所以你代码里的kml变量存的只是加载的对象名字，不是实际的SpatialPolygonsDataFrame，自然没法用空间数据相关的函数处理。

下面给你两种可行的解决方法：

方法一：用get()获取加载后的对象 #

既然load()已经把对象加载到环境里了，我们可以用get()函数通过对象名称去获取实际的对象。修改你的循环代码如下：

polygons <- c('A','B','C','D') 
for (polygon_name in polygons) { 
  kml.path <- paste("~/KML/",polygon_name,".kml.RData",sep="") 
  # 先加载对象，再用get()拿到实际数据
  obj_name <- load(kml.path)
  kml <- get(obj_name)
  proj4string(coordinate.anagog) <- proj4string(kml) 
  points.in.polygons.output <- over(coordinate.anagog, kml)
  # 可选：如果怕环境里变量过多，处理完可以移除加载的对象
  rm(list = obj_name)
}

注意：如果你的每个RData文件里的对象名称都不一样，或者和文件名对应，这个方法没问题；但如果多个文件里的对象重名，后面加载的会覆盖前面的，所以处理完最好用rm()删掉临时对象。

方法二：改用saveRDS()和readRDS()（更推荐） #

这种方法从根源上避免了环境变量的问题，因为saveRDS()是把单个对象保存为文件，readRDS()读取时直接返回这个对象，不需要关心原来的对象名称，非常适合循环场景。

第一步：重新保存你的空间数据（如果还没这么做的话） #

假设你原来的SpatialPolygonsDataFrame对象是A_kml、B_kml等，用saveRDS()保存：

saveRDS(A_kml, file = "~/KML/A.kml.rds")
saveRDS(B_kml, file = "~/KML/B.kml.rds")
# 以此类推处理C和D

第二步：循环读取处理 #

修改你的循环代码为：

polygons <- c('A','B','C','D') 
for (polygon_name in polygons) { 
  kml.path <- paste("~/KML/",polygon_name,".kml.rds",sep="") 
  # 直接读取得到SpatialPolygonsDataFrame对象
  kml <- readRDS(kml.path)
  proj4string(coordinate.anagog) <- proj4string(kml) 
  points.in.polygons.output <- over(coordinate.anagog, kml)
}

这个方法更清晰，也不会在你的工作环境里留下多余的变量，是处理这类循环加载单个对象场景的最佳实践。

内容的提问来源于stack exchange，提问作者Alexandr Fruman