用meteoinfo制作gif时 图像出现叠加而不是更新

cny

源码是参考官网帖子，也添加了cll()，但还是错误
求教

cny

{:lxm_27:} 求教(～￣▽￣)～*

MeteoInfo

参考此脚本：

1. #Set data folders
   
2. obsdir = 'U:/data/micaps'
   
3. mapdir = 'T:/verification/map'
   
4. outdir = 'C:/Temp/test'
   
5. 
   
6. #Set time
   
7. st = datetime.datetime(2015,11,1,0)
   
8. sbjt = st + datetime.timedelta(hours=8)
   
9. 
   
10. #Plot
    
11. figure(figsize=[768,480], newfig=False)
    
12. axesm(position=[0,0,1,1], xyscale=1.2, tickfontsize=12)
    
13. lworld = shaperead(os.path.join(mapdir, 'country1.shp'))
    
14. lchina = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))
    
15. geoshow(lchina)
    
16. geoshow(lworld, edgecolor='k')
    
17. #Set weather list - haze, mist and fog
    
18. weathers = [5,10,11,12,40,41,42,43,44,45,46,47,48,49]
    
19. ls = weatherspec(weathers)
    
20. #Add south China Sea
    
21. lscs = shaperead(os.path.join(mapdir, 'bou1_4l.shp'))
    
22. axesm(position=[0.76,0.09,0.16,0.22], axison=False, xyscale=1.2)
    
23. geoshow(lscs, facecolor=(0,0,255))
    
24. xlim(106, 123)
    
25. ylim(2, 23)
    
26. #Set current plot to 1
    
27. currentplot(1)
    
28. 
    
29. #Create gif animation
    
30. giffn = os.path.join(outdir, 'V_vis_' + st.strftime('%Y%m%d') + '--loop-.gif')
    
31. print giffn
    
32. animation = gifanimation(giffn)
    
33. 
    
34. #Loop
    
35. nn = 0
    
36. bjt = sbjt
    
37. for t in range(1, 25):
    
38.     if nn > 0:
    
39.         cll()
    
40.     bjt = bjt + datetime.timedelta(hours=3)
    
41.     obsfn = os.path.join(obsdir,  str(bjt.year) + '/plot/' + bjt.strftime('%y%m%d%H.000'))
    
42.     if os.path.exists(obsfn):
    
43.         print 'Observation data file: ' + obsfn
    
44.         obsf = addfile_micaps(obsfn)
    
45.         odata = obsf.stationdata('WeatherNow')
    
46.         lobs = scatterm(odata, symbolspec=ls)
    
47.         title('CUACE/Haze-Fog visibility (km) ' + sbjt.strftime('%Y-%m-%d %H:00') + \
    
48.                 ' +' + str(t*3) + ' (' + bjt.strftime('%Y-%m-%d %H:00') + ')', \
    
49.                 bold=False)
    
50.         xlim(68.6, 140.5)
    
51.         ylim(17.3, 54)
    
52.         #Add frame to gif animation
    
53.         gifaddframe(animation)
    
54.         nn += 1
    
55. 
    
56. #Finish gif animation
    
57. animation.finish()
    
58. print 'Finished...'

#Set data folders<br /> obsdir = 'U:/data/micaps'<br /> mapdir = 'T:/verification/map'<br /> outdir = 'C:/Temp/test'<br /> <br /> #Set time<br /> st = datetime.datetime(2015,11,1,0)<br /> sbjt = st + datetime.timedelta(hours=8)<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 /> #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 /> nn = 0<br /> bjt = sbjt<br /> for t in range(1, 25):<br />     if nn > 0:<br />         cll()<br />     bjt = bjt + datetime.timedelta(hours=3)<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 />         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 />         #Add frame to gif animation<br />         gifaddframe(animation)<br />         nn += 1<br /> <br /> #Finish gif animation<br /> animation.finish()<br /> print 'Finished...'

cny

MeteoInfo 发表于 2016-9-6 22:46
参考此脚本：

1. basedir = 'G:\desktop\meteinfo'
   
2. obsdir =os.path.join( basedir, 'micaps/')
   
3. mapdir = os.path.join(basedir, 'MeteoInfo/map')
   
4. outdir = 'G:/desktop/'
   
5. st = datetime.datetime(2015,2,1,2)
   
6. sbjt = st + datetime.timedelta(hours=0)
   
7. #Plot
   
8. figure(figsize=[768,480], newfig=False)
   
9. axesm(position=[0,0,1,1], xyscale=1.2, tickfontsize=12)
   
10. lworld = shaperead(os.path.join(mapdir, 'country1.shp'))
    
11. lchina = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))
    
12. geoshow(lchina)
    
13. geoshow(lworld, edgecolor='k')
    
14. 
    
15. #Add south China Sea
    
16. lscs = shaperead(os.path.join(mapdir, 'bou1_4l.shp'))
    
17. axesm(position=[0.76,0.09,0.16,0.22], axison=False, xyscale=1.2)
    
18. geoshow(lscs, facecolor=(0,0,255))
    
19. xlim(106, 123)
    
20. ylim(2, 23)
    
21. #Set weather list - haze, mist and fog
    
22. weathers = [6,7,31]
    
23. ls = weatherspec(weathers)
    
24. #Set current plot to 1
    
25. currentplot(1)
    
26. giffn = os.path.join(outdir, 'V_vis_' + st.strftime('%Y%m%d') + '--loop-.gif')
    
27. print giffn
    
28. animation = gifanimation(giffn,delay=500)
    
29. nn=0
    
30. bjt= sbjt
    
31. for t in range(1,50):
    
32.     if nn > 0:
    
33.         cll()
    
34.     bjt = bjt + datetime.timedelta(hours=3)
    
35.     obsfn = os.path.join(obsdir,  str(bjt.year) + '/plot/' + bjt.strftime('%y%m%d%H.000'))
    
36.     if os.path.exists(obsfn):
    
37.         print 'Observation data file: ' + obsfn   
    
38.         obsf = addfile_micaps(obsfn)      
    
39.         odata = obsf.stationdata('WeatherNow')
    
40.         lobs = scatterm(odata, symbolspec=ls)
    
41.         title('Dust Event ' + sbjt.strftime('%Y-%m-%d %H:00') + \
    
42.                 ' +' + str(t*3) + ' (' + bjt.strftime('%Y-%m-%d %H:00') + ')', \
    
43.                 bold=False)
    
44.         xlim(68.6, 140.5)
    
45.         ylim(17.3, 54)
    
46.         #colorbar(lfcst, extendrect=False, shrink=0.6)
    
47.         #Add frame to gif animation
    
48.         gifaddframe(animation)
    
49.         nn += 1
    
50. animation.finish()
    
51. print 'Finished...'

复制代码

出来的图跟上一张一样，weather symbol 只是不停叠加，最后就越来越多，越来越密。

不知道错在哪？

MeteoInfo

cny 发表于 2016-9-7 09:58
出来的图跟上一张一样，weather symbol 只是不停叠加，最后就越来越多，越来越密。

不知道错在 ...

你用的什么版本的MeteoInfo？

cny

MeteoInfo 发表于 2016-9-7 10:15
你用的什么版本的MeteoInfo？

Java里的meteoinfo lab跑出来的，版本号是0.8.2

MeteoInfo

cny 发表于 2016-9-7 10:39
Java里的meteoinfo lab跑出来的，版本号是0.8.2

请下载最新版本试试

cny

MeteoInfo 发表于 2016-9-7 11:38
请下载最新版本试试

按您说的下载版本后能跑了，但是跑到nn=137时，提示出现错误，这可怎么办

MeteoInfo

cny 发表于 2016-9-7 19:18
按您说的下载版本后能跑了，但是跑到nn=137时，提示出现错误，这可怎么办

应该是15031714.000数据文件有问题，里面有".00**********"，需要进行修改。

Finnick

MeteoInfo 发表于 2016-9-7 19:34
应该是15031714.000数据文件有问题，里面有".00**********"，需要进行修改。

王老师，参考了汇总帖里王老师的程序，但是出来的gif动图是模糊的，想知道怎么能保存比较清楚的图片啊