问题描述
我是 Spatial data 和 cartogram 库的新手，在尝试重新创建图时遇到一些问题：https://www.r-graph-gallery.com/a-smooth-transition-between-chloropleth-and-cartogram.html
库和数据
library(tidyverse)
library(maptools)
library(cartogram)
library(viridis)
library(sf)
data("wrld_simpl")
afr_cartogram = wrld_simpl[wrld_simpl$REGION==2,]
在此之后，我遇到了一些错误：比如 st_transform ..... 在使用 sf lib 进行了一些谷歌搜索后我修复了它。
afr_sf <- st_as_sf(afr_cartogram)
afr_sf_proj = st_transform(afr_sf,3857)
afr_plot <- cartogram::cartogram(afr_sf_proj,"POP2005",itermax =7)
问题：现在，在这一步之后，我无法像在演示网站中那样重新创建代码，因为我没有列 {{1} } 在我的数据中。
group从哪里可以获得 ggplot() +
geom_polygon(data = afr_plot,aes(fill = POP2005/1000000,x = LON,y = LAT,group = group),size=0,alpha=0.9) +
theme_void()
列？？？
网站中使用的代码：
group解决方法
group 列在这些行中生成
afr_cartogram_df <- tidy(afr_cartogram) %>%
left_join(afr_cartogram@data,by = ("id" = "ISO3"))
afr_df <- tidy(afr) %>%
left_join(afr@data,by = c("id" = "ISO3"))
通过包 tidy 中的 broom 函数，该函数未附加在您的代码中！
使用 broom 附加 library(broom) 或从其命名空间调用 tidy()，如下所示：broom::tidy(...)。
代码中的“数据部分”应如下所示：
data(wrld_simpl)
afr <- wrld_simpl[wrld_simpl$REGION==2,]
afr_cartogram <- wrld_simpl[wrld_simpl$REGION == 2,]
afr_sf <- st_as_sf(afr_cartogram)
afr_sf_proj <- st_transform(afr_sf,3857)
afr_plot <- cartogram_cont(afr_sf_proj,"POP2005",itermax =7)
afr_cartogram_df <- broom::tidy(afr_cartogram) %>%
left_join(afr_cartogram@data,by=c("id" = "ISO3"))
afr_df <- broom::tidy(afr) %>%
left_join(afr@data,by=c("id" = "ISO3"))
随后的 ggplot 代码工作正常：
