[R-es] maps

[Karel L.]

Hola Matías:

La verdad es que tener los puntos ayudaría ...

¿Es esto lo que quieres hacer? (última línea, extraes las coordenadas y 
pintas las etiquetas).

library(maps)
library(mapdata)
library(maptools)
library(sp)

### Creación de datos de muestra
## Poly base
mapSw<-map("worldHires","Sweden",
            xlim=c(12,16.7),ylim=c(55.3,57),col="gray90",
            fill=TRUE,mar = c(4.1, 4.1, par("mar")[3], 0.1))
# 
http://stackoverflow.com/questions/26062280/converting-a-map-object-to-a-spatialpolygon-object
IDs <- sapply(strsplit(mapSw$names, ":"), function(x) x[1])
mapSw <- map2SpatialPolygons(mapSw, IDs=IDs,
                            proj4string=CRS("+proj=longlat +datum=WGS84"))
## Puntos aleatorios
points<-spsample(mapSw,n=20,type="random")
points <- SpatialPointsDataFrame(points, data.frame(id=1:20))

### Plots
plot(mapSw)
plot(points, add=TRUE, pch=19, col="red")
# map.scale(12.4,57, ratio=FALSE,relwidth=0.2,cex=0.8)
box()
## Añadir las etiquetas a los puntos
text(coordinates(points)[,1], coordinates(points)[,2], points en data$id)

Un saludo,


El 19/03/16 a las 02:05, Matias Ledesma escribió:
> Hola a todos!
> Estoy incursionando en la elaboraci�n de mapas con R....y no he encuentro la manera de agregar etiquetas a las estaciones muestreadas (points). Alguien tiene idea de como hacerlo o me puede indicar donde buscarlo?
> library(maps)library(mapdata)
> map("worldHires","Sweden",xlim=c(12,16.7),ylim=c(55.3,57),col="gray90",fill=TRUE,mar = c(4.1, 4.1, par("mar")[3], 0.1))dat<-read.delim("koord.R.txt")points(dat$Long,dat$Lat,pch=19,col="red",cex=1)map.scale(12.4,57,ratio=FALSE,relwidth=0.2,cex=0.8)box()
> Desde ya muchas gracias!SaludosMatias
>
>   		 	   		
> 	[[alternative HTML version deleted]]
>
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es


	[[alternative HTML version deleted]]