matplotlib删除地图投影上及等值线及风场.docVIP

matplotlib删除地图投影上的等值线及风场 【前言】最近在编写一个气象应用程序，用来显示某一时刻某一地区的气温等值线和风场，程序主要用到了第三方库matplotlib及Basemap。在编写的过程中发现，如果不进行擦除操作直接绘制新的等值线或风场，新的等值线（风场）会与原来的等值线（风场）叠加在一起，而绘制的等值线及风场没有单独的remove方法，所以如果想要擦除已经绘制的等值线就要将地图重新投影一遍，如果地图投影精度高一点，整个投影过程就会特别漫长。通过对等值线及风场的返回结果进行研究，我找到了一个不必重新投影地图就可将等值线及风场擦除的方法。 一、matplotlib及Basemap matplotlib是Python常用的数据绘制包。它基于numpy的数组运算功能，可以轻易的画出各种统计图形，如散点图，条行图，饼图，等值线图等。Basemap是Matplotlib的一个子包，负责地图绘制。在数据可视化过程中，我们可以将数据在地图上画出来。 利用matplotlib及Basemap画图的基本步骤是： 创建一个figure实例 在figure里创建Axes容器实例 在Axes容器内创建Basemap实例进行地图投影 调用Basemap实例的contour及barbs方法进行在地图上绘图 二、不进行擦除操作直接绘制新的等值线或风场效果演示 为了方便演示，数据为我自己手中的数据，其中values,x1,y1为等值线的数据及坐标，huvalues,hvvalues,x2,y2为风场的数据及坐标，数据的格式及获得方法就不做过多解释 1. 绘制地图投影 import matplotlib.pyplot as plt from mpl_toolkits.basemap import Basemap fig = plt.figure() ax = fig.add_axes([0,0,1,1]) m = Basemap(projection= laea, lon_0= 90, lat_0= 40, width= \ height= 8000000,resolution=l,ax=ax) m.drawcoastlines(color=tan) plt.show() 我们将得到下图演示的结果 2. 绘制等值线 import matplotlib.pyplot as plt from mpl_toolkits.basemap import Basemap fig = plt.figure() ax = fig.add_axes([0,0,1,1]) m = Basemap(projection= laea, lon_0= 90, lat_0= 40, width= \ height= 8000000,resolution=l,ax=ax) m.drawcoastlines(color=tan) c= m.contour(x1,y1,values,15,linewidths=1.5) plt.show() 得到结果如下： 3. 继续绘制风场 import matplotlib.pyplot as plt from mpl_toolkits.basemap import Basemap fig = plt.figure() ax = fig.add_axes([0,0,1,1]) m = Basemap(projection= laea, lon_0= 90, lat_0= 40, width= \ height= 8000000,resolution=l,ax=ax) m.drawcoastlines(color=tan) c= m.contour(x1,y1,values,15,linewidths=1.5) b= m.barbs(x2,y2,huvalues*2.5,hvvalues*2.5) plt.show() 得到结果如下： 由此我们可以看出，如果在绘制新的等值线或者风场前不进行擦除操作，所有的图像都会叠加在一起 三、通过重新进行地图投影进行擦除操作 1. 绘制地图投影 import matplotlib.pyplot as plt from mpl_toolkits.basemap import Basemap fig = plt.figure() ax = fig.add_axes([0,0,1,1]) m = Basemap(projection= laea, lon_0= 90, lat_0= 40, width= \ height= 8000