请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: MeteoInfo

MeteoInfoLab脚本示例:风场矢量图

[复制链接]

新浪微博达人勋

发表于 2016-12-1 10:26:41 | 显示全部楼层
MeteoInfo 发表于 2016-11-30 22:44
你可以在百度网盘上下载:http://pan.baidu.com/s/1kUENi4n

王老师,下载了您的最新版meteoinfo,成功出图,这个时间变量(datetime.datetime)提取数组真心好用,也不用担心时间设置错误的问题了,点赞!
1111.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-12-3 17:04:30 | 显示全部楼层
MeteoInfo 发表于 2016-11-22 23:36
eta坐标差值为等压面目前还没有现成的函数。

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

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

新浪微博达人勋

 楼主| 发表于 2016-12-3 18:38:49 | 显示全部楼层
本帖最后由 MeteoInfo 于 2016-12-3 18:52 编辑
半颗麦芽糖 发表于 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
(出处: 气象家园)

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

新浪微博达人勋

发表于 2016-12-3 19:10:25 | 显示全部楼层
不错,谢谢分享!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-12-3 21:01:31 | 显示全部楼层
MeteoInfo 发表于 2016-12-3 18:38
这个不用专门的函数,按等间隔取数组就可以了。比如每4个数取一个:
u = u[::4,::4]

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

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

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

新浪微博达人勋

发表于 2016-12-3 21:31:36 | 显示全部楼层
对了,王老师,我用欧洲中心的nc格式资料按照这个脚本画并不存在这样的问题,还有我的模式结果在meteoinfo桌面软件中画的时候也不存在这样的问题。所以我在想是不是投影或者坐标的问题?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

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

数据肯定是有投影的,没看出来有什么问题呀?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-12-4 00:11:35 | 显示全部楼层
MeteoInfo 发表于 2016-12-3 22:30
数据肯定是有投影的,没看出来有什么问题呀?

王老师,青藏高原上下纬度在25到45之间,但是这图是20到40之间,纬度存在存在,或者说纵坐标存在问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-12-4 10:25:43 | 显示全部楼层
半颗麦芽糖 发表于 2016-12-4 00:11
王老师,青藏高原上下纬度在25到45之间,但是这图是20到40之间,纬度存在存在,或者说纵坐标存在问题

在有投影的情况下经纬线不是直线,是弯曲的。你可以把经纬线画出来看看:
axesm(projinfo=f.proj, gridline=True)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-12-4 10:36:45 | 显示全部楼层
本帖最后由 半颗麦芽糖 于 2016-12-4 10:55 编辑
MeteoInfo 发表于 2016-12-4 10:25
在有投影的情况下经纬线不是直线,是弯曲的。你可以把经纬线画出来看看:
axesm(projinfo=f.proj, gridl ...


王老师,我按照您说的加了,果然是我搞错了,他采用了兰伯特投影,所以等纬度线并不是直线,老师我想问一下,如果我想把图改成正常经纬度可否实现,因为是用来作对比的,加了投影还是有点差异。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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