我有一个光栅文件,它由三个独立的2D-numpy数组(数据、经度、经度)和一个shapefile (.shp)组成,它们是从一个NetCDF文件导入的。我想使用shapefile作为掩码来计算掩码中所有数据点的平均值、中位数和标准差。在Python中做这件事最好的地方是什么? This is how it looks like on the map import matplotlib.pyplot as pltplume = geopandas.read_file(shapefile)ax = fig.add_sub