为什么我不能在我的 Matplotlib/Cartopy 地图上添加图像？答案

【问题标题】：Why can't I add an image on top of my Matplotlib/Cartopy map?为什么我不能在我的 Matplotlib/Cartopy 地图上添加图像？
【发布时间】：2018-08-28 20:20:00
【问题描述】：

我正在开展一个项目，该项目涉及使用 Cartopy 在地图上绘制数据。

到目前为止一切正常，但我一直在重构代码以使程序的其他部分可以调用不同的函数。所以为此，我有一个函数可以将我的背景添加到地图中，另一个函数可以在指定的纬度/经度处添加一个地标。显然，我想要背景上方的地标，但我似乎无法让它工作。

对于背景，我希望能够使用 Cartopy 库存图片或网络地图图块。无论哪种方式，问题都是相同的，因此我将 Cartopy 背景用于当前测试目的。这是那个函数：

def custom_background(self, source_point):
    cartmap = self.plot
    source_point = source_point.split(" ")
    source_point = (float(source_point[0]), float(source_point[1]))
    dx = 2.5
    dy = 5
    pad = 0.5
    lon_min, lon_max = source_point[0]-dx, source_point[0]+dx
    lat_min, lat_max = source_point[1]-dy, source_point[1]+dy
    area = 4*dx*dy
    zoom = self.get_zoom(area) ##only relevant when using a map tile
    cartmap.set_extent([lat_min-pad, lat_max+pad, lon_min-pad, lon_max+pad])
    #~ cartmap.add_image(self.tile, zoom)
    cartmap.add_feature(cartopy.feature.LAND, zorder=1)
    return cartmap

这里是地标功能：

def add_point_icon(self, x, y, cartmap):
    src_point = np.array(Image.open('icons/icon63.png'))
    im = OffsetImage(src_point, zoom=1, alpha=1.0, zorder=3)
    ab = AnnotationBbox(im, (x,y), xycoords='data', frameon=False)
    cartmap.add_artist(ab)

这两个都像这样一个接一个地调用：

    cartmap = self.custom_background(mysrc)
    #~ cartmap=self.plot
    self.add_point_icon(x1, y1, cartmap)

结果：

如果我按原样运行代码，这就是地图的外观：

如果我将其更改为（即绕过绘制背景的函数）：

    #~ cartmap = self.custom_background(mysrc)
    cartmap=self.plot
    self.add_point_icon(x1, y1, cartmap)

然后我得到：

为什么我不能让红色的“加号”符号显示在地图顶部？我尝试设置不同对象的“zorder”参数，但它似乎没有做任何事情。我现在完全不知所措。任何帮助将不胜感激，谢谢。

编辑：也许我还应该包括创建子图的行：

def __init__(self, mylevs):
    self.fig, self.header, self.footer, self.plot, self.legend = 
        self.create_spec()

def create_spec(self):
    """Define layout of figure"""

    #left column: header, footer, plt
    fig = plt.figure(figsize=(12,10))
    layout = 1
    if layout == 1: #Default
        widths = [8,1]
        heights = [2, 10, 3]
        column_border = 0.75
        pad = 0.1
        colorbar_width = 0.05
        spec = gridspec.GridSpec(ncols=1, nrows=3, width_ratios = [1], height_ratios=heights, left=0.1, right = column_border)
        #right column: colorbar
        spec2 = gridspec.GridSpec(ncols=1, nrows=1, width_ratios = [1], height_ratios=[1], left=column_border+pad, right=column_border+pad+colorbar_width)

        header = plt.subplot(spec[0,0])
        footer = plt.subplot(spec[2,0])
        plot = plt.subplot(spec[1,0], projection=cimgt.OSM().crs)
        legend = plt.subplot(spec2[0,0])

    return fig, header, footer, plot, legend

【问题讨论】：

我想知道是否正在绘制标志，但在您制作的地图范围之外？这会使它看起来没有被绘制。要快速检查这一点，您可以将行 cartmap.set_extent([...]) 替换为 cartmap.set_global()。如果该标志出现在地图上的其他位置，那么您需要将标志的(x, y) 点值转换为地图的投影（在add_point_icon 方法中执行此操作的好地方）。
就是这样。不是我需要转换点，而是我需要更改在创建轴时使用的投影。

标签： python matplotlib gis cartopy

【解决方案1】：

问题在于，在从地图图块切换回 Cartopy 库存背景时，我忘记将投影从 ccrs.OSM() 切换回 ccrs.PlateCarree()。在地图瓦片投影中，地标被绘制在查看窗口之外。

【讨论】：