То, что я хотел бы иметь, это сгенерировать график упаковки по кругу со значениями, похожими на показанные здесь:1004 *
Однако я хотел бы иметь возможность построить это в пределах границ области на карте, такой как область округа, используя листовки (или другие пакеты) в R.
Я пыталсязапустить код, предоставленный с сайта в их примере, чтобы сгенерировать их сюжет, но я не могу заставить его работать по какой-то причине.
Когда я пытаюсь сделать это самостоятельно в отображенной области, это выглядит примерно так, как показано ниже.Спасибо за помощь.
packages <-c("tidyverse","charlatan" "ggmap", "reshape2", "raster", "GISTools", "rgdal", "sp", "maps", "mapdata", "rgeos", "maptools", "rstan", "leaflet")
lapply(packages, require, character.only = TRUE)
us <- getData('GADM', country='USA', level=2)
texas <- subset(us,NAME_1=="Texas")
Dallas <- subset(texas, NAME_2=="Dallas")
leaflet(Dallas) %>%
addPolygons() %>%
addTiles()
bb <- bbox(Dallas)
minx <- bb[1]
miny <- bb[2]
maxx <- bb[3]
maxy <- bb[4]
set.seed(23)
LL1k <- data.frame(ch_position(n = 7, bbox = c(minx, miny, maxx, maxy)))
LL1kcol <- as.data.frame(t(LL1k))
colnames(LL1kcol) <- c("lon", "lat")
leaflet(Dallas) %>%
addPolygons() %>%
addTiles() %>%
addCircles(lat = LL1kcol$lat, lng = LL1kcol$lon, radius = iris$Sepal.Length*500, color = iris$Species)
