在 R 或 Python 中绘制带有邮政编码的彩色地图

Question

我有一些美国人口统计和公司统计数据。
我想绘制一个州或较小地区（例如城市）的邮政编码区域。每个区域都将用特定于该区域的颜色和/或文本进行注释。输出将类似于 http://maps.huge.info/ 但 a) 带有注释文本； b) pdf 输出； c) 可使用 R 或 Python 编写脚本。

有没有任何软件包和代码可以让我做到这一点？

Answer 1

根据您的应用程序，一个很长的方法可能是使用类似这样的内容：

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

映射您的数据。如果这不是您想要的，您可以从 census.gov 获取原始邮政编码形状文件并手动执行，这非常痛苦。

另外，如果您还没有看到它，这是与类似数据交互的一种巧妙方式，并且可能会提供一些指示：

http://benfry.com/zipdecode/

Answer 2

看看 IBM 这个出色的在线可视化工具
http://manyeyes.alphaworks.ibm.com/manyeyes/

编辑< /strong> 仅供参考，ManyEyes 在其某些可视化中使用 Prefuse 可视化工具包。尽管它是一个基于 java 的框架，但他们也为 Web 提供了 Flash/ActionScript 工具。

Answer 3

在 R 中，有很多方法可以做到这一点（请参阅空间视图 ）；其中许多依赖于“maps”包。

• 查看这个美国 2004 年大选的精彩示例。它最终看起来像这样：

• 这是一个稍微丑陋的模型示例 使用带有“lattice”的“maps”包。

• Andrew Gelman 制作了一些非常漂亮的图，如下所示。例如，请参阅这篇关于红色州的博客文章/blue states 和此后续行动发布。

• 这是一个非常简单的示例使用“gmaps”包，显示了按州划分的每 10 万人因谋杀被捕的逮捕情况地图：

  <前><代码>要求（gmap）
  数据（美国逮捕）
  附上（美国逮捕）
  grid.newpage()
  grid.frame(名称=“地图”)
  grid.pack(“地图”,USALevelPlot(states=rownames(USArrests),levels=Murder,col.fun=reds),height=unit(1,'null'))
  grid.pack("地图",gradientLegendGrob(at=分位数(谋杀),col.fun=reds),side="底部",height=unit(.2,'npc'))
  分离（美国逮捕）
  

Answer 4

有人可能会为您提供更直接的信息，但我发现 O'Reilly 的“R 中的数据混搭”非常有趣……在某种程度上，它是房屋止赎拍卖的空间映射。

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

Answer 5

在 Python 中，您可以使用美国人口普查中的 shapefile 以及 basemap 包。以下是根据人口填写州的示例。

Answer 6

R中有一系列丰富而复杂的包来绘制、分析以及与GIS相关的其他功能。一个起点是 空间数据 上的 CRAN 任务视图：
这是一个复杂且有时神秘的世界，需要一些努力才能理解。

如果您正在寻找一个免费的、功能强大的地图应用程序，我可以建议您：

MapWindow (mapwindow.com)

Answer 7

TechCrunch Trends 的 Daniel Levine 使用 R 中的 maps 包做了很多好事。他的网站上也有可用的代码。

如果您准备学习一门（类似 Java 的）新语言，Paul 建议研究一下 Processing（Ben Fry 曾用它来制作 zipdecode），这也是一个不错的建议。

Answer 8

我假设您想要静态地图。

_{（来源：eduardoleoni.com）}

1）获取形状文件zip 边界和 state census.gov 边界：

2) 使用我在此 所以问题。

例如（假设您的地图子目录中有马里兰州形状文件）：

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()