2

我是Spatial data& cartogramlib 的新手，在尝试从以下位置重新创建绘图时遇到了一些问题：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.....我在使用sflib 进行了一些谷歌搜索后修复了它。

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)

问题：现在在这一步之后，我无法重新创建演示网站中的代码，因为我的数据中没有列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专栏？？？

网站中使用的代码：

data(wrld_simpl)
afr=wrld_simpl[wrld_simpl$REGION==2,]

afr_cartogram <- cartogram(afr, "POP2005", itermax=7)

# Transform these 2 objects in dataframe, plotable with ggplot2
afr_cartogram_df <- tidy(afr_cartogram) %>% left_join(. , afr_cartogram@data, by=c("id"="ISO3")) 
afr_df <- tidy(afr) %>% left_join(. , afr@data, by=c("id"="ISO3")) 
 
# And using the advices of chart #331 we can custom it to get a better result:
ggplot() +
  geom_polygon(data = afr_df, aes(fill = POP2005/1000000, x = long, y = lat, group = group) , size=0, alpha=0.9) +
  theme_void() +
  scale_fill_viridis(name="Population (M)", breaks=c(1,50,100, 140), guide = guide_legend( keyheight = unit(3, units = "mm"), keywidth=unit(12, units = "mm"), label.position = "bottom", title.position = 'top', nrow=1)) +
  labs( title = "Africa", subtitle="Population per country in 2005" ) +
  ylim(-35,35) +
  theme(
    text = element_text(color = "#22211d"), 
    plot.background = element_rect(fill = "#f5f5f4", color = NA), 
    panel.background = element_rect(fill = "#f5f5f4", color = NA), 
    legend.background = element_rect(fill = "#f5f5f4", color = NA),
    plot.title = element_text(size= 22, hjust=0.5, color = "#4e4d47", margin = margin(b = -0.1, t = 0.4, l = 2, unit = "cm")),
    plot.subtitle = element_text(size= 13, hjust=0.5, color = "#4e4d47", margin = margin(b = -0.1, t = 0.4, l = 2, unit = "cm")),
    legend.position = c(0.2, 0.26)
  ) +
  coord_map()
4

1 回答 1

1

这些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代码可以正常工作：

在此处输入图像描述

于 2021-05-17T13:53:53.217 回答