MeteoInfo脚本示例：台风

随缘 · 发表于 2012-9-26 19:34:07

王老师，如果能和wrf模式运算结果wrfout直接衔接的话真是太好了，比如像ncl那样，有官网的脚本示例

MeteoInfo · 发表于 2012-9-26 22:29:10

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

登录/注册后可看大图

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

MarkerType是在MeteoInfoC.Drawing命名空间里，在前面加上：
from MeteoInfoC.Drawing import *

MeteoInfo · 发表于 2012-9-26 22:39:02

随缘发表于 2012-9-26 19:34

登录/注册后可看大图

王老师，如果能和wrf模式运算结果wrfout直接衔接的话真是太好了，比如像ncl那样，有官网的脚本示例

能具体说说“直接衔接”指的什么吗？

随缘 · 发表于 2012-9-27 08:27:50

MeteoInfo 发表于 2012-9-26 22:39

登录/注册后可看大图

能具体说说“直接衔接”指的什么吗？

哦，我的意思也就是像ncl那样开设一个官网的脚本例子，里面有针对wrf模式的部分，大多数人还对meteoinfo的脚本不是很熟悉，通过一些常用的降水啊，台风路径啊这样脚本示例，可以让一些水平较低的同学很快运用meteoinfo，对其脚本多一些了解

MeteoInfo · 发表于 2012-9-27 08:58:09

随缘发表于 2012-9-27 08:27

登录/注册后可看大图

哦，我的意思也就是像ncl那样开设一个官网的脚本例子，里面有针对wrf模式的部分，大多数人还对meteoinfo的 ...

很好的建议，一直也想在www.meteothinker.com上做类似的例子。

良辰 · 发表于 2012-9-27 09:43:59

本帖最后由良辰于 2012-9-27 09:51 编辑

王老师，您看我画出来的为什么和您画的不一样那？标志是字母a。不是台风的标志！呜呜呜
还有批处理的没有实现，代码我给您发出来，麻烦你看一下！
批处理代码如下：

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.MarkerType = MarkerType.Character
      legend.FontName = "weather_symbol"
      legend.CharIndex = 170
      legend.Size = 12

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!'

良辰 · 发表于 2012-9-27 09:52:28

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.MarkerType = MarkerType.Character
legend.FontName = "weather_symbol"
legend.CharIndex = 170
legend.Size = 12
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!'

复制代码

MeteoInfo · 发表于 2012-9-27 10:07:09

良辰发表于 2012-9-27 09:43

登录/注册后可看大图

王老师，您看我画出来的为什么和您画的不一样那？标志是字母a。不是台风的标志！呜呜呜
还有批处理的没有实 ...

你的计算机里没有Weather这个字体吗？（可以在word里用插入符合功能看看字体）

请贴出全部代码及错误信息。

良辰 · 发表于 2012-9-27 10:13:56

#--------------------------------------------------------
# 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 *
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)
#---- Read typhoon data file
print 'Read typhoon data file...'
dir = 'F:\\taifeng\\his_typhoon_track'
list = os.listdir(dir) #列出目录下的所有文件和目录
for line in list:
fn = os.path.join(dir,line)
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.MarkerType = MarkerType.Character
legend.FontName = "weather"
legend.CharIndex = 170
legend.Size = 12
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!'

复制代码

错误提示：
Create typhoon point layer...
Create typhoon polyline layer...
Read typhoon data file...
Traceback (most recent call last):
File "<string>", line 48, in <module>
发现的我计算机里真没有weather这个字体！

MeteoInfo · 发表于 2012-9-27 10:30:07

良辰发表于 2012-9-27 10:13

登录/注册后可看大图

错误提示：
Create typhoon point layer...
Create typhoon polyline layer...

安装MeteoInfo应该会注册上这个字体的。你可以在MeteoInfo的安装目录下找到WeatherSymbol.ttf这个字体文件，把它复制到Windows的字体目录中即可（通常是C:\Windows\Font），如果还有问题请在网上搜索如何在Windows中安装字体。

写MeteoInfo脚本需要了解Python的一些基本语法，这个论坛里也有Python板块，可以去看看，补充一些基本知识，循环体里的代码一定要正确缩进。

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