为什么我不能在Matplotlib / Cartopy地图上添加图像?

时间:2018-08-28 20:20:01

标签: python matplotlib gis cartopy

我正在从事一个项目,该项目涉及使用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)

结果:

如果我按原样运行代码,则地图的外观如下:

map with background

如果我将其更改为(即绕过绘制背景的函数):

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

然后我得到:

enter image description here

为什么我不能在地图顶部显示红色的“加号”?我尝试设置不同对象的“ 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

1 个答案:

答案 0 :(得分:2)

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