MeteoInfoLab脚本示例：gif动画图

MeteoInfo

在MeteoInfoLab脚本中生成gif动画图要用gifanimation()方法创建一个gif动画：

    gifanimation(filename, repeat=0, delay=1000)

方法中第一个参数是gif动画图片的文件名，后面两个参数是可选的，repeat是动画的重复次数，0指的无限次重复；delay指的动画帧之间的时间间隔，单位是毫秒，delay=1000及时间间隔为1秒。


然后在脚本的循环中给gif动画添加帧，用gifaddframe()方法，方法中的参数为之前创建的gif动画变量。

最后用gif动画变量（此例中为animation）的finish()方法结束动画创建过程，同时会完成gif动画图片的写入。

这个示例脚本中的内容很丰富，演示了CUACE-Haze/Fog模式预报的能见度和站点观测的雾-霾天气符号的叠加，从而直观的判断模式的预报效果。

脚本程序：

1. #Set data folders
   
2. fcstdir = 'W:/CUACE_out/haze_54'
   
3. obsdir = 'U:/data/micaps'
   
4. mapdir = 'T:/verification/map'
   
5. outdir = 'T:/verification/cases/cuace_54km/result'
   
6. 
   
7. #Set time
   
8. st = datetime.datetime(2015,11,1,0)
   
9. sbjt = st + datetime.timedelta(hours=8)
   
10. 
    
11. #Get forecasting data files
    
12. fcstfn = os.path.join(fcstdir, str(st.year) + '/' + st.strftime('%Y%m') + '/' + \
    
13.     st.strftime('%Y%m%d%H') + '/produce/haze_post_aero_' + st.strftime('%Y%m%d%H') + '.ctl')
    
14. print 'Forecasting data file: ' + fcstfn
    
15. fcstf = addfile(fcstfn)
    
16. fcstf.bigendian(True)
    
17. 
    
18. #Plot
    
19. figure(figsize=[768,480], newfig=False)
    
20. axesm(position=[0,0,1,1], xyscale=1.2, tickfontsize=12)
    
21. lworld = shaperead(os.path.join(mapdir, 'country1.shp'))
    
22. lchina = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))
    
23. geoshow(lchina)
    
24. geoshow(lworld, edgecolor='k')
    
25. levs = [0.1, 0.5, 1, 5, 10, 20, 30, 50]
    
26. cols = [(192,0,0),(255,69,0),(255,105,180),(255,128,0),(255,192,128),(51,255,255), \
    
27.     (153,255,153),(204,255,204),(255,255,255)]
    
28. #Set weather list - haze, mist and fog
    
29. weathers = [5,10,11,12,40,41,42,43,44,45,46,47,48,49]
    
30. ls = weatherspec(weathers)
    
31. #Add south China Sea
    
32. lscs = shaperead(os.path.join(mapdir, 'bou1_4l.shp'))
    
33. axesm(position=[0.76,0.09,0.16,0.22], axison=False, xyscale=1.2)
    
34. geoshow(lscs, facecolor=(0,0,255))
    
35. xlim(106, 123)
    
36. ylim(2, 23)
    
37. #Set current plot to 1
    
38. currentplot(1)
    
39. 
    
40. #Create gif animation
    
41. giffn = os.path.join(outdir, 'V_vis_' + st.strftime('%Y%m%d') + '--loop-.gif')
    
42. print giffn
    
43. animation = gifanimation(giffn)
    
44. 
    
45. #Loop
    
46. tnum = fcstf.timenum()
    
47. nn = 0
    
48. for t in range(1, 25):
    
49.     if nn > 0:
    
50.         cll()
    
51.         cll()
    
52.     tt = fcstf.gettime(t)
    
53.     bjt = tt + datetime.timedelta(hours=8)
    
54.     obsfn = os.path.join(obsdir,  str(bjt.year) + '/plot/' + bjt.strftime('%y%m%d%H.000'))
    
55.     if os.path.exists(obsfn):
    
56.         print 'Observation data file: ' + obsfn   
    
57.         obsf = addfile_micaps(obsfn)
    
58.         fdata = fcstf['VIS'][t,:,:]
    
59.         lfcst = contourfm(fdata, levs, colors=cols, proj=fcstf.proj)
    
60.         odata = obsf.stationdata('WeatherNow')
    
61.         lobs = scatterm(odata, symbolspec=ls)
    
62.         title('CUACE/Haze-Fog visibility (km) ' + sbjt.strftime('%Y-%m-%d %H:00') + \
    
63.                 ' +' + str(t*3) + ' (' + bjt.strftime('%Y-%m-%d %H:00') + ')', \
    
64.                 bold=False)
    
65.         xlim(68.6, 140.5)
    
66.         ylim(17.3, 54)
    
67.         colorbar(lfcst, extendrect=False, shrink=0.6)
    
68.         #Add frame to gif animation
    
69.         gifaddframe(animation)
    
70.         nn += 1
    
71. 
    
72. #Finish gif animation
    
73. animation.finish()
    
74. print 'Finished...'

#Set data folders<br /> fcstdir = 'W:/CUACE_out/haze_54'<br /> obsdir = 'U:/data/micaps'<br /> mapdir = 'T:/verification/map'<br /> outdir = 'T:/verification/cases/cuace_54km/result'<br /> <br /> #Set time<br /> st = datetime.datetime(2015,11,1,0)<br /> sbjt = st + datetime.timedelta(hours=8)<br /> <br /> #Get forecasting data files<br /> fcstfn = os.path.join(fcstdir, str(st.year) + '/' + st.strftime('%Y%m') + '/' + \<br />     st.strftime('%Y%m%d%H') + '/produce/haze_post_aero_' + st.strftime('%Y%m%d%H') + '.ctl')<br /> print 'Forecasting data file: ' + fcstfn<br /> fcstf = addfile(fcstfn)<br /> fcstf.bigendian(True)<br /> <br /> #Plot<br /> figure(figsize=[768,480], newfig=False)<br /> axesm(position=[0,0,1,1], xyscale=1.2, tickfontsize=12)<br /> lworld = shaperead(os.path.join(mapdir, 'country1.shp'))<br /> lchina = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))<br /> geoshow(lchina)<br /> geoshow(lworld, edgecolor='k')<br /> levs = [0.1, 0.5, 1, 5, 10, 20, 30, 50]<br /> cols = [(192,0,0),(255,69,0),(255,105,180),(255,128,0),(255,192,128),(51,255,255), \<br />     (153,255,153),(204,255,204),(255,255,255)]<br /> #Set weather list - haze, mist and fog<br /> weathers = [5,10,11,12,40,41,42,43,44,45,46,47,48,49]<br /> ls = weatherspec(weathers)<br /> #Add south China Sea<br /> lscs = shaperead(os.path.join(mapdir, 'bou1_4l.shp'))<br /> axesm(position=[0.76,0.09,0.16,0.22], axison=False, xyscale=1.2)<br /> geoshow(lscs, facecolor=(0,0,255))<br /> xlim(106, 123)<br /> ylim(2, 23)<br /> #Set current plot to 1<br /> currentplot(1)<br /> <br /> #Create gif animation<br /> giffn = os.path.join(outdir, 'V_vis_' + st.strftime('%Y%m%d') + '--loop-.gif')<br /> print giffn<br /> animation = gifanimation(giffn)<br /> <br /> #Loop<br /> tnum = fcstf.timenum()<br /> nn = 0<br /> for t in range(1, 25):<br />     if nn > 0:<br />         cll()<br />         cll()<br />     tt = fcstf.gettime(t)<br />     bjt = tt + datetime.timedelta(hours=8)<br />     obsfn = os.path.join(obsdir,  str(bjt.year) + '/plot/' + bjt.strftime('%y%m%d%H.000'))<br />     if os.path.exists(obsfn):<br />         print 'Observation data file: ' + obsfn    <br />         obsf = addfile_micaps(obsfn)<br />         fdata = fcstf['VIS'][t,:,:]<br />         lfcst = contourfm(fdata, levs, colors=cols, proj=fcstf.proj)<br />         odata = obsf.stationdata('WeatherNow')<br />         lobs = scatterm(odata, symbolspec=ls)<br />         title('CUACE/Haze-Fog visibility (km) ' + sbjt.strftime('%Y-%m-%d %H:00') + \<br />                 ' +' + str(t*3) + ' (' + bjt.strftime('%Y-%m-%d %H:00') + ')', \<br />                 bold=False)<br />         xlim(68.6, 140.5)<br />         ylim(17.3, 54)<br />         colorbar(lfcst, extendrect=False, shrink=0.6)<br />         #Add frame to gif animation<br />         gifaddframe(animation)<br />         nn += 1<br /> <br /> #Finish gif animation<br /> animation.finish()<br /> print 'Finished...'

taxueyueming

王老师真厉害，赞一个

梁的丰

王老师了不起...

aas314

好东西 收藏了

锶兎扌斯书

{:eb511:}{:eb511:}

rceclx

很漂亮的动画图，有时间要认真学习！谢谢王老师！

紫宸

王老师，请问最后一张图那个第二个读取的.shp文件是什么呢？在哪里可以下到这个数据

幽居寒舍

厉害了我的王老师，然而，我看不懂

ddd212

{:5_235:}

veelam

厉害了王老师