pcolormesh不绘制点,而是绘制数据点之间的线

时间:2013-04-11 11:37:15

标签: python matplotlib matplotlib-basemap

我需要将雨数据(来自卫星观测)的数据绘制到grib2文件的地图上。最后,我设法通过文本文件将数据加载到numpy数组中,并使用Basemap将其与图片坐标联系起来。但问题是Python不会从数据中放置彩色点,而是倾向于在数据字段中的相邻点之间绘制线条,因此情节看起来很丑陋。 我没有看到问题的根源。

我的代码片段是:

import numpy as np
import matplotlib
matplotlib.use('Agg')
from scipy import *
from pylab import *
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import cm

之后,我读取了所需的数据并生成3个形状大约为(100000,2)的numpy数组,其中包含纬度,经度,度数和每个数据点的值。我用这些命令可视化它:

def joonista(lats,lons,value,nimi,clevs,koloriit):

---------描述重塑的片段(未显示),产生数组“lats”,“lons”,“value”------------

map=Basemap(projection='aea',lat_1=30,lat_2=50,lat_0=45,lon_0=0,llcrnrlon=-30,llcrnrlat=20,urcrnrlon=80,urcrnrlat=53,resolution='l',)

x, y = map(lons, lats)
map.drawcoastlines(linewidth=0.17,color='0.7')
map.drawcountries(linewidth=0.17,color='0.7')
map.drawmeridians(np.arange(-50,60,5),linewidth=0.17,color='0.7',labels=[False,False,False,True])
map.drawparallels(np.arange(-25, 70, 5),linewidth=0.17,color='0.7',labels=[True,False,False,False])
varvid=mpl.colors.ListedColormap(koloriit)
norm = mpl.colors.BoundaryNorm(clevs,varvid.N)
cs = map.pcolormesh(x,y,value,cmap=varvid,norm=norm)
savefig(nimi,dpi=300)
plt.clf()

joonista(latA,lonA,valueA,'h05',[-1,0.00001,0.001,0.01,0.1,0.3,0.5,1,2,3,4,5,6,7,8,9,10,11,12,13],['k','c','#a0fff9','#00b354','#69b300','#97ff03','#C2524D','#FF7500','#b3a900','#fff551','#515bff','#45adff','#da000d','#ff2a36','#ffa0a5','#f003ff','#f778ff','0.5','0.75'])
joonista(latB,lonB,valueB,'h04',[-1,0.0000000000001,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18],['k','c','#a0fff9','#00b354','#69b300','#97ff03','#C2524D','#FF7500','#b3a900','#fff551','#515bff','#45adff','#da000d','#ff2a36','#ffa0a5','#f003ff','#f778ff','0.5','0.75'])

以下是示例图片:

Example

如果我被告知如何解决这个问题,我将不胜感激。

的Aleksei

1 个答案:

答案 0 :(得分:1)

使用Joe Kington建议,我替换了命令

cs=map.pcolormesh(x,y,value,cmap=varvid,norm=norm)

按命令

cs=plt.scatter(x,y,c=value,s=0.6, edgecolors='none',marker=',',cmap=varvid,norm=norm)

可以很好地显示降水分布。

感谢您的帮助!