Skip to content

Instantly share code, notes, and snippets.

@oscarperpinan
Created March 30, 2015 07:33
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save oscarperpinan/504ce14e320d21e1644c to your computer and use it in GitHub Desktop.
Save oscarperpinan/504ce14e320d21e1644c to your computer and use it in GitHub Desktop.
Respuesta a comentario en blog
library(sp)
## Grabar tabla excel en formato CSV usando ; como separador
## Leemos saltando dos primeras líneas
pts <- read.csv2('/tmp/listado-longitud-latitud-municipios-espana.csv', skip = 2)
## Transformamos a SPDF usando columnas de longitud-latitud como coordenadas
munPts <- SpatialPointsDataFrame(pts[, c(5, 4)], pts[,-c(4, 5)],
proj4string = CRS('+proj=longlat +datum=WGS84'))
## Pintamos el resultado para comprobar
spplot(munPts["Altitud"])
## Ahora se puede añadir a otra capa "trellis" existente
p + layer(sp.points(munPts))
@oscarperpinan
Copy link
Author

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment