[英]Eliminate white edges in 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)
这导致此图像:
如您所见,网格单元之间存在白色间隙。 我怎么能摆脱他们?
我知道这是一个老问题,但我想我会为这个问题添加我的解决方案。 当我遇到与你的问题完全相同的问题时,我发现了你的问题,即我的情节中的白线和从-180到180的网格。我的解决方案是使用底图函数addcyclic
from mpl_toolkits.basemap import Basemap, shiftgrid, addcyclic
SSTcyclic, lonCMIP5cyclic = addcyclic(SST, lonCMIP5)
这解决了我的问题。 干杯,特隆德
我曾经遇到过同样的问题。 这个问题很可能是longitude
问题。 确保输入中存在0
和360
。 如果没有,请手动添加它们,并相应地更改mydata
,使它们具有相同的形状。
在我看来,原来的帖子实际上并没有询问0
到360
度之间的白色区域。
我认为OP正在讨论每个颜色方块之间的线条,这与这个bug是一致的:
看来,pdf格式的pcolor图总是包含网格线,但对于像png这样的其他输出格式则不然
以下是开发人员对此问题的看法:
我在gs,xpdf Preview.app中看到结果图像中的网格线,但在Adobe Reader中没有。 当我放大预览时,线条会跳跃一点,并且无论缩放级别如何,屏幕上的宽度始终相同。
在这个例子中绘制的是很多多边形,因此相邻的多边形共享一个具有完全相同坐标的边。 代码填充每个多边形的内部,显然一些渲染算法在多边形之间留下最小宽度线。
所以这是PDF查看器的问题,而不是pcolor
或matplotlib
任何其他方面。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.