Stampa della mappa dei colori con codici postali in R o Python

Question

Ho alcuni dati demografici e firmografici statunitensi.
Vorrei tracciare aree di codice postale in uno stato o una regione più piccola (ad es. Città). Ogni area verrebbe annotata per colore e / o testo specifico per quell'area. L'output sarebbe simile a http://maps.huge.info/ ma a) con testo annotato; b) output pdf; c) programmabile in R o Python.

Esistono pacchetti e codici che mi consentono di farlo?

Answer 1

Suppongo che tu voglia mappe statiche.

_{(fonte: eduardoleoni.com )}

1) Ottieni gli shapefile dei zip e < a href = "http://www.census.gov/geo/www/cob/st2000.html" rel = "nofollow noreferrer"> state confini a census.gov:

2) Utilizza la funzione plot.heat che ho pubblicato in questa domanda SO .

Ad esempio (presuppone che tu abbia gli shapefile maryland nella sottodirectory della mappa):

library(maptools)
##substitute your shapefiles here
state.map <- readShapeSpatial("maps/st24_d00.shp")
zip.map <- readShapeSpatial("maps/zt24_d00.shp")
## this is the variable we will be plotting
zip.map@data$noise <- rnorm(nrow(zip.map@data))
## put the lab point x y locations of the zip codes in the data frame for easy retrieval
labelpos <- data.frame(do.call(rbind, lapply(zip.map@polygons, function(x) x@labpt)))
names(labelpos) <- c("x","y")                        
zip.map@data <- data.frame(zip.map@data, labelpos)
## plot it
png(file="map.png")
## plot colors
plot.heat(zip.map,state.map,z="noise",breaks=c(-Inf,-2,-1,0,1,2,Inf))
## plot text
with(zip.map@data[sample(1:nrow(zip.map@data), 10),] , text(x,y,NAME))
dev.off()

Answer 2

Esistono molti modi per farlo in R (vedi vista spaziale ); molti di questi dipendono dalle " mappe " pacchetto .

• Dai un'occhiata a questo fantastico esempio delle elezioni USA 2004 . Si presenta così:

• Ecco un esempio leggermente brutto di un modello che utilizza le mappe " ; pacchetto con " reticolo " ;.

• Andrew Gelman ha realizzato delle trame molto belle come questa. Vedi, ad esempio, questo post sul blog sugli stati rossi / stati blu e questo seguito postare .

• Ecco un esempio molto semplice utilizzando " gmaps " ; pacchetto , che mostra una mappa degli arresti statali per arresti per 100.000 per omicidio:

  require(gmaps)
  data(USArrests)
  attach(USArrests)
  grid.newpage()
  grid.frame(name="map")
  grid.pack("map",USALevelPlot(states=rownames(USArrests),levels=Murder,col.fun=reds),height=unit(1,'null'))
  grid.pack("map",gradientLegendGrob(at=quantile(Murder),col.fun=reds),side="bottom",height=unit(.2,'npc'))
  detach(USArrests)
  

Answer 3

Qualcuno potrebbe avere qualcosa di più diretto per te, ma ho trovato 'Data Mashups in R' di O'Reilly molto interessante ... in parte, è una mappatura spaziale delle aste di preclusione a casa.

http://oreilly.com/catalog/9780596804770/

Answer 4

In Python, puoi usare gli shapefile dal censimento degli Stati Uniti insieme al pacchetto basemap . Ecco un esempio di compilazione degli stati in base alla popolazione.

Answer 5

Esiste una serie ricca e sofisticata di pacchetti in R per tracciare, fare analisi e altre funzioni relative a GIS. Un punto di partenza è la visualizzazione dell'attività CRAN su Dati spaziali : Questo è un mondo complesso e talvolta arcano e richiede del lavoro per capire.

Se stai cercando un'applicazione di mappatura gratuita e molto funzionale, posso suggerire:

MapWindow (mapwindow.com)

Answer 6

Daniel Levine presso TechCrunch Trends ha fatto cose interessanti con il pacchetto maps in R. Ha anche un codice disponibile sul suo sito.

Il suggerimento di Paul di esaminare Elaborazione - che Ben Fry usava per creare il codice zip - è anche buono, se sei pronto per imparare una nuova lingua (simile a Java).

Answer 7

A seconda dell'applicazione, potrebbe essere molto utile utilizzare qualcosa del genere:

http: // googlemapsmania .blogspot.com / 2006/07 / new-google-maps-us-CAP-mashups.html

Per mappare i tuoi dati. Se non era esattamente quello che volevi, puoi ottenere gli shapefile di codici postali grezzi da census.gov e farlo manualmente, il che è abbastanza doloroso.

Inoltre, se non l'hai visto, questo è un modo semplice per interagire con dati simili e potrebbe offrire alcuni suggerimenti:

http://benfry.com/zipdecode/

Answer 8

Dai un'occhiata a questo eccellente strumento di visualizzazione online di IBM http://manyeyes.alphaworks.ibm.com/manyeyes/

MODIFICA FYI, ManyEyes utilizza il Toolkit di visualizzazione di prefuse per alcuni dei suoi cioè. Anche se si tratta di un framework basato su Java, forniscono anche uno strumento Flash / ActionScript per il Web.