Python：使用 Basemap 在 netCDF 文件中绘制地图

如何解决Python：使用 Basemap 在 netCDF 文件中绘制地图

我正在使用 netCDF 文件绘制地图。但是，我是新手，进展很小，我所做的只是学习了一些教程。在地图上叠加颜色时我遇到了一些困难。这是我的代码（带有下载的 file）：

from mpl_toolkits.basemap import Basemap,cm
from netCDF4 import Dataset
import numpy as np
import matplotlib.pyplot as plt

data = Dataset('filename','r')


lat = data.variables['lat'][:]
lon = data.variables['lon'][:]
time = data.variables['time'][:]
lwe_thickness = data.variables['lwe_thickness'][:]
data.close()


mp = Basemap(projection = 'mill',llcrnrlon=lon.min(),urcrnrlon=lon.max(),llcrnrlat=lat.min(),urcrnrlat=lat.max(),resolution = 'c')

x,y = np.meshgrid(lon,lat)
x,y = mp(lon,lat)

mp.drawcoastlines()
mp.drawstates()
mp.drawcountries()
mp.drawmapboundary()

parallels = np.arange(0.,90,30.)
mp.drawparallels(parallels,labels=[1,0])

meridians = np.arange(-180.,180.,30.)
mp.drawmeridians(meridians,labels=[0,1])

cmesh = mp.pcolormesh(x,y,lwe_thickness,shading='flat',cmap=plt.cm.jet,latlon=True)
cbar = mp.colorbar(cmesh,location='right')

plt.show()