Исключить белые грани в Matplotlib/Basemap pcolor plot

Я рисую данные на карте, используя этот код:

import numpy as np
import matplotlib as mpl
import matplotlib.pyplot as plt
import matplotlib.cm as cm
from mpl_toolkits.basemap import Basemap
from matplotlib.patches import Polygon
from scipy.io import netcdf

ncfile = netcdf.netcdf_file(myfile.nc,'r')
lon = ncfile.variables['longitude'][:]
lat = ncfile.variables['latitude'][:]
data = ncfile.variables['mydata'][:]
ncfile.close()

m = Basemap(projection='nplaea', boundinglat=40, lon_0=270)
m.drawcoastlines(linewidth=.6, zorder=2)
m.drawparallels(np.arange(-80.,81.,20.), zorder=1)
m.drawmeridians(np.arange(-180.,181.,20.), zorder=1)
cNorm = mpl.colors.Normalize(vmin=0, vmax=np.nanmax(data))
cmap = plt.get_cmap('jet')
lons, lats = np.meshgrid(lon, lat)
x, y = m(lons, lats)
datamap = m.pcolor(x, y, data, zorder=0)
datamap.set_norm(cNorm) 
plt.colorbar(datamap, cmap=cmap, norm=cNorm, shrink=0.5)
plt.savefig('figures/map_polar.png', dpi=150, bbox_inches='tight', pad_inches=0.4)

В результате этого изображения: enter image description here

Как вы можете видеть, между ячейками сетки есть белые промежутки. Как я могу избавиться от них?

Ответы

Ответ 1

У меня была одна и та же проблема. Вероятно, проблема в longitude. Убедитесь, что 0 и 360 существуют на входе. Если нет, добавьте их вручную и внесите изменения в mydata соответственно, чтобы они имели одинаковую форму.

Ответ 2

Я знаю, что это старый вопрос, но я думал, что добавлю свое решение этой проблемы. Я нашел ваш вопрос, когда у меня была такая же проблема, как у вас, то есть белая линия в моем сюжете и сетка, идущая от -180 до 180. Решением для меня было использование функции Basemap addcyclic

from mpl_toolkits.basemap import Basemap, shiftgrid, addcyclic
SSTcyclic, lonCMIP5cyclic = addcyclic(SST, lonCMIP5)

Это решило мою проблему. Cheers, Trond

Ответ 3

Мне кажется, что исходное сообщение на самом деле не спрашивает о белой области между градусами 0 и 360.

Я думаю, что OP говорит о линиях между каждым квадратом цвета, который будет соответствовать этой ошибке:

Кажется, что savig pcolor plot в формате pdf всегда включает сетки, что неверно для других выходных форматов, таких как png

Вот что разработчики говорят о проблеме:

Я вижу сетчатые линии в результирующем изображении в gs, xpdf Preview.app, но не в Adobe Reader. Когда я увеличиваю предварительный просмотр, линии немного крутятся, и всегда имеют одинаковую ширину на экране, независимо от уровня масштабирования.

В этом примере получается много полигонов, так что смежные многоугольники разделяют ребро с точно такими же координатами. Код заполняется внутри каждого многоугольника и, по-видимому, некоторые алгоритмы рендеринга оставьте линию минимальной ширины между полигонами.

Так что это проблема с программой просмотра PDF, а не с pcolor или любым другим аспектом matplotlib.