使用不同的时间步以动画形式动画化两个图层

时间:2018-11-30 20:26:13

标签: r ggplot2 gganimate

我正在使用gganimate将动物跟踪数据叠加在猎物可用性上,并且两者都随着时间推移而动画化,以显示动物的运动随着猎物可用性的变化而变化。

跟踪数据集

组成
  • sumDur变量
  • 经度
  • 每小时间隔

猎物数据集

组成
  • 能源价值
  • 经度
  • 每周一次

数据框示例

    head(Tracks)
          sumDur      lon2      lat                 gmt
30641  0.0000000 -170.2978 57.10774 2004-08-18 05:05:00
29380  0.0000000 -170.3264 57.05684 2004-08-18 06:00:00
29381 10.0833333 -170.3565 57.00376 2004-08-18 07:00:00
29382  3.0833333 -170.3859 56.95370 2004-08-18 08:00:00
29383  0.8333333 -170.4147 56.90687 2004-08-18 09:00:00
29384  2.6666667 -170.4429 56.86301 2004-08-18 10:00:00

 head(Prey)
# A tibble: 6 x 4
  TotalEnergyM2   lat  lon2 gmt                
          <dbl> <dbl> <dbl> <dttm>             
1     0          50    160  2004-01-01 00:00:00
2     0          50    162. 2004-01-01 00:00:00
3     0.0000193  50    180. 2004-01-01 00:00:00
4     0          50.3 -151. 2004-01-01 00:00:00
5    38.3        53.9 -158  2004-01-01 00:00:00
6    31.4        53.9 -158. 2004-01-01 00:00:00

目标

我已经能够在静态环境之上成功地对跟踪数据进行动画处理,并且还可以对没有跟踪数据的猎物数据进行动画处理。

我的问题在于尝试同时创建两个动画-我得到了一个动画,但是即使我知道时间重叠,每个动画的时间安排也不正确。

我认为这是由于我的时间列(gmt)的帧长在每个数据帧中不同而引起的;在我的跟踪数据框中,时间步长是一小时,在猎物数据框中,时间步长是一周。不幸的是,我无法验证这一点,因为猎物数据帧太大,无法将其复制到一个小时的时间步长。我一直在尝试一些不同的转换(例如transition_state),但没有成功,也无法在这里或其他地方找到任何解决方案。

如果我没有提供足够的信息或格式已关闭,我们将不胜感激,并且会提前道歉,这是我的第一篇帖子。

示例代码

plot <- ggplot(NULL) +
  geom_raster(data=Prey, aes(x=lon2, y=lat2, fill=TotalEnergyM2) +
  geom_point(data=Tracks, aes(x=lon2, y=lat2, size=sumDur)) +
  scale_fill_viridis_c()+
  transition_time(gmt)+
  ease_aes("linear")+ 
  labs(title="{frame_time}", x="Longitude", y="Latitude")

animate(plot)

如果我注释掉geom_raster或geom_point以便获得单独的动画,或者geom_raster是静态的,则上面的代码可以完美地工作。我已经将两个带有假数据的数据文件的链接附加到了上面的代码中,并重现了我遇到的问题datafiles

1 个答案:

答案 0 :(得分:0)

在我的机器上,您的代码运行没有问题。如果您的datetime的格式正确地设置为日期transition_time(),则将填补数据点之间的所有空白:

plot <- ggplot() + 
    geom_raster(data = Prey, aes(x = lon, y = lat, fill = energy)) +
    scale_fill_viridis_c() +
    geom_point(data = Tracks, aes(x = lon, y = lat, size = sumdur)) +
    transition_time(datetime) +
    ease_aes("linear") +
    labs(title="{frame_time}", x="Longitude", y="Latitude")

animate(plot, nframes = 100)

1

数据

相关问题