当我使用轮廓图在底图上绘制不连续的地理空间数据时,绘图未正确绘制可能是由于其插值。下面给出了数据、图形和示例脚本。
fig,ax = plt.subplots()
m = Basemap(projection='cyl',llcrnrlat=-50,urcrnrlat=50,llcrnrlon=0,urcrnrlon=360,ax=ax)
m.drawcoastlines()
m.contourf(lon,lat,data,latlon=True)
plt.show()
我还尝试使用 pcolormesh 而不是contourf,但由于数据中的 nan 值,它不起作用。
m.pcolormesh(lon,lat,data,ax=ax)
我也尝试了here提到的解决方案,但它给出了以下错误。
ndata = ma.masked_where(np.isnan(data),data)
lat = ma.masked_where(np.isnan(lat),lat)
lon = ma.masked_where(np.isnan(lon),lon)
cs = m.pcolormesh(lon,lat,ndata,ax=ax)
ValueError: x and y arguments to pcolormesh cannot have non-finite values or be of type numpy.ma.core.MaskedArray with masked values
请为我提供正确绘制地理空间数据的解决方案。