更新：已解决，改为经验，CartoPy画色斑图的时候无法同时显示地图和数据

QCD

更新：借助Google找到了解决方法
StackOverFlow上已经有人提过类似的问题（https://stackoverflow.com/questions/45042722/cartopy-matplotlib-contourf-map-overriding-data），是由于两幅图都是同一个类型，使用的zorder是相同的，所以需要手动制定zorder，可以参考matplotlib关于zorder的demo（https://matplotlib.org/examples/pylab_examples/zorder_demo.html）
ax.add_feature(adm2_shapes,zorder=1)
axis = ax.contourf(olon,olat,rain_data_new, levels=levels, cmap='viridis',transform=ccrs.PlateCarree(),zorder=5)

最终结果

最近在学习使用Python进行色斑图，在学习了论坛里多位大神的经验后（如非对称的帖子1和他在Github上的MicapsDraw；平流层的萝卜的maskout方法，也参考过chiqu296的code）成功使用Basemap画出的想要的图。

但是考虑到Basemap势必逐渐被CartoPy取代，MetPy和Iris等的例子里面也多使用CartoPy，因此决定尝试CartoPy，在学习了po_po1的方法后实现了成功剪切，如

但是随后在叠加地图的时候遇到了麻烦，如果使用ax.add_feature来加入用于裁剪的shp文件，就会发生只有地图而没有色斑图的现象，如下图所示

最后发现如果将读入的地图数据的alpha值降低后就可以将2者都显示出来，但是色斑图也被叠加了alpha值处理，如下图所示

不知道社区的大神们有没有遇到过类似的问题，或者有没有其他使用cartopy来画色斑图的成功经验？谢谢！


附上使用的代码：

1. import numpy as np
   
2. from scipy.interpolate import Rbf
   
3. import pandas as pd
   
4. 
   
5. import matplotlib.pyplot as plt
   
6. from matplotlib.pyplot import colormaps as cmap
   
7. 
   
8. import cartopy.crs as ccrs
   
9. import cartopy.feature as cfeat
   
10. 
    
11. from cartopy.io.shapereader import Reader
    
12. from cartopy.feature import ShapelyFeature
    
13. from cartopy.mpl.patch import geos_to_path
    
14. 
    
15. from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter

复制代码

1. <font face="微软雅黑">data = pd.read_csv('../rain05.csv')
   
2. 
   
3. # 插值
   
4. lon = data['Lon']
   
5. lat = data['Lat']
   
6. rain_data = data['PRE_Time_0808']
   
7. 
   
8. # [103,110,24.5,29.3]
   
9. olon = np.linspace(103,110,88)
   
10. olat = np.linspace(24.5,29.3,88)
    
11. olon,olat = np.meshgrid(olon,olat)
    
12. # 插值处理
    
13. func = Rbf(lon, lat, rain_data,function='linear')
    
14. rain_data_new = func(olon, olat)
    
15. 
    
16. 
    
17. 
    
18. levels = [0,10,25,50,100,150,200,300,450]
    
19. 
    
20. 
    
21. fname = '../Micaps-shapefiles/Province.shp'
    
22. fname2 = '../Micaps-shapefiles/Province.shp'
    
23. adm1_shapes = ShapelyFeature(Reader(fname).geometries(),ccrs.PlateCarree(), facecolor='White',edgecolor='black')
    
24. 
    
25. adm2_shapes = ShapelyFeature(Reader(fname2).geometries(),ccrs.PlateCarree(),facecolor='White',edgecolor='black',alpha=0.1,linewidth=2)
    
26. 
    
27. 
    
28. #设定裁剪区域
    
29. reader = Reader(fname)
    
30. provinces = reader.records()
    
31. 
    
32. 
    
33. #Micaps地图中Province.shp中可以使用PAC区分不同省份，同理City.shp和County.shp也可使用该关键字
    
34. us_multipoly, = [province.geometry for province in provinces if province.attributes['PAC'] == 520000]
    
35. main_us_geom = sorted(us_multipoly.geoms, key=lambda geom: geom.area)[0]
    
36. us_path, = geos_to_path(main_us_geom)
    
37. 
    
38. #画图部分
    
39. fig = plt.figure(figsize=(12,9))
    
40. 
    
41. ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())
    
42. plate_carre_data_transform = ccrs.PlateCarree()._as_mpl_transform(ax)
    
43. 
    
44. axis = ax.contourf(olon,olat,rain_data_new, levels=levels, cmap='viridis',transform=ccrs.PlateCarree(),zorder=5)
    
45. 
    
46. for collection in axis.collections:
    
47.     collection.set_clip_path(us_path, plate_carre_data_transform)#设置显示区域
    
48. 
    
49. 
    
50. ax.add_feature(adm2_shapes,zorder=1)
    
51. 
    
52. # ax.set_extent([103,110,24.5,29.3])</font>

复制代码

nickbsb

QCD 发表于 2018-7-12 10:33
地图的问题已经解决了，一个是制定zorder，另一个是加上alpha=1，我没搞懂后面这种为什么也可以起作用

请教下楼主zorder的作用是什么？为什么设置为1和5呢？

QCD

不好意思，代码里漏掉了一部分

1. import numpy as np
   
2. from scipy.interpolate import Rbf
   
3. import pandas as pd
   
4. 
   
5. import matplotlib.pyplot as plt
   
6. from matplotlib.pyplot import colormaps as cmap
   
7. 
   
8. import cartopy.crs as ccrs
   
9. import cartopy.feature as cfeat
   
10. 
    
11. from cartopy.io.shapereader import Reader
    
12. from cartopy.feature import ShapelyFeature
    
13. from cartopy.mpl.patch import geos_to_path
    
14. 
    
15. from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter

import numpy as np<br /> from scipy.interpolate import Rbf<br /> import pandas as pd<br /> <br /> import matplotlib.pyplot as plt<br /> from matplotlib.pyplot import colormaps as cmap<br /> <br /> import cartopy.crs as ccrs<br /> import cartopy.feature as cfeat<br /> <br /> from cartopy.io.shapereader import Reader<br /> from cartopy.feature import ShapelyFeature<br /> from cartopy.mpl.patch import geos_to_path<br /> <br /> from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter

⑨P_bō

你好贵州的同事，我最近也在画色板图

Pilgrim

我想了解插值，线性插值不通用

QCD

⑨P_bō 发表于 2018-7-7 23:27
你好贵州的同事，我最近也在画色板图

卫星遥感的同事？看ID似乎能关联到名字？

⑨P_bō

QCD 发表于 2018-7-8 13:55
卫星遥感的同事？看ID似乎能关联到名字？

哈哈，是的

⑨P_bō

QCD 发表于 2018-7-8 13:55
卫星遥感的同事？看ID似乎能关联到名字？

哈哈，是的

chongzika

楼主现在的代码是已经解决了问题的么？

QCD

chongzika 发表于 2018-7-11 23:15
楼主现在的代码是已经解决了问题的么？

地图的问题已经解决了，一个是制定zorder，另一个是加上alpha=1，我没搞懂后面这种为什么也可以起作用