脚本如何实现MICAPS第一类多要素的叠加显示？【已解决】

孤蓝et · 发表于 2011-11-29 21:08:50

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由孤蓝et 于 2011-12-2 11:24 编辑

1、备注：meteoinfo脚本中对变量名是大小写敏感的，这个大家要注意哦。 myApp.Display("Temperature")
2、问题：meteoinfo自带的Script示例文件只有针对单一变量的绘图，如何实现多个变量叠加显示？
例如：micaps第一类数据，实现温度和风场的叠加显示？这个怎么实现哦？目前我只会用shade之类的绘制，barb绘制风场还没有实现。在baidu博客的meteoinfo相关script中也没有找到相关方法，望斑竹赐教哦？万分感谢。

秀个图：现在需要叠加风场。。。。求方法哟。

孤蓝et · 发表于 2011-11-29 21:11:38

上图绘制的代码，基于meteoinfo目录下的Script提供的示例修改得到如下代码。

# 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 *
#---- 创建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=84")
#---- 按照经纬度范围缩放地图
myApp.ZoomLonLat(76,94,35,50)
#---- 设置屏蔽图层（只绘制中国境内图形）
myApp.SetMaskout("china.shp")
#---- 设置绘图类型为shaded
myApp.SetDrawType("shaded")
#---- 站点数据插值为格点数据的设置
myApp.SetInterpolation(60,140,-20,60,160,160,"IDW_Radius",1,1)
#---- 设置图例文件
myApp.SetLegendScheme("..\\Legend\\tem.lgs")
#---- 设置DefalutLayoutMap（图层显示）
myApp.MapLayout.DefaultLayoutMap.DrawGridLine = True
myApp.MapLayout.DefaultLayoutMap.DrawNeatLine = True
myApp.MapLayout.DefaultLayoutMap.Left = 10
myApp.MapLayout.DefaultLayoutMap.Top = 10
myApp.MapLayout.DefaultLayoutMap.Width = 620
myApp.MapLayout.DefaultLayoutMap.Height = 450
#---- 设置DefaultTitle（图形标题）
myApp.MapLayout.DefaultTitle.SetFont("黑体", 12)
myApp.MapLayout.DefaultTitle.Top = 30
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 = 640
myApp.MapLayout.DefaultLegend.Top = 150
#---- 设置起始结束时间
sTime = DateTime.Parse("2011-8-3 00:00")
eTime = DateTime.Parse("2011-8-3 23:00")
#---- 循环
while sTime <= eTime:
myApp.RemoveDataLayers()
inFile = sTime.ToString("yyMMddHH") + ".000"
print inFile
#---- 打开MICAPS数据文件
myApp.OpenMICAPSData("..\\Sample\\aws\" + inFile)
#---- 生成6小时降水量图层
myApp.Display("Temperature")
print "Display finished"
#---- 调整图层顺序（以避免压盖）
myApp.MoveLayerToTop("bou2_4p.shp")
#myApp.MoveLayerToTop("bou1_4l.shp")
myApp.MoveLayerToTop("res1_4m.shp")
#---- 设置标题名称
title = " 新疆区域温度分布图 (" + sTime.ToString("yyyy-MM-dd HH:00") + ")"
myApp.SetTitle(title)
#---- 设置图例名称
myApp.MapLayout.DefaultLegend.Title = "温度(℃)"
#---- 绘制图形
myApp.MapLayout.PaintGraphics()
#---- 输出图形为文件
outFile = "tem_" + sTime.ToString("yyyyMMddHH") + ".png"
print outFile
myApp.SaveFigure("E:\" + outFile)
#---- 删除6小时降水量图层
if sTime < eTime:
myApp.RemoveLastLayer()
#---- 时间加6小时
sTime = sTime.AddHours(1)

复制代码

MeteoInfo · 发表于 2011-11-29 21:23:09

参考如下代码：

myApp.SetDrawType("barb")
myApp.Display("U","V")

孤蓝et · 发表于 2011-11-29 22:25:10

本帖最后由孤蓝et 于 2011-11-29 22:27 编辑

Micaps第一类数据没有UV分量，所以我打算通过计算得到，然后使用Barb。

# 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 *
.... .....
w_dir = StationData()
w_spd = StationData()

u = StationData()
v = StationData()

w_dir = myApp.GetStationData("WindDirection")
w_spd = myApp.GetStationData("WindSpeed")

u=w_spd*DataMath.sin(w_dir)
v=w_spd*DataMath.cos(w_dir)
.... .....

然后运行脚本提示错误：看提示，貌似stationdata不支持sin函数？
C:\Program Files\MeteoInfo\Script>miscript run_tem.py
11080318.000
Traceback (most recent call last):
File "run_tem.py", line 97, in <module>
AttributeError: 'type' object has no attribute 'sin'

尽头的尽头 · 发表于 2011-11-30 00:19:46

帮顶~~~

MeteoInfo · 发表于 2011-11-30 08:23:56

孤蓝et 发表于 2011-11-29 22:25

登录/注册后可看大图

Micaps第一类数据没有UV分量，所以我打算通过计算得到，然后使用Barb。

# This Python file uses the fo ...

正如你所说，函数名是大小写敏感的，应该是DataMath.Sin , DataMath.Cos

孤蓝et · 发表于 2011-11-30 23:43:15

# 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=84")
#---- 按照经纬度范围缩放地图
myApp.ZoomLonLat(76,94,35,50)
#---- 设置屏蔽图层（只绘制中国境内图形）
myApp.SetMaskout("china.shp")
#---- 站点数据插值为格点数据的设置
myApp.SetInterpolation(60,140,-20,60,160,160,"IDW_Radius",1,1)
#---- 设置图例文件
#myApp.SetLegendScheme("..\\Legend\\tem.lgs")
#---- 设置DefalutLayoutMap（图层显示）
myApp.MapLayout.DefaultLayoutMap.DrawGridLine = True
myApp.MapLayout.DefaultLayoutMap.DrawNeatLine = True
myApp.MapLayout.DefaultLayoutMap.Left = 10
myApp.MapLayout.DefaultLayoutMap.Top = 10
myApp.MapLayout.DefaultLayoutMap.Width = 620
myApp.MapLayout.DefaultLayoutMap.Height = 450
#---- 设置DefaultTitle（图形标题）
myApp.MapLayout.DefaultTitle.SetFont("黑体", 12)
myApp.MapLayout.DefaultTitle.Top = 30
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 = 640
myApp.MapLayout.DefaultLegend.Top = 150
#---- 设置起始结束时间
sTime = DateTime.Parse("2011-8-3 18:00")
eTime = DateTime.Parse("2011-8-3 18:00")
#---- 循环
while sTime <= eTime:
myApp.RemoveDataLayers()
inFile = sTime.ToString("yyMMddHH") + ".000"
print inFile
#---- 打开MICAPS数据文件
myApp.OpenMICAPSData("..\\Sample\\aws\" + inFile)
#---- 设置绘图类型为shaded
#myApp.SetDrawType("Shaded")
#myApp.Display("Temperature")
[color=Red] w_dir = myApp.GetStationData("WindDirection")
w_spd = myApp.GetStationData("WindSpeed")
u=w_spd*DataMath.Sin(w_dir)
v=w_spd*DataMath.Cos(w_dir)
myApp.SetDrawType("Barb")
myApp.Display(u,v)[/color]
print "Display finished"
#---- 调整图层顺序（以避免压盖）
myApp.MoveLayerToTop("bou2_4p.shp")
#myApp.MoveLayerToTop("bou1_4l.shp")
myApp.MoveLayerToTop("res1_4m.shp")
#---- 设置标题名称
title = " 新疆区域温度分布图 (" + sTime.ToString("yyyy-MM-dd HH:00") + ")"
myApp.SetTitle(title)
#---- 设置图例名称
myApp.MapLayout.DefaultLegend.Title = "温度(℃)"
#---- 绘制图形
myApp.MapLayout.PaintGraphics()
#---- 输出图形为文件
outFile = "tem_" + sTime.ToString("yyyyMMddHH") + ".png"
print outFile
myApp.SaveFigure("E:\" + outFile)
#---- 删除6小时降水量图层
if sTime < eTime:
myApp.RemoveLastLayer()
#---- 时间加6小时
sTime = sTime.AddHours(1)

复制代码

为嘛画个简单的barb也这么波折？本想能很轻松的解决，现在又要劳烦斑竹了。唉，不好意思啊。
执行上述脚本报错如下：
C:\Program Files\MeteoInfo\Script>miscript run_tem.py
11080318.000
Traceback (most recent call last):
File "run_tem.py", line 97, in <module>
TypeError: expected str, got StationData

孤蓝et · 发表于 2011-11-30 23:44:56

斑竹要不要直接帮我写好一个例子用micaps第一类数据绘制风场barb撒。。。我被这个纠结了两天了。
愧于使用script不精啊。。。

{:eb343:}{:eb343:}{:eb343:}

MeteoInfo · 发表于 2011-12-1 08:57:48

孤蓝et 发表于 2011-11-30 23:44

登录/注册后可看大图

斑竹要不要直接帮我写好一个例子用micaps第一类数据绘制风场barb撒。。。我被这个纠结了两天了。
愧于使用 ...

这个不是你的问题，我正在更新MIApp类来解决这个问题。

孤蓝et · 发表于 2011-12-1 18:10:29

MeteoInfo 发表于 2011-12-1 08:57

登录/注册后可看大图

这个不是你的问题，我正在更新MIApp类来解决这个问题。

辛苦斑竹了。。。。
继续支持中。

		自动登录	找回密码
密码			立即注册