- 积分
- 1502
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-4-12
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
不知大家有没有用meteoinfo画图,在运行的时候出现这个错误,求指导
#-------------------------------------------------------- # Author: Yaqiang Wang
# Date: 2014-3-25
# Purpose: Calculate average grid data and plot
# Note: Sample
#-----------------------------------------------------------
import clr
from System.Collections.Generic import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data import *
from MeteoInfoC.Data.MeteoData import *
from MeteoInfoC.Legend import *
from MeteoInfoC.Shape import *
#---- Set directories
dataDir = "D:\\MeteoInfo\\lizi\\"
#---- Open data file
mdi = MeteoDataInfo()
mdi.OpenNCData(dataDir + "d03")
#---- Average grid data
tnum = mdi.GetTimeNumber()
print 'Time number: ' + str(tnum)
sumData = GridData()
mdi.LevelIndex = 15
varName = 'QGRAUP'
for i in range(0, tnum):
mdi.TimeIndex = i
gdata = mdi.GetGridData(varName)
if i == 0:
sumData = gdata
else:
sumData = sumData + gdata
aveData = sumData / tnum
#---- Create shaded layer from the grid data
print 'Create shaded layer'
aLS = LegendManage.CreateLegendSchemeFromGridData(aveData, LegendType.GraduatedColor, ShapeTypes.Polygon)
aLayer = DrawMeteoData.CreateShadedLayer(aveData, aLS, varName, "Value")
aLayer.ProjInfo = mdi.ProjInfo
mipy.MapDocument.ActiveMapFrame.AddLayer(aLayer)
mipy.MapDocument.ActiveMapFrame.MoveLayer(aLayer, 0)
print 'Finished!'
|
-
|