MeteoInfoLab脚本示例：风场矢量图

半颗麦芽糖

MeteoInfo 发表于 2016-11-30 22:44
你可以在百度网盘上下载：http://pan.baidu.com/s/1kUENi4n

王老师，下载了您的最新版meteoinfo，成功出图，这个时间变量（datetime.datetime）提取数组真心好用，也不用担心时间设置错误的问题了，点赞！

半颗麦芽糖

MeteoInfo 发表于 2016-11-22 23:36
eta坐标差值为等压面目前还没有现成的函数。

如果风场矢量太密的话，可以先将U, V分量数组抽稀，然后 ...

王老师，这两天一直在考虑这个U/V风场抽稀的问题，找了一圈好像meteoinfo并没有类似grads里面skip的用法，然后想着用气象家园里面说到的grads的re函数来实现，但是也出了很多问题，所以想请问王老师一下meteoinfo脚本能否实现这个，该如何实现？万分感谢

MeteoInfo

半颗麦芽糖 发表于 2016-12-3 17:04
王老师，这两天一直在考虑这个U/V风场抽稀的问题，找了一圈好像meteoinfo并没有类似grads里面skip的用法 ...

这个不用专门的函数，按等间隔取数组就可以了。比如每4个数取一个：
u = u[::4,::4]

参考这个帖子：MeteoInfoLab脚本示例：Hamawari-8 netCDF data
http://bbs.06climate.com/forum.p ... d=37038&fromuid=106
(出处: 气象家园)

在MeteoInfo桌面软件中也有类似skip的功能，参加此贴：关于UV矢量图形
http://bbs.06climate.com/forum.p ... 742&fromuid=106
(出处: 气象家园)

xingkongyao

不错，谢谢分享！

半颗麦芽糖

MeteoInfo 发表于 2016-12-3 18:38
这个不用专门的函数，按等间隔取数组就可以了。比如每4个数取一个：
u = u[::4,::4]

王老师，我根据您的指导对u、v变量进行了抽稀，但是出现了一个问题，地图和数据的位置发生了偏移，经度没有发生变化，但是地图和数据同时向正南方偏了5个纬度，脚本和图片如下，请老师指导

1. 
   
2. <P> f = addfile('E:/gongzuo/3.8.1/test38101.ctl')
   
3. t = datetime.datetime(2014,10,26,18)
   
4. u = f['U'][[t],[500],::4,::4]
   
5. v = f['V'][[t],[500],::4,::4]
   
6. speed = sqrt(u*u+v*v)
   
7. #Plot
   
8. axesm(projinfo=f.proj)
   
9. lchina = shaperead('F:/11/MeteoInfo/map/bou2_4p.shp')
   
10. ltibet = shaperead('F:/11/MeteoInfo/map/DBATP/tibet_ASCII.wmp')
    
11. geoshow(lchina, edgecolor='k')
    
12. geoshow(ltibet, edgecolor='r',size=2)
    
13. layer = quiverm(u,v,speed,5,size=5,proj=f.proj)
    
14. colorbar(layer,fontsize=17)
    
15. #title(t.strftime('%Y-%m-%d %H:00'))
    
16. xlabel(u'2014年10月26日18时(UTC)500hPa风场',fontname='黑体',fontsize=17)
    
17. xticks(fontsize=17)
    
18. yticks(fontsize=17)</P>

复制代码

半颗麦芽糖

对了，王老师，我用欧洲中心的nc格式资料按照这个脚本画并不存在这样的问题，还有我的模式结果在meteoinfo桌面软件中画的时候也不存在这样的问题。所以我在想是不是投影或者坐标的问题？

MeteoInfo

半颗麦芽糖 发表于 2016-12-3 21:01
王老师，我根据您的指导对u、v变量进行了抽稀，但是出现了一个问题，地图和数据的位置发生了偏移，经度没 ...

数据肯定是有投影的，没看出来有什么问题呀？

半颗麦芽糖

MeteoInfo 发表于 2016-12-3 22:30
数据肯定是有投影的，没看出来有什么问题呀？

王老师，青藏高原上下纬度在25到45之间，但是这图是20到40之间，纬度存在存在，或者说纵坐标存在问题

MeteoInfo

半颗麦芽糖 发表于 2016-12-4 00:11
王老师，青藏高原上下纬度在25到45之间，但是这图是20到40之间，纬度存在存在，或者说纵坐标存在问题

在有投影的情况下经纬线不是直线，是弯曲的。你可以把经纬线画出来看看：
axesm(projinfo=f.proj, gridline=True)

半颗麦芽糖

MeteoInfo 发表于 2016-12-4 10:25
在有投影的情况下经纬线不是直线，是弯曲的。你可以把经纬线画出来看看：
axesm(projinfo=f.proj, gridl ...

王老师，我按照您说的加了，果然是我搞错了，他采用了兰伯特投影，所以等纬度线并不是直线，老师我想问一下，如果我想把图改成正常经纬度可否实现，因为是用来作对比的，加了投影还是有点差异。