MeteoInfoLab脚本汇总贴

MeteoInfo

daibingbing 发表于 2017-7-15 09:59
老师，我想问一下有GRIB转ARL数据的的脚本吗？这个汇总贴没找到啊

参考这里：http://www.meteothinker.com/exam ... ctory/grib2arl.html

daibingbing

MeteoInfo 发表于 2017-7-15 15:11
参考这里：http://www.meteothinker.com/examples/meteoinfolab/trajectory/grib2arl.html

谢谢您

走路睡觉的猫WD

太给力了，谢谢分享

Swallow

meteoinfo学习贴

笨笨

{:5_213:}赞，好多python命令能用

♂雨已~○

求教王老师，如何从脚本中获取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...'

效果图：




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

MeteoInfo

♂雨已~○ 发表于 2017-8-14 13:39
求教王老师，如何从脚本中获取HDF文件的时间，并显示在图片上？代码咋个写呢，才开始接触MeteoInfoLab还不 ...

获取第一个时次的时间
f.gettime(0)
可以参考这里：http://bbs.06climate.com/forum.php?mod=viewthread&tid=43434

♂雨已~○

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))

我添加了这些代码实现了，虽然复杂了点，好像也可以了，不过王老师这直接的方法效率更高啊

♂雨已~○

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

请问王老师，这个脚本能加载web地图吗？比如说百度，Google或者高德地图的离线地图文件？如果能，那么该怎么加载呢？

MeteoInfo

♂雨已~○ 发表于 2017-8-14 17:24
请问王老师，这个脚本能加载web地图吗？比如说百度，Google或者高德地图的离线地图文件？如果能，那么该 ...

本贴1楼就有例子，你仔细看看。