我正在尝试将网格拟合到海岸线地图,但无意中绘制了热图的旋转版本。我怀疑这可能是一个投影问题,但不确定如何解决。
领域专家(我从中收到数据)建议我使用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)