用matplotlib移动网格

Question

我正在绘制格式在0到360度之间的数据。我试图在cyl或merc投影上绘制这个，但它只显示从0开始的数据（我想在中心绘制GMT数据，所以需要-180到180的lon网格上的数据） 。如果我移动网格（lon = lon -180）然后显示所有数据，但数据位于错误的位置-180度。

问题：

虽然在正射投影中工作正常。相关代码如下。

lat = np.linspace(90,-90,721)
lon = np.linspace(0,360,1440) 
m = Basemap(projection='cyl',llcrnrlat=-90,urcrnrlat=90,llcrnrlon=0,urcrnrlon=360,resolution='c',)
X, Y = np.meshgrid(lon, lat)
X, Y = m(X, Y)
cs = m.contourf(X,Y,Plot,scale, cmap=cmap)

Answer 1

请尝试：

import numpy as np
from mpl_toolkits.basemap import shiftgrid
from mpl_toolkits.basemap import Basemap

lat = np.linspace(-90, 90, 721)
lon = np.linspace(0, 360, 1440)
Plot, lon = shiftgrid(180., Plot, lon, start=False)  # shiftgrid
m = Basemap(projection='cyl', llcrnrlat=-90, urcrnrlat=90, llcrnrlon=-180, urcrnrlon=180, resolution='c',)
X, Y = np.meshgrid(lon, lat)
X, Y = m(X, Y)
cs = m.contourf(X, Y, Plot, scale, cmap=cmap)

shiftgrid：将全局纬度/经度网格向东或向西移动。

Answer 2

我有一个解决方案（尽管很难看）。通过重新排序数据。

temp = np.zeros((721,1440))
temp[:,0:720] = Plot[:,720:1440]
temp[:,720:1440] = Plot[:,0:720]    
Plot[:]=temp[:]

或者使用np.roll（如果知道要移动多少个网格点）