Базовая карта Matplotlib: удаление океана

Question

У меня есть meshplot, на который я хочу наложить континенты в базовой карте. Я использую этот код:

       from mpl_toolkits.basemap import Basemap
       import matplotlib.pyplot as plt
       m = Basemap(width=12000000,height=9000000,projection='lcc',
            resolution=None,lat_1=45.,lat_2=55,lat_0=50,lon_0=-107.)

       m.drawlsmask(land_color='coral',ocean_color='aqua',lakes=True)
       plt.show()

в отношении this , мои требования противоположны. Я хочу, чтобы континенты располагались поверх сетки или изображения, которое у меня есть, поэтому виден только я sh в районе океана.

Answer 1

Ваши требования:

1. Участок с континентами поверх сетки и / или изображения
2. Видна только сетка / изображение в районе океана

Чтобы получить график, вы должны использовать «zorder» в каждом из задействованных слоев. Данные для построения графика должны быть соответствующим образом преобразованы. Вот код, который вы можете попробовать, и выходной график, который он создает.

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np

m = Basemap(projection='lcc', width=12000000, height=9000000, 
    resolution='c', lat_1=45., lat_2=55, lat_0=50, lon_0=-107.)

# if resolution is None, coastlines wont draw

# draw only land areas with zorder=20
# any other layers with zorder below 20 will be hidden by land areas
m.drawlsmask(land_color='coral', ocean_color='none', lakes=True, zorder=20)
m.drawcoastlines(linewidth=0.3, color='gray', zorder=25)

filename = "small_01.png"  #use your image here
lonmin, lonmax, latmin, latmax = (-130, -40, 35, 45) # set limits of the image

# compute the limits of the image in data coordinates
left, bottom = m (lonmin, latmin)
top, right = m(lonmax, latmax)
image_extent = (left, right, bottom, top)

ax = plt.gca()
# set zorder < 20, to plot the image below land areas
ax.imshow(plt.imread(filename), extent=image_extent, zorder=15)

# plot some meshgrid data
# set zorder above image, but below land
xs = np.linspace(-130, -60, 20)
ys = np.linspace(20, 60, 10)
x2d, y2d = np.meshgrid(xs, ys)

#ax.plot(*m(x2d, y2d), 'ro', zorder=16)  # faster
ax.scatter(*m(x2d, y2d), s=2, zorder=16)

plt.show()

Редактировать 1

Некоторые полезные фрагменты кода:

# This plots shaded relief terrain covering land and sea.
m.shadedrelief(zorder = 25)

# This plots only ocean/sea parts on top.
m.drawlsmask(land_color='none', ocean_color='aqua', zorder=26)