【解决】cartopy 刻度不能显示右端的0度

sam_doggy

1. 
   
2. # ===================================================
   
3. 
   
4. # ===================================================
   
5. from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter
   
6. import cartopy.feature as cfeature
   
7. import cartopy.crs as ccrs
   
8. import matplotlib.pyplot as plt
   
9. scale = '110m'
   
10. fig = plt.figure(figsize=(8, 10))
    
11. ax = plt.axes(projection=ccrs.PlateCarree(central_longitude=180))
    
12. ax.set_global()
    
13. # ===================================================
    
14. 
    
15. land = cfeature.NaturalEarthFeature('physical', 'land', scale, edgecolor='face',
    
16.                                     facecolor=cfeature.COLORS['land'])
    
17. ax.add_feature(land, facecolor='0.75')
    
18. # ===================================================
    
19. 
    
20. ax.coastlines(scale)
    
21. # ===================================================
    
22. 
    
23. ax.set_xticks([0, 60, 120, 180, 240, 300, 360], crs=ccrs.PlateCarree())
    
24. ax.set_yticks([-90, -60, -30, 0, 30, 60, 90], crs=ccrs.PlateCarree())
    
25. 
    
26. lon_formatter = LongitudeFormatter(zero_direction_label=False)
    
27. lat_formatter = LatitudeFormatter()
    
28. ax.xaxis.set_major_formatter(lon_formatter)
    
29. ax.yaxis.set_major_formatter(lat_formatter)
    
30. 
    
31. gl = ax.gridlines()

复制代码

之前一直用basemap,刚准备转用cartoapy

我copy了@阿阿飞飞在 http://bbs.06climate.com/forum.p ... 3601&extra=page%3D1 中的代码，然而右端的0度没显示，进去看了geoaxes的源码，但是没找到问题....

刚刚进官网看示例也是没有右边的tick啊


有替代的解决方法如下：
https://stackoverflow.com/questi ... ertab=votes#tab-top

1. [0, 60, 120, 180, 240, 300, 359.9999999999]

复制代码

为什么阿飞的帖子里是ok 的我就不懂了

edwardli

你可以试试gridliner()这个方法。Cartopy map gridlines and tick labels，https://scitools.org.uk/cartopy/ ... ridliner.xformatter
import matplotlib.pyplot as plt
import matplotlib.ticker as mticker
import cartopy.crs as ccrs

from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER


ax = plt.axes(projection=ccrs.Mercator())
ax.coastlines()

gl = ax.gridlines(crs=ccrs.PlateCarree(), draw_labels=True,
                  linewidth=2, color='gray', alpha=0.5, linestyle='--')
gl.xlabels_top = False
gl.ylabels_left = False
gl.xlines = False
gl.xlocator = mticker.FixedLocator([-180, -45, 0, 45, 180])
gl.xformatter = LONGITUDE_FORMATTER
gl.yformatter = LATITUDE_FORMATTER
gl.xlabel_style = {'size': 15, 'color': 'gray'}
gl.xlabel_style = {'color': 'red', 'weight': 'bold'}

plt.show()

sam_doggy

edwardli 发表于 2019-6-2 05:35
你可以试试gridliner()这个方法。Cartopy map gridlines and tick labels，https://scitools.org.uk/cartop ...

我只有0在两边不显示，180在两边是显示的，我刚刚试了下这个方法，右边还是没有啊。这个应该是因为180有正负，0没有