Устранить белые края в pplot Matplotlib/Basemap

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

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)

Это приводит к этому изображению: введите описание изображения здесь

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

Источник

3 ответа

Решение

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

Источник

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

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

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

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

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

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

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

Источник

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

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

Это решило мою проблему. Ура, Тронд

Источник
Другие вопросы по тегам