Comment superposer des couches de carte (shp et csv) dans R?

Il s'agit d'un crosspost de stackoverflow . Je connais très peu de logiciels SIG et je fais tout mon mapping dans R. Excuses à l'avance si c'est une question trop basique. Disons que j'ai deux fichiers de formes provenant de sources différentes mais avec des attributs différents. Disons, l'un est pour les limites administratives du Texas ( boundaries.shp) et l'autre pour les rivières du Texas ( rivers.shp). J'ai également un troisième dossier towns.csvmontrant les emplacements des villes dans l'état. Après avoir lu les fichiers, je peux surplacer les emplacements des villes sur les limites administratives dans le maptoolspackage:

plot(boundaries); points(towns$lon, towns$lat)

Mais comment puis-je superposer les trois? Il y a sûrement un moyen facile de le faire?

PBSMapping devrait répondre à vos besoins. Il y a un tutoriel chez NCEAS . Le code ci-dessous est adapté de ce tutoriel. Je fais des hypothèses sur vos données btw. Veuillez modifier selon votre situation.

library(PBSmapping)

#prepare towns
pts <- read.csv("towns.csv")
towns <- points(towns$lon, towns$lat)
# read in shapefiles 
rivers <- importShapefile("rivers.shp")
boundaries <- importShapefile("boundaries.shp")


# note that importShapefile reads the .prj file if it exists, but it
# does not adopt the proj4 format used by the above approaches
proj.abbr <- attr(boundaries, "projection") # abbreviated projection info
proj.full <- attr(boundaries, "prj") # full projection info
print(proj.abbr)
# [1] "LL"

# generate map using PBSmapping plotting functions
plotPolys(boundaries, projection=proj.abbr, border="gray",
    xlab="Longitude", ylab="Latitude")
addPoints(towns, pch=20, cex=0.8)
addLines(rivers, col="blue", lwd=2.0)

La façon la plus simple de superposer deux tracés pourrait être d'utiliser l' add = TRUEoption dans plot. Voici un exemple avec des données artificielles

# Load sp package for creating artificial data
library(sp)

# Create sample town points
towns <- data.frame(lon = sample(100), lat = sample(100))
towns <- SpatialPoints(towns)

# Create sample polygon grid
grd <- GridTopology(c(1,1), c(10,10), c(10,10))
polys <- as.SpatialPolygons.GridTopology(grd)

# Plot polygons
plot(polys)

# Add towns (in red colour)
plot(towns, add = TRUE, col = 'red')