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

付亚男

好啊，就是建议换个名字

后面用的人多了，大家都会好奇frykit是什么意思
总不能告诉别人，是炸鸡人的【炸鸡工具箱】哈哈哈

灭火器

付亚男 发表于 2025-3-16 22:01
好啊，就是建议换个名字

后面用的人多了，大家都会好奇frykit是什么意思

被你发现了

尽头的尽头

很不错，我记得天地图好像还给了审图号，这样大家写论文不愁审图号了

junyang517

感谢分享~~~~~~~~~~~~~

astiny

炸鸡老师好强大

独孤酒见

灭火器大神这下应该把关于地图常用到的功能不但实现了，而且是便捷了。大神，在您的主贴中的效果图中怎么实现标注经纬度啊。这种lambert投影的标注我一致没能成功。

灭火器

独孤酒见 发表于 2025-3-20 10:57
灭火器大神这下应该把关于地图常用到的功能不但实现了，而且是便捷了。大神，在您的主贴中的效果图中怎么实 ...

GeoAxes.gridlines 就能标经纬度
https://scitools.org.uk/cartopy/ ... bels/gridliner.html
frykit 里的 set_map_ticks 函数也能，但很不完善，只能标注方形外框的地图。

独孤酒见

灭火器 发表于 2025-3-20 17:24
GeoAxes.gridlines 就能标经纬度
https://scitools.org.uk/cartopy/docs/latest/gallery/gridlines_and_ ...

感谢大神。我看到了你在网页上给的例子了。
https://zhajiman.github.io/post/frykit/
我准备尝试一下。里面的例子就是我想要的标注格式。
我现在经纬度的标注都标注在图里了，没标注在边框上。

感谢大神。不只是高人，还方便广大群众使用。

gdou_lyc

感谢大神提供的工具！！想问问大神，我南海小图里contouf填色是没问题的，但是用quiver画风矢量不显示，想问问大神这边有使用该库包画小图的风矢量的例子吗{:eb302:}

灭火器

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

复制代码

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