MeteoInfoLab脚本示例：Trajectory

MeteoInfo

示例读取HYSPLIT模式输出的气团轨迹数据文件，生成轨迹图层，并显示轨迹各节点的气压图。

脚本程序：

1. from org.meteoinfo.geometry.legend import PointStyle
   
2. fn = os.path.join(migl.get_sample_folder(), 'HYSPLIT', 'tdump')
   
3. f = addfile_hytraj(fn)
   
4. lon = f['lon'][:,:]
   
5. lat = f['lat'][:,:]
   
6. 
   
7. #Plot
   
8. figure(figsize=[526, 489], newfig=False)
   
9. axesm(position=[0.12, 0.3, 0.85, 0.7])
   
10. geoshow('country', edgecolor=(0,0,255), facecolor=(230,230,230))
    
11. cols = makecolors(len(lon))
    
12. levs = arange(0, 101, 5)
    
13. tlayer = plot(lon, lat, levels=levs, isadd=False)
    
14. ls = tlayer.update_legend('unique', 'ID')
    
15. ii = 0
    
16. for lb in ls.getLegendBreaks():
    
17.     lb.setSize(2)
    
18.     lb.setDrawSymbol(True)
    
19.     lb.setSymbolInterval(6)
    
20.     if ii == len(PointStyle.values()):
    
21.         ii = 0
    
22.     lb.setSymbolStyle(PointStyle.values()[ii])
    
23.     ii += 1
    
24. geoshow(tlayer)
    
25. scatter(lon[:,0], lat[:,0], size=6, color='r', marker='S')
    
26. xlim(-92, -55)
    
27. ylim(34, 54)
    
28. yticks(arange(35, 54, 5))
    
29. title('MeteoInfoLab script demo - Trajectory')
    
30. 
    
31. axes(outerposition=[0, 0, 1, 0.3], yreverse=True, xaxistype='time')
    
32. tt = f['time'][:,:]
    
33. data = f['PRESSURE'][:,:]
    
34. plot(tt, data, legend=tlayer.legend())
    
35. xlabel('Time')
    
36. ylabel('hPa')

from org.meteoinfo.geometry.legend import PointStyle<br /> fn = os.path.join(migl.get_sample_folder(), 'HYSPLIT', 'tdump')<br /> f = addfile_hytraj(fn)<br /> lon = f['lon'][:,:]<br /> lat = f['lat'][:,:]<br /> <br /> #Plot<br /> figure(figsize=[526, 489], newfig=False)<br /> axesm(position=[0.12, 0.3, 0.85, 0.7])<br /> geoshow('country', edgecolor=(0,0,255), facecolor=(230,230,230))<br /> cols = makecolors(len(lon))<br /> levs = arange(0, 101, 5)<br /> tlayer = plot(lon, lat, levels=levs, isadd=False)<br /> ls = tlayer.update_legend('unique', 'ID')<br /> ii = 0<br /> for lb in ls.getLegendBreaks():<br />     lb.setSize(2)<br />     lb.setDrawSymbol(True)<br />     lb.setSymbolInterval(6)<br />     if ii == len(PointStyle.values()):<br />         ii = 0<br />     lb.setSymbolStyle(PointStyle.values()[ii])<br />     ii += 1<br /> geoshow(tlayer)<br /> scatter(lon[:,0], lat[:,0], size=6, color='r', marker='S')<br /> xlim(-92, -55)<br /> ylim(34, 54)<br /> yticks(arange(35, 54, 5))<br /> title('MeteoInfoLab script demo - Trajectory')<br /> <br /> axes(outerposition=[0, 0, 1, 0.3], yreverse=True, xaxistype='time')<br /> tt = f['time'][:,:]<br /> data = f['PRESSURE'][:,:]<br /> plot(tt, data, legend=tlayer.legend())<br /> xlabel('Time')<br /> ylabel('hPa')

cfltwy

王老师王老师  和  各位前辈们，我只想要线条，不希望有任何的symbol 标志，怎么修改程序哇？  在线等哇。。。

MeteoInfo

cfltwy 发表于 2017-6-5 12:15
王老师王老师  和  各位前辈们，我只想要线条，不希望有任何的symbol 标志，怎么修改程序哇？  在线等哇。 ...

在
geoshow(tlayer)
语句之前加上：
ls = tlayer.legend()
for lb in ls.getLegendBreaks():
    lb.setDrawSymbol(False)

cfltwy

MeteoInfo 发表于 2017-6-5 12:30
在
geoshow(tlayer)
语句之前加上：

非常感谢王老师的快速回复。  还有一个问题想问您， 因为我下面的时间-高度变化一维图用的是altitude参数， 所以需要将其y轴反转一下，  反转如何实现哇？

MeteoInfo

cfltwy 发表于 2017-6-5 12:40
非常感谢王老师的快速回复。  还有一个问题想问您， 因为我下面的时间-高度变化一维图用的是altitude参数 ...

axes(position=[0.1,0.1,0.8,0.3], yreverse=True, xaxistype='time')
将yreverse参数去掉

cfltwy

MeteoInfo 发表于 2017-6-5 12:57
axes(position=[0.1,0.1,0.8,0.3], yreverse=True, xaxistype='time')
将yreverse参数去掉

  真不应该没注意到 yreverse这个参数。。 王老师，谢谢您！  

wenwen123

厉害了，楼楼

偷影子的人

王老师，我画这个图的时候，只有三条轨迹线，我想把它们设置对比分明的颜色，但是默认的三条轨迹线颜色非常相近，请问您我应该怎么改啊？

偷影子的人

cfltwy 发表于 2017-6-5 12:40
非常感谢王老师的快速回复。  还有一个问题想问您， 因为我下面的时间-高度变化一维图用的是altitude参数 ...

请问你是怎么把时间-高度变化这个图改为altitude参数的？

MeteoInfo

偷影子的人 发表于 2017-9-24 08:53
请问你是怎么把时间-高度变化这个图改为altitude参数的？

data = f.trajvardata(12)
改为：
data = f.trajvardata(13)