我在 www.betydb.org 上看到了美国县级的交互式分区统计图。我想使用 R 重现类似的地图。我只想要地图和工具提示（不是不同缩放级别的所有图块，或者切换地图的能力）
该地图当前创建于ruby，弹出窗口（左下角）查询 MySQL 数据库。写它的程序员已经离开了，我对 Ruby 也不熟悉。
在这里，我将从 csv 文件开始。数据包括州和县名称以及州和县 FIPS。我想绘制Avg_yield.
mydata <- read.csv("https://www.betydb.org/miscanthus_county_avg_yield.csv")
colnames(mydata)
# [1] "OBJECTID" "Join_Count" "TARGET_FID" "COUNTY_NAME" "STATE_NAME" "STATE_FIPS"
# [7] "CNTY_FIPS" "FIPS" "Avg_lat" "Avg_lon" "Avg_yield"
我可以使用以下方法在州一级进行绘图googleVis package
library(googleVis)
p <- gvisGeoChart(data = mydata, locationvar="STATE_NAME", colorvar = 'Avg_yield',
options= list(region="US", displayMode="regions",
resolution="provinces"))
plot(p)
这提供了状态级着色。我的问题是，如何在县级（而不是州级）分辨率下获得类似颜色和工具提示的内容？
The gvisGeoChart帮助（在区域和分辨率下）和谷歌图表文档表明这可能是不可能的，但文档是如此广泛，以至于不清楚我在 R 中的其他选项是什么。
那么，有没有办法获得带有工具提示和县级着色的地图呢？
这是2013年的问题。我不确定是否leaflet当时包裹已经出来了。现在已经2017年年底了，你的任务是有可能完成的。如果您仍然需要执行类似的任务，我想为您留下以下内容。在本例中，数据集中缺少一些县。这些县存在于 USA 多边形数据中，但在mydata。所以我将这些县添加到mydata using setdiff() and bind_rows()。当您绘制传单地图时，您需要指定您的调色板。Avg_yield是一个连续变量。所以你用colorNumeric()。我留下了显示传单地图一部分的屏幕截图。
library(raster)
library(leaflet)
library(tidyverse)
# Get USA polygon data
USA <- getData("GADM", country = "usa", level = 2)
### Get data
mydata <- read.csv("https://www.betydb.org/miscanthus_county_avg_yield.csv",
stringsAsFactors = FALSE) %>%
dplyr::select(COUNTY_NAME, Avg_yield)
### Check counties that exist in USA, but not in mydata
### Create a dummy data frame and bind it with mydata
mydata <- data.frame(COUNTY_NAME = setdiff(USA$NAME_2, mydata$COUNTY_NAME),
Avg_yield = NA,
stringsAsFactors = FALSE) %>%
bind_rows(mydata)
### Create a color palette
mypal <- colorNumeric(palette = "viridis", domain = mydata$Avg_yield)
leaflet() %>%
addProviderTiles("OpenStreetMap.Mapnik") %>%
setView(lat = 39.8283, lng = -98.5795, zoom = 4) %>%
addPolygons(data = USA, stroke = FALSE, smoothFactor = 0.2, fillOpacity = 0.3,
fillColor = ~mypal(mydata$Avg_yield),
popup = paste("Region: ", USA$NAME_2, "<br>",
"Avg_yield: ", mydata$Avg_yield, "<br>")) %>%
addLegend(position = "bottomleft", pal = mypal, values = mydata$Avg_yield,
title = "Avg_yield",
opacity = 1)
