南极立体投影的 cartopy 问题

问题描述

我第一次使用 xarray。我的数据有维度（时间：1000，z_t：60，nlat：384，nlon：320）。我正在尝试绘制特定深度和时间的海洋温度。这在platecarree投影上是成功的，但是当我尝试切换到南极立体投影（使用cartopy）时，我的情节是错误的。南极洲的海岸线与数据的大陆轮廓不匹配。

首先，代码起作用的部分：

import numpy as np
import numpy.ma as ma
import xarray as xr
from matplotlib import pyplot as plt
filename = '/glade/scratch/mb/127kaH11/TEMP_annual_1000yr.nc'
ds = xr.open_dataset(filename)

接下来是我尝试新投影的部分，它似乎无法正常工作：

# try to actually plot on SH projection Now
# This is the map projection we want to plot *onto*
import cartopy.crs as ccrs

#choose a projection
proj = ccrs.southPolarstereo()
#ax.set_global()

ax = plt.subplot(111,projection=proj)

ds.isel(time=999,z_t=33).TEMP.plot.pcolormesh(ax = ax,transform=ccrs.PlateCarree())

ax.set_extent = ([-360,-90,-50],ccrs.PlateCarree())  

ax.coastlines()

我的另一个问题是第二张图片中的白色楔形。有人可以帮我缩小这个差距吗？过去我使用过 add_cyclic_point 但不确定这是否适用于此或如何正确使用它。

非常感谢。

Figure created from first snippet that works

Figure created from second snippet,that is causing problems

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）