在 R 中将数据框添加到空间多边形数据框

Question

我有两个文件，包括covid_provincial.csv和Co_adm1.zip文件。 第一个包括在一个国家/地区确认的 COVID 数量。 第二个是国家边界的shapefile。 我将它们导入 R 并将它们存储到a和b对象中，编写以下代码：

a<-read.csv("the directory of my system which includes covid_provincial.csv")

b<-readOGR(dsn = "the directory (the folder) of my system which includes shape files")

我想在b object 中添加一个名为 Province 的变量。 最后，我想将a合并到b 。

Answer 1

编辑

考虑到@Robert 的建议，最好省略对@data的显式调用。 您可以通过以下方式在 shp 文件中添加一列：

b$Province<-c("provinces", "names", "ecc")

或者像这样，如果你现在a中的省顺序对应于b的元素

b$Province<-a$Province

然后通过合并将a合并到b ：

b <-merge(b,a,by="Province")

这应该有效。 但是，有不同的方法可以做同样的事情。

编辑我添加了一个带有一些数据的示例。 这个例子是不可重现的，但应该让你朝着正确的方向前进：

library(raster)
shp <- shapefile("dati/Limiti_2016_WGS84_Italia/regioni/basilicata_provincie.shp")
plot(shp)

data.frame(shp)

  SHAPE_Leng SHAPE_Area   X_sum
0   593511.8 6593536923 4751953
1   367848.1 3479205085 1371840

创建一个玩具data.frame

df <- data.frame(provincie=c("A","B"),pop=c(10000,20000),covid=c(2000,3500))
shp$provincie <- df$provincie

现在 shp 看起来像这样

  SHAPE_Leng SHAPE_Area   X_sum provincie
0   593511.8 6593536923 4751953         A
1   367848.1 3479205085 1371840         B

最后，我们可以将两者merge

shp <- merge(shp,df,by="province")
 data.frame(shp)

output 是：

  provincie SHAPE_Leng SHAPE_Area   X_sum   pop covid
1         A   593511.8 6593536923 4751953 10000  2000
2         B   367848.1 3479205085 1371840 20000  3500

但shp的 class 被保留

> class(shp)
[1] "SpatialPolygonsDataFrame"
attr(,"package")
[1] "sp"