王老师，求助！ 处理micaps数据中关于时间名的设置

wuwei2163

王老师，你看错了，你说的是第一个错误，我这是第二个错误啊，貌似没有定义文件名吧？？？

MeteoInfo

wuwei2163 发表于 2014-3-18 14:09

                               
登录/注册后可看大图

王老师，你看错了，你说的是第一个错误，我这是第二个错误啊，貌似没有定义文件名吧？？？

图片被气象家园的logo给盖住了，看不清楚。
fileName改为FileName

wuwei2163

王老师，在您的指导下，已经有点熟悉了，下面的这个错误是 定义数据类型时出错了，但不知道怎样改啊

望王老师教我啊

MeteoInfo

wuwei2163 发表于 2014-3-18 14:55

                               
登录/注册后可看大图

王老师，在您的指导下，已经有点熟悉了，下面的这个错误是 定义数据类型时出错了，但不知道怎样改啊

望王 ...

levels数组你是怎么定义赋值的？

wuwei2163

MeteoInfo 发表于 2014-3-18 15:04

                               
登录/注册后可看大图

levels数组你是怎么定义赋值的？

王老师，看你程序上定义的数组是Double,我levels上有7层，感觉定义整型和浮点型都可以的，但总是出错，现在被这个程序搞得头晕脑胀的，我把资料数据以及所有的代码贴出来，麻烦您看看，到底是哪里的问题？？？
# Author: Yaqiang Wang
# Date: 2012-12-11
# Purpose: Convert MICAPS 4 grid data to NetCDF data format
# Note: Sample
#-----------------------------------------------------------
import clr
clr.AddReferenceByPartialName("System.Windows.Forms")
clr.AddReferenceByPartialName("System.Drawing")
from System import *
from System.Windows.Forms import *
from System.Drawing import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data import *
from MeteoInfoC.Data.MeteoData import *

import os.path
import datetime
import time

#---- Set data folder
dataDir = "y:\\ecmwf_thin\\R\\"

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

#---- Create netCDF data info object
outfilepath = "D:\\" + "test.nc"
outDataInfo = NetCDFDataInfo()
outDataInfo.FileName = outfilepath
outDataInfo.MissingValue = 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, 7)
times = []
st = datetime.datetime(1800, 1, 1, 0, 0, 0)
t = datetime.datetime(2014, 03, 17, 8, 0, 0)        
for i in range(0, 9):        
        times.append(t)
        delta = t - st
        hours = (delta.days * 86400 + delta.seconds) / 3600.0
        t = t + datetime.timedelta(hours=3)
        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()

print 'Finished...!'

MeteoInfo

wuwei2163 发表于 2014-3-18 16:36

                               
登录/注册后可看大图

王老师，看你程序上定义的数组是Double,我levels上有7层，感觉定义整型和浮点型都可以的，但总是出错，现 ...

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

#--------------------------------------------------------# Author: Yaqiang Wang<br /> # Date: 2014-3-18<br /> # Purpose: Convert MICAPS 4 grid data to NetCDF data format<br /> # Note: Sample<br /> #-----------------------------------------------------------<br /> import clr<br /> clr.AddReferenceByPartialName("System.Windows.Forms")<br /> clr.AddReferenceByPartialName("System.Drawing")<br /> from System import *<br /> from System.Windows.Forms import *<br /> from System.Drawing import *<br /> clr.AddReference("MeteoInfoC.dll")<br /> from MeteoInfoC import *<br /> from MeteoInfoC.Data import *<br /> from MeteoInfoC.Data.MeteoData import *<br /> <br /> import os.path<br /> import datetime<br /> import time<br /> <br /> #---- Set data folder<br /> dataDir = "D:\\Temp\\test\\R\\"<br /> <br /> #---- Open a MICAPS 4 Data file to get the information<br /> inDataInfo = MeteoDataInfo()<br /> infile = dataDir + "1000\\14031708.000.000"<br /> inDataInfo.OpenMICAPSData(infile)<br /> undef = inDataInfo.MissingValue<br /> lons = inDataInfo.GetX()<br /> lats = inDataInfo.GetY()<br /> lonNum = lons.Length<br /> latNum = lats.Length<br /> ls = [1000, 925, 850, 700, 500, 300, 100]<br /> levelNum = len(ls)<br /> levels = Array.CreateInstance(Double, levelNum)<br /> i = 0<br /> for l in ls:<br />         levels<i> = l;<br />         i += 1<br /> <br /> #---- Create netCDF data info object<br /> outfilepath = dataDir + "test.nc"<br /> outDataInfo = NetCDFDataInfo()<br /> outDataInfo.FileName = outfilepath<br /> outDataInfo.MissingValue = undef<br /> outDataInfo.unlimdimid = 3<br /> <br /> #---- Add dimensions: lon, lat, level, time<br /> lonDim = outDataInfo.AddDimension("lon", lonNum)<br /> latDim = outDataInfo.AddDimension("lat", latNum)<br /> levelDim = outDataInfo.AddDimension("level", levelNum)<br /> timeDim = outDataInfo.AddDimension("time", -1)<br /> <br /> #---- Add variables<br /> outDataInfo.AddVariable("lon", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([lonDim]))<br /> outDataInfo.AddVariable("lat", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([latDim]))<br /> outDataInfo.AddVariable("level", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([levelDim]))<br /> outDataInfo.AddVariable("time", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([timeDim]))<br /> outDataInfo.AddVariable("temadv", NetCDF4.NcType.NC_DOUBLE, Array[Dimension]([timeDim, levelDim, latDim, lonDim]))<br /> #---- Add variable attributes<br /> outDataInfo.AddVariableAttribute("lon", "units", "degrees_east")<br /> outDataInfo.AddVariableAttribute("lon", "long_name", "longitude")<br /> outDataInfo.AddVariableAttribute("lat", "units", "degrees_north")<br /> outDataInfo.AddVariableAttribute("lat", "long_name", "latitude")<br /> outDataInfo.AddVariableAttribute("level", "units", "level");<br /> outDataInfo.AddVariableAttribute("level", "long_name", "level");<br /> outDataInfo.AddVariableAttribute("level", "axis", "z");<br /> outDataInfo.AddVariableAttribute("time", "units", "hours since 1800-1-1 00:00:00")<br /> outDataInfo.AddVariableAttribute("time", "long_name", "time")<br /> outDataInfo.AddVariableAttribute("temadv", "units", "No")<br /> outDataInfo.AddVariableAttribute("temadv", "long_name", "TemperatureAdvection")<br /> outDataInfo.AddVariableAttribute("temadv", "missing_value", undef)<br /> <br /> #---- Add global attributes<br /> outDataInfo.AddGlobalAttribute("title", "Temperature advection field data")<br /> outDataInfo.AddGlobalAttribute("description", "script sample data")<br /> <br /> #---- Create netCDF file<br /> outDataInfo.CreateNCFile(outfilepath)<br /> <br /> #---- Write lon, lat, level, time data<br /> outDataInfo.WriteVar("lon", lons)<br /> outDataInfo.WriteVar("lat", lats)<br /> outDataInfo.WriteVar("level", levels)<br /> <br /> tnum = 6<br /> timeArray = Array.CreateInstance(Double, tnum)<br /> times = []<br /> st = datetime.datetime(1800, 1, 1, 0, 0, 0)<br /> t = datetime.datetime(2014, 3, 17, 8, 0, 0)        <br /> for i in range(0, tnum):        <br />         times.append(t)<br />         delta = t - st<br />         hours = (delta.days * 86400 + delta.seconds) / 3600.0<br />         t = t + datetime.timedelta(hours=3)<br />         timeArray<i> = hours<br />         print timeArray<i><br />         <br /> startArray = Array.CreateInstance(int, 1)<br /> countArray = Array.CreateInstance(int, 1)<br /> startArray[0] = 0<br /> countArray[0] = len(times)<br /> outDataInfo.WriteVara("time", startArray, countArray, timeArray)<br /> <br /> #---- Write data<br /> tidx = 0<br /> st = times[0]<br /> for t in times:<br />         lidx = 0        <br />         for l in ls:<br />                 tt = tidx * 3<br />                 tstr = '.%03d'%tt<br />                 infile = dataDir + str(l) + '\\' + st.strftime("%y%m%d%H") + tstr + tstr<br />                 #print infile<br />                 if os.path.isfile(infile):<br />                         print infile<br />                         inDataInfo.OpenMICAPSData(infile)<br />                         gdata = inDataInfo.GetGridData("var")<br />                         dataArray = gdata.ToOneDimData()        <br />                         startArray = Array.CreateInstance(int, 4)<br />                         countArray = Array.CreateInstance(int, 4)<br />                         startArray[0] = tidx<br />                         startArray[1] = lidx<br />                         startArray[2] = 0<br />                         startArray[3] = 0<br />                         countArray[0] = 1<br />                         countArray[1] = 1<br />                         countArray[2] = latDim.DimLength<br />                         countArray[3] = lonDim.DimLength<br />                         outDataInfo.WriteVara("temadv", startArray, countArray, dataArray)<br />                 lidx += 1<br />         tidx += 1<br /> <br /> #---- Close netCDF file<br /> outDataInfo.CloseNCFile()<br /> <br /> print 'Finished...!'

wuwei2163

MeteoInfo 发表于 2014-3-18 17:15

                               
登录/注册后可看大图

王老师，谢谢您的耐心教导，根据你修改后的程序，可以运行处结果，但是出现了这样的两个错误：
第一个是 ，点击确定后，出现第二个 ，运行出的结果可以看出时间不对，而且结果大小只有2k多.......
论坛上也有人遇到这样的错误，更新了文件后据说可以运行结果，我先是更新个最新的MeteoInfoC.dll文件，结果还是同样的错误，然后我又卸掉了原来的版本，在论坛上下载了最新的1.1.3.0版本，还是不行的，王老师你看这是什么情况？？？