是否可以使用ggrough (https://xvrdm.github.io/ggrough/index.html https://xvrdm.github.io/ggrough/index.html) 对创建的形状进行着色geom_sf（首选）或可能geom_polygon？请参阅此问题了解先前的问题，该问题给出了我想要的情节的外观以及 Z.Lin 的随附答案，该答案修改了程序包以使其与当前版本兼容ggplot2: 无法复制此 ggplot2 图 https://stackoverflow.com/questions/63759571/unable-to-replicate-this-ggplot2-plot?noredirect=1#comment112795312_63759571.
这是使用创建的地图的 MWEgeom_sf我想使用（每个县）来遮蔽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 or geom_polygon;如果没有，添加它们是否容易？
***更新 这是另一个例子，摘自https://ggplot2.tidyverse.org/reference/ggsf.html 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
这会产生：
***更新结束
(Here is an example, created with ggrough, of how I want the shading of the counties to look:
)
这是一次失败的尝试（再次依赖于 Z.Lin 的答案中的代码：无法复制此 ggplot2 图 https://stackoverflow.com/questions/63759571/unable-to-replicate-this-ggplot2-plot?noredirect=1#comment112795312_63759571):
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是优选的。
library(magrittr)
library(ggplot2)
library(ggrough)
Replace parse_rough using trace
trace(ggrough:::parse_rough, edit=TRUE)
在弹出窗口中，粘贴此内容，以便parse_rough将使用parse_sf for GeomSf.
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("GeomArea", "GeomViolin", "GeomSmooth",
"Background")) {
rough_els <- append(rough_els, parse_areas(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))
}
if (geom %in% c("GeomSf")) {
rough_els <- append(rough_els, parse_sf(svg))
}
purrr::map(rough_els, ~purrr::list_modify(.x, geom = geom))
}
创建函数parse_sf.
parse_sf <- function (svg) {
shape <- "path"
keys <- NULL
ggrough:::parse_shape(svg, shape, keys) %>% {
purrr::map(.,
~purrr::list_modify(.x,
points = .x$d,
shape = "path"
))
}
}
创建所需的情节
nc <- sf::st_read(system.file("shape/nc.shp", package = "sf"), quiet = TRUE)
b <- ggplot(nc) +
geom_sf(aes(fill = AREA))
b
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)