我想在 Windows 上使用 Python 使用 NetCDF4 文件中的经度和纬度范围来选择数据

Question

我希望将其转换为 DataFrame 会有所帮助。我计划在生成的数据帧上使用 Between 函数。但转换本身不起作用，即使使用 xarray 也是如此。我还有另一个原因想要将其转换为 DataFrame：使用 statsmodel。

由于我使用的是 Windows，因此无法使用 nctoolkit。

#imports
import requests
import netCDF4 as nc
import gzip
#
# Download the gzipped datafile from GISS and load it into a netCDF Dataset
GISTEMPfile = requests.get('https://data.giss.nasa.gov/pub/gistemp/gistemp250_GHCNv4.nc.gz')
ds = nc.Dataset("dummy_path", mode="r", memory=gzip.decompress(GISTEMPfile.content))
ds.variables
v = ds.variables['tempanomaly']
#
#from pandas import DataFrame
#df = DataFrame(ds)
#
import xarray as xr
temp = xr.open_dataset(ds)
df = temp.to_dataframe()

Answer 1

Xarray 支持使用

sel

方法和切片对象进行基于范围的选择。例如：

ds = xr.open_dataset(...)
ds_region = ds.sel(lon=slice(-85.0, -80.1), lat=slice(-10.0, 10.0))