python-3.x - 从 shapefile 中提取 netCDF4 数据？

Question

在 python 3.6 中，我导入了一个包含全球降水值的 netCDF4 文件。我还导入了一个包含科罗拉多河流域形状的 shapefile。我的目标是只能在我的 shapefile 中读取/提取降水数据。我查找了多个示例，但没有一个真正有帮助。

到目前为止，这是我的代码：

from netCDF4 import Dataset
import numpy as np
import geopandas as gpd

nc = Dataset('filename.nc')
long = nc.variables['lon'][:]
lati = nc.variables['lat'][:]
rainfall = nc.variables['precip'][:]

shapefile=gpd.read_file('filename.shp')

上面的代码没有错误消息。

Answer 1

哦，看，房子里的水文学家！;)

好吧，到目前为止，您还没有对代码做太多事情，您所做的只是将文件读入内存。

当我尝试执行相同的分析（仅使用 grib 文件）时，我找到了一个非常棒的 Python 库，它正是为此目的，称为RasterStats。

它支持使用 ndarray 栅格对象以及大多数 GDAL 支持的栅格文件类型（也必须是 netCDF！），并且它会生成您想要的东西。

有关更多信息，请参阅一本非常简洁的手册，如果您卡在某个地方，请告诉我！