获取基于shapefile掩码的栅格图像值

Question

我有一个栅格文件，它由三个独立的2D-numpy数组(data, lat, lon)和一个shapefile (.shp)组成，它们是从NetCDF文件导入的。我想使用shapefile作为掩码来计算掩码内所有数据点的平均值，中位数和标准差。在Python中做这个最好的是什么?

这是地图上的样子

import matplotlib.pyplot as plt
import geopandas
plume = geopandas.read_file(shapefile)
fig = plt.figure()
ax = fig.add_subplot(111, projection=ccrs.PlateCarree)
f = ax.pcolormesh(lons, lats, np.ma.masked_invalid(data), shading='auto', transform=ccrs.PlateCarree())
plume.boundary.plot(ax=ax)
plt.show()

Answer 1

查看栅格库。用于处理栅格数据的Geopandas (Vectordata)对应项。处理屏蔽以及对数据的计算。