python-3.x - 如何知道切片的 xarray Dataset/DataArray 是否为空？

Question

我正在使用 Python 语言对 Netcdf 文件应用切片和聚合操作。处理此类文件的解决方案之一是使用 Xarray 库。

我还是库功能的新手，所以我想知道 Xarray 对象是否有一些方法来检查切片的 DataSet/DataArray 是否为空，就像 Pandas 一样(在 Pandas 的情况下，可以检查是否数据框/系列通过“空”方法为空)。

我找到的唯一解决方案是始终将 Xarray Dataset/DataArray 转换为 Pandas Dataframe/Series，然后检查它是否为空。

以下是代码片段作为示例:

import xarray as xr

path = 'my_path_to_my_netcdf_file.nc'

Xarray_DataArray = xr.open_dataset(path)

print(Xarray_DataArray)

# this returns something like:

 #     Dimensions:      (lat: 600, lon: 672, time: 37)
 #     Coordinates:
 #     * lat          (lat) float32 -3.9791672 -3.9375012 ... 20.9375 20.979166
 #     * lon          (lon) float32 -60.979168 -60.9375 ... -33.0625 -33.020832
 #     * time         (time) datetime64[ns] 2010-05-19 2010-05-20 ... 2010-06-24
 #     Data variables:
 #       variable_name  (time, lat, lon) float32 dask.array<shape=(37, 600, 672), 
 #         chunksize=(37, 600, 672)>

 # I normally use the 'sel' method to slice the xarray object, like below:

Sliced_Xarray_DataArray = Xarray_DataArray.sel({'lat':slice(-10, -9),
                                                'lon':slice(-170, -169)                  
                                                })


 # but since, Xarray does not possess a proper way to check the slice, I usually have to do the following:

 if Sliced_Xarray_DataArray.to_dataframe().empty():
    print('is empty. Nothing to aggregate')

 else:
    Aggregated_value =  Aggregation_function(Sliced_Xarray_DataArray)

    print('continuing with the analysis')


 #    ... continue

我将不胜感激任何建议。

感谢您抽出宝贵时间，希望能尽快收到您的来信。

您忠诚的，

菲利普·R·莱尔

Answer 1

您可以很容易地检查结果切片中变量的大小是否为 0:

print(Sliced_Xarray_DataArray.time.size)
if Sliced_Xarray_DataArray.time.size == 0:
    print('is empty. Nothing to aggregate')
else:
    print('not empty. Go aggregate')

您的任何坐标变量以及我们的其他变量都可以作为 Sliced_Xarray_DataArray 中的属性访问。所以在你的例子中你可以检查 lat 的大小, lon或 time .