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

【前言】最近在编写一个气象应用程序，用来显示某一时刻某一地区的气温等值线和风场，程序主要用到了第三方库matplotlib及Basemap。在编写的过程当中发现，若是不进行擦除操做直接绘制新的等值线或风场，新的等值线（风场）会与原来的等值线（风场）叠加在一块儿，而绘制的等值线及风场没有单独的remove方法，因此若是想要擦除已经绘制的等值线就要将地图从新投影一遍，若是地图投影精度高一点，整个投影过程就会特别漫长。经过对等值线及风场的返回结果进行研究，我找到了一个没必要从新投影地图就可将等值线及风场擦除的方法。

1、matplotlib及Basemap

matplotlib是Python经常使用的数据绘制包。它基于numpy的数组运算功能，能够轻易的画出各类统计图形，如散点图，条行图，饼图，等值线图等。
Basemap是Matplotlib的一个子包，负责地图绘制。在数据可视化过程当中，咱们能够将数据在地图上画出来。

利用matplotlib及Basemap画图的基本步骤是：

1. 建立一个figure实例
2. 在figure里建立Axes容器实例
3. 在Axes容器内建立Basemap实例进行地图投影
4. 调用Basemap实例的contour及barbs方法进行在地图上绘图

2、不进行擦除操做直接绘制新的等值线或风场效果演示

为了方便演示，数据为我本身手中的数据，其中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= 11000000, \
            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= 11000000, \
            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= 11000000, \
            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()

获得结果以下：

由此咱们能够看出，若是在绘制新的等值线或者风场前不进行擦除操做，全部的图像都会叠加在一块儿

3、经过从新进行地图投影进行擦除操做

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= 11000000, \
            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= 11000000, \
            height= 8000000,resolution='l',ax=ax)
m.drawcoastlines(color='tan')
c= m.contour(x1,y1,values,15,linewidths=1.5)
plt.show()

获得结果以下：

3. 擦除axes从新进行地图投影

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= 11000000, \
            height= 8000000,resolution='l',ax=ax)
m.drawcoastlines(color='tan')
x1,y1=m(lon,lat)
c= m.contour(x1,y1,values,15,linewidths=1.5)
ax.clear()
m = Basemap(projection= 'laea', lon_0= 90, lat_0= 40, width= 11000000, \
            height= 8000000,resolution='l',ax=ax)
m.drawcoastlines(color='tan')
plt.show()

获得下图结果：

4. 绘制风场

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= 11000000, \
            height= 8000000,resolution='l',ax=ax)
m.drawcoastlines(color='tan')
c= m.contour(x1,y1,values,15,linewidths=1.5)
ax.clear()
m = Basemap(projection= 'laea', lon_0= 90, lat_0= 40, width= 11000000, \
            height= 8000000,resolution='l',ax=ax)
m.drawcoastlines(color='tan')
x2,y2= m(hulon,hulat)
b= m.barbs(x2,y2,huvalues*2.5,hvvalues*2.5)
plt.show()

获得结果以下：

此方法虽然能够达到预期的效果，可是咱们的地图投影并未作过改变，从新进行地图投影毫无心义，并且还会占用系统没必要要的资源

4、对等值线实例及风场实例进行研究

1. 等值线

经过dir()命令咱们能够查看建立的等值线实例c的属性方法

print(dir(c))
['_A', '__class__', '__delattr__', '__dict__', '__dir__', '__doc__', '__eq__', '__format__', '__ge__', '__getattribute__', 
'__getstate__', '__gt__', '__hash__', '__init__', '__le__', '__lt__', '__module__', '__ne__', '__new__', '__reduce__', 
'__reduce_ex__', '__repr__', '__setattr__', '__sizeof__', '__str__', '__subclasshook__', '__weakref__', '_add_label',
 '_auto', '_autolev', '_check_xyz', '_contour_args', '_contour_generator', '_contour_level_args', '_corner_mask', 
'_get_allsegs_and_allkinds', '_get_label_clabeltext', '_get_label_text', '_get_lowers_and_uppers', '_initialize_x_y', 
'_levels', '_make_paths', '_process_args', '_process_colors', '_process_levels', '_process_linestyles', '_process_linewidths', 
'_transform', 'add_checker', 'add_label', 'add_label_clabeltext', 'add_label_near', 'allkinds', 'allsegs', 'alpha', 'antialiased', 
'autoscale', 'autoscale_None', 'ax', 'calc_label_rot_and_inline', 'callbacksSM', 'changed', 'check_update', 'clabel', 'cmap',
'collections', 'colorbar', 'colors', 'contour_doc', 'cvalues', 'extend', 'extent', 'filled', 'find_nearest_contour', 'get_alpha', 
'get_array', 'get_clim', 'get_cmap', 'get_label_coords', 'get_label_width', 'get_real_label_width', 'get_text', 'get_transform', 
'hatches', 'labelCValues', 'labelTexts', 'labels', 'layers', 'legend_elements', 'levels', 'linestyles', 'linewidths', 'locate_label',
 'locator', 'logscale', 'monochrome', 'nchunk', 'norm', 'origin', 'pop_label', 'print_label', 'set_alpha', 'set_array', 'set_clim',
 'set_cmap', 'set_label_props', 'set_norm', 'tcolors', 'tlinewidths', 'to_rgba', 'too_close', 'update_dict', 'vmax', 'vmin', 
'zmax', 'zmin']

经过查看，发如今c的属性方法中并无remove()方法，可是咱们发现有一个属性是collections，这里的collections若是与Python的集合类相同，它就应该有remove()方法进删除，因此继续经过dir()命令进行查看

print(dir(c.collections)
['__add__', '__class__', '__contains__', '__delattr__', '__delitem__', '__dict__', '__dir__', '__doc__', '__eq__', 
'__format__', '__ge__', '__getattribute__', '__getitem__', '__getstate__', '__gt__', '__hash__', '__iadd__',
 '__imul__', '__init__', '__iter__', '__le__', '__len__', '__lt__', '__module__', '__mul__', '__ne__', '__new__',
 '__reduce__', '__reduce_ex__', '__repr__', '__reversed__', '__rmul__', '__setattr__', '__setitem__',
 '__setstate__', '__sizeof__', '__str__', '__subclasshook__', '__weakref__', 'append', 'clear', 'copy',
 'count', 'extend', 'index', 'insert', 'pop', 'remove', 'reverse', 'sort', 'type']

经过查看，发现c.collections有remove()方法，因此尝试用c.collections.remove()进行删除，因为remove()方法一次只能删除一个元素，因此要用for循环来进行删除。

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= 11000000, \
            height= 8000000,resolution='l',ax=ax)
m.drawcoastlines(color='tan')
c= m.contour(x1,y1,values,15,linewidths=1.5)
for i in c.collections:
       i.remove()
plt.show()

获得下图结果：

因而可知，此方法可行

2. 风场

经过dir()命令咱们能够查看建立的等值线实例b的属性方法

print(dir(b))
['__add__', '__class__', '__contains__', '__delattr__', '__dir__', '__doc__', '__eq__', '__format__',
 '__ge__', '__getattribute__', '__getitem__', '__getnewargs__', '__gt__', '__hash__', '__init__',
 '__iter__', '__le__', '__len__', '__lt__', '__mul__', '__ne__', '__new__', '__reduce__', '__reduce_ex__', 
'__repr__', '__rmul__', '__setattr__', '__sizeof__', '__str__', '__subclasshook__', 'count', 'index']

能够发现风场的实例b没有remove方法，也没有collections属性，因此决定先经过type()方法看看风场实例是什么类型的

print(type(b))
<class 'tuple'>

能够发现风场实例b是一个元组，因此继续看它的元素有什么属性

print(dir(b[0]))
['_A', '__class__', '__delattr__', '__dict__', '__dir__', '__doc__', '__eq__', '__format__', '__ge__', '__getattribute__',
 '__getstate__', '__gt__', '__hash__', '__init__', '__le__', '__lt__', '__module__', '__ne__', '__new__', '__reduce__', 
'__reduce_ex__', '__repr__', '__setattr__', '__sizeof__', '__str__', '__subclasshook__', '__weakref__', '_agg_filter', 
'_alpha', '_animated', '_antialiaseds', '_axes', '_clipon', '_clippath', '_contains', '_edgecolors', 
'_edgecolors_original', '_facecolors', '_facecolors_original', '_factor', '_find_tails', '_get_bool',
 '_get_value', '_gid', '_hatch', '_is_filled', '_is_stroked', '_label', '_length', '_linestyles', '_linewidths',
 '_make_barbs', '_mouseover', '_offset_position', '_offsets', '_oid', '_path_effects', '_paths', '_picker',
 '_pickradius', '_pivot', '_prepare_points', '_propobservers', '_rasterized', '_remove_method',
 '_set_gc_clip', '_sizes', '_sketch', '_snap', '_stale', '_transOffset', '_transform', '_transformSet',
 '_transforms', '_uniform_offsets', '_url', '_urls', '_visible', 'add_callback', 'add_checker', 'aname',
 'autoscale', 'autoscale_None', 'axes', 'barb_increments', 'barbs_doc', 'callbacksSM', 'changed', 
'check_update', 'clipbox', 'cmap', 'colorbar', 'contains', 'convert_xunits', 'convert_yunits', 'draw', 
'eventson', 'figure', 'fill_empty', 'findobj', 'flip', 'format_cursor_data', 'get_agg_filter', 'get_alpha', 
'get_animated', 'get_array', 'get_axes', 'get_children', 'get_clim', 'get_clip_box', 'get_clip_on',
 'get_clip_path', 'get_cmap', 'get_contains', 'get_cursor_data', 'get_dashes', 'get_datalim', 
'get_edgecolor', 'get_edgecolors', 'get_facecolor', 'get_facecolors', 'get_figure', 'get_fill', 'get_gid', 
'get_hatch', 'get_label', 'get_linestyle', 'get_linestyles', 'get_linewidth', 'get_linewidths',
 'get_offset_position', 'get_offset_transform', 'get_offsets', 'get_path_effects', 'get_paths', 'get_picker',
 'get_pickradius', 'get_rasterized', 'get_sizes', 'get_sketch_params', 'get_snap', 'get_transform', 
'get_transformed_clip_path_and_affine', 'get_transforms', 'get_url', 'get_urls', 'get_visible', 
'get_window_extent', 'get_zorder', 'have_units', 'hitlist', 'is_figure_set', 'is_transform_set', 
'mouseover', 'norm', 'pchanged', 'pick', 'pickable', 'properties', 'remove', 'remove_callback',
 'rounding', 'set', 'set_UVC', 'set_agg_filter', 'set_alpha', 'set_animated', 'set_antialiased', 'set_antialiaseds', 
'set_array', 'set_axes', 'set_clim', 'set_clip_box', 'set_clip_on', 'set_clip_path', 'set_cmap', 'set_color',
 'set_contains', 'set_dashes', 'set_edgecolor', 'set_edgecolors', 'set_facecolor', 'set_facecolors', 'set_figure', 
'set_gid', 'set_hatch', 'set_label', 'set_linestyle', 'set_linestyles', 'set_linewidth', 'set_linewidths', 'set_lw',
 'set_norm', 'set_offset_position', 'set_offsets', 'set_path_effects', 'set_paths', 'set_picker', 'set_pickradius', 
'set_rasterized', 'set_sizes', 'set_sketch_params', 'set_snap', 'set_transform', 'set_url', 'set_urls', 'set_verts',
 'set_verts_and_codes', 'set_visible', 'set_zorder', 'sizes', 'stale', 'stale_callback', 'to_rgba', 'u', 'update',
 'update_dict', 'update_from', 'update_scalarmappable', 'v', 'x', 'y', 'zorder']

能够发现风场实例元组的元素本身就有remove属性，因此经过remove方法尝试一下可不能够进行风场的删除

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= 11000000, \
            height= 8000000,resolution='l',ax=ax)
m.drawcoastlines(color='tan')
b=m.barbs(x2,y2,huvalues*2.5,hvvalues*2.5)
for i in b:
    i.remove()
plt.show()

获得下图结果：

因而可知，此方法可行

5、总结

对于等值线，能够经过等值线实例的collection属性的remove方法进行删除已经绘制等值线；
对于风场，能够经过风场实例元组中个元素的remove方法进行删除已绘制风场；
对于地图上其余绘图方法，也能够经过本文的方法进行一步步的尝试。

---

本文版权归做者和博客园共有，欢迎转载，但未经做者赞成必须保留此段声明，且在文章页面明显位置给出原文链接，不然保留追究法律责任的权利。

想观看Matplotlib教学视频，了解更多Matplotlib实用技巧可关注

微信公众帐号： MatplotlibClass

今日头条号：Matplotlib小讲堂