累积时间

问：

我有一个 nc 文件，其中有 'lon'、'lat'、'time' 和 'foot'（代表脚印）。对于脚变量，我有脚（纬度，纬度，时间）。我做了一个 48 小时的前进轨迹高跷模拟，这意味着以脚为单位的时间包含 49 个不同的“切片”。我现在正试图积累这些时间，这样我就可以拥有一张 2D 地图。我不知道从哪里开始。以下是我的代码：

footprints = "myfile.nc"
print(footprints)
 
with Dataset(footprints) as root:
    time = root.variables['time'][:]
    dates = num2date(time, root.variables['time'].units)
    print(dates[0].strftime('%Y-%m-%d-%H-%M-%S'))
    lats= root.variables['lat'][:]
    lons= root.variables['lon'][:]
    foot= root.variables['foot'][:, :, :]
 
Foot = np.arange(0,49)
for i in Foot:
   ax = plt.axes(projection=ccrs.PlateCarree())
   plt.contourf(lons[:], lats[:], np.squeeze(foot[i,:,:]), cmap = 'YlOrBr', transform=ccrs.PlateCarree())
   ax.set_extent([-178.546, -130.946, 50, 75])   # regional map, x=longitude, y=latitude (xmin, xmax, ymin, ymax)  
   ax.coastlines()
   ax.gridlines(draw_labels=True)
   plt.show()

我不想只映射一个切片/图层，而是将它们全部映射到一个映射中。我应该得到一张地图，但我得到了 49 张。 你能帮忙吗？

当我用xr.open_dataset打开它时，脚如下所示：