MeteoInfo脚本示例：MICAPS第4类数据转为NetCDF

njzqxt

王老师，能不能把10楼转换的 py 脚本文件贴上来，我对照一下呢？谢谢了！

MeteoInfo

njzqxt 发表于 2012-12-19 11:33

                               
登录/注册后可看大图

王老师，能不能把10楼转换的 py 脚本文件贴上来，我对照一下呢？谢谢了！

1. #--------------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2012-12-11
   
4. # Purpose: Convert MICAPS 4 grid data to NetCDF data format
   
5. # Note: Sample
   
6. #-----------------------------------------------------------
   
7. import clr
   
8. clr.AddReferenceByPartialName("System.Windows.Forms")
   
9. clr.AddReferenceByPartialName("System.Drawing")
   
10. from System import *
    
11. from System.Windows.Forms import *
    
12. from System.Drawing import *
    
13. clr.AddReference("MeteoInfoC.dll")
    
14. from MeteoInfoC import *
    
15. from MeteoInfoC.Data import *
    
16. from MeteoInfoC.Data.MeteoData import *
    
17. 
    
18. import os.path
    
19. import datetime
    
20. import time
    
21. 
    
22. #---- Set data folder
    
23. dataDir = "D:\\Temp\\micaps\\ttadv\"
    
24. 
    
25. #---- Open a MICAPS 4 Data file to get the information
    
26. inDataInfo = MeteoDataInfo()
    
27. infile = dataDir + "1000\\12121408.000"
    
28. inDataInfo.OpenMICAPSData(infile)
    
29. undef = inDataInfo.UNDEF
    
30. lons = inDataInfo.GetX()
    
31. lats = inDataInfo.GetY()
    
32. lonNum = lons.Length
    
33. latNum = lats.Length
    
34. ls = [1000, 925, 850, 700, 500, 400, 300, 200, 150, 100]
    
35. levelNum = len(ls)
    
36. levels = Array.CreateInstance(Double, levelNum)
    
37. i = 0
    
38. for l in ls:
    
39.         levels[i] = l;
    
40.         i += 1
    
41. 
    
42. #---- Create netCDF data info object
    
43. outfilepath = "D:\" + "test.nc"
    
44. outDataInfo = NetCDFDataInfo()
    
45. outDataInfo.fileName = outfilepath
    
46. outDataInfo.UNDEF = undef
    
47. outDataInfo.unlimdimid = 3
    
48. 
    
49. #---- Add dimensions: lon, lat, level, time
    
50. lonDim = outDataInfo.AddDimension("lon", lonNum)
    
51. latDim = outDataInfo.AddDimension("lat", latNum)
    
52. levelDim = outDataInfo.AddDimension("level", levelNum)
    
53. timeDim = outDataInfo.AddDimension("time", -1)
    
54. 
    
55. #---- Add variables
    
56. outDataInfo.AddVariable("lon", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([lonDim]))
    
57. outDataInfo.AddVariable("lat", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([latDim]))
    
58. outDataInfo.AddVariable("level", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([levelDim]))
    
59. outDataInfo.AddVariable("time", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([timeDim]))
    
60. outDataInfo.AddVariable("temadv", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([timeDim, levelDim, latDim, lonDim]))
    
61. #---- Add variable attributes
    
62. outDataInfo.AddVariableAttribute("lon", "units", "degrees_east")
    
63. outDataInfo.AddVariableAttribute("lon", "long_name", "longitude")
    
64. outDataInfo.AddVariableAttribute("lat", "units", "degrees_north")
    
65. outDataInfo.AddVariableAttribute("lat", "long_name", "latitude")
    
66. outDataInfo.AddVariableAttribute("level", "units", "level");
    
67. outDataInfo.AddVariableAttribute("level", "long_name", "level");
    
68. outDataInfo.AddVariableAttribute("level", "axis", "z");
    
69. outDataInfo.AddVariableAttribute("time", "units", "hours since 1800-1-1 00:00:00")
    
70. outDataInfo.AddVariableAttribute("time", "long_name", "time")
    
71. outDataInfo.AddVariableAttribute("temadv", "units", "No")
    
72. outDataInfo.AddVariableAttribute("temadv", "long_name", "TemperatureAdvection")
    
73. outDataInfo.AddVariableAttribute("temadv", "missing_value", undef)
    
74. 
    
75. #---- Add global attributes
    
76. outDataInfo.AddGlobalAttribute("title", "Temperature advection field data")
    
77. outDataInfo.AddGlobalAttribute("description", "script sample data")
    
78. 
    
79. #---- Create netCDF file
    
80. outDataInfo.CreateNCFile(outfilepath)
    
81. 
    
82. #---- Write lon, lat, level, time data
    
83. outDataInfo.WriteVar("lon", lons)
    
84. outDataInfo.WriteVar("lat", lats)
    
85. outDataInfo.WriteVar("level", levels)
    
86. 
    
87. timeArray = Array.CreateInstance(Double, 8)
    
88. times = []
    
89. st = datetime.datetime(1800, 1, 1, 0, 0, 0)
    
90. t = datetime.datetime(2012, 12, 12, 8, 0, 0)        
    
91. for i in range(0, 8):        
    
92.         times.append(t)
    
93.         delta = t - st
    
94.         hours = (delta.days * 86400 + delta.seconds) / 3600.0
    
95.         t = t + datetime.timedelta(hours=12)
    
96.         timeArray[i] = hours
    
97.         print timeArray[i]
    
98.         
    
99. startArray = Array.CreateInstance(int, 1)
    
100. countArray = Array.CreateInstance(int, 1)
     
101. startArray[0] = 0
     
102. countArray[0] = len(times)
     
103. outDataInfo.WriteVara("time", startArray, countArray, timeArray)
     
104. 
     
105. #---- Write data
     
106. tidx = 0
     
107. for t in times:
     
108.         lidx = 0
     
109.         for l in ls:
     
110.                 infile = dataDir + str(l) + "\" + t.strftime("%y%m%d%H") + ".000"
     
111.                 print infile
     
112.                 if os.path.isfile(infile):
     
113.                         print infile
     
114.                         inDataInfo.OpenMICAPSData(infile)
     
115.                         gdata = inDataInfo.GetGridData("var")
     
116.                         dataArray = gdata.ToOneDimData()        
     
117.                         startArray = Array.CreateInstance(int, 4)
     
118.                         countArray = Array.CreateInstance(int, 4)
     
119.                         startArray[0] = tidx
     
120.                         startArray[1] = lidx
     
121.                         startArray[2] = 0
     
122.                         startArray[3] = 0
     
123.                         countArray[0] = 1
     
124.                         countArray[1] = 1
     
125.                         countArray[2] = latDim.DimLength
     
126.                         countArray[3] = lonDim.DimLength
     
127.                         outDataInfo.WriteVara("temadv", startArray, countArray, dataArray)
     
128.                 lidx += 1
     
129.         tidx += 1
     
130. 
     
131. #---- Close netCDF file
     
132. outDataInfo.CloseNCFile()
     
133. 
     
134. print 'Finished...!'

复制代码

njzqxt

    谢谢王老师！
    更新了文件后，数据转出来了，不是人品的问题
    上张美丽的像素。

MeteoInfo

njzqxt 发表于 2012-12-20 15:32

                               
登录/注册后可看大图

谢谢王老师！
    更新了文件后，数据转出来了，不是人品的问题
    上张美丽的像素。

人品没问题，运气有点衰呀

njzqxt

不存在。遇上王老师便是最大好运。

njzqxt

    原数据分别来自1000、 925、 850、 700、 500、 400、 300、 200、 150、 100hPa，第层8个数据，共80个数据，转成一个NetCD格式数据。我用GRADS分别作了80张图，与MICAPS原图逐一对照，结果是一模一样。又作某层的时间剖面图（见13楼），也是没有问题。
    作垂直剖面图就出问题，脚本如下，'set lev 1000 100'时能出图，但纵坐标上高度却反了（如图）。
'reinit'
'sdfopen d:\test.nc'
'set t 5'
'set lat 27.5'
'set lon 85 115'
'set lev 100 1000'
'd temadv'
;
    设为'set lev 1000 100'时就出错（如图2），这是什么原因呢？

MeteoInfo

njzqxt 发表于 2012-12-20 18:47

                               
登录/注册后可看大图

原数据分别来自1000、 925、 850、 700、 500、 400、 300、 200、 150、 100hPa，第层8个数据，共80个 ...

我基本不用GrADS，你也许可以试试 set lev 1000 100

卡拉蜂空间

MeteoInfo 发表于 2012-12-18 15:38

                               
登录/注册后可看大图

时间设置错了，改为：
t = datetime.datetime(2012, 12, 12, 8, 0, 0)

麻烦问一下，这个时间设置（2012,12,12,8,0,0 ）分别表示的什么意思？2012年12月。。。然后呢？

卡拉蜂空间

卡拉蜂空间 发表于 2013-4-30 10:54

                               
登录/注册后可看大图

麻烦问一下，这个时间设置（2012,12,12,8,0,0 ）分别表示的什么意思？2012年12月。。。然后呢？

我的文件时1000 925 850 700 500五个文件夹，每个文件夹里有44个数据应该改哪些呀

dataDir = "E:\\micaps\\high-plot\\"

#---- Open a MICAPS 4 Data file to get the information
inDataInfo = MeteoDataInfo()
infile = dataDir + "1000\\12031202.000"
inDataInfo.OpenMICAPSData(infile)
undef = inDataInfo.UNDEF
lons = inDataInfo.GetX()
lats = inDataInfo.GetY()
lonNum = lons.Length
latNum = lats.Length
ls = [1000, 925, 850, 700, 500]
levelNum = len(ls)
levels = Array.CreateInstance(Double, levelNum)
i = 0
for l in ls:
        levels = l;
        i += 1

#---- Create netCDF data info object
outfilepath = dataDir + "test.nc"
outDataInfo = NetCDFDataInfo()
outDataInfo.fileName = outfilepath
outDataInfo.UNDEF = undef
outDataInfo.unlimdimid = 3

#---- Add dimensions: lon, lat, level, time
lonDim = outDataInfo.AddDimension("lon", lonNum)
latDim = outDataInfo.AddDimension("lat", latNum)
levelDim = outDataInfo.AddDimension("level", levelNum)
timeDim = outDataInfo.AddDimension("time", -1)

#---- Add variables
outDataInfo.AddVariable("lon", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([lonDim]))
outDataInfo.AddVariable("lat", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([latDim]))
outDataInfo.AddVariable("level", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([levelDim]))
outDataInfo.AddVariable("time", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([timeDim]))
outDataInfo.AddVariable("temadv", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([timeDim, levelDim, latDim, lonDim]))
#---- Add variable attributes
outDataInfo.AddVariableAttribute("lon", "units", "degrees_east")
outDataInfo.AddVariableAttribute("lon", "long_name", "longitude")
outDataInfo.AddVariableAttribute("lat", "units", "degrees_north")
outDataInfo.AddVariableAttribute("lat", "long_name", "latitude")
outDataInfo.AddVariableAttribute("level", "units", "level");
outDataInfo.AddVariableAttribute("level", "long_name", "level");
outDataInfo.AddVariableAttribute("level", "axis", "z");
outDataInfo.AddVariableAttribute("time", "units", "hours since 1800-1-1 00:00:00")
outDataInfo.AddVariableAttribute("time", "long_name", "time")
outDataInfo.AddVariableAttribute("temadv", "units", "No")
outDataInfo.AddVariableAttribute("temadv", "long_name", "TemperatureAdvection")
outDataInfo.AddVariableAttribute("temadv", "missing_value", undef)

#---- Add global attributes
outDataInfo.AddGlobalAttribute("title", "Temperature advection field data")
outDataInfo.AddGlobalAttribute("description", "script sample data")

#---- Create netCDF file
outDataInfo.CreateNCFile(outfilepath)

#---- Write lon, lat, level, time data
outDataInfo.WriteVar("lon", lons)
outDataInfo.WriteVar("lat", lats)
outDataInfo.WriteVar("level", levels)

timeArray = Array.CreateInstance(Double, 8)
times = []
st = datetime.datetime(1800, 1, 1, 0, 0, 0)
t = datetime.datetime(2012, 12, 8, 8, 0, 0)        
for i in range(0, 8):        
        times.append(t)
        delta = t - st
        hours = (delta.days * 86400 + delta.seconds) / 3600.0
        t = t + datetime.timedelta(hours=12)
        timeArray = hours
        print timeArray
        
startArray = Array.CreateInstance(int, 1)
countArray = Array.CreateInstance(int, 1)
startArray[0] = 0
countArray[0] = len(times)
outDataInfo.WriteVara("time", startArray, countArray, timeArray)

#---- Write data
tidx = 0
for t in times:
        lidx = 0
        for l in ls:
                infile = dataDir + str(l) + "\\" + t.strftime("%y%m%d%H") + ".000"
                #print infile
                if os.path.isfile(infile):
                        print infile
                        inDataInfo.OpenMICAPSData(infile)
                        gdata = inDataInfo.GetGridData("var")
                        dataArray = gdata.ToOneDimData()        
                        startArray = Array.CreateInstance(int, 4)
                        countArray = Array.CreateInstance(int, 4)
                        startArray[0] = tidx
                        startArray[1] = lidx
                        startArray[2] = 0
                        startArray[3] = 0
                        countArray[0] = 1
                        countArray[1] = 1
                        countArray[2] = latDim.DimLength
                        countArray[3] = lonDim.DimLength
                        outDataInfo.WriteVara("temadv", startArray, countArray, dataArray)
                lidx += 1
        tidx += 1

#---- Close netCDF file
outDataInfo.CloseNCFile()

mofangbao

卡拉蜂空间 发表于 2013-4-30 10:54

                               
登录/注册后可看大图

麻烦问一下，这个时间设置（2012,12,12,8,0,0 ）分别表示的什么意思？2012年12月。。。然后呢？

这种不应该来问的吧 不属于MeteoInfo的内容了 是基础的知识