我发布了此question on the GIS stack exchange，但没有看到太多的流量。
我是一名GIS用户，几年来一直使用R进行统计，我很高兴看到很多新的GIS功能被发布（光栅，shapefile，rgdal等）。
我在R中做了很多数据库和表操作，因此能够从shapefile添加和删除属性是一个强大的潜力。
我希望我只是错过了那里的内容，但我似乎无法找到一种在shapefile属性表中添加或删除属性的好方法。
Overflow中的任何人都可以回复我的GIS帖子吗？或者我可以被告知如何将我的列表属性表放入数据帧并返回以替换当前属性表？
自我原帖后的一些进展：
这是在弄清楚如何使用我的属性表appart（.dbf），添加内容，现在我试图将它重新组合起来替换原来的dbf。
>libary(raster); library(rgdal); library(shapefiles)
>shp<-shapefile(Shape) # D.C. area airport polygons
>summary(shp) #Shapefile properties
Object of class SpatialPointsDataFrame
Coordinates:
min max
coords.x1 281314.2 337904.7
coords.x2 4288867.0 4313507.0
Is projected: TRUE
proj4string :
[+proj=utm +zone=18 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0]
Number of points: 4
Data attributes:
ObjectID LOCID NAME FIELD STATE STATE_FIPS ACAIS TOT_ENP TYPE
Min. :134.0 Length:4 Length:4 Length:4 Length:4 Length:4 Length:4 Min. : 271 Length:4
1st Qu.:242.8 Class :character Class :character Class :character Class :character Class :character Class :character 1st Qu.: 3876 Class :character
Median :339.0 Mode :character Mode :character Mode :character Mode :character Mode :character Mode :character Median :3409113 Mode :character
Mean :310.0 Mean :3717251
3rd Qu.:406.2 3rd Qu.:7122488
Max. :428.0
>shp.AT<-read.dbf(gsub(".shp", ".dbf", Shape), header=TRUE) #Read in the attribute table from the .dbf
> shp.AT # First object in the dbf looks like an attribute table!
$dbf
ObjectID LOCID NAME FIELD STATE STATE_FIPS ACAIS TOT_ENP TYPE
1 134 ADW Andrews AFB <NA> Maryland 24 Y 5078 Military
2 279 DCA Washington National <NA> Virginia 51 Y 6813148 <NA>
3 399 HEF Manassas Regional Harry P Davis Field Virginia 51 Y 271 Regional
4 428 IAD Washington Dulles International <NA> Virginia 51 Y 8050506 International
$header
$header$file.version
[1] 3
$header$file.year
[1] 113
$header$file.month
[1] 4
$header$file.day
[1] 12
$header$num.records
[1] 4
$header$header.length
[1] 321
$header$record.length
[1] 148
$header$fields
NAME TYPE LENGTH DECIMAL
1 ObjectID N 10 0
2 LOCID C 5 0
3 NAME C 45 0
4 FIELD C 30 0
5 STATE C 24 0
6 STATE_FIPS C 2 0
7 ACAIS C 1 0
8 TOT_ENP N 11 0
9 TYPE C 20 0
>shp.tab<-as.data.frame(shp.AT[1]) # Grab the first object of the .dbf as a data.frame
> shp.tab # First list object
dbf.ObjectID dbf.LOCID dbf.NAME dbf.FIELD dbf.STATE dbf.STATE_FIPS dbf.ACAIS dbf.TOT_ENP dbf.TYPE
1 134 ADW Andrews AFB <NA> Maryland 24 Y 5078 Military
2 279 DCA Washington National <NA> Virginia 51 Y 6813148 <NA>
3 399 HEF Manassas Regional Harry P Davis Field Virginia 51 Y 271 Regional
4 428 IAD Washington Dulles International <NA> Virginia 51 Y 8050506 International
> shp.tab$NewAT<-1:nrow(shp.tab) # Add my new attribute
> shp.tab # Added my new attribute, now to get this back into my shapefile
dbf.ObjectID dbf.LOCID dbf.NAME dbf.FIELD dbf.STATE dbf.STATE_FIPS dbf.ACAIS dbf.TOT_ENP dbf.TYPE NewAT
1 134 ADW Andrews AFB <NA> Maryland 24 Y 5078 Military 1
2 279 DCA Washington National <NA> Virginia 51 Y 6813148 <NA> 2
3 399 HEF Manassas Regional Harry P Davis Field Virginia 51 Y 271 Regional 3
4 428 IAD Washington Dulles International <NA> Virginia 51 Y 8050506 International 4
>write.dbf(shp.tab, gsub(".shp", ".dbf", Shape)) # Knew this wouldn't work, but demonstrate attempt to write this object as .dbf.
ERROR:
invalid subscript type 'list'
> shp.AT[1]<-shp.tab # Try replacing the old Object[1] with my new table containing the new attribute.
> shp.AT # The table portion fo the shp.AT is gone. No attributes.
$dbf
[1] 134 279 399 428
$header
$header$file.version
[1] 3
$header$file.year
[1] 113
$header$file.month
[1] 4
$header$file.day
[1] 12
$header$num.records
[1] 4
>write.dbf(shp.AT, gsub(".shp", ".dbf", Shape)) # If I go ahead and overwrite anyway...
我的属性已消失，并替换为包含单个字段“dataframe”的属性表。 我递归脚本再次读取属性表。
> shp.tab
dataframe
1 134
2 279
3 399
4 428
所以，我认为我很接近。任何人都可以帮助我重新进入shapefile属性表吗？ 或者，有没有更好的方法？
由于
答案 0 :(得分：5)
我不确定我完全明白你要做什么。看起来你只想在属性表中添加一个新列？如果这是正确的，那么就像对待任何数据帧一样对待它。
library(rgdal)
dsn <- system.file("vectors", package = "rgdal")
shp<-readOGR(dsn = dsn, layer = 'cities')
shp$NewAT<-1:nrow(shp)
这与我系统上的shapefile完美配合。我通常依靠rgdal来读取我的shapefile，使用readOGR（）函数。我很确定你调用的shapefile（）函数也调用了rgdal。
已编辑以添加可重现的数据集。
答案 1 :(得分：0)
Iam使用R进行GIS扩展已经有好几年了，通常与QGIS结合使用。为了处理属性表，我通常将形状另存为CSV，将几何另存为WKT（您可以在QGIS保存对话框中找到这些选项）。在下一步中，我读取了R中的所有csv（shapes）并进行了统计，加入等操作。最后，我将它们写回到HDD并将其加载回QGIS（只需拖放即可，无需导入对话框）并保存它们作为shapefile。
HTH， 乔
答案 2 :(得分：0)
我使用 sf 库和 tidygraph 来处理 shapefile。
library(sf)
Links_Layer <- sf::st_read("D:\\Sample\\Folder", layer = "shapefile", quiet = TRUE)
Links_Layer <- Links_Layer %>% mutate(Attribute_B = Attribute_A^2)
sf::st_write(Links_Layer, "D:\\Sample\\Folder\\updated_shapefile.shp")