消除 Matplotlib/Basemap pcolor 图中的白边答案

【问题标题】：Eliminate white edges in Matplotlib/Basemap pcolor plot消除 Matplotlib/Basemap pcolor 图中的白边
【发布时间】：2012-11-16 09:02:38
【问题描述】：

我正在使用此代码在地图上绘制数据：

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)

这会产生这个图像：

如您所见，网格单元之间存在白色间隙。我怎样才能摆脱它们？

【问题讨论】：

你能把数据文件放在任何地方吗？
不幸的是，没有，但它是规则间隔的网格数据..
无论如何，我无法使用随机生成的（或恒定的）数据重现您的问题...您使用的是什么版本的 matplotlib？也许它是一个已修复的错误？此外，pcolormesh 在这种情况下会更快。尝试使用它而不是pcolor。我怀疑（？）它会解决你的问题，但无论如何它应该更快。
pcolormesh 是否适用于极坐标投影？ (projection='nplaea') 我使用 pcolormesh 得到了一个非常奇怪的结果...（一切都为零）

标签： python numpy matplotlib matplotlib-basemap

【解决方案1】：

我知道这是一个老问题，但我想我会为这个问题添加我的解决方案。当我遇到与您的问题完全相同的问题时，我发现了您的问题，即我的绘图中有一条白线和一个从 -180 到 180 的网格。我的解决方案是使用 Basemap 函数addcyclic

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

这解决了我的问题。干杯，特隆德

【讨论】：

【解决方案2】：

我曾经遇到过同样的问题。问题很可能出在longitude。确保0 和360 都存在于输入中。如果没有，请手动添加它们，然后相应地更改mydata，使它们具有相同的形状。

【讨论】：

我使用的范围是 -180、180。这也应该可以吗？确实是少了180，所以我添加了它，并把最后一列的数据再次复制到数据中。但是，白色条纹仍然存在。

【解决方案3】：

在我看来，原来的帖子实际上并不是在询问 0 和 360 度之间的白色区域。

我认为 OP 正在谈论每个颜色方块之间的线条，这将与 this bug 一致：

似乎将 pcolor 绘图保存为 pdf 格式总是包含网格线，而对于其他输出格式（如 png）则不是这样

以下是开发人员对该问题的看法：

我在 gs、xpdf Preview.app 的结果图像中看到网格线，但没有在 Adobe 阅读器中。当我放大预览时，线条会跳动一点，并且无论缩放级别如何，屏幕上的宽度始终相同。

在这个例子中绘制的是很多多边形，所以相邻的多边形共享一条具有完全相同坐标的边。代码填写每个多边形的内部，显然还有一些渲染算法在多边形之间留下一条最小宽度的线。

所以这是 PDF 查看器的问题，而不是 pcolor 或 matplotlib 的任何其他方面。

【讨论】：