配合 cartopy 使用的中国地图包：frykit

灭火器

项目仓库：https://github.com/ZhaJiMan/frykit
详细介绍：https://zhajiman.github.io/post/frykit/

这个包写了有很久了，最近更新了一波，发到家园宣传一下，有人用才能找出 bug，希望能多得到反馈。之前有个单词拼写错误的 bug 硬是半年才被别人发现……
功能是：

- 自带中国省市县的行政区划数据，无需准备 shapefile 或者 geojson 文件
- 0.7.0 开始自带两套数据，一套高德地图 api 的，一套天地图去年公开的，任选。
- 一个函数快速绘制，快速标注
- 地图白化、指北针、比例尺等小工具

用法类似下面：

1. import matplotlib.pyplot as plt
   
2. import cartopy.crs as ccrs
   
3. import frykit.plot as fplt
   
4. 
   
5. crs = ccrs.PlateCarree()
   
6. ax = plt.axes(projection=crs)
   
7. 
   
8. fplt.add_cn_border(ax)  # 添加国界
   
9. fplt.add_cn_province(ax)  # 添加省界
   
10. fplt.add_cn_city(ax)  # 添加市界
    
11. fplt.add_cn_district(ax)  # 添加县界
    
12. fplt.add_cn_line(ax)  # 添加九段线
    
13. 
    
14. fplt.label_cn_province(ax)  # 添加省名
    
15. fplt.label_cn_city(ax)  # 添加市名
    
16. fplt.label_cn_district(ax)  # 添加县名
    
17. 
    
18. # 用国界裁剪填色图
    
19. cf = ax.contourf(lon, lat, data, transform=crs)
    
20. fplt.clip_by_cn_border(cf)
    
21. 
    
22. fplt.add_compass(ax, 0.1, 0.9)  # 添加指北针
    
23. fplt.add_scale_bar(ax, 0.1, 0.1)  # 添加比例尺

复制代码

效果图如下


最新的 0.7.x 版本安装方式是：

1. pip install frykit[data]

复制代码

conda 包还不会弄，以后有空学学。

灭火器

gdou_lyc 发表于 2025-8-9 10:01
感谢大神提供的工具！！想问问大神，我南海小图里contouf填色是没问题的，但是用quiver画风矢量不显示，想 ...

1. from typing import cast
   
2. 
   
3. import frykit.plot as fplt
   
4. import matplotlib.pyplot as plt
   
5. import numpy as np
   
6. from cartopy.mpl.geoaxes import GeoAxes
   
7. from scipy.ndimage import gaussian_filter
   
8. 
   
9. # 读取数据
   
10. data = fplt.load_test_data()
    
11. X, Y = np.meshgrid(data.lon, data.lat)
    
12. t2m = gaussian_filter(data.t2m - 273.15, sigma=1)
    
13. u10 = gaussian_filter(data.u10, sigma=1)
    
14. v10 = gaussian_filter(data.v10, sigma=1)
    
15. 
    
16. # 设置投影
    
17. map_crs = fplt.CN_AZIMUTHAL_EQUIDISTANT
    
18. data_crs = fplt.PLATE_CARREE
    
19. 
    
20. # 准备地图
    
21. fig = plt.figure(figsize=(8, 5))
    
22. main_ax = fig.add_subplot(projection=map_crs)
    
23. mini_ax = fplt.add_mini_axes(main_ax)
    
24. main_ax = cast(GeoAxes, main_ax)
    
25. mini_ax = cast(GeoAxes, mini_ax)
    
26. main_ax.set_extent((74, 136, 13, 57), crs=data_crs)
    
27. mini_ax.set_extent((105, 122, 2, 25), crs=data_crs)
    
28. for ax in [main_ax, mini_ax]:
    
29.     fplt.add_cn_province(ax, lw=0.25)
    
30.     fplt.add_cn_line(ax)
    
31. 
    
32. # 绘制气温
    
33. for i, ax in enumerate([main_ax, mini_ax]):
    
34.     cf = ax.contourf(
    
35.         X,
    
36.         Y,
    
37.         t2m,
    
38.         levels=np.linspace(-10, 35, 10),
    
39.         cmap="plasma",
    
40.         extend="both",
    
41.         transform=data_crs,
    
42.         transform_first=True,
    
43.     )
    
44.     fplt.clip_by_cn_border(cf)
    
45.     if i == 0:
    
46.         fig.colorbar(cf, ax=ax, label="Temperature (℃)")
    
47. 
    
48. # 用 stride 替代 regrid_shape
    
49. stride = 6
    
50. s = np.s_[::stride, ::stride]
    
51. 
    
52. # 绘制风场
    
53. # 用 inches 单位保证大图小图里的箭头大小相同
    
54. quiver_kwargs = dict(
    
55.     x=X[s],
    
56.     y=Y[s],
    
57.     u=u10[s],
    
58.     v=v10[s],
    
59.     scale=30,
    
60.     scale_units="inches",
    
61.     units="inches",
    
62.     width=0.01,
    
63.     headwidth=4,
    
64.     transform=data_crs,
    
65. )
    
66. main_Q = main_ax.quiver(**quiver_kwargs)
    
67. mini_Q = mini_ax.quiver(**quiver_kwargs)
    
68. fplt.clip_by_cn_border(main_Q)
    
69. fplt.clip_by_cn_border(mini_Q)
    
70. fplt.add_quiver_legend(main_Q, U=10, height=0.12, loc="lower left")
    
71. 
    
72. fplt.savefig("quiver.png")
    
73. plt.close(fig)
    

复制代码

1. matplotlib==3.10.5
   
2. cartopy==0.25.0
   
3. frykit==0.7.5.post1

复制代码

我这里画出来是正常的，就直接画

2417404759

太好了，就是感激像楼主这样的人，之前我使用cartopy时对于它里面的国界定义标注错误不爽，一直以来只画海岸线，要更细还得用shp

御风者

很好的工作，如果方便，可以上传到conda-forge 源中。

安达充

好东西，dd

Lancelot

是炸鸡人，哈哈

一大碗年糕

贫道敬孔

大佬太强了，必须赞，必须顶！！

yqcxm

用过，很好用，没有长江、黄河等湖泊和河流数据

万里路始于足下

炸鸡人大佬我是认得，你的教程真的牛！

hazuki91

插眼，感谢大佬了