底图Lambert Azimuthal等面积图:奇怪的人工制品

时间:2015-08-05 06:30:39

标签: projection matplotlib-basemap

我正在尝试使用以下代码来绘制太平洋半球的网格化数据;然而,正如下面的代码所示,似乎有一个奇怪的神器。我想知道是否有人能够展示我在设置方面出错的地方。谢谢。

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')

output

0 个答案:

没有答案