数据合并

shzljg

我在合并NCEP 1.0*1.0数据的时候遇到了问题。具体看下图。有哪位高手能指导一下。数据是fnl文件。每个时次一个文件。

MeteoInfo

这几天比较忙，一直没顾上测试你的数据，等有空了会仔细看看。

shzljg

好的。麻烦版主了。

MeteoInfo

shzljg 发表于 2012-6-19 11:56

                               
登录/注册后可看大图

好的。麻烦版主了。

改进了netCDF数据合并功能，可以处理这种情况了，需要MeteoInfo最新文件（见置顶帖子）。需要注意的是该数据中的时间维名称并非time（通常情况下都是这个名字），而是initial_time0_hours，所以在合并文件的界面中需要将Time Dimension Name改为正确的时间维名称（如下图所示）。


也可以通过写脚本程序来合并netCDF文件，例如：

1. #--------------------------------------------------------        
   
2. # Author: Yaqiang Wang                                          
   
3. # Date: 2012-5-22                                                
   
4. # Purpose: Join NetCDF data files  
   
5. # Note: Sample                                                   
   
6. #-----------------------------------------------------------     
   
7. import clr
   
8. from System import *
   
9. from System.Collections.Generic import List
   
10. clr.AddReference("MeteoInfoC.dll")
    
11. from MeteoInfoC import *
    
12. from MeteoInfoC.Data import *
    
13. from MeteoInfoC.Data.MeteoData import *
    
14. 
    
15. import os.path
    
16. import datetime
    
17. 
    
18. #---- Set data folders
    
19. dataDir = 'F:\\Temp\\nc\\'
    
20. 
    
21. #---- Join netCDF data files
    
22. infiles = List[str]()
    
23. aTime = datetime.datetime(2007,7,16)
    
24. dt = datetime.timedelta(hours=6)
    
25. for h in range(2):
    
26.         tStr = aTime.strftime('%Y%m%d_%H')
    
27.         afile = dataDir + 'fnl_' + tStr + '_00.nc'
    
28.         print afile
    
29.         if os.path.isfile(afile):
    
30.                 infiles.Add(afile)
    
31.         aTime = aTime + dt
    
32. 
    
33. if infiles.Count > 1:
    
34.         outfile = dataDir + 'fnl_output.nc'
    
35.         NetCDFData.JoinDataFiles_Time(infiles, outfile, 'initial_time0_hours')
    
36. print 'Finished...'

复制代码

shzljg

谢谢版主了。搞好了，我还想问一下，用最新的软件版本打开原来的fnl_20070716_00_00.nc,为什么许多变量应该为26层，但是选择变量后，level仍然为surface.比如打开u_grd_3_isbl_10.

MeteoInfo

shzljg 发表于 2012-6-24 11:06

                               
登录/注册后可看大图

谢谢版主了。搞好了，我还想问一下，用最新的软件版本打开原来的fnl_20070716_00_00.nc,为什么许多变量应该 ...

NetCDF格式本身十分灵活，要让程序自动识别需要遵循某种conventions，比如CF conventions。你的这个数据的维设定比较复杂，不同的变量可能有不同的level维，目前还没有想好怎么自动处理这样的情况，以后有空会仔细研究一下的。

MeteoInfo

shzljg 发表于 2012-6-24 11:06

                               
登录/注册后可看大图

谢谢版主了。搞好了，我还想问一下，用最新的软件版本打开原来的fnl_20070716_00_00.nc,为什么许多变量应该 ...

已经做了改进，可以处理这种情况了，详见此贴：改进NetCDF数据读取
http://bbs.06climate.com/forum.p ... 270&fromuid=106

shzljg

谢谢版主了！