爱气象，爱气象家园！ 
·设为首页 开启辅助访问 新帖子新回复热门贴精华帖

气象家园

 找回密码
 立即注册

QQ登录

只需一步，快速开始

新浪微博登陆

只需一步, 快速开始

搜索
气象家园»论坛 › 精品推荐 › MeteoInfo › MeteoInfo脚本示例：站点数据分析
12下一页
返回列表 发新帖
查看: 14195|回复: 14

MeteoInfo脚本示例：站点数据分析

[复制链接]
MeteoInfo

新浪微博达人勋
发表于 2012-10-27 14:08:39 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 MeteoInfo 于 2012-10-27 16:43 编辑

和格点数据类似，MeteoInfo也可以对站点数据进行各种数学运算和分析。之前讲过利用MICAPS 1数据（站点数据）绘制累计降水量的脚本，这里用新版的脚本程序再说明一下。程序演示了读取两个时次的MICAPS 1数据中的6小时降水量站点数据，并进行相加获得12小时累计降水量的站点数据。为了绘制降水量等值线图需要将站点数据插值为格点数据，在脚本中设置插值方法和参数即可。 需要将软件更新为最新文件（见置顶帖子）。

Image00937.png

脚本代码如下：
  1. # This Python file uses the following encoding: GB2312
  2. #---- 引入类库
  3. import clr
  4. clr.AddReferenceByPartialName("System")
  5. clr.AddReferenceByPartialName("System.Windows.Forms")
  6. clr.AddReferenceByPartialName("System.Drawing")
  7. from System import *
  8. from System.Windows.Forms import *
  9. from System.Drawing import *
  10. clr.AddReference("MeteoInfoC.dll")
  11. from MeteoInfoC import *
  12. from MeteoInfoC.Layout import *
  13. from MeteoInfoC.Data import *
  14. from MeteoInfoC.Data.MeteoData import *
  15. from MeteoInfoC.Legend import *
  16. from MeteoInfoC.Global import *
  17. from MeteoInfoC.Shape import *

  19. #---- 设置路径变量
  20. BaseDir = "C:\\Program Files (x86)\\MeteoInfo\"
  21. MapDir = BaseDir + "Map\"
  22. LegendDir = BaseDir + "Legend\"
  23. DataDir = BaseDir + "Sample\\MICAPS\"

  25. #---- 设置屏蔽图层（只绘制中国境内图形）
  26. mipy.MapDocument.ActiveMapFrame.MapView.MaskOut.SetMaskLayer = True
  27. mipy.MapDocument.ActiveMapFrame.MapView.MaskOut.MaskLayer = "china.shp"

  29. #---- 设置绘图类型为shaded
  30. mipy.MeteoDataset.SetDrawType("shaded")

  32. #---- 站点数据插值为格点数据的设置
  33. gDataSet = GridDataSetting()
  34. gDataSet.DataExtent = Extent(60, 140, -20, 60)
  35. gDataSet.XNum = 160
  36. gDataSet.YNum = 160
  37. interpSet = InterpolationSetting()
  38. interpSet.GridDataSet = gDataSet
  39. interpSet.InterpolationMethod = InterpolationMethods.IDW_Radius
  40. interpSet.Radius = 2
  41. interpSet.MinPointNum = 1
  42. mipy.MeteoDataset.InterpolationSet = interpSet;
  43. mipy.MeteoDataset.UseDefaultInterpSet = True

  45. #---- 设置图例文件
  46. mipy.MeteoDataset.LegendScheme = LegendScheme(ShapeTypes.Polygon)
  47. mipy.MeteoDataset.LegendScheme.ImportFromXMLFile(LegendDir + "rain1.lgs")
  48. mipy.MeteoDataset.UseDefaultLegendScheme = True

  50. #---- 设置起始结束时间
  51. sTime = DateTime.Parse("2010-10-14 14:00")
  52. eTime = DateTime.Parse("2010-10-14 20:00")
  53. aTime = sTime.AddHours(-6)

  55. dataInfo = MeteoDataInfo()
  56. stData = StationData()
  57. i = 1
  58. #---- 循环
  59. while sTime <= eTime:  
  60.   inFile = sTime.ToString("yyMMddHH") + ".000"
  61.   print inFile
  62.   #---- 打开MICAPS数据文件
  63.   dataInfo.OpenMICAPSData(DataDir + inFile)  
  64.   #---- 获取累计降水站点数据
  65.   if i == 1:
  66.           stData = dataInfo.GetStationData("Precipitation6h")
  67.   else:
  68.           stData = stData + dataInfo.GetStationData("Precipitation6h")
  69.   
  70.   #---- 时间加6小时
  71.   sTime = sTime.AddHours(6)
  72.   i = i + 1

  74. #---- 绘制累计降水图层
  75. aLayer = mipy.MeteoDataset.Display(stData)
  76. aLayer.LayerName = "累计降水"
  77. mipy.MapDocument.ActiveMapFrame.UpdateLayerNode(aLayer)
  78. mipy.MapDocument.ActiveMapFrame.MoveLayer(aLayer, 0)

  80. print "Finished!"
复制代码


密码修改失败请联系微信：mofangbao
回复

使用道具 举报

freekiller

新浪微博达人勋
发表于 2012-10-27 18:07:41 | 显示全部楼层
顶               
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

jeydragon

新浪微博达人勋
发表于 2012-10-29 10:49:16 | 显示全部楼层
接着顶起来
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

♂雨已~○

新浪微博达人勋
发表于 2012-10-29 18:46:03 | 显示全部楼层
又学习了，好东西啊
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

weiguangjian

新浪微博达人勋
发表于 2012-11-21 19:46:43 | 显示全部楼层
本帖最后由 weiguangjian 于 2012-11-21 20:16 编辑
  1. # This Python file uses the following encoding: GB2312

  3. #---- 引入类库

  5. import clr

  7. clr.AddReferenceByPartialName("System")

  9. clr.AddReferenceByPartialName("System.Windows.Forms")

  11. clr.AddReferenceByPartialName("System.Drawing")

  13. from System import *

  15. from System.Windows.Forms import *

  17. from System.Drawing import *

  19. clr.AddReference("MeteoInfoC.dll")

  21. from MeteoInfoC import *

  23. from MeteoInfoC.Layout import *

  25. from MeteoInfoC.Data import *

  27. from MeteoInfoC.Data.MeteoData import *

  29. from MeteoInfoC.Legend import *

  31. from MeteoInfoC.Global import *

  33. from MeteoInfoC.Shape import *


  36. #---- 设置路径变量

  38. BaseDir = "C:\\Program Files (x86)\\MeteoInfo\"

  40. MapDir = BaseDir + "Map\"

  42. LegendDir = BaseDir + "Legend\"

  44. DataDir = BaseDir + "Sample\\MICAPS\"


  47. #---- 设置屏蔽图层（只绘制中国境内图形）

  49. mipy.MapDocument.ActiveMapFrame.MapView.MaskOut.SetMaskLayer = True

  51. mipy.MapDocument.ActiveMapFrame.MapView.MaskOut.MaskLayer = "china.shp"


  54. #---- 设置绘图类型为shaded

  56. mipy.MeteoDataset.SetDrawType("shaded")

  58. #---- 站点数据插值为格点数据的设置

  60. gDataSet = GridDataSetting()

  62. gDataSet.DataExtent = Extent(110, 130, 30, 45)

  64. gDataSet.XNum = 130

  66. gDataSet.YNum = 130

  68. interpSet = InterpolationSetting()

  70. interpSet.GridDataSet = gDataSet

  72. interpSet.InterpolationMethod = InterpolationMethods.IDW_Radius

  74. interpSet.Radius = 1

  76. interpSet.MinPointNum = 1

  78. mipy.MeteoDataset.InterpolationSet = interpSet;

  80. mipy.MeteoDataset.UseDefaultInterpSet = True


  83. #---- 设置图例文件

  85. mipy.MeteoDataset.LegendScheme = LegendScheme(ShapeTypes.Polygon)

  87. mipy.MeteoDataset.LegendScheme.ImportFromXMLFile(LegendDir + "rain1.lgs")

  89. mipy.MeteoDataset.UseDefaultLegendScheme = True


  92. #---- 设置起始结束时间

  94. sTime = DateTime.Parse("2012-09-21 14:00")

  96. eTime = DateTime.Parse("2012-09-22 08:00")

  98. aTime = sTime.AddHours(-6)


  101. dataInfo = MeteoDataInfo()

  103. stData = StationData()

  105. i = 1

  107. #---- 循环

  109. while sTime <= eTime:  

  111.   inFile = sTime.ToString("yyMMddHH") + ".sfc"

  113.   print inFile

  115.   #---- 打开MICAPS数据文件

  117.   dataInfo.OpenMICAPSData(DataDir + inFile)  

  119.   #---- 获取累计降水站点数据

  121.   if i == 1:

  123.           stData = dataInfo.GetStationData("Precipitation6h")

  125.   else:

  127.           stData = stData + dataInfo.GetStationData("Precipitation6h")

  129.   

  131.   #---- 时间加6小时

  133.   sTime = sTime.AddHours(6)

  135.   i = i + 1


  138. #---- 绘制累计降水图层

  140. aLayer = mipy.MeteoDataset.Display(stData)

  142. aLayer.LayerName = "累计降水"

  144. mipy.MapDocument.ActiveMapFrame.UpdateLayerNode(aLayer)

  146. mipy.MapDocument.ActiveMapFrame.MoveLayer(aLayer, 0)


  149. print "Finished!"
复制代码
你好！我想用此脚本绘制24小时的累计雨量（2012092108-2012092208），修改执行出现如下错误（这只是出现的第一个错误，后面应该还有错）：
Traceback (most recent call last):
  File "<string>", line 60, in <module>
NameError: name 'GridDataSetting' is not defined
初次接触，除了简单的参数设置外，不知道脚本里面还应做哪些设置。以下是简单设置后的代码：还请赐教，不胜感激！所用数据为micaps一类数据
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

MeteoInfo

新浪微博达人勋
 楼主| 发表于 2012-11-21 20:28:12 | 显示全部楼层
weiguangjian 发表于 2012-11-21 19:46

                               
登录/注册后可看大图

你好！我想用此脚本绘制24小时的累计雨量（2012092108-2012092208），修改执行出现如下错误（这只是出现的第 ...

需要将软件更新为最新文件（见置顶帖子）。
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

weiguangjian

新浪微博达人勋
发表于 2012-11-21 21:03:35 | 显示全部楼层
MeteoInfo 发表于 2012-11-21 20:28

                               
登录/注册后可看大图

需要将软件更新为最新文件（见置顶帖子）。

仍然报错：
Traceback (most recent call last):
  File "<string>", line 87, in <module>
IOError: [Errno 2] Could not find a part of the path 'C:\Program Files (x86)\MeteoInfo\Legend\rain1.lgs'.
我看了，在这个路径没有错啊
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

MeteoInfo

新浪微博达人勋
 楼主| 发表于 2012-11-21 21:28:12 | 显示全部楼层
weiguangjian 发表于 2012-11-21 21:03

                               
登录/注册后可看大图

仍然报错：
Traceback (most recent call last):
  File "", line 87, in

你再仔细查查路径，如果报错肯定是路径有问题。
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

weiguangjian

新浪微博达人勋
发表于 2012-11-21 22:16:02 | 显示全部楼层
MeteoInfo 发表于 2012-11-21 21:28

                               
登录/注册后可看大图

你再仔细查查路径，如果报错肯定是路径有问题。

画出图来了，把BaseDir = "C:\\Program Files (x86)\\MeteoInfo\\"中的（x86）去掉就好了。但是在中国区域外的值怎么没去掉，海上不应该有值的啊？并且没有中国底图呢？！
1.gif

密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

MeteoInfo

新浪微博达人勋
 楼主| 发表于 2012-11-21 23:09:43 | 显示全部楼层
weiguangjian 发表于 2012-11-21 22:16

                               
登录/注册后可看大图

画出图来了，把BaseDir = "C:\\Program Files (x86)\\MeteoInfo\\"中的（x86）去掉就好了。但是在中国区域 ...

你看看MeteoInfo软件的帮助文档以及之前的一些帖子，把软件的基本操作弄清楚。
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

下一页 »
12下一页
返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

小黑屋|帮助中心|清除痕迹

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表，不代表气象家园立场，禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表