MeteoInfoLab怎么用micaps第2类数据画流线？

zyfws

我的代码如下：
f = addfile_micaps(u'D:/Data/09070220.000')
u = f['WindDirection'][:]
v = f['WindSpeed'][:]
axesm()
mlayer = shaperead(r'D:/Program Files (x86)/MeteoInfo/map/country1.shp')
geoshow(mlayer, edgecolor='k')
layer = streamplotm(u, v)
title('Streamline')

结果提示：
indices must be 2 dimensions!
indices must be 2 dimensions!
Traceback (most recent call last):
  File "<iostream>", line 7, in <module>
  File "D:\Program Files (x86)\MeteoInfo\pylib\mipylib\plotlib\miplot.py", line 4865, in streamplotm
    layer = __plot_uvgriddata_m(plot, udata, vdata, None, ls, 'streamplot', isuv, proj=proj, density=density)
  File "D:\Program Files (x86)\MeteoInfo\pylib\mipylib\plotlib\miplot.py", line 4980, in __plot_uvgriddata_m
    layer = DrawMeteoData.createStreamlineLayer(udata.data, vdata.data, density, ls, 'layer', isuv)
AttributeError: 'NoneType' object has no attribute 'data'

请王老师帮看看怎么解决？

MeteoInfo

目前MeteoInfoLab对于micaps2类数据的处理能力还有限，会尽快加强这方面的功能。

zyfws

MeteoInfo 发表于 2017-2-4 18:10
目前MeteoInfoLab对于micaps2类数据的处理能力还有限，会尽快加强这方面的功能。

那请问Java版的MeteoInfo怎么画流线图，软件界面没有流线类型选项，论坛搜到的脚本的实例是C#版的

MeteoInfo

zyfws 发表于 2017-2-6 15:16
那请问Java版的MeteoInfo怎么画流线图，软件界面没有流线类型选项，论坛搜到的脚本的实例是C#版的

两个版本的操作类似。流线分析只能针对格点数据，站点数据必须先插值为格点数据才能进行流线分析。

MeteoInfo

已经升级了MeteoInfo版本为1.4R3，提升了MeteoInfoLab处理micaps2类数据的能力。示例脚本如下：

1. fn = 'D:/Temp/micaps/10120420.000'
   
2. f = addfile_micaps(fn)
   
3. lon = f['Longitude'][:]
   
4. lat = f['Latitude'][:]
   
5. wd = f['WindDirection'][:]
   
6. ws = f['WindSpeed'][:]
   
7. 
   
8. #To grid data
   
9. x = arange(70, 140, 2)
   
10. y = arange(15, 58, 2)
    
11. gwd,gx,gy = griddata((lon, lat), wd, xi=(x, y), method='idw')
    
12. gws,gx,gy = griddata((lon, lat), ws, xi=(x, y), method='idw')
    
13. 
    
14. #Plot
    
15. axesm()
    
16. lworld = shaperead('D:/Temp/map/country1.shp')
    
17. geoshow(lworld, edgecolor='k')
    
18. layer = streamplotm(x, y, gwd, gws, color='b', isuv=False)

fn = 'D:/Temp/micaps/10120420.000'<br /> f = addfile_micaps(fn)<br /> lon = f['Longitude'][:]<br /> lat = f['Latitude'][:]<br /> wd = f['WindDirection'][:]<br /> ws = f['WindSpeed'][:]<br /> <br /> #To grid data<br /> x = arange(70, 140, 2)<br /> y = arange(15, 58, 2)<br /> gwd,gx,gy = griddata((lon, lat), wd, xi=(x, y), method='idw')<br /> gws,gx,gy = griddata((lon, lat), ws, xi=(x, y), method='idw')<br /> <br /> #Plot<br /> axesm()<br /> lworld = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(lworld, edgecolor='k')<br /> layer = streamplotm(x, y, gwd, gws, color='b', isuv=False)

zyfws

MeteoInfo 发表于 2017-2-6 19:39
已经升级了MeteoInfo版本为1.4R3，提升了MeteoInfoLab处理micaps2类数据的能力。示例脚本如下：

很好用，谢谢老师

TL01

过来学习学习