绘制micaps 1类数据出现问题

rhythm9806

我用版主给的script绘制累积雨量，用的数据时MICAPS 1类数据，但是出现问题；
我上传了几个我的资料文件，是国家局FTP上下的，好像和我平时用的稍微有点不一样
版主帮忙看看，问题在哪里
09022417.000 (375.35 KB, 下载次数: 7)

2011-7-8 11:13 上传

点击文件名下载附件

• 关闭
• 附件列表
• 批量上传
  

mofangbao

还不会用python，我out了...

xuminge

micaps3里可以画累积雨量

rhythm9806

我从单位FTP上下的资料是可以用的，我试过，但是没有我要的时间段，所以从国家局下了，但是使用起来就有上面的错误

MeteoInfo

MeteoInfoC.dll中添加标注功能中有一个bug，下载附件并覆盖原来的MeteoInfoC.dll即可。

MeteoInfoC.dll (1.57 MB, 下载次数: 21)

2011-7-8 23:24 上传

点击文件名下载附件

另外MICAPS第一类数据中的6小时降水只在02、08、14、20时有数据（间隔6小时），其它时次所有站点的6小时降水均为无效值，因此在累计降水中只能用02、08、14、20时数据。在MeteoInfo中两个站点数据相加，会通过判断站点号来一一对应进行相加，如果同一站点中有一个数据为无效值，则该站点的累加值也为无效值。

rhythm9806

还是有问题，运行错误信息

MeteoInfo

我这里运行没有问题。

你把程序都贴出来吧，我看看。

rhythm9806

# This Python file uses the following encoding: GB2312
#---- 引入类库
import clr
clr.AddReferenceByPartialName("System")
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.Layout import *
from MeteoInfoC.Data import *

#---- 创建MIApp类的对象
myApp = MIApp()
#---- 添加图层
myApp.OpenLayer(".\\Map\\bou2_4p.shp")
#---- 设置图层渲染
myApp.SetLegendBreak("bou2_4p.shp",0,Color.Yellow,Color.Gray,1,True,False,True)
myApp.OpenLayer(".\\Map\\bou1_4l.shp")
myApp.SetLegendBreak("bou1_4l.shp",0,Color.Blue)
myApp.OpenLayer(".\\Map\\china.shp")
myApp.SetLayerVisible("china.shp", False)
myApp.OpenLayer(".\\Map\\res1_4m.shp")
myApp.SetLegendBreak("res1_4m.shp",0,5,Color.Red,Color.Black,True,True,True)
#---- 获取矢量图层对象，并设置标注
aLayer = myApp.GetVectorLayer("res1_4m.shp")
aLayer.LabelSet.FieldName = "NAME";
aLayer.LabelSet.LabelFont = Font("楷体", 10);
aLayer.LabelSet.Offset = 15;
aLayer.AddLabels();

#---- Lambert投影
myApp.ProjectLayers("+proj=lcc+lat_1=25+lat_2=47+lon_0=105")

#---- 按照经纬度范围缩放地图
myApp.ZoomLonLat(100,125,20,35)

#---- 设置屏蔽图层（只绘制中国境内图形）
myApp.SetMaskout("china.shp")
#---- 设置绘图类型为shaded
myApp.SetDrawType("shaded")
#---- 站点数据插值为格点数据的设置
myApp.SetInterpolation(60,140,-20,60,160,160,"IDW_Radius",1,1)
#---- 设置图例文件
myApp.SetLegendScheme(".\\Legend\\rain1.lgs")

#---- 设置DefalutLayoutMap（图层显示）
myApp.MapLayout.DefaultLayoutMap.DrawGridLine = False
myApp.MapLayout.DefaultLayoutMap.DrawNeatLine = False
myApp.MapLayout.DefaultLayoutMap.Left = 10
myApp.MapLayout.DefaultLayoutMap.Top = 10
myApp.MapLayout.DefaultLayoutMap.Width = 620
myApp.MapLayout.DefaultLayoutMap.Height = 450

#---- 设置DefaultIllustration（插图，这里是南海脚图）
#myApp.MapLayout.DefaultIllustration.Visible = True
#myApp.MapLayout.DefaultIllustration.Left = 40
#myApp.MapLayout.DefaultIllustration.Top = 350
#myApp.MapLayout.DefaultIllustration.Width = 85
#myApp.MapLayout.DefaultIllustration.Height = 109

#---- 设置DefaultTitle（图形标题）
myApp.MapLayout.DefaultTitle.SetFont("黑体", 12)
myApp.MapLayout.DefaultTitle.Top = 40
myApp.MapLayout.DefaultTitle.Left = 80

#---- 设置DefaultLegend（图例）
myApp.MapLayout.DefaultLegend.LegendStyle = LegendStyleEnum.Normal        
myApp.MapLayout.DefaultLegend.DrawNeatLine = False
myApp.MapLayout.DefaultLegend.Title = "降水量（毫米）"
myApp.MapLayout.DefaultLegend.BackColor = Color.White
myApp.MapLayout.DefaultLegend.Font = Font("Arial", 8)
myApp.MapLayout.DefaultLegend.Left = 575
myApp.MapLayout.DefaultLegend.Top = 250

#---- 设置起始结束时间
sTime = DateTime.Parse("2009-02-14 08:00")
eTime = DateTime.Parse("2009-03-08 20:00")
aTime = sTime.AddHours(-6)

stData = StationData()
i = 1
#---- 循环
while sTime <= eTime:
  myApp.RemoveDataLayers()
  inFile = sTime.ToString("yyMMddHH") + ".000"
  print inFile
  #---- 打开MICAPS数据文件
  myApp.OpenMICAPSData(".\\Sample\\MICAPS\\" + inFile)  
  #---- 获取累计降水站点数据
  if i == 1:
      stData = myApp.GetStationData("Precipitation6h")
  else:
      stData = stData + myApp.GetStationData("Precipitation6h")

  #---- 时间加6小时
  sTime = sTime.AddHours(6)
  i = i + 1

#---- 绘制累计降水图层
myApp.Display(stData)
#---- 调整图层顺序（以避免压盖）
myApp.MoveLayerToTop("bou2_4p.shp")  
myApp.MoveLayerToTop("bou1_4l.shp")
myApp.MoveLayerToTop("res1_4m.shp")
#---- 设置标题名称
title = "            全国降水量实况图\n" + aTime.ToString("yyyy年MM月dd日HH时") + " - " + eTime.ToString("yyyy年MM月dd日HH时")
myApp.SetTitle(title)     
#---- 设置图例名称
myApp.MapLayout.DefaultLegend.Title = "降水量（毫米）"
#---- 绘制图形
myApp.MapLayout.PaintGraphics()
#---- 输出图形为文件
outFile = "Prec_" + sTime.ToString("yyyyMMddHH") + "-" + eTime.ToString("yyyyMMddHH") + ".png"
print outFile
myApp.SaveFigure("G:\\" + outFile)   

rhythm9806

就是你的脚本，修改了一下时间

rhythm9806

绘制每日的降水也出现问题，