我正在从事一个项目,该项目涉及使用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
答案 0 :(得分:2)
问题在于从地图图块切换回Cartopy库存背景时,我忘记了将投影从ccrs.OSM()切换回ccrs.PlateCarree()。在地图图块的投影中,地标被绘制在查看窗口的外部。