MeteoInfoLab脚本示例:Trajectory

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

脚本程序:

f = addfile_hytraj('D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo/sample/HYSPLIT/tdump')
tlayer = f.trajlayer()
stlayer = f.trajsplayer()
axesm(position=[0.1,0.4,0.8,0.6])
mlayer = shaperead('D:/Temp/map/country1.shp')
geoshow(mlayer, edgecolor=(0,0,255), facecolor=(230,230,230))
geoshow(tlayer)
ss = makesymbolspec('point', {'marker':'S'})
geoshow(stlayer, symbolspec=ss)
xlim(-92, -55)
ylim(34, 54)
title('MeteoInfoLab script demo - Trajectory')

axes(position=[0.1,0.1,0.8,0.3], yreverse=True, xaxistype='time')
data = f.trajvardata(12)
plot(data, legend=tlayer.legend())
xlabel('Time')
ylabel('hPa')

原文地址:https://www.cnblogs.com/yaqiang/p/4850157.html