绘制风场有压力！

po_po1

目前有两个风场的nc文件，一个是uwnd.nc，一个是vwnd.nc，分别表示两个方向的风场,每个文件都有各自的经纬度坐标（lon/lat）。可是绘图遇到困境了，函数quiver(lon,lat,u,v)没法写

po_po1

po_po1 发表于 2014-2-24 09:51

                               
登录/注册后可看大图

风场绘制出来了，虽然没有添加设置，有点粗糙：

顺便把程序传上来：
from netCDF4 import Dataset
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
u=Dataset(r'C:/Users/Administrator/Desktop/uwnd.2010.nc')
v=Dataset(r'C:/Users/Administrator/Desktop/vwnd.2010.nc')
uwnd=u.variables['uwnd'][0][0][:]
vwnd=v.variables['vwnd'][0][0][:]
m=Basemap(projection='cyl',llcrnrlat=0,urcrnrlat=70,llcrnrlon=50,urcrnrlon=150)
m.drawcountries()
m.drawcoastlines()
lat=u.variables['lat'][:]
lon=u.variables['lon'][:]
index1=np.logical_and(lon>=50,lon<=150);index2=np.logical_and(lat>=0,lat<=70)
lons=lon[index1];lats=lat[index2]
u1=uwnd[index2,:]
u2=u1[:,index1]
v1=vwnd[index2,:]
v2=v1[:,index1]
nx,ny=np.meshgrid(lons,lats)
x,y=m(nx,ny)
m.quiver(x,y,u2,v2)
plt.show()

kongfeng0824

经纬度的步长是否一致？
如果一直可以先读进工作空间，然后在画图

river

用grads很简单就能出来

mofangbao

问题具体一些

topmad

发在python版里面 是用什么工具画  什么资料 具体一些啊

po_po1

topmad 发表于 2014-2-24 08:32

                               
登录/注册后可看大图

发在python版里面 是用什么工具画  什么资料 具体一些啊

当然是python画了，用到的模块：netCDF4，numpy，Basemap
数据资料：
u分量 u=Dataset(r'C:/Users/Administrator/Desktop/uwnd.2010.nc')
获取u分量维度信息：u.dimensions
OrderedDict([(u'lon', <netCDF4.Dimension object at 0x01FEEBA0>), (u'lat', <netCDF4.Dimension object at 0x039ECE70>), (u'level', <netCDF4.Dimension object at 0x039ECEA0>), (u'time', <netCDF4.Dimension object at 0x039ECED0>)])
u分量的变量信息：u.variables
OrderedDict([(u'level', <netCDF4.Variable object at 0x039F00D0>), (u'lat', <netCDF4.Variable object at 0x039F0120>), (u'lon', <netCDF4.Variable object at 0x039F0170>), (u'time', <netCDF4.Variable object at 0x039F01C0>), (u'uwnd', <netCDF4.Variable object at 0x039F0210>)])
变量uwnd具体维度信息:u.variables['uwnd'].dimensions
(u'time', u'level', u'lat', u'lon')
我选取时间（time）、高度（level）都为第一组数据：u=u.variables['uwnd'][0][0][:]   #获得了uwnd的经纬度数据
#############################################################
v分量 v=Dataset(r'C:/Users/Administrator/Desktop/vwnd.2010.nc')
获取v分量维度信息：v.dimensions
OrderedDict([(u'lon', <netCDF4.Dimension object at 0x0396E300>), (u'lat', <netCDF4.Dimension object at 0x0396E330>), (u'level', <netCDF4.Dimension object at 0x0396E2D0>), (u'time', <netCDF4.Dimension object at 0x0396E390>)])
v分量的变量信息：v.variables
OrderedDict([(u'level', <netCDF4.Variable object at 0x039F0260>), (u'lat', <netCDF4.Variable object at 0x039F02B0>), (u'lon', <netCDF4.Variable object at 0x039F0300>), (u'time', <netCDF4.Variable object at 0x039F0350>), (u'vwnd', <netCDF4.Variable object at 0x039F03A0>)])
变量vwnd具体维度信息:v.variables['vwnd'].dimensions
(u'time', u'level', u'lat', u'lon')

我选取时间（time）、高度（level）都为第一组数据：v=v.variables['vwnd'][0][0][:]   #获得了vwnd的经纬度数据

po_po1

river 发表于 2014-2-23 21:53

                               
登录/注册后可看大图

用grads很简单就能出来

对啊，GraDS我试了一下，很简单就出来了。不过我想尝试下python

po_po1

kongfeng0824 发表于 2014-2-23 20:51

                               
登录/注册后可看大图

经纬度的步长是否一致？
如果一直可以先读进工作空间，然后在画图

哦，我试试

po_po1

风场绘制出来了，虽然没有添加设置，有点粗糙：