Wykreślanie mapy kolorów z kodami pocztowymi w R lub Python

Question

Mam kilka danych demograficznych i firmograficznych z USA.
Chciałbym wykreślić obszary kodu pocztowego w stanie lub w mniejszym regionie (np. Mieście). Każdy obszar będzie oznaczony adnotacją koloru i/lub tekstu specyficznego dla tego obszaru. Dane wyjściowe będą podobne do http://maps.huge.info/, ale a) z tekstem opatrzonym adnotacją; b) wydruk pdf; c) skryptowalne w języku R lub Python.

Czy istnieje pakiet i kod, który pozwala mi to zrobić?

Answer 1

jestem zakładając chcesz map statycznych.

alt text http://files.eduardoleoni.com/mapUS.png

1) Get Shapefiles granic zip i state granice na census.gov:

2) Użyj funkcji plot.heat ja pisał w tym SO question.

Na przykład (zakłada masz Shapefiles Maryland w podkatalogu mapie):

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

W zależności od aplikacji, długa droga wokół może być używać coś takiego:

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

mapowania danych. Jeśli nie było to dokładnie to, czego chciałeś, możesz uzyskać surowy shapefiles kodu z pliku census.gov i zrobić to ręcznie, co jest dość uciążliwe.

Ponadto, jeśli nie widziałeś, jest to świetny sposób na interakcję z podobnymi danymi i może zaoferować kilka wskazówek:

http://benfry.com/zipdecode/

Answer 3

Sprawdź to doskonałe narzędzie do wizualizacji online IBM http://manyeyes.alphaworks.ibm.com/manyeyes/

EDIT FYI, ManyEyes używa Prefuse visualization toolkit dla niektórych jej mianowicie. Mimo że jest to środowisko oparte na Javie, udostępnia również narzędzie Flash/ActionScript dla sieci.

Answer 4

Ktoś może mieć dla ciebie coś bardziej bezpośredniego, ale uznałem, że "Data Mashups in R" O'Reilly'ego jest bardzo interesujące ... w części, jest to mapowanie przestrzenne aukcji wykluczenia domu.

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

Answer 5

Jest bogaty i wyrafinowany seria pakietów w R inwestycyjnym, zrobić analizę i inne funkcje związane z GIS. Jednym z miejsc, w którym można zacząć, jest widok zadania CRAN pod numerem Spatial Data: Jest to skomplikowany i niekiedy tajemniczy świat, który wymaga trochę pracy, aby go zrozumieć.

Jeśli szukasz darmowego, bardzo funkcjonalnej aplikacji mapowania, mogę zasugerować:

MapWindow (mapwindow.com)

Answer 6

Daniel Levine at TechCrunch Trends zrobił ładne rzeczy z pakietu w R. maps Ma kod dostępne również na jego stronie.

Sugestia Paula dotycząca przetwarzania, którą Ben Fry użył do stworzenia zipdecodu, jest również dobra, jeśli chcesz nauczyć się nowego języka (podobnego do Java).

Answer 7

Jest wiele sposobów zrobienia tego w R (patrz: spatial view); wiele z tych depend on the "maps" package.

• Sprawdź to cool example of the US 2004 election.Kończy się patrzy tak: alt text http://www.ai.rug.nl/~hedderik/R/US2004/US04Election-PopGraded.png

• Oto nieco brzydki przykładem modelu that uses the "maps" package with "lattice".

• Andrew Gelman wykonane bardzo ładne działki w taki sposób. Zobacz na przykład: this blog post on red states/blue states i this follow up post.

• Oto bardzo prosty przykład using the "gmaps" package, na której przedstawiono mapę aresztowań przez państwowe dla aresztowań za 100,000 za morderstwo:

  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 8

w Pythonie, można użyć Shapefiles ze spisu USA wraz z pakietem basemap. Oto przykład filling in states według populacji.