gpt4 book ai didi

r - 使用 gganimate 时 ggplot geom_polygon 的奇怪行为

转载 作者:行者123 更新时间:2023-12-03 16:50:33 27 4
gpt4 key购买 nike

说我有这些数据。一个用于对州的县进行着色,另一个用于绘制点。

library(tidyverse)
library(gganimate)
devtools::install_github("UrbanInstitute/urbnmapr")
library(urbanmapr)

#1. counties dataset for the shading of the background
data(counties)
#keep only Texas counties
counties <- filter(counties, state_fips==48)

#2. dots dataset for dots over time
dots <- data.frame(group = c(rep(1,3), rep(2,3)),
lat = c(rep(32, 3), rep(33, 3)),
long = c(rep(-100, 3), rep(-99,3)),
year = c(1:3, 1:3))

然后我绘制它们(首先不使用 gganimate ):
ggplot() +
geom_polygon(data = counties, aes(long, lat, group = county_fips, fill = as.numeric(county_fips))) +
scale_fill_gradient(low = "navy", high = "lightskyblue3") +
geom_point(data = dots, aes(long, lat, group=interaction(long, lat), color=year),
size=6, show.legend = FALSE) +
theme_bw() +
scale_color_gradientn(colours = c("red", "yellow", "darkgreen")) +
coord_map() +
labs(subtitle = paste('Year: {frame_time}')) +
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.15, .15)) +
theme(legend.key.size = unit(2,"line"),legend.title=element_text(size=16),
legend.text=element_text(size=14)) +
labs(fill = "abc")

结果如下:

enter image description here

当我尝试使用 gganimate 为点设置动画时,问题就出现了。 :
map <- ggplot() +
geom_polygon(data = counties, aes(long, lat, group = county_fips, fill = as.numeric(county_fips))) +
scale_fill_gradient(low = "navy", high = "lightskyblue3") +
geom_point(data = dots, aes(long, lat, group=interaction(long, lat), color=year),
size=6, show.legend = FALSE) +
theme_bw() +
scale_color_gradientn(colours = c("red", "yellow", "darkgreen")) +
coord_map() +
labs(subtitle = paste('Year: {frame_time}')) +
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.15, .15)) +
theme(legend.key.size = unit(2,"line"),legend.title=element_text(size=16),
legend.text=element_text(size=14)) +
labs(fill = "abc") +
transition_time(year) +
shadow_mark(size=6)

anim_save("output/test.gif", map, end_pause=6, width = 800, height = 800, duration=8)

请注意,除了 transition_time 之外,它完全相同和 shadow_mark线。结果如下:

enter image description here

背景颜色非常不同。可能是颜色颠倒了之类的;我不确定。但是,在我使用更多点和不同阴影值的真实示例中,阴影与实际数据几乎没有相似之处。世界上发生了什么,我该如何解决?

最佳答案

我不确定问题到底是什么,但多边形可能有问题,因为它们往往违反“每个观察一行”的整洁数据规则。此外,作为一般原则,我建议使用 geom_sf()每当您处理地理空间数据时。以下是使用 geom_sf() 重现动画的代码.

理论上,我应该能够用 geom_sf() 绘制点。同样,但我遇到了另一个看起来像 gganimate 错误的问题,所以我不得不解决并提取点的投影坐标并用 geom_point() 绘图。 .

library(tidyverse)
library(gganimate)
library(sf)
#> Linking to GEOS 3.7.2, GDAL 2.4.2, PROJ 5.2.0

# 1. get shape files for texas
library(tidycensus)
options(tigris_use_cache = TRUE)

# download counties for Texas
counties <- get_acs(
state = 48, # omit to download all US counties
geography = "county", year = 2015, geometry = TRUE,
variables = "B19013_001" # need to specify some variable to download, here median income
)
#> Getting data from the 2011-2015 5-year ACS

# 2. dots dataset for dots over time
dots <- data.frame(
group = c(rep(1,3), rep(2,3)),
lat = c(rep(32, 3), rep(33, 3)),
long = c(rep(-100, 3), rep(-99,3)),
year = c(1:3, 1:3)
)

# add geometry column
dots_sf <- dots %>%
mutate(
geometry = st_sfc(
map2(lat, long, function(lat, long) st_point(x = c(long, lat))),
crs = "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"
) %>%
st_transform(crs = 3083) # CRS appropriate for Texas, https://epsg.io/3083
) %>%
st_as_sf()

# and extract transformed coords
dots_sf_transf <- dots_sf %>%
mutate(
long_transf = map_dbl(geometry, function(x) x[[1]]),
lat_transf = map_dbl(geometry, function(x) x[[2]])
)


# 3. static plot
p <- ggplot() +
geom_sf(data = counties, aes(fill = as.numeric(GEOID))) +
scale_fill_gradient(low = "navy", high = "lightskyblue3") +
geom_point(
data = dots_sf_transf,
aes(x = long_transf, y = lat_transf, group = group, color = year),
size = 6, show.legend = FALSE
) +
scale_color_gradientn(colours = c("red", "yellow", "darkgreen")) +
coord_sf(crs = 3083) +
labs(subtitle = "Year: {frame_time}", fill = "abc") +
theme_void() +
theme(legend.position = c(0.15, .15))
p



# animation
p + transition_time(year) + shadow_mark(size=12)



创建于 2019-11-03 由 reprex package (v0.3.0)

关于r - 使用 gganimate 时 ggplot geom_polygon 的奇怪行为,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/58568431/

27 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com