Python绘制CMORPH降水

付亚男

numpy读取数据
proplot绘图

1. """
   
2. Created on Wed Nov 10 19:05:06 2021
   
3. 
   
4. @author: Fu Ya-Nan
   
5. """
   
6. 
   
7. import numpy as np
   
8. import proplot as pplt
   
9. import cartopy.crs as ccrs
   
10. import cartopy.feature as cfeat
    
11. from cartopy.io.shapereader import Reader
    
12. 
    
13. ntime = 2
    
14. nlon = 4948
    
15. nlat = 1649
    
16. fnm = r'F:/cmorph/CMORPH_V1.0_ADJ_8km-30min_2019063023'
    
17. with open(fnm, 'rb') as f:
    
18.     prerate = np.fromfile(f, dtype=np.float32)
    
19. prerate = np.reshape(prerate, (ntime, nlat, nlon))
    
20. preci = prerate[0, :, :]*0.5
    
21. preci[preci<0.0] = 0.0
    
22. 
    
23. lonstart = 0.036378335
    
24. londiff = 0.072756669
    
25. latstart = -59.963614
    
26. latdiff = 0.072771377
    
27. lon = np.arange(nlon)*londiff+lonstart
    
28. lat = np.arange(nlat)*latdiff+latstart
    
29. 
    
30. fig = pplt.figure(refwidth = 5.0)
    
31. ax = fig.subplots(proj = 'cyl')
    
32. cf = ax.contourf(lon[1305:1718], lat[962:1237], preci[962:1237, 1305:1718])
    
33. ax.colorbar(cf)
    
34. ax.format(lonlim = (95, 125), latlim = (10,  30), lonlines = 5, latlines = 5, lonlabels = 'b', latlabels = 'l')
    
35. provinces = cfeat.ShapelyFeature( Reader(r'F:/ngcc/bou2_4m/bou2_4l.shp').geometries(), ccrs.PlateCarree(), edgecolor='black', facecolor='none', )
    
36. river = cfeat.ShapelyFeature( Reader(r'F:/Chinamap-master/cnmap/rivers.shp').geometries(), ccrs.PlateCarree(), edgecolor='lightblue', facecolor='none', )
    
37. ax.add_feature(provinces, linewidth=0.5, zorder=2)
    
38. ax.add_feature(river, linewidth=1.0, zorder=2)
    
39. fig.format(lefttitle = 'CMORPH_20190630_23:00-23:30')

复制代码

gaoyanbo22

请问你的这个数据下载地址有吗？有半小时的吗

gaoyanbo22

preci = prerate[0, :, :]*0.5这里乘以0.5是把降水值处理为一半进行画图的吗

鲸鱼便子

请教一下楼主，这个分辨率半小时的降水单位是 mm/h, 不太理解，是瞬时降水率吗？
如果想要得到多小时的累加降水，应该怎么处理？