利用basemap做潮位变化gif图

葫芦爷儿

为了展示fvcom计算的潮位数据，用basemap来制作gif。做图思路：
读取nc数据 ->设置basemap底图 ->利用matplotlib的ArtistAnimation，循环每一时刻的数据（插值数据 ） -> 保存gif

用到这些东西的原因：
1 fvcom为非结构网格，于是利用matplotlib中的tri工具来插值到结构网格中
2 matplotlib 做动图有三种思路：
   a) 是animation.ArtistAnimation, 将每一张图片生成，存入列表，然后制作动图
   b) 是animation.FuncAnimation,创建函数，刷新数据，然后制作动图
    c) 是控制屏幕刷新时间（不推荐）
注意 linux、mac下可以直接从conda中装imagemagick
       windows要麻烦一些，但是可以通过ffmpe来实现.(下面有参考网址)
3 basemap设置
4 在集群下运行python需要
import matplotlib
matplotlib.use("agg")


下面参考到的网址

• 在集群中使用遇到问题
  

• https://www.cnblogs.com/xionggeclub/p/6146484.html 集群中的字体问题
• https://blog.csdn.net/gdkyxy2013/article/details/79585922 linux下的问题
  

• 调整图像的内容
  

• https://matplotlib.org/gallery/color/colormap_reference.html colormap
• https://stackoverflow.com/questions/3279560/reverse-colormap-in-matplotlib 调换
• https://blog.csdn.net/weixin_41789707/article/details/81035997 不想显示坐标轴
• https://stackoverflow.com/questions/18363987/basemap-how-to-remove-actual-lat-lon-lines-while-keeping-the-ticks-on-the-axis 网格线
• https://stackoverflow.com/questions/14280312/world-map-without-rivers-with-matplotlib-basemap 去除湖泊
• https://blog.csdn.net/pianzang5201/article/details/76070268  加logo
• https://stackoverflow.com/questions/29548704/matplotlib-pyplot-logo-image-in-plot 加logo
• https://www.cnblogs.com/chuangshaogreat/p/7498016.html python format
• https://www.114la.com/other/rgb.htm Colormap  配色参考
• https://matplotlib.org/gallery/images_contours_and_fields/irregulardatagrid.html#sphx-glr-gallery-images-contours-and-fields-irregulardatagrid-py contour等值线参考
• http://www.pythonware.com/products/pil/   这个只有python2.7的内容，建议不要用，
• https://www.bdmweather.com/2018/04/python-m-arcgisimage-basemap-options/  设置一些高清陆地影像
• https://anaconda.org/conda-forge/ffmpeg         windows下安装ffmpe
• https://stackoverflow.com/questions/41699931/how-to-install-imagemagick-with-anaconda-on-windows-64bits    保存动画为GIF和MP4
  

• 调整logo参考的内容
  

• https://blog.csdn.net/weixin_34613450/article/details/80678522 指定像素
• https://matplotlib.org/tutorials/intermediate/imshow_extent.html?highlight=extent 调整extent
  
• http://epanf.scsio.ac.cn/storm/storm.php 参考图片来源
  
  

werewolf

楼主fvcom有没耦合波浪？预报波高，波向等。

werewolf

有效波高。和观测还有一些误差。

方自之

谢谢楼主，当初也是为fvcom，mike21的非结构网格头疼的很久

流羽寻

这个可以用于洪水预报吗

葫芦爷儿

werewolf 发表于 2018-9-3 09:28
有效波高。和观测还有一些误差。

在学习swan，还没有耦合呢。用的holland特征风场驱动的

葫芦爷儿

流羽寻 发表于 2018-9-5 10:54
这个可以用于洪水预报吗

好像水利里用mike21多一些？fvcom非结构网格，也能很好的刻画河流，不过里面的洪水的水动力过程我不懂，你可以看看http://fvcom.smast.umassd.edu 了解一下

葫芦爷儿

关于画logo的补充：
1 之前直接在 fig, ax = plt.subplots()之后调用了ax.imshow()
2 今天，发现在使用cartopy画时，这样图片并不会显示
3 解决办法是，在set_extent()后调用imshow()
这是在git上找到的 https://github.com/SciTools/cartopy/issues/643

wenbinz

兄弟，用python画过单点的风场图吗？就是利用micaps的11类数据绘制风场图？

yijianmingliu

感谢分享，我下载来使用