如何使用Python从NETCDF选择变量的时间范围

Question

我每天都有一个包含单个LAT和LON的点数据。 Vairable（var.shape）的形状是（'time'，'z'），这里z是海洋中的深度。我想在指定的日期范围内（2019-01-01至2020-12-31）和深度范围（100至1500）选择和绘制可变性。我知道如何选择深度范围，但时间不能做同样的时间。我如何选择一定范围的时间？

这里没有工作

import warnings
warnings.filterwarnings('ignore')

import os
import numpy as np
import xarray as xr
import matplotlib.pyplot as plt
import matplotlib as mpl

path_in = "./"
file_ds = xr.open_dataset(path_in + 'ocean.nc')

lon = file_ds.variables['lon'][:]
lat = file_ds.variables['lat'][:]
li = file_ds.variables['depth'][0,:,0,0] 
time = file_ds.variables['time'][:]
var = file_ds.variables['var'][:,:,0,0]

index = np.intersect1d(np.where(li<-100),np.where(li>=-1500))
tindex = np.intersect1d(np.where(time<(['2020-01-01'],dtype='datetime64[ns]'),np.where(time>=['2020-12-31'],dtype='datetime64[ns]'))

var_100_1500 = var[tindex,index] 

fig, ax = plt.subplots(figsize=(15, 9))
fillprof = ax.contourf(time[tindex],li[index],var_100_1500.T,levels=40)

clb=plt.colorbar(fillprof, orientation="vertical", pad=0.02)
clb.ax.tick_params(labelsize=18) 

我尝试了以下代码，但是时间的索引在

  tindex = np.intersect1d(np.where(time<(['2020-01-01'],dtype='datetime64[ns]'),np.where(time>=['2020-12-31'],dtype='datetime64[ns]'))
                                                               ^
SyntaxError: invalid syntax

。

Answer 1

您在tindex分配行的末尾缺少一个“）”。

看起来应该这样：

tindex = np.intersect1d(np.where(time<(['2020-01-01'], dtype='datetime64[ns]'), np.where(time>=['2020-12-31'], dtype='datetime64[ns]')))