近年来,台风(飓风)在世界各地频繁发生,有些台风造成了严重的人身和财产损失。研究台风最基本的条件是弄清台风的路径和各阶段的气象特征。这里我打算用R做台风路径图,包括静态图和动态图。本文是一篇可重复的博文()。感兴趣的读者可以点击门户查看相关数据和代码。
导入相关的扩展包和数据
() () () ()库(gg map)库(maps)库(mapdata) #中文字体设置。# font _ files()show text _ auto(enable=true)font _ add(song ti,song ti . TTC)# import data typoon _ data-read . CSV(。/data/typhoon _ mangkhut.csv)
数据清理
我们需要对数据做进一步的处理,包括转换时间格式和重命名经度。
轨迹-台风_数据%%#转换时间格式突变(time=mdy_hm(time))%%#重命名经度重命名(long=lng)
画静态地图
根据浙江省水利厅的数据,作者绘制了2018年超强台风山竹的路线图。山竹是2018年太平洋台风季节命名的第22个热带气旋。9月初太平洋国际交流日附近形成,然后一路向西;9月16日下午1时登陆中国香港,并发出最高热带气旋警告信号;9月16日下午5时在广东江门市台山登陆,后在中国广西省失踪。山竹超级台风,造成很多地方“停课停工”。据广东省统计,死亡4人,发生广泛经济损失。下图可以清晰显示2018年超级台风山竹的路径图。
tracks _ box-make _ bbox(lon=tracks $ long,lat=tracks$lat,f=0.001)sq _ map-get _ map(location=tracks _ box,maptype=satellite,source=google,zoom=3)p-ggmap(sq_map)主题(text=element _ text(family=Songti))geom _ point(data=tracks,mapping=aes(x=long,y=lat,color=pressure))geom _ line(data=tracks,Mapping=AES (x=long,y=lat,
绘制动态台风路径图
在上述静态路径图的基础上,作者继续绘制动态台风路径图。动态图采用R中gganimate扩展包的功能,以台风发展时间为轴绘制台风路径动态图。在下面的副标题中,作者添加了动态时间显示。需要指出的是,图片中的字体还存在一些问题,与静态图片中的字体不一致;这个动态图中的中文字体需要进一步研究。
ani-ggmap(sq_map)主题(text=element _ text(family=STFangsong,size=18))geom _ point(data=tracks,mapping=aes(x=long,y=lat,Color=pressure))geom _ line(data=tracks,mapping=aes(x=long,y=lat,Color=pressure))# geom _ path(data=tracks,mapping=aes(x=long,y=lat,Color=pressure))scale _ Color _ continuous(name=central pressure(HPA),low=黄色,high
保存动态图表
上图是直接内置在网页中的动态图片。如果读者需要保存动态图片,可以使用下面的命令进行操作。
gif _ ani=animal(ani)anim _ save(mangkhut _ tracking . gif,animation=gif _ ani,path=。/材料/)
摘要
r是非常强大的数据分析和可视化工具,感兴趣的读者可以关注微信官方账号(wliu_2018)更精彩的内容。
最后,感谢上帝。
相关阅读
标签: #山竹台风最新的动态