MeteoInfo脚本示例：生成可见光云图

MeteoInfo · 发表于 2013-8-27 11:07:45

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

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

x

读取MICAPS格式云图数据并生成图片。

脚本代码：

# This Python file uses the following encoding: GB2312
#--------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2013-8-27
# Purpose: Read and plot MICAPS cloud data
# Note: Sample
#-----------------------------------------------------------
#---- 引入类库
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.Data.MapData import *
from MeteoInfoC.Data.MeteoData import *
from MeteoInfoC.Layout import *
from MeteoInfoC.Projections import *
#---- 设置路径变量
baseDir = "C:\\Program Files (x86)\\MeteoInfo\\"
mapDir = baseDir + "Map\\"
legendDir = baseDir + "pal\\"
dataDir = "D:\\Temp\\micaps\\"
#---- 打开图层
bou2Layer = MapDataManage.OpenLayer(mapDir + "bou2_4p.shp")
lb = bou2Layer.LegendScheme.LegendBreaks[0]
lb.Color = Color.Yellow
lb.OutlineColor = Color.Gray
lb.OutlineSize = 1
lb.DrawFill = False
bou1Layer = MapDataManage.OpenLayer(mapDir + "bou1_4l.shp")
lb = bou1Layer.LegendScheme.LegendBreaks[0]
lb.Color = Color.Red
res1Layer = MapDataManage.OpenLayer(mapDir + "res1_4m.shp")
lb = res1Layer.LegendScheme.LegendBreaks[0]
lb.Color = Color.Red
lb.OutlineColor = Color.Black
lb.Size = 5
res1Layer.LabelSet.FieldName = "NAME";
res1Layer.LabelSet.LabelFont = Font("楷体", 10);
res1Layer.LabelSet.YOffset = 15;
res1Layer.LabelSet.LabelColor = Color.Red
res1Layer.AddLabels();
#---- 创建MIApp类的对象
myApp = MIApp()
mapLayout = myApp.MapLayout
layoutMap = mapLayout.ActiveLayoutMap
mapFrame = layoutMap.MapFrame
mapView = mapFrame.MapView
#---- 设置Layout
mapLayout.SetPaperSize(540, 550)
#---- 添加地图图层
mapFrame.AddLayer(bou2Layer)
mapFrame.AddLayer(bou1Layer)
mapFrame.AddLayer(res1Layer)
#---- 设置ActiveLayoutMap（图层显示）
layoutMap.DrawGridLine = False
layoutMap.DrawNeatLine = False
layoutMap.DrawGridLabel = False
layoutMap.DrawGridTickLine = False
layoutMap.Left = 10
layoutMap.Top = 10
layoutMap.Width = 500
layoutMap.Height = 490
#---- 设置图形标题
aText = mapLayout.AddText("Temp",250,50,"黑体",12)
aText.Graphic.Legend.Color = Color.Red
#---- 设置时间
sTime = DateTime.Parse("2013-08-22 16:00")
#---- 设置图例
palfn = legendDir + "I-01.pal"
#---- 设置MeteoDataInfo
mid = MeteoDataInfo()
#---- 设置数据文件
inFile = sTime.ToString("yyMMddHH") + "00.000"
print inFile
#---- 打开MICAPS数据文件
mid.OpenMICAPSData(dataDir + inFile)
#---- 投影
mapView.ProjectLayers(mid.ProjInfo)
#---- 获取云图格点数据
gData = mid.GetGridData("var")
#---- 生成云图图层
cloudLayer = DrawMeteoData.CreateRasterLayer(gData, "Cloud", palfn)
cloudLayer.ProjInfo = mid.ProjInfo
#---- 添加图层
mapFrame.AddLayer(cloudLayer)
print "Display finished"
#---- 调整图层顺序（以避免压盖）
mapFrame.MoveLayerToBottom(cloudLayer)
#---- 缩放至云图范围
mapView.ZoomToExtentEx(cloudLayer.Extent)
#---- 设置标题名称
title = "可见光云图(" + sTime.ToString("yyyy-MM-dd HH:00") + ")"
aText.SetLabelText(title)
#---- 绘制图形
mapLayout.PaintGraphics()
#---- 输出图形为文件
outFile = "Cloud_" + sTime.ToString("yyyyMMddHH") + ".png"
print outFile
#mapLayout.ExportToPicture("E:\\Temp\\" + outFile)
#---- 显示程序窗体（只是为了看效果，自动运行时不需要）
#Application.Run(myApp)
myApp.Show()

# This Python file uses the following encoding: GB2312 
#--------------------------------------------------------         
# Author: Yaqiang Wang                                           
# Date: 2013-8-27                                                 
# Purpose: Read and plot MICAPS cloud data   
# Note: Sample                                                    
#-----------------------------------------------------------   
#---- 引入类库 
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.Data.MapData import * 
from MeteoInfoC.Data.MeteoData import * 
from MeteoInfoC.Layout import * 
from MeteoInfoC.Projections import * 
 
#---- 设置路径变量 
baseDir = "C:\\Program Files (x86)\\MeteoInfo\\" 
mapDir = baseDir + "Map\\" 
legendDir = baseDir + "pal\\" 
dataDir = "D:\\Temp\\micaps\\" 
 
#---- 打开图层 
bou2Layer = MapDataManage.OpenLayer(mapDir + "bou2_4p.shp") 
lb = bou2Layer.LegendScheme.LegendBreaks[0] 
lb.Color = Color.Yellow 
lb.OutlineColor = Color.Gray 
lb.OutlineSize = 1 
lb.DrawFill = False 
bou1Layer = MapDataManage.OpenLayer(mapDir + "bou1_4l.shp") 
lb = bou1Layer.LegendScheme.LegendBreaks[0] 
lb.Color = Color.Red 
res1Layer = MapDataManage.OpenLayer(mapDir + "res1_4m.shp") 
lb = res1Layer.LegendScheme.LegendBreaks[0] 
lb.Color = Color.Red 
lb.OutlineColor = Color.Black 
lb.Size = 5 
res1Layer.LabelSet.FieldName = "NAME"; 
res1Layer.LabelSet.LabelFont = Font("楷体", 10); 
res1Layer.LabelSet.YOffset = 15; 
res1Layer.LabelSet.LabelColor = Color.Red 
res1Layer.AddLabels(); 
 
#---- 创建MIApp类的对象 
myApp = MIApp() 
mapLayout = myApp.MapLayout 
layoutMap = mapLayout.ActiveLayoutMap 
mapFrame = layoutMap.MapFrame 
mapView = mapFrame.MapView 
 
#---- 设置Layout 
mapLayout.SetPaperSize(540, 550) 
 
#---- 添加地图图层 
mapFrame.AddLayer(bou2Layer) 
mapFrame.AddLayer(bou1Layer) 
mapFrame.AddLayer(res1Layer) 
 
#---- 设置ActiveLayoutMap（图层显示） 
layoutMap.DrawGridLine = False 
layoutMap.DrawNeatLine = False 
layoutMap.DrawGridLabel = False 
layoutMap.DrawGridTickLine = False 
layoutMap.Left = 10 
layoutMap.Top = 10 
layoutMap.Width = 500 
layoutMap.Height = 490 
 
#---- 设置图形标题 
aText = mapLayout.AddText("Temp",250,50,"黑体",12) 
aText.Graphic.Legend.Color = Color.Red 
 
#---- 设置时间 
sTime = DateTime.Parse("2013-08-22 16:00") 
 
#---- 设置图例 
palfn = legendDir + "I-01.pal" 
 
#---- 设置MeteoDataInfo 
mid = MeteoDataInfo() 
 
#---- 设置数据文件   
inFile = sTime.ToString("yyMMddHH") + "00.000" 
print inFile 
#---- 打开MICAPS数据文件 
mid.OpenMICAPSData(dataDir + inFile)   
#---- 投影 
mapView.ProjectLayers(mid.ProjInfo) 
#---- 获取云图格点数据 
gData = mid.GetGridData("var") 
#---- 生成云图图层 
cloudLayer = DrawMeteoData.CreateRasterLayer(gData, "Cloud", palfn) 
cloudLayer.ProjInfo = mid.ProjInfo 
#---- 添加图层 
mapFrame.AddLayer(cloudLayer) 
print "Display finished"   
#---- 调整图层顺序（以避免压盖） 
mapFrame.MoveLayerToBottom(cloudLayer)   
#---- 缩放至云图范围 
mapView.ZoomToExtentEx(cloudLayer.Extent) 
#---- 设置标题名称 
title = "可见光云图(" + sTime.ToString("yyyy-MM-dd HH:00") + ")" 
aText.SetLabelText(title) 
#---- 绘制图形 
mapLayout.PaintGraphics() 
#---- 输出图形为文件 
outFile = "Cloud_" + sTime.ToString("yyyyMMddHH") + ".png" 
print outFile 
#mapLayout.ExportToPicture("E:\\Temp\\" + outFile)   
 
#---- 显示程序窗体（只是为了看效果，自动运行时不需要） 
#Application.Run(myApp) 
myApp.Show()

运行结果：

才仔 · 发表于 2013-8-27 15:38:20

大神威武

La_respiration · 发表于 2013-8-27 17:18:56

非常感谢老师

Role · 发表于 2013-8-28 02:00:08

好东西感谢共享

liang · 发表于 2013-8-28 16:09:45

学习一下，收藏备用，多谢分享

寻路人 · 发表于 2013-9-9 12:06:47

大神威武

huangfucheng · 发表于 2014-12-24 15:16:41

老师请问下脚本怎么执行呢

墨墨卡托 · 发表于 2015-6-18 11:23:32

楼主可以分享个MeteoInfo读取、生成FY—2C的可见光、红外、水汽通道的代码吗？

MeteoInfo · 发表于 2015-6-18 11:52:36

墨墨卡托发表于 2015-6-18 11:23
楼主可以分享个MeteoInfo读取、生成FY—2C的可见光、红外、水汽通道的代码吗？

数据是什么？

墨墨卡托 · 发表于 2015-6-18 16:13:19

MeteoInfo 发表于 2015-6-18 11:52
数据是什么？

老师您好，数据虽然不大，但也不小，具体数据我随后发您的邮箱，这里提前写过老师了~

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