我正在尝试使用以下代码来绘制太平洋半球的网格化数据;然而,正如下面的代码所示,似乎有一个奇怪的神器。我想知道是否有人能够展示我在设置方面出错的地方。谢谢。
LON_0 = -160
LAT_0 = 20
plt.close()
fig = plt.figure(figsize=(11.6,11.6))
ax = plt.subplot()
lats = np.linspace(90,-90,181)
lons = np.linspace(180,-180,361)
lons, lats = np.meshgrid(lons,lats)
m = Basemap(width=18000000, height=18000000,
lat_ts=LAT_0, lat_0=LAT_0, lon_0=LON_0,
ax=ax, projection='laea',
resolution='c')
m.drawcoastlines()
data =np.sin(lons*5*math.pi/180.)**2+np.sin(lats*5*math.pi/180.)**2
shiftedLons,shiftedData = m.shiftdata(lons, data,
lon_0=LON_0)
px,py = m(shiftedLons, lats)
cbinfo = m.pcolormesh(px, py, shiftedData, cmap='rainbow')