R语言gganimate疫情数据动态可视化

1047521767

R语言gganimate学习
所需加载包
library(av)
library(ggplot2)
library(gganimate)
library(tidyverse)
library(lubridate)
library(scales)
library(ggrepel)
2 o' P4 \! I; c% G. X! A' ilibrary(cowplot)
( N# j$ A0 i4 k; R数据


ps = ggplot(mydatan, aes(x=reorder(省份, 累计确诊),y=累计确诊, fill=省份,frame=Date)) +
$ o4 Z* j2 y$ s     geom_bar(stat= 'identity', position = 'dodge',show.legend = FALSE) +
6 a' T( E/ R( [- ?9 F    geom_text(aes(label=paste0(累计确诊)),col="black",hjust=-0.2)+
    #theme_bw()+
   #theme(legend.position="none") +
    theme(axis.text.x = element_text(size = 12,angle = 90, hjust = 0.2, vjust = 0.2),legend.position="none") +
" s; h* |: A6 T4 l! ?+ \$ ?& S. B; Y    theme(panel.background=element_rect(fill='transparent'))+
  theme(axis.text.y=element_text(angle=0,colour="black",size=12,hjust=1))+
9 e2 f+ G; z2 ?- S  theme(axis.text.x=element_text(angle=0,colour="white",size=2,hjust=1))+
   theme(panel.grid =element_blank()) +   ## 删去网格线
4 v5 H( z* c( r$ i0 L. W& V# p( Z  theme(axis.text = element_blank()) +   ## 删去所有刻度标签
! t0 X% S! q' z( P  theme(axis.ticks = element_blank()) +  ## 删去所有刻度线
# Here comes the gganimate specific bits
#labs(title = '日期:', x = '省份', y = '累计确诊病例') +
4 g! t3 a( c% d( T#annotate("text",x=0,y=40,label=C,parse=T)+
0 J3 d( A7 ^# q) K$ T8 Fcoord_flip()+
transition_manual(frames=as.Date(Date)) +
, R( e7 D3 @; W$ s#ggdark::dark_theme_bw() + #设置黑色主题
; L  W& l7 E/ k, {$ `labs(title = paste('日期:', '{current_frame}'),x = '', y ='各省累计确诊病例增长(除湖北省外)')+
) H& }& r' N6 `" i' n% h3 Uease_aes('linear')
$ W, y6 T& _2 {6 ]0 Hps
$ [" J" D  ^- [' N" k9 N$ R: W: ?结果展示


视频格式转化，加载BGM
& ~8 s/ c" `1 P#df <- animate(ps, renderer = av_renderer('animation.mp4'),
" \$ o3 K( ~: Q) K# H         #     width = 1280, height = 720, res = 100, fps = 10)#视频制作
+ n" A5 _9 C- h) e$ n7 j+ X# av_encode_video(df, 'output.mp4', framerate = 2,audio ="N.mp3")
% a% z7 M3 p5 m; M! \全国新冠状肺炎26天增长状况
  Y+ U% s/ t+ H/ P& w
! V/ f& ?% J9 P# A3 Spc<-ggplot(data=CNdata_s,aes(x=variable,y=value,fill=variable,frame=Date))+
9 ~$ Y- w9 V2 ^( V) a! k   geom_bar(stat= 'identity', position = 'dodge',show.legend = FALSE,width=0.7) +
   geom_text(aes(label=paste0(value)),col="black",hjust=-0.2)+
    theme(legend.position="none") +
  theme(panel.background=element_rect(fill='transparent'))+
  theme(axis.text.x=element_text(angle=0,colour="black",size=15,hjust=1))+
" E9 h2 v+ K& ]- {2 N  theme(axis.text.y=element_text(angle=0,colour="white",size=2,hjust=1))+
   theme(panel.grid =element_blank()) +   ## 删去网格线
- x$ ]9 H; z) c+ b5 N  theme(axis.text = element_blank()) +   ## 删去所有刻度标签
0 g4 ?: d' X$ {' v- {* T  theme(axis.ticks = element_blank()) +  ## 删去所有刻度线
: `  w0 S5 \$ U% c  #scale_x_continuous(limits = c(0,6))+
# Here comes the gganimate specific bits
: |1 j# _6 N, N" C#labs(title = '日期:', x = '省份', y = '累计确诊病例') +
#annotate("text",x=0,y=40,label=C,parse=T)+
    # coord_flip()+
# z4 K% }2 t- C* d! i     transition_manual(frames=as.Date(Date)) +
    # ggdark::dark_theme_bw() + #设置黑色主题
     labs(title = paste('日期:', '{current_frame}'),x = '全国新冠状肺炎增长', y ='')+
     ease_aes('linear')
pc
8 ]  K7 a  w8 a1 V; c
) t% \' j6 v) H/ o5 X7 N1 V
0 V& P$ L% r7 w' o8 s' I) @9 V. q动态图合并
+ {& g4 ]9 X( wlibrary(magick)
ps_gif <- animate(ps,,width = 720, height = 480)
pc_gif <- animate(pc, width = 360, height = 480)
ps_gifs <- image_read(ps_gif)
pc_gifs <- image_read(pc_gif)
new_gif <- image_append(c(pc_gifs[1], ps_gifs[1]))
for(i in 2:length(pc_gifs)){
3 H: ?$ l& t! R  combined <- image_append(c(pc_gifs, ps_gifs))
. [5 I+ u" G) N& |  new_gif <- c(new_gif, combined)
}
new_gif
结果展示
9 o# U8 n0 b& F- |

. H6 I2 a( k; E- A; D0 ]+ _% [2 }

sjlxdn

111111111111
5 t9 \, I( m# Z! Y5 [* h