如何在R中下载纽约州所有县的数据，以获得传单地图 #来自http://eric.clst.org/Stuff/USGeoJSON 和 # https://en.wikipedia.org/wiki/List_of_United_States_counties_and_county_equivalents 纽约州% addLegend（pal=pal，value=~log10（pop），不透明度=1.0， labFormat=labelFormat（转换=函数（x）轮（10^x）））

以上代码是从中复制的

我知道如何按照代码中的要求下载NY State County数据。（或者换句话说，如何生成NYCountries.geojson文件）

---

我在前两条评论中浏览了这两个网站，但未能从美国的全部数据中对纽约州的数据进行子集划分。

在下载了22 MB的json文件后，我这样做了，它似乎起到了作用

# From http://eric.clst.org/Stuff/USGeoJSON and
# https://en.wikipedia.org/wiki/List_of_United_States_counties_and_county_equivalents
nycounties <- geojsonio::geojson_read("json/nycounties.geojson",
  what = "sp")
# Or use the rgdal equivalent:
# nycounties <- rgdal::readOGR("json/nycounties.geojson", "OGRGeoJSON")

pal <- colorNumeric("viridis", NULL)

leaflet(nycounties) %>%
  addTiles() %>%
  addPolygons(stroke = FALSE, smoothFactor = 0.3, fillOpacity = 1,
    fillColor = ~pal(log10(pop)),
    label = ~paste0(county, ": ", formatC(pop, big.mark = ","))) %>%
  addLegend(pal = pal, values = ~log10(pop), opacity = 1.0,
    labFormat = labelFormat(transform = function(x) round(10^x)))

图书馆（传单）
xy名称（xy）
[1] “地理标识”“州”“县”“名称”“LSAD”“普查地区”
#从维基百科的县列表中，找到Genesse县，
#应该位于纽约州
>xy[grepl（“36037”，xy$GEO_ID），]$STATE
[1] 36
#纽约州应该是36号
纽约市%
addTiles（）%>%
addPolygons（）

---

你从哪里获得数据的？@JorgeLopez查看原始问题中的第一个链接：那么，你知道如何获得数据了吗？

library(leaflet)

xy <- geojsonio::geojson_read("gz_2010_us_050_00_500k.json", what = "sp")

> names(xy)
[1] "GEO_ID"     "STATE"      "COUNTY"     "NAME"       "LSAD"       "CENSUSAREA"

# from Wikipedia list of counties, find Genesse county, 
# which should be located in NY state
> xy[grepl("36037", xy$GEO_ID), ]$STATE
[1] 36

# NY state should be number 36

nyc <- xy[xy$STATE == 36, ]

leaflet(nyc) %>%
  addTiles() %>%
  addPolygons()