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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 17312|回复: 43

MeteoInfoLab脚本示例:CMORPH降水数据

[复制链接]

新浪微博达人勋

发表于 2016-12-31 16:01:23 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 MeteoInfo 于 2017-1-19 09:50 编辑

NOAA CPC的CMORPH降水数据(http://www.cpc.ncep.noaa.gov/pro ... ph_description.html)是二进制格式的格点数据,国家气象信息中心又融合了自动站的降水数据形成了降水融合产品(http://data.cma.cn/data/detail/d ... HOUR_GRID_0.10.html),同样是二进制格点数据。两种数据都有相应的ctl文件供GrADS软件读取处理数据。不过这里示例用MeteoInfoLab的binread()函数读取此类数据,并绘制降水分布图。

NOAA CPC CMORPH降水数据:
  1. fn = 'D:/Temp/binary/CMORPH+MWCOMB_3HRLY-025DEG_20160714'
  2. nvar = 2
  3. ntime = 8
  4. nlat = 480
  5. nlon = 1440
  6. data = binread(fn, [nvar*ntime,nlat,nlon], 'float', byteorder='big_endian')
  7. comb = data[::2,::-1,:]
  8. cpc = data[1::2,::-1,:]
  9. rain = cpc[2,:,:]
  10. lon = arange1(0.125, nlon, 0.25)
  11. lat = arange1(-59.875, nlat, 0.25)
  12. #Plot
  13. axesm()
  14. lworld = shaperead('D:/Temp/Map/country1.shp')
  15. lchina = shaperead('D:/Temp/Map/bou2_4p.shp')
  16. geoshow(lchina, edgecolor='k')
  17. geoshow(lworld, edgecolor='k')
  18. levs = [0.1,0.5,1,2,3,4,5,6,8,10,20,40]
  19. cols = makecolors(len(levs) + 1, cmap='BlAqGrYeOrReVi200')
  20. cols[0] = 'w'
  21. layer = imshowm(lon, lat, rain, levs, colors=cols)
  22. colorbar(layer, label='mm')
  23. title('Precipitation - CMORPH (2016-07-14 06:00)')
  24. #xlim(70, 140)
  25. #ylim(15, 55)


rain_cmorph_global.png

缩放到中国范围:
rain_cmorph.png

国内降水融合数据
  1. fn = 'D:/Temp/grads/SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2016071406.grd'
  2. data = binread(fn, [2, 440, 700], 'float')
  3. crain = data[0,:,:]
  4. crain[crain==-999.0] = nan
  5. gsamp = data[1,:,:]
  6. x = arange1(70.05, 700, 0.1)
  7. y = arange1(15.05, 440, 0.1)
  8. axesm()
  9. lworld = shaperead('D:/Temp/map/country1.shp')
  10. lchina = shaperead('D:/Temp/map/bou2_4p.shp')
  11. geoshow(lworld, edgecolor='k')
  12. geoshow(lchina, edgecolor='k')
  13. levs = [0.1,0.5,1,2,3,4,5,6,8,10,20,40]
  14. cols = makecolors(len(levs) + 1, cmap='BlAqGrYeOrReVi200')
  15. cols[0] = 'w'
  16. layer = imshowm(x, y, crain, levs, colors=cols)
  17. #layer = imshowm(x, y, gsamp)
  18. title('Precipitation - CHN_MERGE (2016-07-14 06:00)')
  19. colorbar(layer, label='mm')
  20. xlim(70, 140)
  21. ylim(15, 55)


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

新浪微博达人勋

发表于 2017-4-17 20:59:32 | 显示全部楼层
我导师给了某年的cmorph daily降水资料,但是全是dat格式的文件,这需要怎么处理呢?用grads画不出来,好慌
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2016-12-31 21:44:52 | 显示全部楼层
支持王老师!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-1-3 22:23:18 | 显示全部楼层
哇,这个真的很有用啊,谢谢分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-1-19 10:44:44 | 显示全部楼层
楼主,我想请问一下,现在从哪可以下到2010年和2011年的数据呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-4-17 23:42:13 | 显示全部楼层
lichangcan 发表于 2017-4-17 20:59
我导师给了某年的cmorph daily降水资料,但是全是dat格式的文件,这需要怎么处理呢?用grads画不出来,好慌

1楼不是有脚本程序吗
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-4-18 14:25:34 | 显示全部楼层
MeteoInfo 发表于 2017-4-17 23:42
1楼不是有脚本程序吗

数据格式不同的,程序操作是一样的吗
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-4-18 14:35:37 | 显示全部楼层
lichangcan 发表于 2017-4-18 14:25
数据格式不同的,程序操作是一样的吗

基本是一样的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-5-5 16:42:46 | 显示全部楼层
王老师,国内降水融合数据(1小时累计量)的时间单位是什么呀?UTC还是北京时呀?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-5-6 14:37:40 | 显示全部楼层
小1龙2虾3 发表于 2017-5-5 16:42
王老师,国内降水融合数据(1小时累计量)的时间单位是什么呀?UTC还是北京时呀?

我也不清楚,你问数据制作者吧。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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