爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 11938|回复: 20

MeteoInfoLab脚本示例:气旋散点图

[复制链接]

新浪微博达人勋

发表于 2015-7-6 21:12:41 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
全球长时间序列中尺度气旋数据(http://cioss.coas.oregonstate.edu/eddies/)有netCDF格式,散点数据类型,只有一个很大的维Nobs = 2590938。尝试读取了部分经纬度和track数据,并绘制散度图。

脚本程序:
  1. fn = 'D:/Temp/nc/tracks.20130125.nc'
  2. f = addfile(fn)
  3. n = 1000
  4. lon = f['lon'][0:n]
  5. lat = f['lat'][0:n]
  6. track = f['track'][0:n]
  7. axesm()
  8. mlayer = shaperead('D:/Temp/map/country1.shp')
  9. geoshow(mlayer, edgecolor=(0,0,255))
  10. layer = scatterm(lon, lat, track)
  11. title('Eddies')
  12. colorbar(layer, orientation='horizontal')
  13. axism([-180, 180, -90, 90])


Image00876.png

评分

参与人数 1金钱 +2 收起 理由
温骸 + 2 很给力!

查看全部评分

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

新浪微博达人勋

发表于 2015-7-7 10:50:18 | 显示全部楼层
谢谢,谢谢!
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2015-7-7 16:15:47 | 显示全部楼层
麻烦问一下,为什么我不能使用addfile这个命令,是不是有什么工具包
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-7-7 17:02:01 | 显示全部楼层
温骸 发表于 2015-7-7 16:15
麻烦问一下,为什么我不能使用addfile这个命令,是不是有什么工具包

你需要看清楚,不是MatLab,是MeteoInfoLab。不同的软件。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-7-7 19:26:39 | 显示全部楼层
MeteoInfo 发表于 2015-7-7 17:02
你需要看清楚,不是MatLab,是MeteoInfoLab。不同的软件。

哈搜的,谢谢老师!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-3-28 08:46:05 | 显示全部楼层
看了一下,都是海洋上气旋/反气旋,请问对于陆地上的气旋/反气旋,该怎么做呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2021-3-12 16:38:20 | 显示全部楼层
尊敬的王老师,我用layer = scatterm(lon, lat, aod)这个语句想画出通过显著性检验的点,但是出来的点是填充了整个范围的点,请教如何解决这个问题,不胜感激!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-3-12 16:42:00 | 显示全部楼层
王老师,我参考您的脚本,编了如下,但画不出来图,如果layer = scatterm(sig)这样,可以出来点,但是点填满了整个范围,请老师指点,感谢~

f = addfile('J:\dbandant\c_bba.nc')
v_lwp = f['Cloud_Water_Path_Liquid_Mean']
v_sig = f['corr_sig']
lon = f['lon'][:]
lat = f['lat'][:]
lwp = v_lwp['35:53','101:134']
sig = v_sig['35:53','101:134']
axesm()
china = shaperead('D:\MeteoInfo_2.3.1\MeteoInfo\map\china.shp')
geoshow(china, edgecolor='lightgray')
procince = shaperead('D:\MeteoInfo_2.3.1\MeteoInfo\map\cn_province.shp')
geoshow(procince, edgecolor='lightgray')
layer = scatterm(lon, lat, sig)
title('LWP')
colorbar(layer)
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-3-12 19:55:56 | 显示全部楼层
晴天 发表于 2021-3-12 16:42
王老师,我参考您的脚本,编了如下,但画不出来图,如果layer = scatterm(sig)这样,可以出来点,但是点填 ...

不太明白你到底想问什么
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-3-13 07:01:32 | 显示全部楼层
MeteoInfo 发表于 2021-3-12 19:55
不太明白你到底想问什么

是这样王老师,我做了云水路径LWP和气溶胶光学厚度在一定范围的相关分析,然后把通过显著性检验的sig点在相关分析图上
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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