微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

南极立体投影的 cartopy 问题

如何解决南极立体投影的 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 举报,一经查实,本站将立刻删除。