使用readOGR()读取的Shapefile空间对象属性值修改失败求助

这个错误是因为你尝试修改Spatial对象子集的@data槽位时，S4对象的特殊结构导致赋值逻辑出了问题。下面给你两种可行的解决方案：

方案1：修正rgdal的赋值方式 #

不要先索引整个Spatial对象再去修改@data，直接对原对象的@data字段使用逻辑索引赋值即可，这样能避开S4子集操作的陷阱：

library(rgdal)
shapefile <- readOGR(dsn="./shapefiles/concelhos-shapefile", layer = "concelhos", use_iconv = TRUE, encoding = "UTF-8")

# 直接操作@data中的列并添加逻辑索引
shapefile@data$NAME_2[shapefile@data$CCA_2 == '0308'] <- "Guimarães"
shapefile@data$NAME_2[shapefile@data$CCA_2 == '0814'] <- "Tavira"

这种方式直接针对属性表的列进行修改，不需要创建临时的Spatial对象子集，能有效避免类型转换错误。

方案2：切换到更现代的sf包（推荐） #

rgdal已经被官方标记为退役状态，现在R中处理空间数据的标准工具是sf包，它的语法更贴近常规的data.frame操作，不需要处理S4槽位，更不容易出错：

# 安装sf包（首次使用需要安装）
install.packages("sf")
library(sf)

# 读取Shapefile，语法更简洁
shapefile_sf <- st_read("./shapefiles/concelhos-shapefile/concelhos.shp", stringsAsFactors = FALSE)

# 修改属性字段，和操作普通data.frame完全一致
shapefile_sf$NAME_2[shapefile_sf$CCA_2 == '0308'] <- "Guimarães"
shapefile_sf$NAME_2[shapefile_sf$CCA_2 == '0814'] <- "Tavira"

sf包不仅解决了S4对象的操作痛点，还支持更多现代GIS功能，后续处理空间数据会更顺畅。

内容的提问来源于stack exchange，提问作者José Luís Borges