我正在使用以下代码在 map 上绘制数据:
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 的网格。我的解决方案是使用 basemap 函数 addcyclic
from mpl_toolkits.basemap import Basemap, shiftgrid, addcyclic
SSTcyclic, lonCMIP5cyclic = addcyclic(SST, lonCMIP5)
这解决了我的问题。干杯,特隆
关于python - 消除 Matplotlib/Basemap pcolor 图中的白边,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/13606304/