MeteoInfo脚本示例：台风

MeteoInfo · 发表于 2012-9-26 13:03:13

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

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

x

本帖最后由 MeteoInfo 于 2012-11-18 16:39 编辑

从台风数据文件中读取每个点的信息，创建台风路径点图层和线图层（包含了相应的属性信息，比如时间、风速等）。数据文件：

198214track.dat (2.28 KB, 下载次数: 210)

2012-9-26 12:59 上传

点击文件名下载附件

代码如下：

#--------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2012-9-26
# Purpose: Read and plot typhoon data
# Note: Sample
#-----------------------------------------------------------
#---- Import clr and classes
import clr
clr.AddReferenceByPartialName("System")
clr.AddReferenceByPartialName("System.Drawing")
from System import *
from System.Drawing import *
from System.Collections.Generic import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data.MapData import *
from MeteoInfoC.Geoprocess import *
from MeteoInfoC.Layer import *
from MeteoInfoC.Shape import *
from MeteoInfoC.Legend import *
import os.path
#---- Create a point and line layers to plot typhoon pathway
print 'Create typhoon point layer...'
pLayer = VectorLayer(ShapeTypes.Point)
pLayer.LayerName = "Typhoon_Point"
pLayer.Visible = True
pLayer.EditAddField("Time", String)
pLayer.EditAddField("Wind", Int32)
pLayer.EditAddField("Status", String)
print 'Create typhoon polyline layer...'
lLayer = VectorLayer(ShapeTypes.Polyline)
lLayer.LayerName = "Typhoon_Line"
lLayer.Visible = True
lLayer.EditAddField("ID", String)
#---- Read typhoon data file
print 'Read typhoon data file...'
fn = "D:\\Temp\\ascii\\198214track.dat"
tf = open(fn)
tf.readline()
aline = tf.readline()
id = aline.split()[1]
tf.readline()
aPolyline = PolylineShape()
plist = List[PointD]()
for aline in tf:
print aline
datalist = aline.split()
lat = float(datalist[1])
lon = float(datalist[2])
t = datalist[3]
wind = int(datalist[4])
stat = datalist[6]
if len(datalist) == 8:
stat = stat + ' ' + datalist[7]
aPS = PointShape()
aPoint = PointD(lon, lat)
aPS.Point = aPoint
plist.Add(aPoint)
shapeNum = pLayer.ShapeNum
if pLayer.EditInsertShape(aPS, shapeNum):
pLayer.EditCellValue("Time", shapeNum, t)
pLayer.EditCellValue("Wind", shapeNum, wind)
pLayer.EditCellValue("Status", shapeNum, stat)
pLayer.UpdateLegendScheme(LegendType.GraduatedColor, "Wind")
for legend in pLayer.LegendScheme.LegendBreaks:
legend.Size = 10
aPolyline.Points = plist
shapeNum = lLayer.ShapeNum
if lLayer.EditInsertShape(aPolyline, shapeNum):
lLayer.EditCellValue("ID", shapeNum, id)
lLayer.UpdateLegendScheme(LegendType.SingleSymbol, "ID")
lLayer.LegendScheme.LegendBreaks[0].Size = 2
mipy.MapDocument.ActiveMapFrame.AddLayer(lLayer)
mipy.MapDocument.ActiveMapFrame.AddLayer(pLayer)
print 'Finished!'

复制代码

表现力还可以再丰富一些，将设置点图层LegendSchem的代码修改如下：

for legend in pLayer.LegendScheme.LegendBreaks:
legend.MarkerType = MarkerType.Character
legend.FontName = "Weather"
legend.CharIndex = 170
legend.Size = 12

复制代码

如果想台风路径线条不同段用不同颜色表示，可以每一段生成一根线条，然后分段渲染。

代码如下：

#--------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2012-11-18
# Purpose: Read and plot typhoon data - polyline with different color
# Note: Sample
#-----------------------------------------------------------
#---- Import clr and classes
import clr
clr.AddReferenceByPartialName("System")
clr.AddReferenceByPartialName("System.Drawing")
from System import *
from System.Drawing import *
from System.Collections.Generic import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data.MapData import *
from MeteoInfoC.Geoprocess import *
from MeteoInfoC.Layer import *
from MeteoInfoC.Shape import *
from MeteoInfoC.Legend import *
from MeteoInfoC.Drawing import *
import os.path
#---- Create a point and line layers to plot typhoon pathway
print 'Create typhoon point layer...'
pLayer = VectorLayer(ShapeTypes.Point)
pLayer.LayerName = "Typhoon_Point"
pLayer.Visible = True
pLayer.EditAddField("Time", String)
pLayer.EditAddField("Wind", Int32)
pLayer.EditAddField("Status", String)
print 'Create typhoon polyline layer...'
lLayer = VectorLayer(ShapeTypes.Polyline)
lLayer.LayerName = "Typhoon_Line"
lLayer.Visible = True
lLayer.EditAddField("ID", String)
lLayer.EditAddField("Time", String)
lLayer.EditAddField("Wind", Int32)
lLayer.EditAddField("Status", String)
#---- Read typhoon data file
print 'Read typhoon data file...'
fn = "D:\\Temp\\ascii\\198214track.dat"
tf = open(fn)
tf.readline()
aline = tf.readline()
id = aline.split()[1]
tf.readline()
plist = List[PointD]()
i = 0
for aline in tf:
print aline
datalist = aline.split()
lat = float(datalist[1])
lon = float(datalist[2])
t = datalist[3]
wind = int(datalist[4])
stat = datalist[6]
if len(datalist) == 8:
stat = stat + ' ' + datalist[7]
aPS = PointShape()
aPoint = PointD(lon, lat)
aPS.Point = aPoint
plist.Add(aPoint)
shapeNum = pLayer.ShapeNum
if pLayer.EditInsertShape(aPS, shapeNum):
pLayer.EditCellValue("Time", shapeNum, t)
pLayer.EditCellValue("Wind", shapeNum, wind)
pLayer.EditCellValue("Status", shapeNum, stat)
if plist.Count > 1:
aPolyline = PolylineShape()
aPolyline.Points = plist
shapeNum = lLayer.ShapeNum
if lLayer.EditInsertShape(aPolyline, shapeNum):
lLayer.EditCellValue("ID", shapeNum, id)
lLayer.EditCellValue("Time", shapeNum, t)
lLayer.EditCellValue("Wind", shapeNum, wind)
lLayer.EditCellValue("Status", shapeNum, stat)
aPoint = PointD(plist[1].X, plist[1].Y)
plist = List[PointD]()
plist.Add(aPoint)
i += 1
pLayer.UpdateLegendScheme(LegendType.GraduatedColor, "Wind")
for legend in pLayer.LegendScheme.LegendBreaks:
legend.MarkerType = MarkerType.Character
legend.FontName = "Weather"
legend.CharIndex = 170
legend.Size = 12
lLayer.UpdateLegendScheme(LegendType.GraduatedColor, "Wind")
for legend in lLayer.LegendScheme.LegendBreaks:
legend.Size = 2
mipy.MapDocument.ActiveMapFrame.AddLayer(lLayer)
mipy.MapDocument.ActiveMapFrame.AddLayer(pLayer)
print 'Finished!'

复制代码

良辰 · 发表于 2012-9-26 16:28:22

非常感谢王老师热情的帮助，我还是第一次用脚本绘图，挺有意思的，但是有点不懂这个语言，好像是Python，没研究过，想问一下王老师，那个设置点图层LegendSchem的代码加在什么位置，万分感激！

MeteoInfo · 发表于 2012-9-26 16:31:38

良辰发表于 2012-9-26 16:28

登录/注册后可看大图

非常感谢王老师热情的帮助，我还是第一次用脚本绘图，挺有意思的，但是有点不懂这个语言，好像是Python，没 ...

是Python的一个分支：IronPython，可以和.Net很好地衔接。第二段代码可以替换第一段代码的69和70行。

良辰 · 发表于 2012-9-26 16:38:02

MeteoInfo 发表于 2012-9-26 16:31

登录/注册后可看大图

是Python的一个分支：IronPython，可以和.Net很好地衔接。第二段代码可以替换第一段代码的69和70行。

谢谢王老师，我如果批处理能做吗？就是画几十条台风轨迹，fn= "F:\\taifeng\\his_typhoon_track\\194806track.dat"
tf = open(fn)，怎么修改，可以把fn定义成数组吗，进行循环打开。希望你先帮帮忙，以后要好好学习Python这个语言了，很不错，主要是metoinfo的脚本处理！

良辰 · 发表于 2012-9-26 16:42:22

F:\未命名.jpg出错了

良辰 · 发表于 2012-9-26 16:42:55

MeteoInfo · 发表于 2012-9-26 16:48:36

良辰发表于 2012-9-26 16:38

登录/注册后可看大图

谢谢王老师，我如果批处理能做吗？就是画几十条台风轨迹，fn= "F:\\taifeng\\his_typhoon_track\\194806t ...

脚本程序很适合来做批处理。你可以把所有的台风文件放在一个目录中，在脚本中可以列出目录中所有的文件，获得一个文件列表，通过循环列表中每个文件就可以了。

dir = 'F:\\taifeng\\his_typhoon_track'
list = os.listdir(dir) #列出目录下的所有文件和目录
for line in list:
fn = os.path.join(dir,line)
.... #例子里的代码，注意代码缩进。

MeteoInfo · 发表于 2012-9-26 16:53:10

良辰发表于 2012-9-26 16:42

登录/注册后可看大图

是代码缩进的问题，你具体在网上查查Python语言代码缩进。

良辰 · 发表于 2012-9-26 16:56:01

MeteoInfo 发表于 2012-9-26 16:53

登录/注册后可看大图

是代码缩进的问题，你具体在网上查查Python语言代码缩进。

王老师，万分感谢，谢谢谢谢

良辰 · 发表于 2012-9-26 17:01:19

良辰发表于 2012-9-26 16:56

登录/注册后可看大图

王老师，万分感谢，谢谢谢谢

它说legend.MarkerType = MarkerType.Character，没有定义，这又是为嘛？

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