python2.7 读取nc文件绘制填色等值线图和风杆图

在 Python 2.7 中，我们可以使用 netCDF4 库来读取 nc 文件，使用 Basemap 库来绘制地图和等值线图，使用 matplotlib 库来绘制填色图和风杆图。

下面是一个简单的示例代码，其中假设 nc 文件中包含了温度和风场数据，我们需要将其绘制在地图上：

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

# 读取 nc 文件
nc_file = 'data.nc'
nc_data = netCDF4.Dataset(nc_file)

# 读取温度和风场数据
temp = nc_data.variables['temp'][:]
u_wind = nc_data.variables['u_wind'][:]
v_wind = nc_data.variables['v_wind'][:]

# 读取经纬度和时间信息
lon = nc_data.variables['lon'][:]
lat = nc_data.variables['lat'][:]
time = nc_data.variables['time'][:]
time_units = nc_data.variables['time'].units

# 绘制地图和等值线图
m = Basemap(projection='mill', llcrnrlon=lon.min(), llcrnrlat=lat.min(),
            urcrnrlon=lon.max(), urcrnrlat=lat.max(), resolution='l')
lon, lat = np.meshgrid(lon, lat)
x, y = m(lon, lat)

plt.figure(figsize=(10, 8))
cs = m.contourf(x, y, temp[0, :, :], cmap='coolwarm')
m.colorbar(cs)
m.drawcoastlines()
m.drawcountries()
m.drawstates()

# 绘制风杆图
skip = (slice(None, None, 3), slice(None, None, 3))
m.barbs(x[skip], y[skip], u_wind[0, skip], v_wind[0, skip], length=5)

plt.title('Temperature and Wind on {}'.format(time[0]))
plt.show()

在这个示例中，我们首先使用 netCDF4 库读取了 nc 文件中的数据，然后使用 Basemap 库绘制了地图和等值线图，使用 matplotlib 库绘制了填色图和风杆图。需要注意的是，我们需要根据 nc 文件中的经纬度信息来设置地图的范围和投影方式，同时需要根据温度和风场数据的形状来设置等值线图和风杆图的位置和大小。

希望这个示例代码能够对你有所帮助！