sto indovinando c'è una soluzione semplice al problema Ho avuto, ma sto avendo qualche problema.Conversione di un oggetto "mappa" ad un oggetto "SpatialPolygon"
Sto cercando di convertire l'map
seguente oggetto:
require(maps)
usa <- map("state")
in un oggetto SpatialPolygon
utilizzando la funzione di map2SpatialPolygons
:
require(maptools)
usa.sp <- map2SpatialPolygons(usa, IDs=usa$names,proj4string=CRS("+proj=longlat"))
continuo a ricevere il seguente errore:
Error in map2SpatialPolygons(usa, IDs = usa$names, proj4string = CRS("+proj=longlat")) :
map and IDs differ in length
Dopo alcuni resea RCH, sembra che gli ID di avere una lunghezza 63 e l'oggetto map
ha lunghezza 169 dopo l'applicazione della funzione .NAmat2xyList(cbind(map$x, map$y))
(per la quale non riesco a trovare la sorgente).
Qualcuno ha qualche idea? Ecco la struttura dell'oggetto usa
mappa:
> str(usa)
List of 4
$ x : num [1:1705] -88.4 -88.1 -88 -87.9 -87.8 ...
$ y : num [1:1705] 30.4 30.4 30.8 30.6 30.3 ...
$ range: num [1:4] -124.7 -67 25.1 49.4
$ names: chr [1:63] "alabama" "arizona" "arkansas" "california" ...
- attr(*, "class")= chr "map"
Quando si utilizza la mappa specificando il database mondiale e si accede ai paesi in base al nome della propria regione, è possibile accedere alle regioni in modo più generale: germania <- mappa ("mondo", regioni = "Germania", fill = TRUE) quanto ho capito solo un insieme limitato di banche dati è definita nel modo in cui hai chiamato esso, la Germania per esempio non lo è. Tuttavia, la tua risposta mi ha salvato la giornata! – Exocom