Проекция базовой карты npstere возвращает пустую карту без данных

Я пытаюсь построить график концентрации морского льда над Северным полюсом, используя полярную стереографическую проекцию с базовой картой. На возвращенной карте выполнено все форматирование проекции, но нет данных. Я могу построить график с цилиндрической проекцией, используя те же данные, но когда я меняю проекцию на «npstere», она остается пустой. Я пробовал разные наборы данных, но карта пуста, независимо от того, какой набор данных я использую.

Вот созданная карта (скриншот, сохраненный график просто пустой)

Blank_map

Вот мой код, я использую Python3.5, Matplotlib 2.2.2 и Basemap 1.0.7.

# load packages
%matplotlib inline
import xarray as xr
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap

# load data
icedata = xr.open_dataset('IceData/G10010_SIBT1850_v1.1_JF.nc', decode_times=False)
seaice_1850 = icedata.seaice_conc_JF[0]

# make map projection
map_fig = Basemap(projection='npstere', boundinglat=50 ,lon_0=0)

# get lat/lon formatted
lon_1d = icedata.longitude.values
lat_1d = icedata.latitude.values
lon_2d, lat_2d = np.meshgrid(lon_1d, lat_1d)

# plot data over map
map_fig.pcolormesh(lon_2d, lat_2d, seaice_1850, cmap='jet')

# make it pretty
map_fig.drawcoastlines()
map_fig.drawparallels(np.arange(-80.,81.,10.))
map_fig.drawmeridians(np.arange(-180.,181.,20.))
plt.colorbar()
plt.show()

person Sarah Thunberg    schedule 12.11.2018    source источник
comment
Устанавливает ли ключевое слово latlon в значение True в pcolormesh (matplotlib.org/basemap/api/basemap_api.html) звонок изменить что-нибудь? Если нет, не могли бы вы предоставить свой файл данных?   -  person Patol75    schedule 13.11.2018
comment
Спасибо! Это сработало, у меня было ощущение, что что-то не так с координатами, но у меня недостаточно опыта, чтобы понять, что не так.   -  person Sarah Thunberg    schedule 13.11.2018


Ответы (1)


Установка ключевого слова latlon в True в вызове pcolormesh (matplotlib.org/basemap/api/basemap_api.html) делает свое дело.

person Patol75    schedule 14.11.2018