cartopy画GMT边框风格的中国地图

灭火器

模仿 GMT 或者说 M_Map 的黑白相间边框风格画中国地图。用到的自制包见 https://github.com/ZhaJiMan/frykit
原理很简单，就是找出刻度的位置，然后在刻度之间画黑白条条。不过目前只支持 PlateCarree 投影的地图，不知道怎么在兰伯特投影那种扇形图里整。

好羡慕隔壁 GMT 的出图效果啊

1. '''
   
2. time: 2024-01-04
   
3. author: 灭火器
   
4. data: https://www.naturalearthdata.com/downloads/50m-raster-data/50m-cross-blend-hypso/
   
5. '''
   
6. import numpy as np
   
7. import matplotlib.pyplot as plt
   
8. import cartopy.crs as ccrs
   
9. import frykit.plot as fplt
   
10. 
    
11. im = plt.imread('./HYP_50M_SR_W/HYP_50M_SR_W.tif')
    
12. 
    
13. crs = ccrs.PlateCarree()
    
14. fig = plt.figure()
    
15. ax = fig.add_subplot(projection=crs)
    
16. ax.imshow(im, extent=[-180, 180, -90, 90], transform=crs)
    
17. fplt.add_cn_province(ax, lw=0.4)
    
18. fplt.add_nine_line(ax, lw=0.4)
    
19. fplt.set_extent_and_ticks(
    
20.     ax=ax,
    
21.     extents=[70, 140, 0, 60],
    
22.     xticks=np.arange(-180, 181, 10),
    
23.     yticks=np.arange(-90, 91, 10)
    
24. )
    
25. ax.tick_params(labelsize='small')
    
26. 
    
27. fplt.gmt_style_frame(ax)
    
28. fplt.add_compass(ax, 0.92, 0.88, size=20, style='circle')
    
29. map_scale = fplt.add_map_scale(ax, 0.05, 0.1)
    
30. map_scale.set_xticks([0, 500, 1000])
    
31. map_scale.set_xticks([250, 750], minor=True)
    
32. fplt.gmt_style_frame(map_scale, width=4)
    
33. 
    
34. fig.savefig('./map.png', dpi=300, bbox_inches='tight')

复制代码

XiaoMaFenJu

确实还挺好看的你别说

独孤酒见

灭火器大神又出牛作了，感谢大神分享。

独孤酒见

灭火器大神又出牛作了，感谢大神分享。
好像漏了个图例。

灭火器

独孤酒见 发表于 2024-1-5 10:30
灭火器大神又出牛作了，感谢大神分享。
好像漏了个图例。

首先我不是大神。图例的话，这个数据是Natural Earth上的浮雕地图，是个RGB数组，不是我用colormap和normalize映射出来的，所以没有colorbar之类的。