也说格点数据插值为站点数据

朝雾

MeteoInfo 发表于 2013-10-25 09:02

                               
登录/注册后可看大图

先自己琢磨吧，有具体问题了再说。

楼主，本人不了解 Python 语言的基本语法。这是GRADS的ctl
dset  rcp45_%y4.dat
title RegCM daily output variables
options little_endian
options template
options  365_day_calendar
undef -1.e34
xdef  27 linear   102  0.5
ydef  31 linear   27 0.5
zdef 1 levels 1008.23
tdef   34310 linear 01jan2006 1dy
vars  4
t2m     0 99 air temperature at 2m (K)
tpr     0 99 total precipitation (mm/day)
t2max   0 99 maximum 2m air temperature (K)
t2min   0 99 minimum 2m air temperature (K)
endvars
时间层次很多，软件脚本可以批量处理吗？
#---- 设置路径变量
BaseDir = "e:\\Program Files\\MeteoInfo\\"
MapDir = "e:\\Program Files\\MeteoInfo\\Map\\"
LegendDir = "e:\\Program Files\\MeteoInfo\\Legend\\"
DataDir = "e:\\Program Files\\MeteoInfo\\Sample\\GrADS\\"
OutDir = "E:\\Temp\\"
#---- 创建MIApp类的对象
myApp = MIApp()
#---- 打开格点数据文件
df = "J:\\CMIP5-data\\RegCM-data1\\data\\rcp45.ctl"
myApp.OpenGrADSData(df)
#---- 获取地面气压格点数据
gData = myApp.GetGridData("t2m")
#---- 设置输入、输出站点数据文件
inf = "J:\\CMIP5-data\\RegCM-data1\\data\\station.csv"
outf = "J:\\CMIP5-data\\RegCM-data1\\data\\stationtest.csv"
#---- 格点数据插值为站点数据，并输出为文件
gData.ToStation(inf, outf)
怎么修改？麻烦楼主指导下，谢谢。

MeteoInfo

朝雾 发表于 2013-10-25 11:12

                               
登录/注册后可看大图

楼主，本人不了解 Python 语言的基本语法。这是GRADS的ctl
dset  rcp45_%y4.dat
title RegCM daily out ...

参考此帖：站点降水数据绘制等值线图新脚本
http://bbs.06climate.com/forum.p ... 695&fromuid=106

里面有时间设置。

朝雾

MeteoInfo 发表于 2013-10-25 11:21

                               
登录/注册后可看大图

参考此帖：站点降水数据绘制等值线图新脚本
http://bbs.06climate.com/forum.php?mod=viewthread&tid=15 ...

好的，谢谢楼主

yanyan83518

最近一直在用meteoinfo，真的是好强大！而且很容易上手，赞

朝雾

MeteoInfo 发表于 2013-10-25 11:21

                               
登录/注册后可看大图

参考此帖：站点降水数据绘制等值线图新脚本
http://bbs.06climate.com/forum.php?mod=viewthread&tid=15 ...

sTime = DateTime.Parse("2006-01-01 00:00")
eTime = DateTime.Parse("2099-12-08 00:00")
while sTime <= eTime:
  inFile = sTime.ToString("yyMMddHH")
  print inFile
  gData = myApp.GetGridData("t2m")
  #---- 设置输入、输出站点数据文件
  inf = "J:\\CMIP5-data\\RegCM-data1\\data\\station.csv"
  outf = "J:\\CMIP5-data\\RegCM-data1\\data\\stationtest.csv"
  #---- 格点数据插值为站点数据，并输出为文件
  gData.ToStation(inf,outf)
  sTime = sTime.AddHours(24)
#Application.Run(myApp)

楼主循环运行后数据会被替代，怎么才能在csv文件中逐列写入，并显示时间？有可以参考的吗？十分感谢

MeteoInfo

朝雾 发表于 2013-10-28 15:07

                               
登录/注册后可看大图

sTime = DateTime.Parse("2006-01-01 00:00")
eTime = DateTime.Parse("2099-12-08 00:00")
while sT ...

outf文件名中加上时间就可以区分开了。

如果想将不同时次插值的结果写入一个文件中需要自己写Python代码读station.csv文件（将所有站点的站号及经纬度读出来），用GridData的ToStation(double x, double y)方法得到每个站点的插值结果（x, y参数分别对应站点经纬度），然后自己写代码去写插值结果文件。

朝雾

MeteoInfo 发表于 2013-10-28 16:09

                               
登录/注册后可看大图

outf文件名中加上时间就可以区分开了。

如果想将不同时次插值的结果写入一个文件中需要自己写Python代 ...

十分感谢楼主的悉心指导，受益匪浅，谢谢！

znk1981

学习中，好东东！！！！

朝雾

MeteoInfo 发表于 2013-10-28 16:09

                               
登录/注册后可看大图

outf文件名中加上时间就可以区分开了。

如果想将不同时次插值的结果写入一个文件中需要自己写Python代 ...

如果想将不同时次插值的结果写入一个文件中需要自己写Python代码读station.csv文件（将所有站点的站号及经纬度读出来），用GridData的ToStation(double x, double y)方法得到每个站点的插值结果（x, y参数分别对应站点经纬度），然后自己写代码去写插值结果文件。

楼主，具体这块有什么可以参考的帖子吗？

MeteoInfo

朝雾 发表于 2013-11-3 20:33

                               
登录/注册后可看大图

如果想将不同时次插值的结果写入一个文件中需要自己写Python代码读station.csv文件（将所有站点的站号及经 ...

帖子你可以在论坛里用搜索功能找。关键你要学习Python读写文件的功能。