求助，MeteoInfoLab脚本改如何写

棉花猪仔

情况是这样的
有站点的U,V分量，想画涡度、散度等值线图。


我的理解是，先把站点的U,V分量插值成格点数据然后计算涡度，散度绘图。


站点数据我是用这样的格式

113.894 22.65 0.0
113.802 22.66 0.5
114.211 22.686 0.0

...



u,v各一个文件  例如u.dat   v.dat

然后我是这样写的fnu = 'D:/uv/u.dat'fnv = 'D:/uv/v.dat'
ncol = numasciicol(fn)
nrow = numasciirow(fn)
ncol2 = numasciicol(fnv)
nrow2 = numasciirow(fnv)
a = asciiread(fn,shape=(nrow,ncol))
a2 = asciiread(fnv,shape=(nrow2,ncol2))
#读取文件lon = a[:,0]
lat = a[:,1]
u = a[:,2]
lon2 = a2[:,0]
lat2 = a2[:,1]
v = a2[:,2]
#转站点格式
uu=asstationdata(u, lon, lat, fill_value=-9999.0)
vv=asstationdata(v, lon2, lat2, fill_value=-9999.0)
x = arange(113.5415067, 114.4700278,0.009)
y = arange(22.33, 22.83716667, 0.009)
#插值prg = uu.griddata((x, y), method='idw', radius=3)
prg2 = vv.griddata((x, y), method='idw', radius=3)

#到这里不知道怎么写了，prg 是PyGridData对象，但是hcurl里的是DimArray
vort = hcurl(prg, prg2)
axesm()
mlayer = shaperead('D:/uv/china.shp')
geoshow(mlayer)
layer = contourfm(prg, 20)
#layer = contourfm(divg, 20)
title('Vorticity')
#title('Divergence')
colorbar(layer, orientation='horizontal', aspect=50)

------------------------------------------------------------------------------------------------------------------------
没学过python脚本，都是看论坛里的示例来画的葫芦，求指教。。。

MeteoInfo

给PyGridData增加了asdimarray()函数。需要下载MeteoInfo最新版：http://yun.baidu.com/share/link?shareid=669776748&uk=51062435

MeteoInfo

明白你的需求了，有时间的时候我会加上PyGridData转为DimArray的功能。

棉花猪仔

非常感谢版主的帮忙!!  有机会必须请你吃饭