MeteoInfoLab脚本示例：台风和热带气旋

MeteoInfo

示例了从一个台风文本文件中读取台风每个节点的经纬度信息并绘制台风路径图，示例数据文件可以从此贴下载：http://bbs.06climate.com/forum.p ... &extra=page%3D1 。

脚本程序：

1. # Read typhoon data file
   
2. fn = 'D:/Temp/ascii/198214track.dat'
   
3. tf = open(fn)
   
4. tf.readline()
   
5. aline = tf.readline()
   
6. id = aline.split()[1]
   
7. tf.readline()
   
8. lon = []
   
9. lat = []
   
10. ws = []
    
11. for aline in tf:
    
12.     print aline
    
13.     datalist = aline.split()
    
14.     lat.append(float(datalist[1]))
    
15.     lon.append(float(datalist[2]))
    
16.     t = datalist[3]
    
17.     ws.append(int(datalist[4]))
    
18.     stat = datalist[6]
    
19.     if len(datalist) == 8:
    
20.         stat = stat + ' ' + datalist[7]
    
21. 
    
22. # Plot
    
23. axesm()
    
24. lworld = shaperead('D:/Temp/map/country1.shp')
    
25. geoshow(lworld, facecolor=[200,200,200])
    
26. plotm(lon, lat, linewidth=2)
    
27. layer = scatterm(lon, lat, ws)
    
28. colorbar(layer, shrink=0.8)
    
29. xlim(110, 140)
    
30. ylim(15, 45)
    
31. title('Typhoon pathway')

# Read typhoon data file<br /> fn = 'D:/Temp/ascii/198214track.dat'<br /> tf = open(fn)<br /> tf.readline()<br /> aline = tf.readline()<br /> id = aline.split()[1]<br /> tf.readline()<br /> lon = []<br /> lat = []<br /> ws = []<br /> for aline in tf:<br />     print aline<br />     datalist = aline.split()<br />     lat.append(float(datalist[1]))<br />     lon.append(float(datalist[2]))<br />     t = datalist[3]<br />     ws.append(int(datalist[4]))<br />     stat = datalist[6]<br />     if len(datalist) == 8:<br />         stat = stat + ' ' + datalist[7]<br /> <br /> # Plot<br /> axesm()<br /> lworld = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(lworld, facecolor=[200,200,200])<br /> plotm(lon, lat, linewidth=2)<br /> layer = scatterm(lon, lat, ws)<br /> colorbar(layer, shrink=0.8)<br /> xlim(110, 140)<br /> ylim(15, 45)<br /> title('Typhoon pathway')

热带气旋的数据文件从此网站上下载：https://climatedataguide.ucar.ed ... one-best-track-data ，数据为netCDF格式，从中读取部分气旋路径的经纬度信息并绘图。

脚本程序：

1. fn = 'D:/Temp/nc/Allstorms.ibtracs_wmo.v03r06.nc'
   
2. f = addfile(fn)
   
3. lons = f['lon_wmo'][:200,:]
   
4. lats = f['lat_wmo'][:200,:]
   
5. data = MIXYListData()
   
6. for i in range(0, lons.dimlen(0)):
   
7.     data.addseries(lons[i,:], lats[i,:])
   
8. 
   
9. # Plot
   
10. axesm()
    
11. mlayer = shaperead('D:/Temp/map/country1.shp')
    
12. geoshow(mlayer)
    
13. layer = plotm(data)
    
14. #layer = scatterm(lons[:,0], lats[:,0], facecolor='b', size=2)
    
15. title('IBTrACS')
    
16. xlim(-180, 180)
    
17. ylim(-90, 90)

fn = 'D:/Temp/nc/Allstorms.ibtracs_wmo.v03r06.nc'<br /> f = addfile(fn)<br /> lons = f['lon_wmo'][:200,:]<br /> lats = f['lat_wmo'][:200,:]<br /> data = MIXYListData()<br /> for i in range(0, lons.dimlen(0)):<br />     data.addseries(lons[i,:], lats[i,:])<br /> <br /> # Plot<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer)<br /> layer = plotm(data)<br /> #layer = scatterm(lons[:,0], lats[:,0], facecolor='b', size=2)<br /> title('IBTrACS')<br /> xlim(-180, 180)<br /> ylim(-90, 90)

当未来

多谢王老师，太赞了

colorlingling

辛苦楼主了,谢谢楼主了

krysunder

收藏，感谢楼主

栗子

画的很漂亮，很强大

rainbow111

努力学习中……多谢老师！

郭小侠V

王老师，您好，我根据您的脚本想画个台风路径图，提示以下错误：



Traceback (most recent call last):
  File "<string>", line 25, in <module>
NameError: name 'axesm' is not defined

MeteoInfo

郭小侠V 发表于 2016-5-11 00:05
王老师，您好，我根据您的脚本想画个台风路径图，提示以下错误：

参考此贴：http://bbs.06climate.com/forum.p ... &extra=page%3D1

郭小侠V

MeteoInfo 发表于 2016-5-11 00:09
参考此贴：http://bbs.06climate.com/forum.php?mod=viewthread&tid=36151&extra=page%3D1

我完全是用您的代码和您的数据，路径也对

Traceback (most recent call last):
  File "<string>", line 25, in <module>
NameError: name 'axesm' is not defined

MeteoInfo

MeteoInfo Java版，MeteoInfoLab