R语言 使用不同的时间步长在gganimate中设置两个层的动画

vwkv1x7d  于 2023-03-10  发布在  其他
关注(0)|答案(1)|浏览(106)

我正在使用gganimate将动物跟踪数据叠加在猎物可用性之上,两者都通过动画显示动物的运动如何随着猎物可用性的变化而变化。
追踪数据集**由

  • sumDur变量
  • 经纬度
  • 每小时一次
    猎物数据集包括
  • 能量值
  • 经纬度
  • 每周一次
    ** Dataframe 示例**
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)的帧长度在每个 Dataframe 中不相同的事实引起的;在我的跟踪 Dataframe 中,时间步长是一小时,而在猎物 Dataframe 中,时间步长是一周。不幸的是,我无法验证这一点,因为猎物 Dataframe 太大,无法将其复制到每小时的时间步长。我尝试了一些不同的转换(例如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问题。

chy5wohz

chy5wohz1#

在我的机器上,您的代码运行没有问题。如果您的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)

数据

相关问题