MeteoInfoLab脚本示例：CMORPH降水数据

MeteoInfo

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)

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

缩放到中国范围：


国内降水融合数据

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)

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

lichangcan

我导师给了某年的cmorph daily降水资料，但是全是dat格式的文件，这需要怎么处理呢？用grads画不出来，好慌

river

支持王老师！

delphi

哇，这个真的很有用啊，谢谢分享

曲梓祎

楼主，我想请问一下，现在从哪可以下到2010年和2011年的数据呢？

MeteoInfo

lichangcan 发表于 2017-4-17 20:59
我导师给了某年的cmorph daily降水资料，但是全是dat格式的文件，这需要怎么处理呢？用grads画不出来，好慌

1楼不是有脚本程序吗

lichangcan

MeteoInfo 发表于 2017-4-17 23:42
1楼不是有脚本程序吗

数据格式不同的，程序操作是一样的吗

MeteoInfo

lichangcan 发表于 2017-4-18 14:25
数据格式不同的，程序操作是一样的吗

基本是一样的

小1龙2虾3

王老师，国内降水融合数据（1小时累计量）的时间单位是什么呀？UTC还是北京时呀？

MeteoInfo

小1龙2虾3 发表于 2017-5-5 16:42
王老师，国内降水融合数据（1小时累计量）的时间单位是什么呀？UTC还是北京时呀？

我也不清楚，你问数据制作者吧。