如何更优雅地绘制/白化中国的地图

clarmy

gaorui123 发表于 2023-4-4 17:22
您好，我使用那个maskout函数进行裁剪时，我使用的ax.pcolormesh而不是ax.contourf。然后就会报错： 'GeoQu ...

maskout 函数应该与 contourf 或 pcolormesh 无关的。

1. import numpy as np
   
2. from cnmaps import get_adm_maps
   
3. import matplotlib.pyplot as plt
   
4. 
   
5. lon = np.linspace(60, 150, 1000)
   
6. lat = np.linspace(0, 60, 1000)
   
7. 
   
8. lons, lats = np.meshgrid(lon, lat)
   
9. data = np.random.random(lons.shape)
   
10. 
    
11. china = get_adm_maps(level="国", record= "first", only_polygon=True, wgs84=True)
    
12. 
    
13. maskout_data = china.maskout(lons, lats, data)
    
14. 
    
15. plt.figure(figsize=(20,8))
    
16. 
    
17. plt.subplot(121)
    
18. plt.pcolormesh(lons, lats, data)
    
19. plt.title("no maskout")
    
20. 
    
21. plt.subplot(122)
    
22. plt.pcolormesh(lons, lats, maskout_data)
    
23. plt.title("maskout")
    
24. plt.show()

复制代码

示例：https://cnmaps.readthedocs.io/zh_CN/latest/content/usage.html

guoaxing

感谢您的经验和学习分享

guoaxing

牛啊，膜拜

曾璨_stella

gaorui123 发表于 2023-4-4 17:22
您好，我使用那个maskout函数进行裁剪时，我使用的ax.pcolormesh而不是ax.contourf。然后就会报错： 'GeoQu ...

把collections改成collections.abc

jinglj2001

我用是SURFER做白化底图文件的（后缀名为BLN)

啸BEST

楼主好，我在调用draw_maps函数绘图时弹出了'MultiPolygon' object is not iterable的报错，无法直接对MULTIPOLYGON格式数据直接进行循环。尝试将for geomestry in map_polygon:修改为for geomestry in map_polygon.geoms后问题似乎得到解决。