从xarray数据集中选择日期列表



我有包含以下信息的xarray数据集:

print(ds)
>><xarray.Dataset>
Dimensions:        (lat: 90, lon: 720, time: 1200)
Coordinates:
* lon            array([-90.      , -89.057592, -88.115183, -87.172775, -86.230366, -85.287958,
-84.34555 , -83.403141, -82.460733, -81.518325, -80.575916, -79.633508,
-78.691099, -77.748691, -76.806283, -75.863874, -74.921466, -73.979058,
-73.036649, -72.094241, -71.151832, -70.209424, -69.267016, -68.324607,
-67.382199, -66.439791, -65.497382, -64.554974, -63.612565, -62.670157,
-61.727749, -60.78534 , -59.842932, -58.900524, -57.958115, -57.015707,
-56.073298, -55.13089 , -54.188482, -53.246073, -52.303665, -51.361257,
-50.418848, -49.47644 , -48.534031, -47.591623, -46.649215, -45.706806,
-44.764398, -43.82199 , -42.879581, -41.937173, -40.994764, -40.052356,
-39.109948, -38.167539, -37.225131, -36.282723, -35.340314, -34.397906,
-33.455497, -32.513089, -31.570681, -30.628272, -29.685864, -28.743455,
-27.801047, -26.858639, -25.91623 , -24.973822, -24.031414, -23.089005,
-22.146597, -21.204188, -20.26178 , -19.319372, -18.376963, -17.434555,
-16.492147, -15.549738, -14.60733 , -13.664921, -12.722513, -11.780105,
-10.837696,  -9.895288,  -8.95288 ,  -8.010471,  -7.068063,  -6.125654,
-5.183246,  -4.240838,  -3.298429,  -2.356021,  -1.413613,  -0.471204,
0.471204,   1.413613,   2.356021,   3.298429,   4.240838,   5.183246,
6.125654,   7.068063,   8.010471,   8.95288 ,   9.895288,  10.837696,
11.780105,  12.722513,  13.664921,  14.60733 ,  15.549738,  16.492147,
17.434555,  18.376963,  19.319372,  20.26178 ,  21.204188,  22.146597,
23.089005,  24.031414,  24.973822,  25.91623 ,  26.858639,  27.801047,
28.743455,  29.685864,  30.628272,  31.570681,  32.513089,  33.455497,
34.397906,  35.340314,  36.282723,  37.225131,  38.167539,  39.109948,
40.052356,  40.994764,  41.937173,  42.879581,  43.82199 ,  44.764398,
45.706806,  46.649215,  47.591623,  48.534031,  49.47644 ,  50.418848,
51.361257,  52.303665,  53.246073,  54.188482,  55.13089 ,  56.073298,
57.015707,  57.958115,  58.900524,  59.842932,  60.78534 ,  61.727749,
62.670157,  63.612565,  64.554974,  65.497382,  66.439791,  67.382199,
68.324607,  69.267016,  70.209424,  71.151832,  72.094241,  73.036649,
73.979058,  74.921466,  75.863874,  76.806283,  77.748691,  78.691099,
79.633508,  80.575916,  81.518325,  82.460733,  83.403141,  84.34555 ,
85.287958,  86.230366,  87.172775,  88.115183,  89.057592,  90.      ])
* lat            array([  0.  ,   1.25,   2.5 , ..., 356.25, 357.5 , 358.75])
* time           array([cftime.DatetimeNoLeap(0001-01-15 12:00:00),
cftime.DatetimeNoLeap(0001-02-14 00:00:00),
cftime.DatetimeNoLeap(0001-03-15 12:00:00), ...,
cftime.DatetimeNoLeap(1200-10-15 12:00:00),
cftime.DatetimeNoLeap(1200-11-15 00:00:00),
cftime.DatetimeNoLeap(1200-12-15 12:00:00)], dtype=object)
Data variables:
tas              (time, lat, lon) float32 ...

现在我想要特定日期的tas值,所以我列出了这样一个列表:

[array(cftime.DatetimeNoLeap(0500-01-14 00:00:00), dtype=object),
array(cftime.DatetimeNoLeap(0341-01-14 00:00:00), dtype=object),
array(cftime.DatetimeNoLeap(1147-01-14 00:00:00), dtype=object),
array(cftime.DatetimeNoLeap(0943-01-14 00:00:00), dtype=object),
array(cftime.DatetimeNoLeap(0011-01-14 00:00:00), dtype=object),
array(cftime.DatetimeNoLeap(0359-01-14 00:00:00), dtype=object),
array(cftime.DatetimeNoLeap(0705-01-14 00:00:00), dtype=object),
array(cftime.DatetimeNoLeap(0494-01-14 00:00:00), dtype=object),
array(cftime.DatetimeNoLeap(0204-01-14 00:00:00), dtype=object),
array(cftime.DatetimeNoLeap(0127-01-14 00:00:00), dtype=object),
array(cftime.DatetimeNoLeap(0386-01-14 00:00:00), dtype=object)]

然后尝试使用这样的sel语句:

ds.sel(time = dates)

但它不起作用。有人能告诉我如何选择这样具体的日期吗。

这种方法适用于我的时间为datetime64格式的情况。

它看起来像是数组中列表中的每个项。您希望每个项目都是该项目。所以你的dates列表应该是这样的:

[cftime.DatetimeNoLeap(0500-01-14 00:00:00),
cftime.DatetimeNoLeap(0341-01-14 00:00:00),
...

这样行吗?

最新更新