如何解决在 r 中创建制图图时出现问题
我是 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
代码工作正常:
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。