Cartopy 下站点资料画图及白化

chiqu296

   之前发过Basemap下画站点数据的帖子，现在转到Cartopy下，所以重新写了一个，白化使用平流层的萝卜的Python完美白化，只是在投影转换那里稍作了改动。

1. import os
   
2. import maskout
   
3. import numpy as np
   
4. import xarray as xr
   
5. from scipy.interpolate import Rbf
   
6. import matplotlib as mpl
   
7. import matplotlib.pyplot as plt
   
8. import cartopy.crs as ccrs
   
9. import cartopy.feature as cfeature
   
10. from cartopy.io.shapereader import Reader
    
11. 
    
12. SHP = './shapefile/china_shp'
    
13. 
    
14. # 数据处理
    
15. ds = xr.open_dataset('./r160.nc')
    
16. da = ds.PRE.sel(month=slice(6, 8)).mean(dim='month').mean(dim='year')
    
17. 
    
18. # 插值
    
19. olon = np.linspace(72, 137, 260)
    
20. olat = np.linspace(15, 55, 120)
    
21. olon, olat = np.meshgrid(olon, olat)
    
22. func = Rbf(ds.station_lon, ds.station_lat, da, function='linear')
    
23. pre_grid = func(olon, olat)
    
24. 
    
25. # 字体
    
26. mpl.rcParams['font.sans-serif'] = ['Times New Roman']
    
27. mpl.rc('font', size=15, weight='normal')
    
28. font = {'family': 'NSimSun','weight': 'normal','size': 15}
    
29. 
    
30. # 颜色
    
31. levels = [0, 0.1, 10, 25, 50, 100, 250, 500]
    
32. cmaps = mpl.colors.ListedColormap(['#FFFFFF', '#A9F18D', '#3DB93D', '#60B8FF', '#0001E2', '#F900FA', '#810045'])
    
33. norm = mpl.colors.BoundaryNorm(levels, 7)
    
34. 
    
35. # 画图
    
36. proj = ccrs.LambertConformal(central_latitude=90, central_longitude=105)
    
37. fig, ax = plt.subplots(figsize=(15, 15), subplot_kw=dict(projection=proj))
    
38. ax.set_extent([80, 130, 13, 55], ccrs.PlateCarree())
    
39. ax.gridlines(linestyle='--')
    
40. 
    
41. ax.add_feature(cfeature.OCEAN.with_scale('50m'))
    
42. ax.add_feature(cfeature.LAND.with_scale('50m'))
    
43. ax.add_feature(cfeature.RIVERS.with_scale('50m'))
    
44. ax.add_feature(cfeature.LAKES.with_scale('50m'))
    
45. 
    
46. ax.add_geometries(Reader(os.path.join(SHP, 'cnhimap.shp')).geometries(),ccrs.PlateCarree(),facecolor='none',edgecolor='k',linewidth=1)
    
47. 
    
48. cf = ax.contourf(olon,olat,pre_grid,levels=levels,cmap=cmaps,norm=norm,zorder=0,transform=ccrs.PlateCarree())
    
49. ax.scatter(ds.station_lon,ds.station_lat,c='k',s=10,marker='o',transform=ccrs.PlateCarree())
    
50. 
    
51. cbar = fig.colorbar(cf,ax=ax,orientation="vertical",aspect=25,pad=0.01,shrink=0.7)
    
52. cbar.set_label('毫米', fontdict=font)
    
53. cbar.set_ticklabels(levels)
    
54. ax.set_title('多年平均夏季降水', fontdict=font)
    
55. 
    

复制代码

muggle

好东西，谢啦

Eegle

很好的白化例子，点赞。思路清晰明了

werewolf

白化是哪部分，没看到

wolfwang

老师，用你的方法，是用的cartopy0.17版本，是不是哪个包不兼容呢？

  File "D:\py\nc\maskout.py", line 32, in shp2clip
    contour.set_clip_path(clip)

UnboundLocalError: local variable 'clip' referenced before assignment

hely8586

学习学习！！

珮瑶瑶

werewolf 发表于 2019-12-30 22:33
白化是哪部分，没看到

我也没有看明白，白化在哪里？

zqf_23

近期要学习用站点资料画全国月平均气温，月降水量图，希望能用上。点我查看给积分不够的朋友的建议，只下了
maskout.py，还没跑通。

zqf_23

用其他下载的地图，运行报错  baowarnings.warn('Tight layout not applied. The left and right margins '   出来的图不对，请问是怎么回事？

zqf_23

('Tight layout not applied. The left and right margins）报错已解决，但为什么颜色填充没有在中国地图内？感觉是个扇形