MeteoInfoLab脚本示例：追踪台风路径

MeteoInfo

利用海平面最低气压可以追踪台风路径（http://bbs.06climate.com/forum.p ... p;extra=&page=1）。从数据中读出海平面气压数组后可以用argmin函数获取最小值所在的网格的序号，然后用unravel_index函数获得最小值在数据各维中的序号，获取最小气压所在位置的经纬度，绘制台风路径图。

1. fn = 'D:/Temp/GrADS/928slp.ctl'
   
2. f = addfile(fn)
   
3. slp = f['slvl']
   
4. lat = slp.dimvalue(1)
   
5. lon = slp.dimvalue(2)
   
6. tt = f.gettimes()
   
7. tlat = []
   
8. tlon = []
   
9. press = []
   
10. for t in range(f.timenum()-1):
    
11.     data = slp[t,:,:]
    
12.     idx = data.argmin()
    
13.     idx = unravel_index(idx, data.shape)
    
14.     tlat.append(lat[idx[0]])
    
15.     tlon.append(lon[idx[1]])
    
16.     press.append(data[idx[0],idx[1]])
    
17. 
    
18. #Plot
    
19. axesm(bgcolor=(204,255,255))
    
20. geoshow('country', edgecolor=[200,200,200], facecolor=(255,251,195))
    
21. plotm(tlon, tlat)
    
22. layer = scatterm(tlon, tlat, press)
    
23. colorbar(layer)
    
24. for t in range(0, f.timenum()-1, 4):
    
25.     text(tlon[t]+0.1, tlat[t], tt[t].strftime('%d:%H'))
    
26. xlim(105, 120)
    
27. ylim(15, 25)
    
28. title('Typoon path tracing')

fn = 'D:/Temp/GrADS/928slp.ctl'<br /> f = addfile(fn)<br /> slp = f['slvl']<br /> lat = slp.dimvalue(1)<br /> lon = slp.dimvalue(2)<br /> tt = f.gettimes()<br /> tlat = []<br /> tlon = []<br /> press = []<br /> for t in range(f.timenum()-1):<br />     data = slp[t,:,:]<br />     idx = data.argmin()<br />     idx = unravel_index(idx, data.shape)<br />     tlat.append(lat[idx[0]])<br />     tlon.append(lon[idx[1]])<br />     press.append(data[idx[0],idx[1]])<br /> <br /> #Plot<br /> axesm(bgcolor=(204,255,255))<br /> geoshow('country', edgecolor=[200,200,200], facecolor=(255,251,195))<br /> plotm(tlon, tlat)<br /> layer = scatterm(tlon, tlat, press)<br /> colorbar(layer)<br /> for t in range(0, f.timenum()-1, 4):<br />     text(tlon[t]+0.1, tlat[t], tt[t].strftime('%d:%H'))<br /> xlim(105, 120)<br /> ylim(15, 25)<br /> title('Typoon path tracing')

小小气象预报员

领教领教啦~最近在学习meteolab，准备把王老师的文章都看了

石青love李鑫淼

王老师，我现在的数据是dat格式的呀，按照您的程序写了，提示错误fn = 'D:/928slp.dat'
f = addfile (fn)
slp = f['slvl']
lat = slp.dimvalue(1)
lon = slp.dimvalue(2)
tt = f.gettimes()
tlat = []
tlon = []
press = []
for t in range(f.timenum()-1):
      data =slp[t,:,:]
      idx = data.argmin()
      idx = unraval_index(idx, data.shape)
      tlat.append(lat[idx[0]])
      tlon.append(lon[idx[1]])
      press.append(data[idx[0],idx[1]])

#plot   
      axesm(bgcolor = (204, 255, 255))
      geoshow('countey', edgecolor=[250, 250, 250], facecolor=(255, 251, 195))
      plotm(tlon, tlat)
      layer = scatterm(tlon, tlat, press)
      colorbar(layer)
      for t in range(0, f.timenum()-1, 4):
       text(tlon[t]+0.1, tlat[t], tt[t].strftime('%d:%H'))
      xlim(105, 120)
      ylim(15, 25)
      title('Tyhpoon path tracing')

   

石青love李鑫淼

Traceback (most recent call last):
  File "<iostream>", line 2, in <module>
  File "D:\MeteoInfo\pylib\mipylib\dataset\midata.py", line 112, in addfile
    meteodata.openData(fname, keepopen)
        at java.lang.NumberFormatException.forInputString(Unknown Source)

        at java.lang.Integer.parseInt(Unknown Source)

        at java.lang.Short.parseShort(Unknown Source)

        at java.lang.Short.parseShort(Unknown Source)

        at org.meteoinfo.data.meteodata.arl.ARLDataInfo.readDataLabel(ARLDataInfo.java:582)

        at org.meteoinfo.data.meteodata.arl.ARLDataInfo.canOpen(ARLDataInfo.java:215)

        at org.meteoinfo.data.meteodata.MeteoDataInfo.openData(MeteoDataInfo.java:462)

        at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)

        at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)

        at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)

        at java.lang.reflect.Method.invoke(Unknown Source)

java.lang.NumberFormatException: java.lang.NumberFormatException: For input string: "<"