将pcolormesh拟合到地图Cartopy

时间:2019-07-10 08:12:58

标签: python matplotlib data-visualization weather cartopy

我正在尝试将网格拟合到海岸线地图,但无意中绘制了热图的旋转版本。我怀疑这可能是一个投影问题,但不确定如何解决。

领域专家(我从中收到数据)建议我使用Lambert Conformal投影,但这似乎可以产生与Plate Carree相同的旋转网格。

import cartopy

# (time, X_grid, Y_grid)
lat = atmospheric.variables['XLAT'][0, :, :]
lon = atmospheric.variables['XLONG'][0, :, :]
extent = [np.min(lon.values), np.max(lon.values), np.min(lat.values), np.max(lat.values)]

projection_crs = ccrs.LambertConformal(central_longitude=np.mean(lon.values), central_latitude=np.mean(lat.values))

ax = plt.axes(projection=projection_crs)
ax.set_extent(extent, crs=projection_crs)

# (time, layer, X_grid, Y_grid)
P_HYD = atmospheric.variables['P_HYD'][0, 0, :, :]
ax.coastlines(resolution='50m', zorder=1)

im = ax.pcolormesh(lon, lat, P_HYD, transform=projection_crs, zorder=0)

enter image description here

0 个答案:

没有答案