Установка pcolormesh для картографирования - PullRequest
0 голосов
/ 10 июля 2019

Я пытаюсь разместить сетку на карте береговой линии, но случайно строю повернутую версию тепловой карты. Я подозреваю, что это может быть проблема проецирования, но я не уверен, как мне решить ее.

Эксперты в области (от которых я получил данные) рекомендовали, чтобы я использовал конформную проекцию Ламберта, но это, похоже, дает ту же повернутую сетку, что и Плато Карри.

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

...