使用NorthPolarStereo绘制带有栅格的网格化数据

Question

我正在从底图移动到Cartopy，并希望绘制覆盖极点的北冰洋的数据。

我决定使用NorthPolarStereo（）投影，很高兴使用pcolormesh或轮廓线f。不幸的是，当我执行以下代码时，我的数据字段没有显示：

import cartopy.crs as ccrso
from netCDF4 import Dataset


def import_envisat_field(year,
                         month):

    data_dir = f'/media/robbie/Seagate Portable Drive/Envisat_thickness/{year}/'
    file = f'ESACCI-SEAICE-L3C-SITHICK-RA2_ENVISAT-NH25KMEASE2-{year}{month}-fv2.0.nc'
    data = Dataset(data_dir+file)

    return(data)

# Import data

data = import_envisat_field("2003","02")

# Make plot

fig = plt.figure(figsize=[10, 5])
ax = plt.axes(projection=ccrs.NorthPolarStereo())

ax.add_feature(cartopy.feature.OCEAN, zorder=0)
ax.add_feature(cartopy.feature.LAND, zorder=1, edgecolor='black')

extent = 2500000

ax.set_extent((-extent,
               extent,
               -extent,
               extent),
              crs=ccrs.NorthPolarStereo())

ax.gridlines()

lon = np.array(data['lon'])
lat = np.array(data['lat'])
field = np.array(data['sea_ice_thickness'])[0]

print(lon.shape,lat.shape,field.shape)
# This print command gives (432, 432) (432, 432) (432, 432)

plt.pcolormesh(lon, lat, field,zorder=2,
             transform=ccrs.NorthPolarStereo())

plt.show()

数据使用底图以一种直接的方式绘制，但是执行上面的代码只是给我一幅很好的北冰洋图，而没有我的数据。

我也尝试过用ax.pcolormesh替换plt.pcolormesh，但是那也不起作用。

Cartopy输出：

具有相同数据的底图输出：

Answer 1

如果您的数据坐标是经度和纬度，则需要使用PlateCarree变换：

plt.pcolormesh(lon, lat, field,zorder=2, transform=ccrs.PlateCarree())

变换描述了数据坐标，并且独立于要绘制的投影。有关更多详细信息，请参见Cartopy文档中的本指南。https://scitools.org.uk/cartopy/docs/latest/tutorials/understanding_transform.html