MeteoInfoLab脚本示例：台风最佳路径

MeteoInfo

演示MeteoInfoLab读取中央台下发的台风最佳路径文件并绘图。数据文件可以在这里下载：http://tcdata.typhoon.gov.cn/zjljsjj_zlhq.html

1. # Read typhoon data file
   
2. fn = 'D:/Temp/ascii/CH2015BST.txt'
   
3. tf = open(fn)
   
4. lons = []
   
5. lats = []
   
6. prss = []
   
7. for line in tf:
   
8.     print line
   
9.     data = line.split()
   
10.     pn = int(data[2])
    
11.     for i in range(pn):
    
12.         line = tf.readline()
    
13.         data = line.split()
    
14.         lat = float(data[2])        
    
15.         lats.append(lat * 0.1)
    
16.         lon = float(data[3])
    
17.         lons.append(lon * 0.1)
    
18.         t = data[0]
    
19.         prs = float(data[4])
    
20.         prss.append(prs)
    
21.     lons.append(nan)
    
22.     lats.append(nan)
    
23.     prss.append(nan)
    
24. 
    
25. # Plot
    
26. axesm()
    
27. lworld = shaperead('D:/Temp/map/country1.shp')
    
28. geoshow(lworld, facecolor=[200,200,200])
    
29. plotm(lons, lats, linewidth=2)
    
30. #geoshow(lats, lons, color='r')
    
31. layer = scatterm(lons, lats, prss)
    
32. colorbar(layer, shrink=0.8)
    
33. xlim(100, 210)
    
34. ylim(0, 60)
    
35. title('Typhoon pathway')

# Read typhoon data file<br /> fn = 'D:/Temp/ascii/CH2015BST.txt'<br /> tf = open(fn)<br /> lons = []<br /> lats = []<br /> prss = []<br /> for line in tf:<br />     print line<br />     data = line.split()<br />     pn = int(data[2])<br />     for i in range(pn):<br />         line = tf.readline()<br />         data = line.split()<br />         lat = float(data[2])        <br />         lats.append(lat * 0.1)<br />         lon = float(data[3])<br />         lons.append(lon * 0.1)<br />         t = data[0]<br />         prs = float(data[4])<br />         prss.append(prs)<br />     lons.append(nan)<br />     lats.append(nan)<br />     prss.append(nan)<br /> <br /> # Plot<br /> axesm()<br /> lworld = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(lworld, facecolor=[200,200,200])<br /> plotm(lons, lats, linewidth=2)<br /> #geoshow(lats, lons, color='r')<br /> layer = scatterm(lons, lats, prss)<br /> colorbar(layer, shrink=0.8)<br /> xlim(100, 210)<br /> ylim(0, 60)<br /> title('Typhoon pathway')

台风路径线条不同颜色表示风速：

1. #Create typhoon layer
   
2. layer = maplayer(shapetype='line')
   
3. 
   
4. # Read typhoon data file
   
5. fn = 'D:/Temp/ascii/CH2015BST.txt'
   
6. tf = open(fn)
   
7. lons = []
   
8. lats = []
   
9. prss = []
   
10. wss = []
    
11. for line in tf:
    
12.     print line
    
13.     data = line.split()
    
14.     pn = int(data[2])
    
15.     for i in range(pn):
    
16.         line = tf.readline()
    
17.         data = line.split()
    
18.         lat = float(data[2])        
    
19.         lats.append(lat * 0.1)
    
20.         lon = float(data[3])
    
21.         lons.append(lon * 0.1)
    
22.         t = data[0]
    
23.         prs = float(data[4])
    
24.         prss.append(prs)
    
25.         ws = float(data[5])
    
26.         wss.append(ws)
    
27.     layer.addshape(lons, lats, z=prss, m=wss)
    
28.     lons = []
    
29.     lats = []
    
30.     prss = []
    
31.     wss = []
    
32. 
    
33. #Set typhoon layer legend
    
34. levs = arange(5, 61, 5)
    
35. cols = makecolors(len(levs) + 1)
    
36. ls = makesymbolspec('line', levels=levs, colors=cols, field='Geometry_M', size=2)
    
37. 
    
38. # Plot
    
39. axesm()
    
40. geoshow('country', facecolor=[200,200,200])
    
41. geoshow(layer, symbolspec=ls)
    
42. colorbar(layer, shrink=0.8, label='Wind speed (m/s)')
    
43. xlim(100, 210)
    
44. ylim(0, 60)
    
45. title('Typhoon pathway with wind speed')

#Create typhoon layer<br /> layer = maplayer(shapetype='line')<br /> <br /> # Read typhoon data file<br /> fn = 'D:/Temp/ascii/CH2015BST.txt'<br /> tf = open(fn)<br /> lons = []<br /> lats = []<br /> prss = []<br /> wss = []<br /> for line in tf:<br />     print line<br />     data = line.split()<br />     pn = int(data[2])<br />     for i in range(pn):<br />         line = tf.readline()<br />         data = line.split()<br />         lat = float(data[2])        <br />         lats.append(lat * 0.1)<br />         lon = float(data[3])<br />         lons.append(lon * 0.1)<br />         t = data[0]<br />         prs = float(data[4])<br />         prss.append(prs)<br />         ws = float(data[5])<br />         wss.append(ws)<br />     layer.addshape(lons, lats, z=prss, m=wss)<br />     lons = []<br />     lats = []<br />     prss = []<br />     wss = []<br /> <br /> #Set typhoon layer legend<br /> levs = arange(5, 61, 5)<br /> cols = makecolors(len(levs) + 1)<br /> ls = makesymbolspec('line', levels=levs, colors=cols, field='Geometry_M', size=2)<br /> <br /> # Plot<br /> axesm()<br /> geoshow('country', facecolor=[200,200,200])<br /> geoshow(layer, symbolspec=ls)<br /> colorbar(layer, shrink=0.8, label='Wind speed (m/s)')<br /> xlim(100, 210)<br /> ylim(0, 60)<br /> title('Typhoon pathway with wind speed')

zenghaode

谢谢王老师，您的效率真高，这么快就解决了。不过我运行您的程序，结果出现了前一个台风结束点与后一个台风开始点用直线连接，这是怎么回事，我没有动程序，该怎么调整呢，谢谢！

MeteoInfo

zenghaode 发表于 2016-9-25 21:34
谢谢王老师，您的效率真高，这么快就解决了。不过我运行您的程序，结果出现了前一个台风结束点与后一个台风 ...

你需要下载最新版本

zenghaode

好的，谢谢您！ 马上下一个，我刚用的是1.3.5版的。另外，提个建议，要是把位置圆点符号换成台风符号，效果可能会更好一些。

zenghaode

我用的1.3.6版，这个是最新的吧，但还是一样的，怎么回事呢？

MeteoInfo

zenghaode 发表于 2016-9-25 22:09
我用的1.3.6版，这个是最新的吧，但还是一样的，怎么回事呢？

1.3.6R2

KIMO23

王老师  您好  不知道您可以分享一下您所使用的台风最佳路径文件档案吗?想知道每个台风如何做区隔，谢谢

MeteoInfo

KIMO23 发表于 2016-9-28 16:39
王老师  您好  不知道您可以分享一下您所使用的台风最佳路径文件档案吗?想知道每个台风如何做区隔，谢谢

在本帖1楼增加了数据下载地址

KIMO23

感谢老师您的网址分享

yinger

lworld = shaperead('D:/Temp/map/country1.shp')

请问这个数据在哪里下载？