2011-12-10 10 views

Respuesta

8

Existen varios mapas que podría usar, dependiendo de qué información/detalle/etc. que necesita, pero para esto el muy agradable wrld_simpl va a hacer muy bien:

library(maptools) 
library(geosphere) 

data(wrld_simpl) 

US_lat = wrld_simpl$LAT[wrld_simpl$NAME == 'United States'] 
US_lon = wrld_simpl$LON[wrld_simpl$NAME == 'United States'] 

SWE_lat = wrld_simpl$LAT[wrld_simpl$NAME == 'Sweden'] 
SWE_lon = wrld_simpl$LON[wrld_simpl$NAME == 'Sweden'] 

points = gcIntermediate(c(US_lon, US_lat), c(SWE_lon, SWE_lat), 100) 

dev.new(width=6, height=4) 
plot(wrld_simpl) 
lines(points, col='red') 

enter image description here

+0

Nice! Lo estoy haciendo con líneas de muchos países a Brasil. Para algunos países después de ~ 120 grados Este (como Australia, Nueva Zelanda, Japón ...) las líneas van hacia el este, se estrechan, van 360 grados hacia el otro lado de la pantalla y luego continúan su trayectoria circular. Cualquier forma de resolver eso? – Rodrigo

Cuestiones relacionadas