我在北极(北极)上使用 ASI 海冰浓度并使用 NorthPolarStereo 进行地图投影,但我找不到合适的转换。实际上,一旦我将其绘制出来,数据就不适合地图了!注意:此 NSIDC 投影应基于 WGS 1984 椭球 (EPSG3413)。+proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
这是代码(x,y,ice 是 lon、lat、海冰浓度):
import cartopy.crs as ccrs
ax = plt.axes(projection=ccrs.NorthPolarStereo(central_longitude=0))
kw = dict(central_latitude=90, central_longitude=-45, true_scale_latitude=70)
cs = ax.pcolormesh(y, x, ice, cmap='cividis', transform=ccrs.Stereographic(**kw))