Joignez le fichier CSV au fichier de formes

Je voulais joindre un fichier CSV à un fichier de formes.

Lorsque j'ai identifié un champ pour faciliter la jointure, par exemple,

map <- spChFIDs(map, as.character(map$ID))

il est revenu

Erreur dans les spChFID (SP, x): les longueurs diffèrent

Quelqu'un peut-il conseiller?

Voici comment vous pouvez le faire dans R en utilisant sp::merge

library(raster)
# read data    
p <- shapefile("path/file.shp")
d <- read.csv("path/file.csv")

# merge on common variable, here called 'key'
m <- merge(p, d, by='key')

# perhaps save as shapefile again
shapefile(m, "path/merged.shp")

Je ne sais pas comment faire avec R pour joindre un fichier de formes avec un csv. Cependant, si cela peut vous aider, vous pouvez également le faire sur QGIS en utilisant le pluggin MMQGIS.

Si vous travaillez dans R, votre fichier de formes et .csv doivent être chargés en tant que data.frame ou data.table. Vous pouvez utiliser join () pour les fusionner.

Avez-vous plus d'informations sur le type de fichier de formes que vous possédez?

J'ai eu la chance de faire ce qui suit lorsque je fusionne des fichiers de formes et des données. Notez que vous devez laisser l'extension .shp dans l'argument "layer". J'ai utilisé "GEOID" où vous mettriez le nom de la "forme" en question.

library("rgdal")
library("data.table")
shapefile = readOGR(dsn = "DIRECTORY WITH SHAPEFILES", layer = "THE ACTUAL SHAPEFILE")
shapefile@data$id = rownames(shapefile@data)
shapefile.points = fortify(shapefile, region = "id")
shapefile.df = join(shapefile.points, shapefile@data, by = "id")
shapefile.df = subset(shapefile.df, select = c(long, lat, group, GEOID))
names(shapefile.df) = c("long", "lat", "group", "GEOID")

Ensuite, vous pouvez fusionner vos données en utilisant quelque chose comme

full.data = join(mydata, shapefile.df, by = "GEOID", type = "full")

Si vous travaillez dans R, la fusion n'est en fait pas nécessaire. Vous pouvez utiliser ggplot pour afficher vos données sur une carte. Voyez, mon exemple. Cette doublure fait des merveilles avec le fichier de forme bangladesh = fortify(p1, region = "ADM1_EN"). Et puis, vous pouvez utiliser

ggplot() + geom_map(data = br2016a, aes(map_id = ADM1_EN, fill = HectareLocal), 
                    map =bangladesh) +
  expand_limits(x = bangladesh$long, y = bangladesh$lat) + 
  coord_fixed(.96) +
  scale_fill_gradient(low="thistle2", high="darkred", 
                      guide="colorbar", na.value="white") +
  labs(title = "Bangladesh Boro Rice (Local variety), 2016", x = element_blank(), y = element_blank(), 
       fill='Hectares') +   theme(legend.position = "bottom") 

Et, la carte que je reçois est ici.