在给定的坐标处添加饼图到cartopy投影

Question

我是一个数据可视化的初学者，甚至更多与cartopy，我知道大多数人我的问题是显而易见的。我正在尝试熟悉cartopy，并成功绘制了文字和点。但我无法实现饼图。

我只想绘制一个特定投影的饼图。但是我真的很迷惑，尽管有关于cartopy的文档。我已经先试试这个：

import cartopy.crs as ccrs 
import matplotlib.pyplot as plt 

ax = plt.axes(projection=ccrs.Robinson()) 
ax.coastlines(resolution='110m') # 110, 50, 10 
ax.stock_img() 

lat, long = 30, 30 # the latitude longitude 
ax.pie(long, lat, [0.25, 0.75], transform=ccrs.PlateCarree()) 

不工作，所以我检查，我发现这个Cartopy coastlines hidden by inset_axes use of Axes.pie，但我不明白的引擎盖下发生了什么，此外似乎仅限于PlateCarre（）。我试图修改它，但我没有设法让它正常工作。

所以我非常简单的问题是如何将几个饼图添加到给定经度和纬度的特定投影？如果你能发展你的答案，你将会非常受欢迎。

Answer 1

您可以使用inset_axes将新轴放置在图中，这将允许承载饼图。 inset_axes的位置由bbox_to_anchor参数决定。要使此参数使用坐标轴的投影坐标（ax），您需要设置bbox_transform=ax.transData。 如果您的坐标位于不同的坐标系中，则需要先使用投影的.transform_point方法将它们转换为使用的坐标。

import cartopy.crs as ccrs 
import matplotlib.pyplot as plt 
from mpl_toolkits.axes_grid1.inset_locator import inset_axes 

ax = plt.axes(projection=ccrs.Robinson()) 
ax.coastlines(resolution='110m') 
ax.stock_img() 


def plot_pie_inset(data,ilon,ilat,ax,width): 
    ax_sub= inset_axes(ax, width=width, height=width, loc=10, 
         bbox_to_anchor=(ilon, ilat), 
         bbox_transform=ax.transData, 
         borderpad=0) 
    wedges,texts= ax_sub.pie(data) 

    ax_sub.set_aspect("equal") 

lon,lat = 90,30 
lonr,latr = ccrs.Robinson().transform_point(lon,lat, ccrs.PlateCarree()) 
plot_pie_inset([0.25, 0.75],lonr,latr,ax,0.5) 


plt.show()