按照给定的空间范围提取站点数据

MeteoInfo

这个script程序展示了如何通过给定经纬度范围来提取站点数据，并将站点数据保存为csv文件（该文件可以在MeteoInfo中使用）。

# This Python file uses the following encoding: GB2312
#---- 引入类库
import clr
clr.AddReferenceByPartialName("System")
from System import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data import *

#---- 设置路径变量
BaseDir = "C:\\Program Files (x86)\\MeteoInfo\\"
DataDir = BaseDir + "Sample\\MICAPS\\"
OutDir = "E:\\temp\\"
#---- 创建MIApp类的对象
myApp = MIApp()

#---- 设置起始结束时间
sTime = DateTime.Parse("2010-10-14 14:00")
eTime = DateTime.Parse("2010-10-14 20:00")

#---- 循环
while sTime <= eTime:
  myApp.RemoveDataLayers()
  inFile = sTime.ToString("yyMMddHH") + ".000"
  print inFile
  #---- 打开MICAPS数据文件
  myApp.OpenMICAPSData(DataDir + inFile)
  #---- 降水
  #---- 获取站点数据
  stData = myApp.GetStationData("Precipitation6h")  
  #---- 根据经纬度范围提取数据（最小经度，最大经度，最小纬度，最大纬度）
  stData = stData.Extract(60,140,-20,60)
  #---- 保存站点数据
  stData.SaveAsCSVFile(OutDir + "Prec_" + sTime.ToString("yyyyMMddHH") + ".csv", "Precipitation")
  #---- 温度
  #---- 获取站点数据
  stData = myApp.GetStationData("Temperature")  
  #---- 根据经纬度范围提取数据（最小经度，最大经度，最小纬度，最大纬度）
  stData = stData.Extract(60,140,-20,60)
  #---- 保存站点数据
  stData.SaveAsCSVFile(OutDir + "Temp_" + sTime.ToString("yyyyMMddHH") + ".csv", "Temperature")

  #---- 时间加6小时
  sTime = sTime.AddHours(6)

mofangbao

谢谢楼主提供哈，这个代码给定的经纬度是不是得是矩形的？有时候需要提取某个省或者区域的站点，如果可以用底图文件作为掩膜，来提取底图范围内的站点那就更强大了，不知道能不能实现啊

MeteoInfo

mofangbao 发表于 2011-7-20 11:42

                               
登录/注册后可看大图

谢谢楼主提供哈，这个代码给定的经纬度是不是得是矩形的？有时候需要提取某个省或者区域的站点，如果可以用 ...

这个也是可以的，用StationData的Maskout方法来做。

MeteoInfo

用多边形图层提取区域内的站点数据。对于格点数据也有类似的方法，区域外的格点值被赋为无效值，和GrADS的maskout函数类似，不过不需要格点的maskout数据，只需要一个多边形图层就可以了。这个例子中用中国多边形图层提取中国区域内的站点数据。

# This Python file uses the following encoding: GB2312
#---- 引入类库
import clr
clr.AddReferenceByPartialName("System")
clr.AddReferenceByPartialName("System.Windows.Forms")
from System import *
from System.Windows.Forms import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data import *

#---- 设置路径变量
BaseDir = "C:\\Program Files (x86)\\MeteoInfo\\"
MapDir = BaseDir + "Map\\"
DataDir = BaseDir + "Sample\\MICAPS\\"
OutDir = "E:\\temp\\"
#---- 创建MIApp类的对象
myApp = MIApp()
myApp.OpenLayer(MapDir + "china.shp")
myApp.Zoom(70,140,10,60)
myApp.SetDrawType("station_point")

#---- 设置起始结束时间
sTime = DateTime.Parse("2010-10-14 20:00")
eTime = DateTime.Parse("2010-10-14 20:00")

#---- 循环
while sTime <= eTime:
  myApp.RemoveDataLayers()
  inFile = sTime.ToString("yyMMddHH") + ".000"
  print inFile
  #---- 打开MICAPS数据文件
  myApp.OpenMICAPSData(DataDir + inFile)
  #---- 降水
  #---- 获取站点数据
  stData = myApp.GetStationData("Precipitation6h")  
  #---- 根据多边形图层提取图层区域内的数据（MapView对象，图层名）
  stData = stData.Maskout(myApp.MapLayout.MapView, "china.shp")
  #---- 保存站点数据
  stData.SaveAsCSVFile(OutDir + "China_Prec_" + sTime.ToString("yyyyMMddHH") + ".csv", "Precipitation")  
  myApp.Display(stData)
  
  #---- 时间加6小时
  sTime = sTime.AddHours(6)

#---- 绘制图形
myApp.MapLayout.PaintGraphics()  
#---- 显示程序窗体
Application.Run(myApp)

mofangbao

MeteoInfo 发表于 2011-7-20 13:36

                               
登录/注册后可看大图

用多边形图层提取区域内的站点数据。对于格点数据也有类似的方法，区域外的格点值被赋为无效值，和GrADS的m ...

太佩服你了!!

ivanfish

学习学习，太需要这要的经验贴了

ivantq

佩服，获取数据的好方法

苏绮

学习了，佩服

kc121

脚本语言放在软件的哪个位置呢？

MeteoInfo

kc121 发表于 2012-2-26 22:47

                               
登录/注册后可看大图

脚本语言放在软件的哪个位置呢？

脚本语言的使用可以看我博客里的具体内容：http://hi.baidu.com/meteoinfo/bl ... fo%20Script/index/2