Question

说我有这些数据。一个用于遮蔽一个州的县，另一个用于绘制点。

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")

这是结果：

当我尝试使用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这几行外，它完全相同。结果如下：

背景颜色非常不同。可能是颜色反转或其他原因。我不确定。但是，在我的实际示例中，为阴影使用更多的点和不同的值，阴影与实际数据几乎没有相似之处。世界上正在发生什么，我该如何解决？

Answer 1

我不确定问题到底是什么，但是多边形可能会出现问题，因为多边形倾向于违反“每观察一行一行”的整齐数据规则。另外，作为一般原则，我建议在处理地理空间数据时使用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)

^{由reprex package（v0.3.0）于2019-11-03创建}

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

1 个答案: