爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: MeteoInfo

MeteoInfoLab脚本汇总贴

  [复制链接]

新浪微博达人勋

 楼主| 发表于 2017-7-15 15:11:06 | 显示全部楼层
daibingbing 发表于 2017-7-15 09:59
老师,我想问一下有GRIB转ARL数据的的脚本吗?这个汇总贴没找到啊

参考这里:http://www.meteothinker.com/exam ... ctory/grib2arl.html
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-7-16 20:52:36 | 显示全部楼层
MeteoInfo 发表于 2017-7-15 15:11
参考这里:http://www.meteothinker.com/examples/meteoinfolab/trajectory/grib2arl.html

谢谢您
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-7-31 17:48:37 | 显示全部楼层
太给力了,谢谢分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-8-2 00:24:27 | 显示全部楼层
meteoinfo学习贴
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-8-11 12:07:54 | 显示全部楼层
{:5_213:}赞,好多python命令能用
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-8-14 13:39:50 | 显示全部楼层
本帖最后由 ♂雨已~○ 于 2017-8-14 14:18 编辑

求教王老师,如何从脚本中获取HDF文件的时间,并显示在图片上?代码咋个写呢,才开始接触MeteoInfoLab还不怎么熟悉,下面是我绘制HDF雷达数据代码:

folder = 'E:/radar/'
fns = '201708121755.dat.hdf'
fn = folder + fns
f = addfile(fn)
vname = 'CompositeReflectivity'
t = f[vname]
DZ= t[0,:,:]
#Plot
axesm()
mlayer1 = shaperead('F:/MeteoInfo_Java_1.4R4_Files/MeteoInfo/map/bou2_4p.shp')
mlayer2=shaperead('F:/MeteoInfo_Java_1.4R4_Files/MeteoInfo/map/airway/JINJIN.shp')
city_layer=shaperead('F:/MeteoInfo_Java_1.4R4_Files/MeteoInfo/map/res1_4m.shp')
geoshow(mlayer1)
geoshow(mlayer2,edgecolor='red')
geoshow(city_layer,facecolor='r',size=4,labelfield='NAME',fontname=u'楷体',fontsize=16,yoffset=15)

levs=[10,15,20,25,30,35,40,45,50,55,60,65,70]
cols=[(255,255,255),(0,162,232),(86,225,250),(3,207,14),(26,152,7),(255,242,0),(217,172,113),(255,147,74),(255,0,0),(204,0,0),(155,0,0),(236,21,236),(130,11,130),(184,108,208)]
layer=imshowm(DZ,levs,colors=cols)
colorbar(layer,shrink=1,aspect=60,orientation='horizontal')
title(vname)
axism([96,113,20,35])

print 'Finish...'

效果图:

1.PNG


我通过上面的代码是绘制出来的反射率是的第一个时次的数据,咋个从数据中获取时间,并在图中加上呢?比如里面的时间有20个层次,而我想画第二个时间层次的时候把时间自动加到图上,这个脚本怎么写呢?这个MeteoInfoLab怎么查询语法使用呢,有很多的使用方法不怎么熟悉,有没有一个使用查询呢?希望王老师解答呢···



密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-8-14 16:25:42 | 显示全部楼层
本帖最后由 MeteoInfo 于 2017-8-14 16:30 编辑
♂雨已~○ 发表于 2017-8-14 13:39
求教王老师,如何从脚本中获取HDF文件的时间,并显示在图片上?代码咋个写呢,才开始接触MeteoInfoLab还不 ...

获取第一个时次的时间
f.gettime(0)
可以参考这里:http://bbs.06climate.com/forum.php?mod=viewthread&tid=43434
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-8-14 16:54:09 | 显示全部楼层
MeteoInfo 发表于 2017-8-14 16:25
获取第一个时次的时间
f.gettime(0)
可以参考这里:http://bbs.06climate.com/forum.php?mod=viewthrea ...

谢谢王老师,我自己摸索出来了,不过有点复杂,也算实现了

  timeStamp=t.dimvalue(0)
  timeStamp *= 1000.0  
  #!print timeStamp  
  timearr = time.localtime(timeStamp)  
  otherStyleTime = time.strftime("%Y-%m-%d %H:%M:%S", timearr)  
  #print otherStyleTime

  title(vname+'\n'+str(otherStyleTime))

我添加了这些代码实现了,虽然复杂了点,好像也可以了,不过王老师这直接的方法效率更高啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-8-14 17:24:41 | 显示全部楼层
MeteoInfo 发表于 2017-8-14 16:25
获取第一个时次的时间
f.gettime(0)
可以参考这里:http://bbs.06climate.com/forum.php?mod=viewthrea ...

请问王老师,这个脚本能加载web地图吗?比如说百度,Google或者高德地图的离线地图文件?如果能,那么该怎么加载呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-8-14 23:00:32 | 显示全部楼层
♂雨已~○ 发表于 2017-8-14 17:24
请问王老师,这个脚本能加载web地图吗?比如说百度,Google或者高德地图的离线地图文件?如果能,那么该 ...

本贴1楼就有例子,你仔细看看。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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