【发布时间】:2021-01-09 20:05:42
【问题描述】:
是否可以使用ggrough (https://xvrdm.github.io/ggrough/index.html) 对geom_sf(首选)或可能geom_polygon 创建的形状进行着色?请参阅此问题以获取先前的问题,该问题给出了我想到的情节的外观以及 Z.Lin 的随附答案,该答案修改了包以使其与ggplot2的当前版本兼容:Unable to replicate this ggplot2 plot。
这是使用 geom_sf 创建的地图的 MWE,我想使用 ggrough 对(每个单独的县)进行着色:
library(tidyverse)
library(magrittr)
library(ggplot2)
library(ggrough)
library(RColorBrewer)
library(tidycensus)
library(viridis)
#install.packages("devtools") # if you have not installed "devtools" package
#devtools::install_github("xvrdm/ggrough")
library(hrbrthemes)
#get nevada shapefile
counties <- get_acs(
geography = "county", year = 2018, geometry = TRUE,
variables = "B19013_001", keep_geo_vars=TRUE
) %>% filter(STATEFP=="32")
counties$GEOID <- as.integer(counties$GEOID)
#############
a <- ggplot() +
geom_sf(data = counties, aes(fill = estimate)) +
scale_fill_viridis(discrete=FALSE, name="", guide=FALSE) +
theme_bw() +
theme(legend.position = c(0.15, .15)) +
theme(plot.subtitle = element_text(hjust = 0.8, vjust=-10, size=30)) +
theme(panel.background = element_rect(fill = 'white')) +
theme(panel.grid = element_blank(),axis.title = element_blank(),
axis.text = element_blank(),axis.ticks = element_blank(),
panel.border = element_blank())+
theme(legend.position = c(0.25, .15), legend.key.size = unit(2,"line"),
legend.title=element_text(size=16),
legend.text=element_text(size=14),
legend.direction = "vertical",
legend.box = "horizontal") +
labs(caption = "")
a
这会产生以下内容:
如何使用ggrough 为这张地图的县着色,或者这不可能?请注意,我认为ggrough 可以处理geom_col, geom_bar, geom_tile, geom_geom_area, geom_ribbon, geom_violin, geom_point, geom_jitter, geom_dotplot, geom_line, and geom_smooth,但我不确定geom_sf 或geom_polygon;如果没有,添加它们是否容易?
***更新 这是另一个例子,取自https://ggplot2.tidyverse.org/reference/ggsf.html:
nc <- sf::st_read(system.file("shape/nc.shp", package = "sf"), quiet = TRUE)
b <- ggplot(nc) +
geom_sf(aes(fill = AREA))
b
这会产生:
***结束更新
(这是一个使用ggrough 创建的示例,说明我希望县的阴影看起来如何:
)
这是一次失败的尝试(再次依赖 Z.Lin 的回答代码:Unable to replicate this ggplot2 plot):
parse_polygons <- function (svg) {
shape <- "polygon" # was "polyline" in ggrough:::parse_areas
keys <- NULL
ggrough:::parse_shape(svg, shape, keys) %>% {
purrr::map(.,
~purrr::list_modify(.x,
points = stringr::str_squish(.x$points) %>%
{stringr::str_glue("M{.}Z")},
shape = "path"))
}
}
trace(ggrough:::parse_rough, edit = TRUE)
# paste the following function into the pop-up window
function (svg, geom) {
rough_els <- list()
if (geom %in% c("GeomCol", "GeomBar", "GeomTile", "Background")) {
rough_els <- append(rough_els, parse_rects(svg))
}
if (geom %in% c("GeomSmooth", "Background")) { # removed GeomArea / GeomViolin from here
rough_els <- append(rough_els, parse_areas(svg))
}
if (geom %in% c("GeomArea", "GeomRibbon", "GeomViolin")) { # new condition here
rough_els <- append(rough_els, parse_polygons(svg))
}
if (geom %in% c("GeomPoint", "GeomJitter", "GeomDotPlot", "Background")) {
rough_els <- append(rough_els, parse_circles(svg))
}
if (geom %in% c("GeomLine", "GeomSmooth", "Background")) {
rough_els <- append(rough_els, parse_lines(svg))
}
if (geom %in% c("Background")) {
rough_els <- append(rough_els, parse_texts(svg))
}
purrr::map(rough_els, ~purrr::list_modify(.x, geom = geom))
}
options <- list(GeomSf=list(fill_style="hachure",
angle_noise=0.5,
gap_noise=0.2,
gap=1.5,
fill_weight=1))
get_rough_chart(a, options)
这会产生错误消息:
Error in `*tmp*`[[i]] : subscript out of bounds
***更新
或者用第二个例子:
options <- list(GeomSf=list(fill_style="hachure",
angle_noise=0.5,
gap_noise=0.2,
gap=1.5,
fill_weight=1))
get_rough_chart(b, options)
同样的错误。
***结束更新。
还请注意,可以使用geom_polygon 创建地图,这也很有趣,尽管geom_sf 是首选。
【问题讨论】:
-
你能发布一个可重现的例子吗?人口普查 API 需要 API 密钥。我已经用ggplot2.tidyverse.org/reference/ggsf.html 的第一个示例尝试了
ggrough,它实际上不会引发错误,而是返回一个空图。 -
是的,我一定会这样做的。感谢您让我知道 Census API 问题 - 我没有考虑过。
-
我刚刚用第二个例子更新了这个问题。我仍然得到同样的错误。有趣的是,你的情节没有,尽管是空的。