MeteoInfoLab脚本示例：多色线条

MeteoInfo

对于线对象的绘制每条线只能是固定的颜色，要实现用不同的颜色绘制一条线需要将该线分隔为多条线来绘制。这里示例读取HYSPLIT输出的气团轨迹数据，生成轨迹图层，在生成轨迹图层的时候每两个节点就生成一条线，并将轨迹节点的高度作为属性加入图层中。利用轨迹高度数据形成色标来显示轨迹，达到一条轨迹不同高度用不同颜色显示的目的。

脚本程序：

1. #Read trajectory data file
   
2. fn = 'E:/Trajectory/RYO-to Yaqiang/05060306'
   
3. f = open(fn)
   
4. line = f.readline()
   
5. mfn = int(line.split()[0])
   
6. for i in range(mfn):
   
7.     f.readline()
   
8. line = f.readline()
   
9. sn = int(line.split()[0])
   
10. for i in range(sn):
    
11.     f.readline()
    
12. f.readline()
    
13. lons = []
    
14. lats = []
    
15. alts = []
    
16. for line in f:
    
17.     dd = line.split()
    
18.     lat = float(dd[9])
    
19.     lon = float(dd[10])
    
20.     alt = float(dd[11])
    
21.     lats.append(lat)
    
22.     lons.append(lon)
    
23.     alts.append(alt)
    
24. alts = array(alts)
    
25. f.close()
    
26. #Create trajectory layer
    
27. layer = MILayer(shapetype='line')
    
28. layer.addfield('Height', 'float')
    
29. lon1 = lons[0]
    
30. lat1 = lats[0]
    
31. for j in range(1, len(lons)):
    
32.     lon = lons[j]
    
33.     lat = lats[j]
    
34.     alt = alts[j]
    
35.     layer.addshape([lon1,lon], [lat1, lat], [alt])
    
36.     lon1 = lon
    
37.     lat1 = lat
    
38. levs = arange(0, 1000, 100)
    
39. cols = makecolors(len(levs) + 1)
    
40. ls = makesymbolspec('line', levels=levs, colors=cols, field='Height', size=2)
    
41. #Plot
    
42. axesm()
    
43. geoshow('country', edgecolor=(204,204,204), facecolor=(250,235,252))
    
44. geoshow(layer, symbolspec=ls)
    
45. scatterm(lons[0], lats[0], size=8, facecolor='r')
    
46. text(lons[0], lats[0]-1, 'Fukushima')
    
47. colorbar(layer)
    
48. title('Multi-color line')
    
49. yticks(arange(30, 51, 5))
    
50. xlim(130, 160)
    
51. ylim(30, 50)

#Read trajectory data file<br /> fn = 'E:/Trajectory/RYO-to Yaqiang/05060306'<br /> f = open(fn)<br /> line = f.readline()<br /> mfn = int(line.split()[0])<br /> for i in range(mfn):<br />     f.readline()<br /> line = f.readline()<br /> sn = int(line.split()[0])<br /> for i in range(sn):<br />     f.readline()<br /> f.readline()<br /> lons = []<br /> lats = []<br /> alts = []<br /> for line in f:<br />     dd = line.split()<br />     lat = float(dd[9])<br />     lon = float(dd[10])<br />     alt = float(dd[11])<br />     lats.append(lat)<br />     lons.append(lon)<br />     alts.append(alt)<br /> alts = array(alts)<br /> f.close()<br /> #Create trajectory layer<br /> layer = MILayer(shapetype='line')<br /> layer.addfield('Height', 'float')<br /> lon1 = lons[0]<br /> lat1 = lats[0]<br /> for j in range(1, len(lons)):<br />     lon = lons[j]<br />     lat = lats[j]<br />     alt = alts[j]<br />     layer.addshape([lon1,lon], [lat1, lat], [alt])<br />     lon1 = lon<br />     lat1 = lat<br /> levs = arange(0, 1000, 100)<br /> cols = makecolors(len(levs) + 1)<br /> ls = makesymbolspec('line', levels=levs, colors=cols, field='Height', size=2)<br /> #Plot<br /> axesm()<br /> geoshow('country', edgecolor=(204,204,204), facecolor=(250,235,252))<br /> geoshow(layer, symbolspec=ls)<br /> scatterm(lons[0], lats[0], size=8, facecolor='r')<br /> text(lons[0], lats[0]-1, 'Fukushima')<br /> colorbar(layer)<br /> title('Multi-color line')<br /> yticks(arange(30, 51, 5))<br /> xlim(130, 160)<br /> ylim(30, 50)

MeteoInfoLab有专门的函数读取HYSPLIT输出的轨迹文件并生成轨迹图层，脚本也可以写成这样：

1. #Read trajectory data file
   
2. fn = 'D:/Temp/HYSPLIT/05060306'
   
3. f = addfile_hytraj(fn)
   
4. layer = f.trajlayer()
   
5. stlayer = f.trajsplayer()
   
6. spoint = stlayer.shapes()[0].getPoint()
   
7. slon = spoint.X
   
8. slat = spoint.Y
   
9. levs = arange(0, 1000, 100)
   
10. cols = makecolors(len(levs) + 1)
    
11. ls = makesymbolspec('line', levels=levs, colors=cols, field='Geometry_Z', size=2)
    
12. 
    
13. #Plot
    
14. axesm()
    
15. geoshow('country', edgecolor=(204,204,204), facecolor=(250,235,252))
    
16. geoshow(layer, symbolspec=ls)
    
17. ss = makesymbolspec('point', {'marker':'S','size':12})
    
18. geoshow(stlayer, symbolspec=ss)
    
19. text(slon, slat-1, 'Fukushima')
    
20. colorbar(layer)
    
21. title('Multi-color line (Geometry_Z field)')
    
22. yticks(arange(30, 51, 5))
    
23. xlim(130, 160)
    
24. ylim(30, 50)

#Read trajectory data file<br /> fn = 'D:/Temp/HYSPLIT/05060306'<br /> f = addfile_hytraj(fn)<br /> layer = f.trajlayer()<br /> stlayer = f.trajsplayer()<br /> spoint = stlayer.shapes()[0].getPoint()<br /> slon = spoint.X<br /> slat = spoint.Y<br /> levs = arange(0, 1000, 100)<br /> cols = makecolors(len(levs) + 1)<br /> ls = makesymbolspec('line', levels=levs, colors=cols, field='Geometry_Z', size=2)<br /> <br /> #Plot<br /> axesm()<br /> geoshow('country', edgecolor=(204,204,204), facecolor=(250,235,252))<br /> geoshow(layer, symbolspec=ls)<br /> ss = makesymbolspec('point', {'marker':'S','size':12})<br /> geoshow(stlayer, symbolspec=ss)<br /> text(slon, slat-1, 'Fukushima')<br /> colorbar(layer)<br /> title('Multi-color line (Geometry_Z field)')<br /> yticks(arange(30, 51, 5))<br /> xlim(130, 160)<br /> ylim(30, 50)

lyysms

诶...突然发现也可以在map里直接画...

aas314

好漂亮  感谢老师

breathfreelylh

老师，请问hysplit输出的气团轨迹数据是在setup run里那个输出的tdmp文件么，是用text打开么

shanie

breathfreelylh 发表于 2016-11-4 20:39
老师，请问hysplit输出的气团轨迹数据是在setup run里那个输出的tdmp文件么，是用text打开么

你会画了么同求

linlinouc

认真学习，感谢分享！

突破的渴望

王老师您好 这个代码要用的话 需要修改哪些地方啊？谢谢王老师~~

lyysms

  File "<iostream>", line 32
    lon = lons<i>
               
SyntaxError: no viable alternative at input '\n'
出现了这个错误...

MeteoInfo

lyysms 发表于 2018-6-23 15:33
File "", line 32
    lon = lons
               

<> 改为 [] ，这是论坛贴代码的bug.

lyysms

谢谢老师回复，但是改完后又出现了这个问题
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
NameError: name 'MILayer' is not defined

lyysms

Traceback (most recent call last):
  File "<iostream>", line 27, in <module>
NameError: name 'MILayer' is not defined